Case-control Genome-wide Association Studies Research Articles

Genome-wide association study for the level of prolificacy in Cameroon’s native goat

ABSTRACT Income from goats highly depends on prolificacy, which is difficult to improve by traditional breeding methods. The study aimed to identify SNP markers for prolificacy, using a case–control genome-wide association study (GWAS) on 111 genotyped Cameroon native goat (CNG) does, based on the 50 K single nucleotide polymorphism (SNP) chip panel. None of the top SNPs reached the significant p-value of 5 × 10−8. The highest p-value was 0.0009. Despite the number of cases being about a quarter of the number of controls, the highest allele frequency of some of the top 20 variants in the cases was indicative of their potential role in the trait. These top variants included the following 15: rs268285661, rs268235169, rs268236449, rs268235135, rs268240394 in Sphingosine-1-phosphate phosphatase 2 (SGPP2) gene, rs268283635 in Solute carrier family 24 member 2 (SLC24A2) gene, rs268251678 in Androgen-induced gene 1 (AIG1) gene, rs268267018, rs268239617, rs268281364, rs268273029, rs268286941, rs268236144, rs268233233 in CEP126 gene and rs268278159, respectively. Our findings indicate that GWAS enable the identification of some loci within genes, with known biological functions and pathways in human being and mice animal model but far-ranging to what was previously hypothesized and tested in goat.

Genetics may affect the risk of undergoing surgery for rhizarthrosis.

Osteoarthritis is a prevalent and severe disease. Involvement of the trapeziometacarpal joint is common and can lead to both pain and disability. Genetics are known to affect the risk of osteoarthritis, but it remains unclear how genetics affect disease trajectories. In this study, we investigated whether the genetic associations of trapeziometacarpal osteoarthritis (rhizarthrosis) vary with the need for surgical treatment. The study was conducted as a case-control genome-wide association study using individuals from the Copenhagen Hospital Biobank pain and degenerative musculoskeletal disease study and the Danish Blood Donor Study (N = 208,342). We identified patients diagnosed with rhizarthrosis and grouped them by treatment status, resulting in two case groups: surgical (N = 1083) and nonsurgical (N = 1888). The case groups were tested against osteoarthritis-free controls in two genome-wide association studies. We then compared variants suggestive of association (p < 10-6) in either of these analyses directly between the treatment groups (surgical vs. nonsurgical rhizarthrosis). We identified 10 variants suggestive of association with either surgical (seven variants) or nonsurgical (three variants) rhizarthrosis. None of the variants reached nominal significance in the opposite treatment group (p ≥ 0.14), and all 10 variants were significantly different between the treatment groups at a false discovery rate of 5%. These results suggest possible differences in the genetic associations of rhizarthrosis depending on surgical treatment. Clinical significance: Uncovering genetic differences between clinically distinct patient groups can reveal biological determinants of disease trajectories.

Genome-wide association study suggests genetic candidate loci of insulin dysregulation in Finnhorses

Equine metabolic syndrome (EMS) is a common welfare problem in horses worldwide. It is characterized by insulin dysregulation (ID), predisposition to laminitis and often obesity. EMS is multifactorial by nature, with both the environment and genetics contributing to the phenotype. Environmental factors, such as feeding and exercise, can be controlled, thus forming the basis for treatment and prevention. Genetic factors, by contrast, are less well-known and not easily controllable. The aim of this study was to identify potential genetic loci influencing ID/EMS in Finnhorses. A single-breed (Finnhorse) case-control genome-wide association study (GWAS) of ID was conducted with controls that included age-appropriate non-ID horses. ID status was determined with an oral sugar test (OST) for fasted horses. Seventy-one Finnhorses participated (n = 34 ID, n = 37 control). DNA samples (hair roots) were genotyped for 65 157 single-nucleotide polymorphisms (SNPs) with the Illumina Equine SNP70 BeadChip, and these data were analysed for association and FST outliers with genomic tools. P-values that exceeded the suggestive threshold (P = 1.00 ×10−5) were found in SNP BIEC2_383954 (P = 3.45 ×10−6) in chromosome 17 and SNP BIEC2_312374 (P = 1.89 ×10−5) in chromosome 15. Hierarchical and Bayesian FST outlier tests also detected these SNPs. Potential candidate genes associated with the ID close to SNP BIEC2_383954, with functions in carbohydrate metabolism, were Arginine and Glutamate Rich 1 (ARGLU1) and Ephrin-B2 (EFNB2).

