Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Genetic Analysis of a Rat Model of Aerobic Capacity and Metabolic Fitness

  • Yu-yu Ren,

    Affiliation Department of Human Genetics, University of Michigan, Ann Arbor, Michigan, United States of America

  • Katherine A. Overmyer,

    Affiliation Department of Internal Medicine, University of Michigan, Ann Arbor, Michigan, United States of America

  • Nathan R. Qi,

    Affiliation Department of Internal Medicine, University of Michigan, Ann Arbor, Michigan, United States of America

  • Mary K. Treutelaar,

    Affiliation Department of Internal Medicine, University of Michigan, Ann Arbor, Michigan, United States of America

  • Lori Heckenkamp,

    Affiliation Department of Anesthesiology, University of Michigan, Ann Arbor, Michigan, United States of America

  • Molly Kalahar,

    Affiliation Department of Anesthesiology, University of Michigan, Ann Arbor, Michigan, United States of America

  • Lauren G. Koch,

    Affiliation Department of Anesthesiology, University of Michigan, Ann Arbor, Michigan, United States of America

  • Steven L. Britton,

    Affiliation Department of Anesthesiology, University of Michigan, Ann Arbor, Michigan, United States of America

  • Charles F. Burant,

    Affiliation Department of Internal Medicine, University of Michigan, Ann Arbor, Michigan, United States of America

  • Jun Z. Li

    junzli@med.umich.edu

    Affiliation Department of Human Genetics, University of Michigan, Ann Arbor, Michigan, United States of America

Correction

24 Mar 2014: The PLOS ONE Staff (2014) Correction: Genetic Analysis of a Rat Model of Aerobic Capacity and Metabolic Fitness. PLOS ONE 9(3): e92345. https://doi.org/10.1371/journal.pone.0092345 View correction

Abstract

Aerobic capacity is a strong predictor of all-cause mortality and can influence many complex traits. To explore the biological basis underlying this connection, we developed via artificial selection two rat lines that diverge for intrinsic (i.e. inborn) aerobic capacity and differ in risk for complex disease traits. Here we conduct the first in-depth pedigree and molecular genetic analysis of these lines, the high capacity runners (HCR) and low capacity runners (LCR). Our results show that both HCR and LCR lines maintain considerable narrow-sense heritability (h2) for the running capacity phenotype over 28 generations (h2 = 0.47 ± 0.02 and 0.43 ± 0.02, respectively). To minimize inbreeding, the lines were maintained by rotational mating. Pedigree records predict that the inbreeding coefficient increases at a rate of <1% per generation, ~37-38% slower than expected for random mating. Genome-wide 10K SNP genotype data for generations 5, 14, and 26 demonstrate substantial genomic evolution: between-line differentiation increased progressively, while within-line diversity deceased. Genome-wide average heterozygosity decreased at a rate of <1% per generation, consistent with pedigree-based predictions and confirming the effectiveness of rotational breeding. Linkage disequilibrium index r2 decreases to 0.3 at ~3 Mb, suggesting that the resolution for mapping quantitative trait loci (QTL) can be as high as 2-3 cM. To establish a test population for QTL mapping, we conducted an HCR-LCR intercross. Running capacity of the F1 population (n=176) was intermediate of the HCR and LCR parentals (28 pairs); and the F2 population (n=645) showed a wider range of phenotypic distribution. Importantly, heritability in the F0-F2 pedigree remained high (h2~0.6). These results suggest that the HCR-LCR lines can serve as a valuable system for studying genomic evolution, and a powerful resource for mapping QTL for a host of characters relevant to human health.

Introduction

Increased intrinsic exercise capacity or aerobic capacity, generally measured as maximal work with a standardized treadmill test, is an excellent predictor of disease risk in humans, with higher capacities associated with enhanced health and resistance to metabolic disease [19]. Peak exercise capacity is a better predictor of mortality than other established risk factors such as hypertension, smoking, and diabetes [7,10]. This frequently observed statistical association suggests a causal connection. However, the biological basis for this connection remains largely unknown [11,12]. Aerobic capacity is a complex phenotype. While some studies document a strong genetic component of maximal oxygen uptake [13,14], other studies report low estimates of heritability [15,16]. In studies involving human subjects, it is often difficult to resolve the relative contributions of genetic and environmental factors to the inter-individual variability of aerobic capacity [13,1721], e.g., to resolve the relative effects of innate endurance from those due to aerobic training.

To understand the genetic and functional basis of aerobic capacity we sought to establish an animal model that allows in-depth analyses of the biology and health impact of innate aerobic capacity. In 1996, a long-term experiment was initiated to create two lines of rats through bi-directional selection for untrained aerobic running capacity [22]. The two lines, termed high capacity runners (HCR) and low capacity runners (LCR), originated from a founder population of 186 genetically heterogeneous rats derived from outcrossing 8 inbred strains (N:NIH stock) [23]. The animals were selected by their performance in run-to-exhaustion tests on a progressively accelerating treadmill, with the highest and lowest runners, one for each sex from each of 13 families, entering into a rotational breeding scheme (see Methods).

One of the original aims for establishing the HCR-LCR lines was to test the hypothesis that artificial selection based on intrinsic aerobic capacity would yield models that also exhibit contrasts in disease risks. This hypothesis has been proven correct: after 28 generations of selection, the HCR and LCR diverged not only for running capacity, but also in other physiological measures, including blood pressure, body mass index, lung capacity, lipid and glucose metabolism [24]. The LCR, relative to the HCR, manifest numerous clinically relevant conditions, including increased susceptibility to cardiac ventricular fibrillation [25] and hepatic steatosis [26]. At the behavioral level the LCR score higher for dysfunctional sleep [27], diminished behavioral strategies for coping with stress [28], and impaired memory and learning [29]. In contrast, the HCR have reduced weight gain [30], increased resistance to the deleterious effects of a high fat diet [31,32], increased capacity for fatty acid oxidation in skeletal muscle [33] and liver [26], and a 28-45% increase of lifespan [34].

A major advantage of the HCR-LCR system is that the pedigree and running phenotype data (n = 11,422) are completely known; and tissue samples for most breeding members (n > 1,500) have been archived. This combination of existing data and reagents, combined with over 70 published physiological studies of the two lines, represents a valuable resource that allows comprehensive analyses of the effects of selection on genomic and phenotypic evolution.

In this study, we carried out a systematic analysis of the running phenotype and related traits over the known pedigree of 0-28 generations. We also collected a genome-wide 10K SNP dataset for a subset of breeding members from three generations (G) (n=142 over G5, G14, and G26), and used these data to examine patterns of genomic evolution in the two lines as they undergo selection. Finally, we performed the first intercross experiments between the HCR and LCR, and analyzed the phenotypic distribution and heritability of the F1 (n = 176) and F2 populations (n = 645). These analyses provided new insights into the genealogical structure, inbreeding patterns, and genetic variability of the two lines, and characterized the intercross animals to assess their suitability as a mapping population for identifying quantitative trait loci (QTL).

Results

Rotational breeding and inbreeding coefficients

The protocols of animal maintenance, phenotyping, and rotational breeding have been described previously [22] (see also Methods). We analyzed the pedigree data for generations 1-28 (Files S1-S2), involving 5,976 HCRs and 5,446 LCRs. For each animal, we calculated its expected inbreeding coefficient (F) by tracing its parental lineage and documenting inbreeding loops. As expected, such pedigree-based estimates of F started to rise at G4-G5 and continued to increase over successive generations (Figure 1). The breeding history included occasional out-of-schedule pairings due to the lack of offspring of a certain sex in a given family or the need to substitute for unproductive mating pairs (see Methods for more details). Despite this, the pattern of F increase in actual pedigrees largely agrees with the expectation assuming perfect adherence to the planned rotation schedule (shown in solid lines in Figure 1). The cyclic rise of inbreeding coefficient every six generations is expected for 13 breeding pairs, due to the inevitable first-cousin paring every half cycle of the rotation [35,36] (further explained in Methods). Importantly, the average increase of estimated F is 0.94% and 0.95% per generation for HCR and LCR, respectively. These predictions are slower than the rate expected under random mating (shown in dotted lines in Figure 1), which increase at 1.51 per generation for both HCR and LCR, starting from the first generation.

