Performance of single nucleotide polymorphisms versus haplotypes for genome-wide association analysis in barley

Genome-wide association studies (GWAS) may benefit from utilizing haplotype information for making marker-phenotype associations. Several rationales for grouping single nucleotide polymorphisms (SNPs) into haplotype blocks exist, but any advantage may depend on such factors as genetic architecture of traits, patterns of linkage disequilibrium in the study population, and marker density. The objective of this study was to explore the utility of haplotypes for GWAS in barley (Hordeum vulgare) to offer a first detailed look at this approach for identifying agronomically important genes in crops. To accomplish this, we used genotype and phenotype data from the Barley Coordinated Agricultural Project and constructed haplotypes using three different methods. Marker-trait associations were tested by the efficient mixed-model association algorithm (EMMA). When QTL were simulated using single SNPs dropped from the marker dataset, a simple sliding window performed as well or better than single SNPs or the more sophisticated methods of blocking SNPs into haplotypes. Moreover, the haplotype analyses performed better 1) when QTL were simulated as polymorphisms that arose subsequent to marker variants, and 2) in analysis of empirical heading date data. These results demonstrate that the information content of haplotypes is dependent on the particular mutational and recombinational history of the QTL and nearby markers. Analysis of the empirical data also confirmed our intuition that the distribution of QTL alleles in nature is often unlike the distribution of marker variants, and hence utilizing haplotype information could capture associations that would elude single SNPs. We recommend routine use of both single SNP and haplotype markers for GWAS to take advantage of the full information content of the genotype data.     

Genome-wide association analysis of meat quality traits in a porcine Large White × Minzhu intercross population

Pork quality is an economically important trait and one of the main selection criteria for breeding in the swine industry. In this genome-wide association study (GWAS), 455 pigs from a porcine Large White × Minzhu intercross population were genotyped using the Illumina PorcineSNP60K Beadchip, and phenotyped for intramuscular fat content (IMF), marbling, moisture, color L*, color a*, color b* and color score in the longissimus muscle (LM). Association tests between each trait and the SNPs were performed via the Genome Wide Rapid Association using the Mixed Model and Regression-Genomic Control (GRAMMAR-GC) approach. From the Ensembl porcine database, SNP annotation was implemented using Sus scrofa Build 9. A total of 45 SNPs showed significant association with one or multiple meat quality traits. Of the 45 SNPs, 36 were located on SSC12. These significantly associated SNPs aligned to or were in close approximation to previously reported quantitative trait loci (QTL) and some were located within introns of previously reported candidate genes. Two haplotype blocks ASGA0100525-ASGA0055225-ALGA0067099-MARC0004712-DIAS0000861, and ASGA0085522-H3GA0056170 were detected in the significant region. The first block contained the genes MYH1, MYH2 and MYH4. A SNP (ASGA0094812) within an intron of the USP43 gene was significantly associated with five meat quality traits. The present results effectively narrowed down the associated regions compared to previous QTL studies and revealed haplotypes and candidate genes on SSC12 for meat quality traits in pigs.     

Tree scanning: a method for using haplotype trees in phenotype/genotype association studies

We use evolutionary trees of haplotypes to study phenotypic associations by exhaustively examining all possible biallelic partitions of the tree, a technique we call tree scanning. If the first scan detects significant associations, additional rounds of tree scanning are used to partition the tree into three or more allelic classes. Two worked examples are presented. The first is a reanalysis of associations between haplotypes at the Alcohol Dehydrogenase locus in Drosophila melanogaster that was previously analyzed using a nested clade analysis, a more complicated technique for using haplotype trees to detect phenotypic associations. Tree scanning and the nested clade analysis yield the same inferences when permutation testing is used with both approaches. The second example is an analysis of associations between variation in various lipid traits and genetic variation at the Apolipoprotein E (APOE) gene in three human populations. Tree scanning successfully identified phenotypic associations expected from previous analyses. Tree scanning for the most part detected more associations and provided a better biological interpretative framework than single SNP analyses. We also show how prior information can be incorporated into the tree scan by starting with the traditional three electrophoretic alleles at APOE. Tree scanning detected genetically determined phenotypic heterogeneity within all three electrophoretic allelic classes. Overall, tree scanning is a simple, powerful, and flexible method for using haplotype trees to detect phenotype/genotype associations at candidate loci.     

