Abstract

Genetic determinants confer variation among inbred mouse strains with respect to the magnitude and pattern of breathing during acute hypoxic challenge. Specifically, inheritance patterns derived from C3H/HeJ (C3) and C57BL/6J (B6) parental strains suggest that differences in hypoxic ventilatory response (HVR) are controlled by as few as two genes. The present study demonstrates that at least one genetic determinant is located on mouse chromosome 9. This genotype-phenotype association was established by phenotyping 52 B6C3F2 (F2) offspring for HVR characteristics. A genome-wide screen was performed using microsatellite DNA markers (n = 176) polymorphic between C3 and B6 mice. By computing log-likelihood values (LOD scores), linkage analysis compared marker genotypes with minute ventilation (V˙e), tidal volume (Vt), and mean inspiratory flow (Vt/Ti, where Ti is inspiratory time) during acute hypoxic challenge (inspired O2 fraction = 0.10, inspired CO2fraction = 0.03 in N2). A putative quantitative trait locus (QTL) positioned in the vicinity of D9Mit207 was significantly associated with hypoxic V˙e (LOD = 4.5), Vt (LOD = 4.0), and Vt/Ti (LOD = 5.1). For each of the three HVR characteristics, the putative QTL explained more than 30% of the phenotypic variation among F2 offspring. In conclusion, this genetic model of differential HVR characteristics demonstrates that a locus ∼33 centimorgans from the centromere on mouse chromosome 9 confers a substantial proportion of the variance inV˙e, Vt, and Vt/Tiduring acute hypoxic challenge.

C3H/HeJ

C57BL/6J

hypoxic hypoventilation

control of breathing

linkage analysis

hypoxic ventilatory sensitivity varies among human subjects and represents a continuum from low to high responsive subgroups (2, 7, 10, 11,34). Hirshman et al. (7) demonstrated a five- to sixfold difference between extreme ventilatory responses in a group of 44 subjects. Although phenotypes in these subjects varied positively with height, weight, and the ventilatory response to hypercapnia, the bimodal attribute of the distribution suggested that a major genetic determinant underlies individual variation in hypoxic ventilatory sensitivity.

Our laboratory has demonstrated a similar continuum among many inbred mouse strains (32). This suggested that the genetic diversity, which regulates variation in hypoxic ventilation, is conserved across species, representing the basis for comparative mapping studies between mouse and human genomes (27). The C3H/HeJ (C3) and C57BL/6J (B6) inbred strains have been the focus of our research because these strains exhibit breathing differences during a broad array of inspirate challenges, including acute hypoxia. Although C3 and B6 mice were comparable with respect to their hypoxic minute ventilation (V˙e), the strains differed in breathing pattern; i.e., C3 mice were characterized by a slow, deep hypoxic ventilatory response (HVR) relative to the rapid, shallow phenotype of B6 mice. Moreover, the HVR profile of first-generation offspring (i.e., B6C3F1/J) of C3 and B6 progenitors represented a third phenotype distinguished by hypoventilation, which was attributable to a marked attenuation in hypoxic tidal volume (Vt). By using quantitative genetic approaches, our studies further showed that variation in HVR characteristics, includingV˙e and Vt, was inherited by second-generation offspring in segregation patterns consistent with two-gene models (31).

The purpose of the present study is to link differential HVR characteristics between C3 and B6 mice to candidate regions of the mouse genome. To achieve this objective, a genome-wide screen was performed using 176 microsatellite markers to analyze DNA samples from 52 B6C3F2 (F2) offspring. The F2offspring were previously characterized with respect to eupneic, hypoxic, and hypercapnic ventilatory responses; a quantitative genetic analysis of these results are reported elsewhere (31). The microsatellite markers represented simple sequence repeat fragments polymorphic between C3 and B6 mice and were evenly distributed across the mouse genome. The hypothesis of the present study suggests that phenotypic differences in hypoxic V˙e and Vt among the parental strains and the B6C3F1/J (F1) offspring are regulated by major genetic determinants. The results demonstrate that a candidate genomic region on mouse chromosome 9 explains a substantial proportion of the phenotypic variation in hypoxic V˙e and Vt and therefore represents a significant quantitative trait locus (QTL) for differential HVR in this genetic model.

