Skip Navigation

Integrative and Comparative Biology 2005 45(3):405-415; doi:10.1093/icb/45.3.405
This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (4)
Right arrow Request Permissions
Google Scholar
Right arrow Articles by Koch, L. G.
Right arrow Articles by Britton, S. L.
Right arrow Search for Related Content
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

The Society for Integrative and Comparative Biology

Divergent Selection for Aerobic Capacity in Rats as a Model for Complex Disease1

Lauren Gerard Koch2,1 and Steven Loyal Britton1
1 Functional Genomics Laboratory, Department of Physical Medicine and Rehabilitation, University of Michigan, 1500 E. Medical Center Dr., Ann Arbor, Michigan 48109-0549


    SYNOPSIS
 TOP
 SYNOPSIS
 INTRODUCTION
 INTRODUCTION
 MODELS OF COMPLEX DISEASE
 EVOLUTION-BASED SELECTION:...
 SELECTION FOR AEROBIC CAPACITY
 USE OF N:NIH STOCK...
 DESIGN OF SELECTION EXPERIMENTS:...
 CURRENT DATA ON RATS...
 References
 
Based upon ideas about evolution, we put forth the argument that the capacity to transfer energy via aerobic metabolism is such a central feature of mammalian biology, that it must also be the primary determinant of complex disease. From this, we hypothesized that artificial selection on low and high capacity for aerobic exercise would create lines that can be used to define the divide between health and disease. In 1996 we began large-scale divergent selection for aerobic treadmill running capacity in a widely heterogeneous stock of rats (N:NIH). By ten generations we developed lines of low capacity runners (LCR) and high capacity runners (HCR) that on average differed by 317%. As a correlated trait, body mass increased at each generation in the LCR while the body mass decreased in the HCR. The lines also separated for key factors of systemic oxygen transport capacity such as maximal oxygen consumption (VO2max), tissue perfusion, capillary density, and oxidative enzyme activity (citrate synthase and B-HAD). We also tested our hypothesis that differences in aerobic energy transfer would produce rats that contrast for risk factors associated with complex disease. Indeed, the lines separated for cardiovascular risk factors including differences in blood pressure, cardiac contractility, visceral adiposity, plasma free fatty acids, and triglycerides. The decrease in aerobic capacity was also associated with low amounts of several proteins required for mitochondrial function.


    INTRODUCTION
 TOP
 SYNOPSIS
 INTRODUCTION
 INTRODUCTION
 MODELS OF COMPLEX DISEASE
 EVOLUTION-BASED SELECTION:...
 SELECTION FOR AEROBIC CAPACITY
 USE OF N:NIH STOCK...
 DESIGN OF SELECTION EXPERIMENTS:...
 CURRENT DATA ON RATS...
 References
 

"Well, complex diseases are really a new frontier for genetics. ...... just like every other area of science we better build a very firm foundation. ......It is a very hard task. Single gene diseases have been worked on in some sense for the entire century. Multi-gene diseases are very much the work of the last 5 years or so. And so complex inheritance is a frontier for the next century."

