SHEsis, a powerful software platform for analyses of linkage disequilibrium, haplotype construction, and genetic association at polymorphism loci

In multiloci-based genetic association studies of complex diseases, a powerful and high efficient tool for analyses of linkage disequilibrium (LD) between markers, haplotype distributions and many chi-square/p values with a large number of samples has been sought for long. In order to achieve the goal of obtaining meaningful results directly from raw data,we developed a robust and user-friendly software platform with a series of tools for analysis in association study with high efficiency. The platform has been well evaluated by several sets of real data.     

An approach to incorporate linkage disequilibrium structure into genomic association analysis

In this study, we propose to use the principal component analysis (PCA) and regression model to incorporate linkage disequilibrium (LD) in genomic association data analysis. To accommodate LD in genomic data and reduce multiple testing, we suggest performing PCA and extracting the PCA score to capture the variation of genomic data, after which regression analysis is used to assess the association of the disease with the principal component score. An empirical analysis result shows that both genotype-basod correlation matrix and haplotype-based LD matrix can produce similar results for PCA. Principal component score seems to be more powerful in detecting genetic association because the principal component score is quantitatively measured and may be able to capture the effect of multiple loci.     

Three novel quantitative trait loci for skin thickness in swine identified by linkage and genome‐wide association studies

Skin is the largest organ in the pig body and plays a key role in protecting the body against pathogens and excessive water loss. Deciphering the genetic basis of swine skin thickness would enrich our knowledge about the skin. To identify the loci for porcine skin thickness, we first performed a genome scan with 194 microsatellite markers in a White Duroc × Erhualian F₂ intercross. We identified three genome‐wide significant QTL on pig chromosomes (SSC) 4, 7 and 15 using linkage analysis. The most significant QTL was found on SSC7 with a small confidence interval of ~5 cM, explaining 23.9 percent of phenotypic variance. Further, we conducted a genome‐wide association study (GWAS) using Illumina PorcineSNP60 Beadchips for the F₂ pedigree and a population of Chinese Sutai pigs. We confirmed significant QTL in the F₂ pedigree and replicated QTL on SSC15 in Chinese Sutai pigs. A meta‐analysis of GWASs on both populations detected a genomic region associated with skin thickness on SSC4. GWAS results were generally consistent with QTL mapping. Identical‐by‐descent analysis defined QTL on SSC7 in a 683‐kb region harboring an interesting candidate gene: HMGA1. On SSC15, the linkage disequilibrium analysis showed a haplotype block of 2.20 Mb that likely harbors the gene responsible for skin thickness. Our findings provide novel insights into the genetic basis of swine skin thickness, which would benefit further understanding of porcine skin function.     

Genome‐wide association studies for hematological traits in swine

Improving immune capacity may increase the profitability of animal production if it enables animals to better cope with infections. Hematological traits play pivotal roles in animal immune capacity and disease resistance. Thus far, few studies have been conducted using a high‐density swine SNP chip panel to unravel the genetic mechanism of the immune capability in domestic animals. In this study, using mixed model‐based single‐locus regression analyses, we carried out genome‐wide association studies, using the Porcine SNP60 BeadChip, for immune responses in piglets for 18 hematological traits (seven leukocyte traits, seven erythrocyte traits, and four platelet traits) after being immunized with classical swine fever vaccine. After adjusting for multiple testing based on permutations, 10, 24, and 77 chromosome‐wise significant SNPs were identified for the leukocyte traits, erythrocyte traits, and platelet traits respectively, of which 10 reached genome‐wise significance level. Among the 53 SNPs for mean platelet volume, 29 are located in a linkage disequilibrium block between 32.77 and 40.59 Mb on SSC6. Four genes of interest are located within the block, providing genetic evidence that this genomic segment may be considered a candidate region relevant to the platelet traits. Other candidate genes of interest for red blood cell, hemoglobin, and red blood cell volume distribution width also have been found near the significant SNPs. Our genome‐wide association study provides a list of significant SNPs and candidate genes that offer valuable information for future dissection of molecular mechanisms regulating hematological traits.     

Impact of polymorphic transposable elements on linkage disequilibrium along chromosomes

