Search

Large-scale gene sequencing studies for complex traits have the potential to identify causal genes with therapeutic implications. We performed gene-based association testing of blood lipid levels with rare (minor allele frequency < 1%) predicted damaging coding variation by using sequence data from >170,000 individuals from multiple ancestries: 97,493 European, 30,025 South Asian, 16,507 African, 16,440 Hispanic/Latino, 10,420 East Asian, and 1,182 Samoan. We identified 35 genes associated with circulating lipid levels; some of these genes have not been previously associated with lipid levels when using rare coding variation from population-based samples. We prioritize 32 genes in array-based genome-wide association study (GWAS) loci based on aggregations of rare coding variants; three (EVI5, SH2B3, and PLIN1) had no prior association of rare coding variants with lipid levels. Most of our associated genes showed evidence of association among multiple ancestries. Finally, we observed an enrichment of gene-based associations for low-density lipoprotein cholesterol drug target genes and for genes closest to GWAS index single-nucleotide polymorphisms (SNPs). Our results demonstrate that gene-based associations can be beneficial for drug target development and provide evidence that the gene closest to the array-based GWAS index SNP is often the functional gene for blood lipid levels., Competing Interests: Declaration of interests The authors declare no competing interests for the present work. P.N. reports investigator-initiated grants from Amgen, Apple, and Boston Scientific; is a scientific advisor to Apple, Blackstone Life Sciences, and Novartis; and has spousal employment at Vertex, all unrelated to the present work. A.V.K. has served as a scientific advisor to Sanofi, Medicines Company, Maze Pharmaceuticals, Navitor Pharmaceuticals, Verve Therapeutics, Amgen, and Color; received speaking fees from Illumina, MedGenome, Amgen, and the Novartis Institute for Biomedical Research; received sponsored research agreements from the Novartis Institute for Biomedical Research and IBM Research; and reports a patent related to a genetic risk predictor (20190017119). C.J.W.’s spouse is employed at Regeneron. L.E.S. is currently an employee of Celgene/Bristol Myers Squibb. Celgene/Bristol Myers Squibb had no role in the funding, design, conduct, and interpretation of this study. M.E.M. receives funding from Regeneron unrelated to this work. E.E.K. has received speaker honoraria from Illumina, Inc and Regeneron Pharmaceuticals. B.M.P. serves on the Steering Committee of the Yale Open Data Access Project funded by Johnson & Johnson. L.A.C. has consulted with the Dyslipidemia Foundation on lipid projects in the Framingham Heart Study. P.T.E. is supported by a grant from Bayer AG to the Broad Institute focused on the genetics and therapeutics of cardiovascular disease. P.T.E. has consulted for Bayer AG, Novartis, MyoKardia, and Quest Diagnostics. S.A.L. receives sponsored research support from Bristol Myers Squibb/Pfizer, Bayer AG, Boehringer Ingelheim, Fitbit, and IBM and has consulted for Bristol Myers Squibb/Pfizer, Bayer AG, and Blackstone Life Sciences. The views expressed in this article are those of the author(s) and not necessarily those of the NHS, the NIHR, or the Department of Health. M.I.M. has served on advisory panels for Pfizer, NovoNordisk, and Zoe Global and has received honoraria from Merck, Pfizer, Novo Nordisk, and Eli Lilly and research funding from Abbvie, Astra Zeneca, Boehringer Ingelheim, Eli Lilly, Janssen, Merck, NovoNordisk, Pfizer, Roche, Sanofi Aventis, Servier, and Takeda. As of June 2019, M.I.M. is an employee of Genentech and a holder of Roche stock. M.E.J. holds shares in Novo Nordisk A/S. H.M.K. is an employee of Regeneron Pharmaceuticals; he owns stock and stock options for Regeneron Pharmaceuticals. M.E.J. has received research grants form Astra Zeneca, Boehringer Ingelheim, Amgen, and Sanofi. S.K. is founder of Verve Therapeutics., (Copyright © 2021 American Society of Human Genetics. All rights reserved.)