thumbnail
Figure 1. Distribution of predicted inbreeding coefficients (F) for generations 0 to 28.

Shown are "violin-plots" for individual generations for HCR (A) and LCR (B). The widths of the ovals indicate the probability density of the data values. The black dots and the thick black lines in the ovals denote the median and the 25-75 percentile range, respectively. The dotted blue line indicate the expected increase in F under random mating given the 13-family breeding scheme, and the solid blue lines indicate the expected F under perfect adherence to the rotational breeding scheme.

Phenotypic response to selection and heritability

For each animal, we collected phenotype data that include maximal running distance, body weight at the time of running trial, and vertical work performed during each run. All animals were tested at 11-12 weeks of age. While both lines were derived from the same base population (indicated in yellow in Figure 2), their running performance gradually diverged over time. After 28 generations, the HCRs and LCRs differ by about 8.3-fold in running distance (9 times of the average within-line standard deviation), compared to ~2.8 fold (range of 298 to 840 meters) among eleven inbred lines commonly used in research [37]. The HCR continue to respond to selection (Figure 2, Table 1-2), with maximal running distance reaching >2000 m, ~2.4 fold higher than the best recorded performance among the inbred lines [37]. The pattern of increase is consistent in both males and females (Figure S1). Body weight increased in LCR and decreased in HCR In the first 12-13 generations, but did not diverge further after G13, stabilizing to a 0.7 to 0.8-fold difference through 28 generations (Figure 3, Table 1-2). In general, females are of lighter weight than males. However, as females tend to run longer, the overall vertical work is near-equivalent between males and females (1.3-fold difference in HCR, 1.1-fold difference in LCR) and larger in HCR than LCR by 6.8 fold.

thumbnail
Figure 2. Distribution of maximal running distance for generations 0 to 28.

Shown are "violin-plots" for individual generations for females and males combined. The yellow oval to the left denotes the founder population (NIH:H, n=153 phenotyped, out of 186), while green and red ovals are for HCR and LCR, respectively. The blue tick marks on the y axis indicate the maximal running distance for eleven inbred lines, which are ordered, from top to bottom, as DA (840m), PVG (718m), AUG (699m), SR (533m), F344 (469m), ACI (450m), LEW (442m), WKY (414m), BUF (373m), MNS (308m) and COP (298m).

MaleFemale
Best Running Distance (meters)Body Weight (grams)Best Running Distance (meters)Body Weight (grams)
GennMeanSDMeanSDnMeanSDMeanSD
157529.0244.3265.727.559564.7224.3181.220.7
282495.3201.9271.329.055652.2263.0174.216.7
374630.4215.1258.229.386721.9248.5175.018.3
464611.3190.2261.425.467784.3270.4176.015.9
586661.6218.0252.926.265812.9270.0171.914.2
679769.9220.0250.628.272916.5271.1169.713.4
767722.6305.7274.534.5661116.2420.6174.417.9
896543.7313.0263.935.489795.8420.8174.115.7
9118600.0272.8257.429.3116889.2381.9172.317.0
10140799.0261.6257.732.6981018.7308.0172.317.5
1192641.2281.3246.230.5110849.3314.0164.517.6
1290756.1348.6251.327.1105883.5528.5166.416.0
13134729.4386.5236.828.61061024.2651.6158.819.1
141121243.2443.2220.932.21101703.6602.7154.916.2
151221667.0284.5221.028.51102023.8424.5152.115.1
161061476.0253.0224.125.11041475.0236.9158.017.2
171461509.2248.9227.428.11371595.6264.7157.417.8
181491540.9241.1220.829.41211658.5256.4156.517.5
191261669.0271.8225.726.21011787.4263.0160.115.4
201351579.2246.0233.726.01131661.3233.2156.615.0
211151834.0306.7233.127.11201934.1291.0160.214.6
221401729.6250.1254.228.21181902.0279.5167.816.0
231221697.3279.0254.031.01051867.1290.6173.516.8
241101649.4226.5258.025.71051834.2264.9173.818.5
251241657.0279.2257.725.21401913.6258.5172.515.2
261681834.9277.7251.823.11292053.1357.0163.512.1
271452000.7364.3246.326.31452295.2392.6163.814.8
281121810.4269.6261.426.11132108.6297.5169.114.0

Table 1. Summary of cohort size, running distance, and body weight by gender and by generation for HCRs.

Total number of animals generated, including animals produced for line maintenance and resource sharing.
CSV
Download CSV
MaleFemale
Best Running Distance (meters)Body Weight (grams)Best Running Distance (meters)Body Weight (grams)
GennMeanSDMeanSDnMeanSDMeanSD
161372.3157.1280.629.064421.1176.7178.317.2
258268.7116.4296.731.268333.8123.6187.512.9
386309.5119.5287.334.281325.396.2186.317.2
472302.689.2297.627.866366.9106.0194.716.0
558252.473.0299.938.492326.1104.9194.719.1
674281.773.0293.831.181333.9108.8195.719.2
778263.0139.2312.729.563356.3150.3203.817.6
870157.793.4313.624.273259.6104.8209.423.9
995171.765.3320.839.3103260.381.8205.520.6
1076167.394.4321.925.674245.088.5213.818.8
11124156.560.9318.638.7121217.772.4207.220.8
1298182.890.6321.436.5104267.9102.8211.921.5
13109150.451.6336.231.5107215.269.5216.916.7
1498205.857.8304.232.6111266.478.4207.220.9
15111240.375.9318.340.4120318.7102.7214.421.4
1663303.269.3314.329.978370.0107.3208.028.0
17127318.375.6316.835.5127367.487.2203.819.4
18127282.658.2323.532.8119374.669.3202.718.0
19123281.378.4327.732.1115363.992.0213.619.1
20113272.364.0326.729.2120330.787.3208.019.4
21112291.390.4336.330.7114356.383.8203.919.2
2299328.963.6335.130.8105389.265.8202.316.6
23103248.338.0339.832.388322.545.7210.817.6
2488238.130.1354.924.384308.947.0214.819.2
25115228.065.2340.633.2110346.573.3205.317.7
26117201.745.1327.050.3123305.349.4206.219.4
27135217.849.6331.631.7119309.364.1206.414.4
28115191.842.9328.231.5111278.345.8208.616.7

Table 2. Summary of cohort size, running distance, and body weight by gender and by generation for LCRs.

Total number of animals generated, including animals produced for line maintenance and resource sharing.
CSV
Download CSV
thumbnail
Figure 3. Distribution of body weight for generations 0 to 28 for HCRs

(A) and LCRs (B). The yellow oval denotes the founder population (NIH:H, n=153). Most distributions are bi-modal, as most females are of lower body weight than most males. The blue tick marks on the y axis indicate the body weight for eleven inbred lines. They are ordered, from top to bottom, as MNS (279g), LEW (265g), SR (254g), BUF (246g), WKY (210g), COP (199g), ACI (193g), F344 (188g), DA (180g), PVG (179g) and AUG (167g).

