Search

Platelets play a key role in thrombosis and hemostasis. Platelet count (PLT) and mean platelet volume (MPV) are highly heritable quantitative traits, with hundreds of genetic signals previously identified, mostly in European ancestry populations. We here utilize whole genome sequencing (WGS) from NHLBI's Trans-Omics for Precision Medicine initiative (TOPMed) in a large multi-ethnic sample to further explore common and rare variation contributing to PLT (n = 61 200) and MPV (n = 23 485). We identified and replicated secondary signals at MPL (rs532784633) and PECAM1 (rs73345162), both more common in African ancestry populations. We also observed rare variation in Mendelian platelet-related disorder genes influencing variation in platelet traits in TOPMed cohorts (not enriched for blood disorders). For example, association of GP9 with lower PLT and higher MPV was partly driven by a pathogenic Bernard-Soulier syndrome variant (rs5030764, p.Asn61Ser), and the signals at TUBB1 and CD36 were partly driven by loss of function variants not annotated as pathogenic in ClinVar (rs199948010 and rs571975065). However, residual signal remained for these gene-based signals after adjusting for lead variants, suggesting that additional variants in Mendelian genes with impacts in general population cohorts remain to be identified. Gene-based signals were also identified at several genome-wide association study identified loci for genes not annotated for Mendelian platelet disorders (PTPRH, TET2, CHEK2), with somatic variation driving the result at TET2. These results highlight the value of WGS in populations of diverse genetic ancestry to identify novel regulatory and coding signals, even for well-studied traits like platelet traits., (© The Author(s) 2021. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Many common and rare variants associated with hematologic traits have been discovered through imputation on large-scale reference panels. However, the majority of genome-wide association studies (GWASs) have been conducted in Europeans, and determining causal variants has proved challenging. We performed a GWAS of total leukocyte, neutrophil, lymphocyte, monocyte, eosinophil, and basophil counts generated from 109,563,748 variants in the autosomes and the X chromosome in the Trans-Omics for Precision Medicine (TOPMed) program, which included data from 61,802 individuals of diverse ancestry. We discovered and replicated 7 leukocyte trait associations, including (1) the association between a chromosome X, pseudo-autosomal region (PAR), noncoding variant located between cytokine receptor genes (CSF2RA and CLRF2) and lower eosinophil count; and (2) associations between single variants found predominantly among African Americans at the S1PR3 (9q22.1) and HBB (11p15.4) loci and monocyte and lymphocyte counts, respectively. We further provide evidence indicating that the newly discovered eosinophil-lowering chromosome X PAR variant might be associated with reduced susceptibility to common allergic diseases such as atopic dermatitis and asthma. Additionally, we found a burden of very rare FLT3 (13q12.2) variants associated with monocyte counts. Together, these results emphasize the utility of whole-genome sequencing in diverse samples in identifying associations missed by European-ancestry-driven GWASs., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2021 American Society of Human Genetics. All rights reserved.)

Genotype-phenotype association studies often combine phenotype data from multiple studies to increase statistical power. Harmonization of the data usually requires substantial effort due to heterogeneity in phenotype definitions, study design, data collection procedures, and data-set organization. Here we describe a centralized system for phenotype harmonization that includes input from phenotype domain and study experts, quality control, documentation, reproducible results, and data-sharing mechanisms. This system was developed for the National Heart, Lung, and Blood Institute's Trans-Omics for Precision Medicine (TOPMed) program, which is generating genomic and other -omics data for more than 80 studies with extensive phenotype data. To date, 63 phenotypes have been harmonized across thousands of participants (recruited in 1948-2012) from up to 17 studies per phenotype. Here we discuss challenges in this undertaking and how they were addressed. The harmonized phenotype data and associated documentation have been submitted to National Institutes of Health data repositories for controlled access by the scientific community. We also provide materials to facilitate future harmonization efforts by the community, which include 1) the software code used to generate the 63 harmonized phenotypes, enabling others to reproduce, modify, or extend these harmonizations to additional studies, and 2) the results of labeling thousands of phenotype variables with controlled vocabulary terms., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Johns Hopkins Bloomberg School of Public Health.)

