Sequence analysis of three canine adipokine genes revealed an association between TNF polymorphisms and obesity in Labrador dogs

Obesity is an emerging health problem in purebred dogs. Due to their crucial role in energy homeostasis control, genes encoding adipokines are considered candidate genes, and their variants may be associated with predisposition to obesity. Searching for polymorphism was carried out in three adipokine genes (TNF, RETN and IL6). The study was performed on 260 dogs, including lean (n = 109), overweight (n = 88) and obese (n = 63) dogs. The largest cohort was represented by Labrador Retrievers (n = 136). Altogether, 24 novel polymorphisms were identified: 12 in TNF (including one missense SNP), eight in RETN (including one missense SNP) and four in IL6. Distributions of five common SNPs (two in TNF, two in RETN and one in IL6) were further analyzed with regard to body condition score. Two SNPs in the non‐coding parts of TNF (c.‐40A>C and c.233+14G>A) were associated with obesity in Labrador dogs. The obtained results showed that the studied adipokine genes are highly polymorphic and two polymorphisms in the TNF gene may be considered as markers predisposing Labrador dogs to obesity.     

Novel SNPs in the <Emphasis Type="Italic">PRDM16</Emphasis> gene and their associations with performance traits in chickens

The PR domain containing 16 (PRDM16) is a member of the Prdm family, and is known to regulate cell differentiation. In the present study, DNA pool sequencing methods were employed to screen genetic variations in the chicken PRDM16 gene. The results revealed four novel single nucleotide polymorphisms (SNPs): NC_006108.2: g.92188G>A, XM_417551: c.1161C>T (Ala/Ala, 387aa), c.1233C>T (Ser/Ser, 411aa) and c.1433G>A (Ser/Asn, 478aa). The BglI polymerase chain reaction–restriction fragment length polymorphism (PCR-RFLP) was used to detect c.1161C>T, while HhaI Forced PCR-RFLP methods were used to detect 1233C>T and c.1433G>A in 964 chickens. The chickens comprised 38 grandparents, 66 F₁ parents and 860 F₂ birds derived from an F₂ resource population of Gushi chickens crossed with Anka broilers. The associations of the polymorphisms in the chicken PRDM16 gene with performance traits were analyzed in the 860 F₂ chickens. The results indicated that the three SNPs were significantly associated with growth, fatness and meat quality traits in the chickens. In particular, the polymorphisms of the missense SNP (c.1433G>A) had positive effects on chicken body weight and body size at different stages. It affected also fatness traits significantly. Comparison of the different genotypes of c.1433G>A showed that the GG genotype favored chicken growth and fatness traits.     

Analysis of non‐synonymous SNPs of the porcine SERPINA6 gene as potential causal variants for a QTL affecting plasma cortisol levels on SSC7

Recently, the SERPINA6 gene encoding corticosteroid‐binding globulin (CBG) has been proposed as a candidate gene for a quantitative trait locus (QTL) affecting cortisol level on pig chromosome 7. The QTL was repeatedly detected in different lines, including a Piétrain × (German Landrace × German Large White) cross (PiF1) and purebred German Landrace (LR). In this study, we investigated whether the known non‐synonymous polymorphisms c.44G>T, c.622C>T, c.770C>T, c.793G>A, c.832G>A and c.919G>A of SERPINA6 are sufficient to explain the QTL in these two populations. Our investigations revealed that SNPs c.44G>T, c.622C>T, c.793G>A and c.919G>A are associated with cortisol level in PiF1 (P < 0.01). Haplotype analysis showed that these associations are largely attributable to differences between a major haplotype carrying SNPs c.793G>A and c.919G>A and a haplotype carrying SNPs c.44G>T and c.622C>T. Furthermore, some SNPs, particularly c.44G>T and c.622C>T and the carrier haplotype, showed association with meat quality traits including pH and conductivity (P < 0.05). In LR, the non‐synonymous SNPs segregate at very low frequency (<5%) and/or show only weak association with cortisol level (SNPs c.832G>A and c.919G>A; P < 0.05). These findings suggest that the non‐synonymous SNPs are not sufficient to explain the QTL across different breeds. Therefore, we examined whether the expression of SERPINA6 is affected by cis‐regulatory polymorphisms in liver, the major organ for CBG production. We found allelic expression imbalance of SERPINA6, which suggests that its expression is indeed affected by genetic variation in cis‐acting elements. This represents candidate causal variation for future studies of the molecular background of the QTL.     

