The GenABEL Project for statistical genomics

Potentially functional variants of INPP5D and EXOSC3 in immunity B cell-related genes are associated with non-small cell lung cancer survival

B cells are adaptive immune cells in the tumor microenvironment and play an important role in tumor development and metastasis. However, the roles of genetic variants of the immunity B cell-related genes in the survival of patients with non-small cell lung cancer (NSCLC) remain unknown. In the present study, we first evaluated associations between 10,776 single nucleotide polymorphisms (SNPs) in 220 immunity B cell-related genes and survival of NSCLC in a discovery dataset of 1,185 patients from the Prostate, Lung, Colorectal and Ovarian (PLCO) Cancer Screening Trial. We found that 369 SNPs were significantly associated with overall survival (OS) of NSCLC in multivariable Cox proportional hazards regression analysis (P ≤ 0.05, Bayesian false discovery probability ≤ 0.80), of which 18 SNPs were validated in another independent genotyping dataset of 984 patients from the Harvard Lung Cancer Susceptibility (HLCS) Study. We then performed linkage disequilibrium (LD) analysis, followed by stepwise analysis with a multivariable Cox regression model. Finally, two independent SNPs, inositol polyphosphate-5-phosphatase D (INPP5D) rs13385922 C>T and exosome component 3 (EXOSC3) rs3208406 A>G, remained significantly associated withNSCLC OS with a combined hazards ratio (HR) of 1.14 (95% confidence interval = 1.06-1.23, P = 2.41×10-4) and 1.20 (95% confidence interval = 1.14-1.28, P = 3.41×10-9), respectively. Furthermore, NSCLC patients with the combination of unfavorable genotypes for these two SNPs were associated with a poor OS (P trend = 0.0002) and disease-specific survival (DSS, P trend < 0.0001) in the PLCO dataset. Expression quantitative trait loci (eQTL) analysis suggested that the INPP5D rs6782875 T allele was significantly correlated with elevated INPP5D mRNA expression levels in normal lung tissues and whole blood samples, while the EXOSC3 rs3208406 G allele was significantly correlated with increased EXOSC3 mRNA expression levels in normal lung tissues. Our data indicated that genetic variants in these immunity B cell-related genes may predict NSCLC survival possibly by influencing the gene expression.

Identifying Regions of the Genome Associated with Conception Rate to the First Service in Holstein Heifers Bred by Artificial Insemination and as Embryo Transfer Recipients

Heifer conception rate to the first service (HCR1) is defined as the number of heifers that become pregnant to the first breeding service compared to the heifers bred. This study aimed to identify loci associated and gene sets enriched for HCR1 for heifers that were bred by artificial insemination (AI, n = 2829) or were embryo transfer (ET, n = 2086) recipients, by completing a genome-wide association analysis and gene set enrichment analysis using SNP data (GSEA-SNP). Three unique loci, containing four positional candidate genes, were associated (p < 1 × 10-5) with HCR1 for ET recipients, while the GSEA-SNP identified four gene sets (NES ≥ 3) and sixty-two leading edge genes (LEGs) enriched for HCR1. While no loci were associated with HCR1 bred by AI, one gene set and twelve LEGs were enriched (NES ≥ 3) for HCR1 with the GSEA-SNP. This included one gene (PKD2) shared between HCR1 AI and ET services. Identifying loci associated or enriched for HCR1 provides an opportunity to use them as genomic selection tools to facilitate the selection of cattle with higher reproductive efficiency, and to better understand embryonic loss.

A Guide to Genome-Wide Association Study Design for Diabetic Retinopathy

Diabetic retinopathy (DR) is the most common microvascular complication related to diabetes. There is evidence that genetics play an important role in DR pathogenesis, but the complexity of the disease makes genetic studies a challenge. This chapter is a practical overview of the basic steps for genome-wide association studies with respect to DR and its associated traits. Also described are approaches that can be adopted in future DR studies. This is intended to serve as a guide for beginners and to provide a framework for further in-depth analysis.

Genome-wide association study for reproduction-related traits in Chinese domestic goose

1. This study measured six reproduction traits in a Sichuan white goose population (209 individuals), including fertility, qualified egg rate, plasma concentrations of progesterone (P), follicle-stimulating hormone (FSH), prolactin (PRL) and oestrogen (E2).2. Whole-genome resequencing data from the same goose population (209 individuals) were used in a genome-wide association study (GWAS) utilising a mixed linear model to investigate the genes and genetic markers associated with reproduction traits. The frequency of the selected SNPs and haplotypes were determined using the Matrix-Assisted Laser Desorption Ionisation Time-Of-Flight Mass Spectrometry (MALDI-TOF MS) method.3. In total, 42 SNPs significantly associated with these traits were identified. A haplotype block was constructed based on five SNPs that were significantly associated with qualified egg rate, with individuals having the haplotype CCTTAAGGAA having the lowest qualified egg rate.4. In conclusion, these results provided potential markers for marker-assisted selection to improve goose reproductive performance and a basis for elucidating the genetics of goose reproduction.

The APOE locus is linked to decline in general cognitive function: 20-years follow-up in the Doetinchem Cohort Study

Cognitive decline is part of the normal aging process. However, some people experience a more rapid decline than others due to environmental and genetic factors. Numerous single nucleotide polymorphisms (SNPs) have been linked to cognitive function, but only a few to cognitive decline. To understand whether cognitive function and cognitive decline are driven by the same mechanisms, we investigated whether 433 SNPs previously linked to cognitive function and 2 SNPs previously linked to cognitive decline are associated with both general cognitive functioning at baseline and general cognitive decline up to 20-years follow-up in the Doetinchem Cohort Study (DCS). The DCS is a longitudinal population-based study that enrolled men and women aged 20-59 years between 1987-1991, with follow-up examinations every 5 years. We used data of rounds 2-6 (1993-2017, n = 2559). General cognitive function was assessed using four cognition tests measuring memory, speed, fluency and flexibility. With these test scores, standardized residuals (adjusted for sex, age and examination round) were calculated for each cognition test at each round and subsequently combined into one general cognitive function measure using principal component analyses. None of the 435 previously identified variants were associated with baseline general cognitive function in the DCS. But rs429358-C, a coding apolipoprotein E (APOE) SNP and one of the variants previously associated with cognitive decline, was associated with general cognitive decline in our study as well (p-value = 1 × 10^-5, Beta = -0.013). These findings suggest that decline of general cognitive function is influenced by other mechanisms than those that are involved in the regulation of general cognitive function.

