Meta‐analysis of genome‐wide association from genomic prediction models

Genome‐wide association (GWA) studies based on GBLUP models are a common practice in animal breeding. However, effect sizes of GWA tests are small, requiring larger sample sizes to enhance power of detection of rare variants. Because of difficulties in increasing sample size in animal populations, one alternative is to implement a meta‐analysis (MA), combining information and results from independent GWA studies. Although this methodology has been used widely in human genetics, implementation in animal breeding has been limited. Thus, we present methods to implement a MA of GWA, describing the proper approach to compute weights derived from multiple genomic evaluations based on animal‐centric GBLUP models. Application to real datasets shows that MA increases power of detection of associations in comparison with population‐level GWA, allowing for population structure and heterogeneity of variance components across populations to be accounted for. Another advantage of MA is that it does not require access to genotype data that is required for a joint analysis. Scripts related to the implementation of this approach, which consider the strength of association as well as the sign, are distributed and thus account for heterogeneity in association phase between QTL and SNPs. Thus, MA of GWA is an attractive alternative to summarizing results from multiple genomic studies, avoiding restrictions with genotype data sharing, definition of fixed effects and different scales of measurement of evaluated traits.     

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.     

A genome‐wide association analysis for carcass traits in a commercial Duroc pig population

We performed a genome‐wide association study to map the genetic determinants of carcass traits in 350 Duroc pigs typed with the Porcine SNP60 BeadChip. Association analyses were carried out using the gemma software. The proportion of phenotypic variance explained by the SNPs ranged between negligible to moderate (= 0.01–0.30) depending on the trait under consideration. At the genome‐wide level, we detected one significant association between backfat thickness between the 3^rd and 4^th ribs and six SNPs mapping to SSC12 (37–40 Mb). We also identified several chromosome‐wide significant associations for ham weight (SSC11: 51–53 Mb, three SNPs; 67–68 Mb, two SNPs), carcass weight (SSC11: 66–68 Mb, two SNPs), backfat thickness between the 3^rd and 4^th ribs (SSC12: 21 Mb, one SNP; 33–40 Mb, 17 SNPs; 51–58 Mb, two SNPs), backfat thickness in the last rib (SSC12: 37 Mb, one SNP; 40–41 Mb, nine SNPs) and lean meat content (SSC13: 34 Mb, three SNPs and SSC16: 45.1 Mb, one SNP; 62–63 Mb, 10 SNPs; 71–75 Mb, nine SNPs). The ham weight trait‐associated region on SSC11 contains two genes (UCHL3 and LMO7) related to muscle development. In addition, the ACACA gene, which encodes an enzyme for the catalysis of fatty acid synthesis, maps to the SSC12 (37–41 Mb) region harbouring trait‐associated regions for backfat thickness traits. Sequencing of these candidate genes may help to uncover the causal mutations responsible for the associations found in the present study.     

Genetic‐based dissection of arsenic accumulation in maize using a genome‐wide association analysis method

Understanding the mechanism of arsenic (As) accumulation in plants is important in reducing As's toxicity to plants and its potential risks to human health. Here, we performed a genome‐wide association study to dissect the genetic basis of the As contents of different maize tissues in Xixian, which was irrigated with As‐rich surface water, and Changge using an association population consisting of 230 representative maize inbred lines. Phenotypic data revealed a wide normal distribution and high repeatability for the As contents in maize tissues. The As concentrations in maize tissues followed the same trend in the two locations: kernels < axes < stems < bracts < leaves. In total, 15, 16 and 15 non‐redundant quantitative trait loci (QTL s) associated with As concentrations were identified (P  ≤ 2.04 × 10^?6) in five tissues from Xixian, Changge, and the combination of the locations, respectively, explaining 9.70%–24.65% of the phenotypic variation for each QTL , on average. Additionally, four QTL s [involving 15 single nucleotide polymorphisms (SNP s)] were detected in the single and the combined locations, indicating that these loci/SNP s might be stable across different environments. The candidate genes associated with these four loci were predicted. In addition, four non‐redundant QTL s (6 SNP s), including a QTL that was detected in multiple locations according to the genome‐wide association study, were found to co‐localize with four previously reported QTL intervals. These results are valuable to understand the genetic architecture of As mechanism in maize and facilitate the genetic improvement of varieties without As toxicity.     

Pathway analysis of genome‐wide association datasets of personality traits