METHODS

Animals.

Male and female F1 mice were purchased from Jackson Laboratory (Bar Harbor, ME) to establish breeding colonies at the Johns Hopkins School of Public Health. Intercross progeny were generated from F1 progenitors and weaned at 4–5 wk of age. From the breeding colonies, randomly selected male F2 offspring (n = 69) were housed in cages of four to six animals for an additional 6–12 wk. Water and mouse chow (Agway Pro-Lab RMH 1000) were provided ad libitum. All animal protocols were reviewed and approved by the Animal Care and Use Committee of the Johns Hopkins School of Public Health.

HVR.

The magnitude and pattern of ventilation were determined by whole body plethysmography that used unrestrained and unanesthetized conditions. The methods and results have been described in detail and are reported elsewhere (31, 32). Briefly, ventilatory function was evaluated during a sequence of acute (3–5 min) inspiratory challenges including hypoxic (inspired O2 fraction = 0.10) admixtures in mild hypercapnia (inspired CO2fraction = 0.03). During hypoxia, the V˙eresponses were comparable between C3 and B6 parental strains but were significantly (P < 0.01) lower in F1 mice. There was a consistent strain difference in breathing frequency (f) between C3 and B6 mice during eupneic, hypoxic, and hypercapnic conditions, and F1 mice resembled the B6 progenitor. In contrast, Vt responses were greatest in C3 mice, and the F1 offspring demonstrated markedly lower Vtresponses during hypoxic conditions compared with B6 mice. These observations, among others, compelled us to characterize the HVR profile of F1 mice as a third phenotype representing hypoxic hypoventilation.

In the present study, cosegregation analysis is used to establish the strength of association between the measurement of hypoxicV˙e and its components. Cosegregation analysis is a powerful genetic tool that resembles aspects of linear regression and correlational analyses by comparing the covariation of two variables. In cosegregation analysis, covariation of two variables suggests that the genes, which regulate phenotypic variation for one variable, are the same as those that regulate variation in the other variable. Although cause and effect relationships cannot be inferred from cosegregation analysis, phenotypes that covary in F2 mice likely originate from the same or a closely linked gene with high probability.

DNA isolation.

After measurements of HVR were made, DNA samples were isolated from 52 randomly selected F2 mice. Each animal was euthanized with pentobarbital sodium to obtain kidney tissue. To isolate high-molecular-weight DNA (22), sample tissue was gently homogenized using a Dounce homogenizer in an isotonic, high-pH solution that included Nonidet P-40. To release genomic DNA, nuclear membranes were lysed in a sodium dodecyl sulfate-proteinase K solution. After incubation at 42°C for 3 h, DNA was separated from other cellular components by extraction using a phenol-chloroform-isoamyl alcohol solution. DNA was precipitated with ice-cold 95% ethanol and 0.15 M potassium chloride. Sample DNA was spooled on a glass rod, rinsed with 70% ethanol, and stored in Tris-EDTA buffer (pH = 8.0). The purity of the sample DNA was determined using spectrophotometry and by calculating the ratio of optical densities (OD) at wavelengths of 260 and 280 nm (i.e., a target OD260-to-OD280 ratio of 1.8).

Acrylamide (6%) denaturing gels were used to separate single strands of amplified DNA according to size. Appropriate PCR controls (i.e., B6, C3, and F1 DNA and no DNA) were used to ensure proper sample loading and to establish the absence of DNA contamination. The genotypes of the F2 offspring were determined by comparing the allelic size of the simple sequence repeat to the known DNA allele size of the parental strains.

Marker selection.