Whole-genome sequencing (WGS), a powerful tool for detecting novel coding and non-coding disease-causing variants, has largely been applied to clinical diagnosis of inherited disorders. Here we leveraged WGS data in up to 62,653 ethnically diverse participants from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program and assessed statistical association of variants with seven red blood cell (RBC) quantitative traits. We discovered 14 single variant-RBC trait associations at 12 genomic loci, which have not been reported previously. Several of the RBC trait-variant associations (RPN1, ELL2, MIDN, HBB, HBA1, PIEZO1, and G6PD) were replicated in independent GWAS datasets imputed to the TOPMed reference panel. Most of these discovered variants are rare/low frequency, and several are observed disproportionately among non-European Ancestry (African, Hispanic/Latino, or East Asian) populations. We identified a 3 bp indel p.Lys2169del (g.88717175_88717177TCT[4]) (common only in the Ashkenazi Jewish population) of PIEZO1, a gene responsible for the Mendelian red cell disorder hereditary xerocytosis (MIM: 194380), associated with higher mean corpuscular hemoglobin concentration (MCHC). In stepwise conditional analysis and in gene-based rare variant aggregated association analysis, we identified several of the variants in HBB, HBA1, TMPRSS6, and G6PD that represent the carrier state for known coding, promoter, or splice site loss-of-function variants that cause inherited RBC disorders. Finally, we applied base and nuclease editing to demonstrate that the sentinel variant rs112097551 (nearest gene RPN1) acts through a cis-regulatory element that exerts long-range control of the gene RUVBL1 which is essential for hematopoiesis. Together, these results demonstrate the utility of WGS in ethnically diverse population-based samples and gene editing for expanding knowledge of the genetic architecture of quantitative hematologic traits and suggest a continuum between complex trait and Mendelian red cell disorders., (Copyright © 2021 American Society of Human Genetics. All rights reserved.)

The Trans-Omics for Precision Medicine (TOPMed) programme seeks to elucidate the genetic architecture and biology of heart, lung, blood and sleep disorders, with the ultimate goal of improving diagnosis, treatment and prevention of these diseases. The initial phases of the programme focused on whole-genome sequencing of individuals with rich phenotypic data and diverse backgrounds. Here we describe the TOPMed goals and design as well as the available resources and early insights obtained from the sequence data. The resources include a variant browser, a genotype imputation server, and genomic and phenotypic data that are available through dbGaP (Database of Genotypes and Phenotypes) 1 . In the first 53,831 TOPMed samples, we detected more than 400 million single-nucleotide and insertion or deletion variants after alignment with the reference genome. Additional previously undescribed variants were detected through assembly of unmapped reads and customized analysis in highly variable loci. Among the more than 400 million detected variants, 97% have frequencies of less than 1% and 46% are singletons that are present in only one individual (53% among unrelated individuals). These rare variants provide insights into mutational processes and recent human evolutionary history. The extensive catalogue of genetic variation in TOPMed studies provides unique opportunities for exploring the contributions of rare and noncoding sequence variants to phenotypic variation. Furthermore, combining TOPMed haplotypes with modern imputation methods improves the power and reach of genome-wide association studies to include variants down to a frequency of approximately 0.01%.

