A model‐based approach to characterize individual inbreeding at both global and local genomic scales

Inbreeding results from the mating of related individuals and may be associated with reduced fitness because it brings together deleterious variants in one individual. In general, inbreeding is estimated with respect to an arbitrary base population consisting of ancestors that are assumed unrelated. We herein propose a model‐based approach to estimate and characterize individual inbreeding at both global and local genomic scales by assuming the individual genome is a mosaic of homozygous‐by‐descent (HBD) and non‐HBD segments. The HBD segments may originate from ancestors tracing back to different periods in the past defining distinct age‐related classes. The lengths of the HBD segments are exponentially distributed with class‐specific parameters reflecting that inbreeding of older origin generates on average shorter stretches of observed homozygous markers. The model is implemented in a hidden Markov model framework that uses marker allele frequencies, genetic distances, genotyping error rates and the sequences of observed genotypes. Note that genotyping errors, low‐fold sequencing or genotype‐by‐sequencing data are easily accommodated under this framework. Based on simulations under the inference model, we show that the genomewide inbreeding coefficients and the parameters of the model are accurately estimated. In addition, when several inbreeding classes are simulated, the model captures them if their ages are sufficiently different. Complementary analyses, either on data sets simulated under more realistic models or on human, dog and sheep real data, illustrate the range of applications of the approach and how it can reveal recent demographic histories among populations (e.g., very recent bottlenecks or founder effects). The method also allows to clearly identify individuals resulting from extreme consanguineous matings.     

From inclusive fitness to fixation probability in homogeneous structured populations

The methods of inclusive fitness provide a powerful analysis of the action of selection on social behaviour. The key component of this analysis is the concept of relatedness R. In infinite populations, a standard method of calculating relatedness coefficients is through coefficients of consanguinity using the notion of genetic identity by descent. In this paper, we show that this approach can also be made to work in finite populations and we assume here that the population has a homogeneous structure, such as an island model. We demonstrate that, under the assumption that genetic effects are small and additive, the resulting formulation of inclusive fitness is equivalent to other significant measures of selection in finite populations, including the change in average allele frequency and fixation probability. The results are illustrated for a model of the evolution of cooperation in a finite island population.     

Native South American genetic structure and prehistory inferred from hierarchical modeling of mtDNA

Genetic diversity in Native South Americans forms a complexpattern at both the continental and local levels. In comparingthe West to the East, there is more variation within groupsand smaller genetic distances between groups. From this pattern,researchers have proposed that there is more variation in theWest and that a larger, more genetically diverse, founding populationentered the West than the East. Here, we question this characterizationof South American genetic variation and its interpretation.Our concern arises because others have inferred regional variationfrom the mean variation within local populations without takinginto account the variation among local populations within thesame region. This failure produces a biased view of the actualvariation in the East. In this study, we analyze the mitochondrial DNA sequence betweenpositions 16040 and 16322 of the Cambridge reference sequence.Our sample represents a total of 886 people from 27 indigenouspopulations from South (22), Central (3), and North America(2). The basic unit of our analyses is nucleotide identity bydescent, which is easily modeled and proportional to nucleotidediversity. We use a forward modeling strategy to fit a seriesof nested models to identity by descent within and between allpairs of local populations. This method provides estimates ofidentity by descent at different levels of population hierarchywithout assuming homogeneity within populations, regions, orcontinents. Our main discovery is that Eastern South America harbors moregenetic variation than has been recognized. We find no evidencethat there is increased identity by descent in the East relativeto the total for South America. By contrast, we discovered thatpopulations in the Western region, as a group, harbor more identityby descent than has been previously recognized, despite thefact that average identity by descent within groups is lower.In this light, there is no need to postulate separate foundingpopulations for the East and the West because the variabilityin the East could serve as a source for the Western gene pools.     

Combinatorial Conflicting Homozygosity (CCH) analysis enables the rapid identification of shared genomic regions in the presence of multiple phenocopies

Background