Gene SNPs and mutations in clinical genetic testing: haplotype-based testing and analysis

Haplotype-based analysis using high-density single nucleotide polymorphism (SNP) markers have gained increasing attention in evaluating candidate genes in various clinical situations. For example, haplotype information is useful for predicting the severity and prognosis of certain genetic disorders. The intragenic cis-interactions between the common polymorphisms and the pathogenic mutations of prion protein (PRNP) and cystic fibrosis transmembrane conductance regulator (CFTR) genes greatly influence the phenotypes and the disease penetrance of hereditary Creutzfeldt-Jakob disease and cystic fibrosis. Merits of haplotype study are more evident in the fine mapping of complex diseases and in identifying genetic variations that influence individual's response to drugs. Knowledge-based approaches and/or linkage analyses using SNP tagged haplotypes are effective tools in detecting genetic associations. For example, haplotype studies in the inflammatory bowel disease susceptibility loci revealed diverse cis and trans gene-gene interactions, which can affect the clinical outcomes. Although currently, we have very limited knowledge on haplotype-phenotypic characterizations of most genes, these examples demonstrate that increased understanding of the clinically relevant haplotypes will provide better results in the diagnosis and possibly in the treatment of both monogenic and polygenic diseases.     

Test for positional candidate genes for body composition on pig chromosome 6

One QTL affecting backfat thickness (BF), intramuscular fat content (IMF) and eye muscle area (MA) was previously localized on porcine chromosome 6 in an F₂ cross between Iberian and Landrace pigs. This work was done to study the effect of two positional candidate genes on these traits: H-FABP and LEPR genes. The QTL mapping analysis was repeated with a regression method using genotypes for seven microsatellites and two PCR-RFLPs in the H-FABP and LEPR genes. H-FABP and LEPR genes were located at 85.4 and 107 cM respectively, by linkage analysis. The effects of the candidate gene polymorphisms were analyzed in two ways. When an animal model was fitted, both genes showed significant effects on fatness traits, the H-FABP polymorphism showed significant effects on IMF and MA, and the LEPR polymorphism on BF and IMF. But when the candidate gene effect was included in a QTL regression analysis these associations were not observed, suggesting that they must not be the causal mutations responsible for the effects found. Differences in the results of both analyses showed the inadequacy of the animal model approach for the evaluation of positional candidate genes in populations with linkage disequilibrium, when the probabilities of the parental origin of the QTL alleles are not included in the model.     

Risk haplotype analysis for bovine paratuberculosis

Paratuberculosis (Johne’s disease), caused by Mycobacterium avium subsp. paratuberculosis, is an important disease for bovines, although its genetic basis is poorly understood. In this study, three candidate genes were typed to study the associations between single nucleotide polymorphisms (SNPs) and paratuberculosis susceptibility (measured in a 1 or 0 form) at the haplotype level. A significant risk haplotype, constructed by a variant allele (C) at the first SNP and a common allele (A) at the second SNP, within the CARD15 gene was detected to trigger genetic effects on paratuberculosis infection in an overdominace manner. Marginally significant haplotypes were identified for the other two genes. The results obtained will provide scientific guidance about the selection and prediction of resistance types in bovines.     

Association of candidate genes with flowering time and water-soluble carbohydrate content in Lolium perenne (L.)

