Benchmarking Relatedness Inference Methods with Genome-Wide Data from Thousands of Relatives

Inferring relatedness from genomic data is an essential component of genetic association studies, population genetics, forensics, and genealogy. While numerous methods exist for inferring relatedness, thorough evaluation of these approaches in real data has been lacking. Here, we report an assessment of 12 state-of-the-art pairwise relatedness inference methods using a data set with 2485 individuals contained in several large pedigrees that span up to six generations. We find that all methods have high accuracy (92–99%) when detecting first- and second-degree relationships, but their accuracy dwindles to \u3c43% for seventh-degree relationships. However, most identical by descent (IBD) segment-based methods inferred seventh-degree relatives correct to within one relatedness degree for \u3e76% of relative pairs. Overall, the most accurate methods are Estimation of Recent Shared Ancestry (ERSA) and approaches that compute total IBD sharing using the output from GERMLINE and Refined IBD to infer relatedness. Combining information from the most accurate methods provides little accuracy improvement, indicating that novel approaches, such as new methods that leverage relatedness signals from multiple samples, are needed to achieve a sizeable jump in performance

Appendix: Wing photomicrographs and averaged yellow reflectance spectra

Photomicrographs of butterfly wings and averaged yellow reflectance spectra. Color-coded circles indicate the location from which the associated reflectance spectra were sampled. N=number of individual specimens measured. Left column: Dorsal wings and dorsal reflectance spectra, Right column: Ventral wings and ventral reflectance spectra. FD, forewing dorsal; FV, forewing ventral; HD, hindwing dorsal; HV, hindwing ventral

Reflectance spectra of Heliconiine butterfly yellow wing pigments

Reflectance spectra of yellow wing colors of Heliconius butterflies and their relatives used in the color space model of birds and butterflies shown in Figure 4 of Bybee et al. Averaged reflectances and location of individual measurements on the wings are shown in Appendix

Data from: UV photoreceptors and UV-yellow wing pigments in Heliconius butterflies allow a color signal to serve both mimicry and intraspecific communication

Mimetic wing coloration evolves in butterflies in the context of predator confusion. Unless butterfly eyes have adaptations for discriminating mimetic color variation, mimicry also carries a risk of confusion for the butterflies themselves. Heliconius butterfly eyes, which express recently duplicated UV opsins, have such an adaptation. To examine bird and butterfly color vision as sources of selection on butterfly coloration we studied yellow wing pigmentation in the tribe Heliconiini. We confirmed using reflectance and mass spectrometry that only Heliconius use 3-hydroxy-DL kynurenine (3-OHK) as a wing pigment. 3-OHK looks yellow to humans but it reflects both UV- and long-wavelength light whereas butterflies in related genera have chemically unknown yellow pigments mostly lacking UV-reflectance. Modeling of these color signals reveals that the two UV photoreceptors of Heliconius are better suited to separating 3-OHK from non-3-OHK spectra compared to the photoreceptors of related genera or birds. The co-occurrence of potentially enhanced UV-vision and a UV-reflecting ‘yellow’ wing pigment could allow unpalatable Heliconius private intraspecific communication in the presence of mimics. Our results are the best available evidence for the correlated evolution of a color signal and color vision. They also suggest that predator visual systems are error-prone in the context of mimicry

UV photoreceptors and UV-yellow wing pigments in Heliconius butterflies allow a color signal to serve both mimicry and intraspecific communication

Mimetic wing coloration evolves in butterflies in the context of predator confusion. Unless butterfly eyes have adaptations for discriminating mimetic color variation, mimicry also carries a risk of confusion for the butterflies themselves. Heliconius butterfly eyes, which express recently duplicated ultraviolet (UV) opsins, have such an adaptation. To examine bird and butterfly color vision as sources of selection on butterfly coloration, we studied yellow wing pigmentation in the tribe Heliconiini. We confirmed, using reflectance and mass spectrometry, that only Heliconius use 3-hydroxy-DL-kynurenine (3-OHK), which looks yellow to humans but reflects both UV- and long-wavelength light, whereas butterflies in related genera have chemically unknown yellow pigments mostly lacking UV reflectance. Modeling of these color signals reveals that the two UV photoreceptors of Heliconius are better suited to separating 3-OHK from non-3-OHK spectra compared with the photoreceptors of related genera or birds. The co-occurrence of potentially enhanced UV vision and a UV-reflecting yellow wing pigment could allow unpalatable Heliconius private intraspecific communication in the presence of mimics. Our results are the best available evidence for the correlated evolution of a color signal and color vision. They also suggest that predator visual systems are error prone in the context of mimicr

Inferring Identical-by-Descent Sharing of Sample Ancestors Promotes High-Resolution Relative Detection

