Direct inference and control of genetic population structure from RNA sequencing data.

RNAseq data can be used to infer genetic variants, yet its use for estimating genetic population structure remains underexplored. Here, we construct a freely available computational tool (RGStraP) to estimate RNAseq-based genetic principal components (RG-PCs) and assess whether RG-PCs can be used to control for population structure in gene expression analyses. Using whole blood samples from understudied Nepalese populations and the Geuvadis study, we show that RG-PCs had comparable results to paired array-based genotypes, with high genotype concordance and high correlations of genetic principal components, capturing subpopulations within the dataset. In differential gene expression analysis, we found that inclusion of RG-PCs as covariates reduced test statistic inflation. Our paper demonstrates that genetic population structure can be directly inferred and controlled for using RNAseq data, thus facilitating improved retrospective and future analyses of transcriptomic data.

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.

Cross-Ancestry Genome-Wide Association Study Defines the Extended CYP2D6 Locus as the Principal Genetic Determinant of Endoxifen Plasma Concentrations.

The therapeutic efficacy of tamoxifen is predominantly mediated by its active metabolites 4-hydroxy-tamoxifen and endoxifen, whose formation is catalyzed by the polymorphic cytochrome P450 2D6 (CYP2D6). Yet, known CYP2D6 polymorphisms only partially determine metabolite concentrations in vivo. We performed the first cross-ancestry genome-wide association study with well-characterized patients of European, Middle-Eastern, and Asian descent (n = 497) to identify genetic factors impacting active and parent metabolite formation. Genome-wide significant variants were functionally evaluated in an independent liver cohort (n = 149) and in silico. Metabolite prediction models were validated in two independent European breast cancer cohorts (n = 287, n = 189). Within a single 1-megabase (Mb) region of chromosome 22q13 encompassing the CYP2D6 gene, 589 variants were significantly associated with tamoxifen metabolite concentrations, particularly endoxifen and metabolic ratio (MR) endoxifen/N-desmethyltamoxifen (minimal P = 5.4E-35 and 2.5E-65, respectively). Previously suggested other loci were not confirmed. Functional analyses revealed 66% of associated, mostly intergenic variants to be significantly correlated with hepatic CYP2D6 activity or expression (ρ = 0.35 to -0.52), and six hotspot regions in the extended 22q13 locus impacting gene regulatory function. Machine learning models based on hotspot variants (n = 12) plus CYP2D6 activity score (AS) increased the explained variability (~ 9%) compared with AS alone, explaining up to 49% (median R2 ) and 72% of the variability in endoxifen and MR endoxifen/N-desmethyltamoxifen, respectively. Our findings suggest that the extended CYP2D6 locus at 22q13 is the principal genetic determinant of endoxifen plasma concentration. Long-distance haplotypes connecting CYP2D6 with adjacent regulatory sites and nongenetic factors may account for the unexplained portion of variability.

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.

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.

Natural resistance to Meningococcal Disease related to CFH loci: Meta-analysis of genome-wide association studies.

Meningococcal disease (MD) remains an important infectious cause of life threatening infection in both industrialized and resource poor countries. Genetic factors influence both occurrence and severity of presentation, but the genes responsible are largely unknown. We performed a genome-wide association study (GWAS) examining 5,440,063 SNPs in 422 Spanish MD patients and 910 controls. We then performed a meta-analysis of the Spanish GWAS with GWAS data from the United Kingdom (combined cohorts: 897 cases and 5,613 controls; 4,898,259 SNPs). The meta-analysis identified strong evidence of association (P-value ≤ 5 × 10-8) in 20 variants located at the CFH gene. SNP rs193053835 showed the most significant protective effect (Odds Ratio (OR) = 0.62, 95% confidence interval (C.I.) = 0.52-0.73; P-value = 9.62 × 10-9). Five other variants had been previously reported to be associated with susceptibility to MD, including the missense SNP rs1065489 (OR = 0.64, 95% C.I.) = 0.55-0.76, P-value = 3.25 × 10-8). Theoretical predictions point to a functional effect of rs1065489, which may be directly responsible for protection against MD. Our study confirms the association of CFH with susceptibility to MD and strengthens the importance of this link in understanding pathogenesis of the disease.

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).

A missense variant in FGD6 confers increased risk of polypoidal choroidal vasculopathy.

Polypoidal choroidal vasculopathy (PCV), a subtype of 'wet' age-related macular degeneration (AMD), constitutes up to 55% of cases of wet AMD in Asian patients. In contrast to the choroidal neovascularization (CNV) subtype, the genetic risk factors for PCV are relatively unknown. Exome sequencing analysis of a Han Chinese cohort followed by replication in four independent cohorts identified a rare c.986A>G (p.Lys329Arg) variant in the FGD6 gene as significantly associated with PCV (P = 2.19 × 10(-16), odds ratio (OR) = 2.12) but not with CNV (P = 0.26, OR = 1.13). The intracellular localization of FGD6-Arg329 is distinct from that of FGD6-Lys329. In vitro, FGD6 could regulate proangiogenic activity, and oxidized phospholipids increased expression of FGD6. FGD6-Arg329 promoted more abnormal vessel development in the mouse retina than FGD6-Lys329. Collectively, our data suggest that oxidized phospholipids and FGD6-Arg329 might act synergistically to increase susceptibility to PCV.