Plasminogen alleles influence susceptibility to invasive aspergillosis

Invasive aspergillosis (IA) is a common and life-threatening infection in immunocompromised individuals. A number of environmental and epidemiologic risk factors for developing IA have been identified. However, genetic factors that affect risk for developing IA have not been clearly identified. We report that host genetic differences influence outcome following establishment of pulmonary aspergillosis in an exogenously immune suppressed mouse model. Computational haplotype-based genetic analysis indicated that genetic variation within the biologically plausible positional candidate gene plasminogen (Plg; Gene ID 18855) correlated with murine outcome. There was a single nonsynonymous coding change (Gly110Ser) where the minor allele was found in all of the susceptible strains, but not in the resistant strains. A nonsynonymous single nucleotide polymorphism (Asp472Asn) was also identified in the human homolog (PLG; Gene ID 5340). An association study within a cohort of 236 allogeneic hematopoietic stem cell transplant (HSCT) recipients revealed that alleles at this SNP significantly affected the risk of developing IA after HSCT. Furthermore, we demonstrated that plasminogen directly binds to Aspergillus fumigatus. We propose that genetic variation within the plasminogen pathway influences the pathogenesis of this invasive fungal infection.     

Association between COMT SNP variation and timidity in Golden and Labrador Retrievers

Timidity in dogs is a trait with high heritability. However, the relevant genetic factors and markers associated with this condition are largely unknown. The function of the catechol‐O‐methyl transferase (COMT) gene has been found to be associated with human fearful or anxious emotions, and the COMT:p.Val158Met polymorphism locus is significantly related to anxious behavior. In the present study, the correlation between timidity and four single nucleotide polymorphism (SNP) variations (C.‐1666C>G c.39A>G, c.216G>A, c.482G>A) of the COMT gene was investigated in dogs. The evaluation was based on the dog courage assessment assay and a genotype and haplotype analysis in Labrador Retrievers (LR) and Golden Retrievers (GR). The principal components analysis factor structure of the courage phenotype was invariant between LR and GR. Sex, breed and age had no statistically significant effect on the timidity of the dogs. All SNP loci detected were in Hardy–Weinberg equilibrium. The c.39A>G locus was removed in the subsequent association analysis due to the significant difference between LR and GR in genotype distributions. Intriguingly, the genotypes and haplotypes of the COMT gene were significantly and highly correlated with the timidity of LR and GR. The A alleles of the COMT:c.216G>A and c.482G>A loci appeared to play a dominant role in the timid behavior of the dogs. This result supports and broadens the warrior/worrier hypothesis and will have important implications for the understanding of the evolution of temperament in dogs. Additionally, the results provide predictive genetic markers for temperament in dogs.     

Association between single nucleotide polymorphisms of TPH1 and TPH2 genes,and depressive disorders