Objective: A recently published genome wide association study of abdominal aortic aneurysms (AAA), based on pooled case control data of European ancestry, identified four new loci for AAA: SMYD2 (top single nucleotide polymorphism [SNP] rs1795061), LINC00540 (rs9316871), PCIF1/MMP9/ZNF335 (rs3827066), and ERG (rs2836411). Of the four, rs1795061 and rs2836411 showed significant heterogeneity across studies and the p value for rs9316871 did not reach the genome wide significance threshold until discovery and replication data were pooled together in that study. The objective of this study was to replicate these newly identified genetic associations for AAA in a US based prospective cohort study, the Atherosclerosis Risk in Communities (ARIC) Study, and a Greece based case control study., Methods: ARIC identified 408 clinically diagnosed AAAs among 8 962 individuals of European ancestry during a median of 22 years of follow up. The Greek case control study included 341 AAAs of European ancestry recruited in a tertiary referral centre and 292 geographically and ethnically matched controls recruited from the same institution. A Cox proportional hazards model was used to analyse the ARIC data and logistic regression to analyse the Greek data., Results: In ARIC, rs9316871 and rs3827066 were significantly associated with AAA risk (HR [p] was 0.77 [.004] and 1.22 [.03], respectively), rs2836411 was associated at borderline significance (1.13 [.08]), whereas rs1795061 was not associated (p = .55). In the Greek case control study, rs1795061 and rs2836411 were significantly associated with AAA (OR [p] was 1.66 [< .001] and 1.29 [.04], respectively), whereas rs9316871 was not (p = .81). Genotyping of rs3827066 did not succeed. In the meta-analysis of the two studies, the association for rs9316871and rs2836411 was statistically significant and consistent between the two studies: p = .02 and .007, respectively., Conclusions: Associations between rs9316871and rs2836411 and AAA risk were replicated in the meta-analysis of the two independent cohorts, providing further support for the importance of these loci in the aetiology of AAA., (Copyright © 2019 European Society for Vascular Surgery. Published by Elsevier B.V. All rights reserved.)

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., Competing Interests: Edwin K Silverman and Michael H Cho have received grant support from GSK, MHC has received consulting fees from Genentech. Scott T. Weiss and Kathleen C. Barnes received royalties from UpToDate. Patrick T. Ellinor is supported by a grant from Bayer AG to the Broad Institute focused on the genetics and therapeutics of cardiovascular diseases, and has also served on advisory boards or consulted for Bayer AG, Quest Diagnostics and Novartis. Steven A Lubitz receives sponsored research support from Bristol Myers Squibb / Pfizer, Bayer HealthCare, and Boehringer Ingelheim, and has consulted for Abbott, Quest Diagnostics, Bristol Myers Squibb / Pfizer. Other authors declared no conflicts of interest.

Importance: Mitochondrial dysfunction is a core component of the aging process and may play a key role in atherosclerotic cardiovascular disease. Mitochondrial DNA copy number (mtDNA-CN), which represents the number of mitochondria per cell and number of mitochondrial genomes per mitochondrion, is an indirect biomarker of mitochondrial function., Objective: To determine whether mtDNA-CN, measured in an easily accessible tissue (buffy coat/circulating leukocytes), can improve risk classification for cardiovascular disease (CVD) and help guide initiation of statin therapy for primary prevention of CVD., Design, Setting, and Participants: Prospective, population-based cohort analysis including 21 870 participants (20 163 free from CVD at baseline) from 3 studies: Cardiovascular Health Study (CHS), Atherosclerosis Risk in Communities Study (ARIC), and Multiethnic Study of Atherosclerosis (MESA). The mean follow-up was 13.5 years. The study included 11 153 participants from ARIC, 4830 from CHS, and 5887 from MESA. Analysis of the data was conducted from March 10, 2014, to January 29, 2017., Exposures: Mitochondrial DNA-CN measured from buffy coat/circulating leukocytes., Main Outcomes and Measures: Incident CVD, which combines coronary heart disease, defined as the first incident myocardial infarction or death owing to coronary heart disease, and stroke, defined as the first nonfatal stroke or death owing to stroke., Results: Of the 21 870 participants, the mean age was 62.4 years (ARIC, 57.9 years; MESA, 62.4 years; and CHS, 72.5 years), and 54.7% of participants were women. The hazard ratios for incident coronary heart disease, stroke, and CVD associated with a 1-SD decrease in mtDNA-CN were 1.29 (95% CI, 1.24-1.33), 1.11 (95% CI, 1.06-1.16), and 1.23 (95% CI, 1.19-1.26). The associations persisted after adjustment for traditional CVD risk factors. Addition of mtDNA-CN to the 2013 American College of Cardiology/American Heart Association Pooled Cohorts Equations for estimating 10-year hard atherosclerosis CVD risk was associated with improved risk classification (continuous net reclassification index, 0.194; 95% CI, 0.130-0.258; P < .001). Mitochondrial DNA-CN further improved sensitivity and specificity for the 2013 American College of Cardiology/American Heart Association recommendations on initiating statin therapy for primary prevention of ASCVD (net 221 individuals appropriately downclassified and net 15 individuals appropriately upclassified)., Conclusions and Relevance: Mitochondrial DNA-CN was independently associated with incident CVD in 3 large prospective studies and may have potential clinical utility in improving CVD risk classification.

