Genome-Wide Association Analysis of Protein-Coding Variants in IgA Nephropathy.

SIGNIFICANCE STATEMENT: Genome-wide association studies have identified nearly 20 IgA nephropathy susceptibility loci. However, most nonsynonymous coding variants, particularly ones that occur rarely or at a low frequency, have not been well investigated. The authors performed a chip-based association study of IgA nephropathy in 8529 patients with the disorder and 23,224 controls. They identified a rare variant in the gene encoding vascular endothelial growth factor A (VEGFA) that was significantly associated with a two-fold increased risk of IgA nephropathy, which was further confirmed by sequencing analysis. They also identified a novel common variant in PKD1L3 that was significantly associated with lower haptoglobin protein levels. This study, which was well-powered to detect low-frequency variants with moderate to large effect sizes, helps expand our understanding of the genetic basis of IgA nephropathy susceptibility. BACKGROUND: Genome-wide association studies have identified nearly 20 susceptibility loci for IgA nephropathy. However, most nonsynonymous coding variants, particularly those occurring rarely or at a low frequency, have not been well investigated. METHODS: We performed a three-stage exome chip-based association study of coding variants in 8529 patients with IgA nephropathy and 23,224 controls, all of Han Chinese ancestry. Sequencing analysis was conducted to investigate rare coding variants that were not covered by the exome chip. We used molecular dynamic simulation to characterize the effects of mutations of VEGFA on the protein's structure and function. We also explored the relationship between the identified variants and the risk of disease progression. RESULTS: We discovered a novel rare nonsynonymous risk variant in VEGFA (odds ratio, 1.97; 95% confidence interval [95% CI], 1.61 to 2.41; P = 3.61×10 -11 ). Further sequencing of VEGFA revealed twice as many carriers of other rare variants in 2148 cases compared with 2732 controls. We also identified a common nonsynonymous risk variant in PKD1L3 (odds ratio, 1.16; 95% CI, 1.11 to 1.21; P = 1.43×10 -11 ), which was associated with lower haptoglobin protein levels. The rare VEGFA mutation could cause a conformational change and increase the binding affinity of VEGFA to its receptors. Furthermore, this variant was associated with the increased risk of kidney disease progression in IgA nephropathy (hazard ratio, 2.99; 95% CI, 1.09 to 8.21; P = 0.03). CONCLUSIONS: Our study identified two novel risk variants for IgA nephropathy in VEGFA and PKD1L3 and helps expand our understanding of the genetic basis of IgA nephropathy susceptibility.

A five-safes approach to a secure and scalable genomics data repository.

Genomic researchers increasingly utilize commercial cloud service providers (CSPs) to manage data and analytics needs. CSPs allow researchers to grow Information Technology (IT) infrastructure on demand to overcome bottlenecks when combining large datasets. However, without adequate security controls, the risk of unauthorized access may be higher for data stored on the cloud. Additionally, regulators are mandating data access patterns and specific security protocols for the storage and use of genomic data. While CSP provides tools for security and regulatory compliance, building the necessary controls required for cloud solutions is not trivial. Research Assets Provisioning and Tracking Online Repository (RAPTOR) by the Genome Institute of Singapore is a cloud-native genomics data repository and analytics platform that implements a "five-safes" framework to provide security and governance controls to data contributors and users, leveraging CSP for sharing and analysis of genomic datasets without the risk of security breaches or running afoul of regulations.

Low frequency variants associated with leukocyte telomere length in the Singapore Chinese population.

The role of low frequency variants associated with telomere length homeostasis in chronic diseases and mortalities is relatively understudied in the East-Asian population. Here we evaluated low frequency variants, including 1,915,154 Asian specific variants, for leukocyte telomere length (LTL) associations among 25,533 Singapore Chinese samples. Three East Asian specific variants in/near POT1, TERF1 and STN1 genes are associated with LTL (Meta-analysis P 2.49×10-14-6.94×10-10). Rs79314063, a missense variant (p.Asp410His) at POT1, shows effect 5.3 fold higher and independent of a previous common index SNP. TERF1 (rs79617270) and STN1 (rs139620151) are linked to LTL-associated common index SNPs at these loci. Rs79617270 is associated with cancer mortality [HR95%CI = 1.544 (1.173, 2.032), PAdj = 0.018] and 4.76% of the association between the rs79617270 and colon cancer is mediated through LTL. Overall, genetically determined LTL is particularly associated with lung adenocarcinoma [HR95%CI = 1.123 (1.051, 1.201), Padj = 0.007]. Ethnicity-specific low frequency variants may affect LTL homeostasis and associate with certain cancers.

Association of Rare CYP39A1 Variants With Exfoliation Syndrome Involving the Anterior Chamber of the Eye.

