#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Use of >100,000 NHLBI Trans-Omics for Precision Medicine (TOPMed) Consortium whole genome sequences improves imputation quality and detection of rare variant associations in admixed African and Hispanic/Latino populations


Authors: Madeline H. Kowalski aff001;  Huijun Qian aff002;  Ziyi Hou aff003;  Jonathan D. Rosen aff001;  Amanda L. Tapia aff001;  Yue Shan aff001;  Deepti Jain aff004;  Maria Argos aff005;  Donna K. Arnett aff006;  Christy Avery aff007;  Kathleen C. Barnes aff008;  Lewis C. Becker aff009;  Stephanie A. Bien aff010;  Joshua C. Bis aff011;  John Blangero aff012;  Eric Boerwinkle aff013;  Donald W. Bowden aff015;  Steve Buyske aff016;  Jianwen Cai aff017;  Michael H. Cho aff018;  Seung Hoan Choi aff020;  Hélène Choquet aff021;  L. Adrienne Cupples aff022;  Mary Cushman aff024;  Michelle Daya aff008;  Paul S. de Vries aff014;  Patrick T. Ellinor aff020;  Nauder Faraday aff009;  Myriam Fornage aff026;  Stacey Gabriel aff027;  Santhi K. Ganesh aff028;  Misa Graff aff007;  Namrata Gupta aff027;  Jiang He aff030;  Susan R. Heckbert aff031;  Bertha Hidalgo aff033;  Chani J. Hodonsky aff007;  Marguerite R. Irvin aff033;  Andrew D. Johnson aff023;  Eric Jorgenson aff021;  Robert Kaplan aff035;  Sharon L. R. Kardia aff036;  Tanika N. Kelly aff030;  Charles Kooperberg aff010;  Jessica A. Lasky-Su aff018;  Ruth J. F. Loos aff037;  Steven A. Lubitz aff020;  Rasika A. Mathias aff009;  Caitlin P. McHugh aff004;  Courtney Montgomery aff039;  Jee-Young Moon aff035;  Alanna C. Morrison aff014;  Nicholette D. Palmer aff015;  Nathan Pankratz aff040;  George J. Papanicolaou aff041;  Juan M. Peralta aff012;  Patricia A. Peyser aff036;  Stephen S. Rich aff042;  Jerome I. Rotter aff043;  Edwin K. Silverman aff018;  Jennifer A. Smith aff044;  Nicholas L. Smith aff031;  Kent D. Taylor aff043;  Timothy A. Thornton aff004;  Hemant K. Tiwari aff046;  Russell P. Tracy aff047;  Tao Wang aff048;  Scott T. Weiss aff018;  Lu-Chen Weng aff020;  Kerri L. Wiggins aff011;  James G. Wilson aff049;  Lisa R. Yanek aff009;  Sebastian Zöllner aff050;  Kari E. North aff007;  Paul L. Auer aff053;  ;  ;  Laura M. Raffield aff054;  Alexander P. Reiner aff031;  Yun Li aff001
Authors place of work: Department of Biostatistics, University of North Carolina, Chapel Hill, North Carolina, United States of America aff001;  Department of Statistics and Operation Research, University of North Carolina, Chapel Hill, North Carolina, United States of America aff002;  Department of Biomedical Informatics, Harvard Medical School, Boston, Massachusetts, United States of America aff003;  Department of Biostatistics, University of Washington, Seattle, Washington, United States of America aff004;  Division of Epidemiology and Biostatistics, University of Illinois at Chicago, Chicago, Illinois, United States of America aff005;  College of Public Health, University of Kentucky, Lexington, Kentucky, United States of America aff006;  Department of Epidemiology, University of North Carolina, Chapel Hill, North Carolina, United States of America aff007;  Department of Medicine, Anschutz Medical Campus, University of Colorado Denver, Aurora, Colorado, United States of America aff008;  GeneSTAR Research Program, Department of Medicine, Johns Hopkins School of Medicine, Baltimore, Maryland, United States of America aff009;  Division of Public Health Sciences, Fred Hutchinson Cancer Research Center, Seattle, Washington, United States of America aff010;  Cardiovascular Health Research Unit, Department of Medicine, University of Washington, Seattle, Washington, United States of America aff011;  Department of Human Genetics and South Texas Diabetes Institute, University of Texas Rio Grande Valley School of Medicine, Brownsville, Texas, United States of America aff012;  Human Genome Sequencing Center, University of Texas Health Science Center at Houston; Baylor College of Medicine, Houston, Texas, United States of America aff013;  Human Genetics Center, Department of Epidemiology, Human Genetics, and Environmental Sciences, School of Public Health, The University of Texas Health Science Center at Houston, Houston, Texas, United States of America aff014;  Department of Biochemistry, Wake Forest School of Medicine, Winston-Salem, North Carolina, United States of America aff015;  Department of Statistics, Rutgers University, Piscataway, New Jersey, United States of America aff016;  Collaborative Studies Coordinating Center, Department of Biostatistics, University of North Carolina, Chapel Hill, North Carolina, United States of America aff017;  Channing Division of Network Medicine, Brigham and Women’s Hospital, Boston, Massachusetts, United States of America aff018;  Department of Medicine, Harvard Medical School, Boston, Massachusetts, United States of America aff019;  Program in Medical and Population Genetics, The Broad Institute of MIT and Harvard, Cambridge, Massachusetts, United States of America aff020;  Division of Research, Kaiser Permanente Northern California, Oakland, California, United States of America aff021;  Department of Biostatistics, Boston University School of Public Health, Boston, Massachusetts, United States of America aff022;  Framingham Heart Study, Framingham, Massachusetts, United States of America aff023;  Departments of Medicine & Pathology, Larner College of Medicine, University of Vermont, Colchester, Vermont, United States of America aff024;  Cardiac Arrhythmia Service and Cardiovascular Research Center, Massachusetts General Hospital, Boston, Massachusetts, United States of America aff025;  School of Public Health, The University of Texas Health Science Center, Houston, Texas, United States of America aff026;  Genomics Platform, Broad Institute, Cambridge, Massachusetts, United States of America aff027;  Department of Internal Medicine, University of Michigan, Ann Arbor, Michigan, United States of America aff028;  Department of Human Genetics, University of Michigan, Ann Arbor, Michigan, United States of America aff029;  Department of Epidemiology, Tulane University School of Public Health and Tropical Medicine, New Orleans, Los Angeles, United States of America aff030;  Department of Epidemiology, University of Washington, Seattle, Washington, United states of America aff031;  Kaiser Permanente Washington Health Research Institute, Kaiser Permanente Washington, Seattle, Washington, United States of America aff032;  Department of Epidemiology, Ryals School of Public Health, University of Alabama at Birmingham, Birmingham, Alabama, United States of America aff033;  Population Sciences Branch, Division of Intramural Research, National Heart, Lung and Blood Institute, Framingham, Massachusetts, United States of America aff034;  Department of Epidemiology & Population Health, Albert Einstein College of Medicine, Bronx, New York, United States of America aff035;  Department of Epidemiology, School of Public Health, University of Michigan, Ann Arbor, Michigan, United States of America aff036;  The Charles Bronfman Institute for Personalized Medicine, Icahn School of Medicine at Mount Sinai, New York, New York, United States of America aff037;  The Mindich Child Health and Development Institute, Icahn School of Medicine at Mount Sinai, New York, New York, United States of America aff038;  Department of Genes and Human Disease, Oklahoma Medical Research Foundation, Oklahoma City, Oklahoma, United States of America aff039;  Department of Laboratory Medicine and Pathology, University of Minnesota, Minneapolis, Minnesota, United States of America aff040;  National Heart, Lung, and Blood Institute, Division of Cardiovascular Sciences, PPSP/EB, NIH, Bethesda, Maryland, United States of America aff041;  Center for Public Health Genomics, Department of Public Health Sciences, University of Virginia, Charlottesville, Virginia, United States of America aff042;  The Institute for Translational Genomics and Population Sciences, Department of Pediatrics, Los Angeles Biomedical Research Institute at Harbor-UCLA Medical Center, Torrance, California, United States of America aff043;  Department of Epidemiology, School of Public Health, University of Michigan, Ann Arbor, Michigan, United States of America aff044;  Seattle Epidemiologic Research and Information Center, Department of Veterans Affairs Office of Research and Development, Seattle, Washington, United States of America aff045;  Department of Biostatistics, Ryals School of Public Health, University of Alabama at Birmingham, Birmingham, Alabama, United States of America aff046;  Departments of Pathology & Laboratory Medicine and Biochemistry, Larrner College of Medicine, University of Vermont, Colchester, Vermont, United States of America aff047;  Department of Epidemiology and Population Health, Albert Einstein College of Medicine, Bronx, New York, United States of America aff048;  Department of Physiology and Biophysics, University of Mississippi Medical Center, Jackson, Mississippi, United States of America aff049;  Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, Michigan, United States of America aff050;  Department of Psychiatry, University of Michigan, Ann Arbor, Michigan, United States of America aff051;  Carolina Center of Genome Sciences, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America aff052;  Zilber School of Public Health, University of Wisconsin-Milwaukee, Milwaukee, Wisconsin, United States of America aff053;  Department of Genetics, University of North Carolina, Chapel Hill, North Carolina, United States of America aff054;  Department of Computer Science, University of North Carolina, Chapel Hill, North Carolina, United States of America aff055
Published in the journal: Use of >100,000 NHLBI Trans-Omics for Precision Medicine (TOPMed) Consortium whole genome sequences improves imputation quality and detection of rare variant associations in admixed African and Hispanic/Latino populations. PLoS Genet 15(12): e32767. doi:10.1371/journal.pgen.1008500
Category: Research Article
doi: https://doi.org/10.1371/journal.pgen.1008500