A common variant mapping to CACNA1A is associated with susceptibility to exfoliation syndrome.

Exfoliation syndrome (XFS) is the most common recognizable cause of open-angle glaucoma worldwide. To better understand the etiology of XFS, we conducted a genome-wide association study (GWAS) of 1,484 cases and 1,188 controls from Japan and followed up the most significant findings in a further 6,901 cases and 20,727 controls from 17 countries across 6 continents. We discovered a genome-wide significant association between a new locus (CACNA1A rs4926244) and increased susceptibility to XFS (odds ratio (OR) = 1.16, P = 3.36 × 10(-11)). Although we also confirmed overwhelming association at the LOXL1 locus, the key SNP marker (LOXL1 rs4886776) demonstrated allelic reversal depending on the ancestry group (Japanese: OR(A allele) = 9.87, P = 2.13 × 10(-217); non-Japanese: OR(A allele) = 0.49, P = 2.35 × 10(-31)). Our findings represent the first genetic locus outside of LOXL1 surpassing genome-wide significance for XFS and provide insight into the biology and pathogenesis of the disease.

Meta-analysis of Genome-Wide Association Studies Identifies Novel Loci Associated With Optic Disc Morphology.

Primary open-angle glaucoma is the most common optic neuropathy and an important cause of irreversible blindness worldwide. The optic nerve head or optic disc is divided in two parts: a central cup (without nerve fibers) surrounded by the neuroretinal rim (containing axons of the retinal ganglion cells). The International Glaucoma Genetics Consortium conducted a meta-analysis of genome-wide association studies consisting of 17,248 individuals of European ancestry and 6,841 individuals of Asian ancestry. The outcomes of the genome-wide association studies were disc area and cup area. These specific measurements describe optic nerve morphology in another way than the vertical cup-disc ratio, which is a clinically used measurement, and may shed light on new glaucoma mechanisms. We identified 10 new loci associated with disc area (CDC42BPA, F5, DIRC3, RARB, ABI3BP, DCAF4L2, ELP4, TMTC2, NR2F2, and HORMAD2) and another 10 new loci associated with cup area (DHRS3, TRIB2, EFEMP1, FLNB, FAM101, DDHD1, ASB7, KPNB1, BCAS3, and TRIOBP). The new genes participate in a number of pathways and future work is likely to identify more functions related to the pathogenesis of glaucoma.

New loci and coding variants confer risk for age-related macular degeneration in East Asians.

Age-related macular degeneration (AMD) is a major cause of blindness, but presents differently in Europeans and Asians. Here, we perform a genome-wide and exome-wide association study on 2,119 patients with exudative AMD and 5,691 controls, with independent replication in 4,226 patients and 10,289 controls, all of East Asian descent, as part of The Genetics of AMD in Asians (GAMA) Consortium. We find a strong association between CETP Asp442Gly (rs2303790), an East Asian-specific mutation, and increased risk of AMD (odds ratio (OR)=1.70, P=5.60 × 10(-22)). The AMD risk allele (442Gly), known to protect from coronary heart disease, increases HDL cholesterol levels by 0.17 mmol l(-1) (P=5.82 × 10(-21)) in East Asians (n=7,102). We also identify three novel AMD loci: C6orf223 Ala231Ala (OR=0.78, P=6.19 × 10(-18)), SLC44A4 Asp47Val (OR=1.27, P=1.08 × 10(-11)) and FGD6 Gln257Arg (OR=0.87, P=2.85 × 10(-8)). Our findings suggest that some of the genetic loci conferring AMD susceptibility in East Asians are shared with Europeans, yet AMD in East Asians may also have a distinct genetic signature.

Variation at HLA-DRB1 is associated with resistance to enteric fever.

Enteric fever affects more than 25 million people annually and results from systemic infection with Salmonella enterica serovar Typhi or Paratyphi pathovars A, B or C(1). We conducted a genome-wide association study of 432 individuals with blood culture-confirmed enteric fever and 2,011 controls from Vietnam. We observed strong association at rs7765379 (odds ratio (OR) for the minor allele = 0.18, P = 4.5 × 10(-10)), a marker mapping to the HLA class II region, in proximity to HLA-DQB1 and HLA-DRB1. We replicated this association in 595 enteric fever cases and 386 controls from Nepal and also in a second independent collection of 151 cases and 668 controls from Vietnam. Imputation-based fine-mapping across the extended MHC region showed that the classical HLA-DRB1*04:05 allele (OR = 0.14, P = 2.60 × 10(-11)) could entirely explain the association at rs7765379, thus implicating HLA-DRB1 as a major contributor to resistance against enteric fever, presumably through antigen presentation.