Importance: Exfoliation syndrome is a systemic disorder characterized by progressive accumulation of abnormal fibrillar protein aggregates manifesting clinically in the anterior chamber of the eye. This disorder is the most commonly known cause of glaucoma and a major cause of irreversible blindness. Objective: To determine if exfoliation syndrome is associated with rare, protein-changing variants predicted to impair protein function. Design, Setting, and Participants: A 2-stage, case-control, whole-exome sequencing association study with a discovery cohort and 2 independently ascertained validation cohorts. Study participants from 14 countries were enrolled between February 1999 and December 2019. The date of last clinical follow-up was December 2019. Affected individuals had exfoliation material on anterior segment structures of at least 1 eye as visualized by slit lamp examination. Unaffected individuals had no signs of exfoliation syndrome. Exposures: Rare, coding-sequence genetic variants predicted to be damaging by bioinformatic algorithms trained to recognize alterations that impair protein function. Main Outcomes and Measures: The primary outcome was the presence of exfoliation syndrome. Exome-wide significance for detected variants was defined as P < 2.5 × 10-6. The secondary outcomes included biochemical enzymatic assays and gene expression analyses. Results: The discovery cohort included 4028 participants with exfoliation syndrome (median age, 78 years [interquartile range, 73-83 years]; 2377 [59.0%] women) and 5638 participants without exfoliation syndrome (median age, 72 years [interquartile range, 65-78 years]; 3159 [56.0%] women). In the discovery cohort, persons with exfoliation syndrome, compared with those without exfoliation syndrome, were significantly more likely to carry damaging CYP39A1 variants (1.3% vs 0.30%, respectively; odds ratio, 3.55 [95% CI, 2.07-6.10]; P = 6.1 × 10-7). This outcome was validated in 2 independent cohorts. The first validation cohort included 2337 individuals with exfoliation syndrome (median age, 74 years; 1132 women; n = 1934 with demographic data) and 2813 individuals without exfoliation syndrome (median age, 72 years; 1287 women; n = 2421 with demographic data). The second validation cohort included 1663 individuals with exfoliation syndrome (median age, 75 years; 587 women; n = 1064 with demographic data) and 3962 individuals without exfoliation syndrome (median age, 74 years; 951 women; n = 1555 with demographic data). Of the individuals from both validation cohorts, 5.2% with exfoliation syndrome carried CYP39A1 damaging alleles vs 3.1% without exfoliation syndrome (odds ratio, 1.82 [95% CI, 1.47-2.26]; P < .001). Biochemical assays classified 34 of 42 damaging CYP39A1 alleles as functionally deficient (median reduction in enzymatic activity compared with wild-type CYP39A1, 94.4% [interquartile range, 78.7%-98.2%] for the 34 deficient variants). CYP39A1 transcript expression was 47% lower (95% CI, 30%-64% lower; P < .001) in ciliary body tissues from individuals with exfoliation syndrome compared with individuals without exfoliation syndrome. Conclusions and Relevance: In this whole-exome sequencing case-control study, presence of exfoliation syndrome was significantly associated with carriage of functionally deficient CYP39A1 sequence variants. Further research is needed to understand the clinical implications of these findings.

Identification of Risk Loci for Parkinson Disease in Asians and Comparison of Risk Between Asians and Europeans: A Genome-Wide Association Study.

Importance: Large-scale genome-wide association studies in the European population have identified 90 risk variants associated with Parkinson disease (PD); however, there are limited studies in the largest population worldwide (ie, Asian). Objectives: To identify novel genome-wide significant loci for PD in Asian individuals and to compare genetic risk between Asian and European cohorts. Design Setting, and Participants: Genome-wide association data generated from PD cases and controls in an Asian population (ie, Singapore/Malaysia, Hong Kong, Taiwan, mainland China, and South Korea) were collected from January 1, 2016, to December 31, 2018, as part of an ongoing study. Results were combined with inverse variance meta-analysis, and replication of top loci in European and Japanese samples was performed. Discovery samples of 31 575 individuals passing quality control of 35 994 recruited were used, with a greater than 90% participation rate. A replication cohort of 1 926 361 European-ancestry and 3509 Japanese samples was analyzed. Parkinson disease was diagnosed using UK Parkinson's Disease Society Brain Bank Criteria. Main Outcomes and Measures: Genotypes of common variants, association with disease status, and polygenic risk scores. Results: Of 31 575 samples identified, 6724 PD cases (mean [SD] age, 64.3 [10] years; age at onset, 58.8 [10.6] years; 3472 [53.2%] men) and 24 851 controls (age, 59.4 [11.4] years; 11 030 [45.0%] men) were analyzed in the discovery study. Eleven genome-wide significant loci were identified; 2 of these loci were novel (SV2C and WBSCR17) and 9 were previously found in Europeans. Replication in European-ancestry and Japanese samples showed robust association for SV2C (rs246814; odds ratio, 1.16; 95% CI, 1.11-1.21; P = 1.17 × 10-10 in meta-analysis of discovery and replication samples) but showed potential genetic heterogeneity at WBSCR17 (rs9638616; I2=67.1%; P = 3.40 × 10-3 for hetereogeneity). Polygenic risk score models including variants at these 11 loci were associated with a significant improvement in area under the curve over the model based on 78 European loci alone (63.1% vs 60.2%; P = 6.81 × 10-12). Conclusions and Relevance: This study identified 2 apparently novel gene loci and found 9 previously identified European loci to be associated with PD in this large, meta-genome-wide association study in a worldwide population of Asian individuals and reports similarities and differences in genetic risk factors between Asian and European individuals in the risk for PD. These findings may lead to improved stratification of Asian patients and controls based on polygenic risk scores. Our findings have potential academic and clinical importance for risk stratification and precision medicine in Asia.