Summary

Most genome-wide association and fine-mapping studies to date have been conducted in individuals of European descent, and genetic studies of populations of Hispanic/Latino and African ancestry are limited. In addition, these populations have more complex linkage disequilibrium structure. In order to better define the genetic architecture of these understudied populations, we leveraged >100,000 phased sequences available from deep-coverage whole genome sequencing through the multi-ethnic NHLBI Trans-Omics for Precision Medicine (TOPMed) program to impute genotypes into admixed African and Hispanic/Latino samples with genome-wide genotyping array data. We demonstrated that using TOPMed sequencing data as the imputation reference panel improves genotype imputation quality in these populations, which subsequently enhanced gene-mapping power for complex traits. For rare variants with minor allele frequency (MAF) < 0.5%, we observed a 2.3- to 6.1-fold increase in the number of well-imputed variants, with 11–34% improvement in average imputation quality, compared to the state-of-the-art 1000 Genomes Project Phase 3 and Haplotype Reference Consortium reference panels. Impressively, even for extremely rare variants with minor allele count <10 (including singletons) in the imputation target samples, average information content rescued was >86%. Subsequent association analyses of TOPMed reference panel-imputed genotype data with hematological traits (hemoglobin (HGB), hematocrit (HCT), and white blood cell count (WBC)) in ~21,600 African-ancestry and ~21,700 Hispanic/Latino individuals identified associations with two rare variants in the HBB gene (rs33930165 with higher WBC [p = 8.8x10-15] in African populations, rs11549407 with lower HGB [p = 1.5x10-12] and HCT [p = 8.8x10-10] in Hispanics/Latinos). By comparison, neither variant would have been genome-wide significant if either 1000 Genomes Project Phase 3 or Haplotype Reference Consortium reference panels had been used for imputation. Our findings highlight the utility of the TOPMed imputation reference panel for identification of novel rare variant associations not previously detected in similarly sized genome-wide studies of under-represented African and Hispanic/Latino populations.

Keywords:

Haplotypes – Alleles – Genome-wide association studies – Hemoglobin – Hematology – Variant genotypes – Genotyping – Consortia

Introduction

Genotype imputation, despite being a standard practice in modern genetic association studies, remains challenging in populations of Hispanic/Latino or African ancestry, particularly for rare variants [16]. One obstacle lies in the lack of appropriate whole genome sequence reference panels for these admixed populations. For individuals of European descent, the relevant haplotypes available have increased by more than 500 times from 120 phased sequences in HapMap2 [7] to more than 64,000 phased sequences in Haplotype Reference Consortium (HRC) [8] reference. However, HRC is predominantly European (other than included 1000 Genomes Project Phase 3 (1000G) SNPs) and includes mostly low-coverage sequencing data (4-8x coverage). The state-of-the-art reference panels for African-ancestry (AA) and Hispanic/Latino cohorts, including the 1000 Genomes Project Phase 3 (1000G) [9] and the Consortium on Asthma among African-ancestry Populations in the Americas (CAAPA) [10], are at least one order of magnitude smaller than HRC. This is especially problematic given the complex LD structure in admixed populations. The NHLBI Trans-Omics for Precision Medicine (TOPMed) Project has recently generated deep-coverage (mean depth 30x) whole genome sequencing (WGS) on more than 50,000 individuals from >26 cohorts and from diverse ancestral backgrounds (notably including ~26% AA and ~10% Hispanic/Latino participants), and now provides an unprecedented opportunity for substantially enhancing imputation quality in under-represented admixed populations and subsequently boosting power for mapping genes and regions underlying complex traits. Here we demonstrate the improvements in rare variant imputation quality in AA and Hispanic/Latino populations using TOPMed as a reference panel versus 1000G and HRC panels, and subsequently identify two low-frequency/rare HBB variant associations with blood cell traits in AA and Hispanic/Latino samples using TOPMed-imputed genotyping array data.

Results and discussion

The cohort and ancestry composition of the TOPMed freeze 5b whole genome sequence reference panel used in our study and the samples with array-based genotyping used for imputation and hematological traits association analyses in self-identified AA and Hispanic/Latino individuals are summarized in S1 and S2 Tables, respectively. We first selected two large U.S. minority cohorts—one AA and one Hispanic/Latino—in order to comprehensively evaluate imputation quality: the Jackson Heart Study (JHS, all AA, n = 3,082) and the Hispanic Community Health Study/Study of Latinos (HCHS/SOL, all Hispanic/Latino, n = 11,887). Both the JHS and HCHS/SOL have external sources of dense genotype data available for comparison. JHS is the largest AA general population cohort sequenced in TOPMed freeze 5b. Therefore, we removed JHS samples from the TOPMed freeze 5b reference panel prior to performing imputation into JHS samples using SNPs genotyped on the Affymetrix 6.0 array, treating the TOPMed freeze 5b calls as true genotypes for evaluation of imputation quality in JHS. HCHS/SOL is the largest and most regionally diverse population-based cohort of Hispanic/Latino individuals living in the US. For HCHS/SOL, we used the entire set of 100,506 phased sequences from TOPMed freeze 5b (including JHS) as reference and performed imputation into 11,887 Hispanic/Latino samples genotyped on the Illumina Omni 2.5 SOL custom array (with high quality genotypes at 2,293,536 markers). As the external source of genotype validation in HCHS/SOL, we used genotypes from the Illumina MEGA array genotyping data (containing >1.7 million multi-ethnic global markers, including low frequency coding variants and ancestry-specific variants) available in the same HCHS/SOL samples to assess imputation quality, evaluating 688,189 imputed markers available on MEGA but not on Omni2.5.

Compared with the 1000G Phase 3 reference panel [9], we were able to increase the number of well-imputed variants from ~28 and ~35 million to ~51 and ~58 million in JHS and HCHS/SOL, respectively (see S3 Table for genome-wide distribution of well-imputed variants). We defined well-imputed variants based on our previous work [1, 2, 4], using MAF-specific estimated R2 thresholds to ensure an average R2 ≥ 0.8 in each imputed cohort separately. For all rare variants with MAF < 0.5%, we observed ~4.2X (2.3X) and ~6.1X (3.3X) increases in the number of well-imputed variants in JHS (HCHS/SOL), compared with 1000G and HRC, respectively. We also observed 22% (11%) and 34% (20%) increases in imputation information content (as measured by average true R2, which is the squared Pearson correlation between imputed and true genotypes) (Fig 1 and S1 Fig, Table 1). For very rare variants with MAF <0.05%, we observed ~22.1X (5.8X) and ~11.8X (10.7X) increases in the number of well-imputed variants, with 6% (5%) and 13% (11%) increases in average true R2, in JHS (HCHS/SOL), compared with 1000G and HRC respectively. Mismatch rates between true and imputed genotypes were low; using the program CalcMatch, the mean concordance for heterozygote individuals (generally the hardest to impute) for Jackson Heart Study is 97.5% for all well imputed variants in Table 1, 96.6% for MAF <0.5%, and 97.6% for MAF < 0.05%. For HCHS/SOL, the mean concordance is 98.2% for all well imputed variants, 92.9% for MAF <0.5%, and 83.8% for MAF < 0.05%. Most well-imputed variants from 1000G and HRC were also included in TOPMed freeze 5b imputation results (S2 Fig).

Comparison of imputation reference panels, for variants with MAF < 1%.
Fig. 1. Comparison of imputation reference panels, for variants with MAF < 1%.
Imputation quality (measured by true R2 [Y-axis]) is plotted with progressively more stringent post-imputation filtering from left to right, with filtering according to estimated R2 (X-axis), for variants with MAF < 1%. Top panels are for the JHS cohort and bottom panels for the HCHS/SOL cohort. Three reference panels are shown: TOPMed (TOPMed freeze 5b), 1000G (the 1000 Genomes Phase 3), and HRC (the Haplotype Reference Consortium).
Tab. 1. Number of well-imputed variants using TOPMed freeze 5b, 1000 Genomes Phase 3 (1000G) and Haplotype Reference Consortium (HRC).
Number of well-imputed variants using TOPMed freeze 5b, 1000 Genomes Phase 3 (1000G) and Haplotype Reference Consortium (HRC).