Genome-wide case-control association study identifies a locus connecting NRF2-ARE regulation and selenium deficiency

Genome-wide case-control association study identifies a locus connecting NRF2-ARE regulation and selenium deficiency

Angular limb deformity associated with TSPAN18, NRG3 and NOVA2 in Rambouillet rams

Angular limb deformity (ALD) affects many species of livestock and companion animals. The mechanisms of ALD development are not well understood, but previous research suggests the involvement of genetic risk factors. A case-control genome-wide association study (GWAS) was conducted with 40 ALD-affected and 302 unaffected Rambouillet rams and 40,945 single nucleotide polymorphisms (SNPs). Forelimbs of 6 ALD-affected rams were examined and diagnosed with osteochondrosis. Genome-wide or chromosome-wide significant SNPs were positioned exonic, intronic or within the 3′UTR of genes TSPAN18, NRG3 and NOVA2, respectively. These genes have previously described roles related to angiogenesis and osteoblast, osteoclast and chondrocyte proliferation and differentiation, which suggests the possibility for their involvement in the pathogenesis of osteochondrosis. Functional consequences of SNPs were evaluated through transcription factor binding site analysis, which predicted binding sites for transcription factors of known importance to bone growth, including SOX6, SOX9 and RUNX2. The identification of genetic risk factors for ALD may help to improve animal welfare and production in Rambouillet, a breed known to be at risk for ALD development. This study proposes genes TSPAN18, NRG3 and NOVA2 as targets for further research towards understanding the etiology of ALD in Rambouillet sheep.

Lymphoma in Border Collies: Genome-Wide Association and Pedigree Analysis.

There has been considerable interest in studying cancer in dogs and its potential as a model system for humans. One area of research has been the search for genetic risk variants in canine lymphoma, which is amongst the most common canine cancers. Previous studies have focused on a limited number of breeds, but none have included Border Collies. The aims of this study were to identify relationships between Border Collie lymphoma cases through an extensive pedigree investigation and to utilise relationship information to conduct genome-wide association study (GWAS) analyses to identify risk regions associated with lymphoma. The expanded pedigree analysis included 83,000 Border Collies, with 71 identified lymphoma cases. The analysis identified affected close relatives, and a common ancestor was identified for 54 cases. For the genomic study, a GWAS was designed to incorporate lymphoma cases, putative "carriers", and controls. A case-control GWAS was also conducted as a comparison. Both analyses showed significant SNPs in regions on chromosomes 18 and 27. Putative top candidate genes from these regions included DLA-79, WNT10B, LMBR1L, KMT2D, and CCNT1.

ADuLT: An efficient and robust time-to-event GWAS

Proportional hazards models have been proposed to analyse time-to-event phenotypes in genome-wide association studies (GWAS). However, little is known about the ability of proportional hazards models to identify genetic associations under different generative models and when ascertainment is present. Here we propose the age-dependent liability threshold (ADuLT) model as an alternative to a Cox regression based GWAS, here represented by SPACox. We compare ADuLT, SPACox, and standard case-control GWAS in simulations under two generative models and with varying degrees of ascertainment as well as in the iPSYCH cohort. We find Cox regression GWAS to be underpowered when cases are strongly ascertained (cases are oversampled by a factor 5), regardless of the generative model used. ADuLT is robust to ascertainment in all simulated scenarios. Then, we analyse four psychiatric disorders in iPSYCH, ADHD, Autism, Depression, and Schizophrenia, with a strong case-ascertainment. Across these psychiatric disorders, ADuLT identifies 20 independent genome-wide significant associations, case-control GWAS finds 17, and SPACox finds 8, which is consistent with simulation results. As more genetic data are being linked to electronic health records, robust GWAS methods that can make use of age-of-onset information will help increase power in analyses for common health outcomes.

