ORIGINAL RESEARCH article

Front. Physiol., 10 September 2026

Sec. Red Blood Cell Physiology

Volume 17 - 2026 | https://doi.org/10.3389/fphys.2026.1908536

Erythroid anion exchanger-1 (SLC4A1) modulates carbonic anhydrase inhibitor acetazolamide on respiration

  • 1. Department of Medicine, Mackay Medical University, New Taipei City, Taiwan

  • 2. Division of Colon and Rectal Surgery, Department of Surgery, MacKay Memorial Hospital, Taipei, Taiwan

  • 3. The Laboratory of Immunogenetics, Department of Medical Research, MacKay Memorial Hospital, New Taipei City, Taiwan

  • 4. Graduate Institute of Biomedical Sciences, College of Medicine, MacKay Medical University, New Taipei City, Taiwan

  • 5. School of Pharmacy, College of Pharmacy, Kaohsiung Medical University, Kaohsiung, Taiwan

  • 6. MacKay Junior College of Medicine, Nursing, and Management, New Taipei City, Taiwan

Abstract

Introduction:

Clinically widely used acetazolamide is a potent carbonic anhydrase (CA) inhibitor. CA accelerates CO2(g)⇌HCO3-(aq) bidirectional conversion. CAII, the major functional isoform, primarily resides inside red blood cells (RBCs). Erythrocyte-specific anion transporter or anion exchanger-1 (AE1) determines HCO3- permeability across the cell membrane and is considered the rate-limiting factor for intraerythrocytic CO2(g)/HCO3(aq) conversion facilitated by CAII. This study aimed to use CA-inhibiting acetazolamide to find out whether AE1 could modulate intraerythrocytic CA catalysis.

Methods:

To explore whether erythroid AE1 could modulate the respiratory effects of acetazolamide, we utilized the GPMur mouse model characteristic of higher erythroid AE1 expression for respiratory measurements using whole-body plethysmography (WBP). The effects of acetazolamide and hypercapnia on murine respiration were compared.

Results:

We found that acetazolamide influenced respiration in the same direction as hypercapnia. Increased erythroid AE1 expression (GPMur) counteracted the impacts of hypercapnia/acetazolamide. Importantly, GPMur/higher AE1 increased respiratory sensitivity and limited responses to acetazolamide.

Conclusion:

Both acetazolamide and hypercapnia drove acidosis. GPMur/increased AE1 reduced the impacts of acidosis with faster respiratory responses. This supports that AE1 and CAII function in concert to facilitate intraerythrocytic CO2/HCO3- conversion for CO2 expiration. Since murine RBC-AE1 could help sensitize respiratory responses to acetazolamide-induced systemic CA inhibition, it warrants future clinical investigation to identify the appropriate dosing of acetazolamide for people with the GP.Mur blood type.

Introduction

Acetazolamide has been used to treat a wide range of illness for over 70 years despite the fact that many of its mechanisms of action are still not completely understood. Acetazolamide inhibits carbonic anhydrases (CA) and reduces the rate of CO2(g)/HCO3(aq) conversion. CA isoforms are ubiquitously present in different tissues (i.e., erythrocytes, eyes, lungs, kidneys, and other gastrointestinal organs as well as choroid plexus, and various brain cells (; )), and thus the inhibiting effects of acetazolamide could be global. Acetazolamide stimulates respiration and remains the first-line prophylactic drug for acute mountain sickness (AMS). Acetazolamide may be used as an add-on medication to treat metabolic alkalosis, particularly diuretic-induced or chloride-resistant cases; it enhances bicarbonate excretion when conventional electrolyte and volume repletion are insufficient or impractical (; ). Acetazolamide also functions as a diuretic, which has been used to lower intraocular pressure (IOP) in glaucoma, to relieve special types of migraine associated with idiopathic intracranial hypertension or CADASIL (cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy) (; ; ; ; ; ; ), and to help decongestion in acute heart failure with volume overload (; ). Currently, acetazolamide is included in the World Health Organization’s list of essential medicines.

Both types I and II of carbonic anhydrases (CAI and CAII) are abundantly present inside human red blood cells. CAII carries out most catalytic activities for CO2/HCO3- conversion inside human and murine erythrocytes (). CAII is the sole catalytic isoform of carbonic anhydrase in murine RBCs (; ); this was demonstrated by CAII-deficient mice that present metabolic and respiratory acidosis (insufficient CO2 excretion) ().