The ability to identify regions of the genome inherited with a dominant trait in one or more families has become increasingly valuable with the wide availability of high throughput sequencing technology. While a number of methods exist for mapping of homozygous variants segregating with recessive traits in consanguineous families, dominant conditions are conventionally analysed by linkage analysis, which requires computationally demanding haplotype reconstruction from marker genotypes and, even using advanced parallel approximation implementations, can take substantial time, particularly for large pedigrees. In addition, linkage analysis lacks sensitivity in the presence of phenocopies (individuals sharing the trait but not the genetic variant responsible). Combinatorial Conflicting Homozygosity (CCH) analysis uses high density biallelic single nucleotide polymorphism (SNP) marker genotypes to identify genetic loci within which consecutive markers are not homozygous for different alleles. This allows inference of identical by descent (IBD) inheritance of a haplotype among a set or subsets of related or unrelated individuals.

Results

A single genome-wide conflicting homozygosity analysis takes <3 seconds and parallelisation permits multiple combinations of subsets of individuals to be analysed quickly. Analysis of unrelated individuals demonstrated that in the absence of IBD inheritance, runs of no CH exceeding 4 cM are not observed. At this threshold, CCH is >97% sensitive and specific for IBD regions within a pedigree exceeding this length and was able to identify the locus responsible for a dominantly inherited kidney disease in a Turkish Cypriot family in which six out 17 affected individuals were phenocopies. It also revealed shared ancestry at the disease-linked locus among affected individuals from two different Cypriot populations.

Conclusions

CCH does not require computationally demanding haplotype reconstruction and can detect regions of shared inheritance of a haplotype among subsets of related or unrelated individuals directly from SNP genotype data. In contrast to parametric linkage allowing for phenocopies, CCH directly provides the exact number and identity of individuals sharing each locus. CCH can also identify regions of shared ancestry among ostensibly unrelated individuals who share a trait. CCH is implemented in Python and is freely available (as source code) from http://sourceforge.net/projects/cchsnp/.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1360-4) contains supplementary material, which is available to authorized users.     

Genetic variation of recent Alu insertions in human populations

The Alu family of intersperesed repeats is comprised of ovr 500,000 members which may be divided into discrete subfamilies based upon mutations held in common between members. Distinct subfamilies of Alu sequences have amplified within the human genome in recent evolutionary history. Several individual Alu family members have amplified so recently in human evolution that they are variable as to presence and absence at specific loci within different human populations. Here, we report on the distribution of six polymorphic Alu insetions in a survey of 563 individuals from 14 human population groups across several continents. Our results indicate that these polymorphic Alu insertions probably have an African origin and that there is a much smaller amount of genetic variation between European populations than that found between other populations groups. Present address: Department of Pathology, Stanley S. Scott Cancer Center, Louisiana State University Medical Center, 1901 Perdido St., New Orleans, LA 70112 Correspondence to: M.A. Batzer     

Marker‐based estimates of relatedness and inbreeding coefficients: an assessment of current methods

Inbreeding (F) of and relatedness (r) between individuals are now routinely calculated from marker data in studies in the fields of quantitative genetics, conservation genetics, forensics, evolution and ecology. Although definable in terms of either correlation coefficient or probability of identity by descent (IBD) relative to a reference, they are better interpreted as correlations in marker‐based analyses because the reference in practice is frequently the current sample or population whose F and r are being estimated. In such situations, negative estimates have a biological meaning, a substantial proportion of the estimates are expected to be negative, and the average estimates are close to zero for r and equivalent to F_IS for F. I show that although current r estimators were developed from the IBD‐based concept of relatedness, some of them conform to the correlation‐based concept of relatedness and some do not. The latter estimators can be modified, however, so that they estimate r as a correlation coefficient. I also show that F and r estimates can be misleading and become biased and marker dependent when a sample containing a high proportion of highly inbred and/or closely related individuals is used as reference. In analyses depending on the comparison between r (or F) estimates and a priori values expected under ideal conditions (e.g. for identifying genealogical relationship), the estimators should be used with caution.     

The probability density of the total IBD length over a single autosome in unilineal relationships

Several authors have studied identity by descent (IBD) by way of a continuous recombination process along a chromosome. Despite its potential uses in, for example, gene mapping or delineation of biological relationships there has been no exact algebraic result given for the probability density function of the IBD proportion in any familial relationship. Other authors have derived algebraic approximations in the case of half-sibs by way of the Poisson clumping heuristic and used computational methods to compute the distribution function of the IBD sharing for unilineal relationships. Here we provide a general numerical method for finding the density of IBD sharing that could be applied to any unilineal relationship and more importantly we derive algebraically an expression for the density for a grandparent-grandchild relationship. Initially we assume that recombination events occur at random along a chromosome, then go on to show how the method could be extended to incorporate a form of genetic interference.     

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.     