Genome-wide association study of abnormal elevation of ALT in patients exposed to atabecestat

BackgroundAtabecestat, a potent brain penetrable BACE1 inhibitor that reduces CSF amyloid beta (Aβ), was developed as an oral treatment for Alzheimer’s disease (AD). Elevated liver enzyme adverse events were reported in three studies although only one case met Hy’s law criteria to predict serious hepatotoxicity.MethodWe performed a case-control genome-wide association study (GWAS) to identify genetic risk variants associated with liver enzyme elevation using 42 cases with alanine transaminase (ALT) above three times the upper limit of normal (ULN) and 141 controls below ULN. Additionally, we performed a GWAS using continuous maximal ALT/ULN (expressed as times the ULN) upon exposure to atabecestat as the outcome measure (n = 285).ResultsNo variant passed the genome-wide significance threshold (p = 5 × 10− 8) in the case-control GWAS. We identified suggestive association signals in genes (NLRP1, SCIMP, and C1QBP) implicated in the inflammatory processes. Among the genes implicated by position mapping using variants suggestively associated (p < 1 × 10− 5) with ALT elevation case-control status, gene sets involved in innate immune response (adjusted p-value = 0.05) and regulation of cytokine production (adjusted p-value = 0.04) were enriched. One genomic region in the intronic region of GABRG3 passed the genome-wide significance threshold in the continuous max(ALT/ULN) GWAS, and this variant was nominally associated with ALT elevation case status (p = 0.009).ConclusionThe suggestive GWAS signals in the case-control GWAS analysis suggest the potential role of inflammation in atabecestat-induced liver enzyme elevation.

Investigating the Genetic Background of Spastic Syndrome in North American Holstein Cattle Based on Heritability, Genome-Wide Association, and Functional Genomic Analyses.

Spastic syndrome is a chronic, progressive disorder of adult cattle characterized by episodes of sudden involuntary muscle contractions or spasms of the extensor and abductor muscles of one or both hind limbs. In this study, a case-control genome-wide association study (GWAS) was performed on an adult Holstein cattle cohort. Based on the 50 K and high-density (HD) SNP panel GWAS, we identified 98 and 522 SNPs, respectively. The most significant genomic regions identified are located on BTA9 at approximately 87 megabase pairs (Mb) and BTA7 between 1 and 20 Mb. Functional analyses of significant SNPs identified genes associated with muscle contraction, neuron growth or regulation, and calcium or sodium ion movement. Two candidate genes (FIG4 and FYN) were identified. FIG4 is ubiquitously expressed in skeletal muscle and FYN is involved with processes such as forebrain development, neurogenesis, locomotion, neurogenesis, synapse development, neuron migration, and the positive regulation of neuron projection development. The CACNA1A gene, which codes for a calcium channel subunit protein in the calcium signaling pathway, seems the most compelling candidate gene, as many calcium ion channel disorders are non-degenerative, and produce spastic phenotypes. These results suggest that spastic syndrome is of polygenic inheritance, with important genomic areas of interest on BTA7 and BTA9.

Association of Risk Variants in the CFH Gene With Elevated Levels of Coagulation and Complement Factors in Idiopathic Multifocal Choroiditis

