1
|
Zhong Y, Zhou X, Pan Z, Zhang J, Pan J. Role of epigenetic regulatory mechanisms in age-related bone homeostasis imbalance. FASEB J 2024; 38:e23642. [PMID: 38690719 DOI: 10.1096/fj.202302665r] [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: 12/24/2023] [Revised: 03/05/2024] [Accepted: 04/22/2024] [Indexed: 05/02/2024]
Abstract
Alterations to the human organism that are brought about by aging are comprehensive and detrimental. Of these, an imbalance in bone homeostasis is a major outward manifestation of aging. In older adults, the decreased osteogenic activity of bone marrow mesenchymal stem cells and the inhibition of bone marrow mesenchymal stem cell differentiation lead to decreased bone mass, increased risk of fracture, and impaired bone injury healing. In the past decades, numerous studies have reported the epigenetic alterations that occur during aging, such as decreased core histones, altered DNA methylation patterns, and abnormalities in noncoding RNAs, which ultimately lead to genomic abnormalities and affect the expression of downstream signaling osteoporosis treatment and promoter of fracture healing in older adults. The current review summarizes the impact of epigenetic regulation mechanisms on age-related bone homeostasis imbalance.
Collapse
Affiliation(s)
- Yunyu Zhong
- State Key Laboratory of Oral Diseases, National Center for Stomatology, National Clinical Research Center for Oral Diseases, Department of Oral and Maxillofacial Surgery, West China Hospital of Stomatology, Sichuan University, Chengdu, Sichuan, China
| | - Xueer Zhou
- State Key Laboratory of Oral Diseases, National Center for Stomatology, National Clinical Research Center for Oral Diseases, Department of Oral and Maxillofacial Surgery, West China Hospital of Stomatology, Sichuan University, Chengdu, Sichuan, China
| | - Zijian Pan
- State Key Laboratory of Oral Diseases, National Center for Stomatology, National Clinical Research Center for Oral Diseases, Department of Oral and Maxillofacial Surgery, West China Hospital of Stomatology, Sichuan University, Chengdu, Sichuan, China
| | - Jiankang Zhang
- State Key Laboratory of Oral Diseases, National Center for Stomatology, National Clinical Research Center for Oral Diseases, Department of Oral and Maxillofacial Surgery, West China Hospital of Stomatology, Sichuan University, Chengdu, Sichuan, China
| | - Jian Pan
- State Key Laboratory of Oral Diseases, National Center for Stomatology, National Clinical Research Center for Oral Diseases, Department of Oral and Maxillofacial Surgery, West China Hospital of Stomatology, Sichuan University, Chengdu, Sichuan, China
| |
Collapse
|
2
|
Wang Y, Huang X, Zhang Q, Cheng C, Qin Z, Lu L, Huang Q. The osteoporosis susceptibility SNP rs188303909 at 2q14.2 regulates EN1 expression by modulating DNA methylation and E2F6 binding. J Mol Med (Berl) 2024; 102:273-284. [PMID: 38153509 DOI: 10.1007/s00109-023-02412-6] [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: 07/19/2023] [Revised: 12/11/2023] [Accepted: 12/19/2023] [Indexed: 12/29/2023]
Abstract
EN1 encodes a homeodomain-containing transcription factor and is a determinant of bone density and fracture. Previous powerful genome-wide association studies (GWASs) have identified multiple single-nucleotide polymorphisms (SNPs) near EN1 at 2q14.2 locus for osteoporosis, but the causal SNPs and functional mechanisms underlying these associations are poorly understood. The target genes regulated by the transcription factor EN1 are also unclear. In this study, we identified rs188303909, a functional CpG-SNP, as a causal SNP for osteoporosis at 2q14.2 through the integration of functional and epigenomic analyses. Functional experiments demonstrated that unmethylated rs188303909 acted as a strong allele-specific distal enhancer to regulate EN1 expression by modifying the binding of transcription factor E2F6, but rs188303909 methylation attenuated the active effect of E2F6 on EN1 expression. Importantly, transcription factor EN1 could differentially bind osteoporosis GWAS lead SNPs rs4869739-T and rs4355801-G to upregulate CCDC170 and COLEC10 expression, thus promoting bone formation. Our study provided a mechanistic insight into expression regulation of the osteoporosis susceptibility gene EN1, which could be a potential therapeutic target for osteoporosis precision medicine. KEY MESSAGES: CpG-SNP rs188303909 is a causal SNP at the osteoporosis susceptibility locus 2q14.2. Rs188303909 distally regulates EN1 expression by modulating DNA methylation and E2F6 binding. EN1 upregulates CCDC170 and COLEC10 expression through osteoporosis GWAS lead SNPs rs4869739 and rs4355801.
Collapse
Affiliation(s)
- Ya Wang
- Hubei Key Laboratory of Genetic Regulation and Integrative Biology, School of Life Sciences, Central China Normal University, Wuhan, Hubei, 430079, China
- Key Laboratory of Laboratory Medicine, Ministry of Education, Zhejiang Provincial Key Laboratory of Medical Genetics, School of Laboratory Medicine and Life Sciences, Wenzhou Medical University, Wenzhou, Zhejiang, 325035, China
| | - Xinyao Huang
- Hubei Key Laboratory of Genetic Regulation and Integrative Biology, School of Life Sciences, Central China Normal University, Wuhan, Hubei, 430079, China
| | - Qiongdan Zhang
- Hubei Key Laboratory of Genetic Regulation and Integrative Biology, School of Life Sciences, Central China Normal University, Wuhan, Hubei, 430079, China
| | - Chen Cheng
- Hubei Key Laboratory of Genetic Regulation and Integrative Biology, School of Life Sciences, Central China Normal University, Wuhan, Hubei, 430079, China
| | - Zixuan Qin
- Hubei Key Laboratory of Genetic Regulation and Integrative Biology, School of Life Sciences, Central China Normal University, Wuhan, Hubei, 430079, China
| | - Li Lu
- Hubei Key Laboratory of Genetic Regulation and Integrative Biology, School of Life Sciences, Central China Normal University, Wuhan, Hubei, 430079, China
| | - Qingyang Huang
- Hubei Key Laboratory of Genetic Regulation and Integrative Biology, School of Life Sciences, Central China Normal University, Wuhan, Hubei, 430079, China.
| |
Collapse
|
3
|
Hatmal MM, Al-Hatamleh MAI, Olaimat AN, Alshaer W, Hasan H, Albakri KA, Alkhafaji E, Issa NN, Al-Holy MA, Abderrahman SM, Abdallah AM, Mohamud R. Immunomodulatory Properties of Human Breast Milk: MicroRNA Contents and Potential Epigenetic Effects. Biomedicines 2022; 10:1219. [PMID: 35740242 PMCID: PMC9219990 DOI: 10.3390/biomedicines10061219] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2022] [Revised: 05/15/2022] [Accepted: 05/17/2022] [Indexed: 02/07/2023] Open
Abstract
Infants who are exclusively breastfed in the first six months of age receive adequate nutrients, achieving optimal immune protection and growth. In addition to the known nutritional components of human breast milk (HBM), i.e., water, carbohydrates, fats and proteins, it is also a rich source of microRNAs, which impact epigenetic mechanisms. This comprehensive work presents an up-to-date overview of the immunomodulatory constituents of HBM, highlighting its content of circulating microRNAs. The epigenetic effects of HBM are discussed, especially those regulated by miRNAs. HBM contains more than 1400 microRNAs. The majority of these microRNAs originate from the lactating gland and are based on the remodeling of cells in the gland during breastfeeding. These miRNAs can affect epigenetic patterns by several mechanisms, including DNA methylation, histone modifications and RNA regulation, which could ultimately result in alterations in gene expressions. Therefore, the unique microRNA profile of HBM, including exosomal microRNAs, is implicated in the regulation of the genes responsible for a variety of immunological and physiological functions, such as FTO, INS, IGF1, NRF2, GLUT1 and FOXP3 genes. Hence, studying the HBM miRNA composition is important for improving the nutritional approaches for pregnancy and infant's early life and preventing diseases that could occur in the future. Interestingly, the composition of miRNAs in HBM is affected by multiple factors, including diet, environmental and genetic factors.
Collapse
Affiliation(s)
- Ma’mon M. Hatmal
- Department of Medical Laboratory Sciences, Faculty of Applied Medical Sciences, The Hashemite University, P.O. Box 330127, Zarqa 13133, Jordan;
| | - Mohammad A. I. Al-Hatamleh
- Department of Immunology, School of Medical Sciences, Universiti Sains Malaysia, Kubang Kerian, Kota Bharu 16150, Malaysia;
| | - Amin N. Olaimat
- Department of Clinical Nutrition and Dietetics, Faculty of Applied Medical Sciences, The Hashemite University, P.O. Box 330127, Zarqa 13133, Jordan; (A.N.O.); (M.A.A.-H.)
| | - Walhan Alshaer
- Cell Therapy Center (CTC), The University of Jordan, Amman 11942, Jordan;
| | - Hanan Hasan
- Department of Pathology, Microbiology and Forensic Medicine, School of Medicine, The University of Jordan, Amman 11942, Jordan;
| | - Khaled A. Albakri
- Faculty of Medicine, The Hashemite University, P.O. Box 330127, Zarqa 13133, Jordan;
| | - Enas Alkhafaji
- Department of Pharmaceutical Sciences, Faculty of Pharmacy, The University of Jordan, Amman 11942, Jordan;
| | - Nada N. Issa
- Department of Medical Laboratory Sciences, Faculty of Applied Medical Sciences, The Hashemite University, P.O. Box 330127, Zarqa 13133, Jordan;
| | - Murad A. Al-Holy
- Department of Clinical Nutrition and Dietetics, Faculty of Applied Medical Sciences, The Hashemite University, P.O. Box 330127, Zarqa 13133, Jordan; (A.N.O.); (M.A.A.-H.)
| | - Salim M. Abderrahman
- Department of Biology and Biotechnology, Faculty of Sciences, The Hashemite University, P.O. Box 330127, Zarqa 13133, Jordan;
| | - Atiyeh M. Abdallah
- Department of Biomedical Sciences, College of Health Sciences, QU Health, Qatar University, Doha 2713, Qatar;
| | - Rohimah Mohamud
- Department of Immunology, School of Medical Sciences, Universiti Sains Malaysia, Kubang Kerian, Kota Bharu 16150, Malaysia;
| |
Collapse
|
4
|
Alsheikh AJ, Wollenhaupt S, King EA, Reeb J, Ghosh S, Stolzenburg LR, Tamim S, Lazar J, Davis JW, Jacob HJ. The landscape of GWAS validation; systematic review identifying 309 validated non-coding variants across 130 human diseases. BMC Med Genomics 2022; 15:74. [PMID: 35365203 PMCID: PMC8973751 DOI: 10.1186/s12920-022-01216-w] [Citation(s) in RCA: 23] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2021] [Accepted: 03/17/2022] [Indexed: 02/08/2023] Open
Abstract
Background The remarkable growth of genome-wide association studies (GWAS) has created a critical need to experimentally validate the disease-associated variants, 90% of which involve non-coding variants. Methods To determine how the field is addressing this urgent need, we performed a comprehensive literature review identifying 36,676 articles. These were reduced to 1454 articles through a set of filters using natural language processing and ontology-based text-mining. This was followed by manual curation and cross-referencing against the GWAS catalog, yielding a final set of 286 articles. Results We identified 309 experimentally validated non-coding GWAS variants, regulating 252 genes across 130 human disease traits. These variants covered a variety of regulatory mechanisms. Interestingly, 70% (215/309) acted through cis-regulatory elements, with the remaining through promoters (22%, 70/309) or non-coding RNAs (8%, 24/309). Several validation approaches were utilized in these studies, including gene expression (n = 272), transcription factor binding (n = 175), reporter assays (n = 171), in vivo models (n = 104), genome editing (n = 96) and chromatin interaction (n = 33). Conclusions This review of the literature is the first to systematically evaluate the status and the landscape of experimentation being used to validate non-coding GWAS-identified variants. Our results clearly underscore the multifaceted approach needed for experimental validation, have practical implications on variant prioritization and considerations of target gene nomination. While the field has a long way to go to validate the thousands of GWAS associations, we show that progress is being made and provide exemplars of validation studies covering a wide variety of mechanisms, target genes, and disease areas. Supplementary Information The online version contains supplementary material available at 10.1186/s12920-022-01216-w.
Collapse
Affiliation(s)
- Ammar J Alsheikh
- Genomics Research Center, AbbVie Inc, North Chicago, Illinois, 60064, USA.
| | - Sabrina Wollenhaupt
- Information Research, AbbVie Deutschland GmbH & Co. KG, 67061, Knollstrasse, Ludwigshafen, Germany
| | - Emily A King
- Genomics Research Center, AbbVie Inc, North Chicago, Illinois, 60064, USA
| | - Jonas Reeb
- Information Research, AbbVie Deutschland GmbH & Co. KG, 67061, Knollstrasse, Ludwigshafen, Germany
| | - Sujana Ghosh
- Genomics Research Center, AbbVie Inc, North Chicago, Illinois, 60064, USA
| | | | - Saleh Tamim
- Genomics Research Center, AbbVie Inc, North Chicago, Illinois, 60064, USA
| | - Jozef Lazar
- Genomics Research Center, AbbVie Inc, North Chicago, Illinois, 60064, USA
| | - J Wade Davis
- Genomics Research Center, AbbVie Inc, North Chicago, Illinois, 60064, USA
| | - Howard J Jacob
- Genomics Research Center, AbbVie Inc, North Chicago, Illinois, 60064, USA
| |
Collapse
|
5
|
Chen Z, Li L, Li Z, Wang X, Han M, Gao Z, Wang M, Hu G, Xie X, Du H, Xie Z, Zhang H. Identification of key serum biomarkers for the diagnosis and metastatic prediction of osteosarcoma by analysis of immune cell infiltration. Cancer Cell Int 2022; 22:78. [PMID: 35151325 PMCID: PMC8841093 DOI: 10.1186/s12935-022-02500-6] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2021] [Accepted: 01/30/2022] [Indexed: 12/15/2022] Open
Abstract
Abstract
Background
The role of circular RNAs (circRNAs) and microRNAs (miRNAs) in osteosarcoma (OS) development has not been fully elucidated. Further, the contribution of the immune response to OS progression is not well defined. However, it is known that circRNAs and miRNAs can serve as biomarkers for the diagnosis, prognosis, and therapy of many cancers. Thus, the aim of this study was to identify novel key serum biomarkers for the diagnosis and metastatic prediction of OS by analysis of immune cell infiltration and associated RNA molecules.
Methods
Human OS differentially expressed circRNAs (DEcircRNAs), differentially expressed miRNAs (DEmiRNAs), and differentially expressed mRNAs (DEmRNAs) were identified by analysis of microarray data downloaded from Gene Expression Omnibus (GEO) datasets. Further, characteristic patterns of OS-infiltrating immune cells were analyzed. On this basis, we identified statistically significant transcription factors. Moreover we performed pathway enrichment analysis, constructed protein–protein interaction networks, and devised competitive endogenous RNA (ceRNA) networks. Biological targets of the ceRNA networks were evaluated and potential OS biomarkers confirmed by RT-qPCR analysis of the patients’ serum.
Results
Seven differentially expressed circRNAs, 166 differentially expressed miRNAs, and 175 differentially expressed mRNAs were identified. An evaluation of cellular OS infiltration identified the highest level of infiltration by M0 macrophages, M2 macrophages, and CD8+ T cells, with M0 macrophages and CD8+ T cells as the most prominent. Significant patterns of tumor-infiltrating immune cells were identified by principal component analysis. Moreover, 185 statistically significant transcription factors were associated with OS. Further, in association with immune cell infiltration, hsa-circ-0010220, hsa-miR-326, hsa-miR-338-3p, and FAM98A were identified as potential novel biomarkers for OS diagnosis. Of these, FAM98A had the most promise as a diagnostic marker for OS and OS metastasis. Most importantly, a novel diagnostic model consisting of these four biomarkers (hsa-circ-0010220, hsa-miR-326, hsa-miR-338-3p, and FAM98A) was established with a 0.928 AUC value.
Conclusions
In summary, potential serum biomarkers for OS diagnosis and metastatic prediction were identified based on an analysis of immune cell infiltration. A novel diagnostic model consisting of these four promising serum biomarkers was established. Taken together, the results of this study provide a new perspective by which to understand immunotherapy of OS.
Collapse
|
6
|
Wang Y, Lu L, Niu Y, Zhang Q, Cheng C, Huang H, Huang X, Huang Q. The osteoporosis risk variant rs9820407 at 3p22.1 acts as an allele-specific enhancer to regulate CTNNB1 expression by long-range chromatin loop formation. Bone 2021; 153:116165. [PMID: 34461284 DOI: 10.1016/j.bone.2021.116165] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 06/08/2021] [Revised: 08/25/2021] [Accepted: 08/25/2021] [Indexed: 11/26/2022]
Abstract
Previous powerful genome-wide association studies (GWASs) and whole-genome sequencing have identified multiple single-nucleotide polymorphisms (SNPs) located over 69 kb upstream of CTNNB1 at 3p22.1 locus associated with osteoporosis. The CTNNB1 gene encodes β-catenin that is an integral part of adherens junctions and the primary mediator of the canonical Wnt signaling pathway. The causal variants and underlying molecular mechanisms of the osteoporosis susceptibility locus 3p22.1 remains unknown. Through comprehensive computational analyses, including expression quantitative trait locus (eQTL), high-throughput chromatin interaction (Hi-C), epigenomic and functional annotation, four enhancer SNPs (rs9820407, rs9878224, rs454690 and rs9832204) were prioritized as potential causal SNPs at 3p22.1 for osteoporosis. Rs9820407 displayed the strongest enhancer activity in dual-luciferase assays. Specifically, the minor rs9820407-A can preferentially bind transcription factor FOXC1, elevate the enhancer activity and increase CTNNB1 expression. The architectural protein CTCF was presumably involved in long-range chromatin interaction between rs9820407 and CTNNB1. Our study provided a mechanistic insight into how noncoding enhancer SNP rs9820407 distally regulates CTNNB1 expression and modulates osteoporosis risk.
Collapse
Affiliation(s)
- Ya Wang
- Hubei Key Laboratory of Genetic Regulation and Integrative Biology, School of Life Sciences, Central China Normal University, Wuhan 430079, China
| | - Li Lu
- Hubei Key Laboratory of Genetic Regulation and Integrative Biology, School of Life Sciences, Central China Normal University, Wuhan 430079, China
| | - Yajing Niu
- Hubei Key Laboratory of Genetic Regulation and Integrative Biology, School of Life Sciences, Central China Normal University, Wuhan 430079, China
| | - Qiongdan Zhang
- Hubei Key Laboratory of Genetic Regulation and Integrative Biology, School of Life Sciences, Central China Normal University, Wuhan 430079, China
| | - Chen Cheng
- Hubei Key Laboratory of Genetic Regulation and Integrative Biology, School of Life Sciences, Central China Normal University, Wuhan 430079, China
| | - Han Huang
- Hubei Key Laboratory of Genetic Regulation and Integrative Biology, School of Life Sciences, Central China Normal University, Wuhan 430079, China
| | - Xinyao Huang
- Hubei Key Laboratory of Genetic Regulation and Integrative Biology, School of Life Sciences, Central China Normal University, Wuhan 430079, China
| | - Qingyang Huang
- Hubei Key Laboratory of Genetic Regulation and Integrative Biology, School of Life Sciences, Central China Normal University, Wuhan 430079, China.
| |
Collapse
|
7
|
Ye W, Wang Y, Hou S, Mei B, Liu X, Huang H, Zhou Q, Niu Y, Chen Y, Zhang M, Huang Q. USF3 modulates osteoporosis risk by targeting WNT16, RANKL, RUNX2, and two GWAS lead SNPs rs2908007 and rs4531631. Hum Mutat 2020; 42:37-49. [PMID: 33058301 DOI: 10.1002/humu.24126] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2020] [Revised: 09/09/2020] [Accepted: 10/05/2020] [Indexed: 01/15/2023]
Abstract
Osteoporotic fractures cause major morbidity and mortality in the aging population. Genome-wide association studies (GWAS) have identified USF3 as the novel susceptibility gene of osteoporosis. However, the functional role in bone metabolism and the target gene of the basic helix-loop-helix transcription factor USF3 are unclear. Here, we show that USF3 enhances osteoblast differentiation and suppresses osteoclastogenesis in cultured human osteoblast-like U-2OS cells. Mechanistic studies revealed that transcription factor USF3 antagonistically interacts with anti-osteogenic TWIST1/TCF12 heterodimer in the WNT16 and RUNX2 promoter, and counteracts CREB1 and JUN/FOS in the RANKL promoter. Importantly, the osteoporosis GWAS variant g.1744A>G (rs2908007A>G) located in the WNT16 promoter confers G-allele-specific transcriptional modulation by USF3, TWIST1/TCF12 and TBX5/TBX15, and USF3 transactivates the osteoclastogenesis suppressor WNT16 promoter activity and antagonizes the repression of WNT16 by TWIST1 and TCF12. The risk G allele of osteoporosis GWAS variant g.3260A>G (rs4531631A>G) in the RANKL promoter facilitates the binding of CREB1 and JUN/FOS and enhances transactivation of the osteoclastogenesis contributor RANKL that is inhibited by USF3. Our findings uncovered the functional mechanisms of osteoporosis novel GWAS-associated gene USF3 and lead single nucleotide polymorphisms rs2908007 and rs4531631 in the regulation of bone formation and resorption.
Collapse
Affiliation(s)
- Weiyuan Ye
- Hubei Key Laboratory of Genetic Regulation and Integrative Biology, School of Life Sciences, Central China Normal University, Wuhan, China
| | - Ya Wang
- Hubei Key Laboratory of Genetic Regulation and Integrative Biology, School of Life Sciences, Central China Normal University, Wuhan, China
| | - Sasa Hou
- Hubei Key Laboratory of Genetic Regulation and Integrative Biology, School of Life Sciences, Central China Normal University, Wuhan, China
| | - Bing Mei
- Hubei Key Laboratory of Genetic Regulation and Integrative Biology, School of Life Sciences, Central China Normal University, Wuhan, China
| | - Xinhong Liu
- Hubei Key Laboratory of Genetic Regulation and Integrative Biology, School of Life Sciences, Central China Normal University, Wuhan, China
| | - Han Huang
- Hubei Key Laboratory of Genetic Regulation and Integrative Biology, School of Life Sciences, Central China Normal University, Wuhan, China
| | - Qian Zhou
- Hubei Key Laboratory of Genetic Regulation and Integrative Biology, School of Life Sciences, Central China Normal University, Wuhan, China
| | - Yajing Niu
- Hubei Key Laboratory of Genetic Regulation and Integrative Biology, School of Life Sciences, Central China Normal University, Wuhan, China
| | - Yuanyuan Chen
- Hubei Key Laboratory of Genetic Regulation and Integrative Biology, School of Life Sciences, Central China Normal University, Wuhan, China
| | - Manling Zhang
- Hubei Key Laboratory of Genetic Regulation and Integrative Biology, School of Life Sciences, Central China Normal University, Wuhan, China
| | - Qingyang Huang
- Hubei Key Laboratory of Genetic Regulation and Integrative Biology, School of Life Sciences, Central China Normal University, Wuhan, China
| |
Collapse
|