Tryptophan catabolites pathway disorders are observed in patients with depression. Moreover, single nucleotide polymorphisms of tryptophan hydroxylase genes may modulate the risk of depression occurrence. The objective of our study was to confirm the association between the presence of polymorphic variants of TPH1 and TPH2 genes, and the development of depressive disorders. Six polymorphisms were selected: c.804‐7C>A (rs10488682), c.‐1668T>A (rs623580), c.803+221C>A (rs1800532), c.‐173A>T (rs1799913)—TPH1, c.‐1449C>A (rs7963803), and c.‐844G>T (rs4570625)—TPH2. A total of 510 DNA samples (230 controls and 280 patients) were genotyped using TaqMan probes. Among the studied polymoorphisms, the G/G genotype and G allele of c.804‐7C>A—TPH1, the T/T homozygote of c.803+221C>A—TPH1, the A/A genotype and A allele of c.1668T>A—TPH1, the G/G homozygote and G allele of c.‐844G>T—TPH2, and the C/A heterozygote and A allele of c.‐1449C>A—TPH2 were associated with the occurrence of depression. However, the T/T homozygote of c.‐1668T>A—TPH1, the G/T heterozygote and T allele of c.‐844G>T—TPH2, and the C/C homozygote and C allele of c.‐1449C>A—TPH2 decreased the risk of development of depressive disorders . Each of the studied polymorphisms modulated the risk of depression for selected genotypes and alleles. These results support the hypothesis regarding the involvement of the pathway in the pathogenesis of depression.     

A CDH23 missense variant in Beauceron dogs with non-syndromic deafness

Congenital coat-colour-related deafness is common among certain canine breeds especially those exhibiting extreme white spotting or merle patterning. We identified a non-syndromic deafness in Beauceron dogs characterised by a bilateral hearing loss in puppies that is not linked to coat colour. Pedigree analysis suggested an autosomal recessive transmission. By combining homozygosity mapping with whole genome sequencing and variant filtering in affected dogs we identified a CDH23:c.700C>T variant. The variant, located in the CHD23 (cadherin related 23) gene, was predicted to induce a CDH23:p.(Pro234Ser) change in the protein. Proline-234 of CDH23 protein is highly conserved across different vertebrate species. In silico tools predicted the CDH23:p.(Pro234Ser) change to be deleterious. CDH23 encodes a calcium-dependent transmembrane glycoprotein localised near the tips of hair-cell stereocilia in the mammalian inner ear. Intact function of these cilia is mandatory for the transformation of the acoustical wave into a neurological signal, leading to sensorineural deafness when impaired. By genotyping a cohort of 90 control Beauceron dogs sampled in France, we found a 3.3% carrier frequency. The CDH23:c.[700C>T] allele is easily detectable with a genetic test to avoid at-risk matings.     

Genetic polymorphisms of XRCC1 gene and susceptibility to gastric cancer in Chinese Han population

Abstract

This study aims to evaluate whether the c.1471G?>?A and c.1686C?>?G genetic polymorphisms of XRCC1 gene influencing gastric cancer susceptibility. A total of 813 subjects with Chinese Han ethnicity were enrolled. Our data suggest that the allele and genotype frequencies are significantly different from gastric cancer patients with cancer-free controls. We find that c.1471G?>?A and c.1686C?>?G genetic polymorphisms statistically increase the risk of gastric cancer. Our findings indicate these two genetic polymorphisms are related with the susceptibility to gastric cancer, and could be used as molecular markers for detecting gastric cancer in Chinese Han ethnicity.     

Associations of polymorphisms in the promoter I of bovine acetyl‐CoA carboxylase‐α gene with beef fatty acid composition