Idiopathic multifocal choroiditis (MFC) is poorly understood, thereby hindering optimal treatment and monitoring of patients. To identify the genes and pathways associated with idiopathic MFC. This was a case-control genome-wide association study (GWAS) and protein study of blood plasma samples conducted from March 2006 to February 2022. This was a multicenter study involving 6 Dutch universities. Participants were grouped into 2 cohorts: cohort 1 consisted of Dutch patients with idiopathic MFC and controls, and cohort 2 consisted of patients with MFC and controls. Plasma samples from patients with idiopathic MFC who had not received treatment were subjected to targeted proteomics. Idiopathic MFC was diagnosed according to the Standardization of Uveitis Nomenclature (SUN) Working Group guidelines for punctate inner choroidopathy and multifocal choroiditis with panuveitis. Data were analyzed from July 2021 to October 2022. Genetic variants associated with idiopathic MFC and risk variants associated with plasma protein concentrations in patients. This study included a total of 4437 participants in cohort 1 (170 [3.8%] Dutch patients with idiopathic MFC and 4267 [96.2%] controls; mean [SD] age, 55 [18] years; 2443 female [55%]) and 1344 participants in cohort 2 (52 [3.9%] patients with MFC and 1292 [96.1%] controls; 737 male [55%]). The primary GWAS association mapped to the CFH gene with genome-wide significance (lead variant the A allele of rs7535263; odds ratio [OR], 0.52; 95% CI, 0.41-0.64; P = 9.3 × 10-9). There was no genome-wide significant association with classical human leukocyte antigen (HLA) alleles (lead classical allele, HLA-A*31:01; P = .002). The association with rs7535263 showed consistent direction of effect in an independent cohort of 52 cases and 1292 control samples (combined meta-analysis OR, 0.58; 95% CI, 0.38-0.77; P = 3.0 × 10-8). In proteomic analysis of 87 patients, the risk allele G of rs7535263 in the CFH gene was strongly associated with increased plasma concentrations of factor H-related (FHR) proteins (eg, FHR-2, likelihood ratio test, adjusted P = 1.1 × 10-3) and proteins involved in platelet activation and the complement cascade. Results suggest that CFH gene variants increase systemic concentrations of key factors of the complement and coagulation cascades, thereby conferring susceptibility to idiopathic MFC. These findings suggest that the complement and coagulation pathways may be key targets for the treatment of idiopathic MFC.

Automatic block-wise genotype-phenotype association detection based on hidden Markov model

BackgroundFor detecting genotype-phenotype association from case–control single nucleotide polymorphism (SNP) data, one class of methods relies on testing each genomic variant site individually. However, this approach ignores the tendency for associated variant sites to be spatially clustered instead of uniformly distributed along the genome. Therefore, a more recent class of methods looks for blocks of influential variant sites. Unfortunately, existing such methods either assume prior knowledge of the blocks, or rely on ad hoc moving windows. A principled method is needed to automatically detect genomic variant blocks which are associated with the phenotype.ResultsIn this paper, we introduce an automatic block-wise Genome-Wide Association Study (GWAS) method based on Hidden Markov model. Using case–control SNP data as input, our method detects the number of blocks associated with the phenotype and the locations of the blocks. Correspondingly, the minor allele of each variate site will be classified as having negative influence, no influence or positive influence on the phenotype. We evaluated our method using both datasets simulated from our model and datasets from a block model different from ours, and compared the performance with other methods. These included both simple methods based on the Fisher’s exact test, applied site-by-site, as well as more complex methods built into the recent Zoom-Focus Algorithm. Across all simulations, our method consistently outperformed the comparisons.ConclusionsWith its demonstrated better performance, we expect our algorithm for detecting influential variant sites may help find more accurate signals across a wide range of case–control GWAS.

Genetic Validation of Psoriasis Phenotyping in UK Biobank Supports the Utility of Self-Reported Data and Composite Definitions for Large Genetic and Epidemiological Studies

Genetic Validation of Psoriasis Phenotyping in UK Biobank Supports the Utility of Self-Reported Data and Composite Definitions for Large Genetic and Epidemiological Studies

Revealing Genetic Variations Associated with Chip-Processing Properties in Potato (Solanum tuberosum L.)