Estimating pairwise relatedness from dominant genetic markers

Knowledge of the genetic relatedness between a pair of individuals is important in many research areas of quantitative genetics, conservation genetics, evolution and ecology. Many estimators have been developed to estimate such pairwise relatedness (r) using codominant markers, such as microsatellites and enzymes. In contrast, only two estimators are proposed to use dominant markers, such as random amplified polymorphic DNAs (RAPDs) and amplified fragment length polymorphisms (AFLPs), in relatedness inference. They are both biased estimators, and their statistical properties and robustness to the sampling errors in allele frequency have not been investigated. In this short paper, I propose two new pairwise relatedness estimators for dominant markers, and compare them in precision, accuracy and robustness to sampling with the two previous estimators using simulations. It was found that the new estimator based on the least squares approach is unbiased when allele frequencies are known or estimated from a sample without correcting for sampling effects. It has, however, a low precision and as a result, an intermediate overall performance among the four estimators in terms of the mean squared deviation (MSD) of estimates from actual values of r. The new estimator based on a similarity index is slightly biased but has generally the lowest MSD among the four estimators compared, regardless of the number of loci, type of actual relationships, allele frequencies known or estimated from samples. Simulations also show that the confidence intervals estimated by bootstrapping are appropriate for different estimators provided that the number of loci used in the estimation is not small.     

The genetic history of the Southern Andes from present-day Mapuche ancestry

Download : Download video (32MB)     

Increased runs of homozygosity in the autosomal genome of Brazilian individuals with neurodevelopmental delay/intellectual disability and/or multiple congenital anomalies investigated by chromosomal microarray analysis

Runs of homozygosity (ROH) in the human genome may be clinically relevant. The aim of this study was to report the frequency of increased ROH of the autosomal genome in individuals with neurodevelopmental delay/intellectual disability and/or multiple congenital anomalies, and to compare these data with a control group. Data consisted of calls of homozygosity from 265 patients and 289 controls. In total, 7.2% (19/265) of the patients showed multiple ROH exceeding 1% of autosomal genome, compared to 1.4% (4/289) in the control group (p=0.0006). Homozygosity ranged from 1.38% to 22.12% among patients, and from 1.53 to 2.40% in the control group. In turn, 1.9% (5/265) of patients presented ROH ≥10Mb in a single chromosome, compared to 0.3% (1/289) of individuals from the control group (p=0.0801). By excluding cases with reported consanguineous parents (15/24), the frequency of increased ROH was 3.4% (9/250) among patients and 1.7% (5/289) in the control group, considering multiple ROH exceeding 1% of the autosome genome and ROH ≥10Mb in a single chromosome together, although not statistically significant (p=0.1873). These results reinforce the importance of investigating ROH, which with complementary diagnostic tests can improve the diagnostic yield for patients with such conditions.     

Sib mating designs for mapping quantitative trait loci

The power to separate the variance of a quantitative trait locus (QTL) from the polygenic variance is determined by the variability of genes identical by descent (IBD) at the QTL. This variability may increase with inbreeding. Selfing, the most extreme form of inbreeding, increases the variability of the IBD value shared by siblings, and thus has a higher efficiency for QTL mapping than random mating. In self-incompatible organisms, sib mating is the closest form of inbreeding. Similar to selfing, sib mating may also increase the power of QTL detection relative to random mating. In this study, we develop an IBD-based method under sib mating designs for QTL mapping. The efficiency of sib mating is then compared with random mating. Monte Carlo simulations show that sib mating designs notably increase the power for QTL detection. When power is intermediate, the power to detect a QTL using full-sib mating is, on average, 7% higher than under random mating. In addition, the IBD-based method proposed in this paper can be used to combine data from multiple families. As a result, the estimated QTL parameters can be applied to a wide statistical inference space relating to the entire reference population. This revised version was published online in July 2006 with corrections to the Cover Date.     

Ubiquity or not ubiquity: That is the question