Children at onset of type 1 diabetes show altered N-glycosylation of plasma proteins and IgG

AIMS/HYPOTHESIS

Individual variation in plasma N-glycosylation has mainly been studied in the context of diabetes complications, and its role in type 1 diabetes onset is largely unknown. Our aims were to undertake a detailed characterisation of the plasma and IgG N-glycomes in patients with recent onset type 1 diabetes, and to evaluate their discriminative potential in risk assessment.

METHODS

In the first part of the study, plasma and IgG N-glycans were chromatographically analysed in a study population from the DanDiabKids registry, comprising 1917 children and adolescents (0.6-19.1 years) who were newly diagnosed with type 1 diabetes. A follow-up study compared the results for 188 of these participants with those for their 244 unaffected siblings. Correlation of N-glycan abundance with the levels and number of various autoantibodies (against IA-2, GAD, ZnT8R, ZnT8W), as well as with sex and age at diagnosis, were estimated by using general linear modelling. A disease predictive model was built using logistic mixed-model elastic net regression, and evaluated using a 10-fold cross-validation.

RESULTS

Our study showed that onset of type 1 diabetes was associated with an increase in the proportion of plasma and IgG high-mannose and bisecting GlcNAc structures, a decrease in monogalactosylation, and an increase in IgG disialylation. ZnT8R autoantibody levels were associated with higher IgG digalactosylated glycan with bisecting GlcNAc. Finally, an increase in the number of autoantibodies (which is a better predictor of progression to overt diabetes than the level of any individual antibody) was accompanied by a decrease in the proportions of some of the highly branched plasma N-glycans. Models including age, sex and N-glycans yielded notable discriminative power between children with type 1 diabetes and their healthy siblings, with AUCs of 0.915 and 0.869 for addition of plasma and IgG N-glycans, respectively.

CONCLUSIONS/INTERPRETATION

We defined N-glycan changes accompanying onset of type 1 diabetes, and developed a predictive model based on N-glycan profiles that could have valuable potential in risk assessment. Increasing the power of tests to identify individuals at risk of disease development would be a considerable asset for type 1 diabetes prevention trials.

Genetic regulation of post-translational modification of two distinct proteins

Post-translational modifications diversify protein functions and dynamically coordinate their signalling networks, influencing most aspects of cell physiology. Nevertheless, their genetic regulation or influence on complex traits is not fully understood. Here, we compare the genetic regulation of the same PTM of two proteins - glycosylation of transferrin and immunoglobulin G (IgG). By performing genome-wide association analysis of transferrin glycosylation, we identify 10 significantly associated loci, 9 of which were not reported previously. Comparing these with IgG glycosylation-associated genes, we note protein-specific associations with genes encoding glycosylation enzymes (transferrin - MGAT5, ST3GAL4, B3GAT1; IgG - MGAT3, ST6GAL1), as well as shared associations (FUT6, FUT8). Colocalisation analyses of the latter suggest that different causal variants in the FUT genes regulate fucosylation of the two proteins. Glycosylation of these proteins is thus genetically regulated by both shared and protein-specific mechanisms.

Linking genetic, morphological, and behavioural divergence between inland island and mainland deer mice

The island syndrome hypothesis (ISH) stipulates that, as a result of local selection pressures and restricted gene flow, individuals from island populations should differ from individuals within mainland populations. Specifically, island populations are predicted to contain individuals that are larger, less aggressive, more sociable, and that invest more in their offspring. To date, tests of the ISH have mainly compared oceanic islands to continental sites, and rarely smaller spatial scales such as inland watersheds. Here, using a novel set of genome-wide SNP markers in wild deer mice (Peromyscus maniculatus) we conducted a genomic assessment of predictions underlying the ISH in an inland riverine island system: analysing island-mainland population structure, and quantifying heritability of phenotypes thought to underlie the ISH. We found clear genomic differentiation between the island and mainland populations and moderate to high marker-based heritability estimates for overall variation in traits previously found to differ in line with the ISH between mainland and island locations. F_ST outlier analyses highlighted 12 loci associated with differentiation between mainland and island populations. Together these results suggest that the island populations examined are on independent evolutionary trajectories, the traits considered have a genetic basis (rather than phenotypic variation being solely due to phenotypic plasticity). Coupled with the previous results showing significant phenotypic differentiation between the island and mainland groups in this system, this study suggests that the ISH can hold even on a small spatial scale.

A multi-omics study of circulating phospholipid markers of blood pressure

High-throughput techniques allow us to measure a wide-range of phospholipids which can provide insight into the mechanisms of hypertension. We aimed to conduct an in-depth multi-omics study of various phospholipids with systolic blood pressure (SBP) and diastolic blood pressure (DBP). The associations of blood pressure and 151 plasma phospholipids measured by electrospray ionization tandem mass spectrometry were performed by linear regression in five European cohorts (n = 2786 in discovery and n = 1185 in replication). We further explored the blood pressure-related phospholipids in Erasmus Rucphen Family (ERF) study by associating them with multiple cardiometabolic traits (linear regression) and predicting incident hypertension (Cox regression). Mendelian Randomization (MR) and phenome-wide association study (Phewas) were also explored to further investigate these association results. We identified six phosphatidylethanolamines (PE 38:3, PE 38:4, PE 38:6, PE 40:4, PE 40:5 and PE 40:6) and two phosphatidylcholines (PC 32:1 and PC 40:5) which together predicted incident hypertension with an area under the ROC curve (AUC) of 0.61. The identified eight phospholipids are strongly associated with triglycerides, obesity related traits (e.g. waist, waist-hip ratio, total fat percentage, body mass index, lipid-lowering medication, and leptin), diabetes related traits (e.g. glucose, insulin resistance and insulin) and prevalent type 2 diabetes. The genetic determinants of these phospholipids also associated with many lipoproteins, heart rate, pulse rate and blood cell counts. No significant association was identified by bi-directional MR approach. We identified eight blood pressure-related circulating phospholipids that have a predictive value for incident hypertension. Our cross-omics analyses show that phospholipid metabolites in the circulation may yield insight into blood pressure regulation and raise a number of testable hypothesis for future research.