There are roughly ~480,000 CAII and ~1 million AE1 protein molecules in each human RBC (). It has been suggested that CAII could directly bind to AE1 in the submembranous region, but there is also evidence against CAII–AE physical interaction (; ). AE1 allows HCO3- efflux or influx following the direction of the anion concentration gradients across the RBC membrane. A major source of blood HCO3- is from CAII-catalyzed hydration of CO2 that takes place inside red blood cells. Because the main components that facilitate CO2(g)/HCO3-(aq) conversion and HCO3-/CO2 transport (i.e., CAII, AE1, deoxygenated hemoglobin (deoxy Hb), and AQP1) are all in proximity of one another in the RBC membrane region, the AE1-centered multi-protein complex has been proposed to form a “CO2-transport metabolon (or hub)” (; ; ; ; ; ). In spite of the controversy about CA–AE1 physical interaction inside RBCs (; ), AE1 (integral membrane protein) and CAII (cytosol protein) share the same substrate HCO3- (which is constantly converted from/to CO2 driven by local concentration gradients inside red blood cells); AE1 and CAII are thus functionally linked or coupled (Figure 1A). In their functional coupling, CAII catalytic turnover (~106/s) far exceeds the capacity of AE1-mediated transport (~50,000–60,000 anions/s) (; ), making AE1 the rate-limiting step. Whether AE1 could modulate CA activities or the effect of acetazolamide remains unclear.

Figure 1

GP.Mur is a clinically important Southeast Asian (SEA) erythrocyte type (with prevalence ranging from 1% to 2% in Malaysia and Indonesia, 9%–18% in Thailand, 6%–9.2% in Vietnam, 7%–10% in the Philippines, 18% in Eastern Taiwan, and 5%–10% in Southern China) (; ; ; ; ; ; ; ; ; ; ; ; ) but rare in non-SEA regions (even northern Asia such as Japan, Korea, and northern China). The human GYP.Mur gene is a naturally occurring recombinant of erythroid glycophorin B (GYPB) and glycophorin A (GYPA) genes. Similar to glycophorin A that promotes AE1 expression, GPMur (glycophorin B-A-B hybrid), as the protein entity of human GP.Mur blood type, also enhances erythroid AE1 levels. Notably, the GPMur blood type is characteristic of 10% or higher erythroid AE1 expression (; ). The higher AE1 expression accelerates transmembrane transport of HCO3- and conceivably drives the conversion between HCO3- and CO2 inside the erythrocytes (). This was supported by the fact that people with the GP.Mur blood type excrete CO2 more efficiently during low- and high-intensity exercise (; ). The phenomenon that GPMur glycophorin B-A-B hybrid protein enhances AE1 expression on human RBCs was verified by a newly established GPMur knock-in (KI) mouse model: human GPMur protein expression in the transgenic mice increased murine erythroid AE1 by ~55% (). Thus, in this study, we adopted the GPMur murine model to explore whether increased RBC AE1 could modulate acetazolamide’s respiratory effects by whole-body plethysmograph or WBP (Figure 1B).

Methods

GYP.Mur humanized knock-in mouse

Since drug-induced respiratory responses could only be studied in mammals and not organoid or cultured cells, this mouse study has been approved by the Institutional Animal Care and Use Committee (IACUC) of Taiwan MacKay Memorial Hospital (MMH-A-S-111-35-R-R) and adhered to the ARRIVE guidelines.

The GYP.Mur humanized KI mouse model (C57BL/6JNarl-Gypaeml(GYPMur)MMH; GPMur KI) was generated by CRISPR/Cas9-mediated DNA cleavage and homologous recombination with a replacing ssDNA on the Gypa allele of C57BL/6JNarl (a sub-strain bred from C57BL/6J; National Center for Biomodels [NCB], Taipei, Taiwan) as described (). For the knock-in strategy, human GYP.Mur cDNA (accession EU338225) was inserted into the 3′ UTR of murine GYPA gene; with this strategy, the expression of human GYP.Mur gene was driven by the promoter of murine erythroid-specific GYPA gene. Knock-in of the GYP.Mur gene was verified by PCR sequencing using genomic DNA extracted from the tail tips of pups. GPMur KI mice produce progenies with equal male/female ratios. GYP.Mur knock-in drives GPMur protein expression on the murine RBC membrane and does not harm or reduce the life expectancy of C57BL6/Narl (). The control or wild-type (WT) mice were age-matched C57BL6/Narl mice purchased from Taiwan-NCB.

Flow cytometric measurement of CAII in mouse RBCs

The RBCs from adult male mice of both strains (C57BL/6JNarl-Gypaeml(GYPMur)MMH [GPMur KI]; C57BL/6JNarl [WT]) were washed with PBS, fixed with 1% paraformaldehyde in PBS at 4 °C for 10 min, and washed with PBS again, followed by permeabilization using 0.05% saponin in PBS/5% fetal bovine serum (FBS; Thermo Fisher Scientific, Waltham, MA, USA; catalog no. 16000044). The fixed and permeabilized RBCs were washed using PBS containing 1% bovine serum albumin (BSA) and incubated with anti-CAII antibody (1:200 dilution; GeneTex, Hsin-Chu, Taiwan) at 4 °C for 1 h, followed by 30 min of secondary antibody labeling with Alexa-Fluor-488-conjugated anti-rabbit Fab (1:200 dilution; Jackson ImmunoResearch, West Grove, PA, USA) at 4 °C in the dark. Intraerythrocytic CAII protein levels were compared based on the intensities of Alexa Fluor 488 fluorescence using FACSCalibur (BD Biosciences). The geometric mean of the WT data was set as 100%, and all data were calculated in percentile (%) relative to the geometric mean of the WT data.

