Genome-wide association study of red blood cell traits in Hispanics/Latinos: The Hispanic Community Health Study/Study of Latinos

Advancements in genetic research by the Hispanic Community Health Study/Study of Latinos: A 10-year retrospective review

The Hispanic Community Health Study/Study of Latinos (HCHS/SOL) is a multicenter, longitudinal cohort study designed to evaluate environmental, lifestyle, and genetic risk factors as they relate to cardiometabolic and other chronic diseases among Hispanic/Latino populations in the United States. Since the study's inception in 2008, as a result of the study's robust genetic measures, HCHS/SOL has facilitated major contributions to the field of genetic research. This 10-year retrospective review highlights the major findings for genotype-phenotype relationships and advancements in statistical methods owing to the HCHS/SOL. Furthermore, we discuss the ethical and societal challenges of genetic research, especially among Hispanic/Latino adults in the United States. Continued genetic research, ancillary study expansion, and consortia collaboration through HCHS/SOL will further drive knowledge and advancements in human genetics research.

Antihypertensive drug targets and breast cancer risk: a two-sample Mendelian randomization study

Findings on the correlation between the use of antihypertensive medication and the risk of breast cancer (BC) have been inconsistent. We performed a two-sample Mendelian randomization (MR) using instrumental variables to proxy changes in gene expressions of antihypertensive medication targets to interrogate this. Genetic instruments for expression of antihypertensive drug target genes were identified with expression quantitative trait loci in blood, which should be associated with systolic blood pressure to proxy for the effect of antihypertensive drug. The association between genetic variants and BC risk were obtained from genome-wide association study summary statistics. The summary-based MR was employed to estimate the drug effects on BC risk. We further performed sensitivity analyses to confirm the discovered MR associations such as assessment of horizontal pleiotropy, colocalization, and multiple tissue enrichment analyses. The overall BC risk was only associated with SLC12A2 gene expression at a Bonferroni-corrected threshold. One standard deviation (SD) decrease of SLC12A2 gene expression in blood was associated with a decrease of 1.12 (95%CI, 0.80-1.58) mmHg of systolic blood pressure, but a 16% increased BC risk (odds ratio, 1.16, 95% confidential interval, 1.06-1.28). This signal was further observed for estrogen receptor positive (ER +) BC (1.17, 1.06-1.28). In addition, one SD decrease in expression of PDE1B in blood was associated with 7% decreased risk of ER + BC (0.93, 0.90-0.97). We detected no evidence of horizontal pleiotropy for these associations and the probability of the causal variants being shared between the gene expression and BC risk was 81.5, 40.5 and 66.8%, respectively. No significant association was observed between other target gene expressions and BC risk. Changes in expression of SLC12A2 and PDE1B mediated possibly via antihypertensive drugs may result in increased and decreased BC risk, respectively.

Red cell distribution width and its polygenic score in relation to mortality and cardiometabolic outcomes

Introduction

Elevated red cell distribution width (RDW) has been associated with a range of health outcomes. This study aims to examine prognostic and etiological roles of RDW levels, both phenotypic and genetic predisposition, in predicting cardiovascular outcomes, diabetes, chronic kidney disease (CKD) and mortality.

Methods

We studied 27,141 middle-aged adults from the Malmö Diet and Cancer study (MDCS) with a mean follow up of 21 years. RDW was measured with a hematology analyzer on whole blood samples. Polygenic scores for RDW (PGS-RDW) were constructed for each participant using genetic data in MDCS and published summary statistics from genome-wide association study of RDW (n = 408,112). Cox proportional hazards regression was used to assess associations between RDW, PGS-RDW and cardiovascular outcomes, diabetes, CKD and mortality, respectively.

Results

