Dense Marker Data Research Articles

Genotyping-by-sequencing and genomic selection applications in hexaploid triticale.

Triticale, a hybrid species between wheat and rye, is one of the newest additions to the plant kingdom with a very short history of improvement. It has very limited genomic resources because of its large and complex genome. Objectives of this study were to generate dense marker data, understand genetic diversity, population structure, linkage disequilibrium (LD), and estimate accuracies of commonly used genomic selection (GS) models on forage yield of triticale. Genotyping-by-sequencing (GBS), using PstI and MspI restriction enzymes for reducing genome complexity, was performed on a triticale diversity panel (n = 289). After filtering for biallelic loci with more than 70% genome coverage, and minor allele frequency (MAF) > 0.05, de novo variant calling identified 16,378 single nucleotide polymorphism (SNP) markers. Sequences of these variants were mapped to wheat and rye reference genomes to infer their homologous groups and chromosome positions. About 45% (7430), and 58% (9500) of the de novo identified SNPs were mapped to the wheat and rye reference genomes, respectively. Interestingly, 28.9% (2151) of the 7430 SNPs were mapped to the D genome of hexaploid wheat, indicating substantial substitution of the R genome with D genome in cultivated triticale. About 27% of marker pairs were in significant LD with an average r2 > 0.18 (P < 0.05). Genome-wide LD declined rapidly to r2 < 0.1 beyond 10 kb physical distance. The three sub-genomes (A, B, and R) showed comparable LD decay patterns. Genetic diversity and population structure analyses identified five distinct clusters. Genotype grouping did not follow prior winter vs spring-type classification. However, one of the clusters was largely dominated by winter triticale. GS accuracies were estimated for forage yield using three commonly used models with different training population sizes and marker densities. GS accuracy increased with increasing training population size while gain in accuracy tended to plateau with marker densities of 2000 SNPs or more. Average GS accuracy was about 0.52, indicating the potential of using GS in triticale forage yield improvement.

A generic hidden Markov model for multiparent populations

A common step in the analysis of multiparent populations (MPPs) is genotype reconstruction: identifying the founder origin of haplotypes from dense marker data. This process often makes use of a probability model for the pattern of founder alleles along chromosomes, including the relative frequency of founder alleles and the probability of exchanges among them, which depend on a model for meiotic recombination and on the mating design for the population. While the precise experimental design used to generate the population may be used to derive a precise characterization of the model for exchanges among founder alleles, this can be tedious, particularly given the great variety of experimental designs that have been proposed. We describe an approximate model that can be applied for a variety of MPPs. We have implemented the approach in the R/qtl2 software, and we illustrate its use in applications to publicly available data on Diversity Outbred and Collaborative Cross mice.

Genomic prediction models trained with historical records enable populating the German ex situ genebank bio-digital resource center of barley (Hordeum sp.) with information on resistances to soilborne barley mosaic viruses

Key messageGenomic prediction with special weight of major genes is a valuable tool to populate bio-digital resource centers.Phenotypic information of crop genetic resources is a prerequisite for an informed selection that aims to broaden the genetic base of the elite breeding pools. We investigated the potential of genomic prediction based on historical screening data of plant responses against the Barley yellow mosaic viruses for populating the bio-digital resource center of barley. Our study includes dense marker data for 3838 accessions of winter barley, and historical screening data of 1751 accessions for Barley yellow mosaic virus (BaYMV) and of 1771 accessions for Barley mild mosaic virus (BaMMV). Linear mixed models were fitted by considering combinations for the effects of genotypes, years, and locations. The best linear unbiased estimations displayed a broad spectrum of plant responses against BaYMV and BaMMV. Prediction abilities, computed as correlations between predictions and observed phenotypes of accessions, were low for the marker-assisted selection approach amounting to 0.42. In contrast, prediction abilities of genomic best linear unbiased predictions were high, with values of 0.62 for BaYMV and 0.64 for BaMMV. Prediction abilities of genomic prediction were improved by up to ~ 5% using W-BLUP, in which more weight is given to markers with significant major effects found by association mapping. Our results outline the utility of historical screening data and W-BLUP model to predict the performance of the non-phenotyped individuals in genebank collections. The presented strategy can be considered as part of the different approaches used in genebank genomics to valorize genetic resources for their usage in disease resistance breeding and research.

