Heterozygosity-fitness correlations in zebra finches: microsatellite markers can be better than their reputation

Numerous studies have reported associations between heterozygosity in microsatellite markers and fitness-related traits (heterozygosity-fitness correlations, HFCs). However, it has often been questioned whether HFCs reflect general inbreeding depression, because a small panel of microsatellite markers does not reflect very well an individual's inbreeding coefficient (F) as calculated from a pedigree. Here, we challenge this prevailing view. Because of chance events during Mendelian segregation, an individual's realized proportion of the genome that is identical by descent (IBD) may substantially deviate from the pedigree-based expectation (i.e. F). This Mendelian noise may result in a weak correlation between F and multi-locus heterozygosity, but this does not imply that multi-locus heterozygosity is a bad estimator of realized IBD. We examined correlations between 11 fitness-related traits measured in up to 1192 captive zebra finches and three measures of inbreeding: (i) heterozygosity across 11 microsatellite markers, (ii) heterozygosity across 1359 single-nucleotide polymorphism (SNP) markers and (iii) F, based on a 5th-generation pedigree. All 11 phenotypic traits showed positive relationships with measures of heterozygosity, especially traits that are most closely related to fitness. Remarkably, the small panel of microsatellite markers produced equally strong HFCs as the large panel of SNP markers. Both marker-based approaches produced stronger correlations with phenotypes than the pedigree-based F, and this did not seem to result from the shortness of our pedigree. We argue that a small panel of microsatellites with high allelic richness may better reflect an individual's realized IBD than previously appreciated, especially in species like the zebra finch, where much of the genome is inherited in large blocks that rarely experience cross-over during meiosis.     

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.     

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.     

What should we weigh to estimate heterozygosity,alleles or loci?

The interest to study the effects of inbreeding in natural populations has increased in the last years. Several microsatellite-derived metrics have recently been developed to infer inbreeding from multilocus heterozygosity data without requiring detailed pedigrees that are difficult to obtain in open populations. Internal relatedness (IR) is currently the most widespread used index and its main attribute is that allele frequency is incorporated into the measure. However, IR underestimates heterozygosity of individuals carrying rare alleles. For example, descendants of immigrants paired with natives (normally more outbred) bearing novel or rare alleles would be considered more homozygous than descendants of native parents. Thus, the analogy between homozygosity and inbreeding that generally is carried out would have no logic in those cases. We propose an alternative index, homozygosity by loci (HL) that avoids such problems by weighing the contribution of each locus to the homozygosity index depending on their allelic variability. Under a wide range of simulated scenarios, we found that our index (HL) correlated better than both IR and uncorrected homozygosity (H(O)), measured as proportion of homozygous loci) with genome-wide homozygosity and inbreeding coefficients in open populations. In these populations, which are likely to prevail in nature, the use of HL instead of IR reduced considerably the sample sizes required to achieve a given statistical power. This is likely to have important consequences on the ability to detect heterozygosity fitness correlations assuming the relationship between genome-wide heterozygosity and fitness traits.