Recombination and selection drive the extent of linkage disequilibrium (LD) among loci and therefore affect the reshuffling of adaptive genetic variation. However, it is poorly known to what extent the enrichment of transposable elements (TEs) in recombinationally‐inert regions reflects their inefficient removal by purifying selection and whether the presence of polymorphic TEs can modify the local recombination rate. In this study, we investigate how TEs and recombination interact at fine scale along chromosomes and possibly support linked selection in natural populations. Whole‐genome sequencing data of 304 individuals from nearby alpine populations of Arabis alpina were used to show that the density of polymorphic TEs is specifically correlated with local LD along chromosomes. Consistent with TEs modifying recombination, the characterization of 28 such LD blocks of up to 5.5 Mb in length revealed strong evidence of selective sweeps at a few loci through either site frequency spectrum or haplotype structure. A majority of these blocks were enriched in genes related to ecologically relevant functions such as responses to cold, salt stress or photoperiodism. In particular, the S‐locus (i.e., supergene responsible for strict outcrossing) was identified in a LD block with high levels of polymorphic TEs and evidence of selection. Another such LD block was enriched in cold‐responding genes and presented evidence of adaptive loci related to photoperiodism and flowering being increasingly linked by polymorphic TEs. These results are consistent with the hypothesis that TEs modify recombination landscapes and thus interact with selection in driving blocks of linked adaptive loci in natural populations.     

Linkage disequilibrium and association studies in higher plants: Present status and future prospects

During the last two decades, DNA-based molecular markers have been extensively utilized for a variety of studies in both plant and animal systems. One of the major uses of these markers is the construction of genome-wide molecular maps and the genetic analysis of simple and complex traits. However, these studies are generally based on linkage analysis in mapping populations, thus placing serious limitations in using molecular markers for genetic analysis in a variety of plant systems. Therefore, alternative approaches have been suggested, and one of these approaches makes use of linkage disequilibrium (LD)-based association analysis. Although this approach of association analysis has already been used for studies on genetics of complex traits (including different diseases) in humans, its use in plants has just started. In the present review, we first define and distinguish between LD and association mapping, and then briefly describe various measures of LD and the two methods of its depiction. We then give a list of different factors that affect LD without discussing them, and also discuss the current issues of LD research in plants. Later, we also describe the various uses of LD in plant genomics research and summarize the present status of LD research in different plant genomes. In the end, we discuss briefly the future prospects of LD research in plants, and give a list of softwares that are useful in LD research, which is available as electronic supplementary material (ESM) Electronic supplementary material Electronic supplementary material is available for this article at and accessible for authorised users.     

Assessment of the functionality of genome‐wide canine SNP arrays and implications for canine disease association studies

Domestic dogs share a wide range of important disease conditions with humans, including cancers, diabetes and epilepsy. Many of these conditions have similar or identical underlying pathologies to their human counterparts and thus dogs represent physiologically relevant natural models of human disorders. Comparative genomic approaches whereby disease genes can be identified in dog diseases and then mapped onto the human genome are now recognized as a valid method and are increasing in popularity. The majority of dog breeds have been created over the past few hundred years and, as a consequence, the dog genome is characterized by extensive linkage disequilibrium (LD), extending usually from hundreds of kilobases to several megabases within a breed, rather than tens of kilobases observed in the human genome. Genome‐wide canine SNP arrays have been developed, and increasing success of using these arrays to map disease loci in dogs is emerging. No equivalent of the human HapMap currently exists for different canine breeds, and the LD structure for such breeds is far less understood than for humans. This study is a dedicated large‐scale assessment of the functionalities (LD and SNP tagging performance) of canine genome‐wide SNP arrays in multiple domestic dog breeds. We have used genotype data from 18 breeds as well as wolves and coyotes genotyped by the Illumina 22K canine SNP array and Affymetrix 50K canine SNP array. As expected, high tagging performance was observed with most of the breeds using both Illumina and Affymetrix arrays when multi‐marker tagging was applied. In contrast, however, large differences in population structure, LD coverage and pairwise tagging performance were found between breeds, suggesting that study designs should be carefully assessed for individual breeds before undertaking genome‐wide association studies (GWAS).     

Beyond pathways: genetic dissection of tocopherol content in maize kernels by combining linkage and association analyses