The objectives of this study were to identify single nucleotide polymorphisms (SNPs) in the promoter I (PI) region of the bovine acetyl‐CoA carboxylase‐α (ACACA) gene and to evaluate the extent to which they were associated with lipid‐related traits. Eight novel SNPs were identified, which were AJ276223:g.2064T>A (SNP1), g.2155C>T (SNP2), g.2203G>T (SNP3), g.2268T>C (SNP4), g.2274G>A (SNP5), g.2340A>G (SNP6), g.2350T>C (SNP7) and g.2370A>G (SNP8). Complete linkage disequilibrium was observed among SNP1, 2, 4, 5, 6 and 8. Phenotypic data were collected from 573 cross‐bred steers with six sire breeds, including Hereford, Angus, Brangus, Beefmaster, Bonsmara and Romosinuano. The genotypes of SNP1/2/4/5/6/8 were significantly associated with adjusted backfat thickness. The genotypes of SNP3 were significantly associated with triacylglycerol (TAG) content and fatty acid composition of longissimus dorsi muscle (LM) in Brangus‐, Romosinuano‐ and Bonsmara‐sired cattle. Cattle with g.2203GG genotype had greater concentrations of TAG, total lipid, total saturated fatty acid and total monounsaturated fatty acid than did cattle with g.2203GT genotype. The genotypes of SNP7 were significantly associated with fatty acid composition of LM. Cattle with genotype g.2350TC had greater amounts of several fatty acids in LM than did cattle with genotype g.2350CC. Our results suggested that the SNPs in the PI region of ACACA gene are associated with variations in the fatty acid contents in LM.     

Cardiac ion channel gene mutations in Greek long QT syndrome patients

The long QT syndrome (LQTS) is an inherited cardiac arrhythmia that may lead to sudden death in the absence of structural heart disease. Mutations in the cardiac potassium and sodium channel genes can be found in approximately 70% of patients with a highly probable clinical diagnosis. In this study, we aimed to genotype and explore the yield of genetic testing of LQTS patients from Greece, for whom there are no collective published data available. We clinically evaluated and genetically screened 17 unrelated patients for mutations in theKCNQ1, KCNH2, SCN5A, KCNE1, andKCNE2 cardiac ion channel genes. Genetic testing was positive in 6 out of 8 patients with a highly probable clinical diagnosis of LQTS and negative for all the other patients. Two patients carriedKCNQ1 mutations (c.580G>C, c.1022C>T), while 4 patients carriedKCNH2 mutations (c.202T>C, c.1714G>A, c.3103delC, c.3136C>T). To the best of our knowledge, the last mentioned mutation (c.3136C>T) is novel. Moreover, 27 single-nucleotide polymorphisms (SNPs) were detected, 5 of which are novel. Our preliminary data indicate a low genetic diversity of the Greek LQTS genetic pool, and are in accordance with international data that genetic testing of the major LQTS genes is efficient in genotyping the majority of patients with a strong clinical diagnosis. Therefore, the transition of an LQTS genetic screening program from research to the diagnostic setting within our ethnic background is feasible.     

A new single nucleotide polymorphism in the rabbit (Oryctolagus cuniculus) myostatin (MSTN) gene is associated with carcass composition traits

This study aimed at the identification of genetic variations in the myostatin (MSTN) gene and testing their effects on carcass quality traits. We comparatively sequenced Giant Grey (GG) and New Zealand White (NZW) rabbits that were founders of a cross‐bred population. Alignment of our sequence data with the GenBank sequence of the rabbit MSTN gene (Ensembl Gene ID ENSOCUG00000012663) identified three single nucleotide polymorphisms (SNPs). The two novel SNPs (c.?125T>C, c.373+234G>A) and one known SNP (c.747+34C>T) were subsequently analysed for linkage with carcass composition traits in 363 F₂ animals of the cross GG × NZW. Significant linkage was found between c.373+234G>A and nine carcass composition traits (P < 0.05). No significant effects were found for c.?125T>C and c.747+34C>T. Because the linked SNP is located in intron 1 and no genetic variation was found in the coding region, further investigations are necessary to understand the functional effect of the c.373+234G>A variant on the variability of the traits.     

SLCO2B1 genetic polymorphisms in a Korean population: pyrosequencing analyses and comprehensive comparison with other populations