Age- and Sex-Related Topological Organization of Human Brain Functional Networks and Their Relationship to Cognition

Age and sex associated with changes in the functional brain network topology and cognition in large population of older adults have been poorly understood. We explored this question further by examining differences in 11 resting-state graph theory measures with respect to age, sex, and their relationships with cognitive performance in 17,127 United Kingdom Biobank participants (mean = 62.83 ± 7.41 years). Age was associated with an overall decrease in the effectiveness of network communication (i.e., integration) and loss of functional specialization (i.e., segregation) of specific brain regions. Sex differences were also observed, with women showing more efficient networks, which were less segregated than in men (FDR adjusted p < 0.05). The age-related changes were also more apparent in men than in women, which suggests that men may be more vulnerable to cognitive decline with age. Interestingly, while network segregation and strength of limbic network were only nominally associated with cognitive performance, the network measures collectively were significantly associated with cognition (FDR adjusted p ≤ 0.002). This may imply that individual measures may be inadequate to capture much of the variance in the neural activity or its output and need further refinement. The complexity of the organization of the functional brain may be shaped by the age and sex of an individual, which ultimately may influence the cognitive performance of older adults. Age and sex stratification may be used to inform clinical neuroscience research to identify older adults at risk of cognitive dysfunction.

Novel genetic variants associated with brain functional networks in 18,445 adults from the UK Biobank

Here, we investigated the genetics of weighted functional brain network graph theory measures from 18,445 participants of the UK Biobank (44–80 years). The eighteen measures studied showed low heritability (mean h²_SNP = 0.12) and were highly genetically correlated. One genome-wide significant locus was associated with strength of somatomotor and limbic networks. These intergenic variants were located near the PAX8 gene on chromosome 2. Gene-based analyses identified five significantly associated genes for five of the network measures, which have been implicated in sleep duration, neuronal differentiation/development, cancer, and susceptibility to neurodegenerative diseases. Further analysis found that somatomotor network strength was phenotypically associated with sleep duration and insomnia. Single nucleotide polymorphism (SNP) and gene level associations with functional network measures were identified, which may help uncover novel biological pathways relevant to human brain functional network integrity and related disorders that affect it.

Effects of Estradiol on Immunoglobulin G Glycosylation: Mapping of the Downstream Signaling Mechanism

Glycans attached to immunoglobulin G (IgG) directly affect this antibody effector functions and regulate inflammation at several levels. The composition of IgG glycome changes significantly with age. In women, the most notable change coincides with the perimenopausal period. Aiming to investigate the effect of estrogen on IgG glycosylation, we analysed IgG and total serum glycomes in 36 healthy premenopausal women enrolled in a randomized controlled trial of the gonadotropin-releasing hormone analogue (GnRHAG) leuprolide acetate to lower gonadal steroids to postmenopausal levels and then randomized to transdermal placebo or estradiol (E2) patch. The suppression of gonadal hormones induced significant changes in the IgG glycome, while E2 supplementation was sufficient to prevent changes. The observed glycan changes suggest that depletion of E2 primarily affects B cell glycosylation, while liver glycosylation stays mostly unchanged. To determine whether previously identified IgG GWAS hits RUNX1, RUNX3, SPINK4, and ELL2 are involved in downstream signaling mechanisms, linking E2 with IgG glycosylation, we used the FreeStyle 293-F transient system expressing IgG antibodies with stably integrated CRISPR/dCas9 expression cassettes for gene up- and downregulation. RUNX3 and SPINK4 upregulation using dCas9-VPR resulted in a decreased IgG galactosylation and, in the case of RUNX3, a concomitant increase in IgG agalactosylation.

Genome-Wide Association Study-Based Identification of SNPs and Haplotypes Associated With Goose Reproductive Performance and Egg Quality

Geese are one of the most economically important waterfowl. However, the low reproductive performance and egg quality of geese hinder the development of the goose industry. The identification and application of genetic markers may improve the accuracy of beneficial trait selection. To identify the genetic markers associated with goose reproductive performance and egg quality traits, we performed a genome-wide association study (GWAS) for body weight at birth (BBW), the number of eggs at 48 weeks of age (EN48), the number of eggs at 60 weeks of age (EN60) and egg yolk color (EYC). The GWAS acquired 2.896 Tb of raw sequencing data with an average depth of 12.44× and identified 9,279,339 SNPs. The results of GWAS showed that 26 SNPs were significantly associated with BBW, EN48, EN60, and EYC. Moreover, five of these SNPs significantly associated with EN48 and EN60 were in a haplotype block on chromosome 35 from 4,512,855 to 4,541,709 bp, oriented to TMEM161A and another five SNPs significantly correlated to EYC were constructed in haplotype block on chromosome 5 from 21,069,009 to 21,363,580, which annotated by TMEM161A, CALCR, TFPI2, and GLP1R. Those genes were enriched in epidermal growth factor-activated receptor activity, regulation of epidermal growth factor receptor signaling pathway. The SNPs, haplotype markers, and candidate genes identified in this study can be used to improve the accuracy of marker-assisted selection for the reproductive performance and egg quality traits of geese. In addition, the candidate genes significantly associated with these traits may provide a foundation for better understanding the mechanisms underlying reproduction and egg quality in geese.

Autozygosity influences cardiometabolic disease-associated traits in the AWI-Gen sub-Saharan African study