Although several genome‐wide association (GWA) studies of human personality have been recently published, genetic variants that are highly associated with certain personality traits remain unknown, due to difficulty reproducing results. To further investigate these genetic variants, we assessed biological pathways using GWA datasets. Pathway analysis using GWA data was performed on 1089 Korean women whose personality traits were measured with the Revised NEO Personality Inventory for the 5‐factor model of personality. A total of 1042 pathways containing 8297 genes were included in our study. Of these, 14 pathways were highly enriched with association signals that were validated in 1490 independent samples. These pathways include association of: Neuroticism with axon guidance [L1 cell adhesion molecule (L1CAM) interactions]; Extraversion with neuronal system and voltage‐gated potassium channels; Agreeableness with L1CAM interaction, neurotransmitter receptor binding and downstream transmission in postsynaptic cells; and Conscientiousness with the interferon‐gamma and platelet‐derived growth factor receptor beta polypeptide pathways. Several genes that contribute to top‐ranked pathways in this study were previously identified in GWA studies or by pathway analysis in schizophrenia or other neuropsychiatric disorders. Here we report the first pathway analysis of all five personality traits. Importantly, our analysis identified novel pathways that contribute to understanding the etiology of personality traits.     

A genome‐wide association study using selective DNA pooling identifies candidate markers for fertility in Holstein cattle

The decline in the reproductive efficiency of dairy cows, especially those with high producing potential, has become a challenging problem. In this study, a selective DNA pooling approach was applied to a cow population whose oocytes were fertilized and cultured to obtain phenotypic records of fertilization rate and blastocyst rate. Using a stringent 5% genome‐wide significance level, 22 and five single nucleotide polymorphisms (SNPs) were found to be associated with fertilization rate and blastocyst rate, respectively. SNPs that showed significant association in selective DNA pooling were further evaluated by individual genotyping. Interestingly, the majority of the SNP associations were confirmed by individual genotyping, testifying to the effectiveness of selective DNA pooling using a high‐density SNP genotyping array. This study is the first application of the selective DNA pooling approach using the BovineSNP50 array in cattle.     

A genome‐wide association study of copy number variations with umbilical hernia in swine

Umbilical hernia (UH) is one of the most common congenital defects in pigs, leading to considerable economic loss and serious animal welfare problems. To test whether copy number variations (CNVs) contribute to pig UH, we performed a case–control genome‐wide CNV association study on 905 pigs from the Duroc, Landrace and Yorkshire breeds using the Porcine SNP60 BeadChip and penncnv algorithm. We first constructed a genomic map comprising 6193 CNVs that pertain to 737 CNV regions. Then, we identified eight CNVs significantly associated with the risk for UH in the three pig breeds. Six of seven significantly associated CNVs were validated using quantitative real‐time PCR. Notably, a rare CNV (CNV14:13030843–13059455) encompassing the NUGGC gene was strongly associated with UH (permutation‐corrected P = 0.0015) in Duroc pigs. This CNV occurred exclusively in seven Duroc UH‐affected individuals. SNPs surrounding the CNV did not show association signals, indicating that rare CNVs may play an important role in complex pig diseases such as UH. The NUGGC gene has been implicated in human omphalocele and inguinal hernia. Our finding supports that CNVs, including the NUGGC CNV, contribute to the pathogenesis of pig UH.     

Population structure and genetic basis of the agronomic traits of upland cotton in China revealed by a genome‐wide association study using high‐density SNPs

Gossypium hirsutum L. represents the largest source of textile fibre, and China is one of the largest cotton‐producing and cotton‐consuming countries in the world. To investigate the genetic architecture of the agronomic traits of upland cotton in China, a diverse and nationwide population containing 503 G. hirsutum accessions was collected for a genome‐wide association study (GWAS) on 16 agronomic traits. The accessions were planted in four places from 2012 to 2013 for phenotyping. The CottonSNP63K array and a published high‐density map based on this array were used for genotyping. The 503 G. hirsutum accessions were divided into three subpopulations based on 11 975 quantified polymorphic single‐nucleotide polymorphisms (SNPs). By comparing the genetic structure and phenotypic variation among three genetic subpopulations, seven geographic distributions and four breeding periods, we found that geographic distribution and breeding period were not the determinants of genetic structure. In addition, no obvious phenotypic differentiations were found among the three subpopulations, even though they had different genetic backgrounds. A total of 324 SNPs and 160 candidate quantitative trait loci (QTL) regions were identified as significantly associated with the 16 agronomic traits. A network was established for multieffects in QTLs and interassociations among traits. Thirty‐eight associated regions had pleiotropic effects controlling more than one trait. One candidate gene, Gh_D08G2376, was speculated to control the lint percentage (LP). This GWAS is the first report using high‐resolution SNPs in upland cotton in China to comprehensively investigate agronomic traits, and it provides a fundamental resource for cotton genetic research and breeding.     