At the time of this study, ∼500 PCR-based markers were available throughout the mouse genome that differed between C3 and B6 strains. Markers were selected to establish 10- to 15-centimorgan (cM) intervals between two loci and to provide coverage of the entire mouse genome with 95% confidence. For F2 offspring, a sample size of 50 mice theoretically excludes 10–15 cM of the genome for each marker using a set of ∼150 well-spaced markers (i.e., given the mouse genome is ∼1,500 cM long). Although this estimate was used in the present study, three to six additional markers were used to saturate genomic regions showing suggestive linkage [log-likelihood values (LOD) > 2.8]. Therefore, the linkage analysis in the present study incorporated ∼95% of the mouse genome with marker intervals of <15 cM. In genomic regions showing suggestive linkage, the size of the gene-containing region was reduced to 1- to 3-cM intervals.

Linkage analysis.

Linkage analysis was performed with the use of the computer software programs Mapmaker-EXP and Mapmaker-QTL (13, 14). The distances between loci map assignments were established using Mapmaker-EXP. Because the order of DNA markers was generally predetermined, the results obtained from Mapmaker-EXP concerning distances between loci were compared with known distances defined for each marker. Interval mapping of phenotypic data was accomplished using Mapmaker-QTL. This analysis incorporated high-density genetic maps to distinguish between plausible QTL and weak effects from distant loci. Subroutines and algorithms for Mapmaker-QTL were based on assumptions that quantitative traits are normally distributed. Therefore, power-transformed data sets were shared between Mapmaker-EXP and Mapmaker-QTL to ascertain linkage and compute LOD scores. Linkage was inferred when the log of the odds ratio (i.e., called a LOD score) in favor of linkage over nonlinkage was 1,000 to 1 (i.e., LOD > 3.3). Alternatively, nonlinkage was deduced when the odds against linkage were 100 to 1 (i.e., LOD < 1.8). The candidate region surrounding putative QTL was examined for genes with biological relevance to hypoxic ventilatory control. Finally, comparative mapping was performed to determine homologies between the mouse and human genomes.

Supplementary data analysis.

Correlation and linear regression analyses were performed to estimate the strength of association between hypoxic V˙e and its components. The paired components of V˙e are either Vt and f or mean inspiratory flow (Vt/Ti, where Ti is inspiratory time) and inspiratory fraction (Ti/Ttot, where Ttot is the total respiratory time). In the latter case, Vt/Ti and Ti/Ttot components are used to identify the “driver” and “timer” elements (18), respectively, as a general model of respiratory control mechanisms. Although the paired components are obvious mathematical correlates of V˙e, dissociation between components acts to accentuate one trait to explain the greatest proportion of genetic variance among individuals.

One- and two-way ANOVAs are used to evaluate the effects of DNA marker loci and to assess the interaction between two genes. A two-tailedt-test was used to compare means between homozygous and heterozygous QTL forms. A 99% confidence level was used to establish statistical significance.

RESULTS

In Fig. 1, cosegregation analysis was used to examine the relationship between hypoxicV˙e responses and either Vt or f components of an individual response. As suggested by correlation coefficients, both Vt and f demonstrated significant (P < 0.0001) positive associations with hypoxicV˙e. The Vt component accounted for ∼72% of the variance in hypoxic V˙e, whereas f explained ∼49% of the variance.

In Fig. 2, cosegregation was used to examine the relationship between hypoxic V˙e responses as a function of both Vt/Ti and Ti/Ttot components. The correlation coefficients suggested that only Vt/Ti responses were significantly (P < 0.0001) associated with hypoxicV˙e; that is, there was relatively little association (P > 0.01) between hypoxic Ti/Ttot andV˙e responses. In this case, the Vt/Ti component accounted for ∼92% of the variance in hypoxic V˙e.