Although tocopherols play an important role in plants and animals, the genetic architecture of tocopherol content in maize kernels has remained largely unknown. In this study, linkage and association analyses were conducted to examine the genetic architecture of tocopherol content in maize kernels. Forty‐one unique quantitative trait loci (QTLs) were identified by linkage mapping in six populations of recombinant inbred lines (RILs). In addition, 32 significant loci were detected via genome‐wide association study (GWAS), 18 of which colocalized with the QTLs identified by linkage mapping. Fine mapping of a major QTL validated the accuracy of GWAS and QTL mapping results and suggested a role for nontocopherol pathway genes in the modulation of natural tocopherol variation. We provided genome‐wide evidence that genes involved in fatty acid metabolism, chlorophyll metabolism and chloroplast function may affect natural variation in tocopherols. These findings were confirmed through mutant analysis of a particular gene from the fatty acid pathway. In addition, the favourable alleles for many of the significant SNPs/QTLs represented rare alleles in natural populations. Together, our results revealed many novel genes that are potentially involved in the variation of tocopherol content in maize kernels. Pyramiding of the favourable alleles of the newly elucidated genes and the well‐known tocopherol pathway genes would greatly improve tocopherol content in maize.     

Genome wide linkage disequilibrium in Chinese asparagus bean (Vigna. unguiculata ssp. sesquipedialis) germplasm: implications for domestication history and genome wide association studies

Association mapping of important traits of crop plants relies on first understanding the extent and patterns of linkage disequilibrium (LD) in the particular germplasm being investigated. We characterize here the genetic diversity, population structure and genome wide LD patterns in a set of asparagus bean (Vigna. unguiculata ssp. sesquipedialis) germplasm from China. A diverse collection of 99 asparagus bean and normal cowpea accessions were genotyped with 1127 expressed sequence tag-derived single nucleotide polymorphism markers (SNPs). The proportion of polymorphic SNPs across the collection was relatively low (39%), with an average number of SNPs per locus of 1.33. Bayesian population structure analysis indicated two subdivisions within the collection sampled that generally represented the 'standard vegetable' type (subgroup SV) and the 'non-standard vegetable' type (subgroup NSV), respectively. Level of LD (r(2)) was higher and extent of LD persisted longer in subgroup SV than in subgroup NSV, whereas LD decayed rapidly (0-2 cM) in both subgroups. LD decay distance varied among chromosomes, with the longest (≈ 5 cM) five times longer than the shortest (≈ 1 cM). Partitioning of LD variance into within- and between-subgroup components coupled with comparative LD decay analysis suggested that linkage group 5, 7 and 10 may have undergone the most intensive epistatic selection toward traits favorable for vegetable use. This work provides a first population genetic insight into domestication history of asparagus bean and demonstrates the feasibility of mapping complex traits by genome wide association study in asparagus bean using a currently available cowpea SNPs marker platform.     

Gene mapping by linkage and association analysis

Genetic analysis is used to map genes, including disease loci, to positions within the human genome. Linkage analysis depends on the co-segregation of a gene (locus) and a phenotype through a pedigree, while association analysis, or linkage disequilibrium mapping, depends on measuring deviation from the random occurrence of alleles in a haplotype in unrelated individuals or nuclear families. Complex computer programs may be used in both forms of analysis. In recent years most interest has focused on identifying genes involved in common, multifactorial diseases. Here I review some current and developing techniques of genetic analysis and give references to where further information can be obtained.     

The extent of linkage disequilibrium in a large cattle population of western Africa and its consequences for association studies

Several previous studies concluded that linkage disequilibrium (LD) in livestock populations from developed countries originated from the impact of strong selection. Here, we assessed the extent of LD in a cattle population from western Africa that was bred in an extensive farming system. The analyses were performed on 363 individuals in a Bos indicus x Bos taurus population using 42 microsatellite markers on BTA04, BTA07 and BTA13. A high level of expected heterozygosity (0.71), a high mean number of alleles per locus (9.7) and a mild shift in Hardy-Weinberg equilibrium were found. Linkage disequilibrium extended over shorter distances than what has been observed in cattle from developed countries. Effective population size was assessed using two methods; both methods produced large values: 1388 when considering heterozygosity (assuming a mutation rate of 10(-3)) and 2344 when considering LD on whole linkage groups (assuming a constant population size over generations). However, analysing the decay of LD as a function of marker spacing indicated a decreasing trend in effective population size over generations. This decrease could be explained by increasing selective pressure and/or by an admixture process. Finally, LD extended over small distances, which suggested that whole-genome scans will require a large number of markers. However, association studies using such populations will be effective.     

