Skip to main content

ORIGINAL RESEARCH article

Front. Physiol., 01 April 2021
Sec. Integrative Physiology
This article is part of the Research Topic Physiological and Pathological Responses to Hypoxia and High Altitude, Volume II View all 22 articles

The Association Between Notching of the Right Ventricular Outflow Tract Flow Velocity Doppler Envelope and Impaired Right Ventricular Function After Acute High-Altitude Exposure

  • 1Institute of Cardiovascular Diseases of PLA, The Second Affiliated Hospital, Third Military Medical University (Army Medical University), Chongqing, China
  • 2Department of Cardiology, The Second Affiliated Hospital, Third Military Medical University (Army Medical University), Chongqing, China
  • 3Department of Geriatric Cardiology, Chinese PLA General Hospital, Beijing, China
  • 4Department of Medical Ultrasonics, The Second Affiliated Hospital, Third Military Medical University (Army Medical University), Chongqing, China

Introduction: Pulmonary artery pressure (PAP) is increased and right ventricular (RV) function is well preserved in healthy subjects upon exposure to high altitude (HA). An increase in PAP may trigger notching of the right ventricular outflow tract Doppler flow velocity envelope (RVOT notch), which is associated with impaired RV function in patients with pulmonary hypertension. However, whether HA exposure can induce RVOT notch formation and the subsequent impact on cardiac function in healthy subjects remains unclear.

Methods: A total of 99 subjects (69 males and 30 females) with a median age of 25 years were enrolled in this study; they traveled from 500 to 4100 m by bus over a 2-day period. All subjects underwent a comprehensive physiological and echocardiographic examination 1 day before ascension at low altitude and 15 ± 3 h after arrival at HA. The RVOT notch was determined by the presence of a notched shape in the RVOT Doppler flow velocity envelope. The systolic PAP (SPAP) was calculated as Bernoulli equation SPAP = 4 × (maximum tricuspid regurgitation velocity)2+5 and mean PAP (mPAP) = 0.61 × SPAP+2. Cardiac output was calculated as stroke volume × heart rate. Pulmonary capillary wedge pressure (PCWP) was calculated as 1.9+1.24 × mitral E/e’. Pulmonary vascular resistance (PVR) was calculated as (mPAP-PCWP)/CO.

Results: After HA exposure, 20 (20.2%) subjects had an RVOT notch [notch (+)], and 79 (79.8%) subjects did not have an RVOT notch [notch (−)]. In the multivariate logistic regression analysis, the SPAP, right ventricular global longitude strain (RV GLS), and tricuspid E/A were independently associated with the RVOT notch. The SPAP, mPAP, PVR, standard deviations of the times to peak systolic strain in the four mid-basal RV segments (RVSD4), peak velocity of the isovolumic contraction period (ICV), and the peak systolic velocity (s’) at the mitral/tricuspid annulus were increased in all subjects. Conversely, the pulse oxygen saturation (SpO2), RV GLS, and tricuspid annulus plane systolic excursion (TAPSE)/SPAP were decreased. However, the increases of SPAP, mPAP, PVR, and RVSD4 and the decreases of SpO2, RV GLS, and TAPSE/SPAP were more pronounced in the notch (+) group than in the notch (−) group. Additionally, increased tricuspid ICV and mitral/tricuspid s’ were found only in the notch (−) group.

Conclusion: HA exposure-induced RVOT notch formation is associated with impaired RV function, including no increase in the tricuspid ICV or s’, reduction of RV deformation, deterioration in RV-pulmonary artery coupling, and RV intraventricular synchrony.

Introduction

Travelers typically flock to high-altitude (HA) areas. However, hypobaric hypoxia at HA stresses the cardiopulmonary system (Parati et al., 2018). Due to HA-induced hypoxic pulmonary vasoconstriction (HPV), the pulmonary artery pressure (PAP) increases (Naeije, 2019), filling of the left ventricle and right ventricle decreases, and the left ventricular (LV)/right ventricular (RV) diastolic function is altered (Maufrais et al., 2017). In contrast, the LV/RV contractile function remains well-preserved (Osculati et al., 2016; Maufrais et al., 2019; Williams et al., 2019; Sareban et al., 2020). According to previous studies, mild pulmonary hypertension (PH) could decrease RV function, and the RV dysfunction progresses regardless of the treatment of lowering pulmonary vascular resistance (PVR) (van de Veerdonk et al., 2011; Huston et al., 2019).

Notching of the right ventricular outflow tract flow velocity Doppler envelope (RVOT notch) has been widely reported in patients with PH. The RVOT notch is caused by the reflected wave propagating to the pulmonary valve prior to the closure of the valve during a systole. The RVOT notch forms in the context of increased PAP and PVR, decreased pulmonary vascular compliance, and the presence of a reflecting site close enough to the pulmonary valve to allow for the velocity of the reflecting wave (Arkles et al., 2011; Ghimire et al., 2016).

The formation of the RVOT notch was mainly attributed to high PAP (Lopez-Candales and Edelman, 2012). As for patients with PH, it has been reported that the presence of an RVOT notch is associated with deteriorated RV function, such as reduced tricuspid annulus plane systolic excursion (TAPSE), RV fractional area change (FAC), peak systolic velocity of the tricuspid annulus (tricuspid s’), and the ratio of stroke volume index to RV end-diastolic area (Arkles et al., 2011; Lopez-Candales and Edelman, 2012). Even though PAP also increased in healthy subjects upon acute HA exposure, the amount of increase only ranged from mild to moderate. It is unclear whether this moderate increase of PAP induced by HPV led to RVOT notch formation. Additionally, since the RV function in healthy subjects upon acute HA exposure was reported to be well preserved; whether the presence of the RVOT notch in healthy subjects upon acute HA exposure is associated with deteriorated RV function is also unclear. Thus, the mechanism of RVOT notch formation and the relationship between RVOT notch and RV function at HA are worthy to investigate.

In a previous study, we had preliminarily reported that compared to the RVOT notch negative population, the correlation between PAP and intraventricular RV dyssynchrony was much greater in people with an RVOT notch. In this study, we tended to formally characterize the RVOT notch and comprehensively unravel the association between HA-induced RVOT notch and cardiac function.

Materials and Methods

Participants and Study Design

