Search

Population genetic analyses of local ancestry tracts routinely assume that the ancestral admixture process is identical for both parents of an individual, an assumption that may be invalid when considering recent admixture. Here we present Parental Admixture Proportion Inference (PAPI), a Bayesian tool for inferring the admixture proportions and admixture times for each parent of a single admixed individual. PAPI analyzes unphased local ancestry tracts and has two components models: a binomial model that exploits the informativeness of homozygous ancestry regions to infer parental admixture proportions, and a hidden Markov model (HMM) that infers admixture times from tract lengths. Crucially, the HMM employs an approximation to the pedigree crossover dynamics that accounts for unobserved within-ancestry recombination, enabling inference of parental admixture times. We compared the accuracy of PAPI’s admixture proportion estimates with those of ANCESTOR in simulated admixed individuals and found that PAPI outperforms ANCESTOR by an average of 46% in a representative set of simulation scenarios, with PAPI’s estimates deviating from the ground truth by 0.047 on average. Moreover, PAPI’s admixture time estimates were strongly correlated with the ground truth in these simulations (R = 0.76), but have an average downward bias of 1.01 generations that is partly attributable to inaccuracies in local ancestry inference. As an illustration of its utility, we ran PAPI on real African Americans from the PAGE study (N = 5, 786) and found strong evidence of assortative mating by ancestry proportion: couples’ ancestry proportions are closer to each other than expected by chance (P < 10−6), and are highly correlated (R = 0.87). We anticipate that PAPI will be useful in studying the population dynamics of admixture and will also be of interest to individuals seeking to learn about their personal genealogies.

Summary The proportion of samples with one or more close relatives in a genetic dataset increases rapidly with sample size, necessitating relatedness modeling and enabling pedigree-based analyses. Despite this, relatives are generally unreported and current inference methods typically detect only the degree of relatedness of sample pairs and not pedigree relationships. We developed CREST, an accurate and fast method that identifies the pedigree relationships of close relatives. CREST utilizes identity by descent (IBD) segments shared between a pair of samples and their mutual relatives, leveraging the fact that sharing rates among these individuals differ across pedigree configurations. Furthermore, CREST exploits the profound differences in sex-specific genetic maps to classify pairs as maternally or paternally related—e.g., paternal half-siblings—using the locations of autosomal IBD segments shared between the pair. In simulated data, CREST correctly classifies 91.5%–100% of grandparent-grandchild (GP) pairs, 80.0%–97.5% of avuncular (AV) pairs, and 75.5%–98.5% of half-siblings (HS) pairs compared to PADRE’s rates of 38.5%–76.0% of GP, 60.5%–92.0% of AV, 73.0%–95.0% of HS pairs. Turning to the real 20,032 sample Generation Scotland (GS) dataset, CREST identified seven pedigrees with incorrect relationship types or maternal/paternal parent sexes, five of which we confirmed as mistakes, and two with uncertain relationships. After correcting these, CREST correctly determines relationship types for 93.5% of GP, 97.7% of AV, and 92.2% of HS pairs that have sufficient mutual relative data; the parent sex in 100% of HS and 99.6% of GP pairs; and it completes this analysis in 2.8 h including IBD detection in eight threads.

Simulations of close relatives and identical by descent (IBD) segments are common in genetic studies, yet most past efforts have utilized sex averaged genetic maps and ignored crossover interference, thus omitting features known to affect the breakpoints of IBD segments. We developed Ped-sim, a method for simulating relatives that can utilize either sex-specific or sex averaged genetic maps and also either a model of crossover interference or the traditional Poisson model for inter-crossover distances. To characterize the impact of previously ignored mechanisms, we simulated data for all four combinations of these factors. We found that modeling crossover interference decreases the standard deviation of pairwise IBD proportions by 10.4% on average in full siblings through second cousins. By contrast, sex-specific maps increase this standard deviation by 4.2% on average, and also impact the number of segments relatives share. Most notably, using sex-specific maps, the number of segments half-siblings share is bimodal; and when combined with interference modeling, the probability that sixth cousins have non-zero IBD sharing ranges from 9.0 to 13.1%, depending on the sexes of the individuals through which they are related. We present new analytical results for the distributions of IBD segments under these models and show they match results from simulations. Finally, we compared IBD sharing rates between simulated and real relatives and find that the combination of sex-specific maps and interference modeling most accurately captures IBD rates in real data. Ped-sim is open source and available from https://github.com/williamslab/ped-sim., Author summary Simulations are ubiquitous throughout statistical genetics in order to generate data with known properties, enabling tests of inference methods and analyses of real world processes in settings where experimental data are challenging to collect. Simulating genetic data for relatives in a pedigree requires the synthesis of chromosomes parents transmit to their children. These chromosomes form as a mosaic of a given parent’s two chromosomes, with the location of switches between the two parental chromosomes known as crossovers. Detailed information about crossover generation based on real data from humans now exists, including the fact that men and women have overall different rates (women produce ∼1.6 times more crossovers) and that real crossovers are subject to interference—whereby crossovers are further apart from one another than expected under a model that selects their locations randomly. Our new method, Ped-sim, can simulate pedigree data using these less commonly modeled crossover features, and we used it to evaluate the impact of sex-specific rates and interference compared to real data. These comparisons show that both factors shape the amount of DNA two relatives share, and that their inclusion in models of crossover better fit data from real relatives.

