Genome‐wide association analysis reveals novel loci for hypoxia adaptability in Tibetan chicken

The Tibetan chicken (TBC), an indigenous chicken breed of the Tibetan Plateau, has adapted to its hypoxic, high‐altitude environment over hundreds of years. The objective of this study was to identify the polymorphisms and genes associated with adaptation to hypoxia in this chicken breed. In the present study, samples were collected during days 18–21 of the incubation period from both surviving chicks and dead embryos, all of which were hatched under hypoxic conditions. A genome‐wide association study was conducted using the Illumina iSelect 60K SNP array with a case–control design, in which the case group consisted of the dead chicken embryos (n = 54) and controls were the surviving chicks (n = 82). Four significant SNPs were detected at the genome‐wide level (P < 0.05), and the results indicated that fork head box G1 (FOXG1) was the main candidate gene. The lead SNP NC_006092.4:g.33368893T>C was confirmed with a polymerase chain reaction‐restriction fragment length polymorphism analysis of 122 cases and 212 controls. A chi‐square test showed a significant association between NC_006092.4:g.33368893T>C and hatchability under hypoxic conditions (P < 0.01). Our results revealed novel polymorphisms and a candidate gene associated with hypoxic adaptation, facilitating further study in this field.     

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.     

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.     

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 search for quantitative trait loci affecting the cortical surface area and thickness of Heschl's gyrus

Heschl's gyrus (HG) is a core region of the auditory cortex whose morphology is highly variable across individuals. This variability has been linked to sound perception ability in both speech and music domains. Previous studies show that variations in morphological features of HG, such as cortical surface area and thickness, are heritable. To identify genetic variants that affect HG morphology, we conducted a genome‐wide association scan (GWAS) meta‐analysis in 3054 healthy individuals using HG surface area and thickness as quantitative traits. None of the single nucleotide polymorphisms (SNPs) showed association P values that would survive correction for multiple testing over the genome. The most significant association was found between right HG area and SNP rs72932726 close to gene DCBLD2 (3q12.1; P = 2.77 × 10^?7). This SNP was also associated with other regions involved in speech processing. The SNP rs333332 within gene KALRN (3q21.2; P = 2.27 × 10^?6) and rs143000161 near gene COBLL1 (2q24.3; P = 2.40 × 10^?6) were associated with the area and thickness of left HG, respectively. Both genes are involved in the development of the nervous system. The SNP rs7062395 close to the X‐linked deafness gene POU3F4 was associated with right HG thickness (Xq21.1; P = 2.38 × 10^?6). This is the first molecular genetic analysis of variability in HG morphology.     

Genome‐wide association study reveals candidate genes associated with body measurement traits in Chinese Wagyu beef cattle

We performed a genome‐wide association study to identify candidate genes for body measurement traits in 463 Wagyu beef cattle typed with the Illumina Bovine HD 770K SNP array. At the genome‐wide level, we detected 18, five and one SNPs associated with hip height, body height and body length respectively. In total, these SNPs are within or near 11 genes, six of which (PENK, XKR4, IMPAD1, PLAG1, CCND2 and SNTG1) have been reported previously and five of which (CSMD3, LAP3, SYN3, FAM19A5 and TIMP3) are novel candidate genes that we found to be associated with body measurement traits. Further exploration of these candidate genes will facilitate genetic improvement in Chinese Wagyu beef cattle.     

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.     

Susceptibility loci for metabolic syndrome and metabolic components identified in Han Chinese: a multi‐stage genome‐wide association study

Metabolic syndrome (MetS), a cluster of metabolic disturbances that increase the risk for cardiovascular disease and diabetes, was because of genetic susceptibility and environmental risk factors. To identify the genetic variants associated with MetS and metabolic components, we conducted a genome‐wide association study followed by replications in totally 12,720 participants from the north, north‐eastern and eastern China. In combined analyses, independent of the top known signal at rs651821 on APOA5, we newly identified a secondary triglyceride‐associated signal at rs180326 on BUD13 (P_combined = 2.4 × 10⁻⁸). Notably, by an integrated analysis of the genotypes and the serum levels of APOA5, BUD13 and triglyceride, we observed that BUD13 was another potential mediator, besides APOA5, of the association between rs651821 and serum triglyceride. rs671 (ALDH2), an east Asian‐specific common variant, was found to be associated with MetS (P_combined = 9.7 × 10⁻²²) in Han Chinese. The effects of rs671 on metabolic components were more prominent in drinkers than in non‐drinkers. The replicated loci provided information on the genetic basis and mechanisms of MetS and metabolic components in Han Chinese.     

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).     