This study was approved by the Clinical Research Ethics Board of Army Medical University (Identification Code 201907501) and was registered at www.chictr.org.cn (ChiCTR-TRC-No.1900025728). In June 2019, we conducted this prospective cohort study on the Qinhai–Tibet Plateau. Han Chinese people who were born and permanently lived at a low altitude (≤500 m) were invited to participate in this study. A total of 111 participants were recruited at low altitudes (Chengdu, China, 500 m) to ascend to HA (Litang, China, 4100 m). Of all 111 subjects, five subjects withdrew from the study prior to ascending, and two subjects returned to low altitude because of severe HA illness. Ultimately, 104 subjects successfully ascended to 4100 m (Litang, China) by bus over a 2-day period. Due to the lack of suitable echocardiographic images in five subjects, the final analysis included 99 subjects. These subjects were healthy adults without HA exposure history in the past 6 months and had no recent medication use. Exclusion criteria included the presence of cardiopulmonary disease, cerebrovascular disease, liver disease, kidney disease, or malignant tumors. All physiological and transthoracic echocardiography examinations were performed 1 day prior to ascension at low altitude and 15 ± 3 h after arrived at HA. This study was conducted in accordance with the Declaration of Helsinki, and all subjects provided written informed consent.

Assessment of Physiological Parameters

Pulse oxygen saturation (SpO2) was recorded using a pulse oximeter (ONYX OR9500, Nonin, Plymouth, MN, United States). Blood pressure was measured with the subject in a sitting position using an electronic sphygmomanometer (Omron HEM-6200, Japan) after a 10-minute rest period. Heart rate (HR) data were recorded using the synchronous electrocardiogram during echocardiography. We measured height and weight with a height and weight scale (RGZ-120, I WISH, China) 1 day before ascension at a low altitude. The body mass index (BMI) was calculated as weight/(height)2.

Transthoracic Echocardiography

Echocardiographic examinations were performed by two experienced sonographers equipped with a 2.5 MHz adult transducer using a CX50 ultrasound system (Philips Ultrasound System, Andover, MA, United States) with the subject in the left lateral decubitus position. The dynamic echocardiographic images, which consisted of three consecutive cardiac cycles, were stored digitally for offline analysis using QLAB 10.5 (Philips Healthcare, Andover, MA, United States) in a blinded fashion. All examinations and measurements were performed according to the recommendations of the American Society of Echocardiography (Rudski et al., 2010; Mitchell et al., 2019), and the value of each echocardiographic parameter was averaged from measurements of three consecutive cardiac cycles.

LV end-diastolic volume, LV end-systolic volume, RV basal transverse diameter (RVD base), and RVD mid were measured. RV end-diastolic and end-systolic areas were determined using the apical four-chamber view by manually tracing the RV endocardium. TAPSE was recorded as the peak excursion of the lateral tricuspid annulus, which was measured with M-mode echocardiography. The maximum tricuspid regurgitation velocity (TRV) was evaluated using continuous-wave Doppler. The LV and RV inflow were measured using the apical four-chamber view of pulsed-wave Doppler. The maximum early (E) and late (A) diastolic velocities were recorded. A pulse tissue Doppler was used to assess the velocities of the mitral and tricuspid annuli on the lateral and septal aspects in early diastole and systole. The peak early diastolic velocity of the annulus (e’), peak systolic velocity of the annulus (s’), and peak velocity of the isovolumic contraction period (ICV) were recorded as the average of the lateral and septal values. Stroke volume (SV), left ventricular ejection fraction (LVEF), and cardiac output (CO) were calculated using LV volume data and HR (Lang et al., 2015). RV FAC was calculated as follows: FAC = (end-diastolic area - end-systolic area)/end-diastolic area × 100. Systolic pulmonary artery pressure (SPAP) was calculated with the following modified Bernoulli equation: SPAP = 4 × (TRV)2+5 (Yock and Popp, 1984). The mean PAP (mPAP) was calculated as follows: mPAP = 0.61 × SPAP+2 (Bossone et al., 2013). Pulmonary capillary wedge pressure (PCWP) was calculated using the following equation: PCWP = 1.9+1.24 × mitral E/e’ (Bossone et al., 2013). PVR was derived according to the equation PVR = (mPAP-PCWP)/CO (Bossone et al., 2013).

To investigate the RVOT blood flow, the pulse Doppler sample volume probe was placed 0.5–1 cm proximally to the pulmonic valve in the parasternal short-axis view. The RVOT acceleration time (ATRVOT) and RVOT ejection time (ETRVOT) were measured. The shapes of the RVOT Doppler flow velocity envelope were categorized as notch (+) (in the presence of a notching pattern) and notch (−) (when no notching pattern was present), as previously described (Arkles et al., 2011). The specific notch relevant parameter mid-systolic flow deceleration time (mid-systolic DT) was measured (Takahama et al., 2017).

2D-STE Analysis

Standard 2D grayscale images (frame rate > 60 fps) in the apical two-, three-, and four-chamber views were used for the LV and RV speckle tracking analyses. Apical, basal, and mid-level parasternal short-axis views of three consecutive beats were also recorded. Using QLAB software v10.5, the operator manually adjusted the region of interest to include all segments of the ventricular myocardium. These were then traced by an algorithm, which automatically calculated the peak strain of each segment on different images.

The LV global longitudinal strain (GLS) was defined as the mean of 17 segments that was calculated on the apical two-, three- and four-chamber images. The LV global circumferential strain (GCS) was defined as the average peak strain in each segment of the basal, mid-, and apical levels using parasternal short-axis images. The RV GLS was defined as the mean of six segments on the RV-focused four-chamber image.

To analyze LV/RV dyssynchrony, the time from the onset of the QRS complexes to the peak strain of each segment on ECG was calculated and correlated to the R-R interval according to Bazett’s formula: corrected interval = measured time to peak strain/(RR interval)1/2. The standard deviation (SD) of the correlated interval for the 12 mid-basal LV segments was termed the TS-12SD (Zoroufian et al., 2014), whereas the SD for the four mid-basal RV segments was termed the RV-SD4 (Badagliacca et al., 2015).

Reproducibility

To determine the intraobserver variability, data from 10 subjects at a low altitude and 10 subjects at a HA were randomly selected and evaluated twice by the same observer with a 1-month interval between evaluations. Interobserver variability was determined by comparing the evaluations conducted by two separate observers blinded to each other. The intra- and interobserver variabilities for the echocardiographic variables are presented in Supplementary Table S1.