Eric S. Lander, DPhil (excerpt from an interview that took place at the "Winding Your Way through DNA" symposium at the University of California San Francisco in 1992) (http://www.accessexcellence.org/RC/CC/lander.html.)


    INTRODUCTION
 TOP
 SYNOPSIS
 INTRODUCTION
 INTRODUCTION
 MODELS OF COMPLEX DISEASE
 EVOLUTION-BASED SELECTION:...
 SELECTION FOR AEROBIC CAPACITY
 USE OF N:NIH STOCK...
 DESIGN OF SELECTION EXPERIMENTS:...
 CURRENT DATA ON RATS...
 References
 
Most common diseases are complex in the sense of being determined by allelic variations in possibly few to hundreds of genes (polygenic) as they interact with variable environments. Diabetes, hypertension, obesity, and coronary artery disease are examples of complex disease and are the major cause of morbidity and mortality, at least in westernized societies. About ten years ago we began selection experiments based upon ideas about evolution to create rat genetic models that divide for low and high health risks. Here we present background information on our approach to selection, the ideas we borrowed from evolution to formulate our hypothesis, and current information about our models.


    MODELS OF COMPLEX DISEASE
 TOP
 SYNOPSIS
 INTRODUCTION
 INTRODUCTION
 MODELS OF COMPLEX DISEASE
 EVOLUTION-BASED SELECTION:...
 SELECTION FOR AEROBIC CAPACITY
 USE OF N:NIH STOCK...
 DESIGN OF SELECTION EXPERIMENTS:...
 CURRENT DATA ON RATS...
 References
 
Current model approaches are less than an ideal representation of complex diseases. Physical maneuvers such as ligation of coronary arteries to model heart disease (Zolotareva and Kogan, 1978Go; Ahn et al., 2004Go) or renal artery occlusion to emulate renovascular hypertension (Goldblatt et al., 1934Go; Pinet et al., 2004Go) more accurately reflect response to injury. Administration of chemical agents to mimic diseases such as streptozotocin-induced diabetes mellitus (Rossini et al., 1977Go) do not reflect the progression of disease, but represent response to chemical injury. Single gene knockout approaches are useful to determine essentiality of a gene's function, but do not represent the polygene condition of complex diseases (Williams and Wagner, 2000Go). Application of gametic mutagens to produce random base changes with the hope that one or more of the mutants will resemble a disease phenotype does not define the combination of alleles that are causative of complex disease. Direct selection for a disease phenotype is also problematic because it is based upon measurable traits or symptoms and not the full complement of underlying mechanisms.

For instance, hypertension is a well-characterized quantitative disease trait, by which selection models for variation in the measure of blood pressure itself seems reasonable (Rapp, 2000Go). Figure 1 panel A demonstrates a hypothetical selection response for extremes in the measure of blood pressure to yield lines low and high for mean blood pressure. Assuming that an average mean arterial pressure is 90 mmHg, selection on the extremes could produce divergent lines that on average demonstrate a mean arterial pressure of 140 mmHg for the high selected line and 75 mmHg for the low selected line after six generations. This approach can be described as "trait modeling" in that the models can then be used to explore mechanism for the trait. The problem is that selection based upon a single measure thought to characterize the disease will not necessarily reproduce the full array of underlying mechanisms and present the full complexity of the disease. This problem is amplified because diseases emerge not as discrete events, but as complexes, such as the cascade represented by metabolic syndrome (Gensini et al., 1998Go) for which hypertension, diabetes, heart failure, and obesity tend to cluster.



View larger version (18K):
[in this window]
[in a new window]
 
FIG. 1. Theoretical data for the development of a selection model to study a complex disease (i.e., hypertension). Panel A demonstrates selection response for measures of blood pressure to produce lines that differ in trait to investigate mechanisms for low and high blood pressure (trait modeling). Whereas, Panel B shows divergent selection response for pressure-diuresis, an assumed mechanism of hypertension to produce lines that differ for trait of blood pressure and other disease-associated traits (mechanistic modeling)

 
From this account, it seems reasonable to build on the idea of "mechanistic modeling" in which selection is based on the mechanism underlying the disease of interest. For example, Guyton and Coleman (1967)Go contend that long-term blood pressure regulation, and thus hypertension, is mediated by the kidney via a process termed pressure diuresis. The basic tenet of pressure diuresis is that an organism has just enough blood pressure to keep the input of salt and water equal to the output (Guyton, 1991Go). If pressure diuresis is indeed the basic mechanism of hypertension, then selection based upon low function for pressure diuresis should yield accurate models that correctly emulate the condition of hypertension. A hypothetical example of such selection is illustrated in Figure 1 panel B. Here, a pressure-diuresis index (PDI) is estimated as the product of mean arterial blood pressure and time (hr) required to excrete one-half of a 4 ml oral NaCl challenge (0.5 mequiv./ml). Presume that a founder population accomplished one-half excretion with an average mean arterial pressure equal to 100 mmHg over a 3 hour period (PDI = 300). After six generations of selection, animals with low pressure-diuresis took 6 hours to excrete the same NaCl load with an average blood pressure of 140 mmHg (PDI = 840). The animals bred for high pressure-diuresis were able to excrete the same NaCl challenge with an average blood pressure of 100 mmHg within 1 hour (PDI = 100). Thus, the high line represents a model for the trait of hypertension as well as having the underlying mechanism of low-pressure diuresis. In contrast, the low line represents a model with a high pressure-diuresis mechanism and low blood pressure. Although conceptually correct (Brand et al., 1991Go; Skarlatos et al., 1994Go; Steele et al., 2000Go), our efforts to yield a selection measure of pressure-diuresis that could be taken to large scale breeding have not been successful. Mechanistic modeling is difficult because, unlike pressure-diuresis for hypertension, not many complex disease examples exist for which an underlying mechanism has been resolved.


    EVOLUTION-BASED SELECTION: ENERGY TRANSFER DEFINES OUR BIOLOGICAL EXISTENCE
 TOP
 SYNOPSIS
 INTRODUCTION
 INTRODUCTION
 MODELS OF COMPLEX DISEASE
 EVOLUTION-BASED SELECTION:...
 SELECTION FOR AEROBIC CAPACITY
 USE OF N:NIH STOCK...
 DESIGN OF SELECTION EXPERIMENTS:...
 CURRENT DATA ON RATS...
 References
 
If evolution determines our current biology, then it seems reasonable to assume that mechanisms of disease originated from our evolutionary history. Ideas about evolution are of course intertwined with consideration of the driving force of life's emergence from the inanimate to the animate (Pross, 2003Go) and of the environments from which natural selection occurred. We reasoned that energy transfer via aerobic metabolism is the most formative feature that typifies mammalian existence at all levels of biological organization; and thus, must also be a primary determinant of disease.

Two billion years of evolution in an oxygen environment resulted in aerobic metabolism as a central feature of mammalian biology (DesMarias, 2000Go; Paytan, 2000Go). Before an oxygen-enriched atmosphere, an organism's capacity for energy transfer was limited by operating within anaerobic pathways to support metabolism and replication. The introduction of oxygen as a final electron acceptor, created an 18-fold greater capacity to transfer energy. It is presumed that this increased free energy transfer afforded by the widened redox potential was permissive for more complexity in the sense that organisms became multicellular (Xiong et al., 2000Go). Obligatory for using oxygen in energy transfer pathways was the co-evolution of enzymes that detoxify the reactive oxygen species formed as by-products (Young and Woodside, 2001Go). Thus, the genetic substrates that mediate oxidation reactions and oxygen detoxification reactions likely form a large foundation for mammalian biology. This is consistent with the view that most diseases are caused by flaws in oxidation/detoxification pathways or defects in mitochondrial structure and function (Young and Woodside, 2001Go; Duchen, 2004Go). This supports the use of maximal oxygen consumption (VO2max) as a clinical reference point and total body aerobic capacity as a predictor of all-cause morbidity and mortality in humans (Myers et al., 2002Go).


    SELECTION FOR AEROBIC CAPACITY
 TOP
 SYNOPSIS
 INTRODUCTION
 INTRODUCTION
 MODELS OF COMPLEX DISEASE
 EVOLUTION-BASED SELECTION:...
 SELECTION FOR AEROBIC CAPACITY
 USE OF N:NIH STOCK...
 DESIGN OF SELECTION EXPERIMENTS:...
 CURRENT DATA ON RATS...
 References
 
The rat was chosen as the model organism for selection for three primary reasons 1) the laboratory rat (Rattus norvegicus) has been a key tool for the study of medicine and pharmacology for human health. A large database of physiological phenotypes for integrated fields such as cardiovascular, neuroscience, and exercise physiology exist in the literature (Jacob and Kwitek, 2002Go). 2) To date, almost all known genes associated with disease in humans have orthologues to the rat genome (Rat genome sequencing project consortium, 2004Go). 3) Over 200 inbred strains of R. norvegicus have been produced from lines and colonies of rats first selectively bred for "disease" alleles (Greenhouse et al., 1990Go). Recent evidence from human studies suggests two genetic substrates as contributors to the aerobic energy capacity phenotype (Bouchard et al., 1999Go, 2000Go). First, there is a set of genes that determine maximal aerobic capacity in the untrained condition (intrinsic). Second, another set of genes dictate the changes in aerobic capacity acquired as a result of exercise training (adaptational). Our aim was to devise a test of treadmill running performance for the rat that most closely estimated intrinsic aerobic capacity.

