- 1Department of Neurosurgery, Xuanwu Hospital, Capital Medical University, Beijing, China
- 2Interventional Neuroradiology Diagnosis and Treatment Center, Xuanwu Hospital, Capital Medical University, Beijing, China
- 3School of Life Science, Beijing Institute of Technology, Beijing, China
- 4School of Medical Imaging, North Sichuan Medical College, Nanchong, China
- 5China International Neuroscience Institute (China-INI), Beijing, China
- 6Neuroendovascular Program, Harvard Medical School, Massachusetts General Hospital, Boston, MA, United States
- 7Department of Radiology, Xuanwu Hospital, Capital Medical University, Beijing, China
Background: Basilar artery fenestration has been proposed as a contributor to ischemic stroke, as unique flow patterns induced by fenestration may be related to thrombus formation or insufficiency. This study aimed to evaluate the hemodynamics of basilar artery fenestration (BAF) using computational fluid dynamics (CFD).
Methods: Patients with BAF and normal vertebrobasilar system were recruited and separately evaluated using CFD. Specific geometric vascular models were reconstructed based on 3D-rotational angiography (3D-RA). Patients were divided into the BAF group and control group (i.e., patients with the normal vertebrobasilar system). Hemodynamic and geometric variables were calculated and compared between groups using Student's t-test or Wilcoxon rank-sum test.
Results: Overall, 24 patients were included, with 12 patients each in the BAF group and the control group. The BAF group had a significantly smaller basilar artery diameter than the control group (3.1 ± 0.51 vs. 3.76 ± 0.4, p = 0.002). Compared to the control group, the BAF group had higher values of maxOSI (median, 0.3 vs. 0.09, p = 0.028), TAWSSG (median, 983.42 vs. 565.39, p = 0.038) in the flow confluence, higher SAR-TAWSSG in bifurcation (median, 70.22 vs. 27.65, p = 0.002) and higher SAR-TAWSSG in basilar artery (median, 48.75 vs. 16.17, p < 0.001) of the vertebrobasilar artery.
Conclusions: This pilot study suggested that hemodynamic differences between BAF and normal vertebrobasilar artery across multiple shear flow parameters. The disturbed flow in the BAF may increase the risk of thrombus formation, plaque instability, and subsequent ischemic cerebrovascular events. These should be confirmed by future studies.
Introduction
Artery fenestration is a rare congenital vascular dysplasia caused by incomplete fusion of primitive embryologic vessels, which are divided into two separate channels and converge at the distal end (1). Fenestration is a developmental vascular anomaly that can occur in multiple cerebral arteries, especially in the anterior communicating artery and basilar artery (2). Basilar artery fenestration (BAF) is the second most common form of fenestration with a prevalence of 0.28–5.26% as reported in post-mortem studies (3), 0.3–0.6% in conventional angiography (4), and 1–2.07% in magnetic resonance angiography (MRA) (5–7). The size of the fenestration largely ranges from 1 mm to 5 mm (8), with a slit-like separation or long segment duplication (9). Given the innate small size and relatively low prevalence rate in imaging studies, BAF may be an easily overlooked cerebrovascular anomaly in clinical practice.
Although BAF is rare, its association with aneurysm formation and ischemic stroke has been reported (2, 10–12). Some cases were related to cerebral infarction and diagnosed as an embolic stroke of undetermined source as these patients usually had no vascular risk factors or other diseases. Therefore, this suggests that BAF is a possible cause of cryptogenic stroke (1, 2, 13). The exact mechanism of fenestration-related infarction remains unknown. Fenestration causes local hemodynamic changes, which may heighten the risk of ischemic stroke or infarction.
Computational fluid dynamics (CFD) is used to simulate flow patterns and analyze the relationship between diseases and hemodynamics. Previous studies have explored the pathophysiological relationship between changes in blood flow and the development of atherosclerosis (14, 15). Some researchers speculate that altered blood flow in the presence of artery fenestration promotes ischemic stroke, but the exact mechanism is not well-defined (16).
In this study, we first simulated flow dynamics of the BAF by using CFD and compared hemodynamic parameters with normal vertebral basilar artery. We hypothesized that there would be hemodynamic differences between the two groups. With this study, we hope to provide researchers and clinicians with updated references for future studies.
Materials and Methods
Patient Selection
The medical ethics committee of the hospital approved this study, and all patients provided written informed consent. From January 2020 to March 2021, 24 patients were included in this study, inclusive of 12 patients with BAF (BAF group) and 12 patients with a normal vertebrobasilar system (control group). The inclusion criteria were as follows: (1) patients diagnosed with BAF by digital subtraction angiography; (2) underwent 3D-rotational angiography (3D-RA) and division of the lumen was observed by virtual artery endoscopy; (3) the quality of images was adequate for CFD analysis; (4) fenestrations were located in the proximal segment of the basilar artery. The exclusion criteria were as follows: (1) BAF with pre-existing stenosis or aneurysm; (2) the quality of images was poor for CFD analysis. Patients who underwent MRA and had normal vertebrobasilar system were enrolled and included in the control group.
Patient Groups and Imaging
To analyze the geometry and hemodynamic characteristics of a BAF, patients were divided into two groups, namely, the BAF group and the control group. Baseline data were collected, such as age, sex, smoking history, posterior circulation infarction (PCI), blood pressure level (systolic and diastolic), low-density lipoprotein (LDL), high-density lipoprotein (HDL), fasting glucose, total cholesterol, and triglycerides. This study focused on the intracranial vertebrobasilar bifurcation and basilar artery (Figure 1). In brief, 5-s rotational angiography was performed using the Innova IGS 630 biplane angiography system (General Electric, America). The protocol refers to our previous published paper (17), wherein we performed a total of 24 ml contrast agents at a flow rate of 3 ml/s, 3 s delay time, and 0.2 mm slice thickness. All patients performed virtual arterial endoscopy and confirmed the presence of fenestration. Geometric variables of interest were calculated, including the bifurcation angle, flow confluence diameter (left and right vertebral arteries merging into a basilar artery), and basilar artery diameter, tortuosity, and length. Two experienced neuroradiologists performed the reconstruction and imaging analysis.
Figure 1. (A) A three-dimensional (3D) geometrical model of basilar artery fenestration (BAF) reconstructed from a patient's 3D-rotational angiography (3D-RA). (B–E) The hemodynamic distribution of BAF, including the following parameters: velocity, time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), and time-averaged wall shear stress gradient (TAWSSG). (F) A 3D geometrical model of a normal vertebrobasilar artery reconstructed from a patient's MRA. (G–J) The hemodynamic distribution (velocity, TAWSS OSI, and TAWSSG) of a normal vertebrobasilar artery.
BAF Modeling
Three-dimensional (3D) geometry and morphometry of vertebrobasilar arteries were extracted from 3D-RA images based on thresholding using the MIMICS software (Materialize, NV, Belgium). The geometries were repaired, cut, smoothed using Geomagic Studio 2012 (3D Systems, North Carolina, USA). Then 3D models were imported into ICEM CFD (ANSYS Inc., Canonsburg, PA, United States) to generate unstructured computational meshes for simulation. All models grid ranged from 1.39 to 1.46 million elements. The wall of the vessel was assumed rigid with no-slip conditions, blood was assumed to be a laminar, incompressible Newtonian fluid. The material blood parameters were set as follows: density, P = 1,066 kg/m3; dynamic viscosity, u = 0.00345 Pa × s. ANSYS CFX 19.2 (ANSYS Inc., Canonsburg, PA, Unites States) was utilized to solve Navier-Stokes equations. The pulsatile flow was imposed at the inlet using the Womersley profile as the boundary condition. The flow waveform was derived from a patient-specific transcranial doppler (TCD) ultrasound of the vertebral artery. Outlet boundaries were set using Murray's law (18). A numerical simulation was carried out for three cardiac cycles, with a time step of 0.01 s; the results of the last cycle were used for the analysis.
Hemodynamic Analysis
Hemodynamic parameters, time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), and TAWSS gradient (TAWSSG) were computed within the flow confluence of the vertebral artery. We further calculated surface area ratio-TAWSS (SAR-TAWSS), surface area ratio-OSI (SAR-OSI), and surface area ratio-TAWSSG (SAR-TAWSSG) in the bifurcation, basilar artery of the vertebrobasilar artery (Figure 1). Nomenclatures analyzed were as follows:
1. SAR-TAWSS in the bifurcation:
( × 100%).
2. SAR-TAWSS in the basilar artery:
( × 100%).
3. SAR-OSI in the bifurcation:
( × 100%).
4. SAR-OSI in the basilar artery:
( × 100%).
5. SAR-TAWSSG in the bifurcation:
( × 100%).
6. SAR-TAWSSG in the basilar artery:
( × 100%).
Statistical Analysis
For baseline characteristics, quantitative variables are described as mean ± SD and qualitative variables are described as numbers and percentages. Given the non-normal distribution, the hemodynamic characteristics are presented as medians and interquartile ranges. Hemodynamic parameters of flow patterns and geometric variables were compared between BAF and control groups. Comparison of categorical variables was performed using the χ2 test (or Fisher's exact test as appropriate). Student's t-test or Wilcoxon rank-sum tests were used to comparing quantitative variables. A two-tailed p < 0.05 was considered significant.
Results
Patient and Fenestration Characteristics
The mean ± SD of the patient age in the BAF and control groups were 55.58 ± 12.24 and 55.25 ± 9.64, respectively (p = 0.942). No significant differences were found in the baseline characteristics between the two groups, including sex distribution (p = 1), smoking history (p = 1), posterior circulation infarction (p = 0.478), blood pressure (systolic, p = 0.839, diastolic, p = 0.708), LDL (p = 0.27), HDL (p = 0.302), fasting glucose (p = 0.22), cholesterol (p = 0.379), and triglyceride (p = 0.818). These data are summarized in Table 1.
BAF Group vs. Control Group
The values of shape, geometric features, and hemodynamics of the BAF and normal vertebrobasilar artery are provided in Table 2. Patients in the BAF group had a smaller basilar artery diameter than that those in the control group (3.1 ± 0.51 vs. 3.76 ± 0.4, p = 0.002). No significant differences were found between the two groups regarding any other geometric or shape parameters, including VA diameter, bifurcation angle, basilar artery curvature, and basilar artery length. When comparing hemodynamic parameters between the two groups, maximum OSI (median, 0.3 vs. 0.09, p = 0.028), TAWSSG (median, 983.42 vs. 565.39, p = 0.038) was higher at the flow confluence of BAF than the normal vertebrobasilar artery. SAR-TAWSSG was higher in bifurcation (median, 70.22 vs. 27.65, p = 0.002) and basilar artery (median, 48.75 vs. 16.17, p < 0.001) of BAF compared with normal vertebrobasilar artery. Data summarizing these are shown in Figure 2.
Figure 2. Boxplots displaying the following hemodynamic parameters: surface area ratio-TAWSSG (SAR-TAWSSG), OSI, and TAWSSG for the patients with BAF and those with a normal vertebrobasilar artery (control group). The P-values were obtained from the Wilcoxon rank-sum test. The BAF group had a significantly higher SAR-TAWSSG value in the bifurcation (A) and basilar artery zones (B) than the control group. The BAF group had significantly higher maximum OSI (C) and TAWSSG values (D) than the control group in flow confluence.
Discussion
This study provides morphometry and flow patterns associated with BAF in vertebrobasilar arteries through CFD. We observed that patients with BAF had higher OSI, TAWSSG, and SAR-TAWSSG values than those of control subjects with normal vertebrobasilar arteries. Furthermore, at the flow confluence, higher OSI and TAWSSG were observed in the proximal segment of patients with BAF. These results verified that fenestration could disturb the hemodynamic patterns in relevant anomalous arteries. Thus, more attention should be paid to this rare congenital vascular disease in clinical settings.
Various CFD parameters are associated with atherosclerosis (14, 19). WSS is used to describe the tangential friction of blood flow on the vascular wall and is known to affect the biological function of endothelial cells (20, 21). Low WSS stimulates atherosclerosis in its early stage, while high WSS can cause endothelial trauma, modulate plaque composition, and reduce cap stability thereby promoting acute thrombus formation upon plaque rupture (22, 23). High WSS is usually accompanied by spatial WSSG in complex vascular geometries. Endothelial cells are more sensitive to WSSG and have different response gradient signals (24). Therefore, in addition to high WSS, WSSG may contribute to alterations in the vulnerable plaque pathobiology. A recent study observed that patients with vertebral artery stenosis had significantly high SAR-TAWSSG and a high incidence of transient ischemic attacks (TIAs) or posterior circulation infarction. Further, high SAR-TAWSSG may induce TIA or posterior circulation infarction (25). In our study, we found that higher SAR-TAWSSG both at the bifurcation and basilar artery zone of patients with BAF compared to that of those patients with normal vertebrobasilar arteries. A previous study reported that BAF is a risk factor for intracranial atherosclerotic plaque, and plaques are mostly located in the proximal and bifurcation segments of BAF (26). High SAR-TAWSSG values were found in the bifurcation of patients with BAF in our study, which may be the cause of unstable plaque formation. The basilar artery is also prone to plaque formation with many perforating arteries. As hemodynamics change, BAF-induced disturbance of flow may further precipitate diffuse diseases in this vessel. Thus, both high SAR-TAWSSG and diffuse vascular diseases may increase the instability of basilar artery plaque, resulting in TIA or infarctions. However, this requires further investigation.
Our study also demonstrates that lower TAWSS is usually located at the outer wall and the proximal as well as bifurcation segments of a BAF (Figure 1). This hemodynamic phenomenon is similar to that in the branches of coronary arteries. According to hydrodynamic theory, WSS in the inner/medial wall is higher than at other components of the vessels. In contrast, the stress in the outer/lateral wall is lower, and its distribution is more irregular (27). We calculated the hemodynamic parameters on the surface of the flow confluence segment of vertebrobasilar arteries. The flow confluence of the BAF converges at the bilateral vertebral artery and then separates into two channels. In a short distance, flow velocity changes rapidly and tend to be characterized by the disturbed flow. Thus, it may create large spatial variations in WSS, which increases the risk of thrombosis formation and subsequent ischemic stroke. BAF-related infarction may be associated with local thrombosis has been reported (13). Further, an artery obstructed by a thrombus or plaque debris lodged is both scenarios that may lead to ischemic stroke. In our study, higher OSI and TAWSSG values at the flow confluence of the BAF were noted, especially in the middle sidewall (Figure 1). Plaque development is associated with low and oscillating WSS, and OSI can quantify the changes in WSS direction over the cardiac cycle. OSI values are associated with changes in endothelial gene expression, lipid accumulation, and inflammatory cell activation. Therefore, a high OSI value is considered an indicator of vascular endothelial dysfunction (28). We speculate that higher OSI and TAWSSG of the BAF can result in multidirectional blood flow disorders, which could easily cause vessel wall dysfunction, promote platelet aggregation, and eventually result in thrombus formation and plaque instability (27). Therefore, we speculate that the flow confluence segments of the BAF are high-risk locations for thrombus formation, plaque development, and rupture. Even though our study excluded BAF with stenosis or aneurysm, a numerically higher rate of posterior circulation infarction was observed in the BAF group (16.6 vs. 0%). Future studies involving more BAF patients with symptoms or serial observation of hemodynamic characteristics and the correlation with cerebrovascular events are warranted.
Limitations
This study has several limitations. First, with the low prevalence of BAF, only 12 patients were included over the study period. Although there is a trend in the results, a small sample may have weakened our conclusion. Second, CFD is a commonly employed method used to obtain hemodynamic parameters but does inevitably involve some assumptions. Modeling assumptions may result in deviation from the actual flow as it would occur in vivo. Finally, a follow-up study with a larger sample size is required to more clearly establish a correlation between BAF, vessel wall changes, and the occurrence of cerebrovascular events.
Conclusions
This pilot study suggests that there may be hemodynamic differences between BAF and normal vertebrobasilar systems. BAF results in high OSI, TAWSSG, and SAR-TAWSSG. Turbulent flow in the BAF might induce thrombus and unstable plaque formation, which increases the risk of ischemic stroke. Combined with high-resolution magnetic resonance imaging for vascular wall imaging, the hemodynamic characteristics of BAF may be more comprehensively evaluated in future studies.
Data Availability Statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding authors.
Author Contributions
LJ and DC contributed to the conception and design of the study. JD drafted the initial manuscript. YM and XT performed the simulation. XB and AD wrote sections of the manuscript. BY, TW, and AP critically reviewed and revised the manuscript. XY, ML, and RY organized the database. All the authors reviewed and approved the final version of the manuscript.
Funding
This work was supported by the Xuanwu Hospital Science Program for Fostering Young Scholars (QNPY2020003) and Beijing Natural Science Foundation (Z190014). The funders have no role in study design, data analysis, and writing for the manuscript.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The authors would like to thank Shukun (Beijing) Network Technology Co., Ltd., for providing additional computing resources and technical support.
References
1. Wu aL X, Zhu J, Cai B. Basilar artery fenestration: an unusual possible cause of ischaemic stroke? BMJ Case Rep. (2018) 2018:bcr2017222910. doi: 10.1136/bcr-2017-222910
2. Gold JJ, Crawford JR. An unusual cause of pediatric stroke secondary to congenital basilar artery fenestration. Case Rep Crit Care. (2013) 2013:627972. doi: 10.1155/2013/627972
3. Wollschlaeger G, Wollschlaeger PB, Lucas FV, Lopez VF. Experience and result with postmortem cerebral angiography performed as routine procedure of the autopsy. Am J Roentgenol Radium Ther Nucl Med. (1967) 101:68–87. doi: 10.2214/ajr.101.1.68
4. Sanders WP, Sorek PA, Mehta BA. Fenestration of intracranial arteries with special attention to associated aneurysms and other anomalies. AJNR Am J Neuroradiol. (1993) 14:675−80.
5. Uchino A, Kato A, Takase Y, Kudo S. Basilar artery fenestrations detected by MR angiography. Radiat Med. (2001) 19:71−4.
6. Uchino A, Saito N, Okada Y, Kozawa E, Nishi N, Mizukoshi W, et al. Fenestrations of the intracranial vertebrobasilar system diagnosed by MR angiography. Neuroradiology. (2012) 54:445–50. doi: 10.1007/s00234-011-0903-x
7. Tanaka M, Kikuchi Y, Ouchi T. Neuroradiological analysis of 23 cases of basilar artery fenestration based on 2280 cases of MR angiographies. Interv Neuroradiol. (2006) 12:39–44. doi: 10.1177/15910199060120S103
8. Finlay HM, Canham PB. The layered fabric of cerebral artery fenestrations. Stroke. (1994) 25:1799–806. doi: 10.1161/01.STR.25.9.1799
9. Kathuria S, Gregg L, Chen J, Gandhi D. Normal cerebral arterial development and variations. Semin Ultrasound CT MR. (2011) 32:242–51. doi: 10.1053/j.sult.2011.02.002
10. Kloska SP, Schlegel PM, Strater R, Niederstadt TU. Causality of pediatric brainstem infarction and basilar artery fenestration? Pediatr Neurol. (2006) 35:436–8. doi: 10.1016/j.pediatrneurol.2006.06.013
11. Palazzo P, Ruff M, Lyerly MJ, Alexandrov AV. Basilar artery thrombus vs. fenestration: a differential diagnostic challenge in acute ischemic stroke. J Neuroimaging. (2014) 24:607–9. doi: 10.1111/jon.12069
12. Trivelato FP, Abud DG, Nakiri GS, de Castro Afonso LH, Ulhoa AC, Manzato LB, et al. Basilar artery fenestration aneurysms: endovascular treatment strategies based on 3D morphology. Clin Neuroradiol. (2016) 26:73–9. doi: 10.1007/s00062-014-0336-0
13. Miyamoto N, Ueno Y, Hira K, Kijima C, Nakajima S, Yamashiro K, et al. Characteristics of clinical symptoms, cerebral images and stroke etiology in vertebro-basilar artery fenestration-related infarction. Brain Sci. (2020) 10. doi: 10.3390/brainsci10040243
14. Malek AM, Alper SL, Izumo S. Hemodynamic shear stress and its role in atherosclerosis. JAMA. (1999) 282:2035–42. doi: 10.1001/jama.282.21.2035
15. Sui B, Gao P, Lin Y, Gao B, Liu L, An J. Blood flow pattern and wall shear stress in the internal carotid arteries of healthy subjects. Acta radiol. (2008) 49:806–14. doi: 10.1080/02841850802068624
16. Jeong SK, Kwak HS, Cho YI. Middle cerebral artery fenestration in patients with cerebral ischemia. J Neurol Sci. (2008) 275:181–4. doi: 10.1016/j.jns.2008.07.037
17. Tong X, Dong J, Zhou G, Zhang X, Wang A, Ji Z, et al. Hemodynamic effects of size and location of basilar artery fenestrations associated to pathological implications. Int J Numer Method Biomed Eng. (2021) 37:e3507. doi: 10.1002/cnm.3507
18. Murray CD. The physiological principle of minimum work applied to the angle of branching of arteries. J Gen Physiol. (1926) 9:835–41. doi: 10.1085/jgp.9.6.835
19. Carallo C, Lucca LF, Ciamei M, Tucci S, de Franceschi MS. Wall shear stress is lower in the carotid artery responsible for a unilateral ischemic stroke. Atherosclerosis. (2006) 185:108–13. doi: 10.1016/j.atherosclerosis.2005.05.019
20. Dolan JM, Kolega J, Meng H. High wall shear stress and spatial gradients in vascular pathology: a review. Ann Biomed Eng. (2013) 41:1411–27. doi: 10.1007/s10439-012-0695-0
21. Meng H, Tutino VM, Xiang J, Siddiqui A. High WSS or low WSS? Complex interactions of hemodynamics with intracranial aneurysm initiation, growth, and rupture: toward a unifying hypothesis. AJNR Am J Neuroradiol. (2014) 35:1254–62. doi: 10.3174/ajnr.A3558
22. Eshtehardi P, Brown AJ, Bhargava A, Costopoulos C, Hung OY, Corban MT, et al. High wall shear stress and high-risk plaque: an emerging concept. Int J Cardiovasc Imaging. (2017) 33:1089–99. doi: 10.1007/s10554-016-1055-1
23. Gijsen FJ, Mastik F, Schaar JA, Schuurbiers JC, van der Giessen WJ, de Feyter PJ, et al. High shear stress induces a strain increase in human coronary plaques over a 6-month period. EuroIntervention. (2011) 7:121–7. doi: 10.4244/EIJV7I1A20
24. Tzima E, Irani-Tehrani M, Kiosses WB, Dejana E, Schultz DA, Engelhardt B, et al. A mechanosensory complex that mediates the endothelial cell response to fluid shear stress. Nature. (2005) 437:426–31. doi: 10.1038/nature03952
25. Feng Y, Liu J, Fan T, Zhang W, Yin X. E Y, et al. Vertebral artery stenoses contribute to the development of diffuse plaques in the basilar artery. Front Bioeng Biotechnol. (2020) 8:168. doi: 10.3389/fbioe.2020.00168
26. Liu L, Zhang XB, Lu S, Liu ZJ, Zhu XJ. Plaque distribution of basilar artery fenestration by 3D high-resolution mr vessel wall imaging. Cell Transplant. (2019) 28:851–5. doi: 10.1177/0963689719843813
27. Soulis JV, Giannoglou GD, Parcharidis GE, Louridas GE. Flow parameters in normal left coronary artery tree. Implication to atherogenesis. Comput Biol Med. (2007) 37:628–36. doi: 10.1016/j.compbiomed.2006.06.006
Keywords: basilar artery fenestration, ischemic stroke, hemodynamic, wall shear stress, computational fluid dynamics
Citation: Dong J, Mei Y, Bai X, Tong X, Dmytriw AA, Yang B, Wang T, Patel AB, Yang X, Li M, Yang R, Chen D and Jiao L (2022) Hemodynamic Differences Between Basilar Artery Fenestration and Normal Vertebrobasilar Artery: A Pilot Study. Front. Neurol. 12:766174. doi: 10.3389/fneur.2021.766174
Received: 28 August 2021; Accepted: 06 December 2021;
Published: 13 January 2022.
Edited by:
Kenneth Liu, Bronson Neuroscience Center, United StatesReviewed by:
Benjamin Yim, Stanford University, United StatesVikas Gupta, Texas Tech University Health Sciences Center El Paso, United States
Copyright © 2022 Dong, Mei, Bai, Tong, Dmytriw, Yang, Wang, Patel, Yang, Li, Yang, Chen and Jiao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Liqun Jiao, bGlxdW5qaWFvJiN4MDAwNDA7c2luYS5jbg==; Duanduan Chen, ZHVhbmR1YW4mI3gwMDA0MDtiaXQuZWR1LmNu
†These authors have contributed equally to this work and share first authorship