PGS-RDW was significantly associated with RDW (Pearson's correlation coefficient = 0.133, p < 0.001). RDW was significantly associated with incidence of stroke (hazard ratio (HR) per 1 standard deviation = 1.06, 95% confidence interval (CI): 1.02-1.10, p = 0.003), atrial fibrillation (HR = 1.09, 95% CI: 1.06-1.12, p < 0.001), heart failure (HR = 1.13, 95% CI: 1.08-1.19, p < 0.001), venous thromboembolism (HR = 1.21, 95% CI: 1.15-1.28, p < 0.001), diabetes (HR = 0.87, 95% CI: 0.84-0.90, p < 0.001), CKD (HR = 1.08, 95% CI: 1.03-1.13, p = 0.004) and all-cause mortality (HR = 1.18, 95% CI: 1.16-1.20, p < 0.001). However, PGS-RDW was significantly associated with incidence of diabetes (HR = 0.96, 95% CI: 0.94-0.99, p = 0.01), but not with any other tested outcomes.

Discussion

RDW is associated with mortality and incidence of cardiovascular diseases, but a significant association between genetically determined RDW and incident cardiovascular diseases were not observed. However, both RDW and PGS-RDW were inversely associated with incidence of diabetes, suggesting a putative causal relationship. The relationship with incidence of diabetes needs to be further studied.

Structural variation across 138,134 samples in the TOPMed consortium

Ever larger Structural Variant (SV) catalogs highlighting the diversity within and between populations help researchers better understand the links between SVs and disease. The identification of SVs from DNA sequence data is non-trivial and requires a balance between comprehensiveness and precision. Here we present a catalog of 355,667 SVs (59.34% novel) across autosomes and the X chromosome (50bp+) from 138,134 individuals in the diverse TOPMed consortium. We describe our methodologies for SV inference resulting in high variant quality and >90% allele concordance compared to long-read de-novo assemblies of well-characterized control samples. We demonstrate utility through significant associations between SVs and important various cardio-metabolic and hemotologic traits. We have identified 690 SV hotspots and deserts and those that potentially impact the regulation of medically relevant genes. This catalog characterizes SVs across multiple populations and will serve as a valuable tool to understand the impact of SV on disease development and progression.

Association of Kidney Comorbidities and Acute Kidney Failure With Unfavorable Outcomes After COVID-19 in Individuals With the Sickle Cell Trait

Importance

Sickle cell trait (SCT), defined as the presence of 1 hemoglobin beta sickle allele (rs334-T) and 1 normal beta allele, is prevalent in millions of people in the US, particularly in individuals of African and Hispanic ancestry. However, the association of SCT with COVID-19 is unclear.

Objective

To assess the association of SCT with the prepandemic health conditions in participants of the Million Veteran Program (MVP) and to assess the severity and sequelae of COVID-19.

Design, Setting, and Participants

COVID-19 clinical data include 2729 persons with SCT, of whom 353 had COVID-19, and 129 848 SCT-negative individuals, of whom 13 488 had COVID-19. Associations between SCT and COVID-19 outcomes were examined using firth regression. Analyses were performed by ancestry and adjusted for sex, age, age squared, and ancestral principal components to account for population stratification. Data for the study were collected between March 2020 and February 2021.

Exposures

The hemoglobin beta S (HbS) allele (rs334-T).

Main Outcomes and Measures

This study evaluated 4 COVID-19 outcomes derived from the World Health Organization severity scale and phenotypes derived from International Classification of Diseases codes in the electronic health records.

Results

Of the 132 577 MVP participants with COVID-19 data, mean (SD) age at the index date was 64.8 (13.1) years. Sickle cell trait was present in 7.8% of individuals of African ancestry and associated with a history of chronic kidney disease, diabetic kidney disease, hypertensive kidney disease, pulmonary embolism, and cerebrovascular disease. Among the 4 clinical outcomes of COVID-19, SCT was associated with an increased COVID-19 mortality in individuals of African ancestry (n = 3749; odds ratio, 1.77; 95% CI, 1.13 to 2.77; P = .01). In the 60 days following COVID-19, SCT was associated with an increased incidence of acute kidney failure. A counterfactual mediation framework estimated that on average, 20.7% (95% CI, -3.8% to 56.0%) of the total effect of SCT on COVID-19 fatalities was due to acute kidney failure.

Conclusions and Relevance

In this genetic association study, SCT was associated with preexisting kidney comorbidities, increased COVID-19 mortality, and kidney morbidity.

A genome-wide association study in a large community-based cohort identifies multiple loci associated with susceptibility to bacterial and viral infections

There is limited data on host-specific genetic determinants of susceptibility to bacterial and viral infections. Genome-wide association studies using large population cohorts can be a first step towards identifying patients prone to infectious diseases and targets for new therapies. Genetic variants associated with clinically relevant entities of bacterial and viral infections (e.g., abdominal infections, respiratory infections, and sepsis) in 337,484 participants of the UK Biobank cohort were explored by genome-wide association analyses. Cases (n = 81,179) were identified based on ICD-10 diagnosis codes of hospital inpatient and death registries. Functional annotation was performed using gene expression (eQTL) data. Fifty-seven unique genome-wide significant loci were found, many of which are novel in the context of infectious diseases. Some of the detected genetic variants were previously reported associated with infectious, inflammatory, autoimmune, and malignant diseases or key components of the immune system (e.g., white blood cells, cytokines). Fine mapping of the HLA region revealed significant associations with HLA-DQA1, HLA-DRB1, and HLA-DRB4 locus alleles. PPP1R14A showed strong colocalization with abdominal infections and gene expression in sigmoid and transverse colon, suggesting causality. Shared significant loci across infections and non-infectious phenotypes in the UK Biobank cohort were found, suggesting associations for example between SNPs identified for abdominal infections and CRP, rheumatoid arthritis, and diabetes mellitus. We report multiple loci associated with bacterial and viral infections. A better understanding of the genetic determinants of bacterial and viral infections can be useful to identify patients at risk and in the development of new drugs.

Whole genome sequence analysis of platelet traits in the NHLBI Trans-Omics for Precision Medicine (TOPMed) initiative

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.

Bayesian model comparison for rare-variant association studies

Whole-genome sequencing studies applied to large populations or biobanks with extensive phenotyping raise new analytic challenges. The need to consider many variants at a locus or group of genes simultaneously and the potential to study many correlated phenotypes with shared genetic architecture provide opportunities for discovery not addressed by the traditional one variant, one phenotype association study. Here, we introduce a Bayesian model comparison approach called MRP (multiple rare variants and phenotypes) for rare-variant association studies that considers correlation, scale, and direction of genetic effects across a group of genetic variants, phenotypes, and studies, requiring only summary statistic data. We apply our method to exome sequencing data (n = 184,698) across 2,019 traits from the UK Biobank, aggregating signals in genes. MRP demonstrates an ability to recover signals such as associations between PCSK9 and LDL cholesterol levels. We additionally find MRP effective in conducting meta-analyses in exome data. Non-biomarker findings include associations between MC1R and red hair color and skin color, IL17RA and monocyte count, and IQGAP2 and mean platelet volume. Finally, we apply MRP in a multi-phenotype setting; after clustering the 35 biomarker phenotypes based on genetic correlation estimates, we find that joint analysis of these phenotypes results in substantial power gains for gene-trait associations, such as in TNFRSF13B in one of the clusters containing diabetes- and lipid-related traits. Overall, we show that the MRP model comparison approach improves upon useful features from widely used meta-analysis approaches for rare-variant association analyses and prioritizes protective modifiers of disease risk.

An Overview of Strategies for Detecting Genotype-Phenotype Associations Across Ancestrally Diverse Populations

Genome-wide association studies (GWAS) have been very successful at identifying genetic variants influencing a large number of traits. Although the great majority of these studies have been performed in European-descent individuals, it has been recognised that including populations with differing ancestries enhances the potential for identifying causal SNPs due to their differing patterns of linkage disequilibrium. However, when individuals from distinct ethnicities are included in a GWAS, it is necessary to implement a number of control steps to ensure that the identified associations are real genotype-phenotype relationships. In this Review, we discuss the analyses that are required when performing multi-ethnic studies, including methods for determining ancestry at the global and local level for sample exclusion, controlling for ancestry in association testing, and post-GWAS interrogation methods such as genomic control and meta-analysis. We hope that this overview provides a primer for those researchers interested in including distinct populations in their studies.

Multi-ethnic genome-wide association analyses of white blood cell and platelet traits in the Population Architecture using Genomics and Epidemiology (PAGE) study

BACKGROUND

Circulating white blood cell and platelet traits are clinically linked to various disease outcomes and differ across individuals and ancestry groups. Genetic factors play an important role in determining these traits and many loci have been identified. However, most of these findings were identified in populations of European ancestry (EA), with African Americans (AA), Hispanics/Latinos (HL), and other races/ethnicities being severely underrepresented.

RESULTS

We performed ancestry-combined and ancestry-specific genome-wide association studies (GWAS) for white blood cell and platelet traits in the ancestrally diverse Population Architecture using Genomics and Epidemiology (PAGE) Study, including 16,201 AA, 21,347 HL, and 27,236 EA participants. We identified six novel findings at suggestive significance (P < 5E-8), which need confirmation, and independent signals at six previously established regions at genome-wide significance (P < 2E-9). We confirmed multiple previously reported genome-wide significant variants in the single variant association analysis and multiple genes using PrediXcan. Evaluation of loci reported from a Euro-centric GWAS indicated attenuation of effect estimates in AA and HL compared to EA populations.

CONCLUSIONS

Our results highlighted the potential to identify ancestry-specific and ancestry-agnostic variants in participants with diverse backgrounds and advocate for continued efforts in improving inclusion of racially/ethnically diverse populations in genetic association studies for complex traits.

The trans-ancestral genomic architecture of glycemic traits

Glycemic traits are used to diagnose and monitor type 2 diabetes and cardiometabolic health. To date, most genetic studies of glycemic traits have focused on individuals of European ancestry. Here we aggregated genome-wide association studies comprising up to 281,416 individuals without diabetes (30% non-European ancestry) for whom fasting glucose, 2-h glucose after an oral glucose challenge, glycated hemoglobin and fasting insulin data were available. Trans-ancestry and single-ancestry meta-analyses identified 242 loci (99 novel; P < 5 × 10-8), 80% of which had no significant evidence of between-ancestry heterogeneity. Analyses restricted to individuals of European ancestry with equivalent sample size would have led to 24 fewer new loci. Compared with single-ancestry analyses, equivalent-sized trans-ancestry fine-mapping reduced the number of estimated variants in 99% credible sets by a median of 37.5%. Genomic-feature, gene-expression and gene-set analyses revealed distinct biological signatures for each trait, highlighting different underlying biological pathways. Our results increase our understanding of diabetes pathophysiology by using trans-ancestry studies for improved power and resolution.

Whole-genome sequencing association analysis of quantitative red blood cell phenotypes: The NHLBI TOPMed program

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.

An Asian-specific MPL genetic variant alters JAK-STAT signaling and influences platelet count in the population

Genomic discovery efforts for hematological traits have been successfully conducted through genome-wide association study on samples of predominantly European ancestry. We sought to conduct unbiased genetic discovery for coding variants that influence hematological traits in a Han Chinese population. A total of 5257 Han Chinese subjects from Beijing, China were included in the discovery cohort and analyzed by an Illumina ExomeChip array. Replication analyses were conducted in 3827 independent Chinese subjects. We analyzed 12 hematological traits and identified 22 exome-wide significant single-nucleotide polymorphisms (SNP)-trait associations with 15 independent SNPs. Our study provides replication for two associations previously reported but not replicated. Further, one association was identified and replicated in the current study, of a coding variant in the myeloproliferative leukemia (MPL) gene, c.793C > T, p.Leu265Phe (L265F) with increased platelet count (β = 20.6 109 cells/l, Pmeta-analysis = 2.6 × 10-13). This variant is observed at ~2% population frequency in East Asians, whereas it has not been reported in gnomAD European or African populations. Functional analysis demonstrated that expression of MPL L265F in Ba/F3 cells resulted in enhanced phosphorylation of Stat3 and ERK1/2 as compared with the reference MPL allele, supporting altered activation of the JAK-STAT signal transduction pathway as the mechanism underlying the novel association between MPL L265F and platelet count.

The Role, Function, and Mechanism of Long Intergenic Noncoding RNA1184 (linc01184) in Colorectal Cancer

Background

Long intergenic noncoding RNA1184 (linc01184) has been recently discovered; however, its role in human diseases is limited to date. The present study is aimed at investigating the expression pattern and mechanism of linc01184 in colorectal cancer (CRC) tumorigenesis.

Methods

The expression of linc01184 in CRC tissues and cell lines was compared with that in normal controls. The functions of linc01184 in CRC cells were identified by overexpression and small interfering RNA (siRNA) approaches in vitro. Meanwhile, the target gene prediction software, luciferase reporter, RNA pull-down, and western blotting assays were used to analyze the oncogenic mechanism.

Results

We found that linc01184 was obviously upregulated in CRC tissues and cells when compared to normal controls, and its upregulation had a positive association with the CRC progression. linc01184 knockdown significantly suppressed CRC cell proliferation and invasion and promoted apoptosis. Besides, linc01184 acted as a competitive endogenous RNA (ceRNA) by directly binding to microRNA-331 (miR-331), and its overexpression resulted in notable increases of human epidermal growth factor receptor 2 (HER2), phosphorylated Ser/Thr kinases (p-Akt), and extracellular regulated protein kinase 1/2 (p-ERK1/2) at posttranscriptional levels in CRC cells, which were antagonized by miR-331.

Conclusions

The findings reveal for the first time that linc01184 is an enhancer for the proliferation and invasion of CRC by functioning as a ceRNA through the linc01184-miR-331-HER2-p-Akt/ERK1/2 pathway regulatory network.

The genetic structure and adaptation of Andean highlanders and Amazonians are influenced by the interplay between geography and culture

Western South America was one of the worldwide cradles of civilization. The well-known Inca Empire was the tip of the iceberg of an evolutionary process that started 11,000 to 14,000 years ago. Genetic data from 18 Peruvian populations reveal the following: 1) The between-population homogenization of the central southern Andes and its differentiation with respect to Amazonian populations of similar latitudes do not extend northward. Instead, longitudinal gene flow between the northern coast of Peru, Andes, and Amazonia accompanied cultural and socioeconomic interactions revealed by archeology. This pattern recapitulates the environmental and cultural differentiation between the fertile north, where altitudes are lower, and the arid south, where the Andes are higher, acting as a genetic barrier between the sharply different environments of the Andes and Amazonia. 2) The genetic homogenization between the populations of the arid Andes is not only due to migrations during the Inca Empire or the subsequent colonial period. It started at least during the earlier expansion of the Wari Empire (600 to 1,000 years before present). 3) This demographic history allowed for cases of positive natural selection in the high and arid Andes vs. the low Amazon tropical forest: in the Andes, a putative enhancer in HAND2-AS1 (heart and neural crest derivatives expressed 2 antisense RNA1, a noncoding gene related to cardiovascular function) and rs269868-C/Ser1067 in DUOX2 (dual oxidase 2, related to thyroid function and innate immunity) genes and, in the Amazon, the gene encoding for the CD45 protein, essential for antigen recognition by T and B lymphocytes in viral-host interaction.