Background: Hypertriglyceridemia has emerged as a critical coronary artery disease (CAD) risk factor. Rare loss-of-function (LoF) variants in apolipoprotein C-III have been reported to reduce triglycerides (TG) and are cardioprotective in American Indians and Europeans. However, there is a lack of data in other Europeans and non-Europeans. Also, whether genetically increased plasma TG due to ApoC-III is causally associated with increased CAD risk is still unclear and inconsistent. The objectives of this study were to verify the cardioprotective role of earlier reported six LoF variants of APOC3 in South Asians and other multi-ethnic cohorts and to evaluate the causal association of TG raising common variants for increasing CAD risk., Methods: We performed gene-centric and Mendelian randomization analyses and evaluated the role of genetic variation encompassing APOC3 for affecting circulating TG and the risk for developing CAD., Results: One rare LoF variant (rs138326449) with a 37% reduction in TG was associated with lowered risk for CAD in Europeans (p = 0.007), but we could not confirm this association in Asian Indians (p = 0.641). Our data could not validate the cardioprotective role of other five LoF variants analysed. A common variant rs5128 in the APOC3 was strongly associated with elevated TG levels showing a p-value 2.8 × 10 - 424 . Measures of plasma ApoC-III in a small subset of Sikhs revealed a 37% increase in ApoC-III concentrations among homozygous mutant carriers than the wild-type carriers of rs5128. A genetically instrumented per 1SD increment of plasma TG level of 15 mg/dL would cause a mild increase (3%) in the risk for CAD (p = 0.042)., Conclusions: Our results highlight the challenges of inclusion of rare variant information in clinical risk assessment and the generalizability of implementation of ApoC-III inhibition for treating atherosclerotic disease. More studies would be needed to confirm whether genetically raised TG and ApoC-III concentrations would increase CAD risk., (© 2021. The Author(s).)

Context: Elevated levels of plasma leucine-rich α-2-glycoprotein 1 (LRG1), a component of transforming growth factor beta signaling, are associated with development and progression of chronic kidney disease in patients with type 2 diabetes (T2D). However, whether this relationship is causal is uncertain., Objectives: To identify genetic variants associated with plasma LRG1 levels and determine whether genetically predicted plasma LRG1 contributes to a rapid decline in kidney function (RDKF) in patients with T2D., Design and Participants: We performed a genome-wide association study of plasma LRG1 among 3694 T2D individuals [1881 (983 Chinese, 420 Malay, and 478 Indian) discovery from Singapore Study of Macro-angiopathy and Micro-vascular Reactivity in Type 2 Diabetes cohort and 1813 (Chinese) validation from Diabetic Nephropathy cohort]. One- sample Mendelian randomization analysis was performed among 1337 T2D Chinese participants with preserved glomerular filtration function [baseline estimated glomerular filtration rate (eGFR) ≥60 mL/min/1.73 m2)]. RDKF was defined as an eGFR decline of 3 mL/min/1.73 m2/year or greater., Results: We identified rs4806985 variant near LRG1 locus robustly associated with plasma LRG1 levels (meta P = 6.66 × 10-16). Among 1337 participants, 344 (26%) developed RDKF, and the rs4806985 variant was associated with higher odds of RDKF (meta odds ratio = 1.23, P = 0.030 adjusted for age and sex). Mendelian randomization analysis provided evidence for a potential causal effect of plasma LRG1 on kidney function decline in T2D (P < 0.05)., Conclusion: We demonstrate that genetically influenced plasma LRG1 increases the risk of RDKF in T2D patients, suggesting plasma LRG1 as a potential treatment target. However, further studies are warranted to elucidate underlying pathways to provide insight into diabetic kidney disease prevention., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Endocrine Society. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Nuclear cataract is the most common type of age-related cataract and a leading cause of blindness worldwide. Age-related nuclear cataract is heritable (h 2  = 0.48), but little is known about specific genetic factors underlying this condition. Here we report findings from the largest to date multi-ethnic meta-analysis of genome-wide association studies (discovery cohort N = 14,151 and replication N = 5299) of the International Cataract Genetics Consortium. We confirmed the known genetic association of CRYAA (rs7278468, P = 2.8 × 10 -16 ) with nuclear cataract and identified five new loci associated with this disease: SOX2-OT (rs9842371, P = 1.7 × 10 -19 ), TMPRSS5 (rs4936279, P = 2.5 × 10 -10 ), LINC01412 (rs16823886, P = 1.3 × 10 -9 ), GLTSCR1 (rs1005911, P = 9.8 × 10 -9 ), and COMMD1 (rs62149908, P = 1.2 × 10 -8 ). The results suggest a strong link of age-related nuclear cataract with congenital cataract and eye development genes, and the importance of common genetic variants in maintaining crystalline lens integrity in the aging eye.