Statistical Analysis

All statistical analyses were performed using SPSS software (Version 22.0, IBM Corp., Armonk, NY, United States). Continuous variables are presented as mean ± standard error of mean if normally distributed or as median and interquartile range if not normally distributed, according to the Kolmogorov-Smirnov test. The paired t-test or Wilcoxon test was used to analyze the differences between paired variables at low altitude and HA, whereas the unpaired variables in the notch (−) and notch (+) groups at low altitude or HA and the unpaired variables in notch (+) subjects with mid-systolic DT < 120 ms and with mid-systolic DT > 120 ms were compared using the independent sample t-test or Mann-Whitney U test. Univariate logistic regression analysis was used to identify the relevant variables associated with RVOT notch. The variables with P-value < 0.05 were included in the multivariate logistic regression analysis. Category variables are presented as a number (percentage) and compared using the chi-square test. Statistical significance was set at P < 0.05.

Results

Subject Characteristics

The 99 subjects consisted of 69 men and 30 women with a median age of 25 years and an average BMI of 22.2 kg/m2. No subject was found to have an RVOT notch at low altitude. Upon HA exposure, an RVOT notch was identified in 20 subjects (20.2%) (Figure 1). There were no significant differences in age, sex, or BMI between the notch (+) and notch (−) groups. Following HA exposure, subjects in all groups had increased systolic blood pressure, diastolic blood pressure, and HR. Both groups also had decreased SpO2 following HA exposure, with the notch (+) group experiencing a lower SpO2 level than the notch (−) group (P = 0.039) (Table 1 and Figure 2).

FIGURE 1
www.frontiersin.org

Figure 1. Occurrence of RVOT notch in healthy subjects at low altitude and upon acute HA exposure. (A) The occurrence of RVOT notch in healthy subjects was 0% at low altitude. Following HA exposure, the occurrence of RVOT notch in healthy subjects was 20.2%. (B) Morphology of the RVOT Doppler flow velocity envelope with or without RVOT notch. (1) RVOT Doppler flow velocity envelope in healthy subjects at low altitude was parabolic in shape. (2) Healthy subjects upon HA exposure with notched RVOT Doppler flow velocity envelope.

TABLE 1
www.frontiersin.org

Table 1. Physiologic parameters in the subjects with or without RVOT notch at low altitude and high altitude.

FIGURE 2
www.frontiersin.org

Figure 2. Differences in SpO2, SPAP, and TAPSE/SPAP between notch (−) and notch (+) groups at low altitude and HA. Different levels of (A) SpO2, (B) SPAP, and (C) TAPSE/SPAP in total subjects, notch (−) subjects and notch (+) subjects upon HA exposure are shown. Differences between notch (−) subjects and notch (+) subjects in (A) SpO2, (B) SPAP, and (C) TAPSE/SPAP were compared at low altitude and HA. Values are presented as mean and standard error of mean. *P < 0.05, **P < 0.01. Abbreviations as in Tables 1, 3.

LV Function

After HA exposure, SV and LV GLS decreased, GCS did not change, and the LVEF, CO, mitral ICV, and TS-12SD increased in all subjects. The mitral s’ in the notch (−) group increased upon HA exposure, but the mitral s’ in the notch (+) group did not change. Following HA exposure, the mitral E, mitral E/A, and mitral E/e’ decreased in all subjects. However, the mitral A increased in all subjects following HA exposure, with a higher value in the notch (+) group than in the notch (−) group. The mitral E/A at HA in the notch (+) group was lower than that in the notch (−) group (Table 2).

TABLE 2
www.frontiersin.org

Table 2. Left ventricular parameters in subjects with or without RVOT notch at low altitude and high altitude.

Pulmonary Hemodynamics and RV Function

Upon HA exposure, TRV, SPAP, mPAP, and PVR increased in all subjects; however, this increase was greater in notch (+) subjects. The ATRVOT, ETRVOT, ATRVOT/ETRVOT, and PCWP decreased in all subjects, except for ETRVOT in notch (+) subjects. The RVD base did not change in all subjects. The RVD mid increased in all subjects, which was majorly observed in the notch (+) group. The RV FAC, TAPSE, and RV GLS decreased in both groups, and the change in RV GLS was greater in the notch (+) group. The tricuspid s’ and tricuspid ICV increased after HA exposure in the notch (−) group (Figure 3). Additionally, the RV GLS decreased in both groups upon HA exposure and was lower in the notch (+) group than in the notch (−) group (Figure 4). The two-dimensional strain of the middle segment in the right ventricular free wall (2DS RVFW mid) and 2DS RVFW base in the notch (+) group at HA were lower than those in the notch (−) group. The TAPSE/SPAP decreased in all subjects after HA exposure, with a greater decrease in the notch (+) group (Figure 2). The RVSD4 increased in all subjects upon HA exposure, and it was more pronounced in the notch (+) group (Figure 4). After HA exposure, the tricuspid A in the notch (+) group was higher than that in the notch (−) group, and the tricuspid E/A in the notch (+) group was lower than that in the notch (−) group (Table 3). Categorizing notch (+) subjects into two groups according to mid-systolic DT (Takahama et al., 2017), we determined the RVD base and RVD mid were larger in mid-systolic DT < 120 ms group than those in mid-systolic DT > 120 ms group (Supplementary Table S2).

FIGURE 3
www.frontiersin.org

Figure 3. Tricuspid pulse Doppler velocity analysis in subjects with and without RVOT notch. Different changes in tricuspid ICV (A) and s’ (B) in notch (−) subjects and notch (+) subjects upon HA exposure. (C) and (D) Left panel: pulse tissue Doppler velocity curves of the tricuspid annulus in notch (−) and notch (+) subjects at low altitude and HA. Right panel: examples of the pulse TDI at tricuspid annulus in a subject without RVOT notch (C) and a subject with RVOT notch (D). TDI, tissue Doppler imaging; *P < 0.05, **P < 0.01. Other abbreviations as in Table 3.

FIGURE 4
www.frontiersin.org