Whole-body plethysmography

From whole-body plethysmography (WBP) pressure measurements, the following ventilatory parameters were calculated: respiratory frequency (Rf, breaths/min), tidal volume (VT, mL), and minute ventilation ( E = Rf × VT, mL/min). In a respiratory cycle, the inspiratory time (Ti, s) and expiratory time (Te, s) can be expressed as duty cycles Ti/Ttot and Te/Ttot, where the total time of one breath or Ttot = Ti + Te. The inspiratory drive and expiratory drive are defined as VT/Ti and VT/Te (mL/s). Neural control of respiration could also be evaluated by the end-inspiratory and end-expiratory pauses (EIP and EEP, ms). The parameters relevant to expiratory kinetics are relaxation time (Tr, s) or the time for box-pressure decay to 36% of the total expiratory pressure and estimated peak expiratory and inspiratory flows (PEFb and PIFb, mL/s).

Two consecutive generations of the GPMur KI and age-matched wild-type (WT) adult male B6 mice (six mice per group × two generations) were tested for respiratory functions using Buxco® WBP FinePointe (Data Sciences International) under unrestrained and conscious conditions at the core facility of the Taiwan Mouse Clinic. To follow the regulation for WBP experimentation using the DSI instrument, each WBP experiment was limited to maximally one test per week and ≤1.5 h per test to minimize the time that a study mouse was put alone in a WBP chamber.

Murine respiratory functions were challenged by hypercapnia and/or acetazolamide. The hypercapnic challenge took place after 45-min stabilization in air flow. For the hypercapnic condition, we used premixed gas containing 5% CO2, 76% N2, and 21% O2. Gas-in of 5% CO2 has been shown to be tolerable for B6J mouse studies (). Each mouse was tested with the same WBP-hypercapnia protocol twice (one with the pre-treatment of acetazolamide and the other without); the two tests were scheduled 1 week apart (see Figure 2A—WBP protocol). To prevent habituation with hypercapnia conditioning or acetazolamide, the animals were not pre-trained in this experimental setting. Each WBP test began instead with 45 min of stabilization in air flow, which allowed individual animals to calm down in their own WBP chambers. As shown in Figure 2A—WBP protocol, after the WBP chamber was injected with air (~0.04% CO2) for 45 min (stabilization phase), pre-mixed high-CO2 gas was injected for 30 min (hypercapnic phase) and changed back to air (recovery phase). In the following week, the same mice were individually weighed and then given acetazolamide (50 mg/kg) by oral gavage before WBP testing. Prior to WBP, their body weights were within the normal range for adult B6 mice (GPMur KI: 29.4 ± 2.0 g; WT: 31.5 ± 2.5 g).

Compared to the typical dosage of acetazolamide for treating human patients (8–50 mg/kg/day in divided doses) (; ), a slightly higher dosage of acetazolamide (40–50 mg/kg) is generally used in rodent respiratory studies to account for human–mouse differences in metabolism and pharmacokinetics (). In our pre-test to determine an appropriate dose for WBP, we examined murine blood gas 1 h after oral gavage with different doses of acetazolamide (0–100 mg/kg BW). We found that the blood pCO2 between WT and the GPMur mice differed most under the dose of 50 mg/kg acetazolamide, and thus this dose was used for WBP. Others have also studied posthypoxic unstable breathing in C57BL/6 mice using 40 mg/kg acetazolamide ().

Statistical analyses

IBM SPSS and OriginPro programs were used for all statistical computations and graphic preparations. Since WBP could be affected by multiple factors (acetazolamide, hypercapnia interval, GPMur phenotype, and body weight), linear mixed model (LMM) was used as the primary confirmation. Descriptive statistics and paired and unpaired t-tests were used for exploratory comparisons to characterize the data. The respiratory data collected at 30–40 min (air), 55–65 min (hypercapnia), and 85–92 min (air/recovery), respectively, were averaged and analyzed using t-tests and LMM. Unpaired t-tests were used to compare GPMur and WT groups. Paired t-tests were used to compare the same mice with and without acetazolamide treatments. In LMM, the hypercapnia-interval factor comprised baseline air and hypercapnia stages; the recovery stage was excluded due to large variations among animals during return to the baseline. We used the main effects model of the LMM to evaluate the overall significance of each factor; this model considers the overall effect of each factor to be independent from the other factors. To account for repeated measures of the same mouse, the ID of each mouse was included as a random effect. LMM in SPSS employs the Satterthwaite method to estimate the degrees of freedom for limited sample sizes. Statistical significance was set at p <0.05.

Results

One-fifth less CAII in murine RBC expressing GPMur