Relatedness inference is an essential component of many genetic analyses and popular in consumer genetic testing. Ramstetter et al. evaluate twelve....., 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 76% 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.

Recombination has essential functions in evolution, meiosis, and breeding. Here, we use the multi-generational pedigree, consisting of 7,165 informative meioses (3,679 female; 3,486 male), and genotyping-by-sequencing (GBS) data from the International Institute of Tropical Agriculture (IITA) to study recombination in cassava (Manihot esculenta). We detected recombination events using SHAPEIT2 and duoHMM, examined the recombination landscape across the 18 chromosomes of cassava and in regions with known introgressed segments from cassava’s wild relative Manihot glaziovii, constructed a genetic map and compared it to an existing map constructed by the International Cassava Genetic Map Consortium (ICGMC), and inspected patterns of recombination placement in male and female meioses to see if there is evidence of sexual dimorphism in crossover distribution and frequency. We found that the placement of crossovers along chromosomes did not vary between the two sexes but that females undergo more meiotic recombination than males. We also observed that introgressions from M. glaziovii decreased recombination in the introgressed region and, in the case of chromosome 4, along the entire length of the chromosome that the introgression is on. We observed a dosage effect on chromosome 1, possibly suggesting the presence of a variant on the M. glaziovii haplotype that leads to lower overall recombination in the introgressed region.

Simulations of close relatives and identical by descent (IBD) segments are common in genetic studies, yet most past efforts have utilized sex averaged genetic maps and ignored crossover interference, thus omitting features known to affect the breakpoints of IBD segments. We developed Ped-sim, a method for simulating relatives that can utilize either sex-specific or sex averaged genetic maps and also either a model of crossover interference or the traditional Poisson model for inter-crossover distances. To characterize the impact of previously ignored mechanisms, we simulated data for all four combinations of these factors. We found that modeling crossover interference decreases the standard deviation of the IBD proportion by 10.4% on average in full siblings through second cousins. By contrast, sex-specific maps increase this standard deviation by 4.2% on average, and also impact the number of segments relatives share. Most notably, using sex-specific maps, the number of segments half-siblings share is bimodal; and when combined with interference modeling, the probability that sixth cousins have non-zero IBD ranges from 9.0 to 13.1%, depending on the sexes of the individuals through which they are related. We present new analytical results for the distributions of IBD segments under these models and show they match results from simulations. Finally, we compared IBD sharing rates between simulated and real relatives and find that the combination of sex-specific maps and interference modeling most accurately captures IBD rates in real data. Ped-sim is open source and available fromhttps://github.com/williamslab/ped-sim.Author summarySimulations are ubiquitous throughout statistical genetics in order to generate data with known properties, enabling tests of inference methods and analyses of real world processes in settings where experimental data are challenging to collect. Simulating genetic data for relatives in a pedigree requires the synthesis of chromosomes parents transmit to their children. These chromosomes form as a mosaic of a given parent’s two chromosomes, with the location of switches between the two parental chromosomes known as crossovers. Detailed information about crossover generation based on real data from humans now exists, including the fact that men and women have overall different rates (women produce ~1.6 times more crossovers) and that real crossovers are subject tointerference—whereby crossovers are further apart from one another than expected under a model that selects their locations randomly. Our new method, Ped-sim, can simulate pedigree data using these less commonly modeled crossover features, and we used it to evaluate the importance of sex-specific rates and interference in real data. These comparisons show that both factors shape the amount of DNA two relatives share identically, and that their inclusion in models of crossover better fit data from real relatives.