A genome‐wide association study identifies a region strongly associated with symmetrical onychomadesis on chromosome 12 in dogs

Symmetrical onychomadesis causes periodic loss of claws in otherwise healthy dogs. Genome‐wide association analysis in 225 Gordon Setters identified a single region associated with symmetrical onychomadesis on chromosome 12 (spanning about 3.3 mb). A meta‐analysis including also English Setters indicated that this genomic region predisposes for symmetrical onychomadesis in English Setters as well. The associated region spans most of the major histocompatibility complex and nearly 1 Mb downstream. Like many other autoimmune diseases, associations of symmetrical onychomadesis with DLA class II alleles have been reported. In this study, no associated markers were revealed within any of the DLA‐DRB1, ‐DQA1 or ‐DQB1 genes, and the odds for symmetrical onychomadesis in the Gordon Setters were much higher, carrying significant single nucleotide polymorphisms compared to the odds of any of the recorded DLA‐DRB1/DQA1/DQB1 haplotypes. We noticed that some of the associated DLA haplotypes were different between the English Setters and the Gordon Setters. Interestingly, associated SNP chip markers showed a more consistent pattern of allelic variants related to cases or controls regardless of breed. In conclusion, the associated genetic markers identified in this study hold the potential to aid in selection of breeding animals to reduce the frequency of symmetrical onychomadesis in the dog.     

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.     

A genome‐wide association study reveals candidate genes for the supernumerary nipple phenotype in sheep (Ovis aries)

Genome‐wide association studies (GWASs) have been widely applied in livestock to identify genes associated with traits of economic interest. Here, we conducted the first GWAS of the supernumerary nipple phenotype in Wadi sheep, a native Chinese sheep breed, based on Ovine Infinium HD SNP BeadChip genotypes in a total of 144 ewes (75 cases with four teats, including two normal and two supernumerary teats, and 69 control cases with two teats). We detected 63 significant SNPs at the chromosome‐wise threshold. Additionally, one candidate region (chr1: 170.723–170.734 Mb) was identified by haplotype‐based association tests, with one SNP (rs413490006) surrounding functional genes BBX and CD47 on chromosome 1 being commonly identified as significant by the two mentioned analyses. Moreover, Gene Ontology enrichment for the significant SNPs identified by the GWAS analysis was functionally clustered into the categories of receptor activity and synaptic membrane. In addition, pathway mapping revealed four promising pathways (Wnt, oxytocin, MAPK and axon guidance) involved in the development of the supernumerary nipple phenotype. Our results provide novel and important insights into the genetic mechanisms underlying the phenotype of supernumerary nipples in mammals, including humans. These findings may be useful for future breeding and genetics in sheep and other livestock.     

A genome‐wide association study of 23 agronomic traits in Chinese wheat landraces

Uncovering the genetic basis of agronomic traits in wheat landraces is important for ensuring global food security via the development of improved varieties. Here, 723 wheat landraces from 10 Chinese agro‐ecological zones were evaluated for 23 agronomic traits in six environments. All accessions could be clustered into five subgroups based on phenotypic data via discriminant function analysis, which was highly consistent with genotypic classification. A genome‐wide association study was conducted for these traits using 52 303 DArT‐seq markers to identify marker‐trait associations and candidate genes. Using both the general linear model and the mixed linear model, 149 significant markers were identified for 21 agronomic traits based on best linear unbiased prediction values. Considering the linkage disequilibrium decay distance in this study, significant markers within 10 cM were combined as a quantitative trait locus (QTL), with a total of 29 QTL identified for 15 traits. Of these, five QTL for heading date, flag leaf width, peduncle length, and thousand kernel weight had been reported previously. Twenty‐five candidate genes associated with significant markers were identified. These included the known vernalization genes VRN‐B1 and vrn‐B3 and the photoperiod response genes Ppd and PRR. Overall, this study should be helpful in elucidating the underlying genetic mechanisms of complex agronomic traits and performing marker‐assisted selection in wheat.     

Crop genome‐wide association study: a harvest of biological relevance