Genebank genomics highlights the diversity of a global barley collection.

Genebanks hold comprehensive collections of cultivars, landraces and crop wild relatives of all major food crops, but their detailed characterization has so far been limited to sparse core sets. The analysis of genome-wide genotyping-by-sequencing data for almost all barley accessions of the German ex situ genebank provides insights into the global population structure of domesticated barley and points out redundancies and coverage gaps in one of the world's major genebanks. Our large sample size and dense marker data afford great power for genome-wide association scans. We detect known and novel loci underlying morphological traits differentiating barley genepools, find evidence for convergent selection for barbless awns in barley and rice and show that a major-effect resistance locus conferring resistance to bymovirus infection has been favored by traditional farmers. This study outlines future directions for genomics-assisted genebank management and the utilization of germplasm collections for linking natural variation to human selection during crop evolution.

Oxytocin receptor gene polymorphisms, attachment, and PTSD: Results from the National Health and Resilience in Veterans Study

The human oxytocin system is implicated in social behavior and stress recovery. Polymorphisms in the oxytocin receptor gene (OXTR) may interact with attachment style to predict stress-related psychopathology like posttraumatic stress disorder (PTSD). The objective of this study was to examine independent and interactive effects of the OXTR single nucleotide polymorphism (SNP) rs53576, which has been associated with stress reactivity, support-seeking, and PTSD in prior studies, and attachment style on risk for PTSD in a nationally representative sample of 2163 European-American (EA) U.S. military veterans who participated in two independent waves of the National Health and Resilience in Veterans Study (NHRVS). Results revealed that insecure attachment style [adjusted odds ratio (OR) = 4.29; p < 0.001] and the interaction of rs53576 and attachment style (OR = 2.58, p = 0.02) were associated with probable lifetime PTSD. Among individuals with the minor A allele, the prevalence of probable PTSD was significantly higher among those with an insecure attachment style (23.9%) than those with a secure attachment style (2.0%), equivalent to an adjusted OR of 10.7. We attempted to replicate these findings by utilizing dense marker data from a genome-wide association study of 2215 high-risk civilians; one OXTR variant, though not rs53576, was associated with PTSD. Exploratory analyses in the veteran sample revealed that the interaction between this variant and attachment style predicting probable PTSD approached statistical significance. Results indicate that polymorphisms in the OXTR gene and attachment style may contribute to vulnerability to PTSD in U.S. military veterans.

Response and inbreeding from a genomic selection experiment in layer chickens.

BackgroundGenomic selection (GS) using estimated breeding values (GS-EBV) based on dense marker data is a promising approach for genetic improvement. A simulation study was undertaken to illustrate the opportunities offered by GS for designing breeding programs. It consisted of a selection program for a sex-limited trait in layer chickens, which was developed by deterministic predictions under different scenarios. Later, one of the possible schemes was implemented in a real population of layer chicken.MethodsIn the simulation, the aim was to double the response to selection per year by reducing the generation interval by 50 %, while maintaining the same rate of inbreeding per year. We found that GS with retraining could achieve the set objectives while requiring 75 % fewer reared birds and 82 % fewer phenotyped birds per year. A multi-trait GS scenario was subsequently implemented in a real population of brown egg laying hens. The population was split into two sub-lines, one was submitted to conventional phenotypic selection, and one was selected based on genomic prediction. At the end of the 3-year experiment, the two sub-lines were compared for multiple performance traits that are relevant for commercial egg production.ResultsBirds that were selected based on genomic prediction outperformed those that were submitted to conventional selection for most of the 16 traits that were included in the index used for selection. However, although the two programs were designed to achieve the same rate of inbreeding per year, the realized inbreeding per year assessed from pedigree was higher in the genomic selected line than in the conventionally selected line.ConclusionsThe results demonstrate that GS is a promising alternative to conventional breeding for genetic improvement of layer chickens.