We describe a candidate gene approach for associating SNPs with variation in flowering time and water-soluble carbohydrate (WSC) content and other quality traits in the temperate forage grass species Lolium perenne. Three analysis methods were used, which took the significant population structure into account. First, a linear mixed model was used enabling a structured association analysis to be incorporated with the nine populations identified in the structure analysis as random variables. Second, a within-population analysis of variance was performed. Third, a tree-scanning method was used, in which haplotype trees were associated with phenotypes on the basis of inferred haplotypes. Analysis of variance within populations identified several associations between WSC, nitrogen (N), and dry matter digestibility with allelic variants within an alkaline invertase candidate gene LpcAI. These associations were only detected in material harvested in one of the two years. By contrast, consistent associations between the L. perenne homolog (LpHD1) of the rice photoperiod control gene HD1 and flowering time were identified. One SNP, in the immediate upstream region of the LpHD1 coding sequence (C-4443-A), was significant in the linear mixed model. Within-population analysis of variance and tree-scanning analysis confirmed and extended this result to the 2118 polymorphisms in some of the populations. The merits of the tree-scanning method are compared to the single SNP analysis. The potential usefulness of the 4443 SNP in marker-assisted selection is currently being evaluated in test crosses of genotypes from this work with turf-grass varieties.     

Haplotype‐based genotyping‐by‐sequencing in oat genome research

In a de novo genotyping‐by‐sequencing (GBS) analysis of short, 64‐base tag‐level haplotypes in 4657 accessions of cultivated oat, we discovered 164741 tag‐level (TL) genetic variants containing 241224 SNPs. From this, the marker density of an oat consensus map was increased by the addition of more than 70000 loci. The mapped TL genotypes of a 635‐line diversity panel were used to infer chromosome‐level (CL) haplotype maps. These maps revealed differences in the number and size of haplotype blocks, as well as differences in haplotype diversity between chromosomes and subsets of the diversity panel. We then explored potential benefits of SNP vs. TL vs. CL GBS variants for mapping, high‐resolution genome analysis and genomic selection in oats. A combined genome‐wide association study (GWAS) of heading date from multiple locations using both TL haplotypes and individual SNP markers identified 184 significant associations. A comparative GWAS using TL haplotypes, CL haplotype blocks and their combinations demonstrated the superiority of using TL haplotype markers. Using a principal component‐based genome‐wide scan, genomic regions containing signatures of selection were identified. These regions may contain genes that are responsible for the local adaptation of oats to Northern American conditions. Genomic selection for heading date using TL haplotypes or SNP markers gave comparable and promising prediction accuracies of up to r = 0.74. Genomic selection carried out in an independent calibration and test population for heading date gave promising prediction accuracies that ranged between r = 0.42 and 0.67. In conclusion, TL haplotype GBS‐derived markers facilitate genome analysis and genomic selection in oat.     

复杂疾病的遗传易感基因区域的精细定位

以单核苷酸多态性(Single-nucleotide polymorphism, SNP)为遗传标记, 采用全基因组关联研究(Genome-wide association studies, GWAS)的策略, 已经在660多种疾病(或性状)中发现了3800多个遗传易感基因区域。但是, 其中最显著关联的遗传变异或致病性的遗传变异位点及其生物学功能并不完全清楚。这些位点的鉴定有助于阐明复杂疾病的生物学机制, 以及发现新的疾病标记物。后GWAS时代的主要任务之一就是通过精细定位研究找到复杂疾病易感基因区域内最显著关联的易感位点或致病性的易感位点并阐明其生物学功能。针对常见变异, 可通过推断或重测序增加SNP密度, 寻找最显著关联的SNP位点, 并通过功能元件分析、表达数量性状位点(Expression quantitative trait locus, eQTL)分析和单体型分析等方法寻找功能性的SNP位点和易感基因。针对罕见变异, 则可采用重测序、罕见单体型分析、家系分析和负荷检验等方法进行精细定位。文章对这些策略和所面临的问题进行了综述。     

Genome-wide association study-based identification genes influencing agronomic traits in rice (Oryza sativa L.)