Population-Matched Transcriptome Prediction Increases TWAS Discovery and Replication Rate

Most genome-wide association studies (GWAS) and transcriptome-wide association studies (TWAS) focus on European populations; however, these results cannot always be accurately applied to non-European populations due to genetic architecture differences. Using GWAS summary statistics in the Population Architecture using Genomics and Epidemiology study, which comprises ∼50,000 Hispanic/Latinos, African Americans, Asians, Native Hawaiians, and Native Americans, we perform TWAS to determine gene-trait associations. We compared results using three transcriptome prediction models derived from Multi-Ethnic Study of Atherosclerosis populations: the African American and Hispanic/Latino (AFHI) model, the European (EUR) model, and the African American, Hispanic/Latino, and European (ALL) model. We identified 240 unique significant trait-associated genes. We found more significant, colocalized genes that replicate in larger cohorts when applying the AFHI model than the EUR or ALL model. Thus, TWAS with population-matched transcriptome models have more power for discovery and replication, demonstrating the need for more transcriptome studies in diverse populations.

•

TWAS mechanistically extends GWAS findings in diverse populations

•

Population-matched transcriptome models detect more replicable associations

•

Colocalization shows GWAS variants likely act through gene expression regulation

•

More GWAS and transcriptome modeling in diverse populations are needed