Even for extremely rare variants with sample minor allele count (MAC) <10 (including cohort singleton variants in the target JHS cohort), average information content rescued (again measured by true R2) was >86%. For example, out of the 8.67 million singleton variants discovered in JHS by TOPMed WGS, 72% (6.24 million) can be well-imputed using Affymetrix 6.0 genotypes and using TOPMed freeze 5b (without JHS individuals) as reference, with an average true R2 of 0.92 (Table 2). Singletons within JHS are defined as variants with MAC = 1 among the JHS samples but which are present in multiple copies in the reference panel. Specifically, the average reference MAC is 29.3 before post-imputation quality control (QC) and 31.0 after QC, with all variants having a MAC>5 in the overall reference panel. Imputation quality is similarly high when examining extremely rare MAC variants in the reference panel, and even higher, as expected, with higher MAC variants within the JHS sample (S4 and S5 Tables). Similar observations hold true for HCHS/SOL, with slightly lower imputation quality (S6 and S7 Tables). Compared to JHS African Americans, the lower imputation quality in HCHS/SOL Hispanic/Latino individuals is likely attributable to multiple reasons, including (1) the more complex LD structure among Hispanic/Latino individuals due to the admixture of three ancestral populations; (2) the availability of a much smaller subset of rare variants for quality evaluation through MEGA array genotyping in HCHS/SOL (in contrast to the availability of nearly all segregating variants in JHS through high-coverage sequencing); and (3) the smaller number of relevant haplotypes in the TOPMed freeze 5b reference (~26% self-identified AAs compared to ~10% self-identified Hispanics/Latinos). Imputation quality for rare and low-frequency variants that are estimated to be well imputed in Table 1 is further stratified by regional background in HCHS/SOL and displayed in S8 Table. We note that greater numbers of AA and Hispanic/Latino individuals will be included in future releases of sequencing datasets from TOPMed, which we anticipate will further improve imputation quality; inclusion of JHS itself in imputation for other AA cohorts would also improve imputation quality.

Tab. 2. Imputation quality for rare variants (minor allele count< = 10) in the Jackson Heart Study (JHS).
Imputation quality for rare variants (minor allele count< = 10) in the Jackson Heart Study (JHS).

Encouraged by these substantial gains in information content for low-frequency and rare variants, we proceeded with imputation in several additional AA and Hispanic/Latino data sets with array-based genotyping (S1 and S9 Tables), followed by association analyses with quantitative blood cell traits to evaluate the power of TOPMed freeze 5b-based imputation in minorities for discovery of genetic variants underlying complex human traits. We specifically chose hematological traits for several reasons. First, these traits are important intermediate clinical phenotypes for a variety of cardiovascular, hematologic, oncologic, immunologic, and infectious diseases [11]. Second, these traits have family-based heritability estimates in the range of 40–65% [12, 13], and have been highly fruitful for gene-mapping with >2,700 common and rare variants identified, though primarily in individuals of European ancestry [1419]. Third, these traits remain under-studied in admixed AA and Hispanic/Latino populations, despite evidence for the existence of variants with distinct genetic architecture in AAs and Hispanics/Latinos [2022]. For example, while hundreds of variants identified in genome-wide association studies (GWAS) of WBC in individuals of European descent explain only ~7% of array heritability, the African specific Duffy null variant DARC rs2814778 alone accounts for 15–20% of population-level WBC variability in AAs [23]. Finally, we have previously successfully leveraged deep-coverage exome sequencing-based imputation using resources from the Exome Sequencing Project for more powerful mapping of genes and regions associated with hematological traits in AAs [1]. Hemoglobin level (HGB), hematocrit (HCT), and WBC were chosen for our primary phenotypic analysis because these traits are available in the largest sample size among the AA and Hispanics/Latinos included in our discovery cohorts.

Our imputation sample used for discovery blood cell trait association analyses included eight cohorts (21,513 AAs and 21,689 Hispanics/Latinos) (S1 Table). These discovery samples do not overlap with individuals sequenced as part of TOPMed freeze 5b (S2 Table). We used the full set of 100,506 phased sequences from TOPMed freeze 5b (including JHS) as the imputation reference panel. We then carried out AA- and Hispanic/Latino-stratified association analyses with quantitative HGB, HCT, and total WBC separately in each cohort genotyping array data set, accounting for ancestry and relatedness. The genome-wide association results for each imputed cohort data set were then meta-analyzed within each ancestry group. S3S8 Figs show the Manhattan plots from ethnic-specific meta-analyses for each trait. QQ plots (S9S14 Figs) show no obvious early departure, with genomic control lambda ranging from 1.008 to 1.044, indicating minimal global inflation of test statistics. For replication of any novel associations identified in the imputation-based discovery analysis, we utilized WGS genotype data and hematological trait data from the non-overlapping set of AA individuals within TOPMed freeze 5b (S10 Table) (see Methods for details).

We first evaluated association statistics for variants previously associated with HGB, HCT, or WBC count in AA and Hispanic/Latino populations (summarized in S11 Table). We assembled a list of 24 AA and 13 Hispanic/Latino previously identified autosomal signals from prior published GWAS or exome-based studies [1, 19, 20, 2430]. Our lists excluded variants reported in multi-ethnic cohorts or meta-analysis including individuals of non-AA or non- Hispanic/Latino ancestry to guard against the scenario that the reported signals were driven predominantly by individuals of European or Asian ancestry. Among the previously reported 24 AA and 13 Hispanic/Latino variants, all but five (four SNPs and a 3.8 kb deletion variant esv2676630) passed variant quality-control filters in TOPMed freeze 5b and were subsequently well-imputed in our target AA and Hispanic/Latino data sets with a stringent post-imputation R2 filter of >0.8 (detailed in S12 Table). Among the 31 known HGB, HCT, or WBC count associations testable with TOPMed freeze 5b, our imputed/discovery cohorts confirmed 84% of these previously reported findings with a consistent direction of effect, using a stringent genome-wide significant threshold of p<5x10-8. Using more lenient p-value thresholds, we could replicate 94% (p<5x10-6) and 100% (p<0.05) of the previously reported findings with the same direction of effect. While these results help confirm the overall validity of our hematological trait association results, it is important to note for these comparisons that many of the samples included in the current TOPMed freeze 5b imputed genome-wide association analysis were also used in the publications originally reporting associations in AA and Hispanic/Latino individuals.

Our ancestry-stratified imputation-based discovery meta-analysis revealed two blood cell trait associations that have not been previously reported, at a genome-wide significant threshold of p<5×10−9 in Hispanics/Latinos and p<1x10-9 in AA populations, based on appropriate significance thresholds for whole genome sequencing analysis [31]. One signal was revealed in each ancestry group: hemoglobin subunit beta (HBB) missense (p.Glu7Lys) variant rs33930165 (gb38:11:5227003:C:T) associated with increased WBC in AAs (β = 0.35 and p = 8.8x10-15, adjusting for SNP rs2814778 and removing potential minor allele homozygotes) (Table 3), and HBB stop-gain (p.Gln40Ter) variant rs11549407 (gb38:11:5226774:G:A) associated with lower HGB and HCT in Hispanics/Latinos (β = -1.92, p = 1.5x10-12; β = -1.66, p = 8.8x10-10). Both variants were either low frequency or rare: the HBB missense variant rs33930165 (hemoglobin C variant) has a MAF of 1.14% among the imputed AA discovery samples and is even rarer in non-AA individuals (absent in Europeans in 1000G); the stop gain variant rs11549407 has a MAF of 0.03% (MAC ~ 15) among the imputed Hispanics/Latinos and is monomorphic among the AAs. Both variants are classified as pathogenic in ClinVar. Both variants were well imputed with R2 ranging from 0.831 to 0.994 and 0.862 to 0.999 in the contributing AA and Hispanic/Latino cohorts, respectively (Table 3). Due to the low allele frequency of these variants in AAs and Hispanics/Latinos and even lower frequency in individuals of European descent, both variants were imputed with lower quality using other reference panels (S13 and S14 Tables): the missense variant HBB rs33930165 had R2 as low as 0.127 and 0.456 using 1000G and HRC, respectively, as references; the HBB stop-gain variant rs11549407 was not available in the 1000G reference panel and had R2 as low as 0.413 using HRC as the reference panel. Carrying the 1000G and HRC imputed genotypes forward to association analyses with hematological traits in the subset of our target imputation cohorts where the variants were well imputed (R2 > 0.8), we observed none of the p-values exceeded genome-wide significance threshold. This explains why these variants were not detectable at a genome-wide significant level using previously available imputation reference panels, with obvious implications for other complex trait association studies in ancestrally diverse study populations.