Objective: Previous studies have identified common genetic variants in 4 chromosomal regions that together account for 14% to 15% of the variance in circulating levels of protein C. To further characterize the genetic architecture of protein C, we obtained denser coverage at some loci, extended investigation of protein C to low-frequency and rare variants, and searched for new associations in genes known to influence protein C., Approach and Results: Genetic associations with protein C antigen level were evaluated in ≤10 778 European and 3190 black participants aged 45 to 64 years. Analyses included >26 million autosomal variants available after imputation to the 1000 Genomes reference panel along with additional low-frequency and rare variants directly genotyped using the Illumina ITMAT-Broad-CARe chip and Illumina HumanExome BeadChip. Genome-wide significant associations ( P <5×10 -8 ) were found for common variants in the GCKR , PROC , BAZ1B , and PROCR-EDEM2 regions in whites and PROC and PROCR-EDEM2 regions in blacks, confirming earlier findings. In a novel finding, the low-density lipoprotein cholesterol-lowering allele of rs12740374, located in the CELSR2-PSRC1-SORT1 region, was associated with lower protein C level in both whites and blacks, reaching genome-wide significance in a meta-analysis combining results from both groups ( P =1.4×10 -9 ). To further investigate a possible link between lipid metabolism and protein C level, we conducted Mendelian randomization analyses using 185 lipid-related genetic variants as instrumental variables. The results indicated that triglycerides, and possibly low-density lipoprotein cholesterol, influence protein C levels., Conclusions: Discovery of variants influencing circulating protein C levels in the CELSR2-PSRC1-SORT1 region may indicate a novel genetic link between lipoprotein metabolism and hemostasis., (© 2017 American Heart Association, Inc.)

Background: C-reactive protein (CRP) is a biomarker of inflammation. Genome-wide association studies (GWAS) have identified single-nucleotide polymorphisms (SNPs) associated with CRP concentrations and inflammation-related traits such as cardiovascular disease, type 2 diabetes mellitus, and obesity. We aimed to replicate previous CRP-SNP associations, assess whether these associations generalize to additional race/ethnicity groups, and evaluate inflammation-related SNPs for a potentially pleiotropic association with CRP., Methods and Results: We selected and analyzed 16 CRP-associated and 250 inflammation-related GWAS SNPs among 40 473 African American, American Indian, Asian/Pacific Islander, European American, and Hispanic participants from 7 studies collaborating in the Population Architecture using Genomics and Epidemiology (PAGE) study. Fixed-effect meta-analyses combined study-specific race/ethnicity-stratified linear regression estimates to evaluate the association between each SNP and high-sensitivity CRP. Overall, 18 SNPs in 8 loci were significantly associated with CRP (Bonferroni-corrected P<3.1×10(-3) for replication, P<2.0×10(-4) for pleiotropy): Seven of these were specific to European Americans, while 9 additionally generalized to African Americans (1), Hispanics (5), or both (3); 1 SNP was seen only in African Americans and Hispanics. Two SNPs in the CELSR2/PSRC1/SORT1 locus showed a potentially novel association with CRP: rs599839 (P=2.0×10(-6)) and rs646776 (P=3.1×10(-5))., Conclusions: We replicated 16 SNP-CRP associations, 10 of which generalized to African Americans and/or Hispanics. We also identified potentially novel pleiotropic associations with CRP for two SNPs previously associated with coronary artery disease and/or low-density lipoprotein-cholesterol. These findings demonstrate the benefit of evaluating genotype-phenotype associations in multiple race/ethnicity groups and looking for pleiotropic relationships among SNPs previously associated with related phenotypes.

A genome screen to identify genes influencing the age at Parkinson disease (PD) onset was completed using 276 families without parkin mutations. Significant evidence of linkage to chromosome 2p near the PARK3 locus (logarithm of odds [lod] = 4.8) was observed. Evidence of linkage was also detected to chromosomes 1q (lod = 3.0) and 8q (lod = 2.6). These data suggest that the genes influencing age at PD onset likely differ from those that contribute to PD susceptibility.