As shown in Fig. 3, LOD scores were plotted as a function of the relative marker distance for hypoxicV˙e, Vt, and Vt/Tiresponses. On the basis of peak LOD scores of 4.5, 4.0, and 5.1 for hypoxic V˙e, Vt, and Vt/Ti responses, respectively, a putative QTL was established for a candidate region in close proximity toD9Mit207 on mouse chromosome 9. For each of the three HVR characteristics, the putative QTL explained ∼30% or more of the phenotypic variation among F2 offspring.

In Fig. 4, individual responses for F2 mice were plotted for homozygous and heterozygous genotypes at the D9Mit207 locus of chromosome 9. The allelic frequencies were not significantly different from Mendelian proportions (χdf=22 = 1.04; P > 0.05). For each HVR characteristic, the response distributions of the heterozygotes were significantly (P < 0.01) lower compared with the two other homozygous groups. In addition, response distributions of the homozygotes with alleles derived from the C3 progenitor spanned the distributions of the other two groups.

Hypoxic V˙e (left), Vt(middle), and Vt/Ti(right) responses for 52 F2 mice are depicted as a function of genotype at the D9Mit207 locus. For each hypoxic response, the heterozygotes (shaded symbols) are significantly (P < 0.01) lower compared with the other two groups. In addition, the variation for C3-like homozygotes (open symbols) spans the range of responses extended by the other two groups. Means ± SE for each genotypic group (i.e., B, B6-like homozygotes; H, heterozygotes; A, C3-like homozygotes) are shown to the right. *P < 0.01.

In Table 1, peak LOD scores were reported for the corresponding hypoxic ventilatory traits. The genome-wide screen produced one other putative QTL (LOD score = 3.3) on mouse chromosome 3 in the vicinity of D3Mit17 for hypoxic Vt/Ti responses. This QTL was shown to account for ∼26% of the variance in the hypoxic Vt/Ti responses among the F2 mice. Two other suggestive QTL (LOD scores = 2.8) were shown to link HVR response variation to mouse chromosomes 2 and 6.

In Fig. 5, the allelic interaction between D9Mit207 and D3Mit17 loci was illustrated for the hypoxic Vt/Ti responses in F2 mice. The allelic frequencies were not significantly different from Mendelian proportions (χdf=82= 7.24; P > 0.05). In addition, a two-way ANOVA suggested that there was a significant effect of bothD9Mit207 (F statisticdf=2 = 10.36; P < 0.001) and D3Mit17 (Fstatisticdf=2 = 6.92; P < 0.01) loci in determining variation in hypoxic Vt/Tiresponses among F2 mice.

Hypoxic Vt/Ti responses are depicted as a function of genotypes at two loci represented byD9Mit207 and D3Mit17 markers. TheD3Mit17 genotype interacts with the D9Mit207locus and influences hypoxic Vt/Ti responses in C3-like homozygotes at the D9Mit207 locus.

DISCUSSION

The present study demonstrates a putative QTL on mouse chromosome 9 that confers variation in HVR responses among F2 mice derived from C3 and B6 progenitors. Collectively, the evidence suggests that the QTL influences hypoxic V˙e by regulating the Vt component of the response (Fig. 1, left). In addition, individual variation in hypoxic Vt/Tiresponses tightly cosegregates with hypoxic V˙eresponses (Fig. 2, left). Furthermore, Vt/Ti emerges to be the HVR characteristic that demonstrates the strongest linkage relationship with a genomic region proximal to D9Mit207, which is a marker ∼33 cM from the centromere on mouse chromosome 9 (Fig. 3). This genotype-phenotype association suggests that the QTL regulates individual variation in hypoxic V˙e by altering a mechanism identified with the neural drive component of breathing (18). In contrast, the respiratory timing components of hypoxic V˙e, including f, Ti, expiratory time, and Ti/Ttot, are not linked to the QTL on mouse chromosome 9 in this segregant offspring class (Table 1).