Since CAII is the main catalytically active isoform of carbonic anhydrase in human and murine RBCs (; ), we first examined whether the expression of human GPMur protein could affect murine erythroid CAII levels in knock-in mice. By flow cytometric quantitation of fixed and permeabilized murine RBCs, we found ~22.4% less CAII in RBCs from GPMur KI than from WT B6 mice (77.6 ± 17.8% [GPMur, n = 11] versus 100 ± 13.6% [WT, n = 8]; **p = 0.009).

GPMur/AE1-sensitized ventilatory responses to acetazolamide/hypercapnia

To visualize how RBC AE1 could modulate acetazolamide’s respiratory effects on GPMur knock-in mice (Figure 1A) (), we plotted the averaged traces from WBP (Figures 2, 3; Supplementary Figures). Acetazolamide and 5% CO2 treatments both aimed to increase blood CO2 and trigger acidosis to activate central and peripheral chemoreceptors (hypercapnic ventilatory responses) (; ). For unrestrained and unanesthetized adult C57BL/6 mice, the reported VT values estimated by Quindry et al. are approximately 9 to 10 mL/kg and Rf 350–468/min (). The ventilatory data of our mice were generally near the upper ends of the reported ranges, presumably because our mice were ~20% heavier and 2 to 3 months older and/or WBP settings could vary among different studies (; ; ; ).

Figure 2

Figure 3

Table 1 summarizes the WBP-respiratory data and presents the initial exploratory analyses in the three gas-in conditions of WBP (baseline air/hypercapnia/recovery air). Acetazolamide and hypercapnia substantially affected almost all respiratory parameters in both WT and GPMur mouse strains. GPMur knock-in lowered VT, although its effects on respiration were mild (Table 1).

Table 1

WBP parameterUnitConditionWT (12)WT+ACZ (12)P-value (WT vs. WT+ACZ)aGPMur (10)GPMur+ACZ (10)P-value (GPMur vs. GPMur+ACZ)aP-value (WT vs. GPMur)P-value (WT+ACZ vs. GPMur+ACZ)
Rf1/minAir333.9 ± 53.1338.0 ± 40.3353.3 ± 64.1387.8 ± 62.6**
Hypercapnia370.4 ± 43.4361.5 ± 17.9356.2 ± 41.3364.4 ± 31.7
Recovery304.5 ± 53.1337.1 ± 24.8*313.5 ± 75.6368.8 ± 46.2*
VTmLAir0.277 ± 0.0210.290 ± 0.0230.257 ± 0.0150.276 ± 0.0280.1**
Hypercapnia0.393 ± 0.0240.441 ± 0.034**0.360 ± 0.0420.415 ± 0.030****0.07
Recovery0.237 ± 0.0230.270 ± 0.025**0.220 ± 0.0250.260 ± 0.030**
EmL/sAir94.6 ± 20.997.7 ± 13.891.9 ± 17.8107.7 ± 24.4*
Hypercapnia144.9 ± 22.3158.8 ± 14.70.06127.2 ± 17.7150.6 ± 17.5*0.051
Recovery73.6 ± 19.491.3 ± 12.4*70.5 ± 20.895.6 ± 15.5**
TesAir0.171 ± 0.0320.128 ± 0.018**0.157 ± 0.0270.110 ± 0.019****
Hypercapnia0.103 ± 0.0130.095 ± 0.0050.050.106 ± 0.0120.094 ± 0.009**
Recovery0.179 ± 0.0220.131 ± 0.016***0.183 ± 0.0510.119 ± 0.011**0.055
Te/TotAir0.733 ± 0.0270.650 ± 0.019***0.726 ± 0.0200.644 ± 0.021***
Hypercapnia0.597 ± 0.0120.565 ± 0.008***0.595 ± 0.0160.562 ± 0.013***
Recovery0.723 ± 0.0440.662 ± 0.029***0.729 ± 0.0150.660 ± 0.031***
TisAir0.061 ± 0.0120.068 ± 0.010*0.059 ± 0.0120.061 ± 0.011
Hypercapnia0.069 ± 0.0070.073 ± 0.0040.070.072 ± 0.0090.073 ± 0.007
Recovery0.069 ± 0.0200.066 ± 0.0040.067 ± 0.0180.061 ± 0.010
Ti/TotAir0.267 ± 0.0270.350 ± 0.019***0.274 ± 0.0200.356 ± 0.021***
Hypercapnia0.403 ± 0.0120.435 ± 0.008***0.405 ± 0.0160.438 ± 0.013***
Recovery0.277 ± 0.0440.338 ± 0.029***0.271 ± 0.0150.340 ± 0.031***
VT/TemL/sAir1.75 ± 0.462.38 ± 0.33**1.74 ± 0.302.63 ± 0.65**
Hypercapnia3.96 ± 0.664.66 ± 0.42**3.49 ± 0.554.45 ± 0.56**
Recovery1.40 ± 0.302.17 ± 0.39***1.33 ± 0.362.25 ± 0.44***
VT/TimL/sAir4.78 ± 0.944.39 ± 0.614.59 ± 0.804.74 ± 0.96
Hypercapnia5.86 ± 0.816.06 ± 0.595.13 ± 0.695.70 ± 0.610.05*
Recovery3.80 ± 1.044.21 ± 0.453.62 ± 0.984.34 ± 0.63*
PAUAir1.04 ± 0.181.19 ± 0.14*0.96 ± 0.141.02 ± 0.12***
Hypercapnia1.20 ± 0.081.23 ± 0.071.10 ± 0.101.09 ± 0.09****
Recovery1.15 ± 0.161.13 ± 0.121.10 ± 0.141.00 ± 0.100.08*
PenhAir0.673 ± 0.1690.876 ± 0.128**0.621 ± 0.1100.714 ± 0.124***
Hypercapnia1.10 ± 0.161.38 ± 0.12***0.96 ± 0.161.15 ± 0.14**0.055***
Recovery0.832 ± 0.2490.836 ± 0.1280.727 ± 0.0930.713 ± 0.133*
TrsAir0.094 ± 0.0180.062 ± 0.007***0.090 ± 0.0110.057 ± 0.009***
Hypercapnia0.048 ± 0.0060.043 ± 0.0030.050.052 ± 0.0060.046 ± 0.006*
Recovery0.093 ± 0.0130.065 ± 0.012***0.095 ± 0.0160.063 ± 0.006**
PEFbmL/sAir4.63 ± 1.214.74 ± 0.724.35 ± 0.864.82 ± 1.10
Hypercapnia7.30 ± 1.429.07 ± 1.06***6.06 ± 1.148.05 ± 0.86****
Recovery3.83 ± 0.874.40 ± 0.680.113.49 ± 1.014.35 ± 0.83*
PIFbmL/sAir8.00 ± 1.457.12 ± 1.020.067.48 ± 1.217.63 ± 1.53
Hypercapnia9.08 ± 1.238.91 ± 0.878.01 ± 0.948.44 ± 0.98*
Recovery6.42 ± 1.726.78 ± 0.646.10 ± 1.516.98 ± 0.980.11
EEPmsAir54.3 ± 15.816.2 ± 12.6***36.9 ± 9.19.9 ± 8.6***
Hypercapnia2.56 ± 1.050.83 ± 0.44***3.92 ± 4.991.58 ± 3.41****
Recovery59.2 ± 22.119.2 ± 10.8***63.6 ± 47.717.0 ± 9.5*
EIPmsAir2.20 ± 0.092.49 ± 0.22**2.28 ± 0.162.40 ± 0.25
Hypercapnia2.22 ± 0.112.59 ± 0.16***2.26 ± 0.172.63 ± 0.20***
Recovery2.33 ± 0.192.36 ± 0.192.38 ± 0.242.43 ± 0.27