Figure 4. Evaluation of RV GLS and dyssynchrony in subjects with and without RVOT notch from low altitude to HA. Different levels of (A) RVSD4 and (B) RV GLS in total subjects, notch (−) subjects, and notch (+) subjects upon HA exposure are shown. (C) and (D) Examples illustrating the 2D-STE measurements of RV GLS and RV dyssynchrony in a subject without RVOT notch (C) and a subject with RVOT notch (D) at HA. *P < 0.05, **P < 0.01. Abbreviations as in Table 3.

TABLE 3
www.frontiersin.org

Table 3. Pulmonary circulation and right ventricular parameters in subjects with or without RVOT notch at low altitude and high altitude.

Determinants of RVOT Notch

In multivariate logistic regression analysis, SPAP (OR, 1.14; 95% CI, 1.05 to 1.23; P = 0.001), RV GLS (OR, 0.74; 95% CI, 0.60 to 0.93; P = 0.011), and tricuspid E/A (OR, 0.13; 95% CI, 0.03 to 0.68; P = 0.016) were found to be independently associated with RVOT notch (Supplementary Table S3).

Discussion

Following rapid elevation gain, an RVOT notch was found in 20 (20.2%) healthy subjects. SpO2, RV deformation, the level of RV-PA coupling, and RV intraventricular synchrony were decreased in all subjects, while the RV afterload, biventricular s’, and ICV were increased in all subjects. However, the increase in the RV afterload and the decrease in SpO2, RV deformation, RV-PA coupling, and RV intraventricular synchrony were more pronounced in subjects with an RVOT notch. Moreover, the mitral s’, tricuspid s’, and tricuspid ICV in subjects with an RVOT notch did not increase upon HA exposure. These results suggest that HA exposure induced the RVOT notch formation, which may be related to impaired RV function.

Association of RVOT Notch, PAP, and SpO2 Upon HA Exposure

In previous studies, the RVOT notch was mainly observed in patients with chronic PH, especially in those with pulmonary vascular disease (Arkles et al., 2011; Takahama et al., 2017). Patients with low arterial compliance, high arterial resistance, and/or a pulmonary embolism proximal to the pulmonary valve may have a pulmonary vascular reflected wave on Doppler, which is seen as a notched Doppler flow velocity contour in the RVOT if the reflected wave rapidly propagates to the pulmonary valve prior to its closure (Torbicki et al., 1999; Arkles et al., 2011). In this study, we found that RVOT notch formation in subjects upon HA exposure is associated with higher SPAP, which is in accordance with the previous finding in PH patients (Lopez-Candales and Edelman, 2012). Upon HA exposure, alveolar hypoxia-induced HPV and inhomogeneous pulmonary vasoconstriction facilitated the ventilation–perfusion matching in the regional lung (Sommer et al., 2008). According to previous studies, SpO2 decreases and SPAP increases with rising altitude (Penaloza and Arias-Stella, 2007). Furthermore, in patients receiving bosentan for PH, SpO2 levels significantly increase with an increase in altitude (Modesti et al., 2006; Mellor et al., 2014). Thus, patients with lower SpO2 are more likely to experience HPV of greater magnitude as an attempt to maintain the ventilation–perfusion balance, which leads to higher SPAP and PVR. Although the presence of an RVOT notch in patients with PH indicates a high SPAP (Kubba et al., 2016), the uneven HPV at HA might induce a reflection site proximal to the pulmonary valve, which leads to the formation of an RVOT notch.

LV Function

Similar to previous studies, our results indicated that LV filling is decreased and LV contractile function is preserved in subjects upon HA exposure (Maufrais et al., 2017, 2019). Although the LV filling is decreased, the HR is increased. The maintained LV function might due to the increase of sympathetic activation. According to a previous study, LV function is not related to RVOT notch formation due to the absence of an RVOT notch in left-heart-disease-associated PH (Kushwaha et al., 2016). However, hypoxia may affect LV function in notch (+) subjects with a lower SpO2. Thus, the increase of mitral s’ in the notch (+) group was blunted in this study. At HA, no difference in TS-12SD was found between the two groups in this study, indicating that upon HA exposure, the SpO2 has little effect on LV dyssynchrony.

Impact of RVOT Notch on RV Function Evaluated by Conventional Echocardiography

In our study, we found that the presence of an RVOT notch indicates impaired RV function in healthy subjects upon HA exposure, which is consistent with previous studies that reported advanced RV dysfunction in PH patients with RVOT notch (Arkles et al., 2011; Lopez-Candales and Edelman, 2012). Previous reports of the changes in RV FAC and TAPSE vary from a decrease (Kurdziel et al., 2017; Netzer et al., 2017), to no change (Maufrais et al., 2017, 2019), and to an increase upon HA exposure (Sareban et al., 2020). These inconsistent results may be due to different exposure times and ascending altitudes, or states of dehydration, all of which have been associated with RV adaptation upon HA exposure. In the present study, we found that the TAPSE and RV FAC decreased upon HA exposure. Although RV contractility may increase in response to the acutely increased RV afterload to maintain the pump function and RV-PA coupling (Naeije et al., 2014), the appropriately increased RV contractility upon HA exposure that is mirrored by load-dependent RV FAC and TAPSE may be under detected due to hypovolemia-induced RV FAC and TAPSE reduction via the Frank-Starling mechanism in patients with moderately increased PAP (Rudski et al., 2010; Najafian et al., 2015; Stembridge et al., 2016; Zhao et al., 2019).

TAPSE/SPAP is noninvasively measured for the evaluation of RV-PA coupling and correlates well with the gold standard multi-beat end-systolic/arterial elastance ratio in patients with PH (Tello et al., 2019; Richter et al., 2020). In this study, the decrease of TAPSE/SPAP in all subjects indicates that the level of RV-PA coupling was decreased in healthy subjects upon HA exposure. The more significant decrease of TAPSE/SPAP in the notch (+) group was attributed to the higher SPAP. The comparable TAPSE between notch (−) and notch (+) groups suggests that the RV contractility in the notch (+) group is maladaptive to the higher RV afterload upon HA exposure. TAPSE/SPAP has been reported to be positively correlated with VO2 peak and workload during cardiopulmonary exercise testing (Martens et al., 2018; Tello et al., 2018). Thus, although TAPSE/SPAP was maintained within a normal range, lower TAPSE/SPAP in the notch (+) subjects upon HA exposure may indicate lower work capacity.