Ancestry-specific associations identified in genome-wide combined-phenotype study of red blood cell traits emphasize benefits of diversity in genomics

BACKGROUND

Quantitative red blood cell (RBC) traits are highly polygenic clinically relevant traits, with approximately 500 reported GWAS loci. The majority of RBC trait GWAS have been performed in European- or East Asian-ancestry populations, despite evidence that rare or ancestry-specific variation contributes substantially to RBC trait heritability. Recently developed combined-phenotype methods which leverage genetic trait correlation to improve statistical power have not yet been applied to these traits. Here we leveraged correlation of seven quantitative RBC traits in performing a combined-phenotype analysis in a multi-ethnic study population.

RESULTS

We used the adaptive sum of powered scores (aSPU) test to assess combined-phenotype associations between ~ 21 million SNPs and seven RBC traits in a multi-ethnic population (maximum n = 67,885 participants; 24% African American, 30% Hispanic/Latino, and 43% European American; 76% female). Thirty-nine loci in our multi-ethnic population contained at least one significant association signal (p < 5E-9), with lead SNPs at nine loci significantly associated with three or more RBC traits. A majority of the lead SNPs were common (MAF > 5%) across all ancestral populations. Nineteen additional independent association signals were identified at seven known loci (HFE, KIT, HBS1L/MYB, CITED2/FILNC1, ABO, HBA1/2, and PLIN4/5). For example, the HBA1/2 locus contained 14 conditionally independent association signals, 11 of which were previously unreported and are specific to African and Amerindian ancestries. One variant in this region was common in all ancestries, but exhibited a narrower LD block in African Americans than European Americans or Hispanics/Latinos. GTEx eQTL analysis of all independent lead SNPs yielded 31 significant associations in relevant tissues, over half of which were not at the gene immediately proximal to the lead SNP.

CONCLUSION

This work identified seven loci containing multiple independent association signals for RBC traits using a combined-phenotype approach, which may improve discovery in genetically correlated traits. Highly complex genetic architecture at the HBA1/2 locus was only revealed by the inclusion of African Americans and Hispanics/Latinos, underscoring the continued importance of expanding large GWAS to include ancestrally diverse populations.

Mosaic chromosome Y loss is associated with alterations in blood cell counts in UK Biobank men

Mosaic loss of Y chromosome (mLOY) is the most frequently detected somatic copy number alteration in leukocytes of men. In this study, we investigate blood cell counts as a potential mechanism linking mLOY to disease risk in 206,353 UK males. Associations between mLOY, detected by genotyping arrays, and blood cell counts were assessed by multivariable linear models adjusted for relevant risk factors. Among the participants, mLOY was detected in 39,809 men. We observed associations between mLOY and reduced erythrocyte count (−0.009 [−0.014, −0.005] × 10¹² cells/L, p = 2.75 × 10⁻⁵) and elevated thrombocyte count (5.523 [4.862, 6.183] × 10⁹ cells/L, p = 2.32 × 10⁻⁶⁰) and leukocyte count (0.218 [0.198, 0.239] × 10⁹ cells/L, p = 9.22 × 10⁻⁹⁵), particularly for neutrophil count (0.174 × [0.158, 0.190]10⁹ cells/L, p = 1.24 × 10⁻⁹⁹) and monocyte count (0.021 [0.018 to 0.024] × 10⁹ cells/L, p = 6.93 × 10⁻⁵⁷), but lymphocyte count was less consistent (0.016 [0.007, 0.025] × 10⁹ cells/L, p = 8.52 × 10⁻⁴). Stratified analyses indicate these associations are independent of the effects of aging and smoking. Our findings provide population-based evidence for associations between mLOY and blood cell counts that should stimulate investigation of the underlying biological mechanisms linking mLOY to cancer and chronic disease risk.

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

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.

