Linked vs unlinked markers: multilocus microsatellite haplotype-sharing as a tool to estimate gene flow and introgression

We have explored the use of multilocus microsatellite haplotypes to study introgression from cultivated (Malus domestica) into wild apple (Malus sylvestris), and to study gene flow among remnant populations of M. sylvestris. A haplotype consisted of alleles at microsatellite loci along one chromosome. As destruction of haplotypes through recombination occurs much faster than loss of alleles due to genetic drift, the lifespan of a multilocus haplotype is much shorter than that of the underlying alleles. When different populations share the same haplotype, this may indicate recent gene flow between populations. Similarly, haplotypes shared between two species would be a strong signal for introgression. As the expected lifespan of a haplotype depends on the strength of the linkage, the length [in centiMorgans (cM)] of the haplotype shared contains information on the number of generations passed. This application of shared haplotypes is distinct from using haplotype-sharing to detect association between markers and a certain trait. We inferred haplotypes for four to eight microsatellite loci on Linkage Group 10 of apple from genotype data using the program phase, and then identified those haplotypes shared between populations and species. Compared with a Bayesian analysis of unlinked microsatellite loci using the program structure, haplotype-sharing detected a partially different set of putative hybrids. Cultivated haplotypes present in M. sylvestris were short (< 1.5 cM), indicating that introgression had taken place many generations ago, except for two Belgian plants that contained a haplotype of 47.1 cM, indicating recent introgression. In the estimation of gene flow, F(ST) based on unlinked loci indicated small (0.032-0.058) but statistically significant differentiation between some populations only. However, various M. sylvestris haplotypes were shared in nearly all pairwise comparisons of populations, and their length indicated recent gene flow. Hence, all Dutch populations should be considered as one conservation unit. The added value of using sharing of multilocus microsatellite haplotypes as a source of population genetic information is discussed.     

High nucleotide diversity and limited linkage disequilibrium in Helicoverpa armigera facilitates the detection of a selective sweep

Insecticides impose extreme selective pressures on populations of target pests and so insecticide resistance loci of these species may provide the footprints of ‘selective sweeps''. To lay the foundation for future genome-wide scans for selective sweeps and inform genome-wide association study designs, we set out to characterize some of the baseline population genomic parameters of one of the most damaging insect pests in agriculture worldwide, Helicoverpa armigera. To this end, we surveyed nine Z-linked loci in three Australian H. armigera populations. We find that estimates of π are in the higher range among other insects and linkage disequilibrium decays over short distances. One of the surveyed loci, a cytochrome P450, shows an unusual haplotype configuration with a divergent allele at high frequency that led us to investigate the possibility of an adaptive introgression around this locus.     

Migration patterns among global populations of the pathogenic fungus Mycosphaerella graminicola

DNA sequences from five nuclear loci and data from three microsatellites were collected from 360 isolates representing 14 globally distributed populations of the plant pathogenic fungus Mycosphaerella graminicola. Haplotype networks were constructed for the five sequence loci and population subdivision was assessed using Hudson's permutation test. Migration estimates were calculated using six regional populations for both the sequence and microsatellite loci. While subdivision was detected among the six regional populations, significant gene flow was indicated among some of the populations. The European and Israeli populations contributed the majority of historical immigrants to the New World. Migration estimates for microsatellite loci were used to infer more recent migration events among specific New World populations. We conclude that gene flow was an important factor in determining the demographic history of Mycosphaerella graminicola.     

Detecting introgressive hybridization between free-ranging domestic dogs and wild wolves (Canis lupus) by admixture linkage disequilibrium analysis