Impact of RVOT Notch on RV Function Evaluated by TDI

Upon HA exposure, the load-independent tricuspid s’ and ICV increased in the context of increased PAP (Abali et al., 2005; Rudski et al., 2010; Negishi et al., 2017), which indicates an increased RV inotropic function in response to the elevated RV afterload (Vonk Noordegraaf et al., 2017). This enhancement of myocardial contractility may be attributed to homeometric autoregulation or sympathetic activation (Rex et al., 2007, 2008). The association of RV contractility with sympathetic activation upon HA exposure is supported by a previous report, which suggests that the tricuspid s’ is correlated with HR (Sareban et al., 2020). However, in the notch (+) group, the tricuspid s’ and ICV did not increase to match the increased RV afterload, leading to a blunted RV contractile response upon HA exposure. This may be due to the fact that the positive inotropic effect of homeometric autoregulation or sympathetic activation was overwhelmed by the negative inotropic effect of hypoxia in notch (+) subjects with lower SpO2 (Silverman et al., 1997).

Impact of RVOT Notch on RV Function Evaluated by 2D-STE

Although the conventional RV function parameters RV FAC and TAPSE were comparable in both groups at HA in this study, the 2D-STE RV function parameter RV GLS was lower in the notch (+) group, especially in the base and mid segments of RVFW. This may be due to the fact that the RV GLS is a more sensitive and accurate assessment for detecting subclinical RV dysfunction than TAPSE and RV FAC (Giusca et al., 2010; Li et al., 2018; Tamulenaite et al., 2018). Even in mild PH, when FAC does not change, the RV GLS begins to decrease (Li et al., 2013). Moreover, the correlation of RV GLS with the gold standard cardiac magnetic resonance-derived RV ejection fraction (RVEF) is better than the correlation of RV FAC and TAPSE with RVEF (Tong et al., 2018; Kavurt et al., 2019); a lower RV GLS may indicate a lower RVEF in notch (+) subjects. These findings are consistent with those of a previous study which reported that the presence of an RVOT notch was related to a reduced RV SV index in patients with PH (Arkles et al., 2011). Intraventricular RV dyssynchrony is an early sign of RV contractile dysfunction. In our previous study, we had found that notch (+) cases in RV dyssynchrony (−) subjects and RV dyssynchrony (+) subjects were comparable. However, in this study, we found that the RVOT notch was associated with higher RVSD4, which is the indicator of RVD. This inconsistency may be attributed to the different ascending modes, exposure times, and ratio of males and females. In the present 2-day period HA exposure study, the subjects with an RVOT notch are more prone to suffer intraventricular RV mechanical dyssynchrony (RVSD4 > 18.7 ms) compared with subjects without an RVOT notch (Badagliacca et al., 2015; Yang et al., 2020). In this study, higher RV dyssynchrony was accompanied by a lower RV GLS. A similar observation was reported in patients after hemodialysis (Unlu et al., 2019).

Impact of Mid-Systolic DT on Cardiac Function

Previous studies reported that PH patients with mid-systolic DT < 120 ms had worse clinical outcomes compared with those subjects with mid-systolic DT > 120 ms (Takahama et al., 2017). In this study, the RVD base and RVD mid in mid-systolic DT < 120 ms group were larger than those in mid-systolic DT > 120 ms group. RV dilation is a sign of RV maladaptation in response to high PAP; however, the outcomes of subjects with mid-systolic DT < 120 ms upon HA exposure still need long-term follow-up studies to investigate.

Several limitations in this study should be noted. For ethical reasons, cardiac catheterization, the gold standard for the measurement of cardiac hemodynamics, was not performed in this study. As 30 females (30.3% of total subjects) were included in this study and the menstrual cycle in female participants may affect heart physiology, menstrual cycle information should be collected, and its association with RVOT notch needs to be investigated. Acute mountain sickness was reported to be associated with reduced LV function, whereas, as this study mainly focused on the association of RVOT notch with RV function, acute mountain sickness cases were not excluded. As the subjects in this study returned to low altitude within two months, the long-term association of an RVOT notch with RV function was not studied. It is not clear whether the subjects with RVOT notch are more likely to develop HA heart disease. Furthermore, the cardiopulmonary exercise test to investigate RV contractile reserves in subjects with an RVOT notch should be performed in future studies as the RV function in these subjects is maladapted to HA exposure. It would be valuable to see whether the RVOT notch would reverse to normal after notch (+) subjects returned to sea level and whether a lower exercise capacity could be observed in the notch (+) subjects at sea level via cardiopulmonary exercise test.

Conclusion

Our study revealed that HA exposure-induced RVOT notch formation may be attributed to high PAP and low SpO2. The presence of an RVOT notch upon HA exposure was associated with impaired RV function as the tricuspid ICV and s’ did not increase, the RV deformation reduced, the RV-PA coupling deteriorated, and RV intraventricular dyssynchrony was observed. Therefore, an RVOT notch may be a potential echocardiographic sign representing impaired RV function at HA. However, whether the impaired RV function in subjects with RVOT notch upon HA exposure persists or is transient still needs further investigation.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by the Clinical Research Ethics Board of Army Medical University. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

LH, FY, and CL contributed to the conception or design of the work. FY, CL, JK, YY, CH, CZ, RR, ZL, and JuY conducted the experiments. FY and CL performed the statistical analyses and drafted the manuscript. SY, SB, JiY, XD, JZ, and HT interpreted the results of the statistical analyses. LH critically revised the manuscript. All authors approved of the final version of the manuscript and agreed to be accountable for all aspects of the work.

Funding

This work was supported by grants from the National Natural Science Foundation of China (Grant No. 81730054); Military Logistics Research Project, PLA (Grant No. BLJ18J007); the Special Health Research Project, Ministry of Health of the P.R. China (Grant No. 201002012); and the PLA Youth Training Project for Medical Science (Grant No. 15QNP062).

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.

Acknowledgments

We thank all participants in this study for their contributions to medical science.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2021.639761/full#supplementary-material

References

