1
|
Allen LH, Fenech M, LeVatte MA, West KP, Wishart DS. Multiomics: Functional Molecular Biomarkers of Micronutrients for Public Health Application. Annu Rev Nutr 2024; 44:125-153. [PMID: 39207879 DOI: 10.1146/annurev-nutr-062322-022751] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/04/2024]
Abstract
Adequate micronutrient intake and status are global public health goals. Vitamin and mineral deficiencies are widespread and known to impair health and survival across the life stages. However, knowledge of molecular effects, metabolic pathways, biological responses to variation in micronutrient nutriture, and abilities to assess populations for micronutrient deficiencies and their pathology remain lacking. Rapidly evolving methodological capabilities in genomics, epigenomics, proteomics, and metabolomics offer unparalleled opportunities for the nutrition research community to link micronutrient exposure to cellular health; discover new, arguably essential micronutrients of microbial origin; and integrate methods of molecular biology, epidemiology, and intervention trials to develop novel approaches to assess and prevent micronutrient deficiencies in populations. In this review article, we offer new terminology to specify nutritional application of multiomic approaches and encourage collaboration across the basic to public health sciences to advance micronutrient deficiency prevention.
Collapse
Affiliation(s)
- Lindsay H Allen
- Western Human Nutrition Research Center, United States Department of Agriculture, Agricultural Research Service, Davis, California, USA
- Department of Nutrition, University of California, Davis, California, USA
| | - Michael Fenech
- Clinical and Health Sciences, University of South Australia, Adelaide, South Australia, Australia
- Genome Health Foundation, North Brighton, South Australia, Australia
| | - Marcia A LeVatte
- Department of Biological Sciences, University of Alberta, Edmonton, Alberta, Canada
| | - Keith P West
- Center for Human Nutrition, Department of International Health, Johns Hopkins Bloomberg School of Public Health, Baltimore, Maryland, USA;
| | - David S Wishart
- Department of Biological Sciences, University of Alberta, Edmonton, Alberta, Canada
- Department of Computing Science, University of Alberta, Edmonton, Alberta, Canada
- Department of Laboratory Medicine and Pathology, University of Alberta, Edmonton, Alberta, Canada
- Faculty of Pharmacy and Pharmaceutical Sciences, University of Alberta, Edmonton, Alberta, Canada
| |
Collapse
|
2
|
Biba DA, Wolf YI, Koonin EV, Rochman ND. Balance between asymmetric allocation and repair of somatic damage in unicellular life forms as an ancient form of r/K selection. Proc Natl Acad Sci U S A 2024; 121:e2400008121. [PMID: 38787879 PMCID: PMC11145259 DOI: 10.1073/pnas.2400008121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/01/2024] [Accepted: 04/24/2024] [Indexed: 05/26/2024] Open
Abstract
Over the course of multiple divisions, cells accumulate diverse nongenetic, somatic damage including misfolded and aggregated proteins and cell wall defects. If the rate of damage accumulation exceeds the rate of dilution through cell growth, a dedicated mitigation strategy is required to prevent eventual population collapse. Strategies for somatic damage control can be divided into two categories, asymmetric allocation and repair, which are not, in principle, mutually exclusive. We explore a mathematical model to identify the optimal strategy, maximizing the total cell number, over a wide range of environmental and physiological conditions. The optimal strategy is primarily determined by extrinsic, damage-independent mortality and the physiological model for damage accumulation that can be either independent (linear) or increasing (exponential) with respect to the prior accumulated damage. Under the linear regime, the optimal strategy is either exclusively repair or asymmetric allocation, whereas under the exponential regime, the optimal strategy is a combination of asymmetry and repair. Repair is preferred when extrinsic mortality is low, whereas at high extrinsic mortality, asymmetric damage allocation becomes the strategy of choice. We hypothesize that at an early stage of life evolution, optimization over repair and asymmetric allocation of somatic damage gave rise to r and K selection strategists.
Collapse
Affiliation(s)
- Dmitry A. Biba
- National Center for Biotechnology Information, National Library of Medicine, NIH, Bethesda, MD20894
- Oak Ridge Institute for Science and Education, Oak Ridge, TN37830
| | - Yuri I. Wolf
- National Center for Biotechnology Information, National Library of Medicine, NIH, Bethesda, MD20894
| | - Eugene V. Koonin
- National Center for Biotechnology Information, National Library of Medicine, NIH, Bethesda, MD20894
| | - Nash D. Rochman
- National Center for Biotechnology Information, National Library of Medicine, NIH, Bethesda, MD20894
- Institute for Implementation Science in Population Health, City University of New York, New York, NY10027
- Department of Epidemiology and Biostatistics, Graduate School of Public Health and Health Policy City, University of New York, New York, NY10027
| |
Collapse
|
3
|
Biba DA, Wolf YI, Koonin EV, Rochman ND. Unicellular life balances asymmetric allocation and repair of somatic damage representing the origin of r/K selection. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.11.21.568103. [PMID: 38076808 PMCID: PMC10705550 DOI: 10.1101/2023.11.21.568103] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/22/2023]
Abstract
Over the course of multiple divisions, cells accumulate diverse non-genetic, somatic damage including misfolded and aggregated proteins and cell wall defects. If the rate of damage accumulation exceeds the rate of dilution through cell growth, a dedicated mitigation strategy is required to prevent eventual population collapse. Strategies for somatic damage control can be divided into two categories, asymmetric allocation and repair, which are not, in principle, mutually exclusive. Through mathematical modelling, we identify the optimal strategy, maximizing the total cell number, over a wide range of environmental and physiological conditions. The optimal strategy is primarily determined by extrinsic (damage-independent) mortality and the physiological model for damage accumulation that can be either independent (linear) or increasing (exponential) with respect to the prior accumulated damage. Under the linear regime, the optimal strategy is either exclusively repair or asymmetric allocation whereas under the exponential regime, the optimal strategy is mixed. Repair is preferred when extrinsic mortality is low, whereas at high extrinsic mortality, asymmetric damage allocation becomes the strategy of choice. We hypothesize that optimization over somatic damage repair and asymmetric allocation in early cellular life forms gave rise to the r and K selection strategies.
Collapse
Affiliation(s)
- Dmitry A. Biba
- National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, MD, USA
- Oak Ridge Institute for Science and Education, Oak Ridge, TN, USA
| | - Yuri I. Wolf
- National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, MD, USA
| | - Eugene V. Koonin
- National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, MD, USA
| | - Nash D. Rochman
- National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, MD, USA
- Institute for Implementation Science in Population Health (ISPH), City University of New York (CUNY), New York, NY, USA
- Department of Epidemiology and Biostatistics, Graduate School of Public Health and Health Policy, City University of New York (CUNY), New York, NY, USA
| |
Collapse
|
4
|
Bai Y, Xie M, Zhu Y. Mechanism underlying Müller cell pyroptosis and its role in the development of proliferative vitreoretinopathy. Clinics (Sao Paulo) 2023; 78:100241. [PMID: 37418795 DOI: 10.1016/j.clinsp.2023.100241] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 11/04/2022] [Revised: 06/08/2023] [Accepted: 06/19/2023] [Indexed: 07/09/2023] Open
Abstract
OBJECTIVES To explore the mechanism underlying Müller Cell Pyroptosis (MCP) and its role in the development of Proliferative Vitreoretinopathy (PVR). METHOD The expression of pyroptosis-related factors, namely, cysteinyl aspartate-specific proteinase (caspase-1), interleukin (IL)-1β, IL-18, and Gasdermin D (GSDMD), was detected by quantitative Reverse Transcription Polymerase Chain Reaction (qRT-PCR) and western blotting at the mRNA and protein levels, respectively, in retinal tissues. Müller and spontaneously Arising Retinal Pigment Epithelia (ARPE)-19 primary cells with GSDMD overexpression or knockdown were cultivated. Western blotting was used to detect the levels of the following pyroptosis-related factors in retinal tissues: caspase-1, IL-1β, IL-18, and GSDMD. Through Cell Adhesion (CA) experiments, the changes in ARPE-19 CA in each group were observed. The migration and invasion of ARPE-19 cells were measured using the Transwell assay. The proliferation of ARPE-19 cells was measured with a Cell Counting Kit 8 (CCK-8) assay. Finally, the expression of the cytokines IL-1β and IL-18 in the ARPE-19 cell culture medium was detected using the Enzyme-Linked Immunosorbent Assay (ELISA). RESULTS Compared with the surrounding normal tissues, the expression of caspase-1, IL-1β, IL-18, and GSDMD at the protein and mRNA levels in the retinal proliferative membrane samples of the patients decreased significantly (p < 0.05). MCP significantly enhanced ARPE-19 CA, migration and invasion, proliferation, and cytokine expression (p < 0.05). CONCLUSIONS MCP can promote the development of PVR lesions.
Collapse
Affiliation(s)
- Yue Bai
- Department of Ophthalmology, The First Affiliated Hospital of Fujian Medical University, Fuzhou, China
| | - Maosong Xie
- Department of Ophthalmology, The First Affiliated Hospital of Fujian Medical University, Fuzhou, China
| | - Yihua Zhu
- Department of Ophthalmology, The First Affiliated Hospital of Fujian Medical University, Fuzhou, China.
| |
Collapse
|
5
|
Duan Y, Yan J, Zhu Y, Zhang C, Tao X, Ji H, Zhang M, Wang X, Wang L. Limited accumulation of high-frequency somatic mutations in a 1700-year-old Osmanthus fragrans tree. TREE PHYSIOLOGY 2022; 42:2040-2049. [PMID: 35640149 DOI: 10.1093/treephys/tpac058] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/10/2022] [Accepted: 05/15/2022] [Indexed: 06/15/2023]
Abstract
Lifespan varies greatly between and within species. Mutation accumulation is considered an important factor explaining this life-history trait. However, direct assessment of somatic mutations in long-lived species is still rare. In this study, we sequenced a 1700-year-old sweet olive tree and analysed the high-frequency somatic mutations accumulated in its six primary branches. We found the lowest per-year mutation accumulation rate in this oldest tree among those studied via the whole-genome sequencing approach. Investigation of mutation profiles suggests that this low rate of high-frequency mutation was unlikely to result from strong purifying selection. More intriguingly, on a per-branching scale, the high-frequency mutation accumulation rate was similar among the long-lived individuals such as oak, wild peach and sweet olive investigated here. We therefore suggest the possibility that the accumulation of high-frequency somatic mutations in very long-lived trees might have an upper boundary due to both the possible limited number of stem cell divisions and the early segregation of the stem cell lineage.
Collapse
Affiliation(s)
- Yifan Duan
- Co-Innovation Center for Sustainable Forestry in Southern China, Key Laboratory of State Forestry and Grassland Administration on Subtropical Forest Biodiversity Conservation, College of Biology and the Environment, Nanjing Forestry University, 159 Longpan Road, Xuanwu District, Nanjing 210037, China
- International Cultivar Registration Center for Osmanthus, Nanjing Forestry University, 159 Longpan Road, Xuanwu District, Nanjing 210037, China
| | - Jiping Yan
- Co-Innovation Center for Sustainable Forestry in Southern China, Key Laboratory of State Forestry and Grassland Administration on Subtropical Forest Biodiversity Conservation, College of Biology and the Environment, Nanjing Forestry University, 159 Longpan Road, Xuanwu District, Nanjing 210037, China
- International Cultivar Registration Center for Osmanthus, Nanjing Forestry University, 159 Longpan Road, Xuanwu District, Nanjing 210037, China
| | - Yue Zhu
- Co-Innovation Center for Sustainable Forestry in Southern China, Key Laboratory of State Forestry and Grassland Administration on Subtropical Forest Biodiversity Conservation, College of Biology and the Environment, Nanjing Forestry University, 159 Longpan Road, Xuanwu District, Nanjing 210037, China
- International Cultivar Registration Center for Osmanthus, Nanjing Forestry University, 159 Longpan Road, Xuanwu District, Nanjing 210037, China
| | - Cheng Zhang
- Co-Innovation Center for Sustainable Forestry in Southern China, Key Laboratory of State Forestry and Grassland Administration on Subtropical Forest Biodiversity Conservation, College of Biology and the Environment, Nanjing Forestry University, 159 Longpan Road, Xuanwu District, Nanjing 210037, China
- International Cultivar Registration Center for Osmanthus, Nanjing Forestry University, 159 Longpan Road, Xuanwu District, Nanjing 210037, China
| | - Xiuhua Tao
- Vegetable and Flowers Research Institute, Jiangxi Academy of Agricultural Sciences, 1738 Liantang Middle Blvd, Nanchang 330200, China
| | - Hongli Ji
- Vegetable and Flowers Research Institute, Jiangxi Academy of Agricultural Sciences, 1738 Liantang Middle Blvd, Nanchang 330200, China
| | - Min Zhang
- Co-Innovation Center for Sustainable Forestry in Southern China, Key Laboratory of State Forestry and Grassland Administration on Subtropical Forest Biodiversity Conservation, College of Biology and the Environment, Nanjing Forestry University, 159 Longpan Road, Xuanwu District, Nanjing 210037, China
- International Cultivar Registration Center for Osmanthus, Nanjing Forestry University, 159 Longpan Road, Xuanwu District, Nanjing 210037, China
| | - Xianrong Wang
- Co-Innovation Center for Sustainable Forestry in Southern China, Key Laboratory of State Forestry and Grassland Administration on Subtropical Forest Biodiversity Conservation, College of Biology and the Environment, Nanjing Forestry University, 159 Longpan Road, Xuanwu District, Nanjing 210037, China
- International Cultivar Registration Center for Osmanthus, Nanjing Forestry University, 159 Longpan Road, Xuanwu District, Nanjing 210037, China
| | - Long Wang
- State Key Laboratory of Pharmaceutical Biotechnology, School of Life Sciences, Nanjing University, 163 Xianlin Avenue, Qixia District. Nanjing 210023, China
| |
Collapse
|
6
|
Bräutigam C, Smerlak M. Diffusion approximations in population genetics and the rate of Muller's ratchet. J Theor Biol 2022; 550:111236. [PMID: 35926567 DOI: 10.1016/j.jtbi.2022.111236] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2021] [Revised: 07/13/2022] [Accepted: 07/25/2022] [Indexed: 10/16/2022]
Abstract
The Wright-Fisher binomial model of allele frequency change is often approximated by a scaling limit in which selection, mutation and drift all decrease at the same 1/N rate. This construction restricts the applicability of the resulting 'Wright-Fisher diffusion equation' to the weak selection, weak mutation regime of evolution. We argue that diffusion approximations of the Wright-Fisher model can be used more generally, for instance in cases where genetic drift is much weaker than selection. One important example of this regime is Muller's ratchet phenomenon, whereby deleterious mutations slowly but irreversibly accumulate through rare stochastic fluctuations. Using a modified diffusion equation we derive improved analytical estimates for the mean click time of the ratchet.
Collapse
Affiliation(s)
- Camila Bräutigam
- Max Planck Institute for Mathematics in the Sciences, Leipzig, Germany
| | - Matteo Smerlak
- Max Planck Institute for Mathematics in the Sciences, Leipzig, Germany.
| |
Collapse
|
7
|
Adamek K, Jones AMP, Torkamaneh D. Accumulation of somatic mutations leads to genetic mosaicism in cannabis. THE PLANT GENOME 2022; 15:e20169. [PMID: 34806848 DOI: 10.1002/tpg2.20169] [Citation(s) in RCA: 16] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/28/2021] [Accepted: 09/22/2021] [Indexed: 05/22/2023]
Abstract
Cannabis (Cannabis sativa L.) is typically propagated using stem cuttings taken from mother plants to produce genetically uniform propagules. However, producers anecdotally report that clonal lines deteriorate over time and eventually produce clones with less vigor and lower cannabinoid levels than the original mother plant. While the cause of this deterioration has not been investigated, one potential contributor is the accumulation of somatic mutations within the plant. To test this, we used deep sequencing of whole genomes (>50×) to compare the variability within an individual cannabis cultivar Honey Banana plant sampled at the bottom, middle, and top. We called over six million sequence variants based on a reference genome and found that the top had the most by a sizable amount. Comparing the variants among the samples uncovered that nearly 600,000 (34%) were unique to the top while the bottom only contained 148,000 (12%), and middle with 77,000 (9%) unique variants. Bioinformatics tools were used to identify mutations in critical cannabinoid-terpene biosynthesis pathways. While none were identified as high impact, four genes contained more than double the average level of nucleotide diversity (π) in or near the gene. Two genes code for essential enzymes required for the cannabinoid pathway while the other two are in the terpene pathways, demonstrating that mutations were accumulating within these pathways and could influence their function. Overall, a measurable number of intraplant genetic diversity was discovered that could impact long-term genetic fidelity of clonal lines and potentially contribute to the observed decline in vigor and cannabinoid content.
Collapse
Affiliation(s)
- Kristian Adamek
- Dep. of Plant Agriculture, Univ. of Guelph, Guelph, ON, N1G 2W1, Canada
| | | | | |
Collapse
|
8
|
Singh LN, Kao SH, Wallace DC. Unlocking the Complexity of Mitochondrial DNA: A Key to Understanding Neurodegenerative Disease Caused by Injury. Cells 2021; 10:cells10123460. [PMID: 34943968 PMCID: PMC8715673 DOI: 10.3390/cells10123460] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2021] [Revised: 12/03/2021] [Accepted: 12/06/2021] [Indexed: 12/12/2022] Open
Abstract
Neurodegenerative disorders that are triggered by injury typically have variable and unpredictable outcomes due to the complex and multifactorial cascade of events following the injury and during recovery. Hence, several factors beyond the initial injury likely contribute to the disease progression and pathology, and among these are genetic factors. Genetics is a recognized factor in determining the outcome of common neurodegenerative diseases. The role of mitochondrial genetics and function in traditional neurodegenerative diseases, such as Alzheimer’s and Parkinson’s diseases, is well-established. Much less is known about mitochondrial genetics, however, regarding neurodegenerative diseases that result from injuries such as traumatic brain injury and ischaemic stroke. We discuss the potential role of mitochondrial DNA genetics in the progression and outcome of injury-related neurodegenerative diseases. We present a guide for understanding mitochondrial genetic variation, along with the nuances of quantifying mitochondrial DNA variation. Evidence supporting a role for mitochondrial DNA as a risk factor for neurodegenerative disease is also reviewed and examined. Further research into the impact of mitochondrial DNA on neurodegenerative disease resulting from injury will likely offer key insights into the genetic factors that determine the outcome of these diseases together with potential targets for treatment.
Collapse
Affiliation(s)
- Larry N. Singh
- Center for Mitochondrial and Epigenomic Medicine, Division of Human Genetics, The Children’s Hospital of Philadelphia, Philadelphia, PA 19104, USA;
- Correspondence:
| | - Shih-Han Kao
- Resuscitation Science Center, The Children’s Hospital of Philadelphia, Philadelphia, PA 19104, USA;
| | - Douglas C. Wallace
- Center for Mitochondrial and Epigenomic Medicine, Division of Human Genetics, The Children’s Hospital of Philadelphia, Philadelphia, PA 19104, USA;
- Department of Pediatrics, Division of Human Genetics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA 19104, USA
| |
Collapse
|
9
|
Mc Auley MT. DNA methylation in genes associated with the evolution of ageing and disease: A critical review. Ageing Res Rev 2021; 72:101488. [PMID: 34662746 DOI: 10.1016/j.arr.2021.101488] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2021] [Revised: 09/30/2021] [Accepted: 10/12/2021] [Indexed: 12/28/2022]
Abstract
Ageing is characterised by a physical decline in biological functioning which results in a progressive risk of mortality with time. As a biological phenomenon, it is underpinned by the dysregulation of a myriad of complex processes. Recently, however, ever-increasing evidence has associated epigenetic mechanisms, such as DNA methylation (DNAm) with age-onset pathologies, including cancer, cardiovascular disease, and Alzheimer's disease. These diseases compromise healthspan. Consequently, there is a medical imperative to understand the link between epigenetic ageing, and healthspan. Evolutionary theory provides a unique way to gain new insights into epigenetic ageing and health. This review will: (1) provide a brief overview of the main evolutionary theories of ageing; (2) discuss recent genetic evidence which has revealed alleles that have pleiotropic effects on fitness at different ages in humans; (3) consider the effects of DNAm on pleiotropic alleles, which are associated with age related disease; (4) discuss how age related DNAm changes resonate with the mutation accumulation, disposable soma and programmed theories of ageing; (5) discuss how DNAm changes associated with caloric restriction intersect with the evolution of ageing; and (6) conclude by discussing how evolutionary theory can be used to inform investigations which quantify age-related DNAm changes which are linked to age onset pathology.
Collapse
Affiliation(s)
- Mark Tomás Mc Auley
- Faculty of Science and Engineering, University of Chester, Exton Park, Chester CH1 4BJ, UK.
| |
Collapse
|
10
|
Abstract
PURPOSE OF REVIEW The aim of this study was to highlight recent evidence on protein requirement and physical activity recommendations in older adults and their interrelationship. RECENT FINDINGS Higher protein intake in older adults is beneficial for better physical function, reduced fracture risk, healthy ageing and lower mortality. However, the association between protein quality and health remains inconsistent. Higher physical activity is associated with a lower risk of physical dependence, frailty, cognitive impairment and mortality. Prospectively increasing or even sustaining physical activity can be protective against adverse health outcomes compared with decreasing physical activity. There is a synergistic protective effect of protein intake and physical activity on health. Data on the bidirectional link between protein intake and physical activity are scarce. Whether increasing physical activity can help to increase energy intake (and thus protein intake) and vice versa require further research. SUMMARY The key role of protein intake and physical activity in optimizing health in older adults is well established. The synergistic effects of these two modifiable lifestyle factors strengthen the need to consider dietary and physical activity recommendations together. Future research may examine the bidirectional link between diet and physical activity, which could potentially be one of the strategies to increase the uptake of recommendations for protein intake and physical activity in older adults.
Collapse
Affiliation(s)
- Suey S Y Yeung
- Department of Medicine and Therapeutics, Faculty of Medicine, The Chinese University of Hong Kong, Hong Kong
| | - Ruth S M Chan
- Department of Applied Biology and Chemical Technology, The Hong Kong Polytechnic University, Hong Kong
| | - Jean Woo
- Department of Medicine and Therapeutics, Faculty of Medicine, The Chinese University of Hong Kong, Hong Kong
- Centre for Nutritional Studies, Faculty of Medicine, The Chinese University of Hong Kong, Hong Kong, China
| |
Collapse
|