It has been demonstrated that oxygen consumption increases linearly as a function of running velocity over a wide range in rats (Gleeson and Baldwin, 1981Go). In accord with this, Brooks and White (1978)Go reported a significant association (R = 0.83) between oxygen consumption and running velocity (range = 14.3 to 43.1 ml/kg/min) in untrained rats tested at a slope of 15%. The test we devised consisted of running each rat on a motorized treadmill up a 15 degree slope with incremental increases in speed until exhausted (Fig. 2). Exhaustion was operationally defined as the third time a rat could no longer keep pace with the speed of the treadmill and remained on the shock grid for two seconds rather than run. The total distance run until exhaustion (meters) was designated as the standard for the estimate of each rat's intrinsic endurance exercise capacity. Each rat was tested over five consecutive days. The single best trial of five was used as the best indicator of capacity determined by intrinsic genetic composition and least indicative of the environmental component of the trait (Barbato et al., 1998Go). This idea of estimating the genetic component from the one best day of performance rather than either the mean or the median of the five trials, for example, has two origins. First, the environment can have an infinite negative influence upon capacity by reducing the distance run to zero. Factors such as subtle differences in housing or daily handling could cause a genetically superior rat to perform below its capacity on a given day. Second, the environment can have only a finite positive influence upon any test of capacity such that, given the best environment, intrinsic capacity is limited by a given genotype. An extension of this logic is that the estimate of breeding value (distance run) has more error for low capacity rats compared to high capacity rats (Britton and Koch, 2001Go).



View larger version (14K):
[in this window]
[in a new window]
 
FIG. 2. The relationship between speed and time (y-axis) to distance run (x-axis) for the treadmill ramped-speed protocol used to estimate aerobic energy capacity in rats. This test was repeated over five days and the best distance run was used as the best estimate of aerobic energy capacity for artificial selection. (Adapted from Koch and Britton, 2001Go)

 
For genetic model development, it is imperative to reduce environmental variance to a minimum; so that, the measure and selection of the trait is primarily for the genetic component. General environmental variance arises from non-localized, more permanent circumstances such as room temperature, light, humidity, diet, and time of day (Falconer and Mackay, 1996Go). Despite experimental control for factors such as housing conditions and manner of testing, other hard to control factors like intrauterine conditions, birthing order, social interactions, and eating patterns may also contribute to general environmental variation among populations. Additionally, special environmental variance for an individual can originate from localized, temporary circumstances such as daily handling (Falconer and Mackay, 1996Go). While artificial selection presumably concentrates genetic components over several generations, it also inadvertently selects for lack of sensitivity to subtle variations in the general environment.

Any performance-based test, whether performed in humans or rats, contains an inherent set of behavioral factors related to age, motivation, tractability, and willingness to perform (Koch et al., 1998Go). Our selection test began when the rats were young adults (10 weeks of age). Each rat was introduced to treadmill exercise over a period of 5 days. The first two days of introduction to treadmill running consisted of simply placing the rat on the treadmill belt that was moving at a velocity of 10 m·min–1 (15° slope) and picking the rat up and moving it forward if it started to slide off the back of the belt. This was done over a period between 1–3 minutes. During introduction days 3–5, the belt speed was gradually increased up to 15 m·min–1 and failure to run caused the rats to slide off of the moving belt and onto a 15 x 15 cm electric shock grid that delivered 1.2 mA of current at 3 Hz. The rats were left on the grid for about 1.5 sec and then moved forward onto the moving belt. With this process, most rats learned to run for five minutes and avoid the mild shock. This amount of exposure to treadmill exercise is likely below that required to produce a significant change in aerobic capacity (Baldwin et al., 1977Go). The ability to achieve this minimal level of treadmill exercise at least once constituted the threshold performance necessary for inclusion in evaluation for maximal running capacity the following week.

A selection experiment for aerobic energy capacity was first done as a test case over three generations in a small outbred population of Sprague Dawley rats (Harlan Sprague-Dawley, Indianapolis, IN) shown in Figure 3 (Koch et al., 1998Go). The founder population (24 females and 28 males) ran to exhaustion by 396 meters with a wide range of scores (149 to 695 meters). The average time to exhaustion for each test was equivalent to about 25 minutes. An event of this duration utilizes aerobic metabolism as the major contributor for energy transfer. A high intensity of breeding selection produced low and high lines for aerobic energy capacity that differed by 70%. We thought this result warranted initiation of a more large-scale selection for treadmill running capacity.



View larger version (13K):
[in this window]
[in a new window]
 