The narrow-sense heritability (h2) for the logarithm of running distance, which measures the proportion of total phenotypic variance explained by additive effects of genes, was calculated for each line separately, and was 0.47 ± 0.02 in HCRs and 0.43 ± 0.03 in LCRs when all 28 generations were considered. To evaluate potential change in heritability over time, we also calculated h2 over four-generation, partially overlapping, intervals and found that while h2 was variable across intervals, it maintained positive values, with no sign of abatement in later generations (Figure S2). The within-line h2 for bodyweight and vertical work are 0.45 ± 0.06 and 0.37 ± 0.02, respectively, for HCRs, and 0.17 ± 0.03 and 0.58 ± 0.02 for LCRs. These results indicate that although LCRs did not show a decrease of running capacity as dramatically as the increase in HCRs, the heritability of running performance was comparable in the two lines. Lower body weight is associated with better running capacity, as shown by the negative correlations between the two phenotypes for both sexes within each line. For HCRs, the spearman correlation is -0.41 ± 0.15 and -0.17 ± 0.16 for males and females, respectively. For LCRs, the correlation is -0.19 ± 0.16 and -0.12 ± 0.13 for males and females, respectively. The fact that HCRs continued to respond to selection and that both lines maintained within-line heritability suggest that causal genetic variants have not been fixed in either line, rather they continue to segregate in both pedigrees.

Increased genomic differentiation between lines

As an initial genetic characterization of the HCR-LCR system, we collected genotype data over a genome-wide panel of ~10K single nucleotide polymorphism (SNP) loci for 142 animals, consisting of 22-25 animals in each of three non-adjacent generations (G5, G14, and G26) in both lines (see Methods). We chose these three generations to profile the long-term genomic changes in the two lines.. We used the average heterozygosity of 61 X chromosome (ChrX) markers to infer sex (Figure S3A), and found no disagreement with the recorded sex of the 142 animals. Documented relatedness was also confirmed by plotting the pairwise proportion of not sharing DNA segments due to identity by descent (P(IBD) = 0, or Z0) against the proportion of sharing one copy IBD (Z1) using the 2,518 SNPs (SNP Panel-2, see Methods). Siblings and non-sibling relatives are separated into distinct clusters (Figure S3B), indicating that the sample identities reflected in the genotype data are consistent with the recorded pedigree. Multidimensional Scaling (MDS) analysis showed that at G5, HCRs and LCRs formed two readily separable clusters (Figure 4). From G5 to G14 and from G14 to G26, between-lines separation increased, indicating a progressively greater differentiation between LCR and HCR.

thumbnail
Figure 4. Progressive genetic differentiation revealed by 10K SNP genotyping data.

Shown is a multidimensional scaling plot (dimensions 1 vs. 2) for 142 genotyped animals in two lines and three generations, as indicated by different symbols, showing that the two lines formed separate clusters at G5, and diverged further in G14 and G26.

To measure the apportionment of total genetic variance into between-line and within-line components we performed an Analysis of Molecular Variance (AMOVA) [38]. The proportion of variance explained by among population difference, as a weighted average over all loci, is increasing over time, from 6.5% at G5, to 15.6% at G14, and to 26.5% at G26.

Decreased genomic diversity within lines

We analyzed genetic diversity at the individual level by calculating the average heterozygosity (H0) across the 2,518 Panel-2 markers for each animal, and averaging within each of the six groups (two lines, at three time points) (Figure 5). At G5, H0 averaged 0.379 in HCR and 0.372 in LCR. At G14 it decreased to 0.338 (-10.8%) in HCR and 0.327 (-12.1%) in LCR. At G26 it decreased further to 0.303 (-10.2% from G14) in HCR and to 0.296 (-9.4%) in LCR. Note that the absolute values of H0 are influenced by the allele frequencies of the genotyped SNP markers, and it is the relative changes of H0 that reflect the altered genomic diversity. The observed rates of decrease are slower than in random mating populations. Using the known numbers of effective breeders at each generation we calculated the expected H0 at each generation assuming random mating, and found that the expected rate of decrease in H0 is on average 1.53 per generation for both HCR and LCR,. The observed rate of decrease, 0.95% and 0.97% per generation for HCR and LCR, respectively, is lower by 37-38%, consistent with the pedigree-based predictions (Figure 1) and confirming that the rotational breeding scheme has successfully reduced inbreeding as predicted [39].

thumbnail
Figure 5. Decrease of average heterozygosity over time in both lines.

Shown is the boxplot of genomewide average heterozygosity for the genotyped animals in two lines and three non-adjacent generations.

As the reduction of H0 over time primarily reflects higher levels of inbreeding, a majority of the increase of homozygote genotypes should be accounted for by the emergence or expansion of long runs of homozygosity (ROH). Using 10,185 SNPs in SNP Panel-1 (see Methods), we found that for HCR, ROH covered an average of 46% of the genome in G5 animals, and this rate increased to 54.8% in G26. For LCR, ROH covered 46.8% of the genome in G5, and 55% in G26. Thus the non-ROH regions shrink by 0.73-0.77% per generation in the two lines, accounting for most of the decrease of H0.

Linkage disequilibrium

We examined linkage disequilibrium (LD) patterns using the Panel-1 SNPs on Chromosome 1 (n = 978). The LD index r2 decays to 0.3 at ~3 Mb in both HCR and LCR (Figure 6). The level of LD is similar between HCR and LCR, showing slightly higher r2 in later generations, and is consistent with those reported for the NIH Heterogeneous Stock [40]. These results also suggest that the resolution of QTL mapping using HCR-LCR can be as high as 2-3 cM, considerably higher than the 20-40 cM resolution of F2 intercross of inbred lines [41].

thumbnail
Figure 6. Linkage disequilibrium (LD) decay over distance in HCR (A) and LCR (B) for chromosome 1.

The LD index r2, averaged for marker pairs falling in discrete distance bins, were plotted against the distance in Mb. LD decays to r=0.3 in about 3 Mb for both HCRs and LCRs.

"F2" intercross of HCR-LCR

As HCR and LCR have evolved separately, a direct between-line comparison of phenotypes and genotypes would incur the effect of population stratification. To create a QTL mapping population with randomized genomes, we performed "F2" intercross experiments using 28 HCR-LCR pairs and obtained 242 F1 animals (176 phenotyped). From the phenotyped F1 population we set up 63 mating pairs and produced 645 F2 animals. The term "F1" or “F2” is applied loosely in this context because our crosses were not based on inbred lines. However, we use F0, F1, and F2 to indicate the generations.

The running phenotype of F1 fell in an intermediate range between that of their HCR and LCR parents, and the F2 animals exhibited larger variations than the F1 rats (Figure 7). This pattern is consistent with the model in which most neutral alleles no longer co-segregate with the phenotypes; and at functional loci (i.e., those responsible for the phenotypic differences between HCRs and LCRs), F1s tend to be heterozygous and F2s carry a wider assortment of genotypes.

thumbnail
Figure 7. Distribution of running performance for animals of the F2 intercross experiment.

Shown are boxplots of the best running distance by generation and by sex.

Importantly, the h2 of the maximal running distance in the F0-F2 pedigree remained high, at 0.60 ± 0.05. The h2 for vertical work is comparable to those of the ancestral lines under selection (0.61 ± 0.05), while the h2 for body weight is lower (0.03 ± 0.04). In addition, we phenotyped a number of other physiological measures in a subset of F2 animals, including heart/body mass ratio (n=380), extensor digitorum longus (EDL) mass/body mass ratio (n=387), and percent body fat (n=490), which showed strong heritability (0.42 ± 0.11, 0.36 ± 0.10, and 0.48 ± 0.10, respectively). These results indicate that the running capacity and related traits are influenced by genetic factors, confirming that it is feasible to use the F2s as a mapping population for QTL.

Discussion

While previous studies have focused on functional, physiological comparisons between HCR and LCR rats, here we conducted the first in-depth pedigree and molecular genetic analysis of the two lines. Phenotypic data over the 28-generation pedigree not only revealed substantial heritability for the running capacity trait, but also showed that the heritability is maintained in later generations, suggesting continued selection. In addition, the strong heritability is recapitulated in the F2 intercross. These findings suggest that causal variants continue to segregate in both HCRs and LCRs and persist in the F2 population. Our study generates the first direct evidence that the trait under selection is highly heritable, providing justification for intercross-based QTL mapping. In addition, the heritability for vertical work, heart/body ratio, EDL/body ratio, and percent body fat suggests that the model can be used for simultaneous QTL mapping for multiple traits.