Impact of Rare and Common Genetic Variants on Diabetes Diagnosis by Hemoglobin A1c in Multi-Ancestry Cohorts: The Trans-Omics for Precision Medicine Program

Hemoglobin A1c (HbA1c) is widely used to diagnose diabetes and assess glycemic control in individuals with diabetes. However, nonglycemic determinants, including genetic variation, may influence how accurately HbA1c reflects underlying glycemia. Analyzing the NHLBI Trans-Omics for Precision Medicine (TOPMed) sequence data in 10,338 individuals from five studies and four ancestries (6,158 Europeans, 3,123 African-Americans, 650 Hispanics, and 407 East Asians), we confirmed five regions associated with HbA1c (GCK in Europeans and African-Americans, HK1 in Europeans and Hispanics, FN3K and/or FN3KRP in Europeans, and G6PD in African-Americans and Hispanics) and we identified an African-ancestry-specific low-frequency variant (rs1039215 in HBG2 and HBE1, minor allele frequency (MAF) = 0.03). The most associated G6PD variant (rs1050828-T, p.Val98Met, MAF = 12% in African-Americans, MAF = 2% in Hispanics) lowered HbA1c (-0.88% in hemizygous males, -0.34% in heterozygous females) and explained 23% of HbA1c variance in African-Americans and 4% in Hispanics. Additionally, we identified a rare distinct G6PD coding variant (rs76723693, p.Leu353Pro, MAF = 0.5%; -0.98% in hemizygous males, -0.46% in heterozygous females) and detected significant association with HbA1c when aggregating rare missense variants in G6PD. We observed similar magnitude and direction of effects for rs1039215 (HBG2) and rs76723693 (G6PD) in the two largest TOPMed African American cohorts, and we replicated the rs76723693 association in the UK Biobank African-ancestry participants. These variants in G6PD and HBG2 were monomorphic in the European and Asian samples. African or Hispanic ancestry individuals carrying G6PD variants may be underdiagnosed for diabetes when screened with HbA1c. Thus, assessment of these variants should be considered for incorporation into precision medicine approaches for diabetes diagnosis.

A Genome-Wide Association Study Identifies Blood Disorder-Related Variants Influencing Hemoglobin A1c With Implications for Glycemic Status in U.S. Hispanics/Latinos

OBJECTIVE

We aimed to identify hemoglobin A1c (HbA1c)-associated genetic variants and examine their implications for glycemic status evaluated by HbA1c in U.S. Hispanics/Latinos with diverse genetic ancestries.

RESEARCH DESIGN AND METHODS

We conducted a genome-wide association study (GWAS) of HbA1c in 9,636 U.S. Hispanics/Latinos without diabetes from the Hispanic Community Health Study/Study of Latinos, followed by a replication among 4,729 U.S. Hispanics/Latinos from three independent studies.

RESULTS

Our GWAS and replication analyses showed 10 previously known and novel loci associated with HbA1c at genome-wide significance levels (P < 5.0 × 10-8). In particular, two African ancestry-specific variants, HBB-rs334 and G6PD-rs1050828, which are causal mutations for sickle cell disease and G6PD deficiency, respectively, had ∼10 times larger effect sizes on HbA1c levels (β = -0.31% [-3.4 mmol/mol]) and -0.35% [-3.8 mmol/mol] per minor allele, respectively) compared with other HbA1c-associated variants (0.03-0.04% [0.3-0.4 mmol/mol] per allele). A novel Amerindian ancestry-specific variant, HBM-rs145546625, was associated with HbA1c and hematologic traits but not with fasting glucose. The prevalence of hyperglycemia (prediabetes and diabetes) defined using fasting glucose or oral glucose tolerance test 2-h glucose was similar between carriers of HBB-rs334 or G6PD-rs1050828 HbA1c-lowering alleles and noncarriers, whereas the prevalence of hyperglycemia defined using HbA1c was significantly lower in carriers than in noncarriers (12.2% vs. 28.4%, P < 0.001). After recalibration of the HbA1c level taking HBB-rs334 and G6PD-rs1050828 into account, the prevalence of hyperglycemia in carriers was similar to noncarriers (31.3% vs. 28.4%, P = 0.28).

CONCLUSIONS

This study in U.S. Hispanics/Latinos found several ancestry-specific alleles associated with HbA1c through erythrocyte-related rather than glycemic-related pathways. The potential influences of these nonglycemic-related variants need to be considered when the HbA1c test is performed.

Genomics of disease risk in globally diverse populations

Risk of disease is multifactorial and can be shaped by socio-economic, demographic, cultural, environmental and genetic factors. Our understanding of the genetic determinants of disease risk has greatly advanced with the advent of genome-wide association studies (GWAS), which detect associations between genetic variants and complex traits or diseases by comparing populations of cases and controls. However, much of this discovery has occurred through GWAS of individuals of European ancestry, with limited representation of other populations, including from Africa, The Americas, Asia and Oceania. Population demography, genetic drift and adaptation to environments over thousands of years have led globally to the diversification of populations. This global genomic diversity can provide new opportunities for discovery and translation into therapies, as well as a better understanding of population disease risk. Large-scale multi-ethnic and representative biobanks and population health resources provide unprecedented opportunities to understand the genetic determinants of disease on a global scale.