ALDsuite: Dense marker MALD using principal components of ancestral linkage disequilibrium.

BackgroundMapping by admixture linkage disequilibrium (MALD) is a whole genome gene mapping method that uses LD from extended blocks of ancestry inherited from parental populations among admixed individuals to map associations for diseases, that vary in prevalence among human populations. The extended LD queried for marker association with ancestry results in a greatly reduced number of comparisons compared to standard genome wide association studies. As ancestral population LD tends to confound the analysis of admixture LD, the earliest algorithms for MALD required marker sets sufficiently sparse to lack significant ancestral LD between markers. However current genotyping technologies routinely provide dense SNP data, which convey more information than sparse sets, if this information can be efficiently used. There are currently no software solutions that offer both local ancestry inference using dense marker data and disease association statistics.ResultsWe present here an R package, ALDsuite, which accounts for local LD using principal components of haplotypes from surrogate ancestral population data, and includes tools for quality control of data, MALD, downstream analysis of results and visualization graphics.ConclusionsALDsuite offers a fast, accurate estimation of global and local ancestry and comes bundled with the tools needed for MALD, from data quality control through mapping of and visualization of disease genes.

Dating Rare Mutations from Small Samples with Dense Marker Data

We present a method for estimating the age of a mutation based on the genetic length of ancestral haplotypes shared between individuals carrying the mutation. The method can be reliably applied to small samples, typical of situations involving rare mutations, and makes effective use of modern high-density SNP data, thus overcoming two of the limitations with existing methods. The method provides age estimates and confidence intervals without the use of asymptotic theory and is applicable to genealogies in which the data are independent or correlated. In the correlated case we estimate the correlation directly from the data, rather than relying on a model for the genealogy. To demonstrate the method's efficacy, we provide simulation results and compare it to other methods. The length data are obtained with a simple procedure, and an R script is available for performing the calculations.

Practical Considerations Regarding the Use of Genotype and Pedigree Data to Model Relatedness in the Context of Genome-Wide Association Studies

Genome-wide association studies of complex traits often are complicated by relatedness among individuals. Ignoring or inappropriately accounting for relatedness often results in inflated type I error rates. Either genotype or pedigree data can be used to estimate relatedness for use in mixed-models when undertaking quantitative trait locus mapping. We performed simulations to investigate methods for controlling type I error and optimizing power considering both full and partial pedigrees and, similarly, both sparse and dense marker coverage; we also examined real data sets. (1) When marker density was low, estimating relatedness by genotype data alone failed to control the type I error rate; (2) this was resolved by combining both genotype and pedigree data. (3) When sufficiently dense marker data were used to estimate relatedness, type I error was well controlled and power increased; however, (4) this was only true when the relatedness was estimated using genotype data that excluded genotypes on the chromosome currently being scanned for a quantitative trait locus.

Identification of Pedigree Relationship from Genome Sharing

Determination of degree of relationship traditionally has been undertaken using genotypic data on individual loci, typically assumed to be independent. With dense marker data as now available, it is possible to identify the regions of the genome shared identical by descent (ibd). This information can be used to determine pedigree relationship (R), e.g., cousins vs. second cousins, and also to distinguish pedigrees that have the same Wright’s relationship (R) such as half-sibs and uncle–nephew. We use simulation to investigate the accuracy with which pedigree relationship can be inferred from genome sharing for uniparental relatives (a common ancestor on only one side of their pedigree), specifically the number, position (whether at chromosome ends), and length of shared regions ibd on each chromosome. Moments of the distribution of the likelihood ratio (including its expectation, the Kullback-Leibler distance) for alternative relationships are estimated for model human genomes, with the ratio of the mean to the SD of the likelihood ratio providing a useful reference point. Two relationships differing in R can be readily distinguished provided at least one has high R, e.g., approximately 98.5% correct assignment of cousins and half-cousins, but only approximately 75% for second cousins once removed and third cousins. Two relationships with the same R can be distinguished only if R is high, e.g., half-sibs and uncle–nephew, with probability of correct assignment being approximately 5/6.