Purpose: To evaluate the roles of known myopia-associated genetic variants for development of myopic macular degeneration (MMD) in individuals with high myopia (HM), using case-control studies from the Consortium of Refractive Error and Myopia (CREAM)., Methods: A candidate gene approach tested 50 myopia-associated loci for association with HM and MMD, using meta-analyses of case-control studies comprising subjects of European and Asian ancestry aged 30 to 80 years from 10 studies. Fifty loci with the strongest associations with myopia were chosen from a previous published GWAS study. Highly myopic (spherical equivalent [SE] ≤ -5.0 diopters [D]) cases with MMD (N = 348), and two sets of controls were enrolled: (1) the first set included 16,275 emmetropes (SE ≤ -0.5 D); and (2) second set included 898 highly myopic subjects (SE ≤ -5.0 D) without MMD. MMD was classified based on the International photographic classification for pathologic myopia (META-PM)., Results: In the first analysis, comprising highly myopic cases with MMD (N = 348) versus emmetropic controls without MMD (N = 16,275), two SNPs were significantly associated with high myopia in adults with HM and MMD: (1) rs10824518 (P = 6.20E-07) in KCNMA1, which is highly expressed in human retinal and scleral tissues; and (2) rs524952 (P = 2.32E-16) near GJD2. In the second analysis, comprising highly myopic cases with MMD (N = 348) versus highly myopic controls without MMD (N = 898), none of the SNPs studied reached Bonferroni-corrected significance., Conclusions: Of the 50 myopia-associated loci, we did not find any variant specifically associated with MMD, but the KCNMA1 and GJD2 loci were significantly associated with HM in highly myopic subjects with MMD, compared to emmetropes., Competing Interests: YLW is an employee of Essilor International, Singapore and received funding in the form of salary. Author KOM is a consultant for Bayer, Nevakah, and Santen, and receiving support in the form of consultancy fees. This does not alter our adherence to PLOS ONE policies on sharing data and materials. There are no patents, products in development or marketed products associated with this research to declare.

Body-fat distribution is a risk factor for adverse cardiovascular health consequences. We analyzed the association of body-fat distribution, assessed by waist-to-hip ratio adjusted for body mass index, with 228,985 predicted coding and splice site variants available on exome arrays in up to 344,369 individuals from five major ancestries (discovery) and 132,177 European-ancestry individuals (validation). We identified 15 common (minor allele frequency, MAF ≥5%) and nine low-frequency or rare (MAF <5%) coding novel variants. Pathway/gene set enrichment analyses identified lipid particle, adiponectin, abnormal white adipose tissue physiology and bone development and morphology as important contributors to fat distribution, while cross-trait associations highlight cardiometabolic traits. In functional follow-up analyses, specifically in Drosophila RNAi-knockdowns, we observed a significant increase in the total body triglyceride levels for two genes (DNAH10 and PLXND1). We implicate novel genes in fat distribution, stressing the importance of interrogating low-frequency and protein-coding variants.

A major challenge in evaluating the contribution of rare variants to complex disease is identifying enough copies of the rare alleles to permit informative statistical analysis. To investigate the contribution of rare variants to the risk of type 2 diabetes (T2D) and related traits, we performed deep whole-genome analysis of 1,034 members of 20 large Mexican-American families with high prevalence of T2D. If rare variants of large effect accounted for much of the diabetes risk in these families, our experiment was powered to detect association. Using gene expression data on 21,677 transcripts for 643 pedigree members, we identified evidence for large-effect rare-variant cis -expression quantitative trait loci that could not be detected in population studies, validating our approach. However, we did not identify any rare variants of large effect associated with T2D, or the related traits of fasting glucose and insulin, suggesting that large-effect rare variants account for only a modest fraction of the genetic risk of these traits in this sample of families. Reliable identification of large-effect rare variants will require larger samples of extended pedigrees or different study designs that further enrich for such variants., Competing Interests: Conflict of interest statement: S.E.L., J. Laramie, and R.G.T. were employees of Complete Genomics during this study. T.M.T. is an employee of Regeneron Pharmaceuticals. D.A. is an employee of Vertex Pharmaceuticals.