The analysis of the effects of autozygosity, measured as the change of the mean value of a trait among offspring of genetic relatives, reveals the existence of directional dominance or overdominance. In this study we detect evidence of the effect of autozygosity in 4 out of 13 cardiometabolic disease-associated traits using data from more than 10,000 sub-Saharan African individuals recruited from Ghana, Burkina Faso, Kenya and South Africa. The effect of autozygosity on these phenotypes is found to be sex-related, with inbreeding having a significant decreasing effect in men but a significant increasing effect in women for several traits (body mass index, subcutaneous adipose tissue, low-density lipoproteins and total cholesterol levels). Overall, the effect of inbreeding depression is more intense in men. Differential effects of inbreeding depression are also observed between study sites with different night-light intensity used as proxy for urban development. These results suggest a directional dominant genetic component mediated by environmental interactions and sex-specific differences in genetic architecture for these traits in the Africa Wits-INDEPTH partnership for Genomic Studies (AWI-Gen) cohort.

Genome-wide associations of human gut microbiome variation and implications for causal inference analyses

Recent population-based1-4 and clinical studies5 have identified a range of factors associated with human gut microbiome variation. Murine quantitative trait loci6, human twin studies7 and microbiome genome-wide association studies1,3,8-12 have provided evidence for genetic contributions to microbiome composition. Despite this, there is still poor overlap in genetic association across human studies. Using appropriate taxon-specific models, along with support from independent cohorts, we show an association between human host genotype and gut microbiome variation. We also suggest that interpretation of applied analyses using genetic associations is complicated by the probable overlap between genetic contributions and heritable components of host environment. Using faecal 16S ribosomal RNA gene sequences and host genotype data from the Flemish Gut Flora Project (n = 2,223) and two German cohorts (FoCus, n = 950; PopGen, n = 717), we identify genetic associations involving multiple microbial traits. Two of these associations achieved a study-level threshold of P = 1.57 × 10-10; an association between Ruminococcus and rs150018970 near RAPGEF1 on chromosome 9, and between Coprococcus and rs561177583 within LINC01787 on chromosome 1. Exploratory analyses were undertaken using 11 other genome-wide associations with strong evidence for association (P < 2.5 × 10-8) and a previously reported signal of association between rs4988235 (MCM6/LCT) and Bifidobacterium. Across these 14 single-nucleotide polymorphisms there was evidence of signal overlap with other genome-wide association studies, including those for age at menarche and cardiometabolic traits. Mendelian randomization analysis was able to estimate associations between microbial traits and disease (including Bifidobacterium and body composition); however, in the absence of clear microbiome-driven effects, caution is needed in interpretation. Overall, this work marks a growing catalogue of genetic associations that will provide insight into the contribution of host genotype to gut microbiome. Despite this, the uncertain origin of association signals will likely complicate future work looking to dissect function or use associations for causal inference analysis.

Dysregulated Antibody, Natural Killer Cell and Immune Mediator Profiles in Autoimmune Thyroid Diseases

The pathogenesis of autoimmune thyroid diseases (AITD) is poorly understood and the association between different immune features and the germline variants involved in AITD are yet unclear. We previously observed systemic depletion of IgG core fucosylation and antennary α1,2 fucosylation in peripheral blood mononuclear cells in AITD, correlated with anti-thyroid peroxidase antibody (TPOAb) levels. Fucose depletion is known to potentiate strong antibody-mediated NK cell activation and enhanced target antigen-expressing cell killing. In autoimmunity, this may translate to autoantibody-mediated immune cell recruitment and attack of self-antigen expressing normal tissues. Hence, we investigated the crosstalk between immune cell traits, secreted proteins, genetic variants and the glycosylation patterns of serum IgG, in a multi-omic and cross-sectional study of 622 individuals from the TwinsUK cohort, 172 of whom were diagnosed with AITD. We observed associations between two genetic variants (rs505922 and rs687621), AITD status, the secretion of Desmoglein-2 protein, and the profile of two IgG N-glycan traits in AITD, but further studies need to be performed to better understand their crosstalk in AITD. On the other side, enhanced afucosylated IgG was positively associated with activatory CD335- CD314+ CD158b+ NK cell subsets. Increased levels of the apoptosis and inflammation markers Caspase-2 and Interleukin-1α positively associated with AITD. Two genetic variants associated with AITD, rs1521 and rs3094228, were also associated with altered expression of the thyrocyte-expressed ligands known to recognize the NK cell immunoreceptors CD314 and CD158b. Our analyses reveal a combination of heightened Fc-active IgG antibodies, effector cells, cytokines and apoptotic signals in AITD, and AITD genetic variants associated with altered expression of thyrocyte-expressed ligands to NK cell immunoreceptors. Together, TPOAb responses, dysregulated immune features, germline variants associated with immunoactivity profiles, are consistent with a positive autoreactive antibody-dependent NK cell-mediated immune response likely drawn to the thyroid gland in AITD.

Pigment Intensity in Dogs is Associated with a Copy Number Variant Upstream of KITLG

Dogs exhibit a wide variety of coat color types, and many genes have been identified that control pigment production, appearance, and distribution. Some breeds, such as the Nova Scotia Duck Tolling Retriever (NSDTR), exhibit variation in pheomelanin pigment intensity that is not explained by known genetic variants. A genome-wide association study comparing light red to dark red in the NSDTR identified a significantly associated region on canine chromosome 15 (CFA 15:23 Mb–38 Mb). Coverage analysis of whole genome sequence data from eight dogs identified a 6 kb copy number variant (CNV) 152 kb upstream of KITLG. Genotyping with digital droplet PCR (ddPCR) confirmed a significant association between an increased copy number with the dark-red coat color in NSDTR (p = 6.1 × 10⁻⁷). The copy number of the CNV was also significantly associated with coat color variation in both eumelanin and pheomelanin-based Poodles (p = 1.5 × 10⁻⁸, 4.0 × 10⁻⁹) and across other breeds. Moreover, the copy number correlated with pigment intensity along the hair shaft in both pheomelanin and eumelanin coats. KITLG plays an important role in melanogenesis, and variants upstream of KITLG have been associated with coat color variation in mice as well as hair color in humans consistent with its role in the domestic dog.