An integrative model of pathway convergence in genetically heterogeneous blast crisis chronic myeloid leukemia.

Targeted therapies against the BCR-ABL1 kinase have revolutionized treatment of chronic phase (CP) chronic myeloid leukemia (CML). In contrast, management of blast crisis (BC) CML remains challenging because BC cells acquire complex molecular alterations that confer stemness features to progenitor populations and resistance to BCR-ABL1 tyrosine kinase inhibitors. Comprehensive models of BC transformation have proved elusive because of the rarity and genetic heterogeneity of BC, but are important for developing biomarkers predicting BC progression and effective therapies. To better understand BC, we performed an integrated multiomics analysis of 74 CP and BC samples using whole-genome and exome sequencing, transcriptome and methylome profiling, and chromatin immunoprecipitation followed by high-throughput sequencing. Employing pathway-based analysis, we found the BC genome was significantly enriched for mutations affecting components of the polycomb repressive complex (PRC) pathway. While transcriptomically, BC progenitors were enriched and depleted for PRC1- and PRC2-related gene sets respectively. By integrating our data sets, we determined that BC progenitors undergo PRC-driven epigenetic reprogramming toward a convergent transcriptomic state. Specifically, PRC2 directs BC DNA hypermethylation, which in turn silences key genes involved in myeloid differentiation and tumor suppressor function via so-called epigenetic switching, whereas PRC1 represses an overlapping and distinct set of genes, including novel BC tumor suppressors. On the basis of these observations, we developed an integrated model of BC that facilitated the identification of combinatorial therapies capable of reversing BC reprogramming (decitabine+PRC1 inhibitors), novel PRC-silenced tumor suppressor genes (NR4A2), and gene expression signatures predictive of disease progression and drug resistance in CP.

Loci for human leukocyte telomere length in the Singaporean Chinese population and trans-ethnic genetic studies.

Genetic factors underlying leukocyte telomere length (LTL) may provide insights into telomere homeostasis, with direct links to disease susceptibility. Genetic evaluation of 23,096 Singaporean Chinese samples identifies 10 genome-wide loci (P < 5 × 10-8). Several of these contain candidate genes (TINF2, PARP1, TERF1, ATM and POT1) with potential roles in telomere biology and DNA repair mechanisms. Meta-analyses with additional 37,505 European individuals reveals six more genome-wide loci, including associations at MPHOSPH6, NKX2-3 and TYMS. We demonstrate that longer LTL associates with protection against respiratory disease mortality [HR = 0.854(0.804-0.906), P = 1.88 × 10-7] in the Singaporean Chinese samples. We further show that the LTL reducing SNP rs7253490 associates with respiratory infections (P = 7.44 × 10-4) although this effect may not be strongly mediated through LTL. Our data expands on the genetic basis of LTL and may indicate on a potential role of LTL in immune competence.

Genome-Wide Analysis of Protein-Coding Variants in Leprosy.

Although genome-wide association studies have greatly advanced our understanding of the contribution of common noncoding variants to leprosy susceptibility, protein-coding variants have not been systematically investigated. We carried out a three-stage genome-wide association study of protein-coding variants in Han Chinese, of whom were 7,048 leprosy patients and 14,398 were healthy control subjects. Seven coding variants of exome-wide significance were discovered, including two rare variants: rs145562243 in NCKIPSD (P = 1.71 × 10-9, odds ratio [OR] = 4.35) and rs149308743 in CARD9 (P = 2.09 × 10-8, OR = 4.75); three low-frequency variants: rs76418789 in IL23R (P = 1.03 × 10-10, OR = 1.36), rs146466242 in FLG (P = 3.39 × 10-12, OR = 1.45), and rs55882956 in TYK2 (P = 1.04 × 10-6, OR = 1.30); and two common variants: rs780668 in SLC29A3 (P = 2.17 × 10-9, OR = 1.14) and rs181206 in IL27 (P = 1.08 × 10-7, OR = 0.83). Discovered protein-coding variants, particularly low-frequency and rare ones, showed involvement of skin barrier and endocytosis/phagocytosis/autophagy, in addition to known innate and adaptive immunity, in the pathogenesis of leprosy, highlighting the merits of protein-coding variant studies for complex diseases.