Meta-analysis of genome-wide association studies identifies novel loci that influence cupping and the glaucomatous process.

Glaucoma is characterized by irreversible optic nerve degeneration and is the most frequent cause of irreversible blindness worldwide. Here, the International Glaucoma Genetics Consortium conducts a meta-analysis of genome-wide association studies of vertical cup-disc ratio (VCDR), an important disease-related optic nerve parameter. In 21,094 individuals of European ancestry and 6,784 individuals of Asian ancestry, we identify 10 new loci associated with variation in VCDR. In a separate risk-score analysis of five case-control studies, Caucasians in the highest quintile have a 2.5-fold increased risk of primary open-angle glaucoma as compared with those in the lowest quintile. This study has more than doubled the known loci associated with optic disc cupping and will allow greater understanding of mechanisms involved in this common blinding condition.

Deep whole-genome sequencing of 100 southeast Asian Malays.

Whole-genome sequencing across multiple samples in a population provides an unprecedented opportunity for comprehensively characterizing the polymorphic variants in the population. Although the 1000 Genomes Project (1KGP) has offered brief insights into the value of population-level sequencing, the low coverage has compromised the ability to confidently detect rare and low-frequency variants. In addition, the composition of populations in the 1KGP is not complete, despite the fact that the study design has been extended to more than 2,500 samples from more than 20 population groups. The Malays are one of the Austronesian groups predominantly present in Southeast Asia and Oceania, and the Singapore Sequencing Malay Project (SSMP) aims to perform deep whole-genome sequencing of 100 healthy Malays. By sequencing at a minimum of 30× coverage, we have illustrated the higher sensitivity at detecting low-frequency and rare variants and the ability to investigate the presence of hotspots of functional mutations. Compared to the low-pass sequencing in the 1KGP, the deeper coverage allows more functional variants to be identified for each person. A comparison of the fidelity of genotype imputation of Malays indicated that a population-specific reference panel, such as the SSMP, outperforms a cosmopolitan panel with larger number of individuals for common SNPs. For lower-frequency (<5%) markers, a larger number of individuals might have to be whole-genome sequenced so that the accuracy currently afforded by the 1KGP can be achieved. The SSMP data are expected to be the benchmark for evaluating the value of deep population-level sequencing versus low-pass sequencing, especially in populations that are poorly represented in population-genetics studies.

Genome-wide association analyses identify three new susceptibility loci for primary angle closure glaucoma.

Primary angle closure glaucoma (PACG) is a major cause of blindness worldwide. We conducted a genome-wide association study including 1,854 PACG cases and 9,608 controls across 5 sample collections in Asia. Replication experiments were conducted in 1,917 PACG cases and 8,943 controls collected from a further 6 sample collections. We report significant associations at three new loci: rs11024102 in PLEKHA7 (per-allele odds ratio (OR)=1.22; P=5.33×10(-12)), rs3753841 in COL11A1 (per-allele OR=1.20; P=9.22×10(-10)) and rs1015213 located between PCMTD1 and ST18 on chromosome 8q (per-allele OR=1.50; P=3.29×10(-9)). Our findings, accumulated across these independent worldwide collections, suggest possible mechanisms explaining the pathogenesis of PACG.

Natural positive selection and north-south genetic diversity in East Asia.

Recent reports have identified a north-south cline in genetic variation in East and South-East Asia, but these studies have not formally explored the basis of these clinical differences. Understanding the origins of these variations may provide valuable insights in tracking down the functional variants in genomic regions identified by genetic association studies. Here we investigate the genetic basis of these differences with genome-wide data from the HapMap, the Human Genome Diversity Project and the Singapore Genome Variation Project. We implemented four bioinformatic measures to discover genomic regions that are considerably differentiated either between two Han Chinese populations in the north and south of China, or across 22 populations in East and South-East Asia. These measures prioritized genomic stretches with: (i) regional differences in the allelic spectrum for SNPs common to the two Han Chinese populations; (ii) differential evidence of positive selection between the two populations as quantified by integrated haplotype score (iHS) and cross-population extended haplotype homozygosity (XP-EHH); (iii) significant correlation between allele frequencies and geographical latitudes of the 22 populations. We also explored the extent of linkage disequilibrium variations in these regions, which is important in combining genetic association studies from North and South Chinese. Two of the regions that emerged are found in HLA class I and II, suggesting that the HLA imputation panel from the HapMap may not be directly applicable to every Chinese sample. This has important implications to autoimmune studies that plan to impute the classical HLA alleles to fine map the SNP association signals.