SLCO2B1, also known as OATP2B1 (Organic Anion Transporter) or OATP-B or SLC21A9, is an organic anion uptake transporter that is encoded by the SLCO2B1 gene. In this study we assessed the frequencies of SLCO2B1 polymorphisms in a Korean population using newly developed pyrosequencing methods and compared their frequencies with those in other ethnic groups. We developed pyrosequencing methods to identify the following six SLCO2B1 non-synonymous polymorphisms: c.1175C > T (rs1621378), c.1457C > T (rs2306168), c.43C > T (rs56837383), c.935G > A (rs12422149), c.601G > A (rs35199625) and c.644A > T (rs72559740). The allele frequencies of these polymorphisms were analyzed in 227 Korean subjects. The allele frequencies of SLCO2B1 polymorphisms in the population tested were as follows: 0.0 for c.1175C > T, c.43C > T and c.644A > T; 0.2687 for c.1457C > T; 0.4273 for c.935G > A; and 0.0727 for c. 601G > A. Even though the allele frequencies of the c.1175C > T and c.1457C > T polymorphisms were comparable to those in Japanese subjects, the frequencies in this Korean population differed from those in other ethnic groups. The developed pyrosequencing methods are rapid and reliable for detecting non-synonymous SLCO2B1 polymorphisms. Large ethnic differences in the frequency of SLCO2B1 genetic polymorphisms were noted among ethnic groups. The SLCO2B1 polymorphisms at c.1175C > T, c.43C > T and c.644A > T were not found in the Korean population while c.1457C > T, c.935G > A and c.601G > A exhibited mostly higher frequencies in Koreans compared with Finnish, Caucasian and African-American populations.     

SNP genotypes of olfactory receptor genes associated with olfactory ability in German Shepherd dogs

To find out the relationship between SNP genotypes of canine olfactory receptor genes and olfactory ability, 28 males and 20 females from German Shepherd dogs in police service were scored by odor detection tests and analyzed using the Beckman GenomeLab SNPstream. The representative 22 SNP loci from the exonic regions of 12 olfactory receptor genes were investigated, and three kinds of odor (human, ice drug and trinitrotoluene) were detected. The results showed that the SNP genotypes at the OR10H1‐like:c.632C>T, OR10H1‐like:c.770A>T, OR2K2‐like:c.518G>A, OR4C11‐like:c.511T>G and OR4C11‐like:c.692G>A loci had a statistically significant effect on the scenting abilities (P < 0.001). The kind of odor influenced the performances of the dogs (P < 0.001). In addition, there were interactions between genotype and the kind of odor at the following loci: OR10H1‐like:c.632C>T, OR10H1‐like:c.770A>T, OR4C11‐like:c.511T>G and OR4C11‐like:c.692G>A (P < 0.001). The dogs with genotype CC at the OR10H1‐like:c.632C>T, genotype AA at the OR10H1‐like:c.770A>T, genotype TT at the OR4C11‐like:c.511T>G and genotype GG at the OR4C11‐like:c.692G>A loci did better at detecting the ice drug. We concluded that there was linkage between certain SNP genotypes and the olfactory ability of dogs and that SNP genotypes might be useful in determining dogs' scenting potential.     

MYLK polymorphism associated with blood eosinophil level among asthmatic patients in a Korean population

The myosin light chain kinase (MYLK) gene encodes both smooth muscle and nonmuscle cell isoforms. Recently, polymorphisms in MYLK have been reported to be associated with several diseases. To examine the genetic effects of polymorphisms on the risk of asthma and related phenotypes, we scrutinized MYLK by re-sequencing/genotyping and statistical analysis in Korean population (n = 1,015). Seventeen common polymorphisms located in or near exons, having pairwise r ² values less than 0.25, were genotyped. Our statistical analysis did not replicate the associations with the risk of asthma and log-transformed total IgE levels observed among African descendant populations. However, two SNPs in intron 16 (+89872C > G and +92263T >C), which were in tight LD (|D′| = 0.99), revealed significant association with log-transformed blood eosinophil level even after correction multiple testing (P = 0.002/P ^corr = 0.01 and P = 0.002/P ^corr = 0.01, respectively). The log-transformed blood eosinophil levels were higher in individuals bearing the minor alleles for +89872C > G and +92263T > C, than in those bearing other allele. In additional subgroup analysis, the genetic effects of both SNPs were much more apparent among asthmatic patients and atopic asthma patients. Among atopic asthma patients, the log-transformed blood eosinophil levels were proportionally increased by gene-dose dependent manner of in both +89872C > G and +92263T > C (P = 0.0002 and P = 0.00007, respectively). These findings suggest that MYLK polymorphisms might be among the genetic factors underlying differential increases of blood eosinophil levels among asthmatic patients. Further biological and/or functional studies are needed to confirm our results.     