Although the potato chip industry is booming, and distinct chip-processing clones have been released over the past 60 years, the genetic architecture of their chip-processing characteristics remains largely unknown. Case-control genome-wide association studies (GWAS) with SolCAP SNP array data for chip-processing clones versus all other market classes in the 393-line potato diversity panel were performed using the GWASpoly R package, enabling detection of significant signals on chromosome 10. Our results were replicated using internal replication of a strata-corrected 190-line panel. Furthermore, the genomic scans employing selective sweep approaches such as the cross-population composite likelihood ratio method (XP-CLR) and PCAdapt redetected the same signals as those in our GWAS. Through applications of four selective sweep approaches, various genetic variants were found across the genome that had been differentially selected. These genomic regions under selection along with transcriptomic data analysis are involved in carbohydrate metabolism-related genes or loci and transcription factors, indicating to be associated with the improvement of chip-processing performance of potato cultivars. Kompetitive allele-specific PCR (KASP) assays were designed for the causal SNPs to use in validating the chip-processing clones. The results could have implications for genomics-assisted breeding of the promising chip-processing cultivars in potato.

Genome-wide Association Study of Methotrexate-Induced Liver Injury in Rheumatoid Arthritis Patients.

Hepatotoxicity is a serious adverse drug reaction related to methotrexate (MTX). However, the cause of drug-induced liver injury (DILI) is still unclear and unpredictable. Genetic risk factors may predispose for MTX-DILI. Therefore, we conducted a nested case-control genome-wide association study to explore genetic risk factors associated with MTX-DILI. Seven international groups contributed blood samples and data of patients with rheumatoid arthritis who used MTX. MTX-DILI was defined as an alanine aminotransferase (ALT) level of at least three times the upper limit of normal (ULN), to increase contrast controls ALT levels did not raise above two times the ULN. Per study site, control subjects and patients with MTX-DILI (ratio 3:1) were matched for age, gender, and duration of MTX use. Patients were genotyped using Illumina GSA MD-24v1-0 and data were imputed using the 1000 Genomes reference panel. Single-nucleotide polymorphisms (SNPs) were analyzed using an additive genetic model, corrected for sex, country, and age. A P-value of ≤ 5 × 10-8 was considered significant, whereas a P-value of ≤ 5 × 10-6 was considered suggestive. A total of 108 MTX-DILI cases and 311 controls were included for association analysis. None of the SNPs were significantly associated with MTX-DILI. However, we found seven suggestive genetic variants associated with MTX-DILI (P-values 7.43 × 10-8 to 4.86 × 10-6 ). Of those, five SNPs were in the intronic protein-coding regions of FTCDNL1, BCOR, FGF14, RBMS3, and PFDN4/DOK5. Investigation of candidates SPATA9 (rs72783407), PLCG2 (rs60427389), RAVER2 (rs72675408), JAK1 (rs72675451), PTPN2 (rs2476601), MTHFR C677T (rs1801133), and into the HLA region did not show significant findings. No genetic variants associated with MTX-DILI were found, whereas suggestive SNPs need further investigation.

Risk factors for, and genetic association with, intestinal atresia in dairy calves.

Intestinal atresia is an under-diagnosed congenital defect in cattle. It results in complete occlusion of the intestinal lumen and, unless surgically corrected, results in death or euthanasia of the affected calf. There is limited information on the incidence of this condition or on risk factors, including predisposing alleles, associated with the defect. In this study, active surveillance of 39 dairy farms over 8 years identified 197 cases of intestinal atresia among 56 454 calves born, an incidence of 0.35%. The majority of cases (83%) had occlusion of the jejunum, although cases with blockage of the colon (14%) or anus (4%) were also identified. The defect was twice as common in male as in female calves (p < 0.0001), and was more common in progeny of older cows than in progeny of first or second lactation cows (p < 0.001). Year and farm of birth were also significantly associated with incidence (p < 0.05). The incidence of intestinal atresia was highest among the progeny of three related Jersey sires, suggesting that a gene for intestinal atresia was segregating within this family. Linkage analysis utilising 28 affected progeny of two half-sib putative carrier sires identified two putative quantitative trait loci associated with the defect, on chromosomes 14 and 26, although no clear candidate genes were identified. There was no evidence of a sire-effect among the progeny of Holstein-Friesian sires. However, a case-control genome-wide association study involving 91 cases and 375 healthy controls, identified 31 SNP in 18 loci as associated with the defect in this breed. These data suggest that intestinal atresia in dairy calves is not a simple Mendelian trait as previously reported but a complex multigenic disorder.

