Two-locus identity probabilities and identity disequilibrium in a partially selfing subdivided population

Measures of association of genes at different loci (linkage disequilibrium) are widely used to determine whether the structure of natural populations is clonal or not, to map genes from population data, or to test for the homogeneity of response of molecular markers to background selection, for example. However, the usual definitions of parameters for gametic associations may not be suitable for all these purposes. In this paper, we derive the recursion equations for one- and two-locus identity probabilities in an infinite island model. We study the role of drift, gene flow, partial selfing and mutation model on the expected association of genes across loci. We define the 'within-subpopulation identity disequilibrium' as the difference between the joint two-locus probability of identity in state and the expected product of one-locus identity probabilities. We evaluate this parameter as a function of recombination rate, effective size, gene flow and selfing rate. Within-subpopulation identity disequilibrium attains maximum values for intermediate immigration rates, whatever the selfing rate. Moreover, identity disequilibrium may be very small, even for high selfing rates. We discuss the implications of these findings for the analysis of data from natural populations.     

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 probabilities of identity states in permutable populations.

Génin and Clerget-Darpoux recently discussed the derivation of the probabilities of identity states for populations in which there was some degree of kinship, primarily to allow the extension of the classical affected-sib-pair method to such populations. It is argued here that their derivation makes certain assumptions that are valid only for some very restricted population models and that are not needed for an appropriate treatment. Here the probabilities of the identity states of two individuals with a given genealogical relationship are specified in terms of the kinship parameters of the underlying population, from which the founders of the individuals'' genealogy have been randomly selected. It is argued that an appropriate representation for a permutable population, one in which gene identity does not depend on the pattern of genes across individuals, requires three parameters. This representation is related to that of Génin and Clerget-Darpoux and to that of Weir.     

The expected probabilities of identity for the circular stepping-stone model

The expected values of the probabilities of identity by descent are derived for the circular stepping-stone model. The results are more easily interpreted than those derived previously.     

Simulation of selected genealogies

Algorithms for generating genealogies with selection conditional on the sample configuration of n genes in one-locus, two-allele haploid and diploid models are presented. Enhanced integro-recursions using the ancestral selection graph, introduced by S. M. Krone and C. Neuhauser (1997, Theor. Popul. Biol. 51, 210-237), which is the non-neutral analogue of the coalescent, enables accessible simulation of the embedded genealogy. A Monte Carlo simulation scheme based on that of R. C. Griffiths and S. Tavaré (1996, Math. Comput. Modelling 23, 141-158), is adopted to consider the estimation of ancestral times under selection. Simulations show that selection alters the expected depth of the conditional ancestral trees, depending on a mutation-selection balance. As a consequence, branch lengths are shown to be an ineffective criterion for detecting the presence of selection. Several examples are given which quantify the effects of selection on the conditional expected time to the most recent common ancestor.     

Gene extinction and allelic origins in complex genealogies

With the increasing emphasis on data analysis in mathematical genetics, problems of parametrizing genealogical structure become of practical importance. A complete specification of the genetic effects of genealogical structure is provided by the probabilities of genetically distinct states of gene identity by descent. Although this provides a direct parametrization for the joint distribution of traits on a set of related individuals, it is an unwieldy tool in the analysis of large and complex genealogies. Probabilities of joint descent of founder genes and likely ancestries of alleles provide alternative characterizations of relationship and have direct application in practical problems. Joint extinction probabilities of founder genes can also be derived as ancestral likelihoods: evolutionarily, the most significant characteristic of a genealogical structure must be its effect on the survival and extinction of genes.     

Intraspecific gene genealogies: trees grafting into networks

Intraspecific gene evolution cannot always be represented by a bifurcating tree. Rather, population genealogies are often multifurcated, descendant genes coexist with persistent ancestors and recombination events produce reticulate relationships. Whereas traditional phylogenetic methods assume bifurcating trees, several networking approaches have recently been developed to estimate intraspecific genealogies that take into account these population-level phenomena.     

Fixation probabilities for the Moran process with three or more strategies: general and coupling results

We study fixation probabilities for the Moran stochastic process for the evolution of a population with three or more types of individuals and frequency-dependent fitnesses. Contrary to the case of populations with two types of individuals, in which fixation probabilities may be calculated by an exact formula, here we must solve a large system of linear equations. We first show that this system always has a unique solution. Other results are upper and lower bounds for the fixation probabilities obtained by coupling the Moran process with three strategies with birth–death processes with only two strategies. We also apply our bounds to the problem of evolution of cooperation in a population with three types of individuals already studied in a deterministic setting by Núñez Rodríguez and Neves (J Math Biol 73:1665–1690, 2016). We argue that cooperators will be fixated in the population with probability arbitrarily close to 1 for a large region of initial conditions and large enough population sizes.

     

Gene genealogies in geographically structured populations.

Population genetics theory has dealt only with the spatial or geographic pattern of degrees of relatedness or genetic similarity separately for each point in time. However, a frequent goal of experimental studies is to infer migration patterns that occurred in the past or over extended periods of time. To fully understand how a present geographic pattern of genetic variation reflects one in the past, it is necessary to build genealogy models that directly relate the two. For the first time, space-time probabilities of identity by descent and coalescence probabilities are formulated and characterized in this article. Formulations for general migration processes are developed and applied to specific types of systems. The results can be used to determine the level of certainty that genes found in present populations are descended from ancient genes in the same population or nearby populations vs. geographically distant populations. Some parameter combinations result in past populations that are quite distant geographically being essentially as likely to contain ancestors of genes at a given population as the past population located at the same place. This has implications for the geographic point of origin of ancestral, "Eve," genes. The results also form the first model for emerging "space-time" molecular genetic data.     

Human hair genealogies and stem cell latency

Background  