Heritability of Horn Size in Thinhorn Sheep

Understanding the genetic basis of fitness-related trait variation has long been of great interest to evolutionary biologists. Secondary sexual characteristics, such as horns in bovids, are particularly intriguing since they can be potentially affected by both natural and sexual selection. Until recently, however, the study of fitness-related quantitative trait variation in wild species has been hampered by a lack of genomic resources, pedigree, and/or phenotype data. Recent innovations in genomic technologies have enabled wildlife researchers to perform marker-based relatedness estimation and acquire adequate loci density, enabling both the “top-down” approach of quantitative genetics and the “bottom-up” approach of association studies to describe the genetic basis of fitness-related traits. Here we combine a cross species application of the OvineHD BeadChip and horn measurements (horn length, base circumference, and volume) from harvested thinhorn sheep to examine the heritability and to perform a genome-wide single-nucleotide polymorphism association study of horn size in the species. Thinhorn sheep are mountain ungulates that reside in the mountainous regions of northwestern North America. Thinhorn sheep males grow massive horns that determine the social rank and mating success. We found horn length, base circumference, and volume to be moderately heritable and two loci to be suggestively associated with horn length.

Genomic Regions Associated with IgE Levels against Culicoides spp. Antigens in Three Horse Breeds

Insect bite hypersensitivity (IBH), which is a cutaneous allergic reaction to antigens from Culicoides spp., is the most prevalent skin disorder in horses. Misdiagnosis is possible, as IBH is usually diagnosed based on clinical signs. Our study is the first to employ IgE levels against several recombinant Culicoides spp. allergens as an objective, independent, and quantitative phenotype to improve the power to detect genetic variants that underlie IBH. Genotypes of 200 Shetland ponies, 127 Icelandic horses, and 223 Belgian Warmblood horses were analyzed while using a mixed model approach. No single-nucleotide polymorphism (SNP) passed the Bonferroni corrected significance threshold, but several regions were identified within and across breeds, which confirmed previously identified regions of interest and, in addition, identifying new regions of interest. Allergen-specific IgE levels are a continuous and objective phenotype that allow for more powerful analyses when compared to a case-control set-up, as more significant associations were obtained. However, the use of a higher density array seems necessary to fully employ the use of IgE levels as a phenotype. While these results still require validation in a large independent dataset, the use of allergen-specific IgE levels showed value as an objective and continuous phenotype that can deepen our understanding of the biology underlying IBH.

Defining the genetic control of human blood plasma N-glycome using genome-wide association study

Glycosylation is a common post-translational modification of proteins. Glycosylation is associated with a number of human diseases. Defining genetic factors altering glycosylation may provide a basis for novel approaches to diagnostic and pharmaceutical applications. Here we report a genome-wide association study of the human blood plasma N-glycome composition in up to 3811 people measured by Ultra Performance Liquid Chromatography (UPLC) technology. Starting with the 36 original traits measured by UPLC, we computed an additional 77 derived traits leading to a total of 113 glycan traits. We studied associations between these traits and genetic polymorphisms located on human autosomes. We discovered and replicated 12 loci. This allowed us to demonstrate an overlap in genetic control between total plasma protein and IgG glycosylation. The majority of revealed loci contained genes that encode enzymes directly involved in glycosylation (FUT3/FUT6, FUT8, B3GAT1, ST6GAL1, B4GALT1, ST3GAL4, MGAT3 and MGAT5) and a known regulator of plasma protein fucosylation (HNF1A). However, we also found loci that could possibly reflect other more complex aspects of glycosylation process. Functional genomic annotation suggested the role of several genes including DERL3, CHCHD10, TMEM121, IGH and IKZF1. The hypotheses we generated may serve as a starting point for further functional studies in this research area.

Protocols, Methods, and Tools for Genome-Wide Association Studies (GWAS) of Dental Traits

Oral health and disease are known to be influenced by complex interactions between environmental (e.g., social and behavioral) factors and innate susceptibility. Although the exact contribution of genomics and other layers of "omics" to oral health is an area of active research, it is well established that the susceptibility to dental caries, periodontal disease, and other oral and craniofacial traits is substantially influenced by the human genome. A comprehensive understanding of these genomic factors is necessary for the realization of precision medicine in the oral health domain. To aid in this direction, the advent and increasing affordability of high-throughput genotyping has enabled the simultaneous interrogation of millions of genetic polymorphisms for association with oral and craniofacial traits. Specifically, genome-wide association studies (GWAS) of dental caries and periodontal disease have provided initial insights into novel loci and biological processes plausibly implicated in these two common, complex, biofilm-mediated diseases. This paper presents a summary of protocols, methods, tools, and pipelines for the conduct of GWAS of dental caries, periodontal disease, and related traits. The protocol begins with the consideration of different traits for both diseases and outlines procedures for genotyping, quality control, adjustment for population stratification, heritability and association analyses, annotation, reporting, and interpretation. Methods and tools available for GWAS are being constantly updated and improved; with this in mind, the presented approaches have been successfully applied in numerous GWAS and meta-analyses among tens of thousands of individuals, including dental traits such as dental caries and periodontal disease. As such, they can serve as a guide or template for future genomic investigations of these and other traits.

Genome-epigenome interactions associated with Myalgic Encephalomyelitis/Chronic Fatigue Syndrome

Myalgic Encephalomyelitis/Chronic Fatigue Syndrome (ME/CFS) is a complex disease of unknown etiology. Multiple studies point to disruptions in immune functioning in ME/CFS patients as well as specific genetic polymorphisms and alterations of the DNA methylome in lymphocytes. However, potential interactions between DNA methylation and genetic background in relation to ME/CFS have not been examined. In this study we explored this association by characterizing the epigenetic (~480 thousand CpG loci) and genetic (~4.3 million SNPs) variation between cohorts of ME/CFS patients and healthy controls. We found significant associations of DNA methylation states in T-lymphocytes at several CpG loci and regions with ME/CFS phenotype. These methylation anomalies are in close proximity to genes involved with immune function and cellular metabolism. Finally, we found significant correlations of genotypes with methylation modifications associated with ME/CFS. The findings from this study highlight the role of epigenetic and genetic interactions in complex diseases, and suggest several genetic and epigenetic elements potentially involved in the mechanisms of disease in ME/CFS.