Rice is one of the most important cereal crops, providing the daily dietary intake for approximately 50% of the global human population. Here, we re-sequenced 259 rice accessions, generating 1371.65 Gb of raw data. Furthermore, we performed genome-wide association studies (GWAS) on 13 agronomic traits using 2.8 million single nucleotide polymorphisms (SNPs) characterized in 259 rice accessions. Phenotypic data and best linear unbiased prediction (BLUP) values of each of the 13 traits over two years of each trait were used for the GWAS. The results showed that 816 SNP signals were significantly associated with the 13 agronomic traits. Then we detected candidate genes related to target traits within 200 kb upstream and downstream of the associated SNP loci, based on linkage disequilibrium (LD) blocks in the whole rice genome. These candidate genes were further identified through haplotype block constructions. This comprehensive study provides a timely and important genomic resource for breeding high yielding rice cultivars.     

Characterization of a QTL region affecting clinical mastitis and protein yield on BTA6

Quantitative trait loci affecting clinical mastitis were detected and fine mapped to a narrow region on bovine chromosome 6 in the Norwegian Red cattle population. The region includes the casein gene cluster and several candidate genes thought to influence clinical mastitis. The most significant results were found for SNPs within the Mucin 7 gene. This gene encodes an antimicrobial peptide and constitutes part of the first line of defence for the mucosal immune system. Detection of long haplotypes extending several Mb may indicate that artificial selection has influenced the haplotype structures in the region. A search for selection sweeps supports this observation and coincides with association results found both by single SNP and haplotype analyses. Our analyses identified haplotypes carrying quantitative trait loci alleles associated with high protein yield and simultaneously fewer incidences of clinical mastitis. The fact that such haplotypes are found in relative high frequencies in Norwegian Red may reflect the combined breeding goal that is characterized by selection for both milk production and disease resistance. The identification of these haplotypes raises the possibility of overcoming the unfavourable genetic correlation between these traits through haplotype-assisted selection.     

HapScope: a software system for automated and visual analysis of functionally annotated haplotypes

We have developed a software analysis package, HapScope, which includes a comprehensive analysis pipeline and a sophisticated visualization tool for analyzing functionally annotated haplotypes. The HapScope analysis pipeline supports: (i) computational haplotype construction with an expectation-maximization or Bayesian statistical algorithm; (ii) SNP classification by protein coding change, homology to model organisms or putative regulatory regions; and (iii) minimum SNP subset selection by either a Brute Force Algorithm or a Greedy Partition Algorithm. The HapScope viewer displays genomic structure with haplotype information in an integrated environment, providing eight alternative views for assessing genetic and functional correlation. It has a user-friendly interface for: (i) haplotype block visualization; (ii) SNP subset selection; (iii) haplotype consolidation with subset SNP markers; (iv) incorporation of both experimentally determined haplotypes and computational results; and (v) data export for additional analysis. Comparison of haplotypes constructed by the statistical algorithms with those determined experimentally shows variation in haplotype prediction accuracies in genomic regions with different levels of nucleotide diversity. We have applied HapScope in analyzing haplotypes for candidate genes and genomic regions with extensive SNP and genotype data. We envision that the systematic approach of integrating functional genomic analysis with population haplotypes, supported by HapScope, will greatly facilitate current genetic disease research.     

Novel strategies to mine alcoholism-related haplotypes and genes by combining existing knowledge framework

High-throughout single nucleotide polymorphism detection technology and the existing knowledge provide strong support for mining the disease-related haplotypes and genes. In this study, first, we apply four kinds of haplotype identification methods (Confidence Intervals, Four Gamete Tests, Solid Spine of LD and fusing method of haplotype block) into high-throughout SNP genotype data to identify blocks, then use cluster analysis to verify the effectiveness of the four methods, and select the alcoholism-related SNP haplotypes through risk analysis. Second, we establish a mapping from haplotypes to alcoholism-related genes. Third, we inquire NCBI SNP and gene databases to locate the blocks and identify the candidate genes. In the end, we make gene function annotation by KEGG, Biocarta, and GO database. We find 159 haplotype blocks, which relate to the alcoholism most possibly on chromosome 1∼22, including 227 haplotypes, of which 102 SNP haplotypes may increase the risk of alcoholism. We get 121 alcoholism-related genes and verify their reliability by the functional annotation of biology. In a word, we not only can handle the SNP data easily, but also can locate the disease-related genes precisely by combining our novel strategies of mining alcoholism-related haplotypes and genes with existing knowledge framework. Supported by the National Natural Science Foundation of China (Grant Nos. 30570424, 60601010 and 30600367), the National High-Tech Research and Development Program of China, (Grant No.2007AA02Z329), the Key Science and Technology Program of Heilongjiang Province(Grant No.GB03C602-4), Natural Science Foundation of Heilongjiang Province (Grant No. F2008-02), Youth Science Foundation of Harbin Medical University (Grant No. 060045) and Science Foundation of Heilongjiang Province Education Department (Grant Nos. 11531113 and 1152hq28).     