Genome-wide association studies (GWAS) have identified >250 loci for body mass index (BMI), implicating pathways related to neuronal biology. Most GWAS loci represent clusters of common, noncoding variants from which pinpointing causal genes remains challenging. Here we combined data from 718,734 individuals to discover rare and low-frequency (minor allele frequency (MAF) < 5%) coding variants associated with BMI. We identified 14 coding variants in 13 genes, of which 8 variants were in genes (ZBTB7B, ACHE, RAPGEF3, RAB21, ZFHX3, ENTPD6, ZFR2 and ZNF169) newly implicated in human obesity, 2 variants were in genes (MC4R and KSR2) previously observed to be mutated in extreme obesity and 2 variants were in GIPR. The effect sizes of rare variants are ~10 times larger than those of common variants, with the largest effect observed in carriers of an MC4R mutation introducing a stop codon (p.Tyr35Ter, MAF = 0.01%), who weighed ~7 kg more than non-carriers. Pathway analyses based on the variants associated with BMI confirm enrichment of neuronal genes and provide new evidence for adipocyte and energy expenditure biology, widening the potential of genetically supported therapeutic targets in obesity.

To investigate the genetic basis of type 2 diabetes (T2D) to high resolution, the GoT2D and T2D-GENES consortia catalogued variation from whole-genome sequencing of 2,657 European individuals and exome sequencing of 12,940 individuals of multiple ancestries. Over 27M SNPs, indels, and structural variants were identified, including 99% of low-frequency (minor allele frequency [MAF] 0.1-5%) non-coding variants in the whole-genome sequenced individuals and 99.7% of low-frequency coding variants in the whole-exome sequenced individuals. Each variant was tested for association with T2D in the sequenced individuals, and, to increase power, most were tested in larger numbers of individuals (>80% of low-frequency coding variants in ~82 K Europeans via the exome chip, and ~90% of low-frequency non-coding variants in ~44 K Europeans via genotype imputation). The variants, genotypes, and association statistics from these analyses provide the largest reference to date of human genetic information relevant to T2D, for use in activities such as T2D-focused genotype imputation, functional characterization of variants or genes, and other novel analyses to detect associations between sequence variation and T2D.

We screened variants on an exome-focused genotyping array in >300,000 participants (replication in >280,000 participants) and identified 444 independent variants in 250 loci significantly associated with total cholesterol (TC), high-density-lipoprotein cholesterol (HDL-C), low-density-lipoprotein cholesterol (LDL-C), and/or triglycerides (TG). At two loci (JAK2 and A1CF), experimental analysis in mice showed lipid changes consistent with the human data. We also found that: (i) beta-thalassemia trait carriers displayed lower TC and were protected from coronary artery disease (CAD); (ii) excluding the CETP locus, there was not a predictable relationship between plasma HDL-C and risk for age-related macular degeneration; (iii) only some mechanisms of lowering LDL-C appeared to increase risk for type 2 diabetes (T2D); and (iv) TG-lowering alleles involved in hepatic production of TG-rich lipoproteins (TM6SF2 and PNPLA3) tracked with higher liver fat, higher risk for T2D, and lower risk for CAD, whereas TG-lowering alleles involved in peripheral lipolysis (LPL and ANGPTL4) had no effect on liver fat but decreased risks for both T2D and CAD.

