In this review, we have summarized our work using combined complex statistical genetics, bioinformatics, and functional genomics to determine the genetic basis of the age-related thymic involution in C57BL/6J X DBA/2J recombinant inbred mice and the parental B6 and D2 mice. We have shown that these mice provided a valuable genetic model that can permit resampling of thymuses from different aged but genetically identical animals and determination of the relative significance of age-associated changes in the thymus. Our results suggest that the quantitative trait loci (QTL) regulating the Con A-induced thymocyte proliferative response were mapped to mouse chromosome Chr 11 (D11Mit51 at 18 cM), a region that harbors the IL-12b gene. The importance of IL-12b in maintaining thymic integrity and function during the aging process was confirmed by a more rapid involution of the thymus in IL-12b knockout (IL-12b-/-) mice compared to wild-type (WT) mice. Functionally, IL-12 provided a strong synergistic effect to augment the IL-7 or IL-2 induced thymocyte proliferative response, especially in both aged WT and IL-12b-/- mice, but not in normal young mice. In contract to the proliferative response, the age-related decline in the total number of thymocytes was determined at different age, and mapped to loci on Chr 9, 62 cM and Chr 10, 32 cM. Using matrix-assisted laser desorption/ionisation-time of flight-mass spectrometry (MALDI-TOF-MS), increased expression of peroxiredoxin was found to be correlated with thymic involution. Our results suggest the possibility to identify the complex molecular network that can be associated with the regulation of thymic involution in aged mice using a high-dimensional functional genomics approach.