Literature
首页医源资料库在线期刊美国临床营养学杂志2004年80卷第6期

Genomic and proteomic approaches for probing the role of vitamin D in health

来源:《美国临床营养学杂志》
摘要:RecentdevelopmentsingenomicsandproteomicshaveprovidednewopportunitiesforustoidentifymoleculartargetsofvitaminDaction。Cantheseapproachesbeusedtohelpclarifywhatconstitutesoptimalserumconcentrationsof25-hydroxyvitaminD。Ipresentanoverviewofhowproteomicand......

点击显示 收起

James C Fleet

1 From the Interdepartmental Nutrition Program and Department of Foods and Nutrition, Purdue University, West Lafayette, IN.

2 Presented at the conference "Vitamin D and Health in the 21st Century: Bone and Beyond," held in Bethesda, MD, October 9–10, 2003.

3 Supported by a grant from the NIH (grant DK54111).

4 Address correspondence to JC Fleet, 700 West State Street, Purdue University, West Lafayette, IN 47907-2059. E-mail: fleetj{at}cfs.purdue.edu.

ABSTRACT

Although we have learned a great deal about vitamin D metabolism and function since it first became apparent that this factor was required for bone health, there are still many gaps in our understanding, at both the basic science (eg, the molecular actions and targets of vitamin D) and applied (eg, what "adequate" vitamin D status means) levels. For example, although the identification of extrarenal 1-hydroxylase activity suggests that autocrine/paracrine actions of 1,25-dihydroxyvitamin D complement the classic endocrine actions of the hormone, the practical implications of this finding are only now being explored. In addition, studies showed that 1,25-dihydroxyvitamin D rapidly activates signal transduction pathways in addition to the classic transcriptional activation pathways that require the vitamin D receptor. These new modes of vitamin D action may be crucial to our understanding of both the traditional calcium-regulating actions of vitamin D and the anticancer actions of this essential mediator. Recent developments in genomics and proteomics have provided new opportunities for us to identify molecular targets of vitamin D action. Cancer researchers have demonstrated that these methods have utility for identifying useful biomarkers of disease states. Can these approaches be used to help clarify what constitutes optimal serum concentrations of 25-hydroxyvitamin D? I present an overview of how proteomic and genomic evaluations of cells, animals, and human subjects have been and can be used to improve our understanding of vitamin D biological processes and the role of vitamin D in health.

Key Words: Proteomics • genomics • biomarker

INTRODUCTION

Traditionally, vitamin D metabolism has been viewed as an endocrine system that responds to changes in serum calcium concentrations (1). Low dietary calcium intake is reflected by a decrease in serum calcium concentrations, which is in turn a signal for the increased production and release of parathyroid hormone (PTH). Among its functions, PTH stimulates renal 1-hydroxylase activity, leading to increased conversion of 25-hydroxyvitamin D3 [25(OH)D3] to 1,25-dihydroxyvitamin D3 [1,25(OH)2D3]. Elevated serum 1,25(OH)2D3 concentrations then stimulate the expression of vitamin D-responsive genes within the primary target tissues that control calcium homeostasis (ie, TRPV6 and calbindin D9k in intestine, osteocalcin and RANKL in bone, and TRPV5 and calbindin D28k in kidney) through activation of the vitamin D receptor (VDR). Although this system is clearly functional and biologically important during periods of calcium stress, it may not be sufficient to explain all of the biological actions of vitamin D.

As several recent reports and other reviews in this conference demonstrated, improved bone health and cancer chemoprevention may be more closely related to changes in serum 25(OH)D3 concentrations than to serum 1,25(OH)2D3 concentrations. For example, increases in serum 25(OH)D3 concentrations were associated with both maximal suppression of PTH (a proresorptive agent) (2–4) and increased efficiency of calcium absorption (5–8). There is also evidence that regular, high-dose, vitamin D supplementation decreased factures rates for common osteoporotic sites (9). This finding and other data suggest that local production of 1,25(OH)2D3, rather than endocrine signaling attributable to renal production, is critical for optimal bone health and cancer prevention (10–12). In addition, as we have come to understand more about the details of the molecular mechanisms controlling VDR function (13, 14), researchers have identified 1,25(OH)2D3 as an activator of various signal transduction pathways leading to the stimulation of protein kinases such as Src kinase, protein kinase C, protein kinase A, and the mitogen-activated protein kinases (15). These examples suggest that strict adherence to the accepted concepts of vitamin D biological mechanisms and actions are not likely to explain fully the health benefits of vitamin D.

In light of these and other recent findings, I contend that the field of vitamin D research would be well served by the use of several new technologies in an attempt to better describe the full range of vitamin D actions. The advances in technologies that permit whole-transcriptome analysis and large-scale proteomic analysis permit us to conduct unbiased evaluations of physiologic states and responses to treatments and interventions. In this review, I briefly demonstrate several instances in which genomic or proteomic analysis has expanded our understanding of vitamin D actions, I present a paradigm for future discovery-based research, and I demonstrate how these and other tools might be used to better define optimal vitamin D status.

GENOMIC AND PROTEOMIC APPROACHES: WHAT DO THEY OFFER TO THE FIELD OF VITAMIN D RESEARCH?

General approach
In the past decade, advances in the areas of genomics and proteomics created a revolution in science. These approaches are reviewed in detail elsewhere (16, 17) but are summarized briefly in Figure 1. Our traditional approach to understanding biological processes has been reductionist. We look for proteins that modulate functions of cells and tissues, and we study the regulation of their production and activity. There is no question that this approach has been very fruitful and will continue to be so. However, while this approach focuses our attention on testable hypotheses, it also provides virtual blinders, ie, we look only for the things we have already been studying. Moving beyond this approach can lead to dramatic advances in our understanding of scientific areas. For example, studies of the known proteins involved in iron metabolism failed to identify the cause of the iron-overload disease hemochromatosis. Only after the gene that is mutated in hemochromatosis was identified through a large-scale sequencing effort (genomics) did researchers find that the mutated gene encoded a protein with features similar to the major histocompatibility complex proteins, rather than a protein that was previously thought to be involved in iron metabolism (18). Similarly, studies have revealed that bone mass and metabolism can be regulated by genes that we would classically associate with obesity, such as leptin (19) and LDL receptor-related protein 5 (20). "Omic" approaches give us the opportunity to see beyond our expectations and discover new gene targets for vitamin D actions and functions.


View larger version (45K):
FIGURE 1.. Comparison of reductionist and genomic/globalist approaches to elucidating biological processes.

 
Studies identifying potential, new, vitamin D-regulated targets with genomics
Several transcript-profiling studies with microarrays have been conducted to elucidate the biological role of 1,25(OH)2D3. Because of the high cost of microarray-based gene expression profiling, most of the studies reported in the literature used very restricted conditions [eg, a single dose and time of treatment with 1,25(OH)2D3] and did not include replicates. This limits the utility of the approach and decreases our confidence in the results. However, more carefully controlled experiments were conducted and addressed some of the experimental design points noted in Figure 2. Figure 2A demonstrates that, depending on the time point examined after 1,25(OH)2D3 treatment and the time course of primary responses to the treatment, the changes in transcript levels could include both primary responses (eg, those resulting from direct, VDR-mediated interactions with gene promoters) and downstream responses that are indirect. Detection of a change in expression does not prove that a gene is a direct vitamin D target gene.


View larger version (33K):
FIGURE 2.. Experimental issues relevant to the use of genomic approaches to study vitamin D biological processes. A: Temporal relationship between vitamin D treatment and the molecular response at the level of the cell. B: Use of multiple independent perturbations to confirm the molecular actions of vitamin D. 1 OHase, 1-hydroxylase.

 
This phenomenon of differential patterns of responses to vitamin D is evident in the work of Lin et al (21). Those authors examined the time course of the response to 100 nmol/L 1,25(OH)2D3 or the vitamin D analog EB1089 in squamous cell carcinoma cells (SCC25) and identified 152 genes as being regulated by vitamin D (89 up and 63 down), with the use of the Affymetrix FL array (Affymetrix, Santa Clara, CA) and a 2.5-fold cutoff for determining a meaningful change in expression. Clustering was performed on the basis of the pattern of expression or the functional classification of the transcripts. Figure 3 shows the diversity of the vitamin D responses in these cells. Even within the genes with documented, functional, vitamin D response elements (VDREs) in their promoters, there was heterogeneity in responses. For example, the CYP24 transcript levels were rapidly increased by vitamin D treatment (significantly increased in 1 h), whereas osteopontin transcript levels increased more slowly in response to treatment (maximal expression at 12 h). This suggests that similar VDREs are differentially regulated, depending on the promoter context (22), but it also demonstrates the difficulty of discerning a direct transcriptional response solely on the basis of a time course, ie, even later responses can be attributable to direct effects. Another interesting finding from that study was that the vitamin D-induced responses in the transcript profile were much more diverse than might have been predicted previously. For example, several transcripts coding for proteins involved in protection from oxidative stress were gradually up-regulated by the vitamin D analog, including glucose-6-phosphate dehydrogenase (generating NADPH), glutathione peroxidase, and selenoprotein P. In addition, the thioredoxin reductase transcript was increased by 1 h after treatment, with peak induction by 6 h. Rapid suppression of transcripts for a variety of signaling peptides (eg, PTH-related protein and galanin) and induction of intracellular cell signaling proteins (eg, Cox-2, phosphoinositide-3-kinase, and p85 subunit) were also observed after treatment. It is not clear which of these responses is primary; none of these genes was previously shown to be regulated by vitamin D or to contain a functional VDRE. However, because 1,25(OH)2D3 promotes cellular differentiation, the up-regulation of some transcripts may represent a vitamin D-induced shift to a more differentiated phenotype. In any case, these data suggest that the traditional approach of examining only the expression of genes controlling cell cycle proteins in an attempt to explain the prodifferentiating action of vitamin D may provide limited information regarding the biological mechanisms of 1,25(OH)2D3 actions in proliferating or cancer cells.


View larger version (20K):
FIGURE 3.. Summary of patterns observed in response to treatment of squamous carcinoma cells with calcitriol or the vitamin D analog EB1089. The symbols represent the 10 distinct groups of up-regulated (5 groups) (black symbols) and down-regulated (5 groups) (white symbols) transcripts. These data demonstrate that there is considerable diversity in the temporal response of transcripts to vitamin D treatment. Adapted from Lin et al (21).

 
To identify direct vitamin D actions, scientists often conduct multiple complementary experiments and compare the results for consistency. In Figure 2B, 3 complementary mouse experiments are illustrated. Treatment of normal mice with 1,25(OH)2D3 would be expected to up-regulate or down-regulate specific genes in a target tissue, whereas the examination of target tissue transcript profiles in mice that lack essential components of the vitamin D signaling system (eg, VDR or 1-hydroxylase) would be expected to demonstrate opposite effects on vitamin D target genes [eg, a gene that is activated with 1,25(OH)2D3 injection would be down-regulated in VDR- or 1-hydroxylase-null mice]. A preliminary attempt at this approach with a small group of animals was recently reported by Li et al (23). By comparing the gene expression profile changes that occurred in kidney after 1,25(OH)2D3 injection with those that resulted from loss of the VDR (wild-type mice compared with VDR knockout mice), those authors identified 95 genes for which the response attributable to vitamin D injection was the opposite of the response attributable to loss of the VDR. Twenty-eight of those transcripts (including 1-hydroxylase mRNA) were up-regulated in VDR-null mice and down-regulated in vitamin D-treated mice, whereas 67 of the transcripts (including 24-hydroxylase mRNA) were down-regulated in VDR-null mice and up-regulated in vitamin D-treated mice. Like the study by Lin et al (21), this study identified many potential vitamin D target genes. Although neither of these studies definitively identified new targets, they narrowed the list of candidates considerably and provided clear guidance for investigators who wish to conduct careful reductionist experiments involving these genes.

Studies identifying new protein complexes necessary for vitamin D actions
Very few studies have taken a proteomics approach to the examination of vitamin D actions in cells. One notable study was recently conducted with keratinocytes by Oda et al (24). Previously, Rachez et al (25) developed an in vitro assay to assess the complex of proteins that associates with the VDR during vitamin D-mediated gene transcription. By using a VDR ligand-binding domain–glutathione S-transferase fusion protein, they were able to identify the VDR-interacting protein (DRIP) complex, a complex of 16 proteins that is essential for vitamin D-mediated gene transcription because of its ability to recruit RNA polymerase II to vitamin D-responsive genes (14, 25). Oda et al (24) used this approach to identify the proteins associated with the VDR in nuclear extracts from proliferating and differentiated keratinocytes. Although they identified a similar complex of proteins interacting with the VDR in proliferating keratinocytes (eg, DRIP complex members and RXR), they found that key members of the complex had changed in differentiated keratinocytes. Proteomic analysis of the proteins associated with the VDR in proliferating and differentiated nuclear extracts showed that at least 5 members of the DRIP complex were lost with differentiation but 2 new proteins, SRC-2 and SRC-3, became prominent members of the complex. That study suggested that the complex mediating VDR-mediated gene expression might not be uniform across vitamin D target tissues or even within the cells of a tissue at different stages of their life spans. This could account for the observed diversity of sensitivity of various cell types/tissues to 1,25(OH)2D3 treatment or the ability of a vitamin D analog to work in one tissue but not another.

At least one other line of vitamin D research might also be improved with a proteomics approach. Specifically, it is now clear that vitamin D stimulates rapid activation of signal transduction pathways, eg, it activates several protein kinases, leading to phosphorylation of various proteins (15). The final targets of the kinases activated by 1,25(OH)2D3 have not yet been identified but it is clear that these targets of phosphorylation are critical for understanding the biological importance of these rapid vitamin D actions. It is apparent that new technologies developed for assessment of the phosphoproteome (26) are likely to identify not only the signal transduction pathways used by 1,25(OH)2D3 but also the terminal proteins whose biological functions are either activated or inhibited through stimulation of these phosphorylation cascades.

Identification of functional biomarkers of vitamin D actions
A major issue in the area of vitamin D research is defining the amount of vitamin D needed for optimal health. Although there are some concerns about the reliability and cross-comparison of various assays for measuring 25(OH)D3 concentrations as an index of vitamin D status (27), the main issue is how 25(OH)D3 concentrations relate to function, ie, bone biological processes related to the risk of osteoporosis or epithelial cell biological processes related to cancer risk. The data from several studies suggest that the cutoff value for adequate serum 25(OH)D3 concentrations may be higher than simply the concentration that prevents rickets (2–4). This idea is based on the association of high vitamin D concentrations with suppression of serum PTH concentrations, a measurement that is being used as a surrogate marker of bone resorption (with the assumption that maximal suppression of PTH is necessary for maximal protection of bone). However, is the PTH concentration an appropriate marker for bone resorption? Or would linking serum 25(OH)D3 concentrations to a more relevant functional endpoint or using a biomarker that is correlated directly with a functional endpoint be better for defining optimal vitamin D status? And what about cancer risk? Should we assume that the serum 25(OH)D3 concentration that is optimal for the protection of bone is optimal for the prevention of cancer? In a perfect world, scientists would be able to directly relate vitamin D status to factors such as fracture incidence, changes in bone density with time, or the development of prostate, breast, or colon cancer. However, such studies would require very large populations and a long study period (as well as being ethically questionable). An alternative might be to correlate vitamin D status with an intermediate endpoint; for bone health, for example, active bone resorption could be measured for a smaller, more controlled, study population. Unfortunately, the means to study active bone resorption (and early functional indices of cancer risk) are currently inadequate for this task. Bone density changes require long periods for reliable observation, and current serum markers are highly variable and thus less reliable. This indicates that new biomarkers of relevant functional endpoints important for bone health and vitamin D biological processes are needed.

Cancer researchers are leading the way toward using the techniques of genomics and proteomics to provide diagnostic markers. For example, Sorlie et al (28) used gene expression profiling of breast tissue biopsies to define the discriminating diagnostic signatures for 6 distinct classes of breast cancer, which suggests that a molecular signature could be used instead of a more subjective histologic evaluation. Although that study and other cancer-profiling studies made great use of gene expression profiles for tumor biopsies as a classification system, the need for a tissue biopsy is a serious limitation to the use of gene expression profiling as a general screening tool for healthy people. Effective biomarkers for assessing healthy people are likely to come from readily available, minimally invasive sampling of blood, blood cells, serum, or urine. The use of serum proteomic profiles, such as that used for ovarian cancer diagnosis, as reported by Petricoin et al (29), may be a more fruitful approach for issues related to nutrient status and health.

In the "omic" approach to identifying and defining biomarkers of disease or of physiologic deficits, the basic idea is to compare the profiles of individuals in 2 well-defined groups (eg, vitamin D-replete and vitamin D-deficient subjects) and then define the changes in the profiles that are correlated best with changes in the condition of interest. For continuous variables such as vitamin D status, the markers may be continuous or they may exhibit breakpoints. Although this approach may provide us with assessment parameters that may be easier to measure than vitamin D metabolite concentrations, without functional correlates this would be no more informative than serum 25(OH)D3 concentrations. Functional assessment could be included to strengthen the relationship; this general scheme is presented in Figure 4. For example, bone resorption could be measured directly through assessment of the release of calcium from bone. Ongoing work in Dr. Connie Weaver's laboratory at Purdue University suggests that a new technique for labeling bones in vivo with small amounts of the natural isotope 41Ca may provide an objective accurate measure of bone resorption (30). 41Ca is a long-lived isotope of calcium that can be produced inexpensively through neutron activation. Because miniscule amounts of 41Ca can be measured with atomic mass spectroscopy, radiologically benign amounts of the isotope can be administered to human subjects to label their bones. After 41Ca has been cleared from soft tissues (100 d), the appearance of 41Ca in the serum or urine is a direct reflection of calcium lost from bone. Because of the long half-life of 41Ca and the relatively low turnover of bone, subjects can be examined repeatedly for >10 y (permitting assessment of multiple treatments or multiple levels of a treatment for the same subject). With the exception of the long study period, this is consistent with work that was previously conducted with animals and high concentrations of 45Ca or tetracycline (31).


View larger version (22K):
FIGURE 4.. Scheme for the identification and cross-validation of serum biomarkers related to vitamin D (VD) status and bone resorption.

 
A discordance in the relationship between the effects of vitamin D status on bone resorption and on the serum proteome could have important physiologic implications. For example, if serum proteins changed as vitamin D status was reduced but this was not accompanied by coordinate changes in bone resorption (41Ca release), then this could be an indication that the vitamin D biomarker is related more to cancer risk (or one of the other potential functions now being proposed for vitamin D, such as diabetes mellitus risk). This could prove to be the basis for future work on the vitamin D–non-bone disease connections.

CONCLUSIONS

The field of vitamin D biology has a long history. However, we are now being confronted with both basic questions related to vitamin D actions and applied issues regarding how to use our fundamental understanding of vitamin D biological processes. Advances in chromatographic techniques have permitted us to better understand vitamin D metabolism, and the revolution in molecular biology has helped us explain the mechanistic basis for vitamin D actions on cells. In this review, I have attempted to explain how the new approaches of genomics and proteomics have the potential to expand our understanding of vitamin D biological processes and the role of vitamin D in human health.

REFERENCES

  1. Holick MF. Vitamin D: a millennium perspective. J Cell Biochem 2003;88:296–307.
  2. Chapuy MC, Preziosi P, Maamer M, et al. Prevalence of vitamin D insufficiency in an adult normal population. Osteoporos Int 1997;7:439–43.
  3. Thomas MK, Lloyd-Jones DM, Thadhani RI, et al. Hypovitaminosis D in medical inpatients. N Engl J Med 1998;338:777–83.
  4. Vieth R, Ladak Y, Walfish PG. Age-related changes in the 25-hydroxyvitamin D versus parathyroid hormone relationship suggest a different reason why older adults require more vitamin D. J Clin Endocrinol Metab 2003;88:185–91.
  5. Barger-Lux MJ, Heaney RP, Lanspa SJ, Healy JC, DeLuca HF. An investigation of sources of variation in calcium absorption efficiency. J Clin Endocrinol Metab 1995;80:406–11.
  6. Heaney RP, Barger-Lux MJ, Dowell MS, Chen TC, Holick MF. Calcium absorptive effects of vitamin D and its major metabolites. J Clin Endocrinol Metab 1997;82:4111–6.
  7. Devine A, Wilson SG, Dick IM, Prince RL. Effects of vitamin D metabolites on intestinal calcium absorption and bone turnover in elderly women. Am J Clin Nutr 2002;75:283–8.
  8. Heaney RP, Dowell MS, Hale CA, Bendich A. Calcium absorption varies within the reference range for serum 25-hydroxyvitamin D. J Am Coll Nutr 2003;22:142–6.
  9. Trivedi DP, Doll R, Khaw KT. Effect of four monthly oral vitamin D3 (cholecalciferol) supplementation on fractures and mortality in men and women living in the community: randomised double blind controlled trial. BMJ 2003;326:469–74.
  10. Zehnder D, Bland R, Williams MC, et al. Extrarenal expression of 25-hydroxyvitamin D3-1-hydroxylase. J Clin Endocrinol Metab 2001;86:888–94.
  11. Whitlatch LW, Young MV, Schwartz GG, et al. 25-Hydroxyvitamin D-1-hydroxylase activity is diminished in human prostate cancer cells and is enhanced by gene transfer. J Steroid Biochem Mol Biol 2002;81:135–40.
  12. Hsu JY, Feldman D, McNeal JE, Peehl DM. Reduced 1-hydroxylase activity in human prostate cancer cells correlates with decreased susceptibility to 25-hydroxyvitamin D3-induced growth inhibition. Cancer Res 2001;61:2852–6.
  13. Haussler MR, Whitfield GK, Haussler CA, et al. The nuclear vitamin D receptor: biological and molecular regulatory properties revealed. J Bone Miner Res 1998;13:325–49.
  14. Rachez C, Freedman LP. Mediator complexes and transcription. Curr Opin Cell Biol 2001;13:274–80.
  15. Farach-Carson MC, Nemere I. Membrane receptors for vitamin D steroid hormones: potential new drug targets. Curr Drug Targets 2003;4:67–76.
  16. Zhu H, Bilgin M, Snyder M. Proteomics. Annu Rev Biochem 2003;72:783–812.
  17. Zhu HJ, Burgess AW. Regulation of transforming growth factor-ß signaling. Mol Cell Biol Res Commun 2001;4:321–30.
  18. Feder JN, Gnirke A, Thomas W, et al. A novel MHC class I-like gene is mutated in patients with hereditary haemochromatosis. Nat Genet 1996;13:399–408.
  19. Takeda S, Karsenty G. Central control of bone formation. J Bone Miner Metab 2001;19:195–8.
  20. Little RD, Carulli JP, Del Mastro RG, et al. A mutation in the LDL receptor-related protein 5 gene results in the autosomal dominant high-bone-mass trait. Am J Hum Genet 2002;70:11–9.
  21. Lin R, Nagai Y, Sladek R, et al. Expression profiling in squamous carcinoma cells reveals pleiotropic effects of vitamin D3 analog EB1089 signaling on cell proliferation, differentiation, and immune system regulation. Mol Endocrinol 2002;16:1243–56.
  22. Toell A, Polly P, Carlberg C. All natural DR3-type vitamin D response elements show a similar functionality in vitro. Biochem J 2000;352:301–9.
  23. Li X, Zheng W, Li YC. Altered gene expression profile in the kidney of vitamin D receptor knockout mice. J Cell Biochem 2003;89:709–19.
  24. Oda Y, Sihlbom C, Chalkley RJ, et al. Two distinct coactivators, DRIP/Mediator and SRC/p160, are differentially involved in VDR transactivation during keratinocyte differentiation. Mol Endocrinol 2003;17:2329–39.
  25. Rachez C, Lemon BD, Suldan Z, et al. Ligand-dependent transcription activation by nuclear receptors requires the DRIP complex. Nature 1999;398:824–8.
  26. Kalume DE, Molina H, Pandey A. Tackling the phosphoproteome: tools and strategies. Curr Opin Chem Biol 2003;7:64–9.
  27. Lips P, Chapuy MC, Dawson-Hughes B, Pols HA, Holick MF. An international comparison of serum 25-hydroxyvitamin D measurements. Osteoporos Int 1999;9:394–7.
  28. Sorlie T, Perou CM, Tibshirani R, et al. Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications. Proc Natl Acad Sci USA 2001;98:10869–74.
  29. Petricoin EF, Ardekani AM, Hitt BA, et al. Use of proteomic patterns in serum to identify ovarian cancer. Lancet 2002;359:572–7.
  30. Jackson GS, Weaver C, Elmore D. Use of accelerator mass spectrometry for studies in nutrition. Nutr Res Rev 2001;14:317–34.
  31. Muhlbauer RC, Fleisch H. A method for continual monitoring of bone resorption in rats: evidence for a diurnal rhythm. Am J Physiol 1990;259:R679–89.

作者: James C Fleet
医学百科App—中西医基础知识学习工具
  • 相关内容
  • 近期更新
  • 热文榜
  • 医学百科App—健康测试工具