FIG. 3. The response to selection for aerobic capacity in a small test population of Sprague-Dawley rats. A 70% difference for average aerobic capacity was produced by three generations of high intensity breeding. This suggests a heritable component for the phenotype of aerobic running capacity based on a treadmill running test. (Adapted from Koch et al., 1998Go)

 

    USE OF N:NIH STOCK AS FOUNDER POPULATION
 TOP
 SYNOPSIS
 INTRODUCTION
 INTRODUCTION
 MODELS OF COMPLEX DISEASE
 EVOLUTION-BASED SELECTION:...
 SELECTION FOR AEROBIC CAPACITY
 USE OF N:NIH STOCK...
 DESIGN OF SELECTION EXPERIMENTS:...
 CURRENT DATA ON RATS...
 References
 
In 1979, Carl Hansen and Karen Spuhler foresaw the need for a genetically heterogeneous rat stock from which selected models could originate (Hansen and Spuhler, 1984Go). At that time, there were about 100 inbred rat strains available; many of which have origin back to an original Wistar Stock. The Wistar stock was made from a relatively small set of breeders as were the Sprague-Dawley and the Long-Evans rats. After phylogenetic consideration, the following eight inbred strains representing the broadest spectrum of laboratory rats available were chosen by Hansen and Spuhler for outcross breeding: Maudsely Reactor (MR/N), Wistar (WN/N), Wistar Kyoto (WKY/N), Marshall 520 (M520), Fischer (F344), A x C 9935 Irish (ACI), Brown Norway (BN/SsN), and Buffalo (BUF/N) (Fig. 4). Since its formation, the N:NIH outbred stock has been maintained as 60 outcrossed families and is available as a research resource from the National Institutes of Health (NIH) Animal Genetic Resource (NIHAGR) located at the NIH (Bethesda, MD). The availability of N:NIH heterogeneous stock of rats is of major benefit for use as a starting population for artificial selection. Because the stock was derived from eight inbred strains, deleterious alleles which may invoke inbreeding depression are already selected against which increases the opportunity for breeding fecundity during selection. At least in theory, one may ultimately be able to determine which causative allelic variants from the eight progenitor strains segregate in a selection process (Flint and Mott, 2001Go).



View larger version (26K):
[in this window]
[in a new window]
 
FIG. 4. Scheme for genetic outcross of 8 inbred rat strains as described by Hanson and Spuhler (1984)Go for the development of the genetically heterogeneous N:NIH stock of rats. N:NIH rats are a useful founder population for selection experiments

 
Large-scale selection for treadmill endurance running capacity was started using a founder population of 96 male and 96 female N:NIH rats (Koch and Britton, 2001Go). Each rat in the founder population was of different parentage, so that initial selection was not among brothers and sisters to obtain wider initial genetic variance (Hartl and Clark, 1988Go). On average (mean ± standard error), the founder population ran to exhaustion by 355 ± 11 m (23 min) with the females running significantly further than the males. The females in the N:NIH founder population ran on average for 380 ± 15 m and the males ran on average for 327 ± 16 m. Figure 5 (upper) is a frequency histogram for running capacity in the N:NIH founder population: females and males combined.



View larger version (14K):
[in this window]
[in a new window]
 
FIG. 5. Frequency distributions for aerobic running capacity in N:NIH rats. The upper figure shows a normal distribution and wide variation (range 35 to 810 meters) in a founder population (n = 168 rats). The lower figure shows results after 10 generations of selection on low and high aerobic running capacity. The magnitude of the response to selection was higher in rats selected for high capacity (+550 meters) in comparison to the low selected (–138 meters). (Figure adapted from Koch and Britton, 2001Go; Wisloff et al., 2005Go)

 

    DESIGN OF SELECTION EXPERIMENTS: ADVANTAGE OF USING DIVERGENT WITH-FAMILY ROTATIONAL BREEDING AND NON-REPLICATE LINES
 TOP
 SYNOPSIS
 INTRODUCTION
 INTRODUCTION
 MODELS OF COMPLEX DISEASE
 EVOLUTION-BASED SELECTION:...
 SELECTION FOR AEROBIC CAPACITY
 USE OF N:NIH STOCK...
 DESIGN OF SELECTION EXPERIMENTS:...
 CURRENT DATA ON RATS...
 References
 
The usual goal of selection is to induce a response attributed from increased frequency of alleles causative of the change in phenotype on which selection criteria is based (Falconer and Mackay, 1996Go). Other sources of response include: 1) sampling errors in estimating the generation means, 2) random changes in gene frequency (i.e., random genetic drift), and 3) unpredictable and subtle changes in the environment. Variation from sampling errors and genetic drift are reduced by first selecting extremes from a large number of animals in the founder population and then using a large number of families for selection at each generation. It is projected that the initial response to selection is greater when the population is smaller and the intensity of selection is high; yet this approach also increases the amount of inbreeding. For relatively short-term breeding programs, inbreeding of about 1% per generation is acceptable to yield retention of genetic variation and thus increase the overall response to selection (Falconer and Mackay, 1996Go). For randomly bred animals with an equal representation of males and females, the rate of inbreeding per generation (DF) is calculated as:

(where N = number of parents).

Making the contributions from each family nearly equal can further reduce the rate of inbreeding. This can be accomplished by selecting a male and female offspring from each family and using them as parents for the next generation. This is termed within-family selection. Kimura and Crow (Kimura and Crow, 1963Go) have evaluated the rates of inbreeding in several rotational breeding programs. Their results demonstrate that rotational within-family selective breeding gives rates of inbreeding quite close to the theoretical rate for minimally low inbreeding for equal representation of families (DF = [1/4N]) (Falconer, 1976Go). Thus, 13 rotational breeding pairs per selected line places the rate of inbreeding at slightly less than 1% per generation (1/[4*26] = 1/104 = 0.96%).