Genome‐wide association analysis identifies potential regulatory genes for eumelanin pigmentation in chicken plumage

Plumage color in chicken is determined by the proportion of eumelanin and pheomelanin pigmentation. As the main ingredient in plumage melanin, eumelanin plays a key role in the dark black, brown and grey coloration. However, very few studies have been performed to identify the related genes and mutations on a genome‐wide scale. Herein, a resource family consisting of one backcross population and two F2 cross populations between a black roster and Yukou Brown I parent stockbreed was constructed for identification of genes related to eumelanin pigmentation. Chickens with eumelanin in their plumage were classified as the case group, and the rest were considered the control group. A genome‐wide association study of this phenotype and genotypes using Affymetrix 600K HD SNP arrays in this F2 family revealed 13 significantly associated SNPs and in 10 separate genes on chromosomes 1, 2, 3 and 5. Based on previous studies in model species, we inferred that genes, including NUAK family kinase 1 (NUAK1) and sonic hedgehog (SHH), may play roles in the development of neural crest cells or melanoblasts during the embryonic period, which may also affect the eumelanin pigmentation. Our results facilitate the understanding of the genetic basis of eumelanin pigmentation in chicken plumage.     

A genome‐wide association study for equine recurrent airway obstruction in European Warmblood horses reveals a suggestive new quantitative trait locus on chromosome 13

Recurrent airway obstruction (RAO), also known as heaves, is an asthma‐like respiratory disease. Its development is strongly influenced by environmental risk factors such as sensitization and exposure to moldy hay, straw bedding and stabling indoors. A hereditary component has been documented in previous studies; however, so far no causative genetic variant that influences the risk of developing RAO has been identified. In this study, we revised an existing dataset and selected 384 horses for genotyping on the Affymetrix high‐density equine SNP array. We performed an allelic case–control genome‐wide association study, which revealed a suggestively significant association on equine chromosome 13 at 32 843 309 bp. This SNP is located in the protein‐coding gene TXNDC11, which is possibly involved in the folding process of the multiprotein complexes DUOX1 and DUOX2. In humans, these proteins are known to take part in regulating the production of H₂O₂ in the respiratory tract epithelium as well as in MUC5AC mucin expression. Therefore, TXNDC11 may be considered a functional candidate gene, and further research is needed to explore its potential role in RAO‐affected horses.     

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 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.     

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.     

A genome‐wide survey reveals abundant rice blast R genes in resistant cultivars

Plant resistance genes (R genes) harbor tremendous allelic diversity, constituting a robust immune system effective against microbial pathogens. Nevertheless, few functional R genes have been identified for even the best‐studied pathosystems. Does this limited repertoire reflect specificity, with most R genes having been defeated by former pests, or do plants harbor a rich diversity of functional R genes, the composite behavior of which is yet to be characterized? Here, we survey 332 NBS‐LRR genes cloned from five resistant Oryza sativa (rice) cultivars for their ability to confer recognition of 12 rice blast isolates when transformed into susceptible cultivars. Our survey reveals that 48.5% of the 132 NBS‐LRR loci tested contain functional rice blast R genes, with most R genes deriving from multi‐copy clades containing especially diversified loci. Each R gene recognized, on average, 2.42 of the 12 isolates screened. The abundant R genes identified in resistant genomes provide extraordinary redundancy in the ability of host genotypes to recognize particular isolates. If the same is true for other pathogens, many extant NBS‐LRR genes retain functionality. Our success at identifying rice blast R genes also validates a highly efficient cloning and screening strategy.     

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.     

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.