A genome‐wide association study suggests new candidate genes for milk production traits in Chinese Holstein cattle

A genome‐wide association study (GWAS) was conducted on 15 milk production traits in Chinese Holstein. The experimental population consisted of 445 cattle, each genotyped by the GGP (GeneSeek genomic profiling)‐BovineLD V3 SNP chip, which had 26 151 public SNPs in its manifest file. After data cleaning, 20 326 SNPs were retained for the GWAS. The phenotypes were estimated breeding values of traits, provided by a public dairy herd improvement program center that had been collected once a month for 3 years. Two statistical models, a fixed‐effect linear regression model and a mixed‐effect linear model, were used to estimate the association effects of SNPs on each of the phenotypes. Genome‐wide significant and suggestive thresholds were set at 2.46E‐06 and 4.95E‐05 respectively. The two statistical models concurrently identified two genome‐wide significant (P < 0.05) SNPs on milk production traits in this Chinese Holstein population. The positional candidate genes, which were the ones closest to these two identified SNPs, were EEF2K (eukaryotic elongation factor 2 kinase) and KLHL1 (kelch like family member 1). These two genes could serve as new candidate genes for milk yield and lactation persistence, yet their roles need to be verified in further function studies.     

Recombination and linkage disequilibrium among mitochondrial genes in structured populations of the gynodioecious plant Silene vulgaris