As genetic datasets increase in size, the fraction of samples with one or more close relatives grows rapidly, resulting in sets of mutually related individuals. We present DRUID-deep relatedness utilizing identity by descent-a method that works by inferring the identical-by-descent (IBD) sharing profile of an ungenotyped ancestor of a set of close relatives. Using this IBD profile, DRUID infers relatedness between unobserved ancestors and more distant relatives, thereby combining information from multiple samples to remove one or more generations between the deep relationships to be identified. DRUID constructs sets of close relatives by detecting full siblings and also uses an approach to identify the aunts/uncles of two or more siblings, recovering 92.2% of real aunts/uncles with zero false positives. In real and simulated data, DRUID correctly infers up to 10.5% more relatives than PADRE when using data from two sets of distantly related siblings, and 10.7%-31.3% more relatives given two sets of siblings and their aunts/uncles. DRUID frequently infers relationships either correctly or within one degree of the truth, with PADRE classifying 43.3%-58.3% of tenth degree relatives in this way compared to 79.6%-96.7% using DRUID

Exome sequencing-based identification of novel type 2 diabetes risk allele loci in the Qatari population.

BACKGROUND:Type 2 diabetes (T2D) susceptibility is influenced by genetic and lifestyle factors. To date, the majority of genetic studies of T2D have been in populations of European and Asian descent. The focus of this study is on genetic variations underlying T2D in Qataris, a population with one of the highest incidences of T2D worldwide. RESULTS:Illumina HiSeq exome sequencing was performed on 864 Qatari subjects (574 T2D cases, 290 controls). Sequence kernel association test (SKAT) gene-based analysis identified an association for low frequency potentially deleterious variants in 6 genes. However, these findings were not replicated by SKAT analysis in an independent cohort of 12,699 exomes, primarly due to the absence of low frequency potentially deleterious variants in 5 of the 6 genes. Interestingly one of the genes identified, catenin beta 1 (CTNNB1, β-catenin), is the key effector of the Wnt pathway and interacts with the nuclear receptor transcription factor 7-like 2 (TCF7L2), variants which are the most strongly associated with risk of developing T2D worldwide. Single variant analysis did not identify any associated variants, suggesting the SKAT association signal was not driven by individual variants. None of the 6 associated genes were among 634 previously described T2D genes. CONCLUSIONS:The observation that genes not previously linked to T2D in prior studies of European and Asian populations are associated with T2D in Qatar provides new insights into the complexity of T2D pathogenesis and emphasizes the importance of understudied populations when assessing genetic variation in the pathogenesis of common disorders

Type 2 Diabetes Risk Allele Loci in the Qatari Population.

The prevalence of type 2 diabetes (T2D) is increasing in the Middle East. However, the genetic risk factors for T2D in the Middle Eastern populations are not known, as the majority of studies of genetic risk for T2D are in Europeans and Asians.All subjects were ≥3 generation Qataris. Cases with T2D (n = 1,124) and controls (n = 590) were randomly recruited and assigned to the 3 known Qatari genetic subpopulations [Bedouin (Q1), Persian/South Asian (Q2) and African (Q3)]. Subjects underwent genotyping for 37 single nucleotide polymorphisms (SNPs) in 29 genes known to be associated with T2D in Europeans and/or Asian populations, and an additional 27 tag SNPs related to these susceptibility loci. Pre-study power analysis suggested that with the known incidence of T2D in adult Qataris (22%), the study population size would be sufficient to detect significant differences if the SNPs were risk factors among Qataris, assuming that the odds ratio (OR) for T2D SNPs in Qatari's is greater than or equal to the SNP with highest known OR in other populations.Haplotype analysis demonstrated that Qatari haplotypes in the region of known T2D risk alleles in Q1 and Q2 genetic subpopulations were similar to European haplotypes. After Benjamini-Hochberg adjustment for multiple testing, only two SNPs (rs7903146 and rs4506565), both associated with transcription factor 7-like 2 (TCF7L2), achieved statistical significance in the whole study population. When T2D subjects and control subjects were assigned to the known 3 Qatari subpopulations, and analyzed individually and with the Q1 and Q2 genetic subpopulations combined, one of these SNPs (rs4506565) was also significant in the admixed group. No other SNPs associated with T2D in all Qataris or individual genetic subpopulations.With the caveats of the power analysis, the European/Asian T2D SNPs do not contribute significantly to the high prevalence of T2D in the Qatari population, suggesting that the genetic risks for T2D are likely different in Qataris compared to Europeans and Asians

Correction: Type 2 Diabetes Risk Allele Loci in the Qatari Population.

[This corrects the article DOI: 10.1371/journal.pone.0156834.]