A prearranged schedule of 13 matings for within-family rotational breeding follows a simple pairing sequence based on assigned family number (1 to 13, F = female, M = male) as shown:


{i1540-7063-045-03-0405-eq2}

When the rotation has completed one entire cycle (i.e., 13 generations), the 1 x 1, 2 x 2, 3 x 3 etc. matings are skipped to avoid sib-matings (Falconer, 1976Go).

Within-family selection is beneficial for two additional reasons. First, it eliminates a large component of environmental variance such as maternal or pre-weaning effects that can be shared if selecting several members from one family. Second, if single-pair matings are to be made, every family contributes two members equally as genetic substrate of the next generation. The contribution from random genetic drift, however, can only be estimated by replicating selected lines. Although it is ideal to utilize replicated lines, the scale of any breeding program is ultimately resource limited. The maintenance of replicate lines operates to decrease fractionally the size of each line to 1/(n + 1) compared with no replicate lines (n = number of replicate lines) (Falconer and Mackay, 1996Go). Our decision was to reduce sampling error and contribution from genetic drift by employing all resources for development of only one low and one high line. The tradeoff was that we have no direct estimate of the contribution of genetic drift.


    CURRENT DATA ON RATS SELECTED FOR AEROBIC CAPACITY
 TOP
 SYNOPSIS
 INTRODUCTION
 INTRODUCTION
 MODELS OF COMPLEX DISEASE
 EVOLUTION-BASED SELECTION:...
 SELECTION FOR AEROBIC CAPACITY
 USE OF N:NIH STOCK...
 DESIGN OF SELECTION EXPERIMENTS:...
 CURRENT DATA ON RATS...
 References
 
The selection response over ten generations from our ongoing experiments is shown in Figure 6 (Koch and Britton, 2001Go). The 13 lowest founder females averaged 205 ± 7 m and the lowest 13 males averaged 167 ± 7 m for aerobic running capacity and were randomly mated. The 13 highest founder females averaged 633 ± 26 m and the 13 highest males averaged 541 ± 21 m for aerobic capacity and were also randomly mated to start the high line. This represents a 3.1-fold difference for the females (633/205) and a 3.2-fold difference for the males (541/167) between the low and high parents selected as progenitors for the selected lines. After ten generations of selection, a total of 2,589 rats were phenotyped for maximal aerobic running capacity. On average the low line members decreased 15 m per generation and the high line members increased 46 m per generation. The overall divergent response was 317% (688 m). At generation 10, the low-selected line averaged 217 ± 10 m (17 min) and the high line averaged 905 ± 22 m (45 min). It is obvious from inspection of the frequency distributions at generation 10 (Fig. 5 lower) that the high line responded more strongly to selection (+550 m) compared to the low line (–138 m). Selection for the low trait character is often of lesser magnitude relative to the high for two primary reasons (Falconer and Mackay, 1996Go). First, selection can be lower simply because the low trait approaches zero or a designated threshold as a limit. Second, it is possible that natural selection has constrained variation in aerobic capacity at the low end and competes against artificial selection to maintain a minimal threshold capacity associated with fitness. In contrast, more genetic variation may exist for high aerobic capacity because this trait has been under less evolutionary pressure (Mousseau and Roff, 1987Go). We named these lines Low Capacity Runners (LCR) and High Capacity Runners (HCR) which best describes the selected trait (Koch and Britton, 2001Go).



View larger version (14K):
[in this window]
[in a new window]
 
FIG. 6. The overall response to selection for low and high aerobic running capacity over 10 generations of selection was 317% difference between lines. On average the high line increased 46 m per generation and the low line decreased 15 m per generation. The mean and SE from a total of 2,589 rats are represented. (Figure adapted from Wisloff et al., 2005Go)

 
A change in an unselected trait produced by selection for another trait is termed a correlated response. Correlated traits are informative and of interest to selection studies for detecting positive and negative associations between traits but can not establish direct cause and effect for three primary reasons (Falconer and Mackay, 1996Go). First, the common property of pleiotropy, a genetic cause of trait differences determined by allelic variants showing a broad spectrum of effects on the expression of many traits, may be undetectable due to the net effect of all the segregating loci. Second, in connection with the process of selection, two traits may be independently influenced by a common selection environment. Third, correlated responses may be directed by evolutionary fitness and thus interconnected to viability and fertility of the lines.

Selection for aerobic running capacity produced a correlated change in body weight in both females and males. In general, the low line (LCR) became heavier by approximately 3 grams per generation and the high line (HCR) became increasingly lighter by an estimate of 1 gram per generation. By generation 10, the LCR females weighed 213 ± 3 g and the HCR females weighed 172 ± 2 g (24% difference). Similarly, the LCR males weighed 310 ± 6 g and the HCR males weighed 255 ± 3 g (22% difference). Multiple linear regression analysis showed that the concomitant changes in body weight across generations accounted for about 7% of the variation in distance run in females and 14–20% in males (Wisloff et al., 2005Go) As such, LCR and HCR lines might serve as contrasting models to determine selected factors that influence body weight and composition.