Abali, G., Tokgozoglu, L., Ozcebe, O. I., Aytemir, K., and Nazli, N. (2005). Which Doppler parameters are load independent? a study in normal volunteers after blood donation. J. Am. Soc. Echocardiogr. 18, 1260–1265. doi: 10.1016/j.echo.2005.06.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Arkles, J. S., Opotowsky, A. R., Ojeda, J., Rogers, F., Liu, T., Prassana, V., et al. (2011). Shape of the right ventricular doppler envelope predicts hemodynamics and right heart function in pulmonary hypertension. Am. J. Respir. Crit. Care Med. 183, 268–276. doi: 10.1164/rccm.201004-0601oc

PubMed Abstract | CrossRef Full Text | Google Scholar

Badagliacca, R., Reali, M., Poscia, R., Pezzuto, B., Papa, S., Mezzapesa, M., et al. (2015). Right intraventricular dyssynchrony in idiopathic, heritable, and anorexigen-induced pulmonary arterial hypertension: clinical impact and reversibility. JACC Cardiovasc. Imaging 8, 642–652. doi: 10.1016/j.jcmg.2015.02.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Bossone, E., D’andrea, A., D’alto, M., Citro, R., Argiento, P., Ferrara, F., et al. (2013). Echocardiography in pulmonary arterial hypertension: from diagnosis to prognosis. J. Am. Soc. Echocardiogr. 26, 1–14. doi: 10.1016/j.echo.2012.10.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghimire, A., Andersen, M. J., Burrowes, L. M., Bouwmeester, J. C., Grant, A. D., Belenkie, I., et al. (2016). The reservoir-wave approach to characterize pulmonary vascular-right ventricular interactions in humans. J. Appl. Physiol. 121, 1348–1353. doi: 10.1152/japplphysiol.00697.2016

PubMed Abstract | CrossRef Full Text | Google Scholar

Giusca, S., Dambrauskaite, V., Scheurwegs, C., D’hooge, J., Claus, P., Herbots, L., et al. (2010). Deformation imaging describes right ventricular function better than longitudinal displacement of the tricuspid ring. Heart 96, 281–288. doi: 10.1136/hrt.2009.171728

PubMed Abstract | CrossRef Full Text | Google Scholar

Huston, J. H., Maron, B. A., French, J., Huang, S., Thayer, T., Farber-Eger, E. H., et al. (2019). Association of mild echocardiographic pulmonary hypertension with mortality and right ventricular function. JAMA Cardiol. 4, 1112–1121. doi: 10.1001/jamacardio.2019.3345

PubMed Abstract | CrossRef Full Text | Google Scholar

Kavurt, A. V., Paç, F. A., Koca, S., Mutlu Mıhçıoðlu, A., and Yiðit, H. (2019). The evaluation of right ventricular systolic function in patients with repaired tetralogy of fallot by conventional echocardiographic methods and speckle tracking echocardiography: compared with the gold standard cardiac mangenetic resonance. Echocardiography 36, 2251–2258. doi: 10.1111/echo.14532

PubMed Abstract | CrossRef Full Text | Google Scholar

Kubba, S., Davila, C. D., and Forfia, P. R. (2016). Methods for evaluating right ventricular function and ventricular-arterial coupling. Prog. Cardiovasc. Dis. 59, 42–51. doi: 10.1016/j.pcad.2016.06.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Kurdziel, M., Wasilewski, J., Gierszewska, K., Kazik, A., Pytel, G., Waclawski, J., et al. (2017). Echocardiographic assessment of right ventricle dimensions and function after exposure to extreme altitude: is an expedition to 8000 m hazardous for right ventricular function? High Alt. Med. Biol. 18, 330–337. doi: 10.1089/ham.2017.0019

PubMed Abstract | CrossRef Full Text | Google Scholar

Kushwaha, S. P., Zhao, Q. H., Liu, Q. Q., Wu, W. H., Wang, L., Yuan, P., et al. (2016). Shape of the pulmonary artery doppler-flow profile predicts the hemodynamics of pulmonary hypertension caused by left-sided heart disease. Clin. Cardiol. 39, 150–156. doi: 10.1002/clc.22493

PubMed Abstract | CrossRef Full Text | Google Scholar

Lang, R. M., Badano, L. P., Mor-Avi, V., Afilalo, J., Armstrong, A., Ernande, L., et al. (2015). Recommendations for cardiac chamber quantification by echocardiography in adults: an update from the American society of echocardiography and the European association of cardiovascular imaging. Eur. Heart J. Cardiovasc. Imaging 16, 233–270. doi: 10.1093/ehjci/jev014

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, A. L., Zhai, Z. G., Zhai, Y. N., Xie, W. M., Wan, J., and Tao, X. C. (2018). The value of speckle-tracking echocardiography in identifying right heart dysfunction in patients with chronic thromboembolic pulmonary hypertension. Int. J. Cardiovasc. Imaging 34, 1895–1904. doi: 10.1007/s10554-018-1423-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y., Xie, M., Wang, X., Lu, Q., and Fu, M. (2013). Right ventricular regional and global systolic function is diminished in patients with pulmonary arterial hypertension: a two-dimensional ultrasound speckle tracking echocardiography study. Int. J. Cardiovasc. Imaging 29, 545–551. doi: 10.1007/s10554-012-0114-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Lopez-Candales, A., and Edelman, K. (2012). Shape of the right ventricular outflow doppler envelope and severity of pulmonary hypertension. Eur. Heart J. Cardiovasc. Imaging 13, 309–316. doi: 10.1093/ejechocard/jer235

PubMed Abstract | CrossRef Full Text | Google Scholar

Martens, P., Verbrugge, F. H., Bertrand, P. B., Verhaert, D., Vandervoort, P., Dupont, M., et al. (2018). Effect of cardiac resynchronization therapy on exercise-induced pulmonary hypertension and right ventricular-arterial coupling. Circ. Cardiovasc. Imaging 11:e007813.

Google Scholar

Maufrais, C., Rupp, T., Bouzat, P., Doucende, G., Verges, S., Nottin, S., et al. (2017). Heart mechanics at high altitude: 6 days on the top of Europe. Eur. Heart J. Cardiovasc. Imaging 18, 1369–1377. doi: 10.1093/ehjci/jew286

PubMed Abstract | CrossRef Full Text | Google Scholar

Maufrais, C., Rupp, T., Bouzat, P., Esteve, F., Nottin, S., Walther, G., et al. (2019). Medex 2015: the key role of cardiac mechanics to maintain biventricular function at high altitude. Exp. Physiol. 104, 667–676. doi: 10.1113/ep087350