The prion‐related protein (testis‐specific) gene (PRNT) is highly polymorphic in Portuguese sheep

The objective of this study was to search for polymorphisms in the ovine prion‐related protein (testis‐specific) gene (PRNT). Sampling included 567 sheep from eight Portuguese breeds. The PRNT gene‐coding region was analyzed by single‐strand conformation polymorphism and sequencing, allowing the identification of the first ovine PRNT polymorphisms, in codons 6, 38, 43 and 48: c.17C>T (p.Ser6Phe, which disrupts a consensus arginine‐X‐X‐serine/threonine motif); c.112G>C (p.Gly38>Arg); c.129T>C and c.144A>G (synonymous) respectively. Polymorphisms in codons 6, 38 and 48 occur simultaneously in 50.6% of the animals, 38.8% presenting as heterozygous. To study the distribution of the polymorphism in codon 43, a restriction fragment length polymorphism analysis was performed. Polymorphic variant c.129C, identified in 89.8% of the animals with 32.8% presented as heterozygous, was considered the wild genotype in Portuguese sheep. Eight different haplotypes which have comparable distribution in all breeds were identified for the PRNT gene. In conclusion, the PRNT coding region is highly polymorphic in sheep, unlike the prion protein 2 dublet gene (PRND), in which we previously found only one synonymous substitution (c.78G>A), in codon 26. The absence or reduced number of PRND heterozygotes (c.78G>A) was significantly associated with three PRNT haplotypes (17C‐112G‐129T‐144A,17CT‐112GC‐129CT‐144AG and 17T‐112C‐129C‐144G), and the only three animals found homozygous at c.78A had the 17C‐112G‐129C‐144A PRNT haplotype. These results constitute evidence of an association between polymorphic variation in PRND and PRNT genes, as has already been observed for PRND and prion protein gene (PRNP).     

Association analysis of v-AKT murine thymoma viral oncogene homolog 1 (AKT1) polymorphisms and type 2 diabetes mellitus in the Korean population

V-AKT murine thymoma viral oncogene homolog 1 (AKT1) is an important downstream target of the insulin-signaling pathway and may be an important regulator of pancreatic beta cell growth. This study investigated the association of theAKT1 gene with susceptibility to type 2 diabetes mellitus and its related traits. By sequencing theAKT1 gene in 24 unrelated individuals, we iden-tified 32 genetic variations including 30 single nucleotide polymorphisms and 2 deletions. For the association analysis, we selected seven single nucleotide polymorphisms (rs10138227, ?726G>A; rs3730358, +12574C>T; rs2494737, +12656T>A; rs2498796, +15761T>C; rs2498799, +19087 A>G; rs2494732, +19789G>A; rs3803304, +19835G>C) based on minor allele frequency (>0.05) and linkage disequilibrium status. The study included 483 type 2 diabetes patients (206 men and 277 women with mean age 64±2.8 years and mean age at onset 56 ± 8.1 years) and 1,138 non-diabetic control subjects (516 men and 622 women with mean age 64 ±2.9 years). Two single nucleotide polymorphisms (rs2498796, +15761T>C and rs2494732, +19789G>A) were found to be associated with risk of type 2 diabetes mellitus, and showed an increased risk of type 2 diabetes mellitus in a recessive model (OR=1.343, 95% CI 1.021–1.765,p=0.035 and OR=1.534, 95% CI 1.058–2.225,p=0.024, respectively). These SNPs were also associated with diabetes-related traits such as levels of fasting blood glucose and hemoglobin A1c. In addition, type 2 diabetes mellitus patients who also have dyslipidemia or high blood pressure showed significant association with single nucleotide polymorphisms in AKT1 when compared with healthy controls. These results indicate that genetic variation in AKT1 influences the development of type 2 diabetes mellitus in the Korean population.     

