Prediction of identity by descent probabilities from marker-haplotypes

The prediction of identity by descent (IBD) probabilities is essential for all methods that map quantitative trait loci (QTL). The IBD probabilities may be predicted from marker genotypes and/or pedigree information. Here, a method is presented that predicts IBD probabilities at a given chromosomal location given data on a haplotype of markers spanning that position. The method is based on a simplification of the coalescence process, and assumes that the number of generations since the base population and effective population size is known, although effective size may be estimated from the data. The probability that two gametes are IBD at a particular locus increases as the number of markers surrounding the locus with identical alleles increases. This effect is more pronounced when effective population size is high. Hence as effective population size increases, the IBD probabilities become more sensitive to the marker data which should favour finer scale mapping of the QTL. The IBD probability prediction method was developed for the situation where the pedigree of the animals was unknown (i.e. all information came from the marker genotypes), and the situation where, say T, generations of unknown pedigree are followed by some generations where pedigree and marker genotypes are known.     

On the genealogy of a duplicated microsatellite

When a microsatellite locus is duplicated in a diploid organism, a single pair of PCR primers may amplify as many as four distinct alleles. To study the evolution of a duplicated microsatellite, we consider a coalescent model with symmetric stepwise mutation. Conditional on the time of duplication and a mutation rate, both in a model of completely unlinked loci and in a model of completely linked loci, we compute the probabilities for a sampled diploid individual to amplify one, two, three, or four distinct alleles with one pair of microsatellite PCR primers. These probabilities are then studied to examine the nature of their dependence on the duplication time and the mutation rate. The mutation rate is observed to have a stronger effect than the duplication time on the four probabilities, and the unlinked and linked cases are seen to behave similarly. Our results can be useful for helping to interpret genetic variation at microsatellite loci in species with a very recent history of gene and genome duplication.     

The Effects of Overdominance on Linkage in a Multilocus System

Computer simulations were performed with overdominant multiple alleles among tightly linked multiple loci under a multiplicative fitness model. The quantity X²/N(n — 1) was introduced as a new measure of linkage disequilibrium which, unlike previously available measures, can be applied to multiple allele models, where N is the sample size, and n is the number of alleles at the locus possessing fewest alleles. Simulations showed that (1) With multiple (three or four) alleles, the approach to stable disequilibrium is slower and the amount of disequilibrium established is weaker than in a two allele system. (2) The number of complementary chromosomes is a function of number of alleles and of population size. (3) As population size increases, the rate of the approach to stable disequilibrium is slower. (4) There is an optimum selection coefficient which minimizes the transient fixation probability of alleles when linkage is present. (5) The absence of linkage disequilibrium is in most cases not a practical method of testing the hypothesis of balancing selection of genetic polymorphisms because it depends strongly on population size in determining linkage disequilibria.     

The inbreeding effective number and the effective number of alleles in a population that varies in size

Consider a population that does not change in size. If it is assumed that there are an infinite number of possible neutral alleles at a locus and u is the probability that a particular gene mutates to some other gene in one generation, the effective number of alleles n_e is computed to be 4N_eu + 1, where N_e is the inbreeding effective population number. It is assumed in this paper that the number of individuals in a monoecious population, or the numbers of males and females in a dioecious population, are states in a finite irreducible Markov chain. In general it is impossible to obtain a single value of n_e. In some cases where the computation of n_e is possible, the results are as follows. When the population is monoecious, N_e is the reciprocal of the asymptotic average, over population sizes, of the probabilities that two gametes uniting to form an individual came from the same individual one generation earlier. In dioecious populations, N_e is the reciprocal of the long-run average of the probabilities that two homologous genes in separate individuals of one generation came from the same individual one generation earlier. Special cases are discussed.     

Genic VERSUS Chromosomal Variation in Natural Populations of DROSOPHILA SUBOBSCURA

Gametic frequencies in one mainland and one island population of D. subobscura were obtained by means of extracting wild chromosomes and subsequently analyzing them for inversions and allozymes. The high degree of cytological heterogeneity which characterizes these populations is not reflected in the genetic data. Two cases of non-random association were observed among eighteen pair-wise comparisons involving gene alleles and inversions to which the locus is linked. In both cases exchange of alleles at the locus is completely suppressed by the inversions. Four cases of linkage disequilibrium were detected among eighteen pairs of loci; two of them could best be explained as transient associations generated by random drift. The results suggest that disequilibria among enzyme loci are not widespread in natural populations—Populations with a lower degree of chromosomal variation are genetically as variable as populations with a higher degree of chromosomal variation. This observation does not support the hypothesis that selection in marginal homokaryotypic populations is for specialized homozygous genotypes.     