Tab. 3. Novel variants detected in TOPMed freeze 5b imputed Hispanic/Latino and African ancestry cohorts, in association analyses with white blood cell count, hemoglobin, and hematocrit.
Novel variants detected in TOPMed freeze 5b imputed Hispanic/Latino and African ancestry cohorts, in association analyses with white blood cell count, hemoglobin, and hematocrit.

Both of our previously unreported genotype-trait associations involve coding variants of HBB, which encodes the beta polypeptide chains in adult hemoglobin. The HBB stop gain (p.Gln40Ter) variant 11:5226774:G:A (rs11549407) is the most common cause of beta zero thalassemia in West Mediterranean countries, particularly among the founder population of Sardinia [32, 33], where the variant has a population allele frequency of ~5%. The Sardinian population is represented in the HRC reference panel (~3500 individuals), which likely contributes to the reasonable imputation quality observed using HRC in most but not all cohorts, in contrast to the absence of this variant in the 1000G reference panel due to very low minor allele count, though imputation quality was clearly improved with the TOPMed freeze 5b reference panel. The p.Gln40Ter mutation is much less prevalent outside of the Western Mediterranean, but has been detected among individuals with beta thalassemia among admixed populations from Central and South America [34, 35], which are geographically and genetically similar to some of the Hispanic/Latino samples included in our imputation-based discovery sample. While the individuals carrying the HBB p.Gln40Ter allele in our unselected population-based Hispanic/Latino sample were all imputed heterozygotes (consistent with “thalassemia minor” and generally considered healthy), there is increasing evidence that silent carriers of beta-thalassemia and sickle cell mutations may be at risk for various health-related conditions [36, 37]. Due to the relatively small number of Hispanic/Latino individuals with blood cell trait data in TOPMed freeze 5b (n~1,080), including only one heterozygote carrier of rs11549407 in those with blood cell traits measured, we were unable to perform a well-powered replication of the association of rs11549407 with HGB and HCT. Moderate anemia is known to occur in some individuals with thalassemia minor, however, concordant with our results [38].

The association of the HBB missense (p.Glu7Lys) variant 11:5227003:C:T or rs33930165 with higher total WBC (β = 0.35, p = 8.8x10-15) among AA was unexpected; rs33930165 has been associated with red blood cell indices such as mean corpuscular hemoglobin concentration [20] but not with white blood cell traits. Because of the higher allele frequency of this variant and also the larger number of AA samples (n = 6,743) in TOPMed freeze 5b, we were able to replicate this HBB rs33930165 association with total WBC in an independent sample (β = 0.27 and p = 4.6x10-4) of AA individuals. By contrast, there was no significant association of the HBB rs33930165 p.Glu7Lys variant with HGB and a modest association with lower HCT in the AA discovery and replication data sets (discovery HCT β = -0.122, p = 0.012; HGB β = 0.110, p = 0.022; replication HCT β = -0.239, p = 0.002; HGB β = -0.009, p = 0.909). The minor allele T of rs33930165 encodes an abnormal form of hemoglobin, Hb C, which in the homozygous state is associated with mild chronic hemolytic anemia and mild to moderate splenomegaly [39]. In our discovery and replication data sets, there were no individuals homozygous for the Hb C variant, nor any compound heterozygotes for Hb S/C (Hb S is sickling form of hemoglobin and individuals homozygous for Hb S have sickle cell disease), which excludes the possibility that the apparently higher WBC is driven by an “inflammatory response” confined to a small number of individuals clinically affected by sickle cell disease or hemoglobin C disease. We next evaluated the association of HBB rs33930165 with circulating number of WBC subtypes, including neutrophils, monocytes, lymphocytes, basophils, and eosinophils. S15 Table shows the results in our AA imputation-based discovery data sets (S16 Table), and TOPMed freeze 5b WGS replication samples (S17 Table), which suggest that the apparent association of HBB rs33930165 with total WBC is mainly driven by an association with higher lymphocyte count, with perhaps a more modest association with higher neutrophil count. Further studies are needed to delineate the putative mechanism of this unexpected association.

Our findings showcase the power of the large, ancestrally diverse TOPMed WGS data set as an imputation reference panel for admixed populations, in terms of both imputation quality and accuracy (especially for rare variants) and subsequent association studies for complex traits. Specifically, we identified two rare variants associated with hematological traits in AA and Hispanic/Latino populations and were able to validate our initial HBB association with WBC in an independent replication sample of sequenced individuals. In our study, we used EAGLE and minimac4 for imputation. We anticipate that the advantages of TOPMed as a reference panel also manifest when using alternative imputation methods. However, making TOPMed available as a reference panel compatible with each imputation method (e.g., corresponding recombination rate information) would be essential. In addition, computing time and memory usage should be taken into consideration as not all existing methods can scale to ~100 million markers in populations containing over thousands of individuals. TOPMed freeze 5b imputation is slightly more computationally intensive than use of the HRC reference panel (and takes nearly eight times longer than 1000G based imputation using the Michigan imputation server). However, we feel this increase in computational time is more than justified by the large number of additional well-imputed variants. We would note that the gains in imputation quality for AA and Hispanic/Latino populations using the TOPMed WGS reference panel likely do not apply to populations poorly represented in TOPMed freeze 5b (such as South Asians); future large-scale sequencing, including in later freezes of TOPMed, will improve imputation quality further across global populations.

Future studies should also evaluate potential increases in statistical power for gene- and region- based tests using TOPMed imputed data. To demonstrate the potential gains, we have performed a targeted analysis of genes previously identified for their association with white blood cell count or hemoglobin/hematocrit levels in exome genotyping arrays or exome sequencing studies. We compared gene-based SKAT test results at these known loci using TOPMed freeze 5b based imputation to gene-based tests performed using 1000G and HRC reference panels. These results are presented in S18S21 Tables and demonstrate that in both African ancestry and Hispanic/Latino populations more previously implicated genes from exome arrays or sequencing based studies were significant using TOPMed freeze 5b as an imputation reference panel versus 1000G phase 3 or HRC imputation. Further exploration of gene- and region-based tests is warranted in future studies, however. We expect the combination of high-quality imputation and higher depth sequencing datasets in larger cohorts of individuals will provide increased power for all rare variant association analyses in diverse populations in the near future.

Methods

Ethics statement

We here performed secondary data analysis on deidentified data only (exempt research). Access to TOPMed data was approved by the University of North Carolina at Chapel Hill Institutional Review board (study 16–2213). All individual studies included in TOPMed were approved by relevant local ethical review boards.

TOPMed 5b sequencing and phasing

The reference panel used for imputation was obtained from deep-coverage whole genome sequences derived from NHLBI’s TOPMed program (www.nhlbiwgs.org), freeze 5b (September 2017). This release included 54,035 non-duplicated, dbGaP released samples, of whom 50,253 have consent to be part of an imputation reference panel. The parent studies that contributed these 50,253 samples are listed in S2 Table. Specific to our analyses, freeze 5b includes 3,082 individuals from the Jackson Heart Study, who were removed from the reference panel for our analysis of imputation quality in this particular cohort. Overall, freeze 5b included 54% European ancestry, 26% AA, 10% Hispanic/Latino, 7% Asian, and 3% other ancestry samples. Detailed sequencing methods used in TOPMed are available at https://www.nhlbiwgs.org/topmed-whole-genome-sequencing-project-freeze-5b-phases-1-and-2. In brief, WGS with mean genome coverage ≥30x was completed at six sequencing centers (New York Genome Center, the Broad Institute of MIT and Harvard, the University of Washington Northwest Genomics Center, Illumina Genomic Services, Macrogen Corp., and Baylor Human Genome Sequencing Center). Sequence data files were transferred from sequencing centers to the TOPMed Informatics Research Center (IRC), where reads were aligned to human genome build GRCh38, using a common pipeline, and joint genotype calling was undertaken. Variants were filtered using a machine learning based support vector machine (SVM) approach, using variants present on genotyping arrays as positive controls and variants with many Mendelian inconsistencies as negative controls. After filtering potentially problematic variant sites, freeze 5b contained ~438 million single nucleotide polymorphisms and ~33 million short insertion-deletion variants. For our imputation analyses, we excluded from the reference panel variants with an overall allele count of 5 or less (leaving 88,062,238 variants in our reference panel, Table 1). Additional sample level quality control (such as detection of sex mismatches, pedigree discrepancies, sample swaps, etc.) was undertaken by the TOPMed Data Coordinating Center (DCC).

Genome-wide genotyping array data sets used for evaluation of imputation quality and/or phenotype association analysis

Hispanic Community Health Study/Study of Latinos (HCHS/SOL)