PubMed Abstract | CrossRef Full Text | Google Scholar

Mellor, A., Boos, C., Holdsworth, D., Begley, J., Hall, D., Lumley, A., et al. (2014). Cardiac biomarkers at high altitude. High Alt. Med. Biol. 15, 452–458. doi: 10.1089/ham.2014.1035

PubMed Abstract | CrossRef Full Text | Google Scholar

Mitchell, C., Rahko, P. S., Blauwet, L. A., Canaday, B., Finstuen, J. A., Foster, M. C., et al. (2019). Guidelines for performing a comprehensive transthoracic echocardiographic examination in adults: recommendations from the american society of echocardiography. J. Am. Soc. Echocardiogr. 32, 1–64. doi: 10.1016/j.echo.2018.06.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Modesti, P. A., Vanni, S., Morabito, M., Modesti, A., Marchetta, M., Gamberi, T., et al. (2006). Role of endothelin-1 in exposure to high altitude: acute mountain sickness and endothelin-1 (ACME-1) study. Circulation 114, 1410–1416. doi: 10.1161/circulationaha.105.605527

PubMed Abstract | CrossRef Full Text | Google Scholar

Naeije, R. (2019). Pulmonary hypertension at high altitude. Eur. Respir. J. 53:1900985.

Google Scholar

Naeije, R., Brimioulle, S., and Dewachter, L. (2014). Biomechanics of the right ventricle in health and disease (2013 grover conference series). Pulm. Circ. 4, 395–406. doi: 10.1086/677354

PubMed Abstract | CrossRef Full Text | Google Scholar

Najafian, J., Taheri, S., Mahaki, B., Molavi, S., Alami, S., Khalesi, S., et al. (2015). Comparing right ventricular function and pulmonary artery pressure before and shortly after hemodialysis in patients with end-stage renal disease. Adv. Biomed. Res. 4:197.

Google Scholar

Negishi, K., Borowski, A. G., Popoviæ, Z. B., Greenberg, N. L., Martin, D. S., Bungo, M. W., et al. (2017). Effect of gravitational gradients on cardiac filling and performance. J. Am. Soc. Echocardiogr. 30, 1180–1188. doi: 10.1016/j.echo.2017.08.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Netzer, N. C., Strohl, K. P., Hogel, J., Gatterer, H., and Schilz, R. (2017). Right ventricle dimensions and function in response to acute hypoxia in healthy human subjects. Acta Physiol. 219, 478–485. doi: 10.1111/apha.12740

PubMed Abstract | CrossRef Full Text | Google Scholar

Osculati, G., Revera, M., Branzi, G., Faini, A., Malfatto, G., Bilo, G., et al. (2016). Effects of hypobaric hypoxia exposure at high altitude on left ventricular twist in healthy subjects: data from HIGHCARE study on mount everest. Eur. Heart J. Cardiovasc. Imaging 17, 635–643. doi: 10.1093/ehjci/jev166

PubMed Abstract | CrossRef Full Text | Google Scholar

Parati, G., Agostoni, P., Basnyat, B., Bilo, G., Brugger, H., Coca, A., et al. (2018). Clinical recommendations for high altitude exposure of individuals with pre-existing cardiovascular conditions: a joint statement by the european society of cardiology, the council on hypertension of the european society of cardiology, the european society of hypertension, the international society of mountain medicine, the italian society of hypertension and the italian society of mountain medicine. Eur. Heart J. 39, 1546–1554.

Google Scholar

Penaloza, D., and Arias-Stella, J. (2007). The heart and pulmonary circulation at high altitudes: healthy highlanders and chronic mountain sickness. Circulation 115, 1132–1146. doi: 10.1161/circulationaha.106.624544

PubMed Abstract | CrossRef Full Text | Google Scholar

Rex, S., Missant, C., Segers, P., and Wouters, P. F. (2007). Thoracic epidural anesthesia impairs the hemodynamic response to acute pulmonary hypertension by deteriorating right ventricular–pulmonary arterial coupling. Crit. Care Med. 35, 222–229. doi: 10.1097/01.ccm.0000250357.35250.a2

CrossRef Full Text | Google Scholar

Rex, S., Missant, C., Segers, P., Rossaint, R., and Wouters, P. F. (2008). Epoprostenol treatment of acute pulmonary hypertension is associated with a paradoxical decrease in right ventricular contractility. Intensive Care Med. 34, 179–189. doi: 10.1007/s00134-007-0831-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Richter, M. J., Peters, D., Ghofrani, H. A., Naeije, R., Roller, F., Sommer, N., et al. (2020). Evaluation and prognostic relevance of right ventricular-arterial coupling in pulmonary hypertension. Am. J. Respir. Crit. Care Med. 201, 116–119. doi: 10.1164/rccm.201906-1195le

PubMed Abstract | CrossRef Full Text | Google Scholar

Rudski, L. G., Lai, W. W., Afilalo, J., Hua, L., Handschumacher, M. D., Chandrasekaran, K., et al. (2010). Guidelines for the echocardiographic assessment of the right heart in adults: a report from the American society of echocardiography endorsed by the european association of echocardiography, a registered branch of the european society of cardiology, and the canadian society of echocardiography. J. Am. Soc. Echocardiogr. 23, 685–713. doi: 10.1016/j.echo.2010.05.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Sareban, M., Perz, T., Macholz, F., Reich, B., Schmidt, P., Fried, S., et al. (2020). Preserved right ventricular function but increased right atrial contractile demand in altitude-induced pulmonary hypertension. Int. J. Cardiovasc. Imaging 36, 1069–1076. doi: 10.1007/s10554-020-01803-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Silverman, H. S., Wei, S., Haigney, M. C., Ocampo, C. J., and Stern, M. D. (1997). Myocyte adaptation to chronic hypoxia and development of tolerance to subsequent acute severe hypoxia. Circ. Res. 80, 699–707. doi: 10.1161/01.res.80.5.699

CrossRef Full Text | Google Scholar

Sommer, N., Dietrich, A., Schermuly, R. T., Ghofrani, H. A., Gudermann, T., Schulz, R., et al. (2008). Regulation of hypoxic pulmonary vasoconstriction: basic mechanisms. Eur. Respir. J. 32, 1639–1651.

Google Scholar