Occasional crossbreeding between free-ranging domestic dogs and wild wolves (Canis lupus) has been detected in some European countries by mitochondrial DNA sequencing and genotyping unlinked microsatellite loci. Maternal and unlinked genomic markers, however, might underestimate the extent of introgressive hybridization, and their impacts on the preservation of wild wolf gene pools. In this study, we genotyped 220 presumed Italian wolves, 85 dogs and 7 known hybrids at 16 microsatellites belonging to four different linkage groups (plus four unlinked microsatellites). Population clustering and individual assignments were performed using a Bayesian procedure implemented in structure 2.1, which models the gametic disequilibrium arising between linked loci during admixtures, aiming to trace hybridization events further back in time and infer the population of origin of chromosomal blocks. Results indicate that (i) linkage disequilibrium was higher in wolves than in dogs; (ii) 11 out of 220 wolves (5.0%) were likely admixed, a proportion that is significantly higher than one admixed genotype in 107 wolves found previously in a study using unlinked markers; (iii) posterior maximum-likelihood estimates of the recombination parameter r revealed that introgression in Italian wolves is not recent, but could have continued for the last 70 (+/- 20) generations, corresponding to approximately 140-210 years. Bayesian clustering showed that, despite some admixture, wolf and dog gene pools remain sharply distinct (the average proportions of membership to wolf and dog clusters were Q(w) = 0.95 and Q(d) = 0.98, respectively), suggesting that hybridization was not frequent, and that introgression in nature is counteracted by behavioural or selective constraints.     

SimPed: a simulation program to generate haplotype and genotype data for pedigree structures

With the widespread availability of SNP genotype data, there is great interest in analyzing pedigree haplotype data. Intermarker linkage disequilibrium for microsatellite markers is usually low due to their physical distance; however, for dense maps of SNP markers, there can be strong linkage disequilibrium between marker loci. Linkage analysis (parametric and nonparametric) and family-based association studies are currently being carried out using dense maps of SNP marker loci. Monte Carlo methods are often used for both linkage and association studies; however, to date there are no programs available which can generate haplotype and/or genotype data consisting of a large number of loci for pedigree structures. SimPed is a program that quickly generates haplotype and/or genotype data for pedigrees of virtually any size and complexity. Marker data either in linkage disequilibrium or equilibrium can be generated for greater than 20,000 diallelic or multiallelic marker loci. Haplotypes and/or genotypes are generated for pedigree structures using specified genetic map distances and haplotype and/or allele frequencies. The simulated data generated by SimPed is useful for a variety of purposes, including evaluating methods that estimate haplotype frequencies for pedigree data, evaluating type I error due to intermarker linkage disequilibrium and estimating empirical p values for linkage and family-based association studies.     

Correction method for null alleles in species with variable microsatellite flanking regions, a case study of the black-lipped pearl oyster Pinctada margaritifera

In bivalves, heterozygote deficiencies and departures from Hardy-Weinberg equilibrium (HWE) in microsatellite analysis are common and mainly attributed to inbreeding, genetic patchiness (Walhund effect), or null alleles. We checked for the occurrence of null alleles at 3 microsatellite loci in 3 populations of black-lipped pearl oyster, Pinctada margaritifera, using a step-by-step method to re-amplify homozygotes and null individuals with redesigned primer pair combinations. After amplification with original primer pairs, the 3 populations exhibited null alleles, absence of structure, and significant departure from HWE for all 3 loci due to heterozygote deficiencies. After 3 re-amplification steps, with modified primer sets, all loci were corrected for null alleles. Once corrected, all populations appeared at HWE, demonstrating that null alleles were responsible for the initial disequilibrium of the populations. Furthermore, analysis from corrected genotypes demonstrates significant genetic differentiation for one population from the other 2.     

Linkage Disequilibrium between Allozymes in Natural Populations of Lodgepole Pine

