Genome-wide Association Loci Research Articles

SLC39A8.p.(Ala391Thr) is associated with poorer cognitive ability: a cross-sectional study of schizophrenia and the general UK population.

The missense SNP NC_000004.12:g.102267552C>T (SLC39A8.p.(Ala391Thr), rs13107325) in SLC39A8, which encodes a zinc transporter, has been linked to schizophrenia and is the likely causal variant for one of the genome-wide association loci associated with the disorder. We tested whether the schizophrenia-risk allele at p.(Ala391Thr) was associated with schizophrenia-related phenotypes, including positive, negative, and disorganised symptoms, cognitive ability, educational attainment, and age of psychosis onset, within three schizophrenia cohorts (combined N=1,232) and, with equivalent phenotypes, in a sample of population controls (UK Biobank, N=355,069). We used regression analyses controlling for age, sex, and population stratification. Within the schizophrenia cohorts, after correction for multiple testing, p.(Ala391Thr) was not significantly associated with any schizophrenia-related phenotypes. In the unaffected participants from the UK Biobank, the schizophrenia-risk allele at p.(Ala391Thr) was associated with significantly poorer cognitive ability and fluid intelligence, a lower probability of obtaining GCSEs or a degree-level qualification, and fewer years in education. There was no association between p.(Ala391Thr) and self-reported psychotic experiences in this cohort. The schizophrenia-risk allele was associated with poorer cognitive ability, but not psychotic experiences, in a volunteer sample drawn from of the general population. To determine whether p.(Ala391Thr) is associated with cognitive phenotypes in people with schizophrenia, and to understand the role of p.(Ala391Thr) in the aetiology of cognitive impairment in schizophrenia, larger independent samples are required.

Genome-wide association and expression quantitative trait loci in cattle reveals common genes regulating mammalian fertility

Most genetic variants associated with fertility in mammals fall in non-coding regions of the genome and it is unclear how these variants affect fertility. Here we use genome-wide association summary statistics for Heifer puberty (pubertal or not at 600 days) from 27,707 Bos indicus, Bos taurus and crossbred cattle; multi-trait GWAS signals from 2119 indicine cattle for four fertility traits, including days to calving, age at first calving, pregnancy status, and foetus age in weeks (assessed by rectal palpation of the foetus); and expression quantitative trait locus for whole blood from 489 indicine cattle, to identify 87 putatively functional genes affecting cattle fertility. Our analysis reveals a significant overlap between the set of cattle and previously reported human fertility-related genes, impling the existence of a shared pool of genes that regulate fertility in mammals. These findings are crucial for developing approaches to improve fertility in cattle and potentially other mammals.

Discerning Endoscopic Severity of Inflammatory Bowel Disease by Scoping the Peripheral Blood Transcriptome

Background and AimsUlcerative colitis (UC) and Crohn’s disease (CD) are chronic inflammatory bowel diseases (IBD) with an incompletely understood etiology and pathogenesis. Identification of suitable drug targets and assessment of disease severity are crucial for optimal management. MethodsUsing RNA-sequencing, we investigated differential gene expression in peripheral blood samples from IBD-patients and non-inflamed controls, analysed pathway enrichment, and identified genes whose expression correlated with endoscopic disease severity. ResultsNeutrophil degranulation emerged as the most significant pathway across all IBD sample types. Signalling by interleukins was prominent in patients with active intestinal inflammation, but also enriched in CD and UC patients without intestinal inflammation. Nevertheless, genes correlated to endoscopic disease severity implicated the primary cilium in CD patients, and translation and focal adhesion in UC patients. Moreover, several of these genes were located in genome-wide associated loci linked to IBD, cholesterol levels, blood cell counts, and levels of markers assessing liver and kidney function. These genes also suggested connections to intestinal epithelial barrier dysfunction, contemporary IBD drug treatment, and new actionable drug targets. A large number of genes associated with endoscopic disease severity corresponded to non-coding RNAs. ConclusionsThis study revealed biological pathways associated with IBD disease state and endoscopic disease severity, thus providing insights into the underlying mechanisms of IBD pathogenesis as well as identifying potential biomarkers and therapies. Peripheral blood might constitute a suitable non-invasive diagnostic sample type, in which gene expression profiles might serve as indicators of ongoing mucosal inflammation, and thus guide personalized treatment decisions.