Unbiased Phenome-Wide Association Studies of Red Cell Distribution Width Identifies Key Associations with Pulmonary Hypertension

Rationale: Red cell distribution width (RDW) is a prognostic factor in many diseases; however, its clinical utility remains limited because the relative value of RDW as a biomarker across disease states has not been established. Objectives: To establish an unbiased RDW disease hierarchy to guide the clinical use of RDW and to assess its relationship to cardiovascular hemodynamic and structural parameters. Methods: We performed phenome-wide association studies for RDW in discovery and replication cohorts derived from a deidentified electronic health record in nonanemic individuals. RDW values obtained within 30 days of echocardiogram or right heart catheterization were tested for association with structural and hemodynamic variables. Results: RDW was associated with 263 phenotypes in both men and women in the discovery cohort (n = 121,530), 48 of which replicated in an independent cohort (n = 2,039). The strongest associations were observed with pulmonary arterial hypertension (odds ratio [OR], 2.1; 95% confidence interval [CI], 1.9-2.3), chronic pulmonary heart disease (OR, 2.0; 95% CI, 1.9-2.2), and congestive heart failure (OR, 1.9; 95% CI, 1.8-2.0); P < 1 × 10-74 for all. By echocardiography, RDW was higher in the setting of right ventricular dysfunction than left ventricular dysfunction (P < 0.001). Measured invasively, mean pulmonary arterial pressure was associated with RDW (21 vs. 33 mm Hg at 25th vs. 75th percentile RDW; P < 1 × 10-7) and remained strongly significant even when controlling for mean pulmonary capillary wedge pressure (21 vs. 29 mm Hg at 25th vs. 75th percentile RDW; P < 1 × 10-7). Conclusions: Among 1,364 coded medical conditions, increased RDW was strongly associated with pulmonary hypertension and heart failure. Hemodynamic and echocardiographic phenotyping confirmed these associations and underscored that the most clinically relevant phenotype associated with RDW was pulmonary hypertension. These hypothesis-generating findings highlight the potential shared pathophysiology of pulmonary hypertension and elevated RDW. Elevated RDW in the absence of anemia should alert clinicians to the potential for underlying cardiopulmonary disease.

Complimentary Methods for Multivariate Genome-Wide Association Study Identify New Susceptibility Genes for Blood Cell Traits

Genome-wide association studies (GWAS) have found hundreds of novel loci associated with full blood count (FBC) phenotypes. However, most of these studies were performed in a single phenotype framework without putting into consideration the clinical relatedness among traits. In this work, in addition to the standard univariate GWAS, we also use two different multivariate methods to perform the first multiple traits GWAS of FBC traits in ∼7000 individuals from the Ugandan General Population Cohort (GPC). We started by performing the standard univariate GWAS approach. We then performed our first multivariate method, in this approach, we tested for marker associations with 15 FBC traits simultaneously in a multivariate mixed model implemented in GEMMA while accounting for the relatedness of individuals and pedigree structures, as well as population substructure. In this analysis, we provide a framework for the combination of multiple phenotypes in multivariate GWAS analysis and show evidence of multi-collinearity whenever the correlation between traits exceeds the correlation coefficient threshold of r ² >=0.75. This approach identifies two known and one novel loci. In the second multivariate method, we applied principal component analysis (PCA) to the same 15 correlated FBC traits. We then tested for marker associations with each PC in univariate linear mixed models implemented in GEMMA. We show that the FBC composite phenotype as assessed by each PC expresses information that is not completely encapsulated by the individual FBC traits, as this approach identifies three known and five novel loci that were not identified using both the standard univariate and multivariate GWAS methods. Across both multivariate methods, we identified six novel loci. As a proof of concept, both multivariate methods also identified known loci, HBB and ITFG3. The two multivariate methods show that multivariate genotype-phenotype methods increase power and identify novel genotype-phenotype associations not found with the standard univariate GWAS in the same dataset.

Varicose veins of lower extremities: Insights from the first large-scale genetic study

Varicose veins of lower extremities (VVs) are a common multifactorial vascular disease. Genetic factors underlying VVs development remain largely unknown. Here we report the first large-scale study of VVs performed on a freely available genetic data of 408,455 European-ancestry individuals. We identified the 12 reliably associated loci that explain 13% of the SNP-based heritability, and prioritized the most likely causal genes CASZ1, PIEZO1, PPP3R1, EBF1, STIM2, HFE, GATA2, NFATC2, and SOX9. VVs-associated variants within these loci exhibited pleiotropic effects on several phenotypes including blood pressure/hypertension and blood cell traits. Gene set enrichment analysis revealed gene categories related to abnormal vasculogenesis. Genetic correlation analysis confirmed known epidemiological associations between VVs and deep venous thrombosis, weight, rough labor, and standing job, and found a genetic overlap with multiple traits that have not been previously suspected to share common genetic background with VVs. These traits included educational attainment, fluid intelligence and prospective memory scores, walking pace (negative correlation with VVs), smoking, height, number of operations, pain, and gonarthrosis (positive correlation with VVs). Finally, Mendelian randomization analysis provided evidence for causal effects of plasma levels of MICB and CD209 proteins, and anthropometric traits such as waist and hip circumference, height, weight, and both fat and fat-free mass. Our results provide novel insight into both VVs genetics and etiology. The revealed genes and proteins can be considered as good candidates for follow-up functional studies and might be of interest as potential drug targets.