WBP data of WT and GPMur mice upon challenge with acetazolamide (ACZ) and hypercapnia, respectively.

Individual WBP data from the three time intervals (1 [30~40 min], 2 [55~65 min], and 3 [85~92 min]) were averaged and shown in mean±SD. WT and GPMur mice were compared by unp1ed two-sample t-tests. aData from the same mice with and without ACZ treatments were compared by p1ed t-tests. *P-values < 0.05, **P-values < 0.01 denoted **, ***P-values < 0.001. The number of mice tested per group is indicated inside the parenthesis.

We used LMM as the confirmatory analyses which accounted for all four factors (hypercapnia, acetazolamide, GPMur expression, and body weight [BW]) and repeated measures of the same mice and found that (1) hypercapnia and acetazolamide were predominant factors influencing WBP-respiratory parameters and (2) GPMur still showed significant lowering effects on VT (Figure 2; Table 2). The latter is consistent with previous findings from the cardiopulmonary exercise testing (CPET) of college athletes: GP.Mur+ athletes persistently showed lower VT throughout CPET (). Differently from GP.Mur+ athletes who showed significantly higher Rf and larger E during CPET, here unrestrained GPMur and WT mice inside WBP chambers did not differ significantly in Rf. This could be due to a lack of strong physical stimulation (e.g., running on a treadmill) and large Rf fluctuation associated with the hypersensitive nature of the B6J strain (lowering the statistical power).

Table 2