The impact of intergenic recombination on the population genetics of plant mitochondrial genomes is unknown. In an effort to study this in the gynodioecious plant Silene vulgaris three-locus PCR/RFLP genotypes (based on the mitochondrial genes atpA, cox1, and cob) were determined for 239 individuals collected from 20 North American populations. Seventeen three-locus PCR/RFLP genotypes were found. Recombination was indicated by observation of each of the four two-locus genotypes possible when the two most common alleles are considered for each of two loci. Based on these common alleles the absolute values of standardized linkage disequilibrium |D'| between pairs of loci range from 0.17 to 0.78. This indicates modest disequilibrium, rather than the maximum value expected in the absence of recombination |D'=1|, or the linkage equilibrium expected if recombination is pervasive (D'=0). Values of D' did not depend on which pair of loci contributed alleles to the analysis. The direction of D' obtained for the common atpA and cox1 alleles was comparable in sign and magnitude to that obtained by examining similar information obtained in a prior study of European samples. All three loci indicated a high degree of population structure (average FST=0.63), which would limit the within-population genetic diversity required for intergenic recombination to create novel genotypes, if most mating is local. Thus, population structure acts as a constraint on the approach to linkage equilibrium.     

Genome‐wide association QTL mapping for teat number in a purebred population of Duroc pigs

Because of increasing litter size in Western pig breeds, additional teats are desirable to increase the capacity for nursing offspring. We applied genome‐wide SNP markers to detect QTL regions that affect teat number in a Duroc population. We phenotyped 1024 animals for total teat number. A total of 36 588 SNPs on autosomes were used in the analysis. The estimated heritability for teat number was 0.34 ± 0.05 on the basis of a genomic relationship matrix constructed from all SNP markers. Using a BayesC method, we identified a total of 18 QTL regions that affected teat number in Duroc pigs; 9 of the 18 regions were newly detected.     

A program to test linkage disequilibrium between loci in autotetraploid species

No software currently implements a test of linkage disequilibrium in autotetraploid species. We propose a program, LD4X that performs a Fisher's exact test between pairs of alleles at two loci. All combinations of alleles from two loci are treated in turn. If two alleles of a pair of loci have a nonrandom distribution, the markers are considered as linked. The program was tested on a set of microsatellite markers in synthetic alfalfa populations.     

Genome‐wide association and epistasis studies unravel the genetic architecture of sudden death syndrome resistance in soybean

Soybean [Glycine max (L.) Merr.] is an economically important crop that is grown worldwide. Sudden death syndrome (SDS), caused by Fusarium virguliforme, is one of the top yield‐limiting diseases in soybean. However, the genetic basis of SDS resistance, especially with respect to epistatic interactions, is still unclear. To better understand the genetic architecture of soybean SDS resistance, genome‐wide association and epistasis studies were performed using a population of 214 germplasm accessions and 31 914 SNPs from the SoySNP50K Illumina Infinium BeadChip. Twelve loci and 12 SNP–SNP interactions associated with SDS resistance were identified at various time points after inoculation. These additive and epistatic loci together explained 24–52% of the phenotypic variance. Disease‐resistant, pathogenesis‐related and chitin‐ and wound‐responsive genes were identified in the proximity of peak SNPs, including stress‐induced receptor‐like kinase gene 1 (SIK1), which is pinpointed by a trait‐associated SNP and encodes a leucine‐rich repeat‐containing protein. We report that the proportion of phenotypic variance explained by identified loci may be considerably improved by taking epistatic effects into account. This study shows the necessity of considering epistatic effects in soybean SDS resistance breeding using marker‐assisted and genomic selection approaches. Based on our findings, we propose a model for soybean root defense against the SDS pathogen. Our results facilitate identification of the molecular mechanism underlying SDS resistance in soybean, and provide a genetic basis for improvement of soybean SDS resistance through breeding strategies based on additive and epistatic effects.     

Analysis of the genetic architecture of maize kernel size traits by combined linkage and association mapping

Kernel size‐related traits are the most direct traits correlating with grain yield. The genetic basis of three kernel traits of maize, kernel length (KL), kernel width (KW) and kernel thickness (KT), was investigated in an association panel and a biparental population. A total of 21 single nucleotide polymorphisms (SNPs) were detected to be most significantly (P < 2.25 × 10^?6) associated with these three traits in the association panel under four environments. Furthermore, 50 quantitative trait loci (QTL) controlling these traits were detected in seven environments in the intermated B73 × Mo17 (IBM) Syn10 doubled haploid (DH) population, of which eight were repetitively identified in at least three environments. Combining the two mapping populations revealed that 56 SNPs (P < 1 × 10^?3) fell within 18 of the QTL confidence intervals. According to the top significant SNPs, stable‐effect SNPs and the co‐localized SNPs by association analysis and linkage mapping, a total of 73 candidate genes were identified, regulating seed development. Additionally, seven miRNAs were found to situate within the linkage disequilibrium (LD) regions of the co‐localized SNPs, of which zma‐miR164e was demonstrated to cleave the mRNAs of Arabidopsis CUC1, CUC2 and NAC6 in vitro. Overexpression of zma‐miR164e resulted in the down‐regulation of these genes above and the failure of seed formation in Arabidopsis pods, with the increased branch number. These findings provide insights into the mechanism of seed development and the improvement of molecular marker‐assisted selection (MAS) for high‐yield breeding in maize.     

Genome-wide linkage disequilibrium in two Japanese beef cattle breeds

There is little knowledge about the degree of linkage disequilibrium (LD) in beef cattle. This study aims to perform a genome-wide search for LD in Japanese Black and Japanese Brown beef cattle and to compare the level of LD between these two breeds. Parameter D' (the LD coefficient) was used as a measure of LD, and LD was tested for significance of allelic associations between syntenic and between non-syntenic marker pairs. Effects of breed, chromosome, genetic map distance and their interactions with D' were tested based on least squares analyses. Both breeds showed high levels of LD, which ranged over several tens of cM and declined as the marker distance increased for syntenic marker pairs. A rapid decline of the D' value was observed between markers that were spaced 5 and 20 cM apart. LD was significant in most cases for marker pairs <40 cM apart but was not significant between non-syntenic loci. The pattern of LD found in these two breeds was similar to that previously published for dairy cattle. The D' value between breeds was not significantly different (P > 0.05), but the interaction between breed and chromosome was highly significant (P < 0.001). Genetic selection seems to have caused the heterogeneity of the D' values among chromosomes within breed. These results indicate that LD mapping is a useful tool for fine-mapping quantitative trait loci of economically important traits in Japanese beef cattle.     

Current opportunities and challenges: genome-wide association studies on pigmentation and skin cancer

Genome-wide association studies (GWAS) have become a widely used approach for genetic association studies of various human traits. A few GWAS have been conducted with the goal of identifying novel loci for pigmentation traits, melanoma, and non-melanoma skin cancer. Nevertheless, the phenotype variation explained by the genetic markers identified so far is limited. In this review, we discuss the GWAS study design and its application in pigmentation and skin cancer research. Furthermore, we summarize recent developments in post-GWAS activities such as meta-analysis, pathway analysis, and risk prediction.