The contribution of SIPA1 and RRP1B germline polymorphisms to breast cancer phenotype,lymph node status and survival in a group of Lithuanian young breast cancer patients

The germline polymorphisms in signal-inducing proliferation-associated protein 1 (SIPA1) and ribosomal RNR processing 1B (RRP1B) might be involved in breast cancer metastasis. The aim of this study was to analyze how SIPA1 and RRP1B polymorphisms contribute to breast cancer phenotype, lymph node status and survival. A group of 100 young, I–II stage breast cancer patients were analyzed for SIPA1 and RRP1B polymorphisms with PCR-RFLP assay. SIPA1 c.2760G?>?A, c.545C?>?T and RRP1B c.436T?>?C polymorphisms were associated with lymph node status, survival and tumor grade, respectively. Our results suggest that SIPA1 and RRP1B germline polymorphisms are important for breast cancer prognosis.     

Studies of type I collagen (<Emphasis Type="Italic">COL1A1</Emphasis>) α1 chain in patients with osteogenesis imperfecta

Nucleotide sequences of exon 51, adjacent intron areas, and regulatory region of the α1 chain of type I collagen (COL1A1) gene were analyzed in 41 patients with osteogenesis imperfecta (OI) from 33 families and their 68 relatives residing at Bashkortostan Republic (BR). Six mutations (four nonsense mutations c.967G>T (p.Gly323X), c.1081C>T (p.Arg361X), c.1243C>T (p.Arg415X), and c.2869C>T (p.Gln957X)) in patients of the Russian origin and two frameshift mutations (c.579delT (p.Gly194ValfsX71), and c.2444delG (p.Gly815AlafsX293)) in patients with OI of Tatar ethnicity as well as 14 single nucleotide polymorphisms in the COL1A1 gene were revealed. Mutations c.967G>T (p.Gly323X) and three alterations in the nucleotide sequence c.544-24C>T, c.643-36delT, and c.957 + 10insA were described for the first time.     

Targeted next-generation sequencing for identifying genes related to horse temperament

It is a fundamental challenge to discover the association of genetic traits with phenotypic traits. In this study, we aimed to identify possible genetic traits related to horse temperament. Based on previous findings, we selected 71 candidate genes related to temperamental trait and examined them in the human and horse reference genomes (hg38 and equCab2, respectively). We found 16 orthologous genes and, by comparing with the human reference genome, 17 homologous genes in the horse reference genome. We designed probes specific for the 33 horse genes. Using the probes, we built sequencing libraries of the genomic DNA samples from eight aggressive and eight docile horses, and sequenced the constructed libraries using the Illumina Hiseq2500 platform. Through the analysis of the targeted exome sequences, we identified single nucleotide polymorphisms (SNPs) in the genes. SNPs could be served as genetic markers to evaluate aggressive or docile levels of horses. To examine whether any genetic variants are associated with horse temperament, we performed genome-wide association study (GWAS) using the SNP data. GWAS analysis identified ten variants (p-value?<0.05) which could be related to horse temperament. We validated the variants using Sanger sequencing. The most significant variants were found in MAOA (c.1164+41T>C) and AR (c.1047+27G>T) genes with 8.09?×?10^?4 p-value. We suggest that the variants might be used to assess horse temperament and to determine superior horses for riding or racing.