Background: Glycated hemoglobin (HbA1c) is used to diagnose type 2 diabetes (T2D) and assess glycemic control in patients with diabetes. Previous genome-wide association studies (GWAS) have identified 18 HbA1c-associated genetic variants. These variants proved to be classifiable by their likely biological action as erythrocytic (also associated with erythrocyte traits) or glycemic (associated with other glucose-related traits). In this study, we tested the hypotheses that, in a very large scale GWAS, we would identify more genetic variants associated with HbA1c and that HbA1c variants implicated in erythrocytic biology would affect the diagnostic accuracy of HbA1c. We therefore expanded the number of HbA1c-associated loci and tested the effect of genetic risk-scores comprised of erythrocytic or glycemic variants on incident diabetes prediction and on prevalent diabetes screening performance. Throughout this multiancestry study, we kept a focus on interancestry differences in HbA1c genetics performance that might influence race-ancestry differences in health outcomes., Methods & Findings: Using genome-wide association meta-analyses in up to 159,940 individuals from 82 cohorts of European, African, East Asian, and South Asian ancestry, we identified 60 common genetic variants associated with HbA1c. We classified variants as implicated in glycemic, erythrocytic, or unclassified biology and tested whether additive genetic scores of erythrocytic variants (GS-E) or glycemic variants (GS-G) were associated with higher T2D incidence in multiethnic longitudinal cohorts (N = 33,241). Nineteen glycemic and 22 erythrocytic variants were associated with HbA1c at genome-wide significance. GS-G was associated with higher T2D risk (incidence OR = 1.05, 95% CI 1.04-1.06, per HbA1c-raising allele, p = 3 × 10-29); whereas GS-E was not (OR = 1.00, 95% CI 0.99-1.01, p = 0.60). In Europeans and Asians, erythrocytic variants in aggregate had only modest effects on the diagnostic accuracy of HbA1c. Yet, in African Americans, the X-linked G6PD G202A variant (T-allele frequency 11%) was associated with an absolute decrease in HbA1c of 0.81%-units (95% CI 0.66-0.96) per allele in hemizygous men, and 0.68%-units (95% CI 0.38-0.97) in homozygous women. The G6PD variant may cause approximately 2% (N = 0.65 million, 95% CI 0.55-0.74) of African American adults with T2D to remain undiagnosed when screened with HbA1c. Limitations include the smaller sample sizes for non-European ancestries and the inability to classify approximately one-third of the variants. Further studies in large multiethnic cohorts with HbA1c, glycemic, and erythrocytic traits are required to better determine the biological action of the unclassified variants., Conclusions: As G6PD deficiency can be clinically silent until illness strikes, we recommend investigation of the possible benefits of screening for the G6PD genotype along with using HbA1c to diagnose T2D in populations of African ancestry or groups where G6PD deficiency is common. Screening with direct glucose measurements, or genetically-informed HbA1c diagnostic thresholds in people with G6PD deficiency, may be required to avoid missed or delayed diagnoses.

Height is a highly heritable, classic polygenic trait with approximately 700 common associated variants identified through genome-wide association studies so far. Here, we report 83 height-associated coding variants with lower minor-allele frequencies (in the range of 0.1-4.8%) and effects of up to 2 centimetres per allele (such as those in IHH, STC2, AR and CRISPLD2), greater than ten times the average effect of common variants. In functional follow-up studies, rare height-increasing alleles of STC2 (giving an increase of 1-2 centimetres per allele) compromised proteolytic inhibition of PAPP-A and increased cleavage of IGFBP-4 in vitro, resulting in higher bioavailability of insulin-like growth factors. These 83 height-associated variants overlap genes that are mutated in monogenic growth disorders and highlight new biological candidates (such as ADAMTS3, IL11RA and NOX4) and pathways (such as proteoglycan and glycosaminoglycan synthesis) involved in growth. Our results demonstrate that sufficiently large sample sizes can uncover rare and low-frequency variants of moderate-to-large effect associated with polygenic human phenotypes, and that these variants implicate relevant genes and pathways.

High blood pressure is a major risk factor for cardiovascular disease and premature death. However, there is limited knowledge on specific causal genes and pathways. To better understand the genetics of blood pressure, we genotyped 242,296 rare, low-frequency and common genetic variants in up to 192,763 individuals and used ∼155,063 samples for independent replication. We identified 30 new blood pressure- or hypertension-associated genetic regions in the general population, including 3 rare missense variants in RBM47, COL21A1 and RRAS with larger effects (>1.5 mm Hg/allele) than common variants. Multiple rare nonsense and missense variant associations were found in A2ML1, and a low-frequency nonsense variant in ENPEP was identified. Our data extend the spectrum of allelic variation underlying blood pressure traits and hypertension, provide new insights into the pathophysiology of hypertension and indicate new targets for clinical intervention., Competing Interests: Conflict of interests N. P. has received financial support from several pharmaceutical companies that manufacture either blood pressure lowering or lipid lowering agents, or both, and consultancy fees. S. K. has received Research Grant-Merck, Bayer, Aegerion; SAB-Catabasis, Regeneron Genetics Center, Merck, Celera; Equity-San Therapeutics, Catabasis; Consulting-Novartis, Aegerion, Bristol Myers-Squibb, Sanofi, AstraZeneca, Alnylam. P. Sever has received research awards from Pfizer Inc. A. Malarstig and M. Uria-Nickelsen are full time employees of Pfizer. D. Reily and M. Hoek are full time employees of Merck and co Inc. M.J. Caulfield is Chief Scientist for Genomics England a UK Government company. The authors declare no competing financial interest.