Although the C3 and B6 parental strains are comparable with respect to hypoxic V˙e, the responses of F1 mice are significantly attenuated relative to both progenitors due to reduced hypoxic Vt and Vt/Ti responses (31). One hypothesis, which may explain the variation in HVR, is supported by two major strain differences between C3 and B6 parental strains: 1) a reduced lung volume and compliance in the B6 progenitor (33) and 2) a muted ventilatory response to chemical inspirate challenge in the C3 progenitor (32). The F1 mice inherit a phenotypic profile that includes a smaller lung volume from the B6 progenitor and a muted chemical response from the C3 progenitor. Therefore, one testable hypothesis suggests that the attenuated HVR of F1 mice is attributable to two major genetic determinants, which are dominant for reduced lung size and diminished chemical responsiveness.

Inheritance studies suggest that the attenuated hypoxic Vtand Vt/Ti response phenotypes of F1mice are conferred to second-generation offspring (31). In addition, similar phenotypic profiles are detectable in selected recombinant inbred strains derived from C3 and B6 progenitors (29). Collectively, these results support the hypothesis that a small number of heritable traits culminate to produce a phenotypic profile representing hypoxic hypoventilation. Recent evidence combining backcross and intercross offspring of C3 and B6 strains suggest that the inheritance patterns for HVR variation are consistent with models that incorporate two unlinked genes (31). The results from the present study demonstrate that, in addition to a locus on mouse chromosome 9, there is another putative QTL on mouse chromosome 3 linked to variation in hypoxic Vt/Ti in this model. This QTL is positioned in a genomic region adjacent to D3Mit17, which is a locus ∼72 cM from the centromere. The D3Mit17 locus differs from the putative QTL assigned to mouse chromosome 3, Itbq1 (located ∼25 cM from the centromere), which determines differences in inspiratory timing at baseline (30). Indeed, our laboratory has shown very little association between hypoxic and baseline ventilatory responses in this genetic model (29,31).

The segregation of hypoxic V˙e, Vt, and Vt/Ti responses of F2 mice based on their genotype at the D9Mit207 locus does not appear to be consistent with dominant or codominant genetic models (15). Although the allelic frequencies are consistent with Mendelian proportions, the average HVR in the heterozygotes are significantly attenuated compared with the average HVR responses of the homozygotes (Fig. 4). Likewise, the variation of the C3-like homozygotes spans the range of the other two groups. It is difficult to explain this segregation pattern based on Mendelian principles. A study by Kingsley et al. (12) constructed a similar linkage map of mouse chromosome 9 using the B6 progenitor. These investigators concluded that there was no evidence of segregation distortion or chromosome rearrangement in the candidate genomic region described in the present study. Collectively, these observations support the hypothesis that a second influential genetic determinant may be interacting with the gene at the D9Mit207 locus. Alternatively, the results are consistent with an underdominance genetic model (36, 37).

To test whether a second gene is interacting with theD9Mit207 locus, hypoxic Vt/Tiresponses are reclassified on the basis of genotypes at both theD9Mit207 and D3Mit17 loci. In Fig. 5, the variation in hypoxic Vt/Ti responses is shown to be independent of the genotype at the D3Mit17 locus for B6-like homozygous and heterozygous genotypes at the D9Mit17locus; that is, the B6-like homozygotes are consistently high, whereas the heterozygotes are consistently low. However, the D3Mit17genotype interacts with the D9Mit207 locus and influences hypoxic Vt/Ti responses in C3-like homozygotes at the D9Mit207 locus. The gene-gene interaction is manifested in B6-like homozygotes at the D3Mit17 locus to produce an attenuated hypoxic Vt/Ti response. Although larger numbers of F2 mice are essential to discern the significance of the potential interaction between the two QTL, the present data suggest that a gene on chromosome 3 modifies the genetic effects on HVR assigned to the QTL on chromosome 9.