A study of associations of the H-FABP genotypes with fat and meat production of pigs

Although the content of intramuscular fat (IMF) influences significantly meat quality, it can be estimated only after the slaughter of animals. Variants of the H-FABP gene were suggested as candidate genes influencing the variability of IMF. The effect of H-FABP - HinfI polymorphism on the content of IMF, backfat thickness, the weight and percentages of major meat parts and of the leg in carcass weight was studied in a group of 97 pigs (46 gilts and 51 barrows) of Large White and Landrace breeds using the test of fattening capacity and carcass value. In the set of experimental animals, the frequencies of genes were H = 0.75 +/- 0.03 and h = 0.25 +/- 0.03. Biometric analyses did not corroborate differences among different H-FABP - HinfI genotypes and all the traits under study. Only in genotypes HH and Hh the differences between least-square means of phenotypic IMF values under study were close to the limit of significance (P = 0.06).     

Polymorphism study of seven SNPs at ADH genes in 15 Chinese populations

It has been shown that the variants of alcohol dehydrogenase (ADH) genes exhibit great diversities among various populations and are associated with susceptibility to alcoholism. To investigate the distribution of SNPs at ADH genes in Chinese populations and the genetic relationship of these groups, we collected 467 individuals from 15 groups distributing widely from north to south in China and genotyped 7 SNPs at ADH genes respectively. The statistic analyses of allele frequencies, estimated haplotype frequencies, pairwise linkage disequilibrium, AMOVA (analysis of molecular variance), pairwise Fst', and cluster analysis indicated (1) that six of these seven SNPs showed great variations in the 15 Chinese populations, and three of them (RsaI, SspI, EcoRI), were confirmed to be informative SNPs. However, the causative SNP ADH1B Arg47His confirmed in case-control studies could not act as significant indicator to distinguish bibulous groups from non-bibulous groups in healthy individuals; (2) haplotypes constructed with ADH SNPs could be used as markers to discern different populations in China, and six-allele haplotype "221211" was the most common one defined in present study; (3) on the basis of SNPs analysis of ADH genes, the 15 populations were grouped into northern groups and southern groups. Moreover, the origin relationship among the populations was indicated according to the results of cluster analysis.     

A Genome-Wide Association Study Identifies Risk Loci to Equine Recurrent Uveitis in German Warmblood Horses

Equine recurrent uveitis (ERU) is a common eye disease affecting up to 3–15% of the horse population. A genome-wide association study (GWAS) using the Illumina equine SNP50 bead chip was performed to identify loci conferring risk to ERU. The sample included a total of 144 German warmblood horses. A GWAS showed a significant single nucleotide polymorphism (SNP) on horse chromosome (ECA) 20 at 49.3 Mb, with IL-17A and IL-17F being the closest genes. This locus explained a fraction of 23% of the phenotypic variance for ERU. A GWAS taking into account the severity of ERU, revealed a SNP on ECA18 nearby to the crystalline gene cluster CRYGA-CRYGF. For both genomic regions on ECA18 and 20, significantly associated haplotypes containing the genome-wide significant SNPs could be demonstrated. In conclusion, our results are indicative for a genetic component regulating the possible critical role of IL-17A and IL-17F in the pathogenesis of ERU. The associated SNP on ECA18 may be indicative for cataract formation in the course of ERU.