Stem cells divide to reproduce themselves and produce differentiated progeny. A fundamental problem in human biology has been the inability to measure how often stem cells divide. Although it is impossible to observe every division directly, one method for counting divisions is to count replication errors; the greater the number of divisions, the greater the numbers of errors. Stem cells with more divisions should produce progeny with more replication errors.     

On the distribution of the general irreversiblen-compartmental model having time-dependent transition probabilities

The cumulant generating function and first two moments are derived for the stochastic distribution of units in a general irreversiblen-compartment model with time-dependent transition probabilities. In this model, a unit in the first compartment can transfer to any one of the remainingn−1 compartments and a unit in the second compartment can transfer to any of the remainingn−2 compartments and so on. In addition, a unit can enter or leave the system through any compartment. The work is related to previous research and a numerical example is given.     

A general method for biological inference: illustrated by the estimation of gene nucleotide transition probabilities

The method of maximum entropy inference developed by Jaynes can be a particularly useful method for obtaining unbiased estimates of biological parameters when the experimental knowledge about a system can be explicitly formulated. Base transition probabilities between genes, though central to evolutionary theory and understanding, present a difficult estimation problem because the ancestral genes are not experimentally accessible. The necessary estimates must therefore be made on the basis of experimental knowledge other than a direct frequency count of base replacements (A leads to C, for example) between contemporary genes. It is shown how maximum entropy inference together with the experimentally observed fact of compositional fidelity in a given gene family can be used to obtain meaningful gene base transition probabilities at each of the three nucleotide positions within codons. Both symmetric and asymmetric transition probabilities are considered. Tables of these probabilities are given for each codon position for the alpha-hemoglobin, beta-hemoglobin, myoglobin, cytochrome c, and the parvalbumin group genes. Tabular values of the average amino acid composition of these five protein families and the average nucleotide composition of their coding genes at varied codon loci are given. It is thus no longer necessary to assume in theories of evolutionary divergence equimolar base ratios A:C:G:U::1:1:1:1 or that each base has an equal chance of mutating to and being fixed as any one of the other three bases.     

Jenti: an efficient tool for mining complex inbred genealogies

SUMMARY: An efficient tool for mining complex inbred genealogies that identify clusters of individuals sharing the same expected amount of relatedness is described. Additionally it allows for the reconstruction of sub-pedigrees suitable for genetic mapping in a systematic way. AVAILABILITY: http://www.jenti.org.     

Estimating genealogies from unlinked marker data: a Bayesian approach

An issue often encountered in statistical genetics is whether, or to what extent, it is possible to estimate the degree to which individuals sampled from a background population are related to each other, on the basis of the available genotype data and some information on the demography of the population. In this article, we consider this question using explicit modelling of the pedigrees and gene flows at unlinked marker loci, but then restricting ourselves to a relatively recent history of the population, that is, considering the genealogy at most some tens of generations backwards in time. As a computational tool we use a Markov chain Monte Carlo numerical integration on the state space of genealogies of the sampled individuals. As illustrations of the method, we consider the question of relatedness at the level of genes/genomes (IBD estimation), using both simulated and real data.     

Estimating genealogies from linked marker data: a Bayesian approach

Background  

Answers to several fundamental questions in statistical genetics would ideally require knowledge of the ancestral pedigree and of the gene flow therein. A few examples of such questions are haplotype estimation, relatedness and relationship estimation, gene mapping by combining pedigree and linkage disequilibrium information, and estimation of population structure.     

Genetic markers, genealogies and biogeographic patterns in the cladocera

Cladoceran crustaceans are an important component of zooplankton in a wide range of freshwater habitats. Although the ecological characteristics of several cladoceran species have been well studied, biogeographical studies have been hampered by problematic taxonomic affiliations. However, recently developed molecular techniques, provide a powerful tool to subject aquatic taxa to comparative analyses. Here we highlight recent molecular approaches in aquatic ecology by presenting a simple method of DNA preparation and PCR amplification of the mitochondrial DNA (16S rDNA) in species from nine different families within the cladocera. On a broad taxonomic scale, sequence analysis of this mtDNA fragment has been used to produce the first molecular based phylogeny of the cladocera. This analysis clustered the cladoceran families in a fashion similar to that suggested by previous systematic classifications. In a more detailed analysis of the family Daphniidae, nuclear randomly amplified polymorphic DNA (RAPD), mitochondrial restriction fragment length polymorphism (RFLP) and morphological analyses were combined to identify species and interspecific hybrids within the Daphnia galeata species complex across 50 lakes in 13 European countries and one lake in Africa. The study revealed interspecific hybridization and backcrossing between some taxa (D. cucullata and D. galeata) to be widespread, and species and hybrids to frequently occur in sympatry. Genetic, as well as morphological information, suggests the occurrence of D. hyalina outside the Holarctic.     

Computation of identity by descent probabilities conditional on DNA markers via a Monte Carlo Markov Chain method

The accurate estimation of the probability of identity by descent (IBD) at loci or genome positions of interest is paramount to the genetic study of quantitative and disease resistance traits. We present a Monte Carlo Markov Chain method to compute IBD probabilities between individuals conditional on DNA markers and on pedigree information. The IBDs can be obtained in a completely general pedigree at any genome position of interest, and all marker and pedigree information available is used. The method can be split into two steps at each iteration. First, phases are sampled using current genotypic configurations of relatives and second, crossover events are simulated conditional on phases. Internal track is kept of all founder origins and crossovers such that the IBD probabilities averaged over replicates are rapidly obtained. We illustrate the method with some examples. First, we show that all pedigree information should be used to obtain line origin probabilities in F2 crosses. Second, the distribution of genetic relationships between half and full sibs is analysed in both simulated data and in real data from an F2 cross in pigs.