Large-scale exome array summary statistics resources for glycemic traits to aid effector gene prioritization.

Genome-wide association studies for glycemic traits have identified hundreds of loci associated with these biomarkers of glucose homeostasis. Despite this success, the challenge remains to link variant associations to genes, and underlying biological pathways. To identify coding variant associations which may pinpoint effector genes at both novel and previously established genome-wide association loci, we performed meta-analyses of exome-array studies for four glycemic traits: glycated hemoglobin (HbA1c, up to 144,060 participants), fasting glucose (FG, up to 129,665 participants), fasting insulin (FI, up to 104,140) and 2hr glucose post-oral glucose challenge (2hGlu, up to 57,878). In addition, we performed network and pathway analyses. Single-variant and gene-based association analyses identified coding variant associations at more than 60 genes, which when combined with other datasets may be useful to nominate effector genes. Network and pathway analyses identified pathways related to insulin secretion, zinc transport and fatty acid metabolism. HbA1c associations were strongly enriched in pathways related to blood cell biology. Our results provided novel glycemic trait associations and highlighted pathways implicated in glycemic regulation. Exome-array summary statistic results are being made available to the scientific community to enable further discoveries.

Expression of the readthrough transcript CiDRE in alveolar macrophages boosts SARS-CoV-2 susceptibility and promotes COVID-19 severity

Lung infection during severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) via the angiotensin-I-converting enzyme 2 (ACE2) receptor induces a cytokine storm. However, the precise mechanisms involved in severe COVID-19 pneumonia are unknown. Here, we showed that interleukin-10 (IL-10) induced the expression of ACE2 in normal alveolar macrophages, causing them to become vectors for SARS-CoV-2. The inhibition of this system in hamster models attenuated SARS-CoV-2 pathogenicity. Genome-wide association and quantitative trait locus analyses identified a IFNAR2-IL10RB readthrough transcript, COVID-19 infectivity-enhancing dual receptor (CiDRE), which was highly expressed in patients harboring COVID-19 risk variants at the IFNAR2 locus. We showed that CiDRE exerted synergistic effects via the IL-10-ACE2 axis in alveolar macrophages and functioned as a decoy receptor for type I interferons. Collectively, our data show that high IL-10 and CiDRE expression are potential risk factors for severe COVID-19. Thus, IL-10R and CiDRE inhibitors might be useful COVID-19 therapies.

Detection of the local adaptive and genome-wide associated loci in southeast Nigerian taro (Colocasia esculenta (L.) Schott) populations

BackgroundTaro has a long history of being consumed and remains orphan and on the hand Nigeria farmers. The role of farmer-driven artificial selection is not negligible to fit landraces to a particular ecological condition. Limited study has been conducted on genome-wide association and no study has been conducted on genome-environment association for clinal adaptation for taro. Therefore, the objective of this study was to detect loci that are associated with environmental variables and phenotype traits and forward input to breeders. The study used 92 geographical referred taro landraces collected from Southeast (SE) Nigeria.ResultsThe result indicates that SE Nigerian taro has untapped phenotype and genetic variability with low admixture. Redundancy analysis indicated that collinear explained SNP variation more than single climatic variable. Overall, the results indicated that no single method exclusively was able to capture population confounding effects better than the others for all six traits. Nevertheless, based on overall model performance, Blink seemed to provide slight advantage over other models and was selected for all subsequent assessment of genome-environment association (GEA) and genome-wide association study (GWAS) models. Genome scan and GEA identified local adapted loci and co-located genes. A total of nine SNP markers associated with environmental variables. Some of the SNP markers (such as S_101024366) co-located with genes which previously reported for climatic adaptation such as astringency, diaminopimelate decarboxylase and MYB transcription factor. Genome-wide association also identified 45, 40 and 34 significant SNP markers associated with studied traits in combined, year 1 and year 2 data sets, respectively. Out of these, five SNP markers (S1_18891752 S3_100795476, S1_100584471 S1_100896936 and S2_10058799) were consistent in two different data sets.ConclusionsThe findings from this study improve our understanding of the genetic control of adaptive and phenotypic traits in Nigerian taro. However, the study suggests further study on identification of local adaptive loci and GWAS through collection of more landraces throughout the country, and across different agro-ecologies.