Pairwise linkage disequilibrium values (D) were estimated for 14 allozyme loci in two natural populations of lodgepole pine (Pinus contorta ssp. latifolia). Maternal multilocus genotypes were inferred from samples of (haploid) megagametophytic seed-endosperms. Coupling/repulsion double heterozygotes were distinguished for closely linked pairs of loci. Assays of seven of the loci in seed embryos allowed estimates of D for these loci in the outcross pollen pool (estimates of outcrossing rates indicate no significant departures from random mating in either population). No disequilibrium was observed between unlinked loci in either maternal genotypes or outcross pollen. However, significant disequilibrium was observed within and between gametes for some allelic combinations of four tightly linked loci; the assumption of random association of gamete types within individuals is thus invalid for some loci in lodgepole pine. Possible causes of the observed D were examined using the noncentrality parameter of the general noncentral chi square distribution. We concluded, from estimates of population size, linkage and measurements of population substructure, that neither drift nor population subdivision was responsible for the significant values of D which were observed and that epistatic selection was the most likely cause of the disequilibrium observed.     

Isolation and characterization of microsatellite markers in musk duck (Biziura lobata: Aves), and their application to other waterfowl species

We isolated and characterized 11 novel microsatellite loci to study paternity in the Australian musk duck (Biziura lobata), using nonradioactive PCR‐based techniques to screen GA and GAAA repeats enriched genomic DNA libraries. Nine of 11 loci showed no evidence of null alleles and were variable (mean H_E = 0.825, mean number of alleles = 9). This set of nine loci is suitable for paternity assignment (exclusion probability for nine unlinked loci = 0.9999). We also demonstrated that many of these loci cross‐amplify in various other waterfowl species.     

Effect of Mating Structure on Variation in Linkage Disequilibrium

Measurement of linkage disequilibrium involves two sampling processes. First, there is the sampling of gametes in the population to form successive generations, and this generates disequilibrium dependent on the effective population size (N_e) and the mating structure. Second, there is sampling of a finite number (n) of individuals to estimate the population disequilibrium.——Two-locus descent measures are used to describe the mating system and are transformed to disequilibrium moments at the final sampling. Approximate eigenvectors for the transition matrix of descent measures are used to obtain formulae for the variance of the observed disequilibria as a function of N_e, mating structure, n, and linkage or recombination parameter.——The variance of disequilibrium is the same for monoecious populations with or without random selfing and for dioecious populations with random pairing for each progeny. With monogamy, the variance is slightly higher, the proportional difference being greater for unlinked loci.     

Isolation and characterization of microsatellite loci from the endangered highlands scrub hypericum (Hypericum cumulicola)

We report the isolation of 19 primer pairs for amplification of polymorphic microsatellite loci for Hypericum cumulicola. These markers were evaluated in 24 individuals from one population; two to four alleles were detected per locus, and observed heterozygosity ranged from 0 to 0.5. Two loci demonstrated significant heterozygote deficiencies, possibly due to null alleles, and significant linkage disequilibrium was found between six pairs of loci. The remaining microsatellite loci will help determine if genetic differentiation is responsible for life‐history differences between natural and anthropogenically disturbed populations of H. cumulicola.     

Identification of Genetic Loci Associated with Quality Traits in Almond via Association Mapping

To design an appropriate association study, we need to understand population structure and the structure of linkage disequilibrium within and among populations as well as in different regions of the genome in an organism. In this study, we have used a total of 98 almond accessions, from five continents located and maintained at the Centro de Investigación y Tecnología Agroalimentaria de Aragón (CITA; Spain), and 40 microsatellite markers. Population structure analysis performed in ‘Structure’ grouped the accessions into two principal groups; the Mediterranean (Western-Europe) and the non-Mediterranean, with K = 3, being the best fit for our data. There was a strong subpopulation structure with linkage disequilibrium decaying with increasing genetic distance resulting in lower levels of linkage disequilibrium between more distant markers. A significant impact of population structure on linkage disequilibrium in the almond cultivar groups was observed. The mean r² value for all intra-chromosomal loci pairs was 0.040, whereas, the r² for the inter-chromosomal loci pairs was 0.036. For analysis of association between the markers and phenotypic traits, five models comprising both general linear models and mixed linear models were selected to test the marker trait associations. The mixed linear model (MLM) approach using co-ancestry values from population structure and kinship estimates (K model) as covariates identified a maximum of 16 significant associations for chemical traits and 12 for physical traits. This study reports for the first time the use of association mapping for determining marker-locus trait associations in a world-wide almond germplasm collection. It is likely that association mapping will have the most immediate and largest impact on the tier of crops such as almond with the greatest economic value.     