Varicose veins of lower extremities (VVs) affect about 30% of adults in developed countries and cause both cosmetic and health problems. A strong body of evidence indicates that heredity plays an important role in the etiology of this condition. However, genetic basis of VVs remains poorly understood. Here, we present the results of the first large-scale genetic study for VVs. We identified genes which are the most likely involved in VVs pathogenesis. We show that VVs are correlated at a genetic level with numerous traits and phenotypes, including those already known from prior epidemiological studies (deep venous thrombosis, body mass index, standing job, etc.) as well as with those that have not been suspected to share common genetic background with VVs (fluid intelligence and prospective memory scores, smoking, walking pace, pain all over the body, and other traits). Finally, using genetic variants as instruments, we demonstrate direct causal effects of the traits related to anthropometry, such as height and weight, and plasma levels of immune-related proteins MICB and CD209. Our study provides novel insight into both VVs genetics and etiology. The revealed genes (CASZ1, PIEZO1, PPP3R1, EBF1, STIM2, HFE, GATA2, NFATC2, and SOX9) and proteins (MICB and CD209) can be considered as good candidates for follow-up functional studies and might be of interest as potential drug targets.

The Future of Genomic Studies Must Be Globally Representative: Perspectives from PAGE

The past decade has seen a technological revolution in human genetics that has empowered population-level investigations into genetic associations with phenotypes. Although these discoveries rely on genetic variation across individuals, association studies have overwhelmingly been performed in populations of European descent. In this review, we describe limitations faced by single-population studies and provide an overview of strategies to improve global representation in existing data sets and future human genomics research via diversity-focused, multiethnic studies. We highlight the successes of individual studies and meta-analysis consortia that have provided unique knowledge. Additionally, we outline the approach taken by the Population Architecture Using Genomics and Epidemiology (PAGE) study to develop best practices for performing genetic epidemiology in multiethnic contexts. Finally, we discuss how limiting investigations to single populations impairs findings in the clinical domain for both rare-variant identification and genetic risk prediction.

Genetic signatures of gene flow and malaria-driven natural selection in sub-Saharan populations of the "endemic Burkitt Lymphoma belt"

Populations in sub-Saharan Africa have historically been exposed to intense selection from chronic infection with falciparum malaria. Interestingly, populations with the highest malaria intensity can be identified by the increased occurrence of endemic Burkitt Lymphoma (eBL), a pediatric cancer that affects populations with intense malaria exposure, in the so called "eBL belt" in sub-Saharan Africa. However, the effects of intense malaria exposure and sub-Saharan populations' genetic histories remain poorly explored. To determine if historical migrations and intense malaria exposure have shaped the genetic composition of the eBL belt populations, we genotyped ~4.3 million SNPs in 1,708 individuals from Ghana and Northern Uganda, located on opposite sides of eBL belt and with ≥ 7 months/year of intense malaria exposure and published evidence of high incidence of BL. Among 35 Ghanaian tribes, we showed a predominantly West-Central African ancestry and genomic footprints of gene flow from Gambian and East African populations. In Uganda, the North West population showed a predominantly Nilotic ancestry, and the North Central population was a mixture of Nilotic and Southern Bantu ancestry, while the Southwest Ugandan population showed a predominant Southern Bantu ancestry. Our results support the hypothesis of diverse ancestral origins of the Ugandan, Kenyan and Tanzanian Great Lakes African populations, reflecting a confluence of Nilotic, Cushitic and Bantu migrations in the last 3000 years. Natural selection analyses suggest, for the first time, a strong positive selection signal in the ATP2B4 gene (rs10900588) in Northern Ugandan populations. These findings provide important baseline genomic data to facilitate disease association studies, including of eBL, in eBL belt populations.

Generalizing polygenic risk scores from Europeans to Hispanics/Latinos

Polygenic risk scores (PRSs) are weighted sums of risk allele counts of single-nucleotide polymorphisms (SNPs) associated with a disease or trait. PRSs are typically constructed based on published results from Genome-Wide Association Studies (GWASs), and the majority of which has been performed in large populations of European ancestry (EA) individuals. Although many genotype-trait associations have generalized across populations, the optimal choice of SNPs and weights for PRSs may differ between populations due to different linkage disequilibrium (LD) and allele frequency patterns. We compare various approaches for PRS construction, using GWAS results from both large EA studies and a smaller study in Hispanics/Latinos: The Hispanic Community Health Study/Study of Latinos (HCHS/SOL, n = 12 , 803 ). We consider multiple approaches for selecting SNPs and for computing SNP weights. We study the performance of the resulting PRSs in an independent study of Hispanics/Latinos from the Women's Health Initiative (WHI, n = 3 , 582 ). We support our investigation with simulation studies of potential genetic architectures in a single locus. We observed that selecting variants based on EA GWASs generally performs well, except for blood pressure trait. However, the use of EA GWASs for weight estimation was suboptimal. Using non-EA GWAS results to estimate weights improved results.

Detecting past and ongoing natural selection among ethnically Tibetan women at high altitude in Nepal

Adaptive evolution in humans has rarely been characterized for its whole set of components, i.e. selective pressure, adaptive phenotype, beneficial alleles and realized fitness differential. We combined approaches for detecting polygenic adaptations and for mapping the genetic bases of physiological and fertility phenotypes in approximately 1000 indigenous ethnically Tibetan women from Nepal, adapted to high altitude. The results of genome-wide association analyses and tests for polygenic adaptations showed evidence of positive selection for alleles associated with more pregnancies and live births and evidence of negative selection for those associated with higher offspring mortality. Lower hemoglobin level did not show clear evidence for polygenic adaptation, despite its strong association with an EPAS1 haplotype carrying selective sweep signals.

Generalization and fine mapping of red blood cell trait genetic associations to multi-ethnic populations: The PAGE Study