Effects on murine WBP parametersBody weightGPMur typeAcetazolamideIncrease to 5% CO2
p-valueEffectsp-valueEffectsp-valueEffectsp-valueEffects
Rf (/min)n.s.n.s.n.s.n.s.
VT (mL)p < 0.001↑p < 0.05↓p < 0.001↑p < 0.001↑
E (mL/min)n.s.n.s.p < 0.01↑p < 0.001↑
Te (s)n.s.n.s.p < 0.001↓p < 0.001↓
Te/Ttotn.s.n.s.p < 0.001↓p < 0.001↓
Ti (s)n.s.n.s.n.s.p < 0.001↑
Ti/Ttotn.s.n.s.p < 0.001↑p < 0.001↑
VT/Ten.s.n.s.p < 0.001↑p < 0.001↑
VT/Tin.s.n.s.n.s.p < 0.001↑
PAUp < 0.01↑p < 0.001↓n.s.p < 0.001↑
Penhp < 0.001↑p < 0.01↓p < 0.001↑p < 0.001↑
Tr (s)n.s.n.s.p < 0.001↓p < 0.001↓
PEFb (mL/s)n.s.n.s.p < 0.001↑p < 0.001↑
PIFb (mL/s)n.s.n.s.n.s.p < 0.001↑
EEP (ms)n.s.n.s.p < 0.01↓p < 0.001↓
EIP (ms)n.s.n.s.p < 0.001↑p < 0.05↑

Assessing the effects of the four variables—body weight, GPMur phenotype, acetazolamide treatment, and CO2 increase to 5% (hypercapnia)—on murine respiratory parameters by linear mixed models (LMM main effects).

P >0.05 (statistically not significant) is denoted n.s. The direction of the effects of a variable is determined from “estimates of fixed effects” in LMM and is indicated by an upward arrow for increasing effects or a downward arrow for decreasing effects.

Expiratory kinetics affected by acidosis and GPMur/AE1 in opposite directions

Expiratory kinetics was characterized by relaxation time (Tr) and peak expiratory/inspiratory flow (PEF/PIF) (Figure 3). Acetazolamide and hypercapnia both triggered acidosis and significantly increased PEF and PIF, indicating a need for high airflow to compensate respiratory adjustments and support acid–base homeostasis. Tr is directly derived from the measured volume curve and is defined as the time from the start of expiration until the time that two-thirds of the tidal volume has been exhaled in a respiratory cycle. In a normal passive exhalation, Tr is essentially the time constant of lung exhalation; prolonged Tr may mean an increase of airway resistance.

The definition of Tr and (Te - Tr) reflects kinetically distinct early and late expiratory phases in a breathing pattern (Figure 1B). PAU compares the time on the early, active part of exhalation (Tr) to the later part of exhalation (Te - Tr) in a breath and is defined as the ratio of (Te - Tr) to Tr. Higher PAU means that the mouse may be in labored breathing, as the animal is spending a longer time in the later than the early active stage of exhalation. Penh is PAU times the ratio of the peak expiratory flow to the peak inspiratory flow (PEF/PIF). PAU is a value of expiratory kinetics, whereas Penh reflects not just expiratory kinetics but also waveforms. Both hypercapnic and acetazolamide challenges increased Penh and PAU; in contrast, GPMur expression persistently lowered Penh and PAU (Table 2, Figures 3A–C). Although whether Penh could be used to assess airway broncoconstriction is controversial (), the observation that GPMur/increased AE1 significantly affected murine expiratory kinetics (Table 2) is reminiscent of similar findings from the human GP.Mur CPET study ().

Finally, we probed into neural control of breathing rhythm. Both hypercapnia and acetazolamide shortened expiratory time (Te) and its fraction in a respiratory cycle (Te/Ttot), while Ti and Ti/Ttot were increased (Supplementary Figure 2). This was driven by the need to increase tidal volume and ventilation upon the onset of acidosis. Both stimulations increased expiratory drive and inspiratory drives (VT/Te and VT/Ti, respectively) in the unrestrained mice (Supplementary Figure 3). Indeed hypercapnia has been reported to trigger active expiration to expand tidal volume by shortening Te and Te/Ttot in unanesthetized rodents ().

Together with shortened Te and lengthened Ti (Supplementary Figure 2), both hypercapnia and acetazolamide substantially shortened end-expiratory pause (EEP: the duration from the end of expiration to the beginning of inspiration) and elongated end-inspiratory pause (EIP: the duration from the end of inspiration to the beginning of expiration) (Figures 3F, G). Although both triggers of acidosis (acetazolamide and 5% CO2) affected EEP and EIP in the same direction, their effects on the neural timing of respiration could be differentiated, as shown in Figures 3F, G. Hypercapnia significantly shortened EEP for the entire 30 min and transiently prolonged EIP. For comparison, acetazolamide shortened EEP and increased EIP substantially, unlike the transient effects of hypercapnia on EIP. Additionally, these data revealed rapid neuronal responses to acetazolamide by respiratory kinetic adjustments; this has implications for the drug mechanisms underlying AMS.

Discussion