Microsatellite loci isolated from the tropical tree Hymenaea courbaril L. (Fabaceae)

We developed 11 microsatellite markers for Hymenaea courbaril for the purpose of studying spatial genetic structure and gene flow. The microsatellite loci were screened in 44 trees from two populations. All loci were polymorphic, exhibiting between two and 16 alleles, and levels of expected heterozygosity from 0.174 to 0.909. Departures from Hardy-Weinberg equilibrium were detected for all loci in one population. The estimated null allele frequency is low or moderate. No locus combinations exhibited linkage disequilibrium.     

Polymorphic microsatellite loci for paternity analysis in the Madagascar paradise flycatcher (Terpsiphone mutata: Aves)

Nine polymorphic microsatellite loci from the Madagascar paradise flycatcher Terpsiphone mutata were isolated using nonradioactive polymerase chain reaction (PCR)‐based techniques to screen an enriched genomic library. Seven polymorphic loci showed no evidence of null alleles and exhibited high levels of variation in 18 unrelated individuals (mean diversity = 0.80, mean number of alleles = 13.6). These loci are therefore suitable for analysis of population structure and paternity (exclusion probability for six unlinked loci = 0.9998).     

Isolation and characterization of polymorphic microsatellite loci in Hemibarbus labeo (Cyprinidae) using PCR‐based isolation of microsatellite arrays (PIMA)

A total of 10 polymorphic microsatellite loci from Hemibarbus labeo were isolated and characterized using an optimized protocol to construct a microsatellite‐enriched genomic library. The analysis of variability was performed in 24 specimens of mainland China. The mean number of alleles across loci was 3.10 ± 1.10 and the level of expected heterozygosity varied from 0.0417 to 0.7482. Frequencies of null alleles of the 10 loci are not significantly greater than zero. No linkage disequilibrium was detected between loci in either population. Five primer pairs cross‐amplify the microsatellites in other species, indicating transportability of the markers within the family Cyprinidae.     

Bias of allele-sharing linkage statistics in the presence of intermarker linkage disequilibrium

Current genome-wide linkage-mapping single-nucleotide polymorphism (SNP) panels with densities of 0.3 cM are likely to have increased intermarker linkage disequilibrium (LD) compared to 5-cM microsatellite panels. The resulting difference in haplotype frequencies versus that predicted may affect multipoint linkage analysis with ungenotyped founders; a common haplotype may be assumed to be rare, leading to inflation of identical-by-descent (IBD) allele-sharing estimates and evidence for linkage. Using data simulated for the Genetic Analysis Workshop 14, we assessed bias in allele-sharing measures and nonparametric linkage (NPL all) and Kong and Cox LOD (KC-LOD) scores in a targeted analysis of regions with and without LD and with and without genes. Using over 100 replicates, we found that if founders were not genotyped, multipoint IBD estimates and delta parameters were modestly inflated and NPL all and KC-LOD scores were biased upwards in the region with LD and no gene; rather than centering on the null, the mean NPL all and KC-LOD scores were 0.51 +/- 0.91 and 0.19 +/- 0.38, respectively. Reduction of LD by dropping markers reduced this upward bias. These trends were not seen in the non-LD region with no gene. In regions with genes (with and without LD), a slight loss in power with dropping markers was suggested. These results indicate that LD should be considered in dense scans; removal of markers in LD may reduce false-positive results although information may also be lost. Methods to address LD in a high-throughput manner are needed for efficient, robust genomic scans with dense SNPs.     

Pleiotropy of FRIGIDA enhances the potential for multivariate adaptation