The LOD plots in Fig. 3 appear to incorporate a secondary peak associated with D9Mit2. Although a two-QTL model is rejected for the genomic region between D9Mit2 andD9Mit208 using Mapmaker-QTL, this broader genomic region cannot be excluded for future consideration. In addition to ∼100 known genes in the 19-cM interval between D9Mit2 andD9Mit208 (17 and 36 cM from the centromere, respectively), there are likely many more novel genes in this candidate genomic region. The current working hypothesis considers genes with biological relevance to differential hypoxic Vt/Ti as likely candidates. More specifically, we postulate that the candidate genes in this region influence HVR variation by modifying the Vt response achieved during acute hypoxic exposure.

The most obvious candidate gene in the genomic region of mouse chromosome 9 is the Drd2 gene, which is located ∼28 cM from the centromere and encodes the dopamine D2 receptor. The effect of D2 receptors on hypoxic ventilatory control has been investigated in many different species, including mice. For example, Olson and Saunders (19) demonstrated that intraperitoneal administration of dopamine and levodopa at relatively high doses caused a depression in acute HVR. In a subsequent study (20), these investigators demonstrated an augmented acute HVR in mice after administration of a dopamine antagonist, droperidol. In general, these results are consistent with the acute effects specific to D2-receptor antagonists such as domperidone. That is, D2-receptor blockade releases the inhibitory effects of dopamine acting via D2 receptors, resulting in an increase in acute hypoxic ventilation. The results of Iturriaga et al. (9) suggested that the release of this modulatory role of dopamine, acting through D2 receptors, resulted from increased excitation of carotid chemosensory activity.

More recently, Huey et al. (8) investigated acute HVR in D2-receptor-deficient homozygous mice on a B6 and 129-hybrid genetic background. Although there were substantial gender effects, these investigators concluded that hypoxic and/or hypercapnic ventilation was significantly greater in D2-deficient mice compared with wild-type littermates. In many cases, genotypic differences in ventilation were derived from a significantly greater Vt response in D2-deficient mice. The consistency between the study of Huey et al. and the results of this study suggest that a DNA polymorphism in the Drd2 gene exists between C3 and B6 mice. This possibility is further supported by results of Cotzias and colleagues (3, 28), showing strain-dependent differences in neurological responses to levodopa between C3 and B6 strains.

Another candidate gene more proximal to the D9Mit207 locus on chromosome 9 is the Chrna5 gene, which is located ∼32 cM from the centromere (4, 9a). Other members of the nicotinic acetylcholine receptor gene family are also located on chromosome 9, including the Chrna3 and Chrnb4 genes. Collectively, the genes encode neuronal nicotinic acetylcholine receptor subunits, which are widely distributed throughout the central and peripheral nervous system. Nicotinic acetylcholine receptor subunits form pentameric channels resulting in heterogeneous expression and functional diversity (5, 23, 35). Pharmacological studies suggested that C3 and B6 mice responded differently to acute nicotine administration in a dose-dependent fashion (1,17). Although the strains did not differ in nicotinic or α-bungartoxin binding in many regions of the central nervous system (16), a dose-dependent increase in ventilatory response to nicotine was greater in B6 compared with C3 mice (17). Whether DNA variation in nicotinic acetylcholine receptor subunits, such as Chrna5 or others on chromosome 9, plays a role in differential HVR between C3 and B6 strains requires future studies.

Considering the observation that the heterozygotes demonstrate an attenuated hypoxic Vt/Ti response relative to both parental strains is not easily understood using simple Mendelian models. One speculative explanation with respect to the present linkage results suggests that underdominance of the heterozygotes is attributable to parental DNA variation in subunit composition of heteromeric channels. Functional deficiencies may emerge, for example, if DNA variation in the Chrna5 gene regulates differential heterologous expression of the α5-subunit. Recent studies suggest that α5 substitution for α7 or in combination with β-subunits produces dysfunctional nicotinic acetylcholine receptor complexes (5, 35). Therefore, a cogitative explanation for underdominance linked to the D9Mit207 locus suggests that variation in the Chrna5 gene modifies heteromeric composition resulting in an amplification of dysfunctional nicotinic acetylcholine receptors. The observation that α7-subunits have been shown in nerve fibers of the carotid body (26) provides neuroanatomic support for this provisional hypothesis.