The HCHS/SOL cohort began in 2006 as a prospective study of Hispanic/Latino populations in the U.S. [4042]. From 2008 to 2011, 16,415 adults were recruited from a random sample of households in four communities (the Bronx, Chicago, Miami, and San Diego). Each Field Center recruited >4,000 participants from diverse socioeconomic groups. Most participants self-identified as having Cuban, Dominican, Puerto Rican, Mexican, Central American, or South American heritage. The cohort has been genotyped both using an Illumina Omni2.5M array (plus 150,000 custom SNP, including ancestry-informative markers, Amerindian population specific variants, previously identified GWAS hits, and other candidate polymorphisms for a total of 2,293,715 SNPs) [43] and using the Illumina Multi-Ethnic Genotyping Array (MEGA) array (containing a total of 1,705,969 SNPs) in efforts from the Population Architecture for Genetic Epidemiology [44] consortium to better assess variation in non-European populations. The MEGA array also includes additional exonic, functional, and clinically-relevant variants. Illumina 2.5M array genotypes were available for 12,802 samples, among whom 11,887 samples also had MEGA array genotypes. The Illumina Omni2.5M array was used for imputation to the TOPMed reference panel, with the MEGA array treated as true genotypes for evaluation of imputation quality. For association analysis, imputation was performed on 11,887 samples after merging Omni2.5M array genotypes and MEGA array genotypes (MEGA genotypes were used for variants in both arrays, which resulted in 2,144,214 variants after quality control). Regional background (for evaluation of stratified imputation quality in S8 Table) was defined using both self-identified background and genetic markers, as described in [43]. For the hematological traits association analysis, 11,588 Hispanic/Latino participants were included.

Women’s Health Initiative

The Women’s Health Initiative (WHI) [45] is a long-term national health study focused heart disease, cancer, and osteoporotic fractures in older women. WHI originally enrolled 161,808 women aged 50–79 between 1993 and 1998 at 40 centers across the US, including both a clinical trial (including three trials for hormone therapy, dietary modification, and calcium/vitamin D) and an observational study arm. The recruitment goal of WHI was to include a socio-demographically diverse population with racial/ethnic minority groups proportionate to the total minority population of US women aged 50–79 years. This goal was achieved; a diverse population, including 26,045 (17%) women from minority populations, was recruited. Two WHI extension studies conducted additional follow-up on consenting women from 2005–2010 and 2010–2015. Genotyping was available on some WHI participants through the WHI SNP Health Association Resource (SHARe) resource, which used the Affymetrix 6.0 array (~906,600 SNPs, 946,000 copy number variation probes) and on other participants through the MEGA array [44]. Imputation and association analysis was performed separately in individuals with Affymetrix only, MEGA only, and both Affymetrix and MEGA data (S1 Table). For variants with both Affymetrix and MEGA genotypes available, MEGA genotypes were used. In total, 4,318 Hispanic/Latino and 8,494 AA women with blood cell traits were included.

UK Biobank

UK Biobank [46] recruited 500,000 people aged between 40–69 years in 2006–2010, establishing a prospective biobank study to understand the risk factors for common diseases such as cancer, heart disease, stroke, diabetes, and dementia). Participants are being followed-up through routine medical and other health-related records from the UK National Health Service. UK Biobank has genotype data on all enrolled participants, as well as extensive baseline questionnaire and physical measures and stored blood and urine samples. Hematological traits were assayed as previously described [14]. Genotyping on custom Axiom arrays and subsequent quality control has been previously described [47]. Samples were included in our analyses if ancestry self-report was “Black Carribean”, “Black African”,” Black or Black British”, “White and Black Carribean”, “White and Black African”, or “Any Other Black Background”. Variants were selected based on call rate exceeding 95%, HWE p-value exceeding 10−8, and MAF exceeding 0.5%. Subsequently, variants in approximate linkage equilibrium were used to generate ten principle components. Samples were excluded if the first principal component exceeded 0.1 and the second principal component exceeded 0.2, to exclude individuals not clustering with most African ancestry individuals. In total, 6,820 AA participants with blood cell traits were included in the analysis.

Genetic Epidemiology Research on Aging (GERA)

The GERA cohort includes over 100,000 adults who are members of the Kaiser Permanente Medical Care Plan, Northern California Region (KPNC) and consented to research on the genetic and environmental factors that affect health and disease, linking together clinical data from electronic health records, survey data on demographic and behavioral factors, and environmental data with genetic data. The GERA cohort was formed by including all self-reported racial and ethnic minority participants with saliva samples (19%); the remaining participants were drawn sequentially and randomly from non-Hispanic White participants (81%). Genotyping was completed as previously described [48] using 4 different custom Affymetrix Axiom arrays with ethnic-specific content to increase genomic coverage. Principal components analysis was used to characterize genetic structure in this multi-ethnic sample, as previously described [49]. Blood cell traits were extracted from medical records. In individuals with multiple measurements, the first visit with complete white blood cell differential (if any) was used for each participant. Otherwise, the first visit was used. In total, 5,783 Hispanic/Latino and 2,246 AA participants with blood cell traits were included in the analysis.

Jackson Heart Study (JHS)

JHS is a population-based study designed to investigate risk factors for cardiovascular disease in African Americans. JHS recruited 5,306 AA participants age 35–84 from urban and rural areas of the three counties (Hinds, Madison and Rankin) that comprise the Jackson, Mississippi metropolitan area from 2000–2004, including a nested family cohort (≥ 21 years old) and some prior participants from the Atherosclerosis Risk in Communities (ARIC) study [50, 51]. Genotyping was performed using an Affymetrix 6.0 array through NHLBI’s Candidate Gene Association Resource (CARe) consortium [52] in 3,029 individuals, with quality control described previously [53]. Due to the greater JHS sample size in TOPMed freeze 5b (n = 3,082), we extracted SNPs genotyped on Affymetrix 6.0 and which passed CARe consortium quality control in the non-duplicated JHS TOPMed sequenced samples included in the imputation reference panel (821,172 variants which passed TOPMed quality controls used for imputation).

Coronary Artery Risk Development in Young Adults (CARDIA)

The CARDIA study is a longitudinal study of cardiovascular disease risk initiated in 1985–86 in 5,115 AA and European ancestry men and women, then aged 18–30 years. The CARDIA sample was recruited at four sites: Birmingham, AL, Chicago, IL, Minneapolis, MN, and Oakland, CA [54, 55]. Similar to JHS, genotyping was performed through the CARe consortium [52, 53] using an Affymetrix 6.0 array. In total, 1,619 AA participants with blood cell traits were included in the analysis.

Atherosclerosis Risk in Communities (ARIC)

The ARIC study was initiated in 1987, when participants were 45–64 years old, recruiting participants age 45–64 years from 4 field centers (Forsyth County, NC; Jackson, MS; northwestern suburbs of Minneapolis, MN; Washington County, MD) in order to study cardiovascular disease and its risk factors [56], including the participants of self-reported AA ancestry included here. Standardized physical examinations and interviewer-administered questionnaires were conducted at baseline (1987–89), three triennial follow-up examinations, a fifth examination in 2011–13, and a sixth exam in 2016–2017. Genotyping was performed through the CARe consortium Affymetrix 6.0 array [52, 53]. In total, 2,392 AA participants with blood cell traits were included in the analysis.

Imputation and post-imputation quality filtering

We first phased individuals from each cohort separately using eagle [57] with default settings. We subsequently performed haplotype-based imputation using minimac4 [58] using phased haplotypes from TOPMed freeze 5b as reference. We used 100,506 TOPMed freeze 5b whole genome sequences as reference for all cohorts except JHS, for which we used 94,342 TOPMed freeze 5b non-JHS sequences. We additionally imputed HCHS/SOL and JHS using 1000 Genomes Phase 3 [9] and HRC [8] reference panels. Post-imputation quality filtering was performed using a R2 threshold specific to each MAF category to ensure average R2 for variants passing threshold was at least 0.8, following our previous work [4, 59]. Restricting to variants passing post-imputation quality control in at least two cohorts resulted in 34.4–35.8 million variants assessed in the AA cohorts and 26.7–27.2 million assessed in the HA cohorts, depending on the exact sample size of the tested trait. Imputation and association analysis included autosomal variants only. We assessed imputation quality (comparing true and estimated average R2) in three selected 3Mb regions: 16-19Mb region (relative to the start of each chromosome) from chromosomes 3, 12, and 20. Example scripts for imputation quality control are available at https://yunliweb.its.unc.edu/topmed5bimputation/index.php.

Hematological traits

HGB, HCT, WBC and differential were measured in both the discovery data sets (S9 and S16 Tables) and a subset of the TOPMed freeze 5b samples (S10 and S17 Tables) using automated clinical hematology analyzers. Prior to association analyses, we excluded extreme outlier values, notably WBC values >200x109/L (as well as WBC subtype count values in these individuals), HCT >60%, and HGB >20g/dL. For longitudinal cohort studies, all values are from the same exam cycle, chosen based on largest available sample size. WBC traits were log transformed due to their skewed distribution. For all traits, we first derived trait residuals adjusting for age, age squared, sex, and principal components/study specific covariates as needed. Trait residuals were then inverse-normalized prior to analysis.