Red blood cell (RBC) traits provide insight into a wide range of physiological states and exhibit moderate to high heritability, making them excellent candidates for genetic studies to inform underlying biologic mechanisms. Previous RBC trait genome-wide association studies were performed primarily in European- or Asian-ancestry populations, missing opportunities to inform understanding of RBC genetic architecture in diverse populations and reduce intervals surrounding putative functional SNPs through fine-mapping. Here, we report the first fine-mapping of six correlated (Pearson's r range: |0.04 - 0.92|) RBC traits in up to 19,036 African Americans and 19,562 Hispanic/Latinos participants of the Population Architecture using Genomics and Epidemiology (PAGE) consortium. Trans-ethnic meta-analysis of race/ethnic- and study-specific estimates for approximately 11,000 SNPs flanking 13 previously identified association signals as well as 150,000 additional array-wide SNPs was performed using inverse-variance meta-analysis after adjusting for study and clinical covariates. Approximately half of previously reported index SNP-RBC trait associations generalized to the trans-ethnic study population (p<1.7x10-4 ); previously unreported independent association signals within the ABO region reinforce the potential for multiple functional variants affecting the same locus. Trans-ethnic fine-mapping did not reveal additional signals at the HFE locus independent of the known functional variants. Finally, we identified a potential novel association in the Hispanic/Latino study population at the HECTD4/RPL6 locus for RBC count (p=1.9x10-7 ). The identification of a previously unknown association, generalization of a large proportion of known association signals, and refinement of known association signals all exemplify the benefits of genetic studies in diverse populations. This article is protected by copyright. All rights reserved.

PheWAS and Beyond: The Landscape of Associations with Medical Diagnoses and Clinical Measures across 38,662 Individuals from Geisinger

Most phenome-wide association studies (PheWASs) to date have used a small to moderate number of SNPs for association with phenotypic data. We performed a large-scale single-cohort PheWAS, using electronic health record (EHR)-derived case-control status for 541 diagnoses using International Classification of Disease version 9 (ICD-9) codes and 25 median clinical laboratory measures. We calculated associations between these diagnoses and traits with ∼630,000 common frequency SNPs with minor allele frequency > 0.01 for 38,662 individuals. In this landscape PheWAS, we explored results within diseases and traits, comparing results to those previously reported in genome-wide association studies (GWASs), as well as previously published PheWASs. We further leveraged the context of functional impact from protein-coding to regulatory regions, providing a deeper interpretation of these associations. The comprehensive nature of this PheWAS allows for novel hypothesis generation, the identification of phenotypes for further study for future phenotypic algorithm development, and identification of cross-phenotype associations.

Common α-globin variants modify hematologic and other clinical phenotypes in sickle cell trait and disease

Co-inheritance of α-thalassemia has a significant protective effect on the severity of complications of sickle cell disease (SCD), including stroke. However, little information exists on the association and interactions for the common African ancestral α-thalassemia mutation (-α3.7 deletion) and β-globin traits (HbS trait [SCT] and HbC trait) on important clinical phenotypes such as red blood cell parameters, anemia, and chronic kidney disease (CKD). In a community-based cohort of 2,916 African Americans from the Jackson Heart Study, we confirmed the expected associations between SCT, HbC trait, and the -α3.7 deletion with lower mean corpuscular volume/mean corpuscular hemoglobin and higher red blood cell count and red cell distribution width. In addition to the recently recognized association of SCT with lower estimated glomerular filtration rate and glycated hemoglobin (HbA1c), we observed a novel association of the -α3.7 deletion with higher HbA1c levels. Co-inheritance of each additional copy of the -α3.7 deletion significantly lowered the risk of anemia and chronic kidney disease among individuals with SCT (P-interaction = 0.031 and 0.019, respectively). Furthermore, co-inheritance of a novel α-globin regulatory variant was associated with normalization of red cell parameters in individuals with the -α3.7 deletion and significantly negated the protective effect of α-thalassemia on stroke in 1,139 patients with sickle cell anemia from the Cooperative Study of Sickle Cell Disease (CSSCD) (P-interaction = 0.0049). Functional assays determined that rs11865131, located in the major alpha-globin enhancer MCS-R2, was the most likely causal variant. These findings suggest that common α- and β-globin variants interact to influence hematologic and clinical phenotypes in African Americans, with potential implications for risk-stratification and counseling of individuals with SCD and SCT.

A standardized framework for representation of ancestry data in genomics studies, with application to the NHGRI-EBI GWAS Catalog

The accurate description of ancestry is essential to interpret, access, and integrate human genomics data, and to ensure that these benefit individuals from all ancestral backgrounds. However, there are no established guidelines for the representation of ancestry information. Here we describe a framework for the accurate and standardized description of sample ancestry, and validate it by application to the NHGRI-EBI GWAS Catalog. We confirm known biases and gaps in diversity, and find that African and Hispanic or Latin American ancestry populations contribute a disproportionately high number of associations. It is our hope that widespread adoption of this framework will lead to improved analysis, interpretation, and integration of human genomics data.

Novel mechanisms of PIEZO1 dysfunction in hereditary xerocytosis

Mutations in PIEZO1 are the primary cause of hereditary xerocytosis, a clinically heterogeneous, dominantly inherited disorder of erythrocyte dehydration. We used next-generation sequencing-based techniques to identify PIEZO1 mutations in individuals from 9 kindreds referred with suspected hereditary xerocytosis (HX) and/or undiagnosed congenital hemolytic anemia. Mutations were primarily found in the highly conserved, COOH-terminal pore-region domain. Several mutations were novel and demonstrated ethnic specificity. We characterized these mutations using genomic-, bioinformatic-, cell biology-, and physiology-based functional assays. For these studies, we created a novel, cell-based in vivo system for study of wild-type and variant PIEZO1 membrane protein expression, trafficking, and electrophysiology in a rigorous manner. Previous reports have indicated HX-associated PIEZO1 variants exhibit a partial gain-of-function phenotype with generation of mechanically activated currents that inactivate more slowly than wild type, indicating that increased cation permeability may lead to dehydration of PIEZO1-mutant HX erythrocytes. In addition to delayed channel inactivation, we found additional alterations in mutant PIEZO1 channel kinetics, differences in response to osmotic stress, and altered membrane protein trafficking, predicting variant alleles that worsen or ameliorate erythrocyte hydration. These results extend the genetic heterogeneity observed in HX and indicate that various pathophysiologic mechanisms contribute to the HX phenotype.