Perspective.

In a 1992 review article honoring Pierre Dejours, Shea and Guz (25) proposed that individual variation in the magnitude and pattern of breathing is likely attributable to the interaction between lung mechanics and chemosensitivity. An individual outcome results in a breathing strategy that optimizes the work of breathing. With acute hypoxic challenge, individual variation in chemosensitivity may emerge to overcome the limits imposed by variation in lung mechanics, and optimal work of breathing may be restored. Powell and colleagues (21) outline specific mechanisms associated with the time course of hypoxic ventilation. These investigators proposed at least three distinct mechanisms operating within the first several minutes to influence acute HVR. The most immediate response is an augmentation of ventilation mediated by peripheral chemoreception and glutamanergic mechanisms in the nucleus tractus solitarius. Short-term potentiation and depression are mechanisms that follow the acute response. The latter two stages involve other respiratory control centers, including respiratory rhythm generation, which are not specific to peripheral chemoreception. Although it is impossible to delineate between the mechanisms described above using our experimental protocol, it is clear that C3 and B6 parental strains vary with respect to respiratory timing, lung mechanics, and ventilatory responsiveness to chemical challenge.

With respect to the present results, a DNA variant at theD9Mit207 locus may be influencing chemosensitivity or central integration, whereas a second DNA variant at theD3Mit17 locus influences lung mechanics. In the B6-like homozygotes at the D9Mit207 locus, HVR is consistently high due to a DNA variant that determines a relatively robust chemical response (see Fig. 5). Likewise, the heterozygotes at theD9Mit207 locus demonstrate consistently low HVR attributable to the dominant effect (i.e., muted chemosensitivity) of an allele derived from the C3 parental strain. The interaction of theD3Mit17 locus (or a gene that regulates lung mechanical differences) emerges to influence the variation in C3-like homozygotes at the D9Mit207 locus. That is, there is an inherent mechanical advantage in individuals with greater lung volumes. The positive phenotypic association between hypoxic ventilatory sensitivity and height observed by Hirshman et al. (7) suggests that interactions between lung mechanics and chemosensitivity affect individual variation in human subjects. Indeed, lung volume is directly related to the cube of body height in humans (24).

In conclusion, a QTL on mouse chromosome 9 determines a substantial proportion of the variance in acute HVR among F2 offspring derived from C3 and B6 parental strains. The genomic region surrounding the D9Mit207 locus incorporates DNA polymorphisms, which regulate differences in hypoxic Vt and Vt/Ti responses between C3, B6, and F1 mice. Phenotypic segregation based on this QTL appears to support an underdominant genetic model in which the heterozygotes consistently demonstrate attenuated HVR. However, our working hypothesis proposes that the attenuated hypoxic Vt and Vt/Ti responses in individual F2mice is regulated by one or more interactive genetic determinants, including a QTL on chromosome 3. Comparative mapping studies suggest that the genomic region encompassing D9Mit207 in the mouse genome is homologous to regions on chromosomes 11 and 15 of the human genome (5). For example, the human DRD2 gene is mapped to 11q22, whereas CHRNA5 is mapped to 15q24 (9a). Future studies are needed, incorporating gene expression and DNA sequence analysis, that search for single nucleotide polymorphisms, which determine individual variation in ventilatory control during hypoxia.

Acknowledgments

The author gratefully acknowledges the technical assistance of Heather Kulaga and the cooperation of the Department of Anesthesiology at the Johns Hopkins School of Medicine.

Footnotes

This study was supported by National Heart, Lung, and Blood Institute Grant HL-53700.

Address for reprint requests and other correspondence: C. G. Tankersley, Division of Physiology, School of Hygiene and Public Health, The Johns Hopkins Univ., 615 N. Wolfe St., Baltimore, MD 21205.

The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked “advertisement” in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.