Association analysis in discovery cohorts

Association analyses were carried out for these variants via EPACTS for all cohorts except for HCHS/SOL, using the q.emmax test to account for relatedness within each cohort. Association tests were performed on inverse normalized residuals (adjusted for age, age squared, sex, and principal components/study specific covariates), further adjusting for kinship matrices constructed in EPACTS using variants with a MAF>1%. Individuals with different starting genotyping platform(s) were also analyzed separately. Inverse-variance weighted meta-analysis were further carried out using GWAMA [60], separately for AAs and Hispanics/Latinos.

Identification and replication of novel associations

To identify putative novel associations, we then filtered out any variant with LD r2 ≥ 0.2 in any ethnic group with any previous reported variant from GWAS, sequencing, or Exome Chip analyses within ±1Mb for a given blood cell trait. We calculated LD in self-reported European ancestry, AA, and Hispanic/Latino individuals from TOPMed freeze 5b. For European and African LD reference panels, we further restricted to individuals with global ancestry estimate ≥0.8. The global ancestry estimates were derived from local ancestry estimates from RFMix [61] using data from the Human Genome Diversity Project (HGDP) [62] as the reference panel with seven populations, namely Sub-Saharan Africa, Central and South Asia, East Asia, Europe, Native America, Oceania, and West Asia and North Africa (Middle East). Global ancestry for each TOPMed individual is defined as the mean local ancestry across all HGDP SNPs. For replication of novel signals, similar to the approach we adopted for the discovery cohorts, we performed association analysis using EPACTS in each contributing cohort and then meta-analyzed with GWAMA.

Supporting information

S1 Fig [pdf]
Comparison of imputation reference panels, for variants with MAF > 1%.

S2 Fig [pdf]
Comparison of well imputed variants included in results from TOPMed (TOPMed freeze 5b), 1000G (the 1000 Genomes Phase 3), and HRC (the Haplotype Reference Consortium).

S3 Fig [pdf]
African ancestry hematocrit analysis Manhattan plot.

S4 Fig [pdf]
African ancestry hemoglobin analysis Manhattan plot.

S5 Fig [pdf]
African ancestry white blood cell count analysis Manhattan plot.

S6 Fig [pdf]
Hispanic/Latino ancestry hematocrit analysis Manhattan plot.

S7 Fig [pdf]
Hispanic/Latino ancestry hemoglobin analysis Manhattan plot.

S8 Fig [pdf]
Hispanic/Latino ancestry white blood cell count analysis Manhattan plot.

S9 Fig [pdf]
African ancestry hematocrit analysis QQ plot.

S10 Fig [pdf]
African ancestry hemoglobin analysis QQ plot.

S11 Fig [pdf]
African ancestry white blood cell count analysis QQ plot.

S12 Fig [pdf]
Hispanic/Latino ancestry hematocrit analysis QQ plot.

S13 Fig [pdf]
Hispanic/Latino ancestry hemoglobin analysis QQ plot.

S14 Fig [pdf]
Hispanic/Latino ancestry white blood cell count analysis QQ plot.

S1 Table [pdf]
Cohorts used for imputation to TOPMed freeze 5b reference panel and subsequent association analysis with hematological traits, including self-identified African ancestry and Hispanic/Latino individuals.

S2 Table [pdf]
Cohorts included in the TOPMed freeze 5b imputation reference panels, with self-reported ancestry.

S3 Table [pdf]
Percentage and number of variants well-imputed with TOPMed freeze5b by chromosome in Jackson Heart Study (JHS) and Hispanic Community Health Study/Study of Latinos (HCHS/SOL).

S4 Table [pdf]
Imputation quality for variants with a minor allele count between 11 and 20 in Jackson Heart Study (JHS).

S5 Table [pdf]
Imputation quality for overall reference panel rare variants (20 or less MAC in TOPMed freeze 5b) in Jackson Heart Study (JHS).

S6 Table [pdf]
Imputation quality for rare variants (20 or less MAC) in Hispanic Community Health Study/Study of Latinos (HCHS/SOL).

S7 Table [pdf]
Imputation quality for rare variants (20 or less MAC) in TOPMed freeze 5b in Hispanic Community Health Study/Study of Latinos (HCHS/SOL).

S8 Table [pdf]
Imputation quality for rare and low frequency variants estimated to be well imputed in stratified by regional background in the Hispanic Community Health Study/Study of Latinos (HCHS/SOL).

S9 Table [pdf]
Demographics, hematological traits, and number of ancestry principal components adjusted for in association analysis models for cohorts imputed to TOPMed freeze 5b reference panel.

S10 Table [pdf]
Demographics, hematological traits, and number of ancestry principal components adjusted for in association analysis for African American cohorts with sequencing and hematological trait data from TOPMed freeze 5b.

S11 Table [pdf]
Overall counts for variants replicated in TOPMed freeze 5b imputed cohorts.

S12 Table [pdf]
Results for previously identified variants in African ancestry and Hispanic/Latino populations in TOPMed freeze 5b imputed samples (included cohorts detailed in and Tables).

S13 Table [pdf]
Imputation of novel variants identified with TOPMed freeze 5b-based imputation using current widely used reference panels from the Haplotype Reference Consortium (HRC) and 1000 Genomes Phase 3, as well as subsequent association analysis results for cohorts where the variants were well-imputed (R2>0.8).

S14 Table [pdf]
Estimated imputation quality for rs33930165 and rs11549407 using 1000G phase 3 and Haplotype Reference Consortium (HRC) as references.

S15 Table [pdf]
Association statistics for the hemoglobin C variant (rs33930165, 11:5227003:C:T) with white blood cell subtypes, adjusting for age, sex, and ancestry principal components.

S16 Table [pdf]
White blood cell subtypes for cohorts imputed to TOPMed freeze 5b reference panel.

S17 Table [pdf]
White blood cell subtypes for African American cohorts with sequencing and hematological trait data from TOPMed freeze 5b.

S18 Table [pdf]
Results for meta-analysis of African ancestry cohorts from sequence kernel association test (SKAT) association results for previously reported genes for hemoglobin (HGB), hematocrit (HCT), or white blood cell count (WBC) using TOPMed freeze 5b, Haplotype Reference Consortium (HRC), and 1000G phase 3 as imputation reference panels.

S19 Table [pdf]
Overall counts for gene results replicated in African ancestry cohorts using TOPMed freeze 5b, 1000G phase 3, and Haplotype Reference Consortium (HRC) as imputation reference panels.

S20 Table [pdf]
Results for meta-analysis of Hispanic/Latino cohorts from sequence kernel association test (SKAT) association results for previously reported genes for hemoglobin (HGB), hematocrit (HCT), or white blood cell count (WBC) using TOPMed freeze 5b, Haplotype Reference Consortium (HRC), and 1000G phase 3 as imputation reference panels.

S21 Table [pdf]
Overall counts for gene results replicated in Hispanic/Latino cohorts using TOPMed freeze 5b, 1000G phase 3, and Haplotype Reference Consortium (HRC) as imputation reference panels.

S1 File [pdf]
TOPMed Banner Authors.

S2 File [pdf]
TOPMed Hematology & Hemostasis Working Group.


Zdroje

1. Auer PL, Johnsen JM, Johnson AD, Logsdon BA, Lange LA, Nalls MA, et al. Imputation of exome sequence variants into population- based samples and blood-cell-trait-associated loci in African Americans: NHLBI GO Exome Sequencing Project. Am J Hum Genet. 2012;91(5):794–808. Epub 2012/10/30. doi: 10.1016/j.ajhg.2012.08.031 23103231.

2. Duan Q, Liu EY, Auer PL, Zhang G, Lange EM, Jun G, et al. Imputation of coding variants in African Americans: better performance using data from the exome sequencing project. Bioinformatics. 2013;29(21):2744–9. Epub 2013/08/21. doi: 10.1093/bioinformatics/btt477 23956302.

3. Lu F-P, Lin K-P, Kuo H-K. Diabetes and the risk of multi-system aging phenotypes: a systematic review and meta-analysis. PloS one. 2009;4(1):e4144. doi: 10.1371/journal.pone.0004144 19127292

4. Liu EY, Buyske S, Aragaki AK, Peters U, Boerwinkle E, Carlson C, et al. Genotype Imputation of MetabochipSNPs Using a Study-Specific Reference Panel of ~4,000 Haplotypes in African Americans From the Women’s Health Initiative. Genet Epidemiol. 2012;36(2):107–17. doi: 10.1002/gepi.21603 22851474

5. Liu EY, Li M, Wang W, Li Y. MaCH-admix: genotype imputation for admixed populations. Genetic epidemiology. 2013;37(1):25–37. Epub 2012/10/18. doi: 10.1002/gepi.21690 23074066.

