BACKGROUND: Genetic variants in complement genes have been associated with a wide range of human disease states, but well-powered genetic association studies of complement activation have not been performed in large multiethnic cohorts. METHODS: We performed medical records-based genome-wide and phenome-wide association studies for plasma C3 and C4 levels among participants of the Electronic Medical Records and Genomics (eMERGE) network. RESULTS: In a GWAS for C3 levels in 3949 individuals, we detected two genome-wide significant loci: chr.1q31.3 (CFH locus; rs3753396-A; β=0.20; 95% CI, 0.14 to 0.25; P=1.52x10-11) and chr.19p13.3 (C3 locus; rs11569470-G; β=0.19; 95% CI, 0.13 to 0.24; P=1.29x10-8). These two loci explained approximately 2% of variance in C3 levels. GWAS for C4 levels involved 3998 individuals and revealed a genome-wide significant locus at chr.6p21.32 (C4 locus; rs3135353-C; β=0.40; 95% CI, 0.34 to 0.45; P=4.58x10-35). This locus explained approximately 13% of variance in C4 levels. The multiallelic copy number variant analysis defined two structural genomic C4 variants with large effect on blood C4 levels: C4-BS (β=-0.36; 95% CI, -0.42 to -0.30; P=2.98x10-22) and C4-AL-BS (β=0.25; 95% CI, 0.21 to 0.29; P=8.11x10-23). Overall, C4 levels were strongly correlated with copy numbers of C4A and C4B genes. In comprehensive phenome-wide association studies involving 102,138 eMERGE participants, we cataloged a full spectrum of autoimmune, cardiometabolic, and kidney diseases genetically related to systemic complement activation. CONCLUSIONS: We discovered genetic determinants of plasma C3 and C4 levels using eMERGE genomic data linked to electronic medical records. Genetic variants regulating C3 and C4 levels have large effects and multiple clinical correlations across the spectrum of complement-related diseases in humans.
Genetic factors play a key role in the pathogenesis of autoimmune diseases, whereas the disease-causing variants remain largely unknown. Herein, we performed an exome-wide association study of systemic sclerosis in a Han Chinese population. In the discovery stage, 527 patients with systemic sclerosis and 5,024 controls were recruited and genotyped. In the validation study, an independent sample set of 479 patients and 1,096 controls were examined. In total, we found that four independent signals reached genome-wide significance. Among them, rs7574865 (Pcombined = 3.87 × 10-12) located within signal transducer and activator of transcription 4 gene was identified previously using samples of European ancestry. Additionally, another signal including three SNPs in linkage disequilibrium might be unreported susceptibility loci located in the epidermis differentiation complex region. Furthermore, two SNPs located within exon 3 of IGHM (rs45471499, Pcombined = 1.15 × 10-9) and upstream of LRP2BP (rs4317244, Pcombined = 4.17 × 10-8) were found. Moreover, rs4317244 was identified as an expression quantitative trait locus for LRP2BP that regulates tight junctions, cell cycle, and apoptosis in endothelial cell lines. Collectively, our results revealed three signals associated with systemic sclerosis in Han Chinese and suggested the importance of LRP2BP in systemic sclerosis pathogenesis. Given the limited sample size and discrepancies between previous results and our study, further studies in multiethnic populations are required for verification.
BACKGROUND: Posterior urethral valves (PUVs) and ureteropelvic junction obstruction (UPJO) are congenital obstructive uropathies that may impair kidney development. OBJECTIVE: To identify genetic variants associated with kidney injury in patients with obstructive uropathy. DESIGN SETTING AND PARTICIPANTS: We included 487 patients born in 1981 or later who underwent pyeloplasty or valve resection before 18 yr of age in the discovery phase, 102 PUV patients in a first replication phase, and 102 in a second replication phase. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: Signs of kidney injury were defined as dialysis, nephrectomy, kidney transplantation, estimated glomerular filtration rate (eGFR) <60 ml/min/1.73 m2, high blood pressure, antihypertensive medication use, proteinuria, and/or one kidney functioning at <45%. We used χ2 tests to calculate p values and odds ratios for >600 000 single-nucleotide polymorphisms (SNPs) in the discovery sample comparing patients with and without signs of kidney injury within 5 yr after surgery. We performed stratified analyses for PUV and UPJO and Kaplan-Meier and Cox regression analyses in the discovery and two replication samples for the associated SNPs, and RNA and protein expression analyses for the associated gene in fetal tissues. RESULTS AND LIMITATIONS: Despite the small and nonhomogeneous sample, we observed suggestive associations for six SNPs in three loci, of which rs6874819 in the CDH12 gene was the most clear (p = 7.5 × 10-7). This SNP also seemed to be associated with time to kidney injury in the PUV discovery and replication samples. RNA expression analyses showed clear CDH12 expression in fetal kidneys, which was confirmed by protein immunolocalization. CONCLUSIONS: This study identified CDH12 as a candidate gene for kidney injury in PUV. PATIENT SUMMARY: We found that variants of the CDH12 gene increase the risk of kidney injury in patients with extra flaps of tissue in the urethra (posterior urethral valves). This is the first report on this gene in this context. Our study provides interesting new information about the pathways involved and important leads for further research for this condition.
Kidney injury 5 years post-surgery in obstructive uropathies
Shigella is a leading cause of moderate-to-severe diarrhea globally and the causative agent of shigellosis and bacillary dysentery. Associated with 80 to 165 million cases of diarrhea and >13% of diarrheal deaths, in many regions, Shigella exposure is ubiquitous while infection is heterogenous. To characterize host-genetic susceptibility to Shigella-associated diarrhea, we performed two independent genome-wide association studies (GWAS) including Bangladeshi infants from the PROVIDE and CBC birth cohorts in Dhaka, Bangladesh. Cases were infants with Shigella-associated diarrhea (n = 143) and controls were infants with no Shigella-associated diarrhea in the first 13 months of life (n = 446). Shigella-associated diarrhea was identified via quantitative PCR (qPCR) threshold cycle (CT ) distributions for the ipaH gene, carried by all four Shigella species and enteroinvasive Escherichia coli Host GWAS were performed under an additive genetic model. A joint analysis identified protective loci on chromosomes 11 (rs582240, within the KRT18P59 pseudogene; P = 6.40 × 10-8; odds ratio [OR], 0.43) and 8 (rs12550437, within the lincRNA RP11-115J16.1; P = 1.49 × 10-7; OR, 0.48). Conditional analyses identified two previously suggestive loci, a protective locus on chromosome 7 (rs10266841, within the 3' untranslated region [UTR] of CYTH3; P conditional = 1.48 × 10-7; OR, 0.44) and a risk-associated locus on chromosome 10 (rs2801847, an intronic variant within MPP7; P conditional = 8.37 × 10-8; OR, 5.51). These loci have all been indirectly linked to bacterial type 3 secretion system (T3SS) activity, its components, and bacterial effectors delivered into host cells. Host genetic factors that may affect bacterial T3SS activity and are associated with the host response to Shigella-associated diarrhea may provide insight into vaccine and drug development efforts for Shigella-associated diarrheal disease.
BACKGROUND: Abdominal aortic aneurysm (AAA) is an important cause of cardiovascular mortality; however, its genetic determinants remain incompletely defined. In total, 10 previously identified risk loci explain a small fraction of AAA heritability. METHODS: We performed a genome-wide association study in the Million Veteran Program testing ≈18 million DNA sequence variants with AAA (7642 cases and 172 172 controls) in veterans of European ancestry with independent replication in up to 4972 cases and 99 858 controls. We then used mendelian randomization to examine the causal effects of blood pressure on AAA. We examined the association of AAA risk variants with aneurysms in the lower extremity, cerebral, and iliac arterial beds, and derived a genome-wide polygenic risk score (PRS) to identify a subset of the population at greater risk for disease. RESULTS: Through a genome-wide association study, we identified 14 novel loci, bringing the total number of known significant AAA loci to 24. In our mendelian randomization analysis, we demonstrate that a genetic increase of 10 mm Hg in diastolic blood pressure (odds ratio, 1.43 [95% CI, 1.24-1.66]; P=1.6×10-6), as opposed to systolic blood pressure (odds ratio, 1.06 [95% CI, 0.97-1.15]; P=0.2), likely has a causal relationship with AAA development. We observed that 19 of 24 AAA risk variants associate with aneurysms in at least 1 other vascular territory. A 29-variant PRS was strongly associated with AAA (odds ratioPRS, 1.26 [95% CI, 1.18-1.36]; PPRS=2.7×10-11 per SD increase in PRS), independent of family history and smoking risk factors (odds ratioPRS+family history+smoking, 1.24 [95% CI, 1.14-1.35]; PPRS=1.27×10-6). Using this PRS, we identified a subset of the population with AAA prevalence greater than that observed in screening trials informing current guidelines. CONCLUSIONS: We identify novel AAA genetic associations with therapeutic implications and identify a subset of the population at significantly increased genetic risk of AAA independent of family history. Our data suggest that extending current screening guidelines to include testing to identify those with high polygenic AAA risk, once the cost of genotyping becomes comparable with that of screening ultrasound, would significantly increase the yield of current screening at reasonable cost.
BACKGROUND: Lung cancer remains the leading cause of cancer death worldwide, with 15% to 20% occurring in never smokers. To assess genetic determinants for prognosis among never smokers, we conducted a genome-wide investigation in the International Lung Cancer Consortium (ILCCO). METHODS: Genomic and clinical data from 1,569 never-smoking patients with lung cancer of European ancestry from 10 ILCCO studies were included. HRs and 95% confidence intervals of overall survival were estimated. We assessed whether the associations were mediated through mRNA expression-based 1,553 normal lung tissues from the lung expression quantitative trait loci (eQTL) dataset and Genotype-Tissue Expression (GTEx). For cross-ethnicity generalization, we assessed the associations in a Japanese study (N = 887). RESULTS: One locus at 13q22.2 was associated with lung adenocarcinoma survival at genome-wide level, with carriers of rs12875562-T allele exhibiting poor prognosis [HR = 1.71 (1.41-2.07), P = 3.60 × 10-8], and altered mRNA expression of LMO7DN in lung tissue (GTEx, P = 9.40 × 10-7; Lung eQTL dataset, P = 0.003). Furthermore, 2 of 11 independent loci that reached the suggestive significance level (P < 10-6) were significant eQTL affecting mRNA expression of nearby genes in lung tissues, including CAPZB at 1p36.13 and UBAC1 at 9q34.3. One locus encoding NWD2/KIAA1239 at 4p14 showed associations in both European [HR = 0.50 (0.38-0.66), P = 6.92 × 10-7] and Japanese populations [HR = 0.79 (0.67-0.94), P = 0.007]. CONCLUSIONS: Based on the largest genomic investigation on the lung cancer prognosis of never smokers to date, we observed that lung cancer prognosis is affected by inherited genetic variants. IMPACT: We identified one locus near LMO7DN at genome-wide level and several potential prognostic genes with cis-effect on mRNA expression. Further functional genomics work is required to understand their role in tumor progression.
Lung adenocarcinoma (survival time) in never smokers
Miscarriage is a common, complex trait affecting ~15% of clinically confirmed pregnancies. Here we present the results of large-scale genetic association analyses with 69,054 cases from five different ancestries for sporadic miscarriage, 750 cases of European ancestry for multiple (≥3) consecutive miscarriage, and up to 359,469 female controls. We identify one genome-wide significant association (rs146350366, minor allele frequency (MAF) 1.2%, P = 3.2 × 10-8, odds ratio (OR) = 1.4) for sporadic miscarriage in our European ancestry meta-analysis and three genome-wide significant associations for multiple consecutive miscarriage (rs7859844, MAF = 6.4%, P = 1.3 × 10-8, OR = 1.7; rs143445068, MAF = 0.8%, P = 5.2 × 10-9, OR = 3.4; rs183453668, MAF = 0.5%, P = 2.8 × 10-8, OR = 3.8). We further investigate the genetic architecture of miscarriage with biobank-scale Mendelian randomization, heritability, and genetic correlation analyses. Our results show that miscarriage etiopathogenesis is partly driven by genetic variation potentially related to placental biology, and illustrate the utility of large-scale biobank data for understanding this pregnancy complication.
Emphysema is an important feature of chronic obstructive pulmonary disease (COPD). Genetic factors likely affect emphysema pathogenesis, but this question has predominantly been studied in those of European ancestry. In this study, we sought to determine genetic components of emphysema severity and characterize the potential function of the associated loci in Korean population. We performed a genome-wide association study (GWAS) on quantitative emphysema in subjects with or without COPD from two Korean COPD cohorts. We investigated the functional consequences of the loci using epigenetic annotation and gene expression data. We also compared our GWAS results with an epigenome-wide association study and previous differential gene expression analysis. In total, 548 subjects (476 [86.9%] male) including 514 COPD patients were evaluated. We identified one genome-wide significant SNP (P < 5.0 × 10-8), rs117084279, near PIBF1. We identified an additional 57 SNPs (P < 5.0 × 10-6) associated with emphysema in all subjects, and 106 SNPs (P < 5.0 × 10-6) in COPD patients. Of these candidate SNPs, 2 (rs12459249, rs11667314) near CYP2A6 were expression quantitative trait loci in lung tissue and a SNP (rs11214944) near NNMT was an expression quantitative trait locus in whole blood. Of note, rs11214944 was in linkage disequilibrium with variants in enhancer histone marks in lung tissue. Several genes near additional SNPs were identified in our previous EWAS study with nominal level of significance. We identified a novel SNP associated with quantitative emphysema on CT. Including the novel SNP, several candidate SNPs in our study may provide clues to the genetic etiology of emphysema in Asian populations. Further research and validation of the loci will help determine the genetic factors for the development of emphysema.
Emphysema
Emphysema in chronic obstructive pulmonary disease
BACKGROUND: Numerous genes have been associated with the three most common allergic diseases (asthma, allergic rhinitis or eczema) but these genes explain only a part of the heritability. In the vast majority of genetic studies, complex phenotypes such as co-morbidity of two of these diseases, have not been considered. This may partly explain missing heritability. OBJECTIVE: To identify genetic variants specifically associated with the co-morbidity of asthma-plus-eczema. METHODS: We first conducted a meta-analysis of four GWAS (Genome-Wide Association Study) of the combined asthma-plus-eczema phenotype (total of 8807 European-ancestry subjects of whom 1208 subjects had both asthma and eczema). To assess whether the association with SNP(s) was specific to the co-morbidity, we also conducted a meta-analysis of homogeneity test of association according to disease status ("asthma-plus-eczema" vs. the presence of only one disease "asthma only or eczema only"). We then used a joint test by combining the two test statistics from the co-morbidity-SNP association and the phenotypic heterogeneity of SNP effect meta-analyses. RESULTS: Seven SNPs were detected for specific association to the asthma-plus-eczema co-morbidity, two with significant and five with suggestive evidence using the joint test after correction for multiple testing. The two significant SNPs are located in the OCA2 gene (Oculocutaneous Albinism II), a new locus never detected for significant evidence of association with any allergic disease. This gene is a promising candidate gene, because of its link to skin and lung diseases, and to epithelial barrier and immune mechanisms. CONCLUSION: Our study underlines the importance of studying sub-phenotypes as co-morbidities to detect new susceptibility genes.
BACKGROUND: Buprenorphine, approved for treating opioid use disorder (OUD), is not equally efficacious for all patients. Candidate gene studies have shown limited success in identifying genetic moderators of buprenorphine treatment response. METHODS: We studied 1616 European-ancestry individuals enrolled in the Million Veteran Program, of whom 1609 had an ICD-9/10 code consistent with OUD, a 180-day buprenorphine treatment exposure, and genome-wide genotype data. We conducted a genome-wide association study (GWAS) of buprenorphine treatment response [defined as having no opioid-positive urine drug screens (UDS) following the first prescription]. We also examined correlates of buprenorphine treatment response in multivariable analyses. RESULTS: Although no variants reached genome-wide significance, 6 loci were nominally significant (p < 1 × 10-5), four of which were located near previously characterized genes: rs756770 (ADAMTSL2), rs11782370 (SLC25A37), rs7205113 (CRISPLD2), and rs13169373 (LINC01947). A higher maximum daily buprenorphine dosage (aOR = 0.98; 95 %CI: 0.97, 0.995), greater number of UDS (aOR = 0.97; 95 %CI: 0.96, 0.99), and history of hepatitis C (HCV) infection (aOR = 0.71; 95 %CI: 0.57, 0.88) were associated with a reduced odds of buprenorphine response. Older age (aOR: 1.01; 95 %CI: 1.000, 1.02) was associated with increased odds of buprenorphine response. CONCLUSIONS: This study had limited statistical power to detect genetic variants associated with a complex human phenotype like buprenorphine treatment response. Meta-analysis of multiple data sets is needed to ensure adequate statistical power for a GWAS of buprenorphine treatment response. The most robust phenotypic predictor of buprenorphine treatment response was intravenous drug use, a proxy for which was HCV infection.
Psychiatric disorders are highly genetically correlated, but little research has been conducted on the genetic differences between disorders. We developed a new method (case-case genome-wide association study; CC-GWAS) to test for differences in allele frequency between cases of two disorders using summary statistics from the respective case-control GWAS, transcending current methods that require individual-level data. Simulations and analytical computations confirm that CC-GWAS is well powered with effective control of type I error. We applied CC-GWAS to publicly available summary statistics for schizophrenia, bipolar disorder, major depressive disorder and five other psychiatric disorders. CC-GWAS identified 196 independent case-case loci, including 72 CC-GWAS-specific loci that were not significant at the genome-wide level in the input case-control summary statistics; two of the CC-GWAS-specific loci implicate the genes KLF6 and KLF16 (from the Krüppel-like family of transcription factors), which have been linked to neurite outgrowth and axon regeneration. CC-GWAS loci replicated convincingly in applications to datasets with independent replication data.
Bipolar disorder vs anorexia nervosa (ordinary least squares (OLS))
Anorexia nervosa vs Tourette's syndrome and other tic disorders (ordinary least squares (OLS))
Schizophrenia vs anorexia nervosa (ordinary least squares (OLS))
Autism spectrum disorder vs Tourette's syndrome and other tic disorders (ordinary least squares (OLS))
Autism spectrum disorder vs obsessive compulsive disorder (ordinary least squares (OLS))
Ulcerative colitis vs rheumatoid arthritis (ordinary least squares (OLS))
Bipolar disorder vs ADHD (ordinary least squares (OLS))
Major depressive disorder vs ADHD (ordinary least squares (OLS))
Bipolar disorder vs major depressive disorder (ordinary least squares (OLS))
Major depressive disorder vs anorexia nervosa (ordinary least squares (OLS))
Schizophrenia vs autism spectrum disorder (ordinary least squares (OLS))
Bipolar disorder vs autism spectrum disorder (ordinary least squares (OLS))
Schizophrenia vs major depressive disorder (ordinary least squares (OLS))
ADHD vs anorexia nervosa (ordinary least squares (OLS))
Schizophrenia vs bipolar disorder (ordinary least squares (OLS))
Bipolar disorder vs obsessive compulsive disorder (ordinary least squares (OLS))
Obsessive compulsive disorder vs Tourette's syndrome and other tic disorders (ordinary least squares (OLS))
Schizophrenia vs Tourette's syndrome and other tic disorders (ordinary least squares (OLS))
Schizophrenia vs ADHD (ordinary least squares (OLS))
ADHD vs autism spectrum disorder (ordinary least squares (OLS))
Bipolar disorder vs Tourette's syndrome and other tic disorders (ordinary least squares (OLS))
Anorexia nervosa vs autism spectrum disorder (ordinary least squares (OLS))
Crohn's disease vs rheumatoid arthritis (ordinary least squares (OLS))
Major depressive disorder vs autism spectrum disorder (ordinary least squares (OLS))
Crohn's disease vs ulcerative colitis (ordinary least squares (OLS))
ADHD vs Tourette's syndrome and other tic disorders (ordinary least squares (OLS))
ADHD vs obsessive compulsive disorder (ordinary least squares (OLS))
In genome-wide association studies, signals associated with rare variants and interactions between genes are hard to detect even when the sample size is in tens of thousands. To overcome these problems, we examine the concept of supervariant. Like the classic concept of the gene, a supervariant is a combination of alleles in multiple loci, but the contributing loci can be anywhere in the genome. We hypothesize that supervariants are easy to detect and the aggregated signals are more stable in their associations with the disease than that from a single nucleoid polymorphism. Using the UK Biobank databases, we develop a ranking and aggregation method for identifying supervariants. Specifically, we examine 9,377 breast cancer cases with 46,861 controls matched by sex and age. In our simulations, the use of supervariants outperforms single-nucleotide polymorphism-based association method in detecting rare variants and signals with interactive structure. In real data analysis, we identify supervariants on Chromosomes 1, 2, 3, 5, 6, 7, 8, 9, 10, 11, 16, and 22 which cover previously reported loci that have associations with breast or other cancers, and several novel loci on Chromosomes 2, 5, 9, and 12. These findings demonstrate the validity of supervariants and its potential of discovering replicable and novel results for complex disease.
BACKGROUND: Type 2 diabetes complications cause a serious emotional and economical burden to patients and healthcare systems globally. Management of both acute and chronic complications of diabetes, which dramatically impair the quality of patients' life, is still an unsolved issue in diabetes care, suggesting a need for early identification of individuals with high risk for developing diabetes complications. METHODS: We performed a genome-wide association study in 601 type 2 diabetes patients after stratifying them according to the presence or absence of four types of diabetes complications: diabetic neuropathy, diabetic nephropathy, macrovascular complications, and ophthalmic complications. RESULTS: The analysis revealed ten novel associations showing genome-wide significance, including rs1132787 (GYPA, OR = 2.71; 95% CI = 2.02-3.64) and diabetic neuropathy, rs2477088 (PDE4DIP, OR = 2.50; 95% CI = 1.87-3.34), rs4852954 (NAT8, OR = 2.27; 95% CI = 2.71-3.01), rs6032 (F5, OR = 2.12; 95% CI = 1.63-2.77), rs6935464 (RPS6KA2, OR = 2.25; 95% CI = 6.69-3.01) and macrovascular complications, rs3095447 (CCDC146, OR = 2.18; 95% CI = 1.66-2.87) and ophthalmic complications. By applying the targeted approach of previously reported susceptibility loci we managed to replicate three associations: MAPK14 (rs3761980, rs80028505) and diabetic neuropathy, APOL1 (rs136161) and diabetic nephropathy. CONCLUSIONS: Together these results provide further evidence for the implication of genetic factors in the development of type 2 diabetes complications and highlight several potential key loci, able to modify the risk of developing these conditions. Moreover, the candidate variant approach proves a strong and consistent effect for multiple variants across different populations.
OBJECTIVES: Nearly 110 susceptibility loci for rheumatoid arthritis (RA) with modest effect sizes have been identified by population-based genetic association studies, suggesting a large number of undiscovered variants behind a highly polygenic genetic architecture of RA. Here, we performed the largest-ever trans-ancestral meta-analysis with the aim to identify new RA loci and to better understand RA biology underlying genetic associations. METHODS: Genome-wide RA association summary statistics in three large case-control collections consisting of 311 292 individuals of Korean, Japanese and European populations were used in an inverse-variance-weighted fixed-effects meta-analysis. Several computational analyses using public omics resources were conducted to prioritise causal variants and genes, RA variant-implicating features (tissues, pathways and transcription factors) and potentially repurposable drugs for RA treatment. RESULTS: We identified 11 new RA susceptibility loci that explained 6.9% and 1.8% of the single-nucleotide polymorphism-based heritability in East Asians and Europeans, respectively, and confirmed 71 known non-human leukocyte antigens (HLA) susceptibility loci, identifying 90 independent association signals. The RA variants were preferentially located in binding sites of various transcription factors and in cell type-specific transcription-activation histone marks that simultaneously highlighted the importance of CD4+ T-cell activation and the potential role of non-immune organs in RA pathogenesis. A total of 615 plausible effector genes, based on gene-based associations, expression-associated variants and chromatin interaction, included targets of drugs approved for RA treatments and potentially repurposable drugs approved for other indications. CONCLUSION: Our findings provide useful insights regarding RA genetic aetiology and variant-driven RA pathogenesis.
Two renal-risk variants in the apolipoprotein L1 gene (APOL1) in African American (AA) deceased donors (DD) are associated with shorter renal allograft survival after transplantation. To identify additional genes contributing to allograft survival, a genome-wide association study was performed in 532 AA DDs. Phenotypic data were obtained from the Scientific Registry of Transplant Recipients. Association and single-nucleotide polymorphism (SNP)-by-APOL1 interaction tests were conducted using death-censored renal allograft survival accounting for relevant covariates. Replication and inverse-variance-weighted meta-analysis were performed using data from 250 AA DD in the Genomics of Transplantation study. Accounting for APOL1, multiple SNPs near the Nudix Hydrolase 7 gene (NUDT7) showed strong independent effects (P = 1.6 × 10-8 -2.2 × 10-8 ). Several SNPs in the Translocation protein SEC63 homolog (SEC63; P = 2 × 10-9 -3.7 × 10-8 ) and plasmacytoma variant translocation 1 (PVT1) genes (P = 4.0 × 10-8 -7 × 10-8 ) modified the effect of APOL1 on allograft survival. SEC63 is expressed in human renal tubule cells and glomeruli, and PVT1 is associated with diabetic kidney disease. Overall, associations were detected for 41 SNPs (P = 2 × 10-9 -5 × 10-8 ) contributing independently or interacting with APOL1 to impact renal allograft survival after transplantation from AA DD. Given the small sample size of the discovery and replication sets, independent validations and functional genomic efforts are needed to validate these results.
Graft survival time in renal transplantation (donor effect) x APOL1 genotype interaction
Graft survival time in renal transplantation (donor effect)
Human skin color is largely determined by genetic factors. Recent GWASs have reported several genetic variants associated with skin color, mostly in European and African populations. In this study, we performed GWAS in 17,019 Korean women to identify genetic variants associated with facial skin color, quantitatively measured as CIELAB color index. We identified variants in three, one, and six genomic loci associated with facial skin color index L∗, a∗, and b∗ values, respectively, and replicated the associations (combined analysis P-value < 5.0 × 10-8). The significant loci included variants in known genes (OCA2 rs74653330, BNC2 rs16935073, rs72620727 near KITLG, and SLC6A17 rs6689641) and to our knowledge previously unreported genes (SCARB1 rs10846744, SYN2 rs12629034, and LINC00486 rs6543678). This is GWAS to elucidate genetic variants of facial skin color in a Korean female population. Further functional characterizations of the investigated genes are warranted to elucidate their contribution to skin pigmentation-related traits.
Importance: The genetic basis of bipolar disorder (BD) in Han Chinese individuals is not fully understood. Objective: To explore the genetic basis of BD in the Han Chinese population. Design, Setting, and Participants: A genome-wide association study (GWAS), followed by independent replication, was conducted to identify BD risk loci in Han Chinese individuals. Individuals with BD were diagnosed based on DSM-IV criteria and had no history of schizophrenia, mental retardation, or substance dependence; individuals without any personal or family history of mental illnesses, including BD, were included as control participants. In total, discovery samples from 1822 patients and 4650 control participants passed quality control for the GWAS analysis. Replication analyses of samples from 958 patients and 2050 control participants were conducted. Summary statistics from the European Psychiatric Genomics Consortium 2 (PGC2) BD GWAS (20 352 cases and 31 358 controls) were used for the trans-ancestry genetic correlation analysis, polygenetic risk score analysis, and meta-analysis to compare BD genetic risk between Han Chinese and European individuals. The study was performed in February 2020. Main Outcomes and Measures: Single-nucleotide variations with P < 5.00 × 10-8 were considered to show genome-wide significance of statistical association. Results: The Han Chinese discovery GWAS sample included 1822 cases (mean [SD] age, 35.43 [14.12] years; 838 [46%] male) and 4650 controls (mean [SD] age, 27.48 [5.97] years; 2465 [53%] male), and the replication sample included 958 cases (mean [SD] age, 37.82 [15.54] years; 412 [43%] male) and 2050 controls (mean [SD] age, 27.50 [6.00] years; 1189 [58%] male). A novel BD risk locus in Han Chinese individuals was found near the gene encoding transmembrane protein 108 (TMEM108, rs9863544; P = 2.49 × 10-8; odds ratio [OR], 0.650; 95% CI, 0.559-0.756), which is required for dendritic spine development and glutamatergic transmission in the dentate gyrus. Trans-ancestry genetic correlation estimation (ρge = 0.652, SE = 0.106; P = 7.30 × 10-10) and polygenetic risk score analyses (maximum liability-scaled Nagelkerke pseudo R2 = 1.27%; P = 1.30 × 10-19) showed evidence of shared BD genetic risk between Han Chinese and European populations, and meta-analysis identified 2 new GWAS risk loci near VRK2 (rs41335055; P = 4.98 × 10-9; OR, 0.849; 95% CI, 0.804-0.897) and RHEBL1 (rs7969091; P = 3.12 × 10-8; OR, 0.932; 95% CI, 0.909-0.956). Conclusions and Relevance: This GWAS study identified several loci and genes involved in the heritable risk of BD, providing insights into its genetic architecture and biological basis.
Autism is heterogeneous multifactorial neurodevelopmental and neuropsychiatric disorder with repetitive and limited behaviors as well as communication deficits. Prevalence of autism in males is predominant than females, but their genetic association is unclear. The study was performed to investigate Y-chromosome haplotypes, significant risk variants and susceptibility genes associated with autistic Saudi young males with autism. Exome genotyping microarray analysis was performed in Saudi young boys with autism (cases, n = 47) and without autism and other genetic or neurodevelopmental disorders (control, n = 43) to identify the functional exonic risk variations among 243,345 exonic variations. The most significant single nucleotide polymorphisms (SNPs) of protein coding associated with autism in Saudi young boys were studied for functional enrichment. Y-chromosome haplotyping analysis of 6 SNPs such as rs1865680, rs2032624, rs2032658, rs2032631, rs9786153 and rs13447352 uncovered the most significant protective (ACGACA p = 2.94 × 10-9) among the controls and the high risk Y-haplotype (GAAGTC p = 6.85 × 10-6) among autistic boys. Exome association study revealed 6 susceptible genes, MCC, AUTS2, VSX1, SETBP1, CNTN3, and PCDH11Y that were known for autistic disorder. The significant predisposed genes with functional variants of Y-chromomere are strongly connected with spermatogenic failure (p = 8.02 × 10-8), azoospermia (p = 6.32 × 10-7), partial chromosome Y deletion (p = 7.66 × 10-6), HDMs demethylate histones pathway (p = 3.55 × 10-4) and immune system diseases (p = 4.11 × 10-3). Y-haplotypes and highly significant pathogenic exonic variants in MCC, AUTS2, VSX1, SETBP1, CNTN3 and PCDH11Y genes are more influential genetic factors for developing autism in boys of Arab origin.
Aiming to uncover a shared genetic basis of abdominal obesity and osteoporosis, we performed a bivariate GWAS meta-analysis of femoral neck BMD (FNK-BMD) and trunk fat mass adjusted by trunk lean mass (TFMadj) in 11,496 subjects from 6 samples, followed by in silico replication in the large-scale UK Biobank (UKB) cohort. A series of functional investigations were conducted on the identified variants. Bivariate GWAS meta-analysis identified two novel pleiotropic loci 12q15 (lead SNP rs73134637, p = 3.45 × 10-7) and 10p14 (lead SNP rs2892347, p = 2.63 × 10-7) that were suggestively associated and that were replicated in the analyses of related traits in the UKB sample (osteoporosis p = 0.06 and 0.02, BMI p = 0.03 and 4.61 × 10-3, N up to 499,520). Cis-eQTL analysis demonstrated that allele C at rs73134637 was positively associated with IFNG expression in whole blood (N = 369, p = 0.04), and allele A at rs11254759 (10p14, p = 9.49 × 10-7) was negatively associated with PRKCQ expression in visceral adipose tissue (N = 313, p = 0.04) and in lymphocytes (N = 117, p = 0.03). As a proof-of-principle experiment, the function of rs11254759, which is 235 kb 5'-upstream from PRKCQ gene, was investigated by the dual-luciferase reporter assay, which clearly showed that the haplotype carrying rs11254759 regulated PRKCQ expression by upregulating PRKCQ promoter activity (p = 4.60 × 10-7) in an allelic specific manner. Mouse model analysis showed that heterozygous PRKCQ deficient mice presented decreased fat mass compared to wild-type control mice (p = 3.30 × 10-3). Mendelian randomization analysis demonstrated that both FNK-BMD and TFMadj were causally associated with fracture risk (p = 1.26 × 10-23 and 1.18 × 10-11). Our findings may provide useful insights into the genetic association between osteoporosis and abdominal obesity.
Femoral neck bone mineral density
Trunk fat mass adjusted for trunk lean mass
Femoral neck bone mineral density and trunk fat mass adjusted by trunk lean mass
How well one does at school is predictive of a wide range of important cognitive, socioeconomic, and health outcomes. The last few years have shown marked advancement in our understanding of the genetic contributions to, and correlations with, academic attainment. However, there exists a gap in our understanding of the specificity of genetic associations with performance in academic subjects during adolescence, a critical developmental period. To address this, the Avon Longitudinal Study of Parents and Children was used to conduct genome-wide association studies of standardised national English (N = 5983), maths (N = 6017) and science (N = 6089) tests. High SNP-based heritabilities (h2SNP) for all subjects were found (41-53%). Further, h2SNP for maths and science remained after removing shared variance between subjects or IQ (N = 3197-5895). One genome-wide significant single nucleotide polymorphism (rs952964, p = 4.86 × 10-8) and four gene-level associations with science attainment (MEF2C, BRINP1, S100A1 and S100A13) were identified. Rs952964 remained significant after removing the variance shared between academic subjects. The findings highlight the benefits of using environmentally homogeneous samples for genetic analyses and indicate that finer-grained phenotyping will help build more specific biological models of variance in learning processes and abilities.
GWAS of Complement system