Age-specific genome-wide association study in glioblastoma identifies increased proportion of 'lower grade glioma'-like features associated with younger age

Glioblastoma (GBM) is the most common malignant brain tumor in the United States. Incidence of GBM increases with age, and younger age-at-diagnosis is significantly associated with improved prognosis. While the relationship between candidate GBM risk SNPs and age-at-diagnosis has been explored, genome-wide association studies (GWAS) have not previously been stratified by age. Potential age-specific genetic effects were assessed in autosomal SNPs for GBM patients using data from four previous GWAS. Using age distribution tertiles (18-53, 54-64, 65+) datasets were analyzed using age-stratified logistic regression to generate p values, odds ratios (OR), and 95% confidence intervals (95%CI), and then combined using meta-analysis. There were 4,512 total GBM cases, and 10,582 controls used for analysis. Significant associations were detected at two previously identified SNPs in 7p11.2 (rs723527 [p54-63 = 1.50x10-9 , OR54-63 = 1.28, 95%CI54-63 = 1.18-1.39; p64+ = 2.14x10-11 , OR64+ = 1.32, 95%CI64+ = 1.21-1.43] and rs11979158 [p54-63 = 6.13x10-8 , OR54-63 = 1.35, 95%CI54-63 = 1.21-1.50; p64+ = 2.18x10-10 , OR64+ = 1.42, 95%CI64+ = 1.27-1.58]) but only in persons >54. There was also a significant association at the previously identified lower grade glioma (LGG) risk locus at 8q24.21 (rs55705857) in persons ages 18-53 (p18-53 = 9.30 × 10-11 , OR18-53 = 1.76, 95%CI18-53 = 1.49-2.10). Within The Cancer Genome Atlas (TCGA) there was higher prevalence of 'LGG'-like tumor characteristics in GBM samples in those 18-53, with IDH1/2 mutation frequency of 15%, as compared to 2.1% [54-63] and 0.8% [64+] (p = 0.0005). Age-specific differences in cancer susceptibility can provide important clues to etiology. The association of a SNP known to confer risk for IDH1/2 mutant glioma and higher prevalence of IDH1/2 mutation within younger individuals 18-53 suggests that more younger individuals may present initially with 'secondary glioblastoma.'

Adrenocortical Challenge Response and Genomic Analyses in Scottish Terriers With Increased Alkaline Phosphate Activity

Scottish terriers (ST) frequently have increased serum alkaline phosphatase (ALP) of the steroid isoform. Many of these also have high serum concentrations of adrenal sex steroids. The study's objective was to determine the cause of increased sex steroids in ST with increased ALP. Adrenal gland suppression and stimulation were compared by low dose dexamethasone (LDDS), human chorionic gonadotropin (HCG) and adrenocorticotropic hormone (ACTH) response tests. Resting plasma pituitary hormones were measured. Steroidogenesis-related mRNA expression was evaluated in six ST with increased ALP, eight dogs of other breeds with pituitary-dependent hyperadrenocorticism (HAC), and seven normal dogs. The genome-wide association of single nucleotide polymorphisms (SNP) with ALP activity was evaluated in 168 ST. ALP (reference interval 8–70 U/L) was high in all ST (1,054 U/L) and HAC (985 U/L) dogs. All HAC dogs and 2/8 ST had increased cortisol post-ACTH administration. All ST and 2/7 Normal dogs had increased sex steroids post-ACTH. ST and Normal dogs had similar post-challenge adrenal steroid profiles following LDDS and HCG. Surprisingly, mRNA of hydroxysteroid 17-beta dehydrogenase 2 (HSD17B2) was lower in ST and Normal dogs than HAC. HSD17B2 facilities metabolism of sex steroids. A SNP region was identified on chromosome 5 in proximity to HSD17B2 that correlated with increased serum ALP. ST in this study with increased ALP had a normal pituitary-adrenal axis in relationship to glucocorticoids and luteinizing hormone. We speculate the identified SNP and HSD17B2 gene may have a role in the pathogenesis of elevated sex steroids and ALP in ST.

Genetic variation within endolysosomal system is associated with late-onset Alzheimer’s disease

Late-onset Alzheimer’s disease is the most common dementia type, yet no treatment exists to stop the neurodegeneration. Evidence from monogenic lysosomal diseases, neuronal pathology and experimental models suggest that autophagic and endolysosomal dysfunction may contribute to neurodegeneration by disrupting the degradation of potentially neurotoxic molecules such as amyloid-β and tau. However, it is uncertain how well the evidence from rare disorders and experimental models capture causal processes in common forms of dementia, including late-onset Alzheimer’s disease. For this reason, we set out to investigate if autophagic and endolysosomal genes were enriched for genetic variants that convey increased risk of Alzheimer’s disease; such a finding would provide population-based support for the endolysosomal hypothesis of neurodegeneration. We quantified the collective genetic associations between the endolysosomal system and Alzheimer’s disease in three genome-wide associations studies (combined n = 62 415). We used the Mergeomics pathway enrichment algorithm that incorporates permutations of the full hierarchical cascade of SNP-gene-pathway to estimate enrichment. We used a previously published collection of 891 autophagic and endolysosomal genes (denoted as AphagEndoLyso, and derived from the Lysoplex sequencing platform) as a proxy for cellular processes related to autophagy, endocytosis and lysosomal function. We also investigated a subset of 142 genes of the 891 that have been implicated in Mendelian diseases (MenDisLyso). We found that both gene sets were enriched for genetic Alzheimer’s associations: an enrichment score 3.67 standard deviations from the null model (P = 0.00012) was detected for AphagEndoLyso, and a score 3.36 standard deviations from the null model (P = 0.00039) was detected for MenDisLyso. The high enrichment score was specific to the AphagEndoLyso gene set (stronger than 99.7% of other tested pathways) and to Alzheimer’s disease (stronger than all other tested diseases). The APOE locus explained most of the MenDisLyso signal (1.16 standard deviations after APOE removal, P = 0.12), but the AphagEndoLyso signal was less affected (3.35 standard deviations after APOE removal, P = 0.00040). Additional sensitivity analyses further indicated that the AphagEndoLyso Gene Set contained an aggregate genetic association that comprised a combination of subtle genetic signals in multiple genes. We also observed an enrichment of Parkinson’s disease signals for MenDisLyso (3.25 standard deviations) and for AphagEndoLyso (3.95 standard deviations from the null model), and a brain-specific pattern of gene expression for AphagEndoLyso in the Gene Tissue Expression Project dataset. These results provide evidence that a diffuse aggregation of genetic perturbations to the autophagy and endolysosomal system may mediate late-onset Alzheimer’s risk in human populations.