6. Vergara C, Parker MM, Franco L, Cho MH, Valencia-Duarte AV, Beaty TH, et al. Genotype imputation performance of three reference panels using African ancestry individuals. Hum Genet. 2018;137(4):281–92. Epub 2018/04/11. doi: 10.1007/s00439-018-1881-4 29637265.

7. The International HapMap Consortium. A second generation human haplotype map of over 3.1 million SNPs. Nature. 2007;449:851–61. doi: 10.1038/nature06258 17943122

8. McCarthy S, Das S, Kretzschmar W, Delaneau O, Wood AR, Teumer A, et al. A reference panel of 64,976 haplotypes for genotype imputation. Nat Genet. 2016;48(10):1279–83. Epub 2016/08/23. doi: 10.1038/ng.3643 27548312.

9. The 1000 Genomes Project Consortium. A global reference for human genetic variation. Nature. 2015;526(7571):68–74. Epub 2015/10/04. doi: 10.1038/nature15393 26432245.

10. Mathias RA, Taub MA, Gignoux CR, Fu W, Musharoff S, O’Connor TD, et al. A continuum of admixture in the Western Hemisphere revealed by the African Diaspora genome. Nature communications. 2016;7:12522. Epub 2016/10/12. doi: 10.1038/ncomms12522 27725671 other authors declare no competing financial interests.

11. Crosslin DR, McDavid A, Weston N, Nelson SC, Zheng X, Hart E, et al. Genetic variants associated with the white blood cell count in 13,923 subjects in the eMERGE Network. Hum Genet. 2012;131(4):639–52. Epub 2011/11/01. doi: 10.1007/s00439-011-1103-9 22037903.

12. Whitfield JB, Martin NG, Rao DC. Genetic and environmental influences on the size and number of cells in the blood. Genetic epidemiology. 1985;2(2):133–44. doi: 10.1002/gepi.1370020204 4054596

13. Garner C, Tatu T, Reittie JE, Littlewood T, Darley J, Cervino S, et al. Genetic influences on F cells and other hematologic variables: a twin heritability study. Blood. 2000;95(1):342–6. Epub 1999/12/23. 10607722.

14. Astle WJ, Elding H, Jiang T, Allen D, Ruklisa D, Mann AL, et al. The Allelic Landscape of Human Blood Cell Trait Variation and Links to Common Complex Disease. Cell. 2016;167(5):1415–29 e19. Epub 2016/11/20. doi: 10.1016/j.cell.2016.10.042 27863252.

15. van der Harst P, Zhang W, Mateo Leach I, Rendon A, Verweij N, Sehmi J, et al. Seventy-five genetic loci influencing the human red blood cell. Nature. 2012;492(7429):369–75. Epub 2012/12/12. doi: 10.1038/nature11677 23222517.

16. Mousas A, Ntritsos G, Chen MH, Song C, Huffman JE, Tzoulaki I, et al. Rare coding variants pinpoint genes that control human hematological traits. PLoS Genet. 2017;13(8):e1006925. Epub 2017/08/09. doi: 10.1371/journal.pgen.1006925 28787443.

17. Chami N, Chen MH, Slater AJ, Eicher JD, Evangelou E, Tajuddin SM, et al. Exome Genotyping Identifies Pleiotropic Variants Associated with Red Blood Cell Traits. Am J Hum Genet. 2016. Epub 2016/06/28. doi: 10.1016/j.ajhg.2016.05.007 27346685.

18. Eicher JD, Chami N, Kacprowski T, Nomura A, Chen MH, Yanek LR, et al. Platelet-Related Variants Identified by Exomechip Meta-analysis in 157,293 Individuals. Am J Hum Genet. 2016. Epub 2016/06/28. doi: 10.1016/j.ajhg.2016.05.005 27346686.

19. Tajuddin SM, Schick UM, Eicher JD, Chami N, Giri A, Brody JA, et al. Large-Scale Exome-wide Association Analysis Identifies Loci for White Blood Cell Traits and Pleiotropy with Immune-Mediated Diseases. Am J Hum Genet. 2016. Epub 2016/06/28. doi: 10.1016/j.ajhg.2016.05.003 27346689.

20. Hodonsky CJ, Jain D, Schick UM, Morrison JV, Brown L, McHugh CP, et al. Genome-wide association study of red blood cell traits in Hispanics/Latinos: The Hispanic Community Health Study/Study of Latinos. PLoS genetics. 2017;13(4):e1006760. Epub 2017/04/30. doi: 10.1371/journal.pgen.1006760 28453575.

21. Group. CCHW. Meta-analysis of rare and common exome chip variants identifies S1PR4 and other loci influencing blood cell traits. Nat Genet. 2016;48(8):867–76. Epub 2016/07/12. doi: 10.1038/ng.3607 27399967.

22. Lo KS, Wilson JG, Lange LA, Folsom AR, Galarneau G, Ganesh SK, et al. Genetic association analysis highlights new loci that modulate hematological trait variation in Caucasians and African Americans. Hum Genet. 2011;129(3):307–17. Epub 2010/12/15. doi: 10.1007/s00439-010-0925-1 21153663.

23. Tournamille C, Colin Y, Cartron JP, Le Van Kim C. Disruption of a GATA motif in the Duffy gene promoter abolishes erythroid gene expression in Duffy-negative individuals. Nat Genet. 1995;10(2):224–8. Epub 1995/06/01. doi: 10.1038/ng0695-224 7663520.

24. Chen Z, Tang H, Qayyum R, Schick UM, Nalls MA, Handsaker R, et al. Genome-wide association analysis of red blood cell traits in African Americans: the COGENT Network. Hum Mol Genet. 2013;22(12):2529–38. Epub 2013/03/01. doi: 10.1093/hmg/ddt087 23446634.

25. Li J, Glessner JT, Zhang H, Hou C, Wei Z, Bradfield JP, et al. GWAS of blood cell traits identifies novel associated loci and epistatic interactions in Caucasian and African-American children. Hum Mol Genet. 2013;22(7):1457–64. Epub 2012/12/25. doi: 10.1093/hmg/dds534 23263863.

26. van Rooij FJA, Qayyum R, Smith AV, Zhou Y, Trompet S, Tanaka T, et al. Genome-wide Trans-ethnic Meta-analysis Identifies Seven Genetic Loci Influencing Erythrocyte Traits and a Role for RBPMS in Erythropoiesis. Am J Hum Genet. 2017;100(1):51–63. Epub 2016/12/27. doi: 10.1016/j.ajhg.2016.11.016 28017375.

27. Jain D, Hodonsky CJ, Schick UM, Morrison JV, Minnerath S, Brown L, et al. Genome-wide association of white blood cell counts in Hispanic/Latino Americans: the Hispanic Community Health Study/Study of Latinos. Hum Mol Genet. 2017;26(6):1193–204. Epub 2017/02/06. doi: 10.1093/hmg/ddx024 28158719.

28. Polfus LM, Khajuria RK, Schick UM, Pankratz N, Pazoki R, Brody JA, et al. Whole-Exome Sequencing Identifies Loci Associated with Blood Cell Traits and Reveals a Role for Alternative GFI1B Splice Variants in Human Hematopoiesis. Am J Hum Genet. 2016;99(2):481–8. Epub 2016/08/04. doi: 10.1016/j.ajhg.2016.06.016 27486782.

29. Keller MF, Reiner AP, Okada Y, van Rooij FJ, Johnson AD, Chen MH, et al. Trans-ethnic meta-analysis of white blood cell phenotypes. Hum Mol Genet. 2014;23(25):6944–60. Epub 2014/08/07. doi: 10.1093/hmg/ddu401 25096241.

30. Reiner AP, Lettre G, Nalls MA, Ganesh SK, Mathias R, Austin MA, et al. Genome-wide association study of white blood cell count in 16,388 African Americans: the continental origins and genetic epidemiology network (COGENT). PLoS genetics. 2011;7(6):e1002108. Epub 2011/07/09. doi: 10.1371/journal.pgen.1002108 21738479.

31. Pulit SL, de With SA, de Bakker PI. Resetting the bar: Statistical significance in whole-genome sequencing-based association studies of global populations. Genetic epidemiology. 2017;41(2):145–51. Epub 2016/12/19. doi: 10.1002/gepi.22032 27990689.

32. Trecartin RF, Liebhaber SA, Chang JC, Lee KY, Kan YW, Furbetta M, et al. beta zero thalassemia in Sardinia is caused by a nonsense mutation. The Journal of clinical investigation. 1981;68(4):1012–7. Epub 1981/10/01. doi: 10.1172/JCI110323 6457059.

33. Rosatelli MC, Dozy A, Faa V, Meloni A, Sardu R, Saba L, et al. Molecular characterization of beta-thalassemia in the Sardinian population. Am J Hum Genet. 1992;50(2):422–6. Epub 1992/02/01. 1734721.