The genetic architecture of common traits, including the number, frequency, and effect sizes of inherited variants that contribute to individual risk, has been long debated. Genome-wide association studies have identified scores of common variants associated with type 2 diabetes, but in aggregate, these explain only a fraction of the heritability of this disease. Here, to test the hypothesis that lower-frequency variants explain much of the remainder, the GoT2D and T2D-GENES consortia performed whole-genome sequencing in 2,657 European individuals with and without diabetes, and exome sequencing in 12,940 individuals from five ancestry groups. To increase statistical power, we expanded the sample size via genotyping and imputation in a further 111,548 subjects. Variants associated with type 2 diabetes after sequencing were overwhelmingly common and most fell within regions previously identified by genome-wide association studies. Comprehensive enumeration of sequence variation is necessary to identify functional alleles that provide important clues to disease pathophysiology, but large-scale sequencing does not support the idea that lower-frequency variants have a major role in predisposition to type 2 diabetes.

There has been limited success in identifying causal variants underlying association signals observed in genome-wide association studies (GWAS). The use of 1000 Genomes Project (1KGP) allows the imputation to estimate the genetic information at untyped variants. However, long stretches of high linkage disequilibrium within the genome prevent us from differentiating between causal variants and perfect surrogates, thus limiting our ability to identify causal variants. Transethnic strategies have been proposed as a possible solution to mitigate this. However, these studies generally rely on imputing genotypes from multiple ancestries from 1KGP but not against population-specific reference panels. Here, we perform the first transethnic fine-mapping study across three Asian cohorts from diverse ancestries at the loci implicated with eye and blood lipid traits, using population-specific reference panels that have been generated by whole-genome sequencing samples from the same ancestry groups. Our study outlines several challenges faced in a fine-mapping exercise where one simply aims to meta-analyse existing GWAS that have been imputed against reference haplotypes from the 1KGP.

Insulin secretion has a crucial role in glucose homeostasis, and failure to secrete sufficient insulin is a hallmark of type 2 diabetes. Genome-wide association studies (GWAS) have identified loci contributing to insulin processing and secretion; however, a substantial fraction of the genetic contribution remains undefined. To examine low-frequency (minor allele frequency (MAF) 0.5-5%) and rare (MAF < 0.5%) nonsynonymous variants, we analyzed exome array data in 8,229 nondiabetic Finnish males using the Illumina HumanExome Beadchip. We identified low-frequency coding variants associated with fasting proinsulin concentrations at the SGSM2 and MADD GWAS loci and three new genes with low-frequency variants associated with fasting proinsulin or insulinogenic index: TBC1D30, KANK1 and PAM. We also show that the interpretation of single-variant and gene-based tests needs to consider the effects of noncoding SNPs both nearby and megabases away. This study demonstrates that exome array genotyping is a valuable approach to identify low-frequency variants that contribute to complex traits.

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.

Copy number variants (CNVs) extend our understanding of the genetic diversity in humans. However, the distribution and characteristics of CNVs in Asian populations remain largely unexplored, especially for rare CNVs that have emerged as important genetic factors for complex traits. In the present study, we performed an in-depth investigation of common and rare CNVs across 8,148 individuals from the three major Asian ethnic groups: Chinese (n = 1,945), Malays (n = 2,399), and Indians (n = 2,217) in Singapore, making this investigation the most comprehensive genome-wide survey of CNVs outside the European-ancestry populations to date. We detected about 16 CNVs per individual and the ratio of loss to gain events is ∼2:1. The majority of the CNVs are of low frequency (<10%), and 40% are rare (<1%). In each population, ∼20% of the CNVs are not previously catalogued in the Database of Genomic Variants (DGV). Contrary to findings from European studies, the common CNVs (>5%) in our populations are not well tagged by SNPs in Illumina 1M and 610K arrays, and most disease-associated common CNVs previously reported in Caucasians are rare in our populations. We also report noticeable population differentiation in the CNV landscape of these Asian populations, with the greatest diversity seen between the Indians and the Chinese., (© 2011 Wiley Periodicals, Inc.)