© 2015 by the American Diabetes Association. Insulin sensitivity, insulin secretion, insulin clearance, and glucose effectiveness exhibit strong genetic components, although few studies have examined their genetic architecture or influence on type 2 diabetes (T2D) risk. We hypothesized that loci affecting variation in these quantitative traits influence T2D. We completed a multicohort genome-wide association study to search for loci influencing T2D-related quantitative traits in 4,176 Mexican Americans. Quantitative traits were measured by the frequently sampled intravenous glucose tolerance test (four cohorts) or euglycemic clamp (three cohorts), and random-effects models were used to test the association between loci and quantitative traits, adjusting for age, sex, and admixture proportions (Discovery). Analysis revealed a significant (P < 5.00 3 1028) association at 11q14.3 (MTNR1B) with acute insulin response. Loci with P < 0.0001 among the quantitative traits were examined for translation to T2D risk in 6,463 T2D case and 9,232 control subjects of Mexican ancestry (Translation). Nonparametric meta- Analysis of the Discovery and Translation cohorts identified significant associations at 6p24 (SLC35B3/TFAP2A) with glucose effectiveness/T2D, 11p15 (KCNQ1) with disposition index/T2D, and 6p22 (CDKAL1) and 11q14 (MTNR1B) with acute insulin response/T2D. These results suggest that T2D and insulin secretion and sensitivity have both shared and distinct genetic factors, potentially delineating genomic components of these quantitative traits that drive the risk for T2D.

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 methods in real data has been lacking. Here, we report an assessment of 11 state-ofthe-art relatedness inference methods using a dataset with 2,485 individuals contained in several large pedigrees that span up to six generations. We find that all methods have high accuracy (∼93%−99%) when reporting first and second degree relationships, but their accuracy dwindles to less than 60% for fifth degree relationships. However, the inferred relationships were correct to within one relatedness degree at a rate of 83%−99% across all methods and considered relationship degrees. Furthermore, most methods infer unrelated individuals correctly at a rate of ∼99%, suggesting a low rate of false positives. Overall, the most accurate methods were ERSA 2.0 and approaches that classify relationships using the IBD segments inferred by Refined IBD and IBDseq. Combining results from the most accurate methods provides little accuracy improvement, indicating that novel approaches for relatedness inference may be needed to achieve a sizeable jump in performance.

Importance: Latino populations have one of the highest prevalences of type 2 diabetes worldwide. Objectives: To investigate the association between rare protein-coding genetic variants and prevalence of type 2 diabetes in a large Latino population and to explore potential molecular and physiological mechanisms for the observed relationships. Design, Setting, and Participants: Whole-exome sequencing was performed on DNA samples from 3756 Mexican and US Latino individuals (1794 with type 2 diabetes and 1962 without diabetes) recruited from 1993 to 2013. One variant was further tested for allele frequency and association with type 2 diabetes in large multiethnic data sets of 14 276 participants and characterized in experimental assays. Main Outcome and Measures: Prevalence of type 2 diabetes. Secondary outcomes included age of onset, body mass index, and effect on protein function. Results: A single rare missense variant (c.1522G>A [p.E508K]) was associated with type 2 diabetes prevalence (odds ratio [OR], 5.48; 95% CI, 2.83-10.61; P = 4.4 × 10−7) in hepatocyte nuclear factor 1-α (HNF1A), the gene responsible for maturity onset diabetes of the young type 3 (MODY3). This variant was observed in 0.36% of participants without type 2 diabetes and 2.1% of participants with it. In multiethnic replication data sets, the p.E508K variant was seen only in Latino patients (n = 1443 with type 2 diabetes and 1673 without it) and was associated with type 2 diabetes (OR, 4.16; 95% CI, 1.75-9.92; P = .0013). In experimental assays, HNF-1A protein encoding the p.E508K mutant demonstrated reduced transactivation activity of its target promoter compared with a wild-type protein. In our data, carriers and noncarriers of the p.E508K mutation with type 2 diabetes had no significant differences in compared clinical characteristics, including age at onset. The mean (SD) age for carriers was 45.3 years (11.2) vs 47.5 years (11.5) for noncarriers (P = .49) and the mean (SD) BMI for carriers was 28.2 (5.5) vs 29.3 (5.3) for noncarriers (P = .19). Conclusions and Relevance: Using whole-exome sequencing, we identified a single low-frequency variant in the MODY3-causing gene HNF1A that is associated with type 2 diabetes in Latino populations and may affect protein function. This finding may have implications for screening and therapeutic modification in this population, but additional studies are required. publishedVersion

Hapi is a new dynamic programming algorithm that ignores uninformative states and state transitions in order to efficiently compute minimum-recombinant and maximum likelihood haplotypes. When applied to a dataset containing 103 families, Hapi performs 3.8 and 320 times faster than state-of-the-art algorithms. Because Hapi infers both minimum-recombinant and maximum likelihood haplotypes and applies to related individuals, the haplotypes it infers are highly accurate over extended genomic distances., National Institutes of Health (U.S.) (NIH grant 5-T90-DK070069), National Institutes of Health (U.S.) (Grant 5-P01-NS055923), National Science Foundation (U.S.) (Graduate Research Fellowship)