Genomic Selection: Status in Different Species and Challenges for Breeding

Technical advances and development in the market for genomic tools have facilitated access to whole-genome data across species. Building-up on the acquired knowledge of the genome sequences, large-scale genotyping has been optimized for broad use, so genotype information can be routinely used to predict genetic merit. Genomic selection (GS) refers to the use of aggregates of estimated marker effects as predictors which allow improved individual differentiation at young age. Realizable benefits of GS are influenced by several factors and vary in quantity and quality between species. General characteristics and challenges of GS in implementation and routine application are described, followed by an overview over the current status of its use, prospects and challenges in important animal species. Genetic gain for a particular trait can be enhanced by shortening of the generation interval, increased selection accuracy and increased selection intensity, with species- and breed-specific relevance of the determinants. Reliable predictions based on genetic marker effects require assembly of a reference for linking of phenotype and genotype data to allow estimation and regular re-estimation. Experiences from dairy breeding have shown that international collaboration can set the course for fast and successful implementation of innovative selection tools, so genomics may significantly impact the structures of future breeding and breeding programmes. Traits of great and increasing importance, which were difficult to improve in the conventional systems, could be emphasized, if continuous availability of high-quality phenotype data can be assured. Equally elaborate strategies for genotyping and phenotyping will allow tailored approaches to balance efficient animal production, sustainability, animal health and welfare in future.

Inferring coancestry in population samples in the presence of linkage disequilibrium.

In both pedigree linkage studies and in population-based association studies there has been much interest in the use of modern dense genetic marker data to infer segments of gene identity by descent (ibd) among individuals not known to be related, to increase power and resolution in localizing genes affecting complex traits. In this article, we present a hidden Markov model (HMM) for ibd among a set of chromosomes and describe methods and software for inference of ibd among the four chromosomes of pairs of individuals, using either phased (haplotypic) or unphased (genotypic) data. The model allows for missing data and typing error, but does not model linkage disequilibrium (LD), because fitting an accurate LD model requires large samples from well-studied populations. However, LD remains a major confounding factor, since LD is itself a reflection of coancestry at the population level. To study the impact of LD, we have developed a novel simulation approach to generate realistic dense marker data for the same set of markers but at varying levels of LD. Using this approach, we present results of a study of the impact of LD on the sensitivity and specificity of our HMM model in estimating segments of ibd among sets of four chromosomes and between genotype pairs. We show that, despite not incorporating LD, our model has been quite successful in detecting segments as small as 10(6) bp (1 Mpb); we present also comparisons with fastIBD which uses an LD model in estimating ibd.

Development of admixture mapping panels for African Americans from commercial high-density SNP arrays