Stembridge, M., Ainslie, P. N., and Shave, R. (2016). Mechanisms underlying reductions in stroke volume at rest and during exercise at high altitude. Eur. J. Sport Sci. 16, 577–584. doi: 10.1080/17461391.2015.1071876

PubMed Abstract | CrossRef Full Text | Google Scholar

Takahama, H., Mccully, R. B., Frantz, R. P., and Kane, G. C. (2017). Unraveling the RV ejection doppler envelope: insight into pulmonary artery hemodynamics and disease severity. JACC Cardiovasc. Imaging 10, 1268–1277. doi: 10.1016/j.jcmg.2016.12.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Tamulenaite, E., Zvirblyte, R., Ereminiene, R., Ziginskiene, E., and Ereminiene, E. (2018). Changes of left and right ventricle mechanics and function in patients with end-stage renal disease undergoing haemodialysis. Medicina 54:87. doi: 10.3390/medicina54050087

PubMed Abstract | CrossRef Full Text | Google Scholar

Tello, K., Axmann, J., Ghofrani, H. A., Naeije, R., Narcin, N., Rieth, A., et al. (2018). Relevance of the TAPSE/PASP ratio in pulmonary arterial hypertension. Int. J. Cardiol. 266, 229–235. doi: 10.1016/j.ijcard.2018.01.053

PubMed Abstract | CrossRef Full Text | Google Scholar

Tello, K., Wan, J., Dalmer, A., Vanderpool, R., Ghofrani, H. A., Naeije, R., et al. (2019). Validation of the tricuspid annular plane systolic excursion/systolic pulmonary artery pressure ratio for the assessment of right ventricular-arterial coupling in severe pulmonary hypertension. Circ. Cardiovasc. Imaging 12:e009047.

Google Scholar

Tong, X., Poon, J., Li, A., Kit, C., Yamada, A., Shiino, K., et al. (2018). Validation of cardiac magnetic resonance tissue tracking in the rapid assessment of RV function: a comparative study to echocardiography. Clin. Radiol. 73, 324 e329–324 e318.

Google Scholar

Torbicki, A., Kurzyna, M., Ciurzynski, M., Pruszczyk, P., Pacho, R., Kuch-Wocial, A., et al. (1999). Proximal pulmonary emboli modify right ventricular ejection pattern. Eur. Respir. J. 13, 616–621. doi: 10.1183/09031936.99.13361699

PubMed Abstract | CrossRef Full Text | Google Scholar

Unlu, S., Pagourelias, E., Sezenoz, B., Sahinarslan, A., Uludag, O. M., Gokalp, G., et al. (2019). Higher ultrafiltration rate is associated with right ventricular mechanical dispersion. Anatol. J. Cardiol. 21, 206–213.

Google Scholar

van de Veerdonk, M. C., Kind, T., Marcus, J. T., Mauritz, G. J., Heymans, M. W., Bogaard, H. J., et al. (2011). Progressive right ventricular dysfunction in patients with pulmonary arterial hypertension responding to therapy. J. Am. Coll. Cardiol. 58, 2511–2519. doi: 10.1016/j.jacc.2011.06.068

PubMed Abstract | CrossRef Full Text | Google Scholar

Vonk Noordegraaf, A., Westerhof, B. E., and Westerhof, N. (2017). The relationship between the right ventricle and its load in pulmonary hypertension. J. Am. Coll. Cardiol. 69, 236–243. doi: 10.1016/j.jacc.2016.10.047

PubMed Abstract | CrossRef Full Text | Google Scholar

Williams, A. M., Ainslie, P. N., Anholm, J. D., Gasho, C., Subedi, P., and Stembridge, M. (2019). Left ventricular twist is augmented in hypoxia by beta1-adrenergic-dependent and beta1-adrenergic-independent factors, without evidence of endocardial dysfunction. Circ. Cardiovasc. Imaging 12:e008455.

Google Scholar

Yang, Y., Liu, C., Tian, J., Ding, X., Yu, S., Bian, S., et al. (2020). Preliminary study of right ventricular dyssynchrony under high-altitude exposure: determinants and impacts. Front. Physiol. 11:703.

Google Scholar

Yock, P. G., and Popp, R. L. (1984). Noninvasive estimation of right ventricular systolic pressure by doppler ultrasound in patients with tricuspid regurgitation. Circulation 70, 657–662. doi: 10.1161/01.cir.70.4.657

CrossRef Full Text | Google Scholar

Zhao, H., Kang, Y., Pickle, J., Wang, J., and Han, Y. (2019). Tricuspid annular plane systolic excursion is dependent on right ventricular volume in addition to function. Echocardiography 36, 1459–1466. doi: 10.1111/echo.14439

PubMed Abstract | CrossRef Full Text | Google Scholar

Zoroufian, A., Razmi, T., Taghavi-Shavazi, M., Lotfi-Tokaldany, M., and Jalali, A. (2014). Evaluation of subclinical left ventricular dysfunction in diabetic patients: longitudinal strain velocities and left ventricular dyssynchrony by two-dimensional speckle tracking echocardiography study. Echocardiography 31, 456–463. doi: 10.1111/echo.12389

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: high altitude, right ventricular function, right ventricular outflow tract, speckle tracking echocardiography, tissue Doppler imaging

Citation: Yuan F, Liu C, Yu S, Bian S, Yang J, Ding X, Zhang J, Tan H, Ke J, Yang Y, He C, Zhang C, Rao R, Liu Z, Yang J and Huang L (2021) The Association Between Notching of the Right Ventricular Outflow Tract Flow Velocity Doppler Envelope and Impaired Right Ventricular Function After Acute High-Altitude Exposure. Front. Physiol. 12:639761. doi: 10.3389/fphys.2021.639761

Received: 09 December 2020; Accepted: 26 February 2021;
Published: 01 April 2021.

Edited by:

Nicolas Voituron, Université Paris 13, France

Reviewed by:

Abdul Hameed, Sheffield Teaching Hospital, United Kingdom
David Cristóbal Andrade, University of Antofagasta, Chile

Copyright © 2021 Yuan, Liu, Yu, Bian, Yang, Ding, Zhang, Tan, Ke, Yang, He, Zhang, Rao, Liu, Yang and Huang. 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: Lan Huang, huanglan260@126.com

These authors have contributed equally to this work

Disclaimer: 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.