34. Perea FJ, Magana MT, Cobian JG, Sanchez-Lopez JY, Chavez ML, Zamudio G, et al. Molecular spectrum of beta-thalassemia in the Mexican population. Blood cells, molecules & diseases. 2004;33(2):150–2. Epub 2004/08/19. doi: 10.1016/j.bcmd.2004.06.001 15315794.

35. Silva AN, Cardoso GL, Cunha DA, Diniz IG, Santos SE, Andrade GB, et al. The Spectrum of beta-Thalassemia Mutations in a Population from the Brazilian Amazon. Hemoglobin. 2016;40(1):20–4. Epub 2015/09/16. doi: 10.3109/03630269.2015.1083443 26372288.

36. Key NS, Connes P, Derebail VK. Negative health implications of sickle cell trait in high income countries: from the football field to the laboratory. British journal of haematology. 2015;170(1):5–14. Epub 2015/03/11. doi: 10.1111/bjh.13363 25754217.

37. Graffeo L, Vitrano A, Scondotto S, Dardanoni G, Pollina Addario WS, Giambona A, et al. beta-Thalassemia heterozygote state detrimentally affects health expectation. European journal of internal medicine. 2018;54:76–80. Epub 2018/06/24. doi: 10.1016/j.ejim.2018.06.009 29934240.

38. Galanello R, Origa R. Beta-thalassemia. Orphanet journal of rare diseases. 2010;5:11-. doi: 10.1186/1750-1172-5-11 20492708.

39. Fairhurst RM, Casella JF. Images in clinical medicine. Homozygous hemoglobin C disease. Q1. 2004;350(26):e24. Epub 2004/06/25. doi: 10.1056/NEJMicm030486 15215497.

40. Sorlie PD, Aviles-Santa LM, Wassertheil-Smoller S, Kaplan RC, Daviglus ML, Giachello AL, et al. Design and implementation of the Hispanic Community Health Study/Study of Latinos. Ann Epidemiol. 2010;20(8):629–41. Epub 2010/07/09. doi: 10.1016/j.annepidem.2010.03.015 20609343.

41. Daviglus ML, Talavera GA, Aviles-Santa ML, Allison M, Cai J, Criqui MH, et al. Prevalence of major cardiovascular risk factors and cardiovascular diseases among Hispanic/Latino individuals of diverse backgrounds in the United States. Jama. 2012;308(17):1775–84. Epub 2012/11/03. doi: 10.1001/jama.2012.14517 23117778.

42. Lavange LM, Kalsbeek WD, Sorlie PD, Aviles-Santa LM, Kaplan RC, Barnhart J, et al. Sample design and cohort selection in the Hispanic Community Health Study/Study of Latinos. Ann Epidemiol. 2010;20(8):642–9. Epub 2010/07/09. doi: 10.1016/j.annepidem.2010.05.006 20609344.

43. Conomos MP, Laurie CA, Stilp AM, Gogarten SM, McHugh CP, Nelson SC, et al. Genetic Diversity and Association Studies in US Hispanic/Latino Populations: Applications in the Hispanic Community Health Study/Study of Latinos. Am J Hum Genet. 2016;98(1):165–84. Epub 2016/01/11. doi: 10.1016/j.ajhg.2015.12.001 26748518.

44. Wojcik G, Graff M, Nishimura KK, Tao R, Haessler J, Gignoux CR, et al. The PAGE Study: How Genetic Diversity Improves Our Understanding of the Architecture of Complex Traits. bioRxiv. 2018:188094. doi: 10.1101/188094

45. The Women’s Health Initiative Study Group. Design of the Women’s Health Initiative clinical trial and observational study. Control Clin Trials. 1998;19(1):61–109. Epub 1998/03/11. doi: 10.1016/s0197-2456(97)00078-0 9492970.

46. UK Biobank. UK Biobank: rationale, design and development of a large-scale prospective resource. 2007. http://www.ukbiobank.ac.uk/resources/.

47. Bycroft C, Freeman C, Petkova D, Band G, Elliott LT, Sharp K, et al. The UK Biobank resource with deep phenotyping and genomic data. Nature. 2018;562(7726):203–9. doi: 10.1038/s41586-018-0579-z 30305743

48. Kvale MN, Hesselson S, Hoffmann TJ, Cao Y, Chan D, Connell S, et al. Genotyping Informatics and Quality Control for 100,000 Subjects in the Genetic Epidemiology Research on Adult Health and Aging (GERA) Cohort. Genetics. 2015;200(4):1051–60. Epub 2015/06/21. doi: 10.1534/genetics.115.178905 26092718.

49. Banda Y, Kvale MN, Hoffmann TJ, Hesselson SE, Ranatunga D, Tang H, et al. Characterizing Race/Ethnicity and Genetic Ancestry for 100,000 Subjects in the Genetic Epidemiology Research on Adult Health and Aging (GERA) Cohort. Genetics. 2015;200(4):1285–95. doi: 10.1534/genetics.115.178616 26092716

50. Taylor HA Jr., Wilson JG, Jones DW, Sarpong DF, Srinivasan A, Garrison RJ, et al. Toward resolution of cardiovascular health disparities in African Americans: design and methods of the Jackson Heart Study. Ethn Dis. 2005;15(4 Suppl 6):S6-4-17. Epub 2005/12/02. 16320381.

51. Wilson JG, Rotimi CN, Ekunwe L, Royal CD, Crump ME, Wyatt SB, et al. Study design for genetic analysis in the Jackson Heart Study. Ethn Dis. 2005;15(4 Suppl 6):S6-30-7. Epub 2005/12/02. 16317983.

52. Musunuru K, Lettre G, Young T, Farlow DN, Pirruccello JP, Ejebe KG, et al. Candidate gene association resource (CARe): design, methods, and proof of concept. Circulation Cardiovascular genetics. 2010;3(3):267–75. Epub 2010/04/20. doi: 10.1161/CIRCGENETICS.109.882696 20400780.

53. Lettre G, Palmer CD, Young T, Ejebe KG, Allayee H, Benjamin EJ, et al. Genome-wide association study of coronary heart disease and its risk factors in 8,090 African Americans: the NHLBI CARe Project. PLoS genetics. 2011;7(2):e1001300. Epub 2011/02/25. doi: 10.1371/journal.pgen.1001300 21347282.

54. Friedman GD, Cutter GR, Donahue RP, Hughes GH, Hulley SB, Jacobs DR Jr., et al. CARDIA: study design, recruitment, and some characteristics of the examined subjects. J Clin Epidemiol. 1988;41(11):1105–16. Epub 1988/01/01. doi: 10.1016/0895-4356(88)90080-7 3204420.

55. Cutter GR, Burke GL, Dyer AR, Friedman GD, Hilner JE, Hughes GH, et al. Cardiovascular risk factors in young adults. The CARDIA baseline monograph. Control Clin Trials. 1991;12(1 Suppl):1S–77S. Epub 1991/02/11. doi: 10.1016/0197-2456(91)90002-4 1851696.

56. The Atherosclerosis Risk in Communities (ARIC) Study: design and objectives. The ARIC investigators. Am J Epidemiol. 1989;129(4):687–702. Epub 1989/04/01. 2646917.

57. Loh PR, Danecek P, Palamara PF, Fuchsberger C, AR Y, KF H, et al. Reference-based phasing using the Haplotype Reference Consortium panel. Nat Genet. 2016;48(11):1443–8. Epub 2016/10/28. doi: 10.1038/ng.3679 27694958.

58. Das S, Forer L, Schonherr S, Sidore C, Locke AE, Kwong A, et al. Next-generation genotype imputation service and methods. Nat Genet. 2016;48(10):1284–7. Epub 2016/08/30. doi: 10.1038/ng.3656 27571263.

59. Duan Q, Liu EY, Croteau-Chonka DC, Mohlke KL, Li Y. A comprehensive SNP and indel imputability database. Bioinformatics. 2013;29(4):528–31. Epub 2013/01/08. doi: 10.1093/bioinformatics/bts724 23292738.

60. Magi R, Lindgren CM, Morris AP. Meta-analysis of sex-specific genome-wide association studies. Genet Epidemiol. 2010;34(8):846–53. Epub 2010/11/26. doi: 10.1002/gepi.20540 21104887.

61. Maples Brian K, Gravel S, Kenny Eimear E, Bustamante Carlos D. RFMix: A Discriminative Modeling Approach for Rapid and Robust Local-Ancestry Inference. The American Journal of Human Genetics. 2013;93(2):278–88. http://dx.doi.org/10.1016/j.ajhg.2013.06.020 23910464

62. Li JZ, Absher DM, Tang H, Southwick AM, Casto AM, Ramachandran S, et al. Worldwide human relationships inferred from genome-wide patterns of variation. Science. 2008;319(5866):1100–4. Epub 2008/02/23. doi: 10.1126/science.1153717 18292342.

Štítky
Genetika Reprodukční medicína

Článek vyšel v časopise

PLOS Genetics


2019 Číslo 12
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy Podcasty Doporučená témata Časopisy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#