BackgroundAdmixture mapping is a powerful approach for identifying genetic variants involved in human disease that exploits the unique genomic structure in recently admixed populations. To use existing published panels of ancestry-informative markers (AIMs) for admixture mapping, markers have to be genotyped de novo for each admixed study sample and samples representing the ancestral parental populations. The increased availability of dense marker data on commercial chips has made it feasible to develop panels wherein the markers need not be predetermined.ResultsWe developed two panels of AIMs (~2,000 markers each) based on the Affymetrix Genome-Wide Human SNP Array 6.0 for admixture mapping with African American samples. These two AIM panels had good map power that was higher than that of a denser panel of ~20,000 random markers as well as other published panels of AIMs. As a test case, we applied the panels in an admixture mapping study of hypertension in African Americans in the Washington, D.C. metropolitan area.ConclusionsDeveloping marker panels for admixture mapping from existing genome-wide genotype data offers two major advantages: (1) no de novo genotyping needs to be done, thereby saving costs, and (2) markers can be filtered for various quality measures and replacement markers (to minimize gaps) can be selected at no additional cost. Panels of carefully selected AIMs have two major advantages over panels of random markers: (1) the map power from sparser panels of AIMs is higher than that of ~10-fold denser panels of random markers, and (2) clusters can be labeled based on information from the parental populations. With current technology, chip-based genome-wide genotyping is less expensive than genotyping ~20,000 random markers. The major advantage of using random markers is the absence of ascertainment effects resulting from the process of selecting markers. The ability to develop marker panels informative for ancestry from SNP chip genotype data provides a fresh opportunity to conduct admixture mapping for disease genes in admixed populations when genome-wide association data exist or are planned.

Comparison of association methods for dense marker data

While data sets based on dense genome scans are becoming increasingly common, there are many theoretical questions that remain unanswered. How can a large number of markers in high linkage disequilibrium (LD) and rare disease variants be simulated efficiently? How should markers in high LD be analyzed: individually or jointly? Are there fast and simple methods to adjust for correlation of tests? What is the power penalty for conservative Bonferroni adjustments? Assuming that association scans are adequately powered, we attempt to answer these questions. Performance of single-point and multipoint tests, and their hybrids, is investigated using two simulation designs. The first simulation design uses theoretically derived LD patterns. The second design uses LD patterns based on real data. For the theoretical simulations we used polychoric correlation as a measure of LD to facilitate simulation of markers in LD and rare disease variants. Based on the simulation results of the two studies, we conclude that statistical tests assuming only additive genotype effects (i.e. Armitage and especially multipoint T(2)) should be used cautiously due to their suboptimal power in certain settings. A false discovery rate (FDR)-adjusted combination of tests for additive, dominant and recessive effects had close to optimal power. However, the common genotypic chi(2) test performed adequately and could be used in lieu of the FDR combination. While some hybrid methods yield (sometimes spectacularly) higher power they are computationally intensive. We also propose an "exact" method to adjust for multiple testing, which yields nominally higher power than the Bonferroni correction.

Estimation of Pairwise Identity by Descent From Dense Genetic Marker Data in a Population Sample of Haplotypes

I present a new approach for calculating probabilities of identity by descent for pairs of haplotypes. The approach is based on a joint hidden Markov model for haplotype frequencies and identity by descent (IBD). This model allows for linkage disequilibrium, and the method can be applied to very dense marker data. The method has high power for detecting IBD tracts of genetic length of 1 cM, with the use of sufficiently dense markers. This enables detection of pairwise IBD between haplotypes from individuals whose most recent common ancestor lived up to 50 generations ago.

Efficient computation of genotype probabilities for loci with many alleles: II. Iterative method for large, complex pedigrees.

An algorithm for computing genotype probabilities for marker loci with many alleles in large, complex pedigrees with missing marker data is presented. The algorithm can also be used to calculate grandparental origin probabilities, which summarize the segregation pattern and are useful for mapping quantitative trait loci. The algorithm is iterative and is based on peeling on alleles instead of the traditional peeling on genotypes. This makes the algorithm more computationally efficient for loci with many alleles. The algorithm is approximate in pedigrees that contain loops, including loops generated by full sibs. The algorithm has no restrictions on pedigree structure or missing marker phenotypes, although together those factors affect the degree of approximation. In livestock pedigrees with dense marker data, the degree of approximation may be minimal. The algorithm can be used with an incomplete penetrance model for marker loci. Thus, it takes into account the possibility of marker scoring errors and helps to identify them. The algorithm provides a computationally feasible method to analyze genetic marker data in large, complex livestock pedigrees.