Identification of six genomic regions associated with pigment intensity in alpacas

Context Currently, there is a growing consumer demand for more ecologically sustainable practices in the textile industry. Fabric dyeing is highly pollutive, and one way to avoid dyeing is to use naturally coloured fibres. Alpacas exhibit a wide range of fleece colours, thereby making them a good source of fibre for sustainable textile production. Our understanding of the colour genetics of alpacas is improving, but there is still no explanation for all the variation seen in alpaca coat colours. Aims To identify a region or regions in the alpaca genome that contribute to differences in pigment intensity. Methods Colorimetric analysis using L*a*b* colourspace of fibre from white and black alpacas, that had been genotyped using the Neogen Australasia alpaca coat colour test, was used to classify the samples into intense and dilute groups for each colour. Illumina short-read genotyping by sequencing of the DNA from these alpacas was used to identify single nucleotide polymorphisms that were subsequently used in a case–control genome-wide association study (GWAS) comparing the extreme dilute and intensely pigmented animals. Key results Intense eumelanic fibre is darker (P = 0.0003), less red (P = 0.004), and more blue (P = 0.001) than is dilute eumelanic fibre. Intense pheomelanic fibre is darker (P = 1 × 10−7), more red (P = 3 × 10−20), and more yellow (P = 2 × 10−6) than is dilute pheomelanic fibre. The GWAS showed six regions of genome-wide significance. After manual inspection of these six regions, the best candidate region was upstream of KITLG, a gene previously associated with pigment intensity in dogs. Conclusions In combination with ASIP genotype, a regulatory mutation in a region upstream of KITLG in alpacas potentially has a major effect on pigment intensity in the species. Implications Successful identification of a marker for pigment intensity will allow breeders to select more precisely for breeding and production animals that will assist them in supplying the desired fibre colours to the textile industry.

Spotlight in Plastic Surgery: January 2023.

Spotlight in Plastic Surgery: January 2023.

Rationale, Design, and the Baseline Characteristics of the RHDGen (The Genetics of Rheumatic Heart Disease) Network Study†.

The genetics of rheumatic heart disease (RHDGen) Network was developed to assist the discovery and validation of genetic variations and biomarkers of risk for rheumatic heart disease (RHD) in continental Africans, as a part of the global fight to control and eradicate rheumatic fever/RHD. Thus, we describe the rationale and design of the RHDGen study, comprising participants from 8 African countries. RHDGen screened potential participants using echocardiography, thereafter enrolling RHD cases and ethnically-matched controls for whom case characteristics were documented. Biological samples were collected for conducting genetic analyses, including a discovery case-control genome-wide association study (GWAS) and a replication trio family study. Additional biological samples were also collected, and processed, for the measurement of biomarker analytes and the biomarker analyses are underway. Participants were enrolled into RHDGen between December 2012 and March 2018. For GWAS, 2548 RHD cases and 2261 controls (3301 women [69%]; mean age [SD], 37 [16.3] years) were available. RHD cases were predominantly Black (66%), Admixed (24%), and other ethnicities (10%). Among RHD cases, 34% were asymptomatic, 26% had prior valve surgery, and 23% had atrial fibrillation. The trio family replication arm included 116 RHD trio probands and 232 parents. RHDGen presents a rare opportunity to identify relevant patterns of genetic factors and biomarkers in Africans that may be associated with differential RHD risk. Furthermore, the RHDGen Network provides a platform for further work on fully elucidating the causes and mechanisms associated with RHD susceptibility and development.

