Abstract
Background α-Globin is expressed in endothelial cells of resistance arteries, where it limits endothelial nitric oxide signaling and enhances α-adrenergic–mediated vasoconstriction. α-Globin gene (HBA) copy number is variable in people of African descent and other populations worldwide. Given the protective effect of nitric oxide in the kidney, we hypothesized that HBA copy number would be associated with kidney disease risk.
Methods Community-dwelling Black Americans aged ≥45 years old were enrolled in a national longitudinal cohort from 2003 through 2007. HBA copy number was measured using droplet digital PCR. The prevalence ratio (PR) of CKD and the relative risk (RR) of incident reduced eGFR were calculated using modified Poisson multivariable regression. The hazard ratio (HR) of incident ESKD was calculated using Cox proportional hazards multivariable regression.
Results Among 9908 participants, HBA copy number varied from 2 to 6. In analyses adjusted for demographic, clinical, and genetic risk factors, a one-copy increase in HBA was associated with 14% greater prevalence of CKD (PR, 1.14; 95% CI, 1.07 to 1.21; P<0.0001). While HBA copy number was not associated with incident reduced eGFR (RR, 1.06; 95% CI, 0.94 to 1.19; P=0.38), the hazard of incident ESKD was 32% higher for each additional copy of HBA (HR, 1.32; 95% CI, 1.09 to 1.61; P=0.005).
Conclusions Increasing HBA copy number was associated with a greater prevalence of CKD and incidence of ESKD in a national longitudinal cohort of Black Americans.
- human genetics
- nitric oxide
- hemodynamics and vascular regulation
- alpha globin
- alpha thalassemia
- kidney disease
- African Americans
- prevalence
- gene dosage
Black Americans develop kidney disease at a younger age than other Americans and are three times more likely to develop end-stage kidney disease (ESKD), even after accounting for socioeconomic factors and comorbid medical conditions.1⇓⇓⇓⇓–6 DNA sequence variants that increase the risk of kidney disease (including one in the hemoglobin β locus [HBB]), sickle cell trait,7,8 and apolipoprotein L-1 (APOL1)9⇓–11 are more common among Black Americans, yet only partly explain the racial disparity in kidney disease. The evaluation of genetic risk factors for diseases common in minority populations in the United States remains a high priority. However, Black Americans and other minority populations are underrepresented in many large-scale longitudinal US population studies, hindering genetic investigations into kidney disease and other common chronic diseases.12,13
Basic science14⇓⇓⇓⇓⇓–20 and clinical studies21⇓–23 have identified nitric oxide signaling as important in protection against and recovery from ischemic or oxidative injury to the kidney; however, to date there have been few clear associations between genetic variants in nitric oxide signaling pathways in the kidney and kidney function or kidney disease.24⇓⇓–27 α-Globin has been proposed to function as a regulator of nitric oxide signaling in the endothelial cells of small resistance arteries.28 Genetic or pharmacologic disruption of α-globin or its stabilizing protein leads to enhanced nitric oxide signaling and decreased vasoconstriction in response to α-adrenergic stimuli,29⇓–31 potentially conferring protection against kidney injury.
α-Globin is expressed by the HBA1 and HBA2 genes that are organized in tandem on chromosome 16. The sequence similarity of the HBA1 and HBA2 genes makes them susceptible to deletions, which are common among individuals from Africa, Southeast Asia, the Pacific Islands, and the Mediterranean.32 The 3.7-kilobase insertion/deletion variant is common among Black Americans. Because of the unique challenges of genotyping this structural variant, it has been omitted from genome-wide association studies. This deletion spans parts of the HBA2 and HBA1 genes and reduces the functional gene copy number by one, leading to a decrease in total α-globin gene (HBA) expression in red cell precursors.33 The quantitative relationship between HBA copy number and α-globin expression in red blood cell (RBC) precursors is the molecular basis for the spectrum of α-thalassemia syndromes.33⇓–35 In light of the newly recognized function of α-globin as a restrictor of endothelial nitric oxide signaling,28 we hypothesized that α-globin gene deletions would be associated with protection against prevalent CKD, incident reduced eGFR, and incident ESKD among Black Americans.
Methods
Study Design
The REasons for Geographic and Racial Differences in Stroke (REGARDS) study is a longitudinal cohort study designed to determine the reasons for racial disparities in stroke and cognitive decline in Black and White Americans aged ≥45 years.36 REGARDS enrolled 30,239 community-dwelling participants from the 48 contiguous United States from 2003 to 2007. A total of 12,514 (41%) participants were Black, and 56% were from states considered to be in the stroke belt, where stroke incidence in the United States is highest. Exclusion criteria included self-reported race other than Black or White, residence in or on the waiting list for a nursing home, active cancer within the past year, or inability to communicate in English.
Baseline variables were collected via standardized computer-assisted telephone interview, self-administered questionnaire, and in-home physical examination.36,37 Trained personnel administered computer-assisted telephone interviews to collect study participant age, sex, region of residence, insurance status, education level, income, and reports of physician-diagnosed comorbid conditions. Trained personnel conducted in-home examinations to measure height, weight, and BP; collect blood and urine specimens; and review medication bottles. Baseline measurements were repeated in some participants during a second in-home visit conducted a median (25th, 75th percentile) of 9.5 (8.7, 9.9) years after baseline.
All participants provided oral and written informed consent. The REGARDS study was approved by participating center Institutional Review Boards. All self-reported Black participants consenting to genetic research were eligible to participate in this study. Participants were excluded if they had hemoglobin SS (n=5), SC (n=1), or CC (n=6); or had ESKD at enrollment (n=91) (Figure 1). This study followed the Strengthening the Reporting of Observational Studies in Epidemiology reporting guideline.
Cohort flow diagram.
Main Exposure Variable
HBA copy number was evaluated as a numeric variable with values of 2, 3, 4, 5, or 6 as determined by droplet digital PCR (ddPCR) analysis of genomic DNA. The ddPCR copy number assay targeted a unique sequence within the 3.7-kb insertion/deletion polymorphism. We determined copy number of this target relative to a reference gene, EIF2C1 (human eukaryotic translation initiation factor 2C, 1). Two-dimensional clusters of droplet counts for target and reference genes were manually gated using Quantasoft (Bio-Rad) per the manufacturer's protocols. Droplet counts, copy number variant values, and 95% CIs for copy number variant were extracted and visualized, and genotype was assigned using custom scripts in the R computing environment without user intervention (see the Supplemental Material).
Outcome Measures
CKD prevalence was defined by an eGFR<60 ml/min per 1.73 m2 or a urine albumin-creatinine ratio ≥30 mg/g or at baseline. The Chronic Kidney Disease Epidemiology Collaboration equation was used to estimate GFR from serum creatinine concentrations. In addition, prevalence of eGFR<60 ml/min per 1.73 m2 and a urine albumin-creatinine ratio ≥30 mg/g at baseline were evaluated as separate outcome measures. Incident reduced eGFR was defined by an eGFR<60 ml/min per 1.73 m2 at the follow-up in-home visit and >40% decline in eGFR from baseline, among those who had eGFR≥60 ml/min per 1.73 m2 at baseline.38 This definition of incident reduced eGFR emphasized the clinically relevant transition across the 60-ml/min per 1.73 m2 threshold and did not examine changes among those who were below that threshold at baseline. The incident reduced eGFR analysis was limited to those participants with available follow-up visit data (Figure 1). Incident ESKD was identified by linkage to United States Renal Data System (USRDS) data through December 31, 2018.
Covariates
Hypertension was defined as one or more of the following: SBP≥140 mm Hg or DBP≥90 mm Hg (recorded as the average of two measurements performed according to a standardized protocol39); self-reported current use of medication to control BP; or two or more antihypertensive medications found on medication bottle review during the baseline in-home study visit. A clinical definition of hypertension was employed as a covariate instead of measured BP because measured BP would be variably affected by antihypertensive use. Hemoglobin, RBC count, and mean cell volume (MCV) values were measured with an automated hematology analyzer and mean cell hemoglobin (MCH), mean corpuscular hemoglobin concentration (MCHC), and red cell distribution width–coefficient of variation were calculated using standard formulae.40 Diabetes mellitus was defined by a fasting glucose level ≥126 mg/dl, a random glucose ≥200 mg/dl, or use of antidiabetic medication. Body mass index was calculated from measured height and weight. Smoking status was self-reported and categorized as never, past, or current smoker. Region was defined as three geographic areas: stroke belt buckle (coastal plains of North Carolina, South Carolina, and Georgia), stroke belt (the rest of North Carolina, South Carolina, and Georgia and the entire states of Tennessee, Alabama, Mississippi, Arkansas, and Louisiana), and stroke non-belt (the remaining contiguous United States).36,41 Age, sex, race, health insurance (yes or no), highest education level obtained (less than high school, high school, some college, college or more), and annual income (≤$20K, $20–34K, $35–74K, ≥$75K) were self-reported. Sickle cell trait and hemoglobin C genotypes were measured by a TaqMan single nucleotide polymorphism (SNP) genotyping assay (Applied Biosystems/Thermo Fisher Scientific).8 APOL1 high-risk genotype was defined as the presence of two renal risk alleles compared with a reference of less than two high-risk alleles (Supplemental Material).42,43
Statistical Analysis
The analytic plan and outcome measure definitions were prespecified and approved by the REGARDS Cohort Study Executive Committee in December of 2018.
Unadjusted Analyses of Participant Characteristics by HBA Copy Number
For continuous measures of participant characteristics, medians and 25th and 75th percentiles were calculated by HBA copy number. Differences between copy number groups were assessed by Kruskal–Wallis nonparametric ANOVA. Categoric variables were calculated as percentages within each category and differences were assessed by chi-squared tests of association.
Adjusted Analyses of HBA Copy Number and Kidney Disease Outcome Measures
The primary outcome measure was the presence or absence of CKD at baseline (CKD prevalence) and the main exposure variable was HBA copy number. In addition, secondary outcome measures for incident reduced eGFR and incident ESKD were evaluated. Thirteen known risk factors for kidney disease were included in multivariable regression models for CKD prevalence, incident reduced eGFR, and incident ESKD (see ‘Covariates’ section above). A modified Poisson regression model with robust variance estimation was used to estimate the prevalence ratio (PR) and 95% confidence interval (95% CI) for CKD, the separate outcomes of eGFR<60 ml/min per 1.73 m2 or urine albumin-creatinine ratio ≥30 mg/g at baseline, and the relative risk (RR) and 95% CI for incident reduced eGFR. These models estimated the difference in log-transformed PR or RR for each additional HBA allele. The use of Poisson regression with robust variance estimation allows for estimating the RR directly (reported here as the PR for prevalent CKD and RR for incident reduced eGFR), rather than the odds ratio that arises from logistic regression. The RR is regarded as a more clinically meaningful and transparent assessment of risk than the odds ratio in cohort studies.44 Incident ESKD hazard ratio (HR) and 95% CI were estimated with a multivariable Cox proportional hazards model for the time between baseline and date of ESKD onset from the USRDS. Baseline CKD and eGFR were not included as covariates in the incident ESKD model because they are intermediate variables (not independent variables) on the causal pathway to ESKD. We performed a time-to-event analysis including the time from baseline until (1) a participant dies; (2) is registered in USRDS for ESKD; or (3) follow-up stops. The timing of events for each genotype, unadjusted for covariates, is provided in Supplemental Figure 1.
Missing Data Considerations
Missing data for the outcome measures and explanatory variables were rare (<0.5%) with some exceptions, e.g., complete blood count results (Table 1). A multiple imputation approach was employed in all multivariable regression models (Supplemental Material).45
Clinical and demographic characteristics according to HBA copy number
Sensitivity Analyses and Tests for Interaction
To address the possibility of association by population stratification, we performed a sensitivity analysis including the first ten principal components of ancestry among 7641 (77%) participants for whom Infinium Expanded Multi-Ethnic Genotyping Array data were available.
Prespecified tests for interaction between HBA copy number and sickle cell trait on the outcomes of CKD prevalence, incident reduced eGFR, and incident ESKD were performed for the fully adjusted multivariable models. Additionally, prespecified tests for interaction between HBA copy number and each of the following variables, APOL1, age, sex, and hypertension, were performed on the outcome of CKD prevalence for the fully adjusted multivariable model. A post hoc sensitivity analysis was performed in which HBA copy number was analyzed as a categoric variable instead of as a quantitative variable in the analysis of prevalent CKD and incident ESKD outcomes.
Relationships between the Structural Variant and Flanking SNPs
To exclude the possibility that the associations with kidney disease risk measured at the structural variant in the HBA1 and HBA2 loci were due to linkage disequilibrium with flanking SNPs, we analyzed SNP genotypes available on 8841 participants from the Infinium Expanded Multi-Ethnic Genotyping Array. Pairwise linkage disequilibrium was calculated between the −3.7-kb structural variant and each of 13,549 SNPs in the first 1 Mb of chromosome 16 using the r2 statistic in the R (v.3.6.2) package genetics. Associations between individual SNPs and CKD prevalence were calculated using PLINK (v.1.9) in a regression model that included age, sex, and the first four principal components of ancestry. The threshold for significance for a 1-Mb segment of the genome was estimated to be P<1.5 × 10−4 on the basis of the threshold for genome-wide significance (P<5 × 10−8) divided by the fraction of the genome (1 Mb of 3097 Mb) analyzed.46
Population Preventable Fraction
To estimate the fraction of kidney disease prevented by the protective genotypes of HBA copy numbers 2 or 3 among Black Americans, we calculated the population preventable fraction and the upper and lower bounds of the 95% CI for the fully adjusted models (further methods in the Supplemental Material).
Results
HBA Copy Number Variation
HBA copy number was variable among the 9918 participants who met inclusion and exclusion criteria (Figure 1): 393 participants (4%) had two HBA copies, 2744 (28%) had three HBA copies, 6668 (67%) had four HBA copies, 101 (1%) had five HBA copies, and two (<1%) participants had six HBA copies (Table 1). CKD prevalence increased with HBA copy number. CKD prevalence was 20% among those with two HBA copies; 25% among those with three copies; 25% among those with four copies; and 35% among those with ≥5 copies (P=0.008; Table 1). RBC parameters differed according to HBA copy number (P<0.0001): a hypochromic, microcytic anemia was observed in those with two copies of HBA.
HBA Copy Number and Kidney Disease Outcome Measures
Prevalent CKD
A total of 2476 (25%) participants had CKD at baseline. In an unadjusted analysis, each additional copy of HBA was associated with an 8% greater prevalence of CKD (PR, 1.08; 95% CI, 1.02 to 1.15; P=0.01). After adjusting for 13 CKD risk factors selected a priori, each additional copy of HBA was associated with a 14% greater prevalence of CKD (PR, 1.14; 95% CI, 1.07 to 1.21; P<0.0001; Table 2 ). In separate fully adjusted models analyzing the outcomes of prevalent eGFR<60 ml/min per 1.73 m2 or urine albumin-creatinine ratio ≥30 mg/g, each additional HBA copy increased the prevalence of reduced eGFR (PR, 1.22; 95% CI, 1.11 to 1.34; P<0.0001) and increased the prevalence of elevated urine albumin-creatinine ratio (PR, 1.09; 95% CI, 1.01 to 1.17; P=0.02; Table 2).
Association of HBA gene copy number with CKD prevalence—analysis adjusted for all listed covariates
Incident Reduced eGFR
Data were available for 3733 participants who were evaluated at a second in-home visit (Figure 1). A total of 655 (17.5%) developed incident reduced eGFR over a median (25th, 75th percentile) of 9.5 (8.7, 9.9) years. In unadjusted analysis there was no association between HBA copy number and incident reduced eGFR (RR, 1.01; 95% CI, 0.90 to 1.15; P=0.83).
After adjusting for relevant risk factors, there was no association between HBA copy number and incident reduced eGFR (RR, 1.06; 95% CI, 0.94 to 1.19; P=0.38; Table 3).
Association of HBA gene copy number with incident reduced eGFR—analysis adjusted for all listed covariates
Incident ESKD
Out of 9905 eligible participants for the ESKD analysis, 342 participants developed ESKD over a median (25th, 75th percentile) of 10.7 (5.5, 14.0) years of follow-up. In an unadjusted analysis there was no significant association between HBA copy number and the hazard of incident ESKD (HR, 1.07; 95% CI, 0.88 to 1.29; P=0.50; unadjusted survival curves are reported in Supplemental Figure 1). In contrast, an analysis adjusted for 13 risk factors selected a priori found that each additional copy of HBA was associated with a 32% increase in the hazard of incident ESKD (HR, 1.32; 95% CI, 1.09 to 1.61; P=0.005; Table 4).
Association of HBA gene copy number with incident ESKD– analysis adjusted for all listed covariates
Sensitivity Analyses and Tests of Interaction
The point estimates of the associations between HBA copy number and the outcomes of CKD prevalence, incident reduced eGFR, and incident ESKD were not substantially changed and remained significant for both CKD prevalence and incident ESKD after adjustment for the first ten principal components of ancestry (Supplemental Table 1). There was no interaction between HBA copy number and sickle cell trait for the outcomes of CKD prevalence, incident reduced eGFR, or incident ESKD (Supplemental Table 2), nor was there an interaction between HBA copy number and each of age, sex, hypertension, or APOL1 on the outcome of CKD prevalence (Supplemental Table 3). The P value of each crossproduct interaction term was ≥0.20 in the fully adjusted models. When HBA copy number was analyzed as a categoric variable instead of as a quantitative variable, the risks associated with each copy number genotype for prevalent CKD and incident ESKD were similar to risks determined from prespecified models that analyzed copy number as a quantitative variable; however, the risk estimates for the categoric analysis of the 5–6 copy category lacked precision due to the small number of observations for these genotypes (Supplemental Tables 4 and 5).
Relationships between the HBA Structural Variant and Nearby Sequence Variants
To assess whether associations with HBA copy number could be attributed to sequence variants flanking the HBA1 and HBA2 genes, we examined pairwise linkage disequilibrium between the −3.7-kb structural variant and 13,549 SNPs within the first 1 Mb of the chromosome. Pairwise linkage disequilibrium was weak between the structural variant and each SNP; the maximum r2 value was 0.25 and most SNPs had r2 values close to zero (Supplemental Figure 2, A and B). Each SNP was tested for association with CKD prevalence. None of the SNPs reached the threshold for significance for a 1-Mb genome region (Supplemental Figure 2, C and D), suggesting that the association between the structural variant and kidney disease risk was not due to linkage disequilibrium with SNPs in flanking gene regions.
Population Preventable Fraction
To assess the fraction of kidney disease prevented by the protective two- and three-copy HBA genotypes, we estimated the population preventable fraction. The two- and three-copy HBA genotypes prevent 4.3% (95% CI, 2.0 to 6.5) of CKD prevalence and 11.4% (95% CI, 2.2 to 20.6) of incident ESKD among Black Americans (Supplemental Material).
Discussion
In this national longitudinal cohort study of Black Americans, each additional HBA copy was associated with a 14% greater prevalence of CKD and a 32% greater risk of incident ESKD after adjustment for established genetic, biomedical, demographic, and social risk factors. In addition, there was a significant association between HBA copy number and each component of the prevalent CKD outcome, low eGFR and elevated urine albumin-creatinine ratio. Kidney disease risk increased with higher HBA copy number, suggesting that the level of HBA gene expression, rather than a specific allele, was responsible for this genetic association. These findings support the hypothesis that higher HBA copy number is associated with greater risk of kidney disease among Black Americans; in other words, the common 3.7-kb deletion is protective. Moreover, the protective two- and three-copy HBA genotypes (i.e., homozygotes and heterozygotes for the 3.7-kb gene deletion) were estimated to prevent approximately 4% of prevalent CKD and 11% of incident ESKD among Black Americans at the population level.
HBA copy number was associated with both a greater prevalence of CKD at baseline and greater incidence of ESKD over time; HBA copy number was not associated with incident reduced eGFR. To be included in the incident reduced eGFR analysis, participants had to undergo a second in-home visit. Participants in the group providing second in-home visit data were healthier (less smoking, hypertension, and diabetes mellitus at baseline) and had a higher socioeconomic status (higher education and income) at baseline than those who did not take part in the second in-home visit (see Supplemental Table 6). Together, the reduced disease severity, smaller sample size available for analysis (n=3733), and the stringent prespecified definition of reduced eGFR (>40% decrease among those with eGFR>60 ml/min per 1.73 m2 at baseline) may have reduced the power to detect an effect of HBA copy number on this secondary end point.
We performed ddPCR to quantify HBA copy number, which allowed a variety of underlying genotypes (−a/−a, −a/aa, aa/aa, aa/aaa, and aaa/aaa) to be summarized as an integer copy number (2, 3, 4, 5, 6) that reflects the number of functional HBA gene copies per genome. We evaluated a prespecified additive relationship between HBA copy number and kidney disease risk. Although the specific relationship between HBA copy number, HBA transcript abundance, and α-globin protein function in the kidney has yet to be defined, hematologic observations demonstrate a quantitative relationship between HBA copy number and α-globin protein abundance.35,47,48 Alternative models, in which HBA copy number was analyzed as a categoric variable instead of as a quantitative variable, further supported the conclusion that CKD and ESKD risk increased with HBA copy number (Supplemental Tables 4 and 5). Genetic associations can arise from population structure, i.e., underlying differences in genetic ancestry between those who develop kidney disease and those who do not. To address this, we performed a prespecified sensitivity analysis that included principal components of ancestry and confirmed that the association between HBA copy number and CKD prevalence was not due to population structure (Supplemental Table 1). The association with HBA copy number was unlikely due to linkage disequilibrium with sequence variants in nearby genes, given the established relationship between HBA copy number and functional gene expression; the rapid decline of linkage disequilibrium away from the HBA loci49; and the additive mode of inheritance (i.e., a dose-response relationship between HBA copy number and CKD risk) which transcends any single HBA allele. Nevertheless, we measured pairwise linkage disequilibrium between the structural variant and SNPs directly and measured association between each SNP and CKD prevalence. We found no SNPs that could explain the observed association between HBA copy number and kidney disease risk (Supplemental Figure 2). Changes in the level of HBA gene expression determined by functional gene copy number remains the most plausible genetic mechanism responsible for the associations with kidney disease risk.
α-Thalassemia is known to modify the risk of sickle cell disease nephropathy50,51; therefore, we considered whether the reduction in HBA copy number mitigated the effect of sickle cell trait on kidney disease risk. We found the association between HBA copy number and kidney disease risk to be independent of sickle cell trait, and there was no interaction between HBA copy number and sickle cell trait on kidney disease outcomes in the prespecified analysis. A previous study reported that HBA genotype modified CKD risk among those with sickle cell trait, but found no association between HBA and CKD in the general population.52 This study differs substantially from the prior report in the following ways: (1) ddPCR was used to quantify HBA copy number from 2 to 6; (2) a prespecified log-linear effect of HBA copy number was employed; (3) this study population was more than three times larger, with older individuals and a higher prevalence of CKD; and (4) additional established social, demographic, and clinical factors were included in the current analysis. Together, the precise estimation of copy number, increased sample size, and additional covariates may have increased the power to detect an association between HBA copy number and kidney disease risk in a general population sample of Black Americans.
There are several potential mechanisms through which HBA copy number variation could modify kidney disease risk among Black Americans. One such mechanism is through decreased erythroid expression of HBA leading to anemia, a risk factor for kidney disease. However, HBA alleles associated with protection against kidney disease tended to be associated with lower, not higher, hemoglobin levels. Moreover, the associations between HBA copy number and kidney disease risk remained significant when adjusted for hemoglobin level. We considered whether anemia associated with CKD could have contributed to overestimation of the effect of HBA copy number in the CKD prevalence model. In a post hoc sensitivity analysis omitting hemoglobin, the association between HBA copy number and CKD prevalence remained (Supplemental Table 7). Red cell microcytosis associated with lower HBA copy numbers could potentially confer protection against kidney disease through improved blood rheology53; however, this mechanism does not explain the increased risk associated with higher HBA copy numbers (5, 6) that have normal MCV. Thus, a clear hematologic mechanism explaining the association between HBA copy number and kidney disease risk is not apparent.
α-Globin is expressed not only in red cell precursors, but also in the endothelium of resistance arteries where it regulates nitric oxide signaling.28 In human omental arteries, α-globin binds with eNOS at the myoendothelial junction where it regulates vasoconstriction to phenylephrine in a NOS-dependent fashion.54(preprint) In mice, pharmacologic disruption of the α-globin/eNOS complex in endothelial cells increased nitric oxide release and led to a fall in BP in both normal and hypertensive animals, but had no effect in endothelial eNOS–knockout mice.29,30 Genetic deletion of Hba or α-hemoglobin stabilizing protein (Ahsp) in mice altered small artery vasoreactivity by reducing vasoconstriction to α-adrenergic stimuli in a NOS-dependent manner.28,31 Individuals with HBA deletions had enhanced tissue perfusion after transient arterial occlusion or in response to temperature changes, consistent with increased endothelial nitric oxide release.55,56 Together, these human and mouse studies identify variation in α-globin expression to be associated with changes in vascular nitric oxide signaling and vasoconstriction to α-1–adrenergic stimuli. We found HBA copy number to be independently associated with each component of CKD, low eGFR or elevated urine albumin-creatine ratio. We speculate that decreased endothelial HBA expression in individuals with decreased HBA copy number could lead to increases in endothelial nitric oxide signaling in human arteries and decreased vasoconstriction in response to sympathetic nervous system activity. This could confer protection against kidney disease in a manner consistent with the protective role of increased nitric oxide signaling in experimental and human clinical studies of kidney injury.15,16,21,22,57 Further investigation is necessary to understand how reduced HBA expression affects renal function, response to injury, and progression of kidney disease.
This study has strengths and limitations. The REGARDS study is one of the largest available cohorts of Black Americans and has clearly defined kidney disease outcome measures and data on social, biomedical, and genetic factors. This allows adjustment for established kidney disease risk factors when estimating risk. We prespecified the analysis plan to reduce the chance of spurious associations. We employed a robust and quantitative ddPCR method to genotype the insertion/deletion polymorphism and the gene triplication in the HBA loci. We evaluated genetic risk factors for CKD specific to Black individuals, including sickle cell trait and high-risk APOL1 status, and determined that HBA copy number was associated with kidney disease independently of these factors. This analysis focused on a Black American population, in which the 3.7-kb gene deletion is common. A limitation of this study is the lack of an independent replication cohort. Although the REGARDS study provides one of the largest longitudinal cohorts with which to study kidney disease in Black Americans, the novel associations reported here would be further strengthened by replication in similarly powered cohorts of Black Americans and in other populations in which HBA deletions are common.
The major genetic variants that determine kidney disease risk among Black Americans arose in Africa as mutations conferring resistance to infectious diseases: HBA and HBB variants (like those that cause α-thalassemia and sickle cell trait) confer protection against malaria whereas APOL1 variants protect against trypanosomiasis, also known as African sleeping sickness.58⇓–60 In contrast to the variant alleles of HBB (such as sickle cell trait) and APOL1, which confer protection against infectious diseases at the cost of increased kidney disease risk, the 3.7-kb deletion in HBA is associated with protection against both malaria and kidney disease (Supplemental Figure 3).
In conclusion, we report that higher HBA copy number was independently associated with greater CKD prevalence and ESKD incidence after accounting for known clinical, demographic, and genetic risk factors in this national longitudinal study of Black Americans. The high frequency of the HBA gene deletion found in Black Americans may act to reduce the overall burden of kidney disease in this population.
Disclosures
M. Cushman reports Scientific Advisor or Membership: International Society on Thrombosis and Haemostasis—journal editor, Research and Practice in Thrombosis and Haemostasis. O.M. Gutierrez discloses receiving grant funding and consulting fees from Akebia Therapeutics; grant funding and consulting fees from Amgen; grant funding from GlaxoSmithKline; consulting fees from QED Therapeutics; and Honoraria from Akebia, Amgen, Ardelyx, AstraZeneca, and Reata. J.B. Kopp reports a Patent (Vpr antibodies) with the National Institutes of Health; Scientific Advisor or Membership as a Board Member of the National Kidney Foundation serving the National Capitol Area, NIH Foundation for the Advancement of Education in the Sciences; Editorial Boards of the American Journal of Nephrology and American Journal of Physiology—Renal; and Other Interests/Relationships with the National Institute of Diabetes and Digestive and Kidney Diseases via Confidentiality disclosure agreement with Intermune, Ionis, and Third Rock; Collaboration research and development agreement with Sanofi; and Research collaboration agreement with Astra Zeneca, Genentech, Merck, and Third Rock. B.T. Mott reports Consultancy Agreements with Geminus Therapeutics; Ownership Interest in Geminus Therapeutics; and Patents and Inventions with Chinook Therapeutics and Geminus Therapeutics. R.P. Naik reports Consultancy Agreements with Elsevier; and Research Funding from Rigel. L.H. Pecker reports Scientific Advisor or Membership as Sickle Cell Reproductive Health Education Directive. N.A. Zakai reports Honoraria from the American Society of Hematology and the International Society on Thrombosis and Haemostasis; and Scientific Advisor or Membership via Editorial Board for Research and Practice of Thrombosis and Haemostasis, Various committees for the Hemophilia and Thrombosis Research Society, and the American Society of Hematology. All remaining authors have nothing to disclose.
Funding
This is an ancillary study supported by cooperative agreement U01 NS041588 cofunded by the National Institute of Neurological Disorders and Stroke (NINDS) and the National Institute on Aging (NIA), National Institutes of Health, Department of Health and Human Service. This research was supported in part by the Divisions of Intramural Research, National Institute of Allergy and Infectious Diseases (NIAID) (project AI001150) (A.P. Ruhl, H.C. Ackerman) and the National Heart, Lung, and Blood Institute (NHLBI) (project HL006196) (A.P. Ruhl, Y. Yang, H.C. Ackerman). This work was also funded in part by the National Cancer Institute (NCI) Intramural Research Program under contract HHSN26120080001E (C.A. Winkler), the National Institute of Diabetes and Digestive and Kidney Diseases (NIDDK) Z01 DK04312 (J.B. Kopp), and the NHLBI grants K08HL12510 (R.P. Naik) and K08HL096841 (N.A. Zakai).
Acknowledgments
The authors thank the investigators, staff, and participants of the Reasons for Geographic and Racial Differences in Stroke (REGARDS) study for their valuable contributions. A full list of participating REGARDS investigators and institutions can be found at http://www.regardsstudy.org. The content is solely the responsibility of the authors and does not necessarily represent the official views of the NINDS,NIA,NIAID,NCI,NIDDK, or NHLBI. The content of this publication does not necessarily reflect the view or policy of the Department of Health and Human Services, nor does mention of trade names, commercial products, or organizations imply endorsement by the government. Some of the data reported here have been supplied by the United States Renal Data System (USRDS). The interpretation and reporting of these data are the responsibility of the authors and in no way should be seen as an official policy or interpretation of the US government.
Dr. A. Parker Ruhl and Dr. Hans C. Ackerman conceived and designed the study and drafted the first version of the manuscript. Dr. Yu Yang, Dr. Bryan T. Mott, Dr. Cheryl A. Winkler, and Dr. Jeffrey B. Kopp performed sample analyses. Dr. Leslie A. Lange performed the principal components of ancestry analysis. Dr. Amit Patki performed the linkage disequilibrium analysis. Dr. Neal Jeffries designed and implemented the statistical analyses. All authors contributed to data interpretation and revised and approved the manuscript.
Supplemental Material
This article contains the following supplemental material online at http://jasn.asnjournals.org/lookup/suppl/doi:10.1681/ASN.2021050653/-/DCSupplemental.
I. Supplemental Figure and Tables
Supplemental Table 1. Association of HBA copy number with CKD prevalence, incident reduced eGFR, and incident ESKD—fully adjusted models including ten principal components of ancestry.
Supplemental Table 2. Prespecified test for interaction between HBA copy number and SCT on the outcomes of CKD prevalence, incident reduced eGFR, and incident ESKD in fully adjusted models.
Supplemental Table 3. Prespecified tests for interaction between each of age, sex, hypertension, or APOL1 and HBA on the outcome of CKD prevalence in fully adjusted models.
Supplemental Table 4. Post-hoc sensitivity analysis of the association of categorical HBA copy number with CKD prevalence.
Supplemental Table 5. Post-hoc sensitivity analysis of the association of categorical HBA copy number with ESKD incidence.
Supplemental Table 6. Comparison of participants with and without second visit data.
Supplemental Table 7. Post-hoc sensitivity analysis of the association of HBA copy number with CKD prevalence when hemoglobin is omitted from the model.
Supplemental Figure 1. Unadjusted survival curves of ESKD onset by genotype.
Supplemental Figure 2. Linkage disequilibrium and association analysis of sequence variants in the first 1 Mb of chromosome 16 flanking HBA1 and HBA2.
Supplemental Figure 3. Comparison of prevalence ratios for HBA, APOL1, and HBB genetic risk factors for CKD.
II. Estimation of population preventable fraction of HBA copy number on kidney disease.
III. Additional Methods.
a. HBA Genotyping Methods
b. APOL1 Genotyping Methods
c. ESKD Time-to-event Analysis Methods
d. Multiple Imputation Procedure
e. Diagnostic Modeling Description
f. Assessment of the Missing at Random Assumption
IV. Supplemental References.
Footnotes
Published online ahead of print. Publication date available at www.jasn.org.
- Copyright © 2022 by the American Society of Nephrology