Brain Proteome-Wide and Transcriptome-Wide Asso-ciation Studies, Bayesian Colocalization, and Mendelian Randomization Analyses Reveal Causal Genes of Parkinson's Disease.

How genome-wide associated loci confer risk for Parkinson's disease is unclear. We aim to reveal causal genes through effects on brain proteins to provide new pathogenesis insights for Parkinson's disease. Proteome-wide and transcriptome-wide associations were determined by functional summary-based imputation leveraging data from genome-wide association summary (56 306 Europeans, 1.4 million controls), brain proteomes (528 cases from 2 separate data sets), and transcriptome (452 cases), followed by Mendelian randomization, Bayesian colocalization, cell-type-specific and brain regional expression, and drug-gene interaction analyses. As a result, genetically regulated protein abundances of 11 genes were associated with Parkinson's disease. Five genes (CD38, GPNMB, TMEM175, RAB7L1, and HIP1R) were colocalized. Four genes (GPNMB, SEC23IP, CD38, and DGKQ) demonstrated Mendelian randomized correlations (p < 8.10 × 10-5). Higher GPNMB level (1.47, 1.28-1.68) and lower CD38 level (0.319, 0.24-0.43) were causally associated with higher risk of Parkinson's disease, consistent with transcriptomic evaluations. CD38 and GPNMB were preferentially enriched in astrocytes and oligodendrocyte precursor cells, respectively. And CD38 and GPNMB were suggested to be the targets of many oncological drugs from Drug-Gene Interaction database. In conclusion, utilizing multidimensional data, GPNMB and CD38 were prioritized as the causal genes of Parkinson's disease, crucial for mechanistic and therapeutic investigations.

13.2 Identification and Characterization of Genome-Wide Associated Loci in OCD

OCD is a heritable disorder, but recent genome-wide association studies (GWAS) had limited success to identify replicated OCD susceptibility loci at scale. On behalf of the working group in the Psychiatric Genomics Consortium focused on Tourette’s disorder and OCD, we report results from the largest OCD case-control sample to date.

Inverse and Concordant Mucosal Pathway Gene Expressions in Inflamed and Non-Inflamed Ulcerative Colitis Patients: Potential Relevance to Aetiology and Pathogenesis.

Ulcerative colitis (UC) arises from a complex interplay between host and environmental factors, but with a largely unsolved pathophysiology. The pathophysiology was outlined by RNA-sequencing of mucosal biopsies from non-inflamed and inflamed colon of UC patients (14 and 17, respectively), and from 27 patients without intestinal inflammation. Genes differentially expressed (DE), or present in enriched gene sets, were investigated using statistical text analysis of functional protein information. Compared with controls, inflamed and non-inflamed UC mucosa displayed 9360 and 52 DE genes, respectively. Seventy-three non-pseudogenes were DE relative to both gender and inflammation. Mitochondrial processes were downregulated in inflamed and upregulated in non-inflamed UC mucosa, whereas angiogenesis and endoplasmic reticulum (ER) stress were upregulated in both tissue states. Immune responses were upregulated in inflamed mucosa, whereas the non-inflamed UC mucosa presented both up- and downregulated gene sets. DE and enriched genes overlapped with genes present in inflammatory bowel disease genome-wide associated loci (p = 1.43 × 10−18), especially regarding immune responses, respiratory chain, angiogenesis, ER stress, and steroid hormone metabolism. Apart from confirming established pathophysiological mechanisms of immune cells, our study provides evidence for involvement of less described pathways (e.g., respiratory chain, ER stress, fatty-acid oxidation, steroid hormone metabolism and angiogenesis).