Sex-specific glioma genome-wide association study identifies new risk locus at 3p21.31 in females, and finds sex-differences in risk at 8q24.21

Incidence of glioma is approximately 50% higher in males. Previous analyses have examined exposures related to sex hormones in women as potential protective factors for these tumors, with inconsistent results. Previous glioma genome-wide association studies (GWAS) have not stratified by sex. Potential sex-specific genetic effects were assessed in autosomal SNPs and sex chromosome variants for all glioma, GBM and non-GBM patients using data from four previous glioma GWAS. Datasets were analyzed using sex-stratified logistic regression models and combined using meta-analysis. There were 4,831 male cases, 5,216 male controls, 3,206 female cases and 5,470 female controls. A significant association was detected at rs11979158 (7p11.2) in males only. Association at rs55705857 (8q24.21) was stronger in females than in males. A large region on 3p21.31 was identified with significant association in females only. The identified differences in effect of risk variants do not fully explain the observed incidence difference in glioma by sex.

Using genome wide association studies to identify common QTL regions in three different genetic backgrounds based on Iberian pig breed

One of the major limitation for the application of QTL results in pig breeding and QTN identification has been the limited number of QTL effects validated in different animal material. The aim of the current work was to validate QTL regions through joint and specific genome wide association and haplotype analyses for growth, fatness and premier cut weights in three different genetic backgrounds, backcrosses based on Iberian pigs, which has a major role in the analysis due to its high productive relevance. The results revealed nine common QTL regions, three segregating in all three backcrosses on SSC1, 0–3 Mb, for body weight, on SSC2, 3–9 Mb, for loin bone-in weight, and on SSC7, 3 Mb, for shoulder weight, and six segregating in two of the three backcrosses, on SSC2, SSC4, SSC6 and SSC10 for backfat thickness, shoulder and ham weights. Besides, 18 QTL regions were specifically identified in one of the three backcrosses, five identified only in BC_LD, seven in BC_DU and six in BC_PI. Beyond identifying and validating QTL, candidate genes and gene variants within the most interesting regions have been explored using functional annotation, gene expression data and SNP identification from RNA-Seq data. The results allowed us to propose a promising list of candidate mutations, those identified in PDE10A, DHCR7, MFN2 and CCNY genes located within the common QTL regions and those identified near ssc-mir-103-1 considered PANK3 regulators to be further analysed.

Genomic analysis of morphometric traits in bighorn sheep using the Ovine Infinium® HD SNP BeadChip

Elucidating the genetic basis of fitness-related traits is a major goal of molecular ecology. Traits subject to sexual selection are particularly interesting, as non-random mate choice should deplete genetic variation and thereby their evolutionary benefits. We examined the genetic basis of three sexually selected morphometric traits in bighorn sheep (Ovis canadensis): horn length, horn base circumference, and body mass. These traits are of specific concern in bighorn sheep as artificial selection through trophy hunting opposes sexual selection. Specifically, horn size determines trophy status and, in most North American jurisdictions, if an individual can be legally harvested. Using between 7,994–9,552 phenotypic measures from the long-term individual-based study at Ram Mountain (Alberta, Canada), we first showed that all three traits are heritable (h² = 0.15–0.23). We then conducted a genome-wide association study (GWAS) utilizing a set of 3,777 SNPs typed in 76 individuals using the Ovine Infinium^® HD SNP BeadChip. We found suggestive association for body mass at a single locus (OAR9_91647990). The absence of strong associations with SNPs suggests that the traits are likely polygenic. These results represent a step forward for characterizing the genetic architecture of fitness related traits in sexually dimorphic ungulates.

Multivariate discovery and replication of five novel loci associated with Immunoglobulin G N-glycosylation

Joint modeling of a number of phenotypes using multivariate methods has often been neglected in genome-wide association studies and if used, replication has not been sought. Modern omics technologies allow characterization of functional phenomena using a large number of related phenotype measures, which can benefit from such joint analysis. Here, we report a multivariate genome-wide association studies of 23 immunoglobulin G (IgG) N-glycosylation phenotypes. In the discovery cohort, our multi-phenotype method uncovers ten genome-wide significant loci, of which five are novel (IGH, ELL2, HLA-B-C, AZI1, FUT6-FUT3). We convincingly replicate all novel loci via multivariate tests. We show that IgG N-glycosylation loci are strongly enriched for genes expressed in the immune system, in particular antibody-producing cells and B lymphocytes. We empirically demonstrate the efficacy of multivariate methods to discover novel, reproducible pleiotropic effects.Multivariate analysis methods can uncover the relationship between phenotypic measures characterised by modern omic techniques. Here the authors conduct a multivariate GWAS on IgG N-glycosylation phenotypes and identify 5 novel loci enriched in immune system genes.

Limited influence of germline genetic variation on all-cause mortality in women with early onset breast cancer: evidence from gene-based tests, single-marker regression, and whole-genome prediction