The whole-body plethysmography study examined the different impacts of erythroid AE1 with acidosis triggers (acetazolamide and hypercapnia) on respiration. We found that acetazolamide elicited earlier ventilatory responses in GPMur knock-in than WT mice. The lower tidal volume in GPMur mice (Figure 2; Supplementary Figure 1), consistent with GP.Mur+ people, may reflect more efficient CO2 excretion per breath due to enhanced erythroid AE1-mediated HCO3- transport. Ventilation with smaller tidal volumes potentially reduces the work of breathing. In human athletes, this is compensated by increased respiratory rates, resulting in higher overall minute ventilation (). By increasing blood CO2 and acidosis (verified by blood gas measurements), the two stimuli (hypercapnia and acetazolamide) substantially enlarged minute ventilation (). Notably, acetazolamide stimulation in both strains of the B6 mice manifested in persistent VT expansion and not in transient bursts of the respiratory rates (Figure 2; Supplementary Figure 1). This verifies others’ findings that acetazolamide predominantly expands tidal volume to increase minute ventilation ().

In these WBP experiments, the effects of GPMur/increased AE1 were opposite to the effects of acidosis by acetazolamide or hypercapnia (Table 2). Both acetazolamide and hypercapnia increased VT, minute ventilation ( E), inspiratory time (Ti), inspiratory duty cycle (Ti/Ttot), expiratory drive (VT/Te), estimated peak expiratory flow (PEFb), and end inspiratory pause (EIP); both also decreased expiratory time (Te), expiratory duty cycle (Te/Ttot), relaxation time (Tr), and end expiratory pause (EEP). In contrast, GPMur mice showed lowered tidal volume (VT) and distinct expiratory kinetics. The altered expiratory kinetics in GPMur mice likely reflects enhanced CO2 excretion rates at the alveolar–capillary interface. Faster CO2 unloading in pulmonary capillaries is a consequence of accelerated erythroid CO2/HCO3- conversion facilitated by AE1 transport; these alter sensing of central and peripheral chemoreceptors that control breathing patterns.

Respiratory control involves peripheral (RBC and carotid body chemoreceptors) and central (brainstem chemoreceptors) components. The present WBP study focuses on erythroid HCO3- transport. It has been reported that H+-sensing central chemoreceptors in the brainstem respiratory center contribute 70%–80% of the hypercapnic ventilatory responses (HVR), and peripheral carotid chemoreceptors contribute only 20%–30% HVR (). However, the contribution of RBC had not been clear. Here the faster ventilatory responses of GPMur mice demonstrated the RBC-AE1 contribution to respiration (opposite direction to that of acidosis). Thus, the effect of RBC-AE1 through CO2/HCO3- conversion should be viewed as an integral part of systemic respiratory control.

From the perspective of molecular mechanism, these findings support that AE1-mediated HCO3- transport, rather than CAII catalysis, is the rate-limiting step in erythroid AE1–CA functional coupling, despite ~22% less CAII in murine GPMur than WT RBCs. This interpretation is grounded in the hypothesis of “erythrocyte CO2-transport metabolon,” in which AE1, GPA/GPB (GPMur), Rh/RhAG, AQP1, deoxyhemoglobin (deoxyHb), and CAII are physically or indirectly associated with one another in close proximity on the RBC membrane (; ). The interaction between AE1 complexes and cytosol-abundant deoxyHb and CAII is transient. Erythroid AE1 and CAII share the substrate HCO3-, and yet erythroid CAII catalytic turnover (~106/s) far exceeds AE1 transport capacity (~104–105/s) (); therefore, AE1-mediated HCO3- membrane transport is the rate-limiting step in their functional coupling (). In this study, increased erythroid AE1 drove faster ventilatory responses despite reduced CAII in murine GPMur RBCs (Figure 1); this again points to AE1 (not CAII) as the rate-limiting factor in blood CO2 metabolism. Nonetheless, an alternative explanation should be considered. Because acetazolamide-mediated CA inhibition is systemic and beyond RBC, any strain/species differences in CA tissue distribution and acetazolamide pharmacokinetics might also contribute to the indirect effects of acetazolamide or compensatory mechanisms. We favor AE1 being rate-limiting because it is the only explanation reconciling the faster ventilatory responses in GPMur mice and GP.Mur+ individuals (; ).

Besides AE1 and CAII, Rhesus-associated glycoprotein (RhAG), AQP1, and UT-B on the erythrocyte membrane also participate in the transport of CO2 and H2O. RhAG is a gas channel that primarily transports NH3 over CO2 (; ). Erythroid AQP1 is responsible for ~85% of RBC osmotic water permeability and ~60% of RBC CO2 permeability (the latter likely through the dual pathways of AQP1) (; ; ). The urea transporter UT-B (SLC14A1; the Kidd/Jk antigen) mainly transports urea, but it also transports H2O and NH3; the physiological significance for the transport of the latter two substrates remains unclear (; ; ; ). Among these three erythroid transporters, the expression of UT-B (~14,000 copies per RBC) is less than 1/10 of the levels of RhAG (~100,000–200,000 copies) or AQP1 (~200,000 copies) (; ). RhAG is physically associated with tetrameric AE1 macrocomplexes (; ). AQP1 transiently interacts with AE1 on the RBC membrane, as revealed by cryo-EM and FLIM-FRET (; ). Conceivably, GP.Mur-enhanced AE1 expression may remodel the RBC membrane that, in principle, could alter the abundance or organization of other transporters/channels. Indeed we previously found fewer Rh/RhAG complexes in GP.Mur+/+ human RBCs (). Their different expressions and potential co-regulation on the RBC membrane are directions for future work.