ZEB2 Shapes the Epigenetic Landscape of Atherosclerosis.

Smooth muscle cells (SMCs) transition into a number of different phenotypes during atherosclerosis, including those that resemble fibroblasts and chondrocytes, and make up the majority of cells in the atherosclerotic plaque. To better understand the epigenetic and transcriptional mechanisms that mediate these cell state changes, and how they relate to risk for coronary artery disease (CAD), we have investigated the causality and function of transcription factors at genome-wide associated loci. We used CRISPR-Cas 9 genome and epigenome editing to identify the causal gene and cells for a complex CAD genome-wide association study signal at 2q22.3. Single-cell epigenetic and transcriptomic profiling in murine models and human coronary artery smooth muscle cells were used to understand the cellular and molecular mechanism by which this CAD risk gene exerts its function. CRISPR-Cas 9 genome and epigenome editing showed that the complex CAD genetic signals within a genomic region at 2q22.3 lie within smooth muscle long-distance enhancers for ZEB2, a transcription factor extensively studied in the context of epithelial mesenchymal transition in development of cancer. Zeb2 regulates SMC phenotypic transition through chromatin remodeling that obviates accessibility and disrupts both Notch and transforming growth factor β signaling, thus altering the epigenetic trajectory of SMC transitions. SMC-specific loss of Zeb2 resulted in an inability of transitioning SMCs to turn off contractile programing and take on a fibroblast-like phenotype, but accelerated the formation of chondromyocytes, mirroring features of high-risk atherosclerotic plaques in human coronary arteries. These studies identify ZEB2 as a new CAD genome-wide association study gene that affects features of plaque vulnerability through direct effects on the epigenome, providing a new therapeutic approach to target vascular disease.

Abstract 9564: Zeb2 Shapes the Epigenetic Landscape of Atherosclerosis and Modulates the Risk of Myocardial Infarction

Smooth muscle cells (SMC) transition into a number of different phenotypes during atherosclerosis, including those that resemble fibroblasts and chondrocytes, and make up the majority of cells in atherosclerotic plaque. To better understand the epigenetic and transcriptional mechanisms that mediate these cell state changes, and how they relate to risk for coronary artery disease (CAD), we have investigated the causality and function of transcription factors (TFs) at genome wide associated loci. Employing in vitro CRISPR-Cas 9 genome and epigenome editing we show that complex genetic signals within a gene desert at 2q22.3 lie within novel smooth muscle long-distance enhancers for ZEB2 , a TF extensively studied in the context of epithelial mesenchymal transition (EMT) in development and cancer. Single-cell epigenetic and transcriptomic profiling demonstrated that Zeb2 regulates SMC phenotypic transition through chromatin remodeling that obviates accessibility and disrupts both Notch and TGFβ signaling, altering the epigenetic trajectory of SMC transitions. These changes resulted in an inability for transition SMC to turn off contractile programing and take on a fibroblast-like phenotype, and accelerated chondromyocyte differentiation, mirroring features of high-risk atherosclerotic plaques in human coronary arteries. Parallel epigenetic alterations are observed in human coronary artery smooth muscle cells, with ZEB binding motifs preferentially enriched for human genetic signals for coronary artery disease, affirming its role as a master epigenetic regulator of atherosclerosis pathology.

Type 2 diabetes is associated with the MTNR1B gene, a genetic bridge between circadian rhythm and glucose metabolism, in a Turkish population.