Corneal curvature (CC) is a key determinant of major eye diseases, such as keratoconus, myopia and corneal astigmatism. No prior studies have discovered the genes for CC. Here we report the findings from four genome-wide association studies of CC in 10 008 samples from three population groups in Singapore. Our discovery phase surveyed 2867 Chinese and 3072 Malays, allowing us to identify two loci that were associated with CC variation: FRAP1 on chromosome 1p36.2 and PDGFRA on chromosome 4q12. These findings were subsequently replicated in a validation study involving an additional 2953 Asian Indians and a further collection of 1116 Chinese children. The effect sizes of the identified variants were consistent across all four cohorts, with seven single nucleotide polymorphisms (SNPs) in FRAP1 (lead SNP: rs17036350, meta P-value = 4.06 × 10(-13)) and six SNPs in PDGFRA (lead SNP: rs2114039, meta P-value = 1.33 × 10(-9)) attaining genome-wide significance in the SNP-based meta-analysis of the four studies. This is the first genome-wide survey of CC variation and we have identified two implicated loci in three genetically diverse Asian populations, suggesting the presence of common genetic etiology across multiple populations.

Genome-wide association studies have discovered and confirmed a large number of loci that are implicated with disease susceptibility and severity. Polymorphisms that emerged from these studies are mostly indirectly associated to the phenotype, and the natural progression is to identify the causal variants that are functionally responsible for these association signals. Long stretches of high linkage disequilibrium (LD) benefitted the initial discovery phase in a genome-wide scan, allowing commercial genotyping products with imperfect coverage to detect genomic regions genuinely associated with the phenotype. However, regions of high LD confound the fine-mapping phase, as markers that are perfectly correlated to the causal variants display similar evidence of phenotypic association, hampering the process of differentiating the functional polymorphisms from neighboring surrogates. Here, we explore the potential of integrating information across different populations for narrowing the candidate region that a causal variant resides in, and compare the efficacy of this process of trans-population fine-mapping with the extent of variation in patterns of LD between the populations. In addition, we explore two different strategies for pooling data across multiple populations for the purpose of prioritizing the rankings of the causal variants. Our results clearly establish the benefits of trans-population analysis in reducing the number of possible candidates for the causal variants, particularly in genomic regions displaying strong evidence of inter-population LD variation. Directly integrating the statistical evidence by summing the test statistics outperforms the standard meta-analytic procedure. These findings have direct relevance to the design and analysis of ongoing fine-mapping studies., (© 2010 Wiley-Liss, Inc.)

Research on the role of copy number variations (CNVs) in the genetic risk of diseases in Asian populations has been hampered by a relative lack of reference CNV maps for Asian populations outside the East Asians. In this article, we report the population characteristics of CNVs in Chinese, Malay, and Asian Indian populations in Singapore. Using the Illumina Human 1M Beadchip array, we identify 1,174 CNV loci in these populations that corroborated with findings when the same samples were typed on the Affymetrix 6.0 platform. We identify 441 novel loci not previously reported in the Database of Genomic Variations (DGV). We observe a considerable number of loci that span all three populations and were previously unreported, as well as population-specific loci that are quite common in the respective populations. From this we observe the distribution of CNVs in the Asian Indian population to be considerably different from the Chinese and Malay populations. About half of the deletion loci and three-quarters of duplication loci overlap UCSC genes. Tens of loci show population differentiation and overlap with genes previously known to be associated with genetic risk of diseases. One of these loci is the CYP2A6 deletion, previously linked to reduced susceptibility to lung cancer., ((c) 2010 Wiley-Liss, Inc.)

The Singapore Genome Variation Project (SGVP) provides a publicly available resource of 1.6 million single nucleotide polymorphisms (SNPs) genotyped in 268 individuals from the Chinese, Malay, and Indian population groups in Southeast Asia. This online database catalogs information and summaries on genotype and phased haplotype data, including allele frequencies, assessment of linkage disequilibrium (LD), and recombination rates in a format similar to the International HapMap Project. Here, we introduce this resource and describe the analysis of human genomic variation upon agglomerating data from the HapMap and the Human Genome Diversity Project, providing useful insights into the population structure of the three major population groups in Asia. In addition, this resource also surveyed across the genome for variation in regional patterns of LD between the HapMap and SGVP populations, and for signatures of positive natural selection using two well-established metrics: iHS and XP-EHH. The raw and processed genetic data, together with all population genetic summaries, are publicly available for download and browsing through a web browser modeled with the Generic Genome Browser.