PURPOSE

Women diagnosed with breast cancer have heterogeneous survival outcomes that cannot be fully explained by known prognostic factors, and germline variation is a plausible but unconfirmed risk factor.

METHODS

We used three approaches to test the hypothesis that germline variation drives some differences in survival: mortality loci identification, tumor aggressiveness loci identification, and whole-genome prediction. The 2954 study participants were women diagnosed with breast cancer before age 50, with a median follow-up of 15 years who were genotyped on an exome array. We first searched for loci in gene regions that were associated with all-cause mortality. We next searched for loci in gene regions associated with five histopathological characteristics related to tumor aggressiveness. Last, we also predicted 10-year all-cause mortality on a subset of 1903 participants (3,245,343 variants after imputation) using whole-genome prediction methods.

RESULTS

No risk loci for mortality or tumor aggressiveness were identified. This null result persisted when restricting to women with estrogen receptor-positive tumors, when examining suggestive loci in an independent study, and when restricting to previously published risk loci. Additionally, the whole-genome prediction model also found no evidence to support an association.

CONCLUSION

Despite multiple complementary approaches, our study found no evidence that mortality in women with early onset breast cancer is influenced by germline variation.

Association study reveals novel risk loci for sporadic inclusion body myositis

BACKGROUND AND PURPOSE

The aim was to identify potential genetic risk factors associated with sporadic inclusion body myositis (sIBM).

METHODS

An association based case-control approach was utilized on whole exome sequencing data of 30 Finnish sIBM patients and a control cohort (n = 193). A separate Italian cohort of sIBM patients (n = 12) was used for evaluation of the results.

RESULTS

Seven single nucleotide polymorphisms were identified in five genes that have a considerably higher observed frequency in Finnish sIBM patients compared to the control population, and the previous association of the genetic human leukocyte antigen region was confirmed.

CONCLUSIONS

All seven identified variants could individually or in combination increase the susceptibility for sIBM.

Estimating Disequilibrium Coefficients

Gametic phase disequilibrium is the nonrandom association of alleles within gametes. Linkage disequilibrium (LD) describes the special case of deviation from independence between alleles at two linked genetic loci. Estimation of allelic LD requires knowledge of haplotypes. Genotype-based LD measures dispense with the haplotype estimation step and avoid bias in LD estimation. In this chapter, the most important measures for allelic and genotypic LD are introduced. The use of software packages for LD estimation is illustrated.

Deciphering the regulation of porcine genes influencing growth, fatness and yield-related traits through genetical genomics

Genetical genomics approaches aim at identifying quantitative trait loci for molecular traits, also known as intermediate phenotypes, such as gene expression, that could link variation in genetic information to physiological traits. In the current study, an expression GWAS has been carried out on an experimental Iberian × Landrace backcross in order to identify the genomic regions regulating the gene expression of those genes whose expression is correlated with growth, fat deposition, and premium cut yield measures in pig. The analyses were conducted exploiting Porcine 60K SNP BeadChip genotypes and Porcine Expression Microarray data hybridized on mRNA from Longissimus dorsi muscle. In order to focus the analysis on productive traits and reduce the number of analyses, only those probesets whose expression showed significant correlation with at least one of the seven phenotypes of interest were selected for the eGWAS. A total of 63 eQTL regions were identified with effects on 36 different transcripts. Those eQTLs overlapping with phenotypic QTLs on SSC4, SSC9, SSC13, and SSC17 chromosomes previously detected in the same animal material were further analyzed. Moreover, candidate genes and SNPs were analyzed. Among the most promising results, a long non-coding RNA, ALDBSSCG0000001928, was identified, whose expression is correlated with premium cut yield. Association analysis and in silico sequence domain annotation support TXNRD3 polymorphisms as candidate to regulate ALDBSSCG0000001928 expression, which can be involved in the transcriptional regulation of surrounding genes, affecting productive and meat quality traits.

Non-additive genome-wide association scan reveals a new gene associated with habitual coffee consumption

Coffee is one of the most consumed beverages world-wide and one of the primary sources of caffeine intake. Given its important health and economic impact, the underlying genetics of its consumption has been widely studied. Despite these efforts, much has still to be uncovered. In particular, the use of non-additive genetic models may uncover new information about the genetic variants driving coffee consumption. We have conducted a genome-wide association study in two Italian populations using additive, recessive and dominant models for analysis. This has uncovered a significant association in the PDSS2 gene under the recessive model that has been replicated in an independent cohort from the Netherlands (ERF). The identified gene has been shown to negatively regulate the expression of the caffeine metabolism genes and can thus be linked to coffee consumption. Further bioinformatics analysis of eQTL and histone marks from Roadmap data has evidenced a possible role of the identified SNPs in regulating PDSS2 gene expression through enhancers present in its intron. Our results highlight a novel gene which regulates coffee consumption by regulating the expression of the genes linked to caffeine metabolism. Further studies will be needed to clarify the biological mechanism which links PDSS2 and coffee consumption.

The GenABEL Project for statistical genomics

Development of free/libre open source software is usually done by a community of people with an interest in the tool. For scientific software, however, this is less often the case. Most scientific software is written by only a few authors, often a student working on a thesis. Once the paper describing the tool has been published, the tool is no longer developed further and is left to its own device. Here we describe the broad, multidisciplinary community we formed around a set of tools for statistical genomics. The GenABEL project for statistical omics actively promotes open interdisciplinary development of statistical methodology and its implementation in efficient and user-friendly software under an open source licence. The software tools developed withing the project collectively make up the GenABEL suite, which currently consists of eleven tools. The open framework of the project actively encourages involvement of the community in all stages, from formulation of methodological ideas to application of software to specific data sets. A web forum is used to channel user questions and discussions, further promoting the use of the GenABEL suite. Developer discussions take place on a dedicated mailing list, and development is further supported by robust development practices including use of public version control, code review and continuous integration. Use of this open science model attracts contributions from users and developers outside the “core team”, facilitating agile statistical omics methodology development and fast dissemination.