We observed continued response to selection in HCR during 20-28 generations. This observation is notable because it can be interpreted in two possible scenarios. The first is that the aerobic running capacity may be influenced by many interacting QTLs, and as variants in some loci become fixed under selection, the previously hidden phenotypic effects of other variants can be "released" and come under selection, thereby fueling prolonged responsiveness. This agrees with previous observations in similar systems that long-term selection did not exhaust the genetic variation for the selected trait due to the graduate shifts in the capacitors of cryptic genetic variation [42,43]. The second scenario is that the trait may be governed by many QTL of small effects, hence the strength of selection (~20% of animals become breeders in each generation (see Table 3)) may not have effectively driven the rapid changes of allele frequencies. The two scenarios are not mutually exclusive, and the observation means that not all causal alleles have been differentially fixed in the two lines. Therefore the F2 mapping approach needs to consider the possibility that the causal variants may be segregating within one or both lines.

HCRLCR
GenNNbNo. Male (Nm)No. Female (Nf)NNbNo. Male (Nm)No. Female (Nf)
0-261313-261313
1116281414125261313
2137281414126261313
3160281414167301515
4131261313138261313
5151281414150261313
6151301515155261313
7133301515141261313
8185281315143281414
9234301416198271413
10238291316150271314
11202341618245281414
12195301416202261313
13240341618216281414
14222301515209281414
15232321616231301515
16210381919141301515
17283422121254381919
18270462323246381919
19227422022238452223
20248442222233502525
21235502525226391920
22258482424204502525
23227482424191522626
24215442222172502525
25264502525225442222
26297442222240502525
27290---254---
28225---226---
Total59769674764915446920459461

Table 3. Number of phenotyped animals (N) by line and generation, and the number of animals chosen as effective breeders (Nb), separately shown for male (Nm) and female (Nf).

All the offspring generated including backup lines.
CSV
Download CSV

HCRs exhibited accelerated improvement of running capacity during G12-G15 (Figure 2). To identify the cause(s) of this acceleration we examined factors such as diet, running protocols, the breeding schedule, and "Operator", i.e., the experimenter or a team of experimenters assessing the running phenotype. The average litter size (i.e., fecundity) in the recorded HCR pedigree was not changed significantly during this period (Figure S4A), hence there was no noticeable change in fertility or the strength of selection (i.e., the fraction of animals chosen as breeders). There was no systematic correlation between litter size and inbreeding coefficient of the offspring (not shown); and there was no documented change in diet, running apparatus, or running protocol. The breeding schedules for the two lines were closely synchronized across all 28 generations (Figure S5A). The pedigree-based prediction of F was increasing in both HCR and LCR as expected (Figure 1). However, a more detailed retrospective analysis of the breeding records found three factors having changes during the G12-G15 period. The first is Operator: a team supervised by Operator 3 performed the running tests during G7-G13, while a team supervised by Operator 4 performed the tests during G14-G15. The second factor is the number of animals in the pedigree with no entries for running data which accrued mostly from rats that "refused" to run. The number of non-compliant rats in both lines gradually increased during G7-G13 in both lines, dropped immediately at G14, and remained low for most of subsequent generations (Figure S4B). Despite presumed standardization of the running protocol, the loss of running data may be Operator dependent in the sense that "refusal to run" is a subjective measure. Third, the fraction of mating pairs that were out-of-schedule increased in G12-G14 in HCR, and dropped after G15 (Figure S4C). The simplest interpretation of these co-occurrences is that Operator 3 subjectively determined that a large number of animals refused to run. Those who did run showed no improvement over G6-G13. With Operator 4, nearly all animals were able to run, and ran better than previous generations. While plausible, this simple scenario does not explain all the observations. First, despite being kept and tested under the same conditions as the HCRs, the LCRs exhibited no comparable acceleration or deceleration in running capacity. Second, the acceleration in HCRs began in G12-G13 with the unexplained emergence in some families of one or two exceptional runners, whose running distance were often more than twice as long as that of their siblings (Figure S6). The performance of these runners could not be easily explained by Operator. Partly because the exceptional runners tended to be selected as breeders, such improved performance spread wider across the cohort in G14-G15 and gradually became the norm after G16. However, there was not a clear-cut Mendelian segregation pattern in these generations: the pairing of two exceptional runners often still produced mediocre offspring. Among HCR mating pairs in G12-G15 there were 13 out-of-schedule pairs, which did not produce more exceptional runners than on-schedule pairs (not shown). The location of animal facility changed between G15 and G16 for both lines, but this change took place after the acceleration had started. Despite these complications, heritability estimates for HCR, when calculated for three-generation sections of the pedigree and shifted by one generation, did not show dramatic changes over the generations (Figure S5B).

The accelerated improvement of running capacity in HCRs during G12-G15 could also reflect genetic changes. However, emergence of a single high-impact de novo mutation is unlikely, as prodigious running capacity arose in multiple families concurrently. Such a pattern, however, is compatible with a scenario in which causal "high" alleles in multiple genes interact in a non-linear fashion. Various combinations of the high alleles could have undergone gradual enrichment and in G12-G15, began to manifest as improved phenotype when the most favored combinations were formed. Future studies, including linkage analysis of these intermediate generations, are needed to characterize the genetic changes accompanying the apparent varying tempo of trait evolution.

In the F2 generation of the intercross, the running distance distribution is wider than in F1, but did not reach the full range seen in F0 animals. The fact that none of the F2 animals performed as well as their HCR grandparents, and very few performed as poorly as the LCR grandparents, strongly suggests that multiple genetic loci are involved. The Castle-Wright estimator of the effective number of QTLs is calculated as 4-10 using our F2 data [44,45]. Caution should be taken as the calculation is based on simplifying assumptions such as unlinked loci of equal effects that have no interaction. Only the actual linkage or association studies can reveal the number and impact of QTL underlying the trait in question.

The HCR-LCR system was initiated in 1996 [22] and reached G28 in 2011. During this time, the two lines have diverged in innate endurance running capacity and showed marked differences in body type and metabolic traits. The HCR animals show a lower weight gain than LCR, in both young and adult rat, and this can be partly accounted for by higher spontaneous activity and lower fuel economy during activity in HCRs [32]. The two lines also diverged for many health indicators, with HCR showing a relative resistance to obesity, higher insulin sensitivity, lower blood pressure, improved lipid parameters, and enhanced longevity [30,31,34]. These phenotypes are of immense public health interest, as prevalence of diabetes, cardiovascular disorders, obesity, and metabolic syndrome is rising at an alarming rate and account for a major portion of disease burden worldwide [46]. The model system used in this study is ideally suited for elucidating the fundamental biology of metabolic health. Understanding the genetic architecture and molecular underpinnings of the remarkable HCR-LCR differences has the potential to provide new insights into the relationship between exercise capacity and metabolic health in humans.

Taken as a whole, the results presented above suggest that the HCR-LCR system is well-suited to serve as a novel model system for studying genome evolution under sustained selection and for dissecting the functional and genetic basis of polygenic traits. The model exhibits large phenotypic divergence, sustained heritability for a wide range of cardiovascular and metabolic traits, and maintained outbred character. The complete pedigree is known, with running phenotype for all animals already collected, and tissue sample for most breeders archived. Compared to inbred line-based gene mapping, our system offers some additional advantages. First, while the F2 generation could be subjected to conventional linkage mapping [47], the two lines have accumulated ~60 generations of historical recombination (~30 as the NIH Heterogeneous Stock, 28 generations of divergent selection, and F2 intercross). Consequently, animals in both lines carry fine-grained genomic mosaics of eight "ancestral" inbred strains, with LD structure on the order of 3 Mb, allowing for greater resolution in association analysis [4850]. Second, our system has maintained genetic diversity through rotational breeding, such that networks of interacting QTLs may have evolved jointly under selection, making the system particularly suitable for detection of interaction QTL [51]. Combining QTL mapping with the wealth of existing knowledge of the HCR-LCR system is expected to allow the identification and prioritization of high quality candidate genes that will shed insight into the biology of oxidative capacity and metabolic fitness.