It seems reasonable to predict that two-way selection for treadmill running capacity would produce lines that differ in the capacity for oxygen utilization as governed by the cardiac output and the ability of peripheral tissues to extract oxygen will be at least one of the predominant factors selected, although, other factors such as skeletal mechanics, neuromuscular coordination, and heat dissipation may also be involved. At generation 7, the VO2max (maximal oxygen consumption) was 21% higher in HCR compared to the LCR rats (Henderson et al., 2002Go). For evaluation, oxygen transport system was considered as four sequential components of oxygen utilization: 1) pulmonary ventilation, 2) alveolar-capillary diffusion, 3) convective transport from blood to tissues, and 4) blood–to-tissue transfer. A 30% difference in oxygen transfer at the tissue level was the major factor contributing to the difference in VO2 max between the LCR and HCR rats (Henderson et al., 2002Go). Although the total capillary and fiber number were similar, the fiber area in the HCR was 37% less (Howlett et al., 2003Go). As a result, the number of capillaries per unit area of muscle was 32% more in the HCR compared to the LCR rats. Both citrate sythase and ß–hydroxyacyl-CoA dehydrognease enzyme activities were about 40% higher in the HCR relative to LCR while phosphofructokinase was lower (Table 1).


View this table:
[in this window]
[in a new window]
 
TABLE 1. Skeletal muscle capillarity and muscle enzyme activities from the medial gastrocnemius muscle of rat selected lines for low aerobic running capacity (LCR) and high aerobic running capacity (HCR). Data from Howlett et al.,2003

 
High blood pressure, depressed cardiac contractility, body weight, and elevated blood lipid levels are standard estimates of cardiovascular disease risk. Abdominal aortic blood pressure measured via telemetry was higher in the LCR during the day (105 ± 13 vs. 89 ± 8 mm Hg), at night (98 ± 3 vs. 91 ± 7 mm Hg), and for the combined 24-hour period (102 ± 6 vs. 90 ± 7 mm Hg) (Fig. 7A). Stroke volume was 47% less in a working Langendorff preparation and field-stimulated isolated left ventricular myocytes shortened 22% less in LCR relative the HCR (Fig. 7B, 7C). Visceral adiposity relative to body weight was 63% higher, plasma free fatty acids 95% higher, and plasma triglycerides 167% higher in LCR compared to HCR rats (Fig. 2D, E, F). These data suggest increased cardiovascular risk factors associate with low aerobic capacity (Hussain et al., 2001Go; Wisloff et al., 2005Go).



View larger version (17K):
[in this window]
[in a new window]
 
FIG. 7. Rats selected on the basis of low versus high aerobic running capacity differ in several indices of cardiovascular risk factors linked to complex disease. A. Mean blood pressures recorded via telemetry in LCR (n = 10) over a 24-hour period was higher than in the HCR (n = 10). B. Cardiac contractility as assessed by isolated working heart preparation was lower for stroke volume (µl per gram heart weight) in LCR hearts. C. Cardiomyocytes isolated from LCR rats demonstrated lesser percent shortening when field stimulated. D. Visceral adiposity (% per body weight) was greater in LCR, E. Free fatty acid levels (mEQ/L), and F. Plasma triglycerides (mg/dL) were elevated in LCR by comparison to HCR. Values are Mean ± SE. (Figure adapted from Wisloff et al., 2005Go)

 
In view of the lower aerobic capacity we hypothesized that LCR have compromised mitochondrial oxidative function relative to the HCR rats. To test this, we measured the cellular content of proteins required for mitochondrial biogenesis and function (Mootha et al., 2003Go) in soleus muscle, which is composed largely of highly oxidative type I fibers. As Figure 8 reveals, the content peroxisome proliferative activated receptor, gamma, coactivator 1, alpha (PGC-1{alpha}), peroxisome proliferative activated receptor, gamma (PPAR-{gamma}), ubiquinol-cytochrome c oxidoreductase core 2 subunit (UQCRC2), cytochrome c oxidase subunit I (COXI), uncoupling protein 2 (UCP2), and ATP synthase H+ transporting mitochondrial F1 complex (F1-ATP synthase) was markedly reduced in the LCR rats by comparison to the HCR (Wisloff et al., 2005Go). The uniform decline in these proteins is consistent with the hypothesis that reduced aerobic metabolism plays a causal role in the development of the differences between the LCR and HCR rats. Although our experiments do not prove a direct cause and effect relationship, they are consonant with the view that impaired regulation of oxidative pathways in mitochondria may be a common factor linking reduced total-body aerobic capacity to complex disease.



View larger version (26K):
[in this window]
[in a new window]
 
FIG. 8. Proteins known to be integral for mitochondrial function were more highly expressed in soleus skeletal muscles from HCR relative to LCR. A. Representative Western blots from LCR and HCR rats for six proteins: 1) peroxisome proliferative activated receptor {gamma} (PPAR-), 2) PPAR-{gamma} coactivator 1 {alpha} (PGC-1{alpha}), 3) ubiquinol-cytochrome c oxidoreductase core 2 subunit (UQCRC2), 4) cytochrome c oxidase subunit I (COXI), 5) uncoupling protein 2 (UCP2), and 6) ATP synthase H+-transporting mitochondrial F1 complex (F1-ATP synthase). B. Mean values (± 1 standard deviation) for expression of each of the six proteins (n = 6 LCR and 6 HCR muscles). Y axis units are densitometric measures for each protein relative to the {alpha}-actin signal. (Figure adapted from Wisloff et al., 2005Go)

 


    ACKNOWLEDGMENTS
 
We appreciate the collaborative efforts with U. Wisloff, O. Ellingsen, S. Najjar, S. Swoap, P. D. Wagner, R. Howlett, and N.C. Gonzalez, This work was supported by grants from the United States Public Health Service, National Institutes of Health, (Heart, Lung and Blood Institute HL 64270 and National Center for Research Resources RR 17718) to LGK and SLB.


    FOOTNOTES
 
1 From the Symposium Selection Experiments as a Tool in Evolutionary and Comparative Physiology: Insights into Complex Traits, at the Annual Meeting of the Society for Integrative and Comparative Biology, 5–8 January 2004, at New Orleans, Louisiana, USA. Back