With the advent of rapid genotyping and next‐generation sequencing technologies, genome‐wide association study (GWAS) has become a routine strategy for decoding genotype–phenotype associations in many species. More than 1000 such studies over the last decade have revealed substantial genotype–phenotype associations in crops and provided unparalleled opportunities to probe functional genomics. Beyond the many ‘hits’ obtained, this review summarizes recent efforts to increase our understanding of the genetic architecture of complex traits by focusing on non‐main effects including epistasis, pleiotropy, and phenotypic plasticity. We also discuss how these achievements and the remaining gaps in our knowledge will guide future studies. Synthetic association is highlighted as leading to false causality, which is prevalent but largely underestimated. Furthermore, validation evidence is appealing for future GWAS, especially in the context of emerging genome‐editing technologies.     

Genome‐wide association study for pigmentation traits in Chinese Holstein population

With the Illumina BovineSNP50K BeadChip, we performed a genome‐wide association study (GWAS) for two pigmentation traits in a Chinese Holstein population: proportion of black (PB) and teat colour (TC). A case–control design was used. Cases were the cows with PB <0.30 (n = 129) and TC <2 points (n = 140); controls were those with PB >0.90 (n = 58) and TC >4 points (n = 281). The RM test of roadtrips (version 1.2) was applied to detect SNPs for the two traits with 42 883 and 42 741 SNPs respectively. A total of nine and 12 genome‐wide significant (P < 0.05) SNPs associated with PB and TC respectively were identified. Of these, two SNPs for PB were located within the KIT and IGFBP7 genes, and the other four SNPs were 23~212 kb away from the PDGFRA gene on BTA6; nine SNPs associated with TC were located within or 21~78.8 kb away from known genes on chromosomes 4, 11, 22, 23 and 24. By combing through our GWAS results and the biological functions of the genes, we suggest that the KIT, IGFBP7, PDGFRA, MITF, ING3 and WNT16 genes are promising candidates for PB and TC in Holstein cattle, providing a basis for further investigation on the genetic mechanism of pigmentation formation.     

Metabolome‐genome‐wide association study dissects genetic architecture for generating natural variation in rice secondary metabolism

Plants produce structurally diverse secondary (specialized) metabolites to increase their fitness for survival under adverse environments. Several bioactive compounds for new drugs have been identified through screening of plant extracts. In this study, genome‐wide association studies (GWAS) were conducted to investigate the genetic architecture behind the natural variation of rice secondary metabolites. GWAS using the metabolome data of 175 rice accessions successfully identified 323 associations among 143 single nucleotide polymorphisms (SNPs) and 89 metabolites. The data analysis highlighted that levels of many metabolites are tightly associated with a small number of strong quantitative trait loci (QTLs). The tight association may be a mechanism generating strains with distinct metabolic composition through the crossing of two different strains. The results indicate that one plant species produces more diverse phytochemicals than previously expected, and plants still contain many useful compounds for human applications.     

Non‐replication of genome‐wide‐based associations of efficient food conversion in dairy cows

Animal growth relative to food energy input is of key importance to agricultural production. Several recent studies highlighted genetic markers associated with food conversion efficiency in beef cattle, and there is now a requirement to validate these associations in additional populations and to assess their potential utility for selecting animals with enhanced food‐use efficiency. The current analysis tested a population of dairy cattle using 138 DNA markers previously associated with food intake and growth in a whole‐genome association analysis of beef animals. Although seven markers showed point‐wise significance at P < 0.05, none of the single‐nucleotide polymorphisms tested were significantly associated with food conversion efficiency after correction for multiple testing. These data do not support the involvement of this subset of previously implicated markers in the food conversion efficiency of the physiologically distinct New Zealand Holstein‐Friesian dairy breed.     

Metabolite‐based genome‐wide association study enables dissection of the flavonoid decoration pathway of wheat kernels

The marriage of metabolomic approaches with genetic design has proven a powerful tool in dissecting diversity in the metabolome and has additionally enhanced our understanding of complex traits. That said, such studies have rarely been carried out in wheat. In this study, we detected 805 metabolites from wheat kernels and profiled their relative contents among 182 wheat accessions, conducting a metabolite‐based genome‐wide association study (mGWAS) utilizing 14 646 previously described polymorphic SNP markers. A total of 1098 mGWAS associations were detected with large effects, within which 26 candidate genes were tentatively designated for 42 loci. Enzymatic assay of two candidates indicated they could catalyse glucosylation and subsequent malonylation of various flavonoids and thereby the major flavonoid decoration pathway of wheat kernel was dissected. Moreover, numerous high‐confidence genes associated with metabolite contents have been provided, as well as more subdivided metabolite networks which are yet to be explored within our data. These combined efforts presented the first step towards realizing metabolomics‐associated breeding of wheat.