Development and characterization of 11 microsatellite markers in the rock sparrow,Petronia petronia

Eleven polymorphic microsatellite loci were isolated in the passeriform Petronia petronia using the enrichment protocol FIASCO (fast isolation by AFLP of sequences containing repeats). We detected three to 13 alleles per locus in 25 specimens collected from an Italian population. The level of expected heterozygosity ranged from 0.439 to 0.856. One locus is sex linked to the Z chromosome. The total exclusionary probabilities using these loci for the first and the second parent were 0.978 and 0.999, respectively. These are the first microsatellite loci characterized from the rock sparrow that can be used for estimating population structure and for large‐scale parentage analysis.     

元江鲤种群遗传多样性

选择12对微卫星标记检测了于2011年采集自元江(红河上游中国江段)5个样点192尾鲤的群体遗传多样性.共检测到201个等位基因,每个位点等位基因2-27个.各群体各位点平均等位基因(NA)12.25-14.67个,平均有效等位基因(NE)8.28-9.73个,平均观察杂合度(Ho)o.7765-0.8037,平均期望杂合度(HE)0.7761-0.8080,平均多态信息含量(PIC)0.7534-0.7843.元江鲤种群192个个体各位点NA、NE、Ho、HE、PIC分别为16.50、11.26、0.7927、0.8049、0.7966,种群遗传多样性水平高.元江鲤群体之间遗传分化小,可作为一个种群管理单元进行管理.增殖放流要防止遗传多样性丧失.     

The Evolution of the Y Chromosome with X-Y Recombination

A theoretical population genetic model is developed to explore the consequences of X-Y recombination in the evolution of sex chromosome polymorphism. The model incorporates one sex-determining locus and one locus subject to natural selection. Both loci have two alleles, and the rate of classical meiotic recombination between the loci is r. The alleles at the sex-determining locus specify whether the chromosome is X or Y, and the alleles at the selected locus are arbitrarily labeled A and a. Natural selection is modeled as a process of differential viabilities. The system can be expressed in terms of three recurrence equations, one for the frequency of A on the X-bearing gametes produced by females, one for each of the frequency of A on the X- and Y-bearing gametes produced by males. Several special cases are examined, including X chromosome dominance and symmetric selection. Unusual equilibria are found with the two sexes having very different allele frequencies at the selected locus. A significant finding is that the allowance of recombination results in a much greater opportunity for polymorphism of the Y chromosome. Tighter linkage results in a greater likelihood for equilibria with a large difference between the sex chromosomes in allele frequency.     

Derivation of single-locus relationship coefficients conditional on marker information

The coefficient of relationship is defined as the correlation between the additive genetic values of two individuals. This coefficient can be defined specifically for a single quantitative trait locus (QTL) and may deviate considerably from the overall expectation if it is taken conditional on information from linked marker loci. Conditional halfsib correlations are derived under a simple genetic model with a biallelic QTL linked to a biallelic marker locus. The conditional relationship coefficients are shown to depend on the recombination rate between the marker and the QTL and the population frequency of the marker alleles, but not on parameters of the QTL, i.e. number and frequency of QTL alleles, degree of dominance etc., nor on the (usually unknown) QTL genotype of the sire. Extensions to less simplified cases (multiple alleles at the marker locus and the QTL, two marker loci flanking the QTL) are given. For arbitrary pedigrees, conditional relationship coefficients can also be derived from the conditional gametic covariance matrix suggested by Fernando and Grossman (1989). The connection of these two approaches is discussed. The conditional relationship coefficient can be used for marker-assisted genetic evaluation as well as for the detection of QTL and the estimation of their effects.     

Hitchhiking and associative overdominance at a microsatellite locus

The possible effects of a selected locus on a closely linked microsatellite locus are discussed and analyzed in terms of coalescent theory and models of the mutation process. Background selection caused by recurrent deleterious mutations will reduce the variance of allele size at a microsatellite locus. The occasional substitution of advantageous alleles (genetic hitchhiking) will also reduce the variance, but a high mutation rate at a microsatellite locus can restore the variance relatively rapidly. Overdominance at the selected locus will increase the variance at the microsatellite locus and create partitioning of the variation in allele size among gametes carrying one or the other of the overdominant alleles. These results suggest that neutral microsatellite loci can provide indicators of selective processes at closely linked loci.      