Genome-wide association study of Parkinson's disease in East Asians.

Genome-wide association studies (GWAS) on Parkinson's disease (PD) have mostly been done in Europeans and Japanese. No study has been done in Han Chinese, which make up nearly a fifth of the world population. We conducted the first Han Chinese GWAS analysing a total of 22,729 subjects (5,125 PD cases and 17,604 controls) from Singapore, Hong Kong, Malaysia, Korea, mainland China and Taiwan. We performed imputation, merging and logistic regression analyses of 2,402,394 SNPs passing quality control filters in 779 PD cases, 13,227 controls, adjusted for the first three principal components. 90 SNPs with association P < 10-4 were validated in 9 additional sample collections and the results were combined using fixed-effects inverse-variance meta-analysis. We observed strong associations reaching genome-wide significance at SNCA, LRRK2 and MCCC1, confirming their important roles in both European and Asian PD. We also identified significant (P < 0.05) associations at 5 loci (DLG2, SIPA1L2, STK39, VPS13C and RIT2), and observed the same direction of associations at 9 other loci including BST1 and PARK16. Allelic heterogeneity was observed at LRRK2 while European risk SNPs at 6 other loci including MAPT and GBA-SYT11 were non-polymorphic or very rare in our cohort. Overall, we replicate associations at SNCA, LRRK2, MCCC1 and 14 other European PD loci but did not identify Asian-specific loci with large effects (OR > 1.45) on PD risk. Our results also demonstrate some differences in the genetic contribution to PD between Europeans and Asians. Further pan-ethnic meta-analysis with European GWAS cohorts may unravel new PD loci.

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

BACKGROUND: Extranodal natural killer T-cell lymphoma (NKTCL), nasal type, is a rare and aggressive malignancy that occurs predominantly in Asian and Latin American populations. Although Epstein-Barr virus infection is a known risk factor, other risk factors and the pathogenesis of NKTCL are not well understood. We aimed to identify common genetic variants affecting individual risk of NKTCL. METHODS: We did a genome-wide association study of 189 patients with extranodal NKTCL, nasal type (WHO classification criteria; cases) and 957 controls from Guangdong province, southern China. We validated our findings in four independent case-control series, including 75 cases from Guangdong province and 296 controls from Hong Kong, 65 cases and 983 controls from Guangdong province, 125 cases and 1110 controls from Beijing (northern China), and 60 cases and 2476 controls from Singapore. We used imputation and conditional logistic regression analyses to fine-map the associations. We also did a meta-analysis of the replication series and of the entire dataset. FINDINGS: Associations exceeding the genome-wide significance threshold (p<5 × 10(-8)) were seen at 51 single-nucleotide polymorphisms (SNPs) mapping to the class II MHC region on chromosome 6, with rs9277378 (located in HLA-DPB1) having the strongest association with NKTCL susceptibility (p=4·21 × 10(-19), odds ratio [OR] 1·84 [95% CI 1·61-2·11] in meta-analysis of entire dataset). Imputation-based fine-mapping across the class II MHC region suggests that four aminoacid residues (Gly84-Gly85-Pro86-Met87) in near-complete linkage disequilibrium at the edge of the peptide-binding groove of HLA-DPB1 could account for most of the association between the rs9277378*A risk allele and NKTCL susceptibility (OR 2·38, p value for haplotype 2·32 × 10(-14)). This association is distinct from MHC associations with Epstein-Barr virus infection. INTERPRETATION: To our knowledge, this is the first time that a genetic variant conferring an NKTCL risk is noted at genome-wide significance. This finding underlines the importance of HLA-DP antigen presentation in the pathogenesis of NKTCL. FUNDING: Top-Notch Young Talents Program of China, Special Support Program of Guangdong, Specialized Research Fund for the Doctoral Program of Higher Education (20110171120099), Program for New Century Excellent Talents in University (NCET-11-0529), National Medical Research Council of Singapore (TCR12DEC005), Tanoto Foundation Professorship in Medical Oncology, New Century Foundation Limited, Ling Foundation, Singapore National Cancer Centre Research Fund, and the US National Institutes of Health (1R01AR062886, 5U01GM092691-04, and 1R01AR063759-01A1).