Materials and Methods

Ethics statement

This study was approved by the University Committee on Use and Care of Animals, Ann Arbor, Michigan (Approval Numbers: #08905 and #03797). The proposed animal use procedures are in compliance with University guidelines, and State and Federal regulations.

Rotational breeding scheme

In practice, each line contains at least 13 mating pairs through all generations. From each of the family produced, one male and one female are selected as breeders for the subsequent generation. For HCR, the male and female with the greatest running distance are selected, whereas in LCR, those with the lowest distance are selected. The breeders are paired between different families to avoid brother-sister mating, and the pairings rotate in successive generations to minimize inbreeding [39]. When the 13th rotation is reached, same-family mating is skipped, and the pairings are reiterated starting again in the same way as rotation 1 (Figure S7A). Sometimes, if a particular mating fails, or if a family lacks animal of one sex, substitute mating is attempted involving a male from another family (Table 2). In some cases, one male is mated to two females. After G12, female HCR with extremely low body weights were not selected to be breeders in order to avoid reduced fecundity [52]. During G9-G13 in HCR there were also three cross-generation matings, whose offspring were incorporated into subsequent generations. Further, occasionally additional pairs are bred to generate experimental cohorts for study by us or for sharing with collaborators, and the progenies in these "analytical families" are not used for maintaining the lines, and are not counted in our calculation of the expected inbreeding levels (see below). They are, however, used to calculate heritability and the distribution of trait values.

One inevitable consequence of this breeding scheme is the mating of first-cousins at every half interval (Figure S7B). For example, at G7, every breeding pair, such as 1M-7F (a male from Family 1 and a female from Family 7) involves first cousins, because they are from 1F-7M and 7F-13M matings respectively, in G6, in which 7M and 7F are siblings. This results in a 6% spike in inbreeding values in G8 (Figure 1), and such a cyclic pattern continues in subsequent generations, resulting in spikes at G14, G20, and G26. The actual pedigree deviates from a perfectly executed breeding scheme due to the inclusion of substitute breeders, and the resulting estimates of inbreeding coefficient from the actual pedigree depart slightly from the expectations (Figure 1). As some breeding pairs were assembled to generate offspring for research use rather than line propagation, the "effective breeders" are those that contribute offspring who are also used as breeders, and do not include those whose offspring were used only for research (Table 3).

Running phenotype

Eleven week old animals are subjected to run-to-exhaustion tests without prior training, except for brief sessions of treadmill education during the week prior to the tests. The purpose of such education sessions is to familiarize the rats to the experimenters and the testing equipment and to ensure that each rat has the ability to achieve a minimal level of continual running for 5 minutes at least once, which constitutes the threshold performance necessary for inclusion in the actual running tests the following week. During education, the rats learn to keep running in order to avoid a mild shock (1.2 mA of current at 3 Hz) induced by the electrified grid located at the back of the treadmill. For all sessions the treadmill is set at a 15-degree upward slope.

During the run-to-exhaustion test, each rat was evaluated on five consecutive days (Mon-Fri) for G0-G16 and on three alternating days (Mon-Wed-Fri) for G17-28. Each trial starts at a velocity of 10 m/min, which increases by 1 m/min every 2 min until the rat reaches exhaustion. Exhaustion point is defined as the third time a rat can no longer keep pace with the treadmill and remains on the shock grid for two seconds rather than resuming running. At this point the rat is removed from the treadmill and weighed. For each rat, the best distance out of the multiple trials is taken as the best estimate of its intrinsic capacity, and used as the criterion for breeder selection. The vertical work during each trial is estimated using the equation:

work = (running distance) x (body weight) x (sin[15°]) x (9.8m/s2)/1000 in which the unit for work is joule (J=kg•m2/s2). Unit for running distance and body weight is meter and gram, respectively.

Phenotype distribution

To display phenotype distribution we produced violin plots (as shown in Figures 1, 2, and 3) using the vioplot function from the Vioplot package in R. To calculate the Spearman's rank correlation coefficient (ρ) between maximal running distance and body weight, separately for two sexes within each line, we used the cor and cor.test functions in R. The ρ values were calculated for each generation, and averaged over G1-G28.

Inbreeding coefficient and Heritability

We calculated the inbreeding coefficient (F) for each animal in the pedigree using the calcInbreeding function from the pedigree package in R [53]. The pedigree for earlier generations of NIH:H animals, i.e., those that preceded the founders of our lines, was not available. We are therefore limited to calculate the increase of inbreeding coefficients from those of the founders, effectively assuming they were unrelated, while in fact they were related according to the (unknown) breeding patterns in the preceding generations. To calculate the expected F under random mating, we used the equation Fn+1 = Fn+(Nf+Nm)/(8*Nf*Nm)-Fn(Nf+Nm)/(8*Nf*Nm), where the Fn and Fn+1 are the inbreeding coefficients at the n-th and (n+1)-th generation, respectively, and Nf and Nm are the numbers of male and female breeders at the n-th generation (Nf = Nm = 13 in every generation for a 13-family breeding scheme). To calculate the expected F under perfect adherence to the rotational scheme, we generated an idealized pedigree of 13 mating pairs of exact mating patterns as intended, and used the pedigree package to calculate F for every member of the pedigree.

To calculate the narrow-sense heritability (h2), we applied the variance and covariance component models as implemented in SOLAR version 4.3.1 [54]. We estimated h2 for maximal running distance both over the entire pedigree and for four-generation intervals (with one and three generation overlap) to assess the h2 variation over time. We also estimated h2 for body weight and work for the G0-G28 pedigree, and for additional traits for the F2 intercross. For the G0-G28 analysis, we included sex and operator as covariates; and for the F2 intercross we included sex and batch because the breeding was performed in two batches, containing 154 and 491 F2 animals, respectively.

Genotyping and data processing

DNA from 22-25 breeders from both lines in three non-adjacent generations (G5, G14, and G26, n=142) was extracted from frozen liver tissue, and genotyped across 10,846 SNP loci using the Affymetrix Rat Mapping 10K GeneChip. Attempts to extract DNA from generations earlier than G5 revealed that many samples in G0 and G4 were degraded. We therefore chose G5 as the earliest generation in our analysis due to its assured DNA quality. In assessing the quality of SNP markers we removed 28 duplicate SNPs, 496 SNPs with genotype missing rate >10%, and 137 SNPs with Hardy-Weinberg Equilibrium test p < 0.001. These steps led to 10,185 SNPs that formed the "Panel-1" markers. As some analyses require a reduced set of SNPs without rare variants and without strong linkage disequilibrium, we removed from Panel-1 an addition set of 7,284 SNPs selected by trimming SNP pairs in linkage disequilibrium with r2 value >0.05 (in windows of 10 SNPs, sliding by 2 SNPs each time), and 572 SNPs with minor allele frequency (MAF) <5%. After these steps, 2,518 SNPs remained and formed the "Panel-2" markers. The Panel-2 markers were used in calculations of IBD, AMOVA, and genome-wide average heterozygosity. Pairwise Identity-by-State (IBS) matrix was estimated in PLINK [55] using the -genome command and Panel-2 markers. Multidimensional scaling analysis of the IBS matrix was performed in R [56]. Analysis of molecular variance (AMOVA), as implemented in the program Arlequin, was used to calculate the within- and among-group differentiation [57].

We assessed the accuracy of recorded sex for each genotyped animal by calculating the average heterozygosity of the X chromosome SNPs. Male and female animals are confirmed by non-overlapping distributions of ChrX heterozygosity values. We confirmed known sibling pairs among the genotyped animals by plotting pairwise Z0 vs Z1 values in R. Z0 and Z1 values were determined in PLINK using the -genome command.

Runs of homozygosity

Using the 10,185 Panel-1 markers, we identified long runs of homozygosity (ROH), in PLINK using the -homozyg command. We defined ROHs as genomic segments with at least 4 homozygous markers and having a density of at least 1 SNP per 500Kb. Total ROH length in each animal was obtained by summing over all ROH and also reported as the fraction of the rat genome (2.75 Gb).

Genomewide average heterozygosity

Using the Panel-2 markers, we calculated the average heterozygosity in PLINK using the -hardy command, and compared across genotyped lines and generations using boxplots. The expected heterozygosity values were calculated using the equation Hn+1 = Hn(1-((Nf+Nm)/(8*Nf*Nm))), where the Hn and Hn+1 are the heterozygosity at the n-th and (n+1)-th generation, respectively, and Nf and Nm are the numbers of male and female breeders at the n-th generation, respectively (Nf = Nm = 13 in every generation for a 13-family breeding scheme) (Table 3).

LD calculation

Pairwise measurements of LD (r2) were calculated for marker pairs within 5 Mb on Chromosome 1 using Haploview [58]. Chromosome 1 was chosen as a representaive autosome. To show the relationship between r2 and inter-marker distance, we calculated average r2 values for groups of marker pairs falling in discrete bins of inter-marker distance, in 500Kb increments, and plotted the values for G5, G14, and G26 in both HCR and LCR (as shown in Figure 6).

“F2” intercross

We performed the F2 intercross in two batches. For the first batch, we randomly selected 4 males and 4 females from G26 of each line to form 8 HCR-LCR reciprocal pairs, which generated 79 F1 rats, from which 20 males and 20 females were randomly selected to form pairs between different F1 families (i.e., avoiding brother-sister mating). This generated 154 F2 rats. For the second batch, we selected 9 males and 9 females from G28 of each line to form 18 mating pairs, which generated 163 F1 rats, from which 97 were phenotyped. Out of the 97 F1 animals, 40 males and 40 females were selected across the 18 families with equal representation between HCR/LCR parentage, such that we had 4 combinations (HCR-mom male with HCR-mom female, HCR-mom male with LCR-mom female, LCR-mom male with HCR-mom female, and LCR-mom male with LCR-mom female) and generated 491 F2 rats. The two batches together yielded 645 F2 rats (Figure 7). Phenotyping for running performance followed the same protocols as described above. Additionally, for F2 animals in the second batch we measured lean mass, fat mass, fluid mass, fasting blood glucose, heart mass, and EDL muscle mass at 16-20 weeks of age. Body composition was determined via NMR using Bunter Optics Minispec LF90 II. Blood glucose after a 4-hour fast was determined using Accu-Check Aviva meter. At time of dissection, heart and EDL muscles were weighed immediately upon harvesting.

Supporting Information

Figure S1.

Distribution of maximal running distance for generations 0 to 28 for males (A) and females (B). Shown are "violin-plots" for individual generations for females and males separately. The blue tick marks on the y axis indicate the maximal running distance for eleven inbred lines, which are ordered, from top to bottom for males (A) as DA (968m), PVG (754m), SR (615m), AUG (594m), ACI (447m), LEW (405m), WKY (387m), BUF (355m), F344 (332m), MNS (302m) and COP (262m), and for females (B) as AUG (805m), DA (712m), PVG (682m), F344 (606m), LEW (479m), ACI (453m), WKY (441m), SR (409m), BUF (391m), COP (333m), and MNS (315m).

https://doi.org/10.1371/journal.pone.0077588.s001

(TIFF)

Figure S2.

Narrow-sense heritability of running capacity remained positive over time. Shown are heritability (h2) estimates and standard errors for maximal running distance in four-generation intervals that overlap by one-generation for HCR (A) and LCR (B).

https://doi.org/10.1371/journal.pone.0077588.s002

(TIFF)

Figure S3.

Sample quality assessment in 10K SNP genotype data. (A) Average heterozygosity over 61 SNPs on the X chromosome for 142 genotyped animals, ordered by line and generation. Males and females fall in two non-overlapping clusters, indicating that the observed X-chromosome heterozygosity is consistent with the reported sex of the animals. (B) Scatter plot of Z0 (genomic proportion that a pair of animals share 0 allele identical-by-descent [IBD]) versus Z1 (proportion that a pair of animals share 1 allele IBD), showing that known sib pairs form a separate cluster than more distant relatives.

https://doi.org/10.1371/journal.pone.0077588.s003

(TIFF)

Figure S4.

Assessment of litter size, missing phenotype, and rotational breeding schedule over the generations. (A) Average litter size over G1-G28 for the HCR animals in the recorded pedigree, accompanied by standard deviations among all the families for each generation. No significant change in litter size is observed over the course of selection. (B) Number of animals in each generation without recorded running phenotype for HCR (blue) and LCR (red). The increase in missing phenotypes between G7-G13 overlaps operator 3 (shown as the orange bar below). (C) Percent of mating pairs out-of-schedule (off-rotation) per generation for HCR (green) and LCR (red).

https://doi.org/10.1371/journal.pone.0077588.s004

(TIFF)

Figure S5.

Between-line synchrony and short-term heritability estimates in HCR. (A) Dates of birth (x-axis) for G1-G28 animals (y-axis). The horizontal bars indicate the range of birth dates, and the dots indicate the average. The close match between HCR and LCR shows that the two lines are synchronized. (B) HCR narrow-sense heritability for adjacent 3-generation intervals, showing no significant change over time.

https://doi.org/10.1371/journal.pone.0077588.s005

(TIFF)

Figure S6.

Running performance of mid-parents and offspring for on/off-rotation mating types for G12 (A), G13 (B), G14 (C), and G15 (D). Shown are the scatterplots of mid-parent running distance (x-axis) versus male and female offspring running distance (y-axis) for on-rotation mating (indicated by XX.0), off-rotation mating due to mother (XX.1), off-rotation mating due to father (XX.2), off-rotation mating due to both mother and father (XX.3), and either father or mother from the enrichment cohort (XX.4). Data points above the solid diagonal line indicate offspring with better running performance than their mid-parent.

https://doi.org/10.1371/journal.pone.0077588.s006

(TIFF)

Figure S7.

Rotational breeding scheme. (A) Mate pairing matrix for breeding rotations 1 through 12, where females from families 1 through 13 are designated in columns (from left to right), while males are designated in rows (from top to bottom). Male-female pairs are formed differently in successive generations as indicated by the rotation numbers in the matrix. (B) An example to show how, at the mid-cycle through rotational breeding among 13 families, every breeding pair is a first-cousin mating. The example shown is rotation 7 between Family 1 and Family 7. A male from family 1 (1M) to be mated with a female from family 7 (7F) are both offspring of breeding members of family 7 in rotation 6, thus making them first cousins.

https://doi.org/10.1371/journal.pone.0077588.s007

(TIFF)

File S1.

HCR pedigree in G1-G28 and select phenotypes.

https://doi.org/10.1371/journal.pone.0077588.s008

(XLSX)

File S2.

LCR pedigree in G1-G28 and select phenotypes.

https://doi.org/10.1371/journal.pone.0077588.s009

(XLSX)

Acknowledgments

We acknowledge Thomas Hancock, Christopher Chermside-Scabbo, and Abigail Burant, for assistance in phenotyping F1 and F2 intercrosses. We also acknowledge Frank Albert (Princeton University) and Alexander Cagan (Max Planck Institute for Evolution Anthropology in Leipzig, Germany), for advice and stimulating discussions.

Author Contributions

Conceived and designed the experiments: YR KAO NRQ MKT LH MK LGK SLB CFB JZL. Performed the experiments: YR KAO NRQ MKT LH MK LGK SLB. Analyzed the data: YR LGK SLB CFB JZL. Contributed reagents/materials/analysis tools: YR KAO NRQ MKT LH MK LGK SLB. Wrote the manuscript: YR LGK SLB CFB JZL.

References

  1. 1. Blair SN, Kampert JB, Kohl HW 3rd, Barlow CE, Macera CA et al. (1996) Influences of cardiorespiratory fitness and other precursors on cardiovascular disease and all-cause mortality in men and women. JAMA 276: 205-210. doi:https://doi.org/10.1001/jama.276.3.205. PubMed: 8667564.
  2. 2. Blair SN, Kohl HW 3rd, Paffenbarger RS Jr., Clark DG, Cooper KH et al. (1989) Physical fitness and all-cause mortality. A prospective study of healthy men and women. JAMA 262: 2395-2401. doi:https://doi.org/10.1001/jama.1989.03430170057028. PubMed: 2795824.
  3. 3. Church TS, Cheng YJ, Earnest CP, Barlow CE, Gibbons LW et al. (2004) Exercise capacity and body composition as predictors of mortality among men with diabetes. Diabetes Care 27: 83-88. doi:https://doi.org/10.2337/diacare.27.1.83. PubMed: 14693971.
  4. 4. Gulati M, Pandey DK, Arnsdorf MF, Lauderdale DS, Thisted RA et al. (2003) Exercise capacity and the risk of death in women: the St James Women Take Heart Project. Circulation 108: 1554-1559. doi:https://doi.org/10.1161/01.CIR.0000091080.57509.E9. PubMed: 12975254.
  5. 5. Kokkinos P, Myers J, Kokkinos JP, Pittaras A, Narayan P et al. (2008) Exercise capacity and mortality in black and white men. Circulation 117: 614-622. doi:https://doi.org/10.1161/CIRCULATIONAHA.107.734764. PubMed: 18212278.
  6. 6. Lakka TA, Venäläinen JM, Rauramaa R, Salonen R, Tuomilehto J et al. (1994) Relation of leisure-time physical activity and cardiorespiratory fitness to the risk of acute myocardial infarction. N Engl J Med 330: 1549-1554. doi:https://doi.org/10.1056/NEJM199406023302201. PubMed: 8177243.
  7. 7. Myers J, Prakash M, Froelicher V, Do D, Partington S et al. (2002) Exercise capacity and mortality among men referred for exercise testing. N Engl J Med 346: 793-801. doi:https://doi.org/10.1056/NEJMoa011858. PubMed: 11893790.
  8. 8. Sandvik L, Erikssen J, Thaulow E, Erikssen G, Mundal R et al. (1993) Physical fitness as a predictor of mortality among healthy, middle-aged Norwegian men. N Engl J Med 328: 533-537. doi:https://doi.org/10.1056/NEJM199302253280803. PubMed: 8426620.
  9. 9. Sui X, LaMonte MJ, Laditka JN, Hardin JW, Chase N et al. (2007) Cardiorespiratory fitness and adiposity as mortality predictors in older adults. JAMA 298: 2507-2516. doi:https://doi.org/10.1001/jama.298.21.2507. PubMed: 18056904.
  10. 10. Leeper NJ, Myers J, Zhou M, Nead KT, Syed A et al. (2012) Exercise capacity is the strongest predictor of mortality in patients with peripheral arterial disease. J Vasc Surg, 57: 728–33. PubMed: 23044259.
  11. 11. Pedersen BK, Saltin B (2006) Evidence for prescribing exercise as therapy in chronic disease. Scand J Med Sci Sports 16 Suppl 1: 3-63. doi:https://doi.org/10.1111/j.1600-0838.2006.00520.x. PubMed: 16451303.
  12. 12. Bray MS, Hagberg JM, Pérusse L, Rankinen T, Roth SM et al. (2009) The human gene map for performance and health-related fitness phenotypes: the 2006-2007 update. Med Sci Sports Exerc 41: 35-73. doi:https://doi.org/10.1249/01.mss.0000353371.87524.43. PubMed: 19123262.
  13. 13. Bouchard C, Rankinen T, Chagnon YC, Rice T, Pérusse L et al. (2000) Genomic scan for maximal oxygen uptake and its response to training in the HERITAGE Family Study. J Appl Physiol 88: 551-559. PubMed: 10658022.
  14. 14. Klissouras V (1971) Heritability of adaptive variation. J Appl Physiol 31: 338-344. PubMed: 5165174.
  15. 15. Engström LM, Fischbein S (1977) Physical capacity in twins. Acta Genet Med Gemellol (Roma) 26: 159-165. PubMed: 563658.
  16. 16. Howald H (1976) Ultrastructure and biochemical function of skeletal muscle in twins. Ann Hum Biol 3: 455-462. doi:https://doi.org/10.1080/03014467600001721. PubMed: 988780.
  17. 17. Bray MS (2000) Genomics, genes, and environmental interaction: the role of exercise. J Appl Physiol 88: 788-792. PubMed: 10658051.
  18. 18. Kaprio J, Koskenvuo M, Sarna S (1981) Cigarette smoking, use of alcohol, and leisure-time physical activity among same-sexed adult male twins. Prog Clin Biol Res 69: 37–46 Pt C: 37-46 PubMed: 7198250.
  19. 19. Pérusse L, Tremblay A, Leblanc C, Bouchard C (1989) Genetic and environmental influences on level of habitual physical activity and exercise participation. Am J Epidemiol 129: 1012-1022. PubMed: 2705422.
  20. 20. Prud'homme D, Bouchard C, Leblanc C, Landry F, Fontaine E (1984) Sensitivity of maximal aerobic power to training is genotype-dependent. Med Sci Sports Exerc 16: 489-493. doi:https://doi.org/10.1249/00005768-198410000-00012. PubMed: 6542620.
  21. 21. Simoneau JA, Lortie G, Boulay MR, Marcotte M, Thibault MC et al. (1986) Inheritance of human skeletal muscle and anaerobic capacity adaptation to high-intensity intermittent training. Int J Sports Med 7: 167-171. doi:https://doi.org/10.1055/s-2008-1025756. PubMed: 3733313.
  22. 22. Koch LG, Britton SL (2001) Artificial selection for intrinsic aerobic endurance running capacity in rats. Physiol Genomics 5: 45-52. PubMed: 11161005.
  23. 23. Hansen C, Spuhler K (1984) Development of the National Institutes of Health genetically heterogeneous rat stock. Alcohol Clin Exp Res 8: 477-479. doi:https://doi.org/10.1111/j.1530-0277.1984.tb05706.x. PubMed: 6391259.
  24. 24. Koch LG, Britton SL, Wisløff U (2012) A rat model system to study complex disease risks, fitness, aging, and longevity. Trends Cardiovasc Med 22: 29-34. doi:https://doi.org/10.1016/j.tcm.2012.06.007. PubMed: 22867966.
  25. 25. Lujan HL, Britton SL, Koch LG, DiCarlo SE (2006) Reduced susceptibility to ventricular tachyarrhythmias in rats selectively bred for high aerobic capacity. Am J Physiol Heart Circ Physiol 291: H2933-H2941. doi:https://doi.org/10.1152/ajpheart.00514.2006. PubMed: 16891405.
  26. 26. Thyfault JP, Rector RS, Uptergrove GM, Borengasser SJ, Morris EM et al. (2009) Rats selectively bred for low aerobic capacity have reduced hepatic mitochondrial oxidative capacity and susceptibility to hepatic steatosis and injury. J Physiol 587: 1805-1816. doi:https://doi.org/10.1113/jphysiol.2009.169060. PubMed: 19237421.
  27. 27. Muncey AR, Saulles AR, Koch LG, Britton SL, Baghdoyan HA et al. (2010) Disrupted sleep and delayed recovery from chronic peripheral neuropathy are distinct phenotypes in a rat model of metabolic syndrome. Anesthesiology 113: 1176-1185. doi:https://doi.org/10.1097/ALN.0b013e3181f56248. PubMed: 20938334.
  28. 28. Burghardt PR, Flagel SB, Burghardt KJ, Britton SL, Gerard-Koch L et al. (2011) Risk-assessment and coping strategies segregate with divergent intrinsic aerobic capacity in rats. Neuropsychopharmacology 36: 390-401. doi:https://doi.org/10.1038/npp.2010.144. PubMed: 20927049.
  29. 29. Wikgren J, Mertikas GG, Raussi P, Tirkkonen R, Äyräväinen L et al. (2012) Selective breeding for endurance running capacity affects cognitive but not motor learning in rats. Physiol Behav 106: 95-100. doi:https://doi.org/10.1016/j.physbeh.2012.01.011. PubMed: 22285210.
  30. 30. Wisløff U, Najjar SM, Ellingsen O, Haram PM, Swoap S et al. (2005) Cardiovascular risk factors emerge after artificial selection for low aerobic capacity. Science 307: 418-420. doi:https://doi.org/10.1126/science.1108177. PubMed: 15662013.
  31. 31. Noland RC, Thyfault JP, Henes ST, Whitfield BR, Woodlief TL et al. (2007) Artificial selection for high-capacity endurance running is protective against high-fat diet-induced insulin resistance. Am J Physiol Endocrinol Metab 293: E31-E41. doi:https://doi.org/10.1152/ajpendo.00500.2006. PubMed: 17341547.
  32. 32. Novak CM, Escande C, Burghardt PR, Zhang M, Barbosa MT et al. (2010) Spontaneous activity, economy of activity, and resistance to diet-induced obesity in rats bred for high intrinsic aerobic capacity. Horm Behav 58: 355-367. doi:https://doi.org/10.1016/j.yhbeh.2010.03.013. PubMed: 20350549.
  33. 33. Lessard SJ, Rivas DA, Chen ZP, van Denderen BJ, Watt MJ et al. (2009) Impaired skeletal muscle beta-adrenergic activation and lipolysis are associated with whole-body insulin resistance in rats bred for low intrinsic exercise capacity. Endocrinology 150: 4883-4891. doi:https://doi.org/10.1210/en.2009-0158. PubMed: 19819977.
  34. 34. Koch LG, Kemi OJ, Qi N, Leng SX, Bijma P et al. (2011) Intrinsic aerobic capacity sets a divide for aging and longevity. Circ Res 109: 1162-1172. doi:https://doi.org/10.1161/CIRCRESAHA.111.253807. PubMed: 21921265.
  35. 35. Farid A, Makarechian M, Strobeck C (1987) Inbreeding under a Cyclical Mating System. Theor Appl Genet 73: 506-515. doi:https://doi.org/10.1007/BF00289187.
  36. 36. Windig JJ, Kaal L (2008) An effective rotational mating scheme for inbreeding reduction in captive populations illustrated by the rare sheep breed Kempisch Heideschaap. Animals 2: 1733-1741. PubMed: 22444078.
  37. 37. Barbato JC, Koch LG, Darvish A, Cicila GT, Metting PJ et al. (1998) Spectrum of aerobic endurance running performance in eleven inbred strains of rats. J Appl Physiol 85: 530-536. PubMed: 9688730.
  38. 38. Excoffier L, Smouse PE, Quattro JM (1992) Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics 131: 479-491. PubMed: 1644282.
  39. 39. Nomura T, Yonezawa K (1996) A comparison of four systems of group mating for avoiding inbreeding. Genet Sel Evol 28: 141-159. doi:https://doi.org/10.1186/1297-9686-28-2-141.
  40. 40. Johannesson M, Lopez-Aumatell R, Stridh P, Diez M, Tuncel J et al. (2009) A resource for the simultaneous high-resolution mapping of multiple quantitative trait loci in rats: the NIH heterogeneous stock. Genome Res 19: 150-158. PubMed: 18971309.
  41. 41. Flint J, Valdar W, Shifman S, Mott R (2005) Strategies for mapping and cloning quantitative trait genes in rodents. Nat Rev Genet 6: 271-286. doi:https://doi.org/10.1038/ni1170. PubMed: 15803197.
  42. 42. Dworkin I, Palsson A, Birdsall K, Gibson G (2003) Evidence that Egfr contributes to cryptic genetic variation for photoreceptor determination in natural populations of Drosophila melanogaster. Curr Biol 13: 1888-1893. doi:https://doi.org/10.1016/j.cub.2003.10.001. PubMed: 14588245.
  43. 43. Queitsch C, Sangster TA, Lindquist S (2002) Hsp90 as a capacitor of phenotypic variation. Nature 417: 618-624. doi:https://doi.org/10.1038/nature749. PubMed: 12050657.
  44. 44. Castle WE (1921) An Improved Method of Estimating the Number of Genetic Factors Concerned in Cases of Blending Inheritance. Science 54: 223. doi:https://doi.org/10.1126/science.54.1393.223. PubMed: 17792870.
  45. 45. Wright S (1968) Evolution and the genetics of populations; a treatise. Chicago:: University of Chicago Press. v. p.
  46. 46. James PT, Rigby N, Leach R (2004) The obesity epidemic, metabolic syndrome and future prevention strategies. Eur J Cardiovasc Prev Rehabil 11: 3-8. doi:https://doi.org/10.1097/01.hjr.0000114707.27531.48. PubMed: 15167200.
  47. 47. Albert FW, Carlborg O, Plyusnina I, Besnier F, Hedwig D et al. (2009) Genetic architecture of tameness in a rat model of animal domestication. Genetics 182: 541-554. doi:https://doi.org/10.1534/genetics.109.102186. PubMed: 19363126.
  48. 48. Mott R, Talbot CJ, Turri MG, Collins AC, Flint J (2000) A method for fine mapping quantitative trait loci in outbred animal stocks. Proc Natl Acad Sci U S A 97: 12649-12654. doi:https://doi.org/10.1073/pnas.230304397. PubMed: 11050180.
  49. 49. Valdar W, Flint J, Mott R (2006) Simulating the collaborative cross: power of quantitative trait loci detection and mapping resolution in large sets of recombinant inbred strains of mice. Genetics 172: 1783-1797. PubMed: 16361245.
  50. 50. Valdar WS, Flint J, Mott R (2003) QTL fine-mapping with recombinant-inbred heterogeneous stocks and in vitro heterogeneous stocks. Mamm Genome 14: 830-838. doi:https://doi.org/10.1007/s00335-003-3021-1. PubMed: 14724734.
  51. 51. Carlborg O, Jacobsson L, Ahgren P, Siegel P, Andersson L (2006) Epistasis and the release of genetic variation during long-term selection. Nat Genet 38: 418-420. doi:https://doi.org/10.1038/ng1761. PubMed: 16532011.
  52. 52. Bünger L, Hill WG (1999) Inbred lines of mice derived from long-term divergent selection on fat content and body weight. Mamm Genome 10: 645-648. doi:https://doi.org/10.1007/s003359901063. PubMed: 10341104.
  53. 53. Bastiaansen JW, Coster A, Calus MP, van Arendonk JA, Bovenhuis H (2012) Long-term response to genomic selection: effects of estimation method and reference population structure for different genetic architectures. Genet Sel Evol 44: 3. doi:https://doi.org/10.1186/1297-9686-44-3. PubMed: 22273519.
  54. 54. Almasy L, Blangero J (1998) Multipoint quantitative-trait linkage analysis in general pedigrees. Am J Hum Genet 62: 1198-1211. doi:https://doi.org/10.1086/301844. PubMed: 9545414.
  55. 55. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR et al. (2007) PLINK: A tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 81: 559-575. doi:https://doi.org/10.1086/519795. PubMed: 17701901.
  56. 56. R Development Core Team (2010) R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing.
  57. 57. Excoffier L, Laval G, Schneider S (2005) Arlequin (version 3.0): an integrated software package for population genetics data analysis. Evol Bioinform Online 1: 47-50. PubMed: 19325852.
  58. 58. Barrett JC, Fry B, Maller J, Daly MJ (2005) Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 21: 263-265. doi:https://doi.org/10.1093/bioinformatics/bth457. PubMed: 15297300.