Identifying Signatures of Selection in Genetic Time Series

Both genetic drift and natural selection cause the frequencies of alleles in a population to vary over time. Discriminating between these two evolutionary forces, based on a time series of samples from a population, remains an outstanding problem with increasing relevance to modern data sets. Even in the idealized situation when the sampled locus is independent of all other loci, this problem is difficult to solve, especially when the size of the population from which the samples are drawn is unknown. A standard χ²-based likelihood-ratio test was previously proposed to address this problem. Here we show that the χ²-test of selection substantially underestimates the probability of type I error, leading to more false positives than indicated by its P-value, especially at stringent P-values. We introduce two methods to correct this bias. The empirical likelihood-ratio test (ELRT) rejects neutrality when the likelihood-ratio statistic falls in the tail of the empirical distribution obtained under the most likely neutral population size. The frequency increment test (FIT) rejects neutrality if the distribution of normalized allele-frequency increments exhibits a mean that deviates significantly from zero. We characterize the statistical power of these two tests for selection, and we apply them to three experimental data sets. We demonstrate that both ELRT and FIT have power to detect selection in practical parameter regimes, such as those encountered in microbial evolution experiments. Our analysis applies to a single diallelic locus, assumed independent of all other loci, which is most relevant to full-genome selection scans in sexual organisms, and also to evolution experiments in asexual organisms as long as clonal interference is weak. Different techniques will be required to detect selection in time series of cosegregating linked loci.     

The number of heterozygous loci between two randomly chosen completely linked sequences of loci in two subdivided population models

The stationary probability distribution of the number of heterozygous loci in two randomly chosen sequences of completely linked infinite alleles loci, with mutation at each locus, is found in the island model for within and between islands. Results for an infinite site model are found as a limit. A single charge state locus is also studied in the island model and distributions found for the charge difference between two genes. Similar results are derived for a stepping stone model.     

On the Theory of Partially Inbreeding Finite Populations. I. Partial Selfing

Some stochastic theory is developed for monoecious populations of size N in which there are probabilities beta and 1 - beta of reproduction by selfing and by random mating. It is assumed that beta much greater than N-1. Expressions are derived for the inbreeding coefficient of one random individual and the coefficient of kinship of two random separate individuals at time t. The mean and between-lines variance of the fraction of copies of a locus that are identical in two random separate individuals in an equilibrium population are obtained under the assumption that there is an infinite number of possible alleles. It is found that the theory for random mating populations holds if the effective population number is Ne = N'/(1 + FIS), where FIS is the inbreeding coefficient at equilibrium when N is infinite and N' is the reciprocal of the probability that two gametes contributing to random separate adults come from the same parent. When there is a binomial distribution of successful gametes emanating from each adult, N' = N. An approximation to the probability that an allele A survives if it is originally present in one AA heterozygote is found to be 2(N'/N)(FISS1 + (1 - FIS)S2), where S1 and S2 are the selective advantages of AA and AA in comparison with AA. In the last section it is shown that if there is partial full sib mating and binomial offspring distributions Ne = N/(1 + 3FIS).     

The Coalescent Process in Models with Selection and Recombination

The statistical properties of the process describing the genealogical history of a random sample of genes at a selectively neutral locus which is linked to a locus at which natural selection operates are investigated. It is found that the equations describing this process are simple modifications of the equations describing the process assuming that the two loci are completely linked. Thus, the statistical properties of the genealogical process for a random sample at a neutral locus linked to a locus with selection follow from the results obtained for the selected locus. Sequence data from the alcohol dehydrogenase (Adh) region of Drosophila melanogaster are examined and compared to predictions based on the theory. It is found that the spatial distribution of nucleotide differences between Fast and Slow alleles of Adh is very similar to the spatial distribution predicted if balancing selection operates to maintain the allozyme variation at the Adh locus. The spatial distribution of nucleotide differences between different Slow alleles of Adh do not match the predictions of this simple model very well.     

De novo developed microsatellite markers in gill parasites of the genus Dactylogyrus (Monogenea): Revealing the phylogeographic pattern of population structure in the generalist parasite Dactylogyrus vistulae