The nature of population structure in eukaryotic microbes has been the subject of intense debate, but until recently the tools to test these hypotheses were either problematic (e.g., allozymes that cannot detect all genetic changes) or beyond financial and technological limits of most laboratories (e.g., high throughput sequencing). In a recent issue of Molecular Ecology, Craig et al. (2019) use a genomic approach to investigate the population structure of a model alga, the chlorophyte Chlamydomonas reinhardtii (Figure 1). Using high throughput sequencing, read mapping, and variant calling, they detected strong signals of differentiation at a continental scale, while local patterns of admixture were complex. Population genomic techniques such as these have not been used extensively in studies of microbial eukaryotes and the fields of conservation genetics and evolution stand to benefit vastly from the adoption of these techniques to studies of diverse protist lineages.     

A Likelihood-Based Trait-Model-Free Approach for Linkage Detection of Binary Trait

Summary .  Trait-model-free (or "allele-sharing") approach to linkage analysis is a popular tool in genetic mapping of complex traits, because of the absence of explicit assumptions about the underlying mode of inheritance of the trait. The likelihood framework introduced by Kong and Cox (1997,  American Journal of Human Genetics   61, 1179–1188) allows calculation of accurate p-values and LOD scores to test for linkage between a genomic region and a trait. Their method relies on the specification of a model for the trait-dependent segregation of marker alleles at a genomic region linked to the trait. Here we propose a new such model that is motivated by the desire to extract as much information as possible from extended pedigrees containing data from individuals related over several generations. However, our model is also applicable to smaller pedigrees, and has some attractive features compared with existing models ( Kong and Cox, 1997 ), including the fact that it incorporates information on both affected and unaffected individuals. We illustrate the proposed model on simulated and real data, and compare its performance with the existing approach ( Kong and Cox, 1997 ). The proposed approach is implemented in the program lm_ibdtests within the framework of MORGAN 2.8 ( http://www.stat.washington.edu/thompson/Genepi/MORGAN/Morgan.shtml ).     

Haplotype kinship for three populations of the Goettingen minipig

To overcome limitations of diversity measures applied to livestock breeds marker based estimations of kinship within and between populations were proposed. This concept was extended from the single locus consideration to chromosomal segments of a given length in Morgan. Algorithms for the derivation of haplotype kinship were suggested and the behaviour of marker based haplotype kinship was investigated theoretically. In the present study the results of the first practical application of this concept are presented. Full sib pairs of three sub-populations of the Goettingen minipig were genotyped for six chromosome segments. After haplotype reconstruction the haplotypes were compared and mean haplotype kinships were estimated within and between populations. Based on haplotype kinships a distance measure is proposed which is approximatively linear with the number of generations since fission. The haplotype kinship distances, the respective standard errors and the pedigree-based expected values are presented and are shown to reflect the true population history better than distances based on single-locus kinships. However the marker estimated haplotype kinship reveals variable among segments. This leads to high standard errors of the respective distances. Possible reasons for this phenomenon are discussed and a pedigree-based approach to correct for identical haplotypes which are not identical by descent is proposed.     

相关个体基因型联合概率分布及在身份鉴定中的应用

从联合父系基因概率出发,得出处在同一代的多个个体的联合基因型概率,讨论两种符合我国国情的家谱图,得到同一家族内第m代独生子女之间的联合基因型概率,相应的方法可用来求多个家族、多代独生子女之间的联合基因型概率．列举了两个案例来说明相关个体的联合基因型概率在身份鉴定中的应用．     

Modeling the effect of inbreeding among founders in linkage analysis

In this paper, we present a unified mathematical model for linkage analysis that allows for inbreeding among founders in all families. The identical by descent (IBD) configuration of each pedigree is modeled as a Markov process containing two parameters; the inverse inbreeding and kinship coefficient and a rate parameter proportional to the inverse expected length of chromosome segments shared IBD by two different founder haplotypes. We use hidden Markov models and define a forward-backward algorithm for computing the conditional IBD-distribution given marker data, thereby extending the multipoint method of Lander and Green [1987. Construction of multilocus genetic maps in humans, Proc. Natl. Acad. Sci. USA 84, 2363-2367] to situations where founders are inbred. Our methodology is valid for arbitrary pedigree structures. Simulation and theoretical approximations for nonparametric linkage (NPL) analysis based on affected sib pairs reveal that NPL scores are inflated and type 1 errors increased when the inbreeding coefficient or rate parameter is underestimated. When the parents are genotyped, we present a general way of modifying the score function to drastically reduce this effect.