Finally, the GPMur knock-in model recapitulates the direction of the human GP.Mur phenotype, namely: (1) human GPMur protein expression increases murine erythroid AE1 in mice and men (; ; ), (2) here the mice reproduce ventilatory patterns (lower tidal volume and altered expiratory kinetics) documented in GP.Mur+ people during exercise (; ), and (3) both GPMur+ mice and men have a higher blood pressure that is related to nitric oxide metabolism (; ). The first two warrant further translational–clinical investigation (e.g., blood gas assessment of acetazolamide dose responses) to verify whether GP.Mur+ individuals could achieve equivalent respiratory effects at lower doses.

Limitation

First, this is a murine WBP study using GPMur knock-in mice characteristic of ~55% higher erythroid AE1 expression. The magnitudes of the observed respiratory effects may not be translated directly to humans (10%–20% higher erythroid AE1in GP.Mur+ people). Second, murine RBCs lack CAI isoform that is present in human RBCs. Although the role of CAI in RBC is unclear, whether erythroid AE1–CA functional coupling could be affected by different CA isoforms is also unknown. Third, CA is ubiquitously expressed, while AE1 expression is restricted to the erythroid system.

It is thus unknown whether AE1 through circulating red blood cells could remotely influence carbonic anhydrases in other systems, such as the brainstem respiratory center, peripheral chemoreceptors, and/or kidneys. As acetazolamide systemically inhibits CA, the present WBP study design could not separate the drug effects on erythrocytes versus on other organs. Lastly, WBP inferred chemoreceptor and pulmonary responses from ventilatory-pressure measurements rather than direct measurements of chemoreceptor activity or alveolar–capillary CO2 exchange. Considering that respiratory mechanics and chemoreflex set points may vary between mice and men, how erythroid AE1 influences respiratory neuronal/chemoreceptor responses remains to be explored.

Conclusion

From this murine WBP study of conscious, unrestrained breathing, we found that increased RBC-AE1 expression (the GPMur phenotype) enhanced acetazolamide sensitivity in respiration. This implies that GP.Mur+ individuals may achieve therapeutic respiratory effects at lower acetazolamide doses and thus warrants future clinical investigation.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.

Ethics statement

The animal study was approved by MacKay Memorial Hospital IACUC. The study was conducted in accordance with the local legislation and institutional requirements.

Author contributions

C-TL: Conceptualization, Supervision, Writing – review & editing. J-HL: Formal analysis, Investigation, Writing – review & editing. TH: Formal analysis, Investigation, Writing – original draft. H-LL: Formal analysis, Investigation, Writing – review & editing. S-WW: Conceptualization, Writing – review & editing. KH: Conceptualization, Formal analysis, Funding acquisition, Methodology, Supervision, Validation, Writing – original draft, Writing – review & editing.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This work was funded by grants from Taiwan National Science and Technology Council (NSTC 113-2320-B-195-004; NSTC114-2320-B-195-003) and from MacKay Memorial Hospital to K.H. (MMH 113-67; MMH 113-95; MMH 114-25; MMH 115-030; MMH 115-090).

Acknowledgments

We thank P. L. Chen, L. Y. Chen, Taiwan National Center for Biomodels (NCB), and Taiwan Mouse Clinic for technical support in professional colony management and mouse experiments as well as Dr. Eri R. Swenson for critical reading of this manuscript.

Conflict of interest

The GPMur Knock-in mouse model is patented. The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2026.1908536/full#supplementary-material

References

Summary

Keywords

acetazolamide, anion exchanger-1 (AE1; band 3), carbonic anhydrase (CA), GP.Mur blood type (Mi.III), hypercapnia, whole body plethysmography (WBP)

Citation

Lai C-T, Lin J-H, Hu T, Lee H-L, Wang S-W and Hsu K (2026) Erythroid anion exchanger-1 (SLC4A1) modulates carbonic anhydrase inhibitor acetazolamide on respiration. Front. Physiol. 17:1908536. doi: 10.3389/fphys.2026.1908536

Received

13 June 2026

Revised

12 July 2026

Accepted

27 July 2026

Published

10 September 2026

Volume

17 - 2026

Edited by

Tim McMahon, Duke University, United States

Reviewed by

Doris Cerecedo, National Polytechnic Institute (IPN), Mexico

Saroj Tripathy, All India Institute of Medical Sciences, Deoghar (AIIMS Deoghar), India

Updates

Copyright

*Correspondence: Kate Hsu,

†ORCID: Kate Hsu, orcid.org/0000-0003-1009-0997

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics