Research Article | Open Access
Systems Biology Genetic Approach Identifies Serotonin Pathway as a Possible Target for Obstructive Sleep Apnea: Results from a Literature Search Review
Rationale. Overall validity of existing genetic biomarkers in the diagnosis of obstructive sleep apnea (OSA) remains unclear. The objective of this systematic genetic study is to identify “novel” biomarkers for OSA using systems biology approach. Methods. Candidate genes for OSA were extracted from PubMed, MEDLINE, and Embase search engines and DisGeNET database. The gene ontology (GO) analyses and candidate genes prioritization were performed using Enrichr tool. Genes pertaining to the top 10 pathways were extracted and used for Ingenuity Pathway Analysis. Results. In total, we have identified 153 genes. The top 10 pathways associated with OSA include (i) serotonin receptor interaction, (ii) pathways in cancer, (iii) AGE-RAGE signaling in diabetes, (iv) infectious diseases, (v) serotonergic synapse, (vi) inflammatory bowel disease, (vii) HIF-1 signaling pathway, (viii) PI3-AKT signaling pathway, (ix) regulation lipolysis in adipocytes, and (x) rheumatoid arthritis. After removing the overlapping genes, we have identified 23 candidate genes, out of which >30% of the genes were related to the genes involved in the serotonin pathway. Among these 4 serotonin receptors SLC6A4, HTR2C, HTR2A, and HTR1B were strongly associated with OSA. Conclusions. This preliminary report identifies several potential candidate genes associated with OSA and also describes the possible regulatory mechanisms.
According to the Centers for Disease Control and Prevention, sleep disorders now constitute a public health epidemic. The National Healthy Sleep Awareness Project estimated that at least 25 million American adults are affected by OSA for the year 2014 . The underlying pathophysiology and physiological mechanisms of OSA are multifactorial, may vary considerably between individuals, and are not clearly understood . It is estimated that majority of the patients with sleep apnea are undiagnosed  increasing evidence supporting an association between OSA and obesity , hypertension , metabolic syndrome , type 2 diabetes (T2DM)  and cardiovascular mortality , increased risk of psychiatric disorders , and traffic and occupational accidents . Given the potentially serious consequences of OSA, timely diagnosis and recognition, risk stratification, and appropriate treatment are of utmost importance .
In-lab polysomnography (PSG) is the “gold standard” test for the diagnosis of OSA . However, PSG is resource-intensive and cumbersome and requires patients to spend the night under observation in a foreign (hospital) environment. Although home-based sleep monitoring devices overcome these problems, they are less accurate, are prone to technical breakdown, fail to capture other forms of sleep disorders, and underestimate the severity of OSA . Recently, several sleep based psychometric questionnaires are available to screen individuals at low to high OSA risk. However, the lack of high quality evidence from large scale epidemiological studies hinders their applicability in clinical practice . Also, questionnaires did not allow reliable differentiation between different phenotypes of OSA; that is, they were unable to identify individuals who are susceptible to metabolic, cardiovascular, and neurologic effects of OSA. Also, sleep based questionnaires may not be suitable to the general population because they include questions related to sleepiness and not all patients, even those with severe OSA, report sleepiness. For instance, Wisconsin Sleep Cohort Study found that only 37% of patients with severe OSA (Apnea-Hypopnea Index (AHI) score ≥ 30 events/h) are aware of questions related to daytime sleepiness .
Blood biomarkers may offer alternative approach for screening and identifying individuals at risk of developing OSA and related outcomes. Biomarker is defined as a characteristic that is objectively measured and evaluated as an indicator of normal biologic processes, pathogenic processes, or pharmacologic response to a therapeutic intervention. In the case of OSA, the candidate biomarkers should be useful for the following: diagnosis, assessing disease burden and severity, and evaluating response to treatment. Several epidemiological studies suggest that the circulatory levels of natriuretic peptide (BNP) or N-terminal-pro-BNP , leptin, C-reactive protein (CRP) [16, 17], TNFα , vaspin , ghrelin , and interleukin-6  are associated with OSA. To ensure diagnostic accuracy, a biomarker should yield both high sensitivity and specificity to identify individuals at high risk of developing OSA. Such diagnostic performance would limit the necessity of PSG, an expensive and laborious modality, at least in some patients. Additional aspects should include low cost and facility in use, as well as the ability to simultaneously evaluate multiple pathogenic pathways. However, current knowledge of biomarkers for OSA is limited by data collection techniques, disease complexity, and potential confounding factors.
Advances in computational approaches and experimental omics methods that allow the simultaneous analysis of multidimensional data such as DNA, RNA, and proteins in a single analysis have made these systems biology approaches feasible for biomarker discovery . Systems biology is a scientific discipline that endeavors to quantify all of the molecular elements of a biological system to assess their interactions and to integrate that information into graphical network models . Therefore, the objective of this report is to (a) identify “novel” biochemical markers and associated pathways of OSA using systems biology approach, (b) assess the association of the selected candidate genes with OSA comorbid conditions, and (c) study potential interaction of the candidate genes in a biological network context.
2.1. Literature Mining and Identification of Candidate Genes
Literature search was conducted in PubMed, MEDLINE, and Embase search engines for studies published for “genes, polymorphisms, and biomarkers” using the terms “Apnea, Obstructive Sleep” or “Apneas, Obstructive Sleep” or “Obstructive Sleep Apnea” or “Obstructive Sleep Apneas” or “Obstructive Sleep Apnea Syndrome” or “OSAHS” or “Sleep Apnea Hypopnea Syndrome” or “Sleep Apneas, Obstructive” or “Sleep Apnea Syndrome, Obstructive” or “Syndrome, Obstructive Sleep Apnea” or “Syndrome, Sleep Apnea, Obstructive” or “Syndrome, Upper Airway Resistance, Sleep Apnea” or “Upper Airway Resistance Sleep Apnea Syndrome” for obstructive sleep apnea (MeSH® ID: D020181). Furthermore, we queried the DisGeNET database (http://www.disgenet.org/web/DisGeNET/menu)  that integrates information from four repositories (search term: umls:C0520679 for OSA): (a) Online Mendelian Inheritance in Man (OMIM), (b) UniProt/SwissProt (UniProt), (c) Pharmacogenomics Knowledge Base (PHARMGKB), and (d) Comparative Toxicogenomics Database (CTD). DisGeNET can be accessed through the Cytoscape 2.8.3, a platform for complex network analysis. The identified genes were then evaluated to verify authenticity and to remove redundancies. Genes with unknown functions were excluded from the analysis (<5%).
2.2. Gene-Set Enrichment Analysis and Candidate Gene Prioritization
2.3. Network Analysis by Ingenuity Pathway Analysis (IPA)
Finally, the shortlisted genes were submitted to Ingenuity Pathway Analysis (IPA 4.0, Ingenuity Systems Inc., https://www.ingenuity.com/) for mapping to canonical pathways and identifying upstream regulators using a set of criteria: genes and endogenous chemicals, direct and indirect interactions, maximum molecules per network, networks per analysis, humans as the selected species, and all tissues and primary cells. The resulting networks were scored based on the fold change provided by Cuffdiff as log2 (fold change) for each gene. The obtained values correspond to Fisher’s exact test, with the null hypothesis that the molecules within the networks are connected based on chance.
3. Results: Candidate Gene Identification
Literature mining for the search terms as described in the methodology resulted in 2200 abstracts. After carefully screening the abstracts we got 550 genes related to OSA and from that we have shortlisted 153 genes (removal of redundancies, duplicates, unknown functions, etc.), for further analysis. The analysis of statistically overrepresented pathways in the shortlist of genes revealed 28 canonical pathway maps (by ToppGene suite) with confidence level (adjusted for false discovery rate (FDR)). The top 10 pathways and their corresponding genes revealed by Enrichr tool are given in Figures 1(a) and 1(b). This includes (i) serotonin receptor interaction, (ii) pathways in cancer, (iii) AGE-RAGE signaling in diabetes, (iv) infectious diseases, (v) serotonergic synapse, (vi) inflammatory bowel disease, (vii) HIF-1 signaling pathway, (viii) PI3-AKT signaling pathway, (ix) regulation lipolysis in adipocytes, and (x) rheumatoid arthritis. After merging all the pathways, “23” unique genes remained for the core analysis (Box 1). The set of unique genes implicated in OSA showed maximal association with OSA comorbid conditions (Figure 2) such as anxiety disorders (), obesity (), obsessive-compulsive disorder (), dyslipidemia (), social anxiety disorder (), T2DM (), blood pressure (), eating disorders (), depression (), and metabolic syndrome () (Figure 3). Among twenty-three genes identified, 9 (36.0%) genes are involved in serotonin receptor mediated pathways. Therefore, for the subsequent network analysis we focused only on serotonin pathways.
For network analysis Figure 3 shows the Ingenuity Pathway Analysis of the selected candidate genes and their corresponding molecular interactions. The significant network associations include (a) carbohydrate metabolism, small molecule biochemistry, and psychological disorders; (b) behavior, cell signaling, and neurological disease; and (c) drug metabolism, endocrine system development and function, and lipid metabolism. In the predicted merged network generated by IPA algorithm, four serotonin receptors, HTR1B, SLC6A4, HTR2A, and HTR2C, were found to be the potential candidates regulating the OSA. The expression of our candidate genes could be modulated mainly by 3 molecules: 5-hydroxy tryptamine (for HTR1B and SLC6A4), PRL (HTR2A and HTR2C), and HSP90AB1 (for SL6CA4). The important regulators of these genes were insulin (), melatonin (), glucocorticoid (), lithium (), and RAB1A ().
The overall validity of existing biomarkers in the diagnosis of obstructive sleep apnea (OSA) remains unclear. Although the existing adipokines and inflammatory markers seem to have a favorable potential to become a good biomarker to identify OSA they are not specific for the disease. In this study, we have performed the gene mining analyses and found 23 candidate genes that are likely associated with OSA. Functional enrichment analysis suggested that most of these genes are directly involved in the above-mentioned biological processes (Figure 2), which are highly relevant for OSA. Several pathways found in our study have already been discussed by many authors with inconclusive evidence [26–33] before as well as their involvement in the pathogenesis of OSA. For instance, Diefenbach et al.  showed the positive association between EDN1 variant Lys198Asn genotype and OSA. However, plasma levels of EDN1 were not associated with the OSA severity.
The disease pathway analysis showed that the selected candidate genes were also associated with the OSA comorbid conditions such as anxiety disorders, obesity, and dyslipidemia. For the subsequent analysis, we focused only on serotonin pathway because majority of serotonin genes make up this pathway. Dysfunction of the serotoninergic system has long been suspected in sleep disorders and respiratory diseases [35, 36]. 5-HT acts as both a neurotransmitter and a neuromodulator in the human central nervous system (CNS). Serotonin pathway affects food intake, sleep, anxiety, sexual behavior, and mood . Approximately 2% of the body’s serotonin is stored in central nervous system (CNS) and the rest of the body’s serotonin is in the gut and stored peripherally, where it operates as a peripheral hormone, affecting vasoconstriction, intestinal motility, primary hemostasis, liver repair, and the control of the T-cell-mediated immune system . Emerging evidence from animal studies showed that both peripheral and central 5-HT and its receptors were involved in sleep regulation . The oxygen deprivation and hyperventilation which are characteristic of OSA can be mimicked and induced in panic attack which might derange the serotonin pathway mediators. Recently, 5-HT system has been shown to be involved in circadian regulation , glucose and lipid metabolism, and adipocyte differentiation [41, 42]. Serotonergic neurons have a significant influence on sleep/wake cycles due to their multiple connections throughout the cortex, basal forebrain, limbic system, and brainstem areas. The 5-HT systems function predominantly to promote wakefulness and inhibit rapid eye movement (REM) sleep. The recent finding suggests that serotonergic connections have effects on numerous CNS processes, and dysfunction in this system may be implicated in respiratory pathology such as OSA . Serotonin is thought to cause its effects through these membrane-bound receptors . In our analyses, we have identified 4 receptors: SLC6A4, , , and which associated with OSA.
The genomewide association studies demonstrated the association of 5-HTR 2A/2C polymorphisms with OSA. Furthermore, by regulating the magnitude and duration of serotonergic responses, the SLC6A4 (i.e., 5-HTT) is pivotal to the fine-tuning of brain serotonergic neurotransmission and of the peripheral actions of 5-HT .
It is clear that there is an intricate network of interacting genes regulating sleep and the derangement in them leads to OSA. In this preliminary report, we identified 23 novel biomarkers especially from the serotonin pathway which might associate with OSA. Consequently, we demonstrated network-level analysis of potential OSA candidate genes and the plausible regulatory mechanisms. Limited evidence from the clinical trials showed that 5-HT has a critical role influencing respiratory control during sleep. 5-HT acts as a potent central ventilatory stimulant and serves to maintain upper airway patency and maintain eucapnia via the chemoreceptor properties. Relative reductions in 5-HT may lead to the development and worsening of OSA. Further biological experiments are warranted to study the role of peripheral 5-HT and its receptors on the pathogenesis of OSA.
Parts of this study were presented at the 13th Annual World Congress on Insulin Resistance, Diabetes and Cardiovascular Disease (WCIRDC), Los Angeles, California, 19–21 Nov 2015.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Ram Jagannathan contributed to conceptualization of hypothesis, selected genes, analyzed the data, interpreted the results, wrote the manuscript, and reviewed/edited the manuscript. Azizi Seixas interpreted the results, wrote the manuscript, and reviewed/edited the manuscript; David St-Jules, April Rogers, and Lu Hu reviewed/edited the manuscript; Girardin Jean-Louis and Mary Ann Sevick contributed to discussion and reviewed/edited the manuscript.
Mary Ann Sevick is supported by NIH-K24-NR012226 and NIH-R01DK100492.
- P. E. Peppard, T. Young, J. H. Barnet, M. Palta, E. W. Hagen, and K. M. Hla, “Increased prevalence of sleep-disordered breathing in adults,” The American Journal of Epidemiology, vol. 177, no. 9, pp. 1006–1014, 2013.
- D. J. Eckert and A. Malhotra, “Pathophysiology of adult obstructive sleep apnea,” Proceedings of the American Thoracic Society, vol. 5, no. 2, pp. 144–153, 2008.
- T. Young, L. Evans, L. Finn, and M. Palta, “Estimation of the clinically diagnosed proportion of sleep apnea syndrome in middle-aged men and women,” Sleep, vol. 20, no. 9, pp. 705-706, 1997.
- A. R. Schwartz, S. P. Patil, A. M. Laffan, V. Polotsky, H. Schneider, and P. L. Smith, “Obesity and obstructive sleep apnea: pathogenic mechanisms and therapeutic approaches,” Proceedings of the American Thoracic Society, vol. 5, no. 2, pp. 185–192, 2008.
- J. M. Marin, A. Agusti, I. Villar et al., “Association between treated and untreated obstructive sleep apnea and risk of hypertension,” The Journal of the American Medical Association, vol. 307, no. 20, pp. 2169–2176, 2012.
- S. Xu, Y. Wan, M. Xu et al., “The association between obstructive sleep apnea and metabolic syndrome: a systematic review and meta-analysis,” BMC Pulmonary Medicine, vol. 15, no. 1, article 105, 2015.
- X. Wang, Y. Bi, Q. Zhang, and F. Pan, “Obstructive sleep apnoea and the risk of type 2 diabetes: A meta-analysis of prospective cohort studies,” Respirology, vol. 18, no. 1, pp. 140–146, 2013.
- X. Wang, Y. Ouyang, Z. Wang, G. Zhao, L. Liu, and Y. Bi, “Obstructive sleep apnea and risk of cardiovascular disease and all-cause mortality: a meta-analysis of prospective cohort studies,” International Journal of Cardiology, vol. 169, no. 3, pp. 207–214, 2013.
- A. Sharafkhaneh, N. Giray, P. Richardson, T. Young, and M. Hirshkowitz, “Association of psychiatric disorders and sleep apnea in a large cohort,” Sleep, vol. 28, no. 11, pp. 1405–1411, 2005.
- S. Garbarino, O. Guglielmi, A. Sanna, G. L. Mancardi, and N. Magnavita, “Risk of occupational accidents in workers with obstructive sleep apnea: Systematic review and meta-Analysis,” Sleep, vol. 39, no. 6, pp. 1211–1218, 2016.
- L. E. Costa, C. H. G. Uchôa, R. R. Harmon, L. A. Bortolotto, G. Lorenzi-Filho, and L. F. Drager, “Potential underdiagnosis of obstructive sleep apnoea in the cardiology outpatient setting,” Heart, vol. 101, no. 16, pp. 1288–1292, 2015.
- A. Qaseem, P. Dallas, D. K. Owens, M. Starkey, J.-E. C. Holty, and P. Shekelle, “Diagnosis of obstructive sleep apnea in adults: A clinical practice guideline from the American College of Physicians,” Annals of Internal Medicine, vol. 161, no. 3, pp. 210–220, 2014.
- M. Ahmed, N. P. Patel, and I. Rosen, “Portable monitors in the diagnosis of obstructive sleep apnea,” Chest, vol. 132, no. 5, pp. 1672–1677, 2007.
- T. Young, L. Finn, P. E. Peppard et al., “Sleep disordered breathing and mortality: eighteen-year follow-up of the wisconsin sleep cohort,” Sleep, vol. 31, no. 8, pp. 1071–1078, 2008.
- M. Ljunggren, B. Lindahl, J. Theorell-Haglöw, and E. Lindberg, “Association between obstructive sleep apnea and elevated levels of type B natriuretic peptide in a community-based sample of women,” Sleep, vol. 35, no. 11, pp. 1521–1527, 2012.
- A. S. M. Shamsuzzaman, M. Winnicki, P. Lanfranchi et al., “Elevated C-reactive protein in patients with obstructive sleep apnea,” Circulation, vol. 105, no. 21, pp. 2462–2464, 2002.
- T. Yokoe, K. Minoguchi, H. Matsuo et al., “Elevated levels of C-reactive protein and interleukin-6 in patients with obstructive sleep apnea syndrome are decreased by nasal continuous positive airway pressure,” Circulation, vol. 107, no. 8, pp. 1129–1134, 2003.
- D. Gozal, L. D. Serpero, L. Kheirandish-Gozal, O. S. Capdevila, A. Khalyfa, and R. Tauman, “Sleep measures and morning plasma TNF-α levels in children with sleep-disordered breathing,” Sleep, vol. 33, no. 3, pp. 319–325, 2010.
- M. Kiskac, M. Zorlu, M. E. Akkoyunlu et al., “Vaspin and lipocalin-2 levels in severe obsructive sleep apnea,” Journal of Thoracic Disease, vol. 6, no. 6, pp. 720–725, 2014.
- A. Ursavas, Y. O. Ilcol, N. Nalci, M. Karadag, and E. Ege, “Ghrelin, leptin, adiponectin, and resistin levels in sleep apnea syndrome: Role of obesity,” Annals of Thoracic Medicine, vol. 5, no. 3, pp. 161–165, 2010.
- L. Hood, J. R. Heath, M. E. Phelps, and B. Lin, “Systems biology and new technologies enable predictive and preventative medicine,” Science, vol. 306, no. 5696, pp. 640–643, 2004.
- J. Piñero, N. Queralt-Rosinach, À. Bravo et al., “DisGeNET: a discovery platform for the dynamical exploration of human diseases and their genes,” Database, vol. 2015, Article ID bav028, 2015.
- J. Chen, E. E. Bardes, B. J. Aronow, and A. G. Jegga, “ToppGene Suite for gene list enrichment analysis and candidate gene prioritization,” Nucleic Acids Research, vol. 37, no. 2, pp. W305–W311, 2009.
- E. Y. Chen, C. M. Tan, Y. Kou et al., “Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool,” BMC Bioinformatics, vol. 14, article 128, 2013.
- M. V. Kuleshov, M. R. Jones, A. D. Rouillard et al., “Enrichr: a comprehensive gene set enrichment analysis web server 2016 update,” Nucleic Acids Research, vol. 44, no. W1, pp. W90–W97, 2016.
- V. S. Rotenberg, “REM sleep function and brain monoamine regulation: an application of the search activity concept,” Sleep and Sleep Disorders: A Neuropsychopharmacological Approach, pp. 27–35, 2006.
- E. K. Weir and A. Olschewski, “Role of ion channels in acute and chronic responses of the pulmonary vasculature to hypoxia,” Cardiovascular Research, vol. 71, no. 4, pp. 630–641, 2006.
- W. Wu, Z. Li, T. Tang, J. Wu, F. Liu, and L. Gu, “5-HTR2A and IL-6 polymorphisms and obstructive sleep apnea-hypopnea syndrome,” Biomedical Reports, vol. 4, no. 2, pp. 203–208, 2016.
- S. C. Veasey, “Serotonin agonists and antagonists in obstructive sleep apnea: Therapeutic potential,” American Journal of Respiratory Medicine, vol. 2, no. 1, pp. 21–29, 2003.
- I. A. Harsch, P. C. Konturek, C. Koebnick et al., “Leptin and ghrelin levels in patients with obstructive sleep apnoea: effect of CPAP treatment,” European Respiratory Journal, vol. 22, no. 2, pp. 251–257, 2003.
- K. Karkoulias, D. Lykouras, F. Sampsonas et al., “The role of endothelin-1 in obstructive sleep apnea syndrome and pulmonary arterial hypertension: pathogenesis and endothelin-1 antagonists,” Current Medicinal Chemistry, vol. 17, no. 11, pp. 1059–1066, 2010.
- H. Xu, J. Guan, H. Yi et al., “Elevated low-density lipoprotein cholesterol is independently associated with obstructive sleep apnea: evidence from a large-scale cross-sectional study,” Sleep & breathing = Schlaf & Atmung, vol. 20, no. 2, pp. 627–634, 2016.
- A. E. Mirrakhimov and A. M. Ali, “Pathobiology of Obstructive Sleep Apnea-Related Dyslipidemia: Focus on the Liver,” ISRN Cardiology, vol. 2013, Article ID 687069, 5 pages, 2013.
- K. Diefenbach, K. Kretschmer, S. Bauer et al., “Endothelin-1 gene variant Lys198Asn and plasma endothelin level in obstructive sleep apnea,” Cardiology, vol. 112, no. 1, pp. 62–68, 2008.
- S. Sood, J. L. Morrison, H. Liu, and R. L. Horner, “Role of endogenous serotonin in modulating genioglossus muscle activity in awake and sleeping rats,” American Journal of Respiratory and Critical Care Medicine, vol. 172, no. 10, pp. 1338–1347, 2005.
- M. Yilmaz, Y. A. Bayazit, T. U. Ciftci et al., “Association of serotonin transporter gene polymorphism with obstructive sleep apnea syndrome,” Laryngoscope, vol. 115, no. 5, pp. 832–836, 2005.
- I. Lucki, “The spectrum of behaviors influenced by serotonin,” Biological Psychiatry, vol. 44, no. 3, pp. 151–162, 1998.
- H. Watanabe, M. T. Rose, and H. Aso, “Role of peripheral serotonin in glucose and lipid metabolism,” Current Opinion in Lipidology, vol. 22, no. 3, pp. 186–191, 2011.
- D. W. Carleij and M. Rachilovacki, “Role of peripheral serotonin in the regulation of central sleep apneas in rats,” Chest, vol. 115, no. 5, pp. 1397–1401, 1999.
- C. M. Ciarleglio, H. E. S. Resuehr, and D. G. McMahon, “Interactions of the serotonin and circadian systems: Nature and nurture in rhythms and blues,” Neuroscience, vol. 197, pp. 8–16, 2011.
- H. Watanabe, D. Akasaka, H. Ogasawara et al., “Peripheral serotonin enhances lipid metabolism by accelerating bile acid turnover,” Endocrinology, vol. 151, no. 10, pp. 4776–4786, 2010.
- E. Hajduch, F. Rencurel, A. Balendran, I. H. Batty, C. P. Downes, and H. S. Hundal, “Serotonin (5-hydroxytryptamine), a novel regulator of glucose transport in rat skeletal muscle,” Journal of Biological Chemistry, vol. 274, no. 19, pp. 13563–13568, 1999.
- G. B. Richerson, “Serotonergic neurons as carbon dioxide sensors that maintain pH homeostasis,” Nature Reviews Neuroscience, vol. 5, no. 6, pp. 449–461, 2004.
- D. Hoyer, J. P. Hannon, and G. R. Martin, “Molecular, pharmacological and functional diversity of 5-HT receptors,” Pharmacology Biochemistry and Behavior, vol. 71, no. 4, pp. 533–554, 2002.
- K.-P. Lesch, D. Bengel, A. Heils et al., “Association of anxiety-related traits with a polymorphism in the serotonin transporter gene regulatory region,” Science, vol. 274, no. 5292, pp. 1527–1531, 1996.
Copyright © 2017 Ram Jagannathan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.