Type 2 diabetes (T2D) is a complicated public health problem in Turkey as well as worldwide. Genome-wide approaches have been guiding in very challenging situations, such as the elucidation of genetic variations underlying complex diseases such as T2D. Despite intensive studies worldwide, few studies have determined the genetic susceptibility to T2D in Turkish populations. In this study, we investigated the effect of genes that are strongly associated with T2D in genome-wide association (GWA) studies, including MTNR1B, CDKAL1, THADA, ADAMTS9 and ENPP1, on T2D and its characteristic traits in a Turkish population. In 824 nonobese individuals (454 T2D patients and 370 healthy individuals), prominent variants of these GWA genes were genotyped by real-time PCR using the LightSNiP Genotyping Assay System. The SNP rs1387153 C/T, which is located 28kb upstream of the MTNR1B gene, was significantly associated with T2D and fasting blood glucose levels (P < 0.05). The intronic SNP rs10830963 C/G in the MTNR1B gene was not associated with T2D, but it was associated with fasting blood glucose, HbA1C and LDL levels (P < 0.05). The other important GWA loci investigated in our study were not found to be associated with T2D or its traits. Only the SNP rs1044498 (A/C variation) in the ENPP1 gene was determined to be related to fasting blood glucose (P < 0.05). Our study suggests, consistent with the literature, that the MTNR1B locus, which has a prominent role in glucose regulation, is associated with T2D development by affecting blood glucose levels in our population.

The impact of non-additive genetic associations on age-related complex diseases

Genome-wide association studies (GWAS) are not fully comprehensive, as current strategies typically test only the additive model, exclude the X chromosome, and use only one reference panel for genotype imputation. We implement an extensive GWAS strategy, GUIDANCE, which improves genotype imputation by using multiple reference panels and includes the analysis of the X chromosome and non-additive models to test for association. We apply this methodology to 62,281 subjects across 22 age-related diseases and identify 94 genome-wide associated loci, including 26 previously unreported. Moreover, we observe that 27.7% of the 94 loci are missed if we use standard imputation strategies with a single reference panel, such as HRC, and only test the additive model. Among the new findings, we identify three novel low-frequency recessive variants with odds ratios larger than 4, which need at least a three-fold larger sample size to be detected under the additive model. This study highlights the benefits of applying innovative strategies to better uncover the genetic architecture of complex diseases.

Single nucleotide variations in ZBTB46 are associated with post-thrombolytic parenchymal haematoma.

Haemorrhagic transformation is a complication of recombinant tissue-plasminogen activator treatment. The most severe form, parenchymal haematoma, can result in neurological deterioration, disability, and death. Our objective was to identify single nucleotide variations associated with a risk of parenchymal haematoma following thrombolytic therapy in patients with acute ischaemic stroke. A fixed-effect genome-wide meta-analysis was performed combining two-stage genome-wide association studies (n = 1904). The discovery stage (three cohorts) comprised 1324 ischaemic stroke individuals, 5.4% of whom had a parenchymal haematoma. Genetic variants yielding a P-value < 0.05 1 × 10-5 were analysed in the validation stage (six cohorts), formed by 580 ischaemic stroke patients with 12.1% haemorrhagic events. All participants received recombinant tissue-plasminogen activator; cases were parenchymal haematoma type 1 or 2 as defined by the European Cooperative Acute Stroke Study (ECASS) criteria. Genome-wide significant findings (P < 5 × 10-8) were characterized by in silico functional annotation, gene expression, and DNA regulatory elements. We analysed 7 989 272 single nucleotide polymorphisms and identified a genome-wide association locus on chromosome 20 in the discovery cohort; functional annotation indicated that the ZBTB46 gene was driving the association for chromosome 20. The top single nucleotide polymorphism was rs76484331 in the ZBTB46 gene [P = 2.49 × 10-8; odds ratio (OR): 11.21; 95% confidence interval (CI): 4.82-26.55]. In the replication cohort (n = 580), the rs76484331 polymorphism was associated with parenchymal haematoma (P = 0.01), and the overall association after meta-analysis increased (P = 1.61 × 10-8; OR: 5.84; 95% CI: 3.16-10.76). ZBTB46 codes the zinc finger and BTB domain-containing protein 46 that acts as a transcription factor. In silico studies indicated that ZBTB46 is expressed in brain tissue by neurons and endothelial cells. Moreover, rs76484331 interacts with the promoter sites located at 20q13. In conclusion, we identified single nucleotide variants in the ZBTB46 gene associated with a higher risk of parenchymal haematoma following recombinant tissue-plasminogen activator treatment.