An evolutionary response to selection requires genetic variation; however, even if it exists, then the genetic details of the variation can constrain adaptation. In the simplest case, unlinked loci and uncorrelated phenotypes respond directly to multivariate selection and permit unrestricted paths to adaptive peaks. By contrast, ‘antagonistic’ pleiotropic loci may constrain adaptation by affecting variation of many traits and limiting the direction of trait correlations to vectors that are not favoured by selection. However, certain pleiotropic configurations may improve the conditions for adaptive evolution. Here, we present evidence that the Arabidopsis thaliana gene FRI (FRIGIDA) exhibits ‘adaptive’ pleiotropy, producing trait correlations along an axis that results in two adaptive strategies. Derived, low expression FRI alleles confer a ‘drought escape’ strategy owing to fast growth, low water use efficiency and early flowering. By contrast, a dehydration avoidance strategy is conferred by the ancestral phenotype of late flowering, slow growth and efficient water use during photosynthesis. The dehydration avoidant phenotype was recovered when genotypes with null FRI alleles were transformed with functional alleles. Our findings indicate that the well-documented effects of FRI on phenology result from differences in physiology, not only a simple developmental switch.     

Test for interaction between two unlinked loci

Despite the growing consensus on the importance of testing gene-gene interactions in genetic studies of complex diseases, the effect of gene-gene interactions has often been defined as a deviance from genetic additive effects, which is essentially treated as a residual term in genetic analysis and leads to low power in detecting the presence of interacting effects. To what extent the definition of gene-gene interaction at population level reflects the genes' biochemical or physiological interaction remains a mystery. In this article, we introduce a novel definition and a new measure of gene-gene interaction between two unlinked loci (or genes). We developed a general theory for studying linkage disequilibrium (LD) patterns in disease population under two-locus disease models. The properties of using the LD measure in a disease population as a function of the measure of gene-gene interaction between two unlinked loci were also investigated. We examined how interaction between two loci creates LD in a disease population and showed that the mathematical formulation of the new definition for gene-gene interaction between two loci was similar to that of the LD between two loci. This finding motived us to develop an LD-based statistic to detect gene-gene interaction between two unlinked loci. The null distribution and type I error rates of the LD-based statistic for testing gene-gene interaction were validated using extensive simulation studies. We found that the new test statistic was more powerful than the traditional logistic regression under three two-locus disease models and demonstrated that the power of the test statistic depends on the measure of gene-gene interaction. We also investigated the impact of using tagging SNPs for testing interaction on the power to detect interaction between two unlinked loci. Finally, to evaluate the performance of our new method, we applied the LD-based statistic to two published data sets. Our results showed that the P values of the LD-based statistic were smaller than those obtained by other approaches, including logistic regression models.     

Mendel's First Law: partisan interests and the parliament of genes

Mendel’s First Law requires explanation because of the possibility of ‘meiotic drivers’, genes that distort fair segregation for selfish gain. The suppression of drive, and the restoration of fair segregation, is often attributed to genes at loci unlinked to the drive locus—such genes cannot benefit from drive but do suffer its associated fitness costs. However, selection can also favour suppressors at loci linked to the drive locus, raising the question of whether suppression of drive usually comes from linked or unlinked loci. Here, I study linked and unlinked suppression in a two-locus model with initial stable polymorphism at the drive locus. I find that the invasion rate of suppressors is a decreasing function of the recombination fraction between the drive and suppressor loci. Surprisingly, the relative likelihood of unlinked vs. linked suppression increases with the strength of drive and is insensitive to the fitness costs of the driver allele. I find that the chromosomal position of the driver influences how rapidly it is suppressed, with a driver in the middle of a chromosome suppressed more rapidly than a driver near the tip. When drive is strong, only a small number of chromosomes are required for suppression usually to derive from unlinked loci. In contrast, when drive is weak, and especially when suppressor alleles are associated with fitness costs, suppression will usually come from linked loci unless the genome comprises many chromosomes.Subject terms: Evolutionary genetics, Population genetics