2 E-mail: lgkoch{at}med.umich.edu Back


    References
 TOP
 SYNOPSIS
 INTRODUCTION
 INTRODUCTION
 MODELS OF COMPLEX DISEASE
 EVOLUTION-BASED SELECTION:...
 SELECTION FOR AEROBIC CAPACITY
 USE OF N:NIH STOCK...
 DESIGN OF SELECTION EXPERIMENTS:...
 CURRENT DATA ON RATS...
 References
 
Ahn, J., J. Varagic, M. Slama, D. Susic, and E. D. Frohlich. 2004. Cardiac structural and functional responses to salt loading in SHR. Am. J. Physiol, 287:H767-H772.

Baldwin, K. M., D. A. Cooke, and W. G. Cheadle. 1977. Time course adaptations in cardiac and skeletal muscle in different running program. J. Appl. Physiol, 42:267-272.[Abstract/Free Full Text]

Barbato, J. C., L. G. Koch, A. Darvish, G. T. Cicila, P. J. Metting, and S. L. Britton. 1998. Spectrum of aerobic endurance running performance in eleven inbred strains of rats. J. Appl. Physiol, 85:530-536.[Abstract/Free Full Text]

Bouchard, C., P. An, T. Rice, J. S. Skinner, J. G. Wilmore, J. Gagnon, L. Perusse, A. S. Leon, and D. C. Rao. 1999. Familial aggregation of VO2 max response to exercise training: Results from the Heritage family study. J. Appl. Physiol, 87:1003-1008.[Abstract/Free Full Text]

Bouchard, C., T. Rankinen, Y. C. Chagono, T. Rice, L. Perusse, J. Gagnon, I. Borecki, P. An, A. S. Leon, J. S. Skinner, J. H. Wilmore, M. Province, and D. C. Rao. 2000. Genomic scan for maximal oxygen uptake and its response to training in the HERITAGE family study. J. Appl. Physiol, 88:551-559.[Abstract/Free Full Text]

Brand, P. H., K. B. Coyne, K. A. Kostrzewski, D. Shier, P. J. Metting, and S. L. Britton. 1991. Pressure diuresis and autonomic function in conscious dogs. Am. J. Physiol, 261:R802-R810.

Britton, S. L., and L. G. Koch. 2001. Animal genetic models for complex traits of physical capacity. Exerc. Sport. Sci. Rev, 29:7-14.[CrossRef][Medline]

Brooks, G. A., and T. P. White. 1978. Determination of metabolic and heart rate responses of rats to treadmill exercise. J. Appl. Physiol, 45:1009-1015.[Free Full Text]

DesMarias, D. J. 2000. When did photosynthesis emerge on earth? Science, 289:1703-1705.[Free Full Text]

Duchen, M. R. 2004. Mitochondria in health and disease: Perspectives on a new mitochondrial biology. Mol. Aspects Med, 25:365-451.[Medline]

Falconer, D. S. 1976. Genetic aspects of breeding methods. In UFAW Handbook on the care and management of laboratory animals, pp. 7–26. Churchill Livingstone, Edinburgh.

Falconer, D. S., and T. F. C. Mackay. 1996. Introduction to quantitative genetics. 4th ed. Addison Wesley Longman, Ltd., Essex, England.

Flint, J., and R. Mott. 2001. Finding the molecular basis of quantitative traits: Successes and pitfalls. Nat. Rev. Genet, 2:437-445.[CrossRef][Web of Science][Medline]

Gensini, G. F., M. Comeglio, and A. Calella. 1998. Classical risk factors and emerging elements in the risk profile for coronary artery disease. Eur. Heart J, 19:A53-A61.

Gleeson, T. T., and K. M. Baldwin. 1981. Cardiovascular response to treadmill exercise in untrained rats. J. Appl. Physiol, 50:1206-1211.[Abstract/Free Full Text]

Goldblatt, H., J. Lynch, R. J. Hanzal, and W. W. Summerville. 1934. The production of persistent elevation of systolic blood pressure by means of renal ischemia. J. Exp. Med, 59:347-380.[Abstract]

Greenhouse, D. D., M. F. W. Festing, S. Hasan, and A. L. Cohen. 1990. Inbred strains of rats. In H. J. Hedrich (ed.), Genetic monitoring of inbred strains of rats, pp. 410–480. Gustav Fisher, Stuttgart.

Guyton, A. C. 1991. Blood pressure control—special role of the kidneys and body fluids. Science, 252:1813-1816.[Abstract/Free Full Text]

Guyton, A. C., and T. G. Coleman. 1967. Long-term regulation of the circulation: Interrelationships with body fluid volumes. In E. B. and A. C. Guyton (eds.), W. B. Physical bases of circulatory transport, regulation, and exchange, pp. 179–201. Saunders Co., Philadelphia, Pennsylvania.

Hansen, C., and K. Spuhler. 1984. Development of the National Institutes of Health genetically heterogeneous rat stock. Alcohol. Clin. Exp. Res, 8:477-479.[Web of Science][Medline]

Hartl, D. L., and A. G. Clark. 1988. Principles of population genetics. Sinauer Associates, Inc., Sunderland, Massachusettes.

Henderson, K. K., H. Wagner, F. Favret, S. L. Britton, L. G. Koch, P. D. Wagner, and N. C. Gonzalez. 2002. Determinants of maximal O2 uptake in rats selectively bred for endurance running capacity. J. Appl. Physiol, 93:1265-1274.[Abstract/Free Full Text]

Howlett, R. A., N. C. Gonzalez, H. E. Wagner, Z. Fu, S. L. Britton, L. G. Koch, and P. D. Wagner. 2003. Selected contribution: Skeletal muscle capillarity and enzyme activity in rats selectively bred for running endurance. J. Appl. Physiol, 94:1682-1688.[Abstract/Free Full Text]