Approaches using microsatellite markers are considered the gold standard for modern population genetic studies. However, although they have found application in research into various platyhelminth taxa, they remained substantially underutilized in the study of monogeneans. In the present study, a newly developed set of 24 microsatellite markers was used to investigate the genetic diversity of the generalist monogenean species Dactylogyrus vistulae. The analyzed parasite specimens were collected from 13 cyprinoid species from 11 sites in the Apennine and Balkan peninsulas. A total of 159 specimens were genotyped at each of the loci and the number of alleles per locus ranged from 2 to 16, with a mean number of 6.958 alleles per locus. Exceptionally high genetic diversity was observed among D. vistulae individuals in the southern Balkans (mean N _A per locus = 3.917), suggesting that generalist D. vistulae expanded from the south to the north in the Balkans and later into central Europe. The initial clustering analysis divided all investigated specimens into three major clusters; however, the results of the subsequent analyses revealed the existence of various subpopulations, suggesting that the population structure of D. vistulae is associated with the diversification of their cyprinoid hosts. In addition, the partition of the parasite population was observed in regions of the sympatric occurrence of two host species, indicating that these hosts may represent a barrier for gene flow, even for generalist parasite species.     

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.     

Characterization of microsatellite loci isolated in Mountain Plover (Charadrius montanus)

Primers for 15 microsatellite loci were developed for Mountain Plover, a species whose distribution and abundance have been reduced drastically in the past 30 years. In a screen of 126 individuals collected from four breeding locales across the species’ range, levels of polymorphism ranged from two to 13 alleles per locus. No two loci were found to be linked, although one locus revealed significant departures from Hardy–Weinberg equilibrium. These microsatellite loci can be used in population genetic studies, ultimately aiding in management efforts for Mountain Plover. Additionally, these markers can potentially be used in studies investigating the mating system of Mountain Plover.     

Genetic variability in European <Emphasis Type="Italic">Leptestheria dahalacensis</Emphasis> (Rüppel, 1837) (Crustacea,Branchiopoda, Spinicaudata)

The genetic variability of the gonochoric Leptestheria dahalacensis (Rüppel, 1837) was studied through the analysis of mitochondrial and nuclear (microsatellite loci) markers in eight Italian and two Central European populations. Mitochondrial data exhibited a low variability, as only six mitotypes were scored: five in Italy and one for both Central European samples, with a very low number of substitutions. All analysed microsatellite loci were variable, with 3–5 alleles per locus and 1–4 alleles per population. All populations were at the Hardy–Weinberg equilibrium, with the exceptions of two samples for locus ldAC-16, due to heterozygote excess, and of four populations for locus ldAC-11, probably linked to the presence of null alleles. A substantial population structuring was found between Central European and Italian samples for both utilized markers. This observation may be explained by isolation by distance and/or recent isolation events. On the other hand, the absence of a clear inter-pond variability in Italian sample comparisons may be ascribed to high dispersal ability in the short range. Handling editor: K. Martens     

Genetic Diversity of Some Pear Cultivars and Genotypes Using Simple Sequence Repeat (SSR) Markers

The genetic diversity and relationships among 47 pear cultivars and genotypes (Pyrus spp.), including 4 Japanese pears (Pyrus pyrifolia), 40 European pears (Pyrus communis), 1 Chinese pear (Pyrus bretschneideri) as well as 2 wild relatives (Pyrus salicifolia and Pyrus mazandaranica) were studied using 28 microsatellite primer pairs. A total of 174 alleles were produced at the 28 SSR loci with their sizes ranging from 81 to 290?bp. The number of observed alleles for each locus ranged from 3 (TsuENH014 and TsuENH046) to 12 (NB103a), with an average of 6.21 alleles per locus. In some SSR loci, more than two alleles were amplified in some cultivars and genotypes, suggesting that duplication has occurred in those accessions. This information suggests that at least two genomic regions exist for these loci in the pear genome. The observed heterozygosity (H _o) values of amplified loci ranged from 0.17 (TsuENH006) to 0.97 (NB103a). Shannon's information index (I) value was observed to be highest (2.14) in the NB103a locus, while the TsuENH006 locus had the lowest value with an average of 1.37 among SSR loci. The Dice genetic similarity coefficient ranged from 0.29 (??Nijisseiki?? and P. mazandaranica) to 0.91 (??Chojuro?? and ??Nijisseiki??) among samples. UPGMA cluster analysis showed two major groups corresponding to the Japanese and European pears.