Abstract 14588: A Large-Scale Genome-Wide Association Study of Angiographically Determined Burden of Coronary Atherosclerosis in a Genetically Diverse Population

Introduction: Few large-scale genome wide association studies (GWAS) of angiographically derived burden of coronary atherosclerosis have been reported to date. The yield of this approach relative to case-control GWAS of coronary artery disease (CAD) remains unclear. Methods: We used the Veterans Affairs Clinical Assessment, Reporting, and Tracking Program, a national quality and safety organization for invasive cardiac procedures, to conduct the largest ever GWAS of burden of coronary atherosclerosis involving participants of the Million Veteran Program (MVP). Obstructive disease of a native major vessel was defined as the presence of at least one obstruction &gt;50% of the luminal diameter, or a prior revascularization procedure, involving the main artery and/or its major branches. Non-obstructive disease was defined as one or more vessels having one or more obstructions of &gt;20% but none more than 50%. We classified an individual’s extent of disease as normal, non-obstructive, one vessel, two vessel, or three vessel / left main CAD. GWAS was performed stratified by racial/ethnic group using linear regression in REGENIE, with age, sex, and genetic principal components as covariates. Stratified results were then combined through meta-analysis using GWAMA. Results: The extent of disease was quantified in 49,213 non-Hispanic White, 11,326 non-Hispanic Black, and 4,058 Hispanic MVP participants, a majority of which were subjects with clinical CAD. Forty-three loci reached genome wide significance including 27 among Whites, four among Blacks at SORT1 , LPL , COL4A2 , and LDLR , two among Hispanics at CDKN2B-AS1 and NRG4 , and an additional 13 through meta-analysis of all three strata. Two loci, HHEX and NRG4 , have not been previously established as susceptibility loci CAD. Conclusions: Large-scale GWAS of burden of coronary atherosclerosis complements discovery efforts of case-control GWAS of CAD. Despite substantially smaller sample sizes, this quantitative approach recapitulated many established loci for CAD, uncovered novel loci, and identified several known susceptibility loci in admixed Black and Hispanic populations. The latter required a remarkably larger sample size to achieve through GWAS of clinically defined events alone.

A genome-wide association study to investigate genetic loci associated with primary glaucoma in American Cocker Spaniels.

To identify genetic associations with primary glaucoma (PG) in American Cocker Spaniels using a genome-wide association study (GWAS). A nationwide ambidirectional case-control cohort study was performed in American Cocker Spaniels that had an ophthalmic examination performed by a veterinarian. Ninety-four dogs with PG (cases) and 111 dogs without glaucoma (controls) met phenotypic criteria and had a blood sample collected after receiving informed owner consent. Genomic DNA was extracted from whole blood samples and genotyped (CanineHD BeadChip, Illumina Inc). A case-control GWAS using a linear mixed model was performed, and 3 significance thresholds were calculated (1) using a Bonferroni correction on all single nucleotide polymorphisms (SNPs) included in the GWAS, (2) using a Bonferroni correction on only the unlinked SNPs from a pruned data set, and (3) using 10,000 random phenotype permutations. Following genotype data quality control, 89 cases and 93 controls were included in the GWAS. We identified an association on canine chromosome (CFA10); however, it did not reach statistical significance. Potential candidate genes within the surrounding linkage disequilibrium interval include coiled-coil domain containing 85A (CCDC85A) and extracellular growth factor containing fibulin extracellular matrix protein 1 (EFEMP1). Primary glaucoma in the American Cocker Spaniel is a complex heterogeneous disease that may be influenced by a locus on CFA10. The candidate genes CCDC85A and EFEMP1 within the identified linkage disequilibrium interval have been shown to be involved in human open-angle glaucoma.