Hussain, S. O., J. C. Barbato, L. G. Koch, P. J. Metting, and S. L. Britton. 2001. Cardiac function in rats selectively bred for low (LCR) and high (HCR) capacity running. Am. J. Physiol, 281:R1787-R1791.

Jacob, H. J., and A. E. Kwitek. 2002. Rat genetics: Attaching physiology and pharmacology to the genome. Nat. Rev. Genet, 3:33-42.[CrossRef][Web of Science][Medline]

Kimura, M., and J. F. Crow. 1963. On the maximum avoidance of inbreeding. Genet. Res, 4:399-415.[Medline]

Koch, L. G., and S. L. Britton. 2001. Artificial selection for intrinsic aerobic endurance running capacity in rats. Physiol. Genomics, 5:45-52.[Abstract/Free Full Text]

Koch, L. G., T. A. Meredith, T. D. Fraker, P. J. Metting, and S. L. Britton. 1998. Heritability of treadmill running endurance in rats. Am. J. Physiol, 275:R1455-R1460.

Mootha, V. K., C. M. Lindgren, K. F. Eriksson, A. Subramanian, S. Sihag, J. Lehar, P. Puigserver, E. Carlsson, M. Ridderstrale, E. Laurila, N. Houstis, M. J. Daly, N. Patterson, J. P. Mesirov, T. R. Golub, P. Tamayo, B. Spiegelman, E. S. Lander, J. N. Hirschhorn, D. Altshuler, and L. C. Groop. 2003. PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat. Genet, 34:267-273.[CrossRef][Web of Science][Medline]

Mousseau, T. A., and D. A. Roff. 1987. Natural selection and the heritability of fitness components. Heredity, 59:181-197.

Myers, J., M. Prakash, V. Froelicher, D. Do, S. Partington, and J. E. Atwood. 2002. Exercise capacity and mortality among men referred for exercise testing. N. Engl. J. Med, 346:793-801.[Abstract/Free Full Text]

Paytan, A. 2000. Sulfate clues for the early history of atmospheric oxygen. Science, 288:626-627.[Free Full Text]

Pinet, F., F. Poirier, S. Fuchs, P. Tharaux, M. Caron, P. Carvol, J. Michel, and R. Joubert-Caron. 2004. Troponin T as a marker of diffenentiation revealed by proteomic analysis in renal arterioles. FASEB J, 18:585-586.[Abstract/Free Full Text]

Pross, A. 2003. The driving force for life's emergence: Kinetic and thermodynamic considerations. J. Theor. Biol, 220:393-406.[CrossRef][Web of Science][Medline]

Rapp, J. P. 2000. Genetic analysis of inherited hypertension in the rat. Physiol. Rev, 80:135-172.[Abstract/Free Full Text]

Rat genome sequencing project consortium., 2004. Genome sequence of the Brown Norway rat yields insights into mammalian evolution. Nature, 428:493-521.[CrossRef][Medline]

Rossini, A., A. A. Like, W. L. Chick, M. C. Appel, and G. F. Cahill Jr. 1977. Studies of streptozotocin-induced insulitis and diabetes. Proc. Natl. Acad. Sci. U.S.A, 74:2485-9.[Abstract/Free Full Text]

Skarlatos, S., P. H. Brand, P. J. Metting, and S. L. Britton. 1994. Spontaneous changes in arterial blood pressure and renal interstitial hydrostatic pressure in conscious rats. J. Physiol. 481, 3:743-752.

Steele, J. E., L. G. Koch, and P. H. Brand. 2000. State-dependent expression of pressure diuresis in conscious rats. Proc. Soc. Exp. Biol. Med, 224:109-15.[Abstract/Free Full Text]

Williams, R. S., and P. D. Wagner. 2000. Transgenic animals in integrative biology: Approaches and interpretations of outcome. J. Appl. Physiol, 88:1119-26.[Abstract/Free Full Text]

Wisloff, U., S. M. Najjar, O. Ellingsen, P. M. Haram, S. Swoap, Q. Al-Share, M. Fernstrom, K. Rezaei, S. J. Lee, L. G. Koch, and S. L. Britton. 2005. Cardiovascular risk factors from artificial selection for low aerobic capacity in rats. Science, 307:418-420.[Abstract/Free Full Text]

Xiong, J., W. M. Fischer, K. Inouw, M. Nakahara, and C. E. Bauer. 2000. Molecular evidence for the early evolution of photosynthesis. Science, 289:1724-1730.[Abstract/Free Full Text]

Young, I. S., and J. V. Woodside. 2001. Antioxidants in health and disease. J. Clin. Pathol, 54:176-186.[Abstract/Free Full Text]

Zolotareva, A. G., and M. E. Kogan. 1978. Production of experimental occlusive myocardial infarction in mice. Cor. Vasa, 20:308-14.[Web of Science][Medline]


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?


This article has been cited by other articles:


Home page
Integr. Comp. Biol.Home page
J. G. Swallow and T. Garland Jr.
Selection Experiments as a Tool in Evolutionary and Comparative Physiology: Insights into Complex Traits--an Introduction to the Symposium
Integr. Comp. Biol., June 1, 2005; 45(3): 387 - 390.
[Full Text] [PDF]


Home page
Integr. Comp. Biol.Home page
J. G. Swallow, J. S. Rhodes, and T. Garland Jr.
Phenotypic and Evolutionary Plasticity of Organ Masses in Response to Voluntary Exercise in House Mice
Integr. Comp. Biol., June 1, 2005; 45(3): 426 - 437.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (4)
Right arrow Request Permissions
Google Scholar
Right arrow Articles by Koch, L. G.
Right arrow Articles by Britton, S. L.
Right arrow Search for Related Content
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?