Integrative analysis of liver-specific non-coding regulatory SNPs associated with the risk of coronary artery disease

Genetic factors underlying coronary artery disease (CAD) have been widely studied using genome-wide association studies (GWASs). However, the functional understanding of the CAD loci has been limited by the fact that a majority of GWAS variants are located within non-coding regions with no functional role. High cholesterol and dysregulation of the liver metabolism such as non-alcoholic fatty liver disease confer an increased risk of CAD. Here, we studied the function of non-coding single-nucleotide polymorphisms in CAD GWAS loci located within liver-specific enhancer elements by identifying their potential target genes using liver cis-eQTL analysis and promoter Capture Hi-C in HepG2 cells. Altogether, 734 target genes were identified of which 121 exhibited correlations to liver-related traits. To identify potentially causal regulatory SNPs, the allele-specific enhancer activity was analyzed by (1) sequence-based computational predictions, (2) quantification of allele-specific transcription factor binding, and (3) STARR-seq massively parallel reporter assay. Altogether, our analysis identified 1,277 unique SNPs that display allele-specific regulatory activity. Among these, susceptibility enhancers near important cholesterol homeostasis genes (APOB, APOC1, APOE, and LIPA) were identified, suggesting that altered gene regulatory activity could represent another way by which genetic variation regulates serum lipoprotein levels. Using CRISPR-based perturbation, we demonstrate how the deletion/activation of a single enhancer leads to changes in the expression of many target genes located in a shared chromatin interaction domain. Our integrative genomics approach represents a comprehensive effort in identifying putative causal regulatory regions and target genes that could predispose to clinical manifestation of CAD by affecting liver function.

The overlap of genetic susceptibility to schizophrenia and cardiometabolic disease can be used to identify metabolically different groups of individuals

Understanding why individuals with severe mental illness (Schizophrenia, Bipolar Disorder and Major Depressive Disorder) have increased risk of cardiometabolic disease (including obesity, type 2 diabetes and cardiovascular disease), and identifying those at highest risk of cardiometabolic disease are important priority areas for researchers. For individuals with European ancestry we explored whether genetic variation could identify sub-groups with different metabolic profiles. Loci associated with schizophrenia, bipolar disorder and major depressive disorder from previous genome-wide association studies and loci that were also implicated in cardiometabolic processes and diseases were selected. In the IMPROVE study (a high cardiovascular risk sample) and UK Biobank (general population sample) multidimensional scaling was applied to genetic variants implicated in both psychiatric and cardiometabolic disorders. Visual inspection of the resulting plots used to identify distinct clusters. Differences between these clusters were assessed using chi-squared and Kruskall-Wallis tests. In IMPROVE, genetic loci associated with both schizophrenia and cardiometabolic disease (but not bipolar disorder or major depressive disorder) identified three groups of individuals with distinct metabolic profiles. This grouping was replicated within UK Biobank, with somewhat less distinction between metabolic profiles. This work focused on individuals of European ancestry and is unlikely to apply to more genetically diverse populations. Overall, this study provides proof of concept that common biology underlying mental and physical illness may help to stratify subsets of individuals with different cardiometabolic profiles.

Genetic variation in GC and CYP2R1 affects 25-hydroxyvitamin D concentration and skeletal parameters: A genome-wide association study in 24-month-old Finnish children

Vitamin D is important for normal skeletal homeostasis, especially in growing children. There are no previous genome-wide association (GWA) studies exploring genetic factors that influence vitamin D metabolism in early childhood. We performed a GWA study on serum 25-hydroxyvitamin D (25(OH)D) and response to supplementation in 761 healthy term-born Finnish 24-month-old children, who participated in a randomized clinical trial comparing effects of 10 μg and 30 μg of daily vitamin D supplementation from age 2 weeks to 24 months. Using the Illumina Infinium Global Screening Array, which has been optimized for imputation, a total of 686085 markers were genotyped across the genome. Serum 25(OH)D was measured at the end of the intervention at 24 months of age. Skeletal parameters reflecting bone strength were determined at the distal tibia at 24 months using peripheral quantitative computed tomography (pQCT) (data available for 648 children). For 25(OH)D, two strong GWA signals were identified, localizing to GC (Vitamin D binding protein) and CYP2R1 (Vitamin D 25-hydroxylase) genes. The GWA locus comprising the GC gene also associated with response to supplementation. Further evidence for the importance of these two genes was obtained by comparing association signals to gene expression data from the Genotype-Tissue Expression project and performing colocalization analyses. Through the identification of haplotypes associated with low or high 25(OH)D concentrations we used a Mendelian randomization approach to show that haplotypes associating with low 25(OH)D were also associated with low pQCT parameters in the 24-month-old children. In this first GWA study on 25(OH)D in this age group we show that already at the age of 24 months genetic variation influences 25(OH)D concentrations and determines response to supplementation, with genome-wide significant associations with GC and CYP2R1. Also, the dual association between haplotypes, 25(OH)D and pQCT parameters gives support for vertical pleiotropy mediated by 25(OH)D.

Investigating APOE, APP-A\u03b2 metabolism genes and Alzheimer\u2019s disease GWAS hits in brain small vessel ischemic disease

Alzheimer’s disease and small vessel ischemic disease frequently co-exist in the aging brain. However, pathogenic links between these 2 disorders are yet to be identified. Therefore we used Taqman genotyping, exome and RNA sequencing to investigate Alzheimer’s disease known pathogenic variants and pathways: APOE ε4 allele, APP-Aβ metabolism and late-onset Alzheimer’s disease main genome-wide association loci (APOE, BIN1, CD33, MS4A6A, CD2AP, PICALM, CLU, CR1, EPHA1, ABCA7) in 96 early-onset small vessel ischemic disease Caucasian patients and 368 elderly neuropathologically proven controls (HEX database) and in a mouse model of cerebral hypoperfusion. Only a minority of patients (29%) carried APOE ε4 allele. We did not detect any pathogenic mutation in APP, PSEN1 and PSEN2 and report a burden of truncating mutations in APP-Aß degradation genes. The single-variant association test identified 3 common variants with a likely protective effect on small vessel ischemic disease (0.54>OR > 0.32, adj. p-value <0.05) (EPHA1 p.M900V and p.V160A and CD33 p.A14V). Moreover, 5/17 APP-Aß catabolism genes were significantly upregulated (LogFC > 1, adj. p-val<0.05) together with Apoe, Ms4a cluster and Cd33 during brain hypoperfusion and their overexpression correlated with the ischemic lesion size. Finally, the detection of Aβ oligomers in the hypoperfused hippocampus supported the link between brain ischemia and Alzheimer’s disease pathology.

Genetic risk of extranodal natural killer T-cell lymphoma: a genome-wide association study in multiple populations.

Extranodal natural killer T-cell lymphoma (NKTCL; nasal type) is an aggressive malignancy with a particularly high prevalence in Asian and Latin American populations. Epstein-Barr virus infection has a role in the pathogenesis of NKTCL, and HLA-DPB1 variants are risk factors for the disease. We aimed to identify additional novel genetic variants affecting risk of NKTCL. We did a genome-wide association study of NKTCL in multiple populations from east Asia. We recruited a discovery cohort of 700 cases with NKTCL and 7752 controls without NKTCL of Han Chinese ancestry from 19 centres in southern, central, and northern regions of China, and four independent replication samples including 717 cases and 12 650 controls. Three of these independent samples (451 cases and 5301 controls) were from eight centres in the same regions of southern, central, and northern China, and the fourth (266 cases and 7349 controls) was from 11 centres in Hong Kong, Taiwan, Singapore, and South Korea. All cases had primary NKTCL that was confirmed histopathologically, and matching with controls was based on geographical region and self-reported ancestry. Logistic regression analysis was done independently by geographical regions, followed by fixed-effect meta-analyses, to identify susceptibility loci. Bioinformatic approaches, including expression quantitative trait loci, binding motif and transcriptome analyses, and biological experiments were done to fine-map and explore the functional relevance of genome-wide association loci to the development of NKTCL. Genetic data were gathered between Jan 1, 2008, and Jan 23, 2019. Meta-analysis of all samples (a total of 1417 cases and 20 402 controls) identified two novel loci significantly associated with NKTCL: IL18RAP on 2q12.1 (rs13015714; p=2·83 × 10-16; odds ratio 1·39 [95% CI 1·28-1·50]) and HLA-DRB1 on 6p21.3 (rs9271588; 9·35 × 10-26 1·53 [1·41-1·65]). Fine-mapping and experimental analyses showed that rs1420106 at the promoter of IL18RAP was highly correlated with rs13015714, and the rs1420106-A risk variant had an upregulatory effect on IL18RAP expression. Cell growth assays in two NKTCL cell lines (YT and SNK-6 cells) showed that knockdown of IL18RAP inhibited cell proliferation by cell cycle arrest in NKTCL cells. Haplotype association analysis showed that haplotype 47F-67I was associated with reduced risk of NKTCL, whereas 47Y-67L was associated with increased risk of NKTCL. These two positions are component parts of the peptide-binding pocket 7 (P7) of the HLA-DR heterodimer, suggesting that these alterations might account for the association at HLA-DRB1, independent of the previously reported HLA-DPB1 variants. Our findings provide new insights into the development of NKTCL by showing the importance of inflammation and immune regulation through the IL18-IL18RAP axis and antigen presentation involving HLA-DRB1, which might help to identify potential therapeutic targets. Taken in combination with additional genetic and other risk factors, our results could potentially be used to stratify people at high risk of NKTCL for targeted prevention. Guangdong Innovative and Entrepreneurial Research Team Program, National Natural Science Foundation of China, National Program for Support of Top-Notch Young Professionals, Chang Jiang Scholars Program, Singapore Ministry of Health's National Medical Research Council, Tanoto Foundation, National Research Foundation Singapore, Chang Gung Memorial Hospital, Recruitment Program for Young Professionals of China, First Affiliated Hospital and Army Medical University, US National Institutes of Health, and US National Cancer Institute.

Genetic variation in GC and CYP2R1 affects 25-hydroxyvitamin D concentration and skeletal parameters: A genome-wide association study in 24-month-old Finnish children.

Vitamin D is important for normal skeletal homeostasis, especially in growing children. There are no previous genome-wide association (GWA) studies exploring genetic factors that influence vitamin D metabolism in early childhood. We performed a GWA study on serum 25-hydroxyvitamin D (25(OH)D) and response to supplementation in 761 healthy term-born Finnish 24-month-old children, who participated in a randomized clinical trial comparing effects of 10 μg and 30 μg of daily vitamin D supplementation from age 2 weeks to 24 months. Using the Illumina Infinium Global Screening Array, which has been optimized for imputation, a total of 686085 markers were genotyped across the genome. Serum 25(OH)D was measured at the end of the intervention at 24 months of age. Skeletal parameters reflecting bone strength were determined at the distal tibia at 24 months using peripheral quantitative computed tomography (pQCT) (data available for 648 children). For 25(OH)D, two strong GWA signals were identified, localizing to GC (Vitamin D binding protein) and CYP2R1 (Vitamin D 25-hydroxylase) genes. The GWA locus comprising the GC gene also associated with response to supplementation. Further evidence for the importance of these two genes was obtained by comparing association signals to gene expression data from the Genotype-Tissue Expression project and performing colocalization analyses. Through the identification of haplotypes associated with low or high 25(OH)D concentrations we used a Mendelian randomization approach to show that haplotypes associating with low 25(OH)D were also associated with low pQCT parameters in the 24-month-old children. In this first GWA study on 25(OH)D in this age group we show that already at the age of 24 months genetic variation influences 25(OH)D concentrations and determines response to supplementation, with genome-wide significant associations with GC and CYP2R1. Also, the dual association between haplotypes, 25(OH)D and pQCT parameters gives support for vertical pleiotropy mediated by 25(OH)D.