A protocol for the assessment of diastolic function using Doppler ultrasonography in a preclinical model of pulmonary hypertension due to left heart disease is presented.
Method Article
A protocol for the assessment of diastolic function using Doppler ultrasonography in a preclinical model of pulmonary hypertension due to left heart disease is presented.
Diastolic dysfunction and heart failure with preserved ejection fraction (HFpEF) are significant contributors to pulmonary hypertension (PH) due to left heart disease. Standard in vivo evaluation of relaxation abnormalities in preclinical models includes two-dimensional (2D) echocardiography with Doppler assessment of transmitral flow and tissue relaxation, which is complex, expensive, and requires sophisticated echocardiographic equipment. Here, a valuable surrogate method to evaluate diastolic dysfunction in PH due to left heart disease is demonstrated in a well-established rodent model using pulsed wave Doppler ultrasound without 2D echocardiography. Diastolic transmitral flow pattern is identified from the apical window and is correlated with a simultaneously recorded ECG tracing. Indices of diastolic function are collected, and the identification of these indices is described in this protocol. This low-cost, readily implemented technique identifies relevant markers of diastolic dysfunction associated with PH, cardiovascular disease, and HFpEF. Moreover, it reproducibly detects its progression over time, when 2D echocardiography is not preferred or available. The primary limitations of this method are related to the absence of visualization of cardiac structures, the potential for inaccuracies due to the angle dependence of the Doppler signal, as well as the animal's body habitus, which can be mitigated by adjusting the location of the probe and optimizing its alignment with blood flow.
Pulmonary hypertension (PH) is a deadly disease without a cure that is characterized by increased pulmonary blood pressure1. Left heart disease (LHD) is by far the most common cause of PH worldwide and is associated with high cost and significant morbidity2,3. Cardiac culprits of PH-LHD include multiple different entities that result in either systolic or diastolic dysfunction, such as heart failure, which is a multifaceted clinical syndrome that continues to be among the leading causes of mortality worldwide4. Among patients with heart failure, HFpEF accounts for over half of the incident-based hospital admissions4,5. Diastolic dysfunction has been identified as one of the most important precursors of, and a fundamental feature in the interplay between HFpEF and PH due to LHD (PH-LHD)6,7. Due to the lack of suitable preclinical models that recapitulate most clinical and subclinical features of PH and HFpEF, the understanding of the mechanisms underlying these diseases is incomplete7.
Diastolic dysfunction is a hallmark of PH-LHD and can be assessed by 2D and Doppler echocardiography8. Two-dimensional and Doppler methods for the assessment of left ventricular diastolic function include among other parameters the peak early diastolic E-wave velocity, peak late diastolic A-wave, mitral valve A duration (time interval from A-wave onset to end of A wave at zero baseline), mitral valve E/A ratio, E-wave deceleration time (DT, time interval from peak E-wave extrapolated to the zero-velocity baseline), tissue Doppler imaging (TDI) of the septal mitral valve annulus (e', echocardiographic parameter of tissue relaxation of the mitral valve annulus), and IVRT (isovolumic relaxation time, time between the closure of the aortic valve and the opening of the mitral valve). Echocardiographic markers of impaired relaxation include a combination of reduced or reversed E/A ratio in the setting of A-wave augmentation and E-wave reduction, prolonged IVRT, reduced E-wave deceleration time, delayed tissue Doppler indices (prolonged e' velocity, increased average septal E/e'), and a lower absolute value of global longitudinal strain (GLS)9. Studies comparing echocardiographic and invasive hemodynamic measurements have demonstrated that echocardiographic measurements and invasive measurements are comparable and are a reliable method for evaluating cardiac hemodynamics10. However, 2D echocardiography is complex, expensive, and requires sophisticated echocardiographic equipment and personnel. These factors are a major limitation in following the progression of the disease in animals. Consequently, investigators either rely on two measures (start and end of study) for longitudinal studies or single measurements in a cross-sectional study design.
Amongst these parameters of diastolic dysfunction, E, A, E/A, A duration, DT, and IVRT can be assessed by pulsed wave Doppler echocardiography, obviating the need for 2D echocardiography. In this article, a protocol for the Doppler assessment of diastolic function is demonstrated in the well-established ZSF-1 rat model of PH-LHD11. A total of 25 (10 genetically obese and 15 genetically lean) male ZSF1 rats were longitudinally followed and evaluated sonographically for the development of diastolic dysfunction for sixteen weeks between 8 and 24 weeks of age. The same echocardiographic parameters were evaluated across the two cohorts. All animals were evaluated at baseline (8 weeks of age). Doppler evidence of diastolic dysfunction was manifested by 16 weeks of age, when the interim measurements were taken. End-point measurements were taken at 24 weeks of age. While the investigator performing data collection was unable to be blinded due to the obvious weight gain in the obese cohort, the investigator performing data analysis was blinded to the animals' phenotype. While this readily implementable method offers a route for longitudinal assessment of disease progression in rodents, limitations of this technique include the angle dependence of a good-quality Doppler signal, which may be challenging to obtain depending on the probe alignment and animal size.
All procedures were conducted in accordance with the Johns Hopkins University Animal Care and Use Committee (ACUC) approval for protocol RA21M273. No animals were exposed to discomfort for data collection or the purposes of the measurements. The reagents and the equipment used are listed in the Table of Materials.
1. Preparation of the equipment and the animal
2. Obtaining an ECG tracing
3. Setup of computer
4. Performing Doppler measurements
5. Obtaining data from Doppler ultrasonography
6. Ending the experiment
This method allows longitudinal assessment of diastolic function in PH-LHD. For the most accurate assessment of diastolic function, antegrade and retrograde flow patterns related to the opening and closing of the mitral and aortic valve need to be visualized within the same cardiac cycle over multiple cycles (Figure 1). The closure of the aortic valve marks the beginning of diastole, the earliest phase of which is the isovolumic relaxation, preceding the opening of the mitral valve and the subsequent commencement of the biphasic transmitral flow during early and late diastole. Diastole ends with the cessation of blood flow and the closure of the mitral valve. These flow pattern changes are shown in Figure 1 where the systolic transaortic flow is depicted below the zero-flow baseline, isovolumic relaxation is represented by zero flow between the closure of the aortic valve and the opening of the mitral valve, the biphasic diastolic transmitral flow is depicted as an M-shaped pattern above the baseline, followed by zero flow between the closure of the mitral valve and the opening of the aortic valve. In the absence of the simultaneous detection of the aortic or transmitral flow by spectral Doppler, isovolumic contraction and relaxation times cannot be determined.
In addition, there are multiple indices that can be determined using this model from the biphasic flow pattern (Figure 2), such as: E (normal value for rats: 72.5 cm/s ± 13.9 cm/s), A (normal value for rats: 46.9 cm/s ±15.4 cm/s), E/A ratio (normal value for rats: 1. 7 ± 0.6), and E deceleration time (DT) indicating the deceleration time of the early peak filling velocity (normal value for rats: 48.8 ms ± 13.9 ms)12. However, while these parameters are frequently utilized in humans, they have proven unreliable in rodent models13,14. In the context of this article the IVRT (normal value for rats: 25.5 ±6.3 ms), which is the time interval between the end of systole marked by the cessation of aortic flow, and the beginning of diastole, marked by the commencement of the biphasic transmitral flow, is of special interest for PH-LHD as an indicator for diastolic dysfunction in rodents and has been proven to be reliable for diastolic dysfunction in these models13,14. This diastolic dysfunction is indicated by an initially reversed, then normalized and increased E/A ratio, shortened DT, and prolonged IVRT15. Here, IVRT is a valuable descriptor of left ventricular diastolic function and reflects the speed of relaxation during the earliest, isovolumic phase of diastole. IVRT prolongation suggests an impaired, slower rate of relaxation. This has been found to be the earliest echocardiographic manifestation of PH-LHD in a HFpEF model16 (Figure 3). Additionally, end-point measurements, both 2D and Doppler echo on each animal at the conclusion of the experiment, were performed, and it was found that the measurements obtained were comparable (2D echo 22.64 ms ± 0.90 ms vs. Doppler 23.61 ms ± 1.29 ms). These IVRTs are also consistent with the data reported in the literature, obtained via 2D echo12.
In this model (ZSF1 rats), there is a prolongation of the IVRT over a period of multiple weeks, which can be detected as early as 16 weeks of age, and that continues to increase in obese animals over time, as demonstrated by using Doppler echocardiography only (Figure 3).

Figure 1: Representative Doppler image. (A) This figure shows a representative tracing of the Doppler signal measuring blood flow across the mitral valve (top) and the R wave of the ECG (below). (B) This figure shows a representative tracing of the Doppler signal measuring blood flow across the tricuspid valve (top) and the R wave of the ECG (below). Abbreviations: ES: early start; EPV: early peak velocity; EEAS: early end, atrial start; APV: atrial peak velocity; AE: atrial end; IVCTE: isovolumic contraction time end; IVRT: isovolumic relaxation time start. Please click here to view a larger version of this figure.

Figure 2: Doppler Indices for Mitral valve flow. Depicted here are multiple indices (mean +/- standard deviation) that can be determined from the biphasic flow pattern across the mitral valve in young (8-week-old) male healthy (lean ZSF1) rats: E (n = 14, 71.71 cm/s ± 27.14 cm/s), A (n = 14, 51.25 cm/s ± 24.69 cm/s), E/A ratio (n = 13, 1.364 ± 0.33), and E deceleration time (DT, n = 10, 18.35 ms ± 6.69 ms). Please click here to view a larger version of this figure.

Figure 3: Longitudinal determination of IVRT. This figure shows that the isovolumic relaxation time (IVRT) increases in the ZSF1 rat model of PH-LHD during aging. The values were obtained at the following timepoints: baseline (8 weeks, lean: 21.30 ms ± 3.12 ms, obese:23.61 ms ± 3.33 ms), interim (16 weeks, lean: 26.59 ms ± 5.78 ms, obese: 31.77 ms ± 5.99 ms), and final (24 weeks of age, lean: 23.59 ms ± 3.87 ms, obese: 35.52 ms ± 5.99 ms). Blue represents lean control male rats (n = 10, 9, and 9 at baseline, interim, and final measurements), and red represents ZSF1 obese male rats (n =15, 15, and 10 at baseline, interim, and final measurements) that develop PH-LHD by 24 weeks of age. P<0.0001 by two-way ANOVA. Please click here to view a larger version of this figure.
In this article, a surrogate method to evaluate diastolic dysfunction in PH-LHD is demonstrated in a well-established rodent model using pulsed wave Doppler ultrasound without 2D echocardiography. Diastolic transmitral flow pattern is identified from the apical window and is correlated with a simultaneously recorded ECG tracing, which allows for the determination of indices of diastolic function. This low-cost technique identifies markers of diastolic function and dysfunction longitudinally, avoiding the need for costly 2D echocardiography.
Diastolic dysfunction is a key feature of PH-LHD and is characterized by increased passive ventricular stiffness and impaired active relaxation. Decreased diastolic compliance of the left ventricle leads to increased left ventricular end-diastolic pressure (LVEDP) and pulmonary capillary wedge pressure (PCWP) at rest (≥15 mmHg) and during exercise (≥25 mmHg), which, while diagnostic of HFpEF and a precursor for PH-LHD, are challenging to determine without invasive techniques or advanced imaging modalities17,18. The current gold standard for obtaining these values is cardiac catheterization; however, this is not always feasible (e.g., longitudinal measurements) or affordable (high cost of surgery or core services) in preclinical models of PH-LHD19. Given the difficulty in obtaining those values directly, it is important to use additional markers of diastolic dysfunction that can be obtained non-invasively, like the peak rate of early filling (E), peak rate filling rate during atrial systole (A), E/A ratio, E deceleration time, and IVRT11. This can be accomplished using echocardiographic modalities, including image-guided pulse wave Doppler echocardiography, which is considered the standard noninvasive method to evaluate diastolic function both in the clinical and in the preclinical setting; however, this too is an expensive technique, especially in rodents9. Diastolic transmitral flow characteristics and their temporal relation to tissue relaxation and cardiac events can be assessed from the apical 4-chamber view using 2D echocardiography, by aligning the ultrasound beam with the blood flow, and placing the sample volume at the level of leaflet coaptation20. A major limitation of this technique is the fact that E, E/A, and IVRT are load and heart rate-dependent indices of diastolic function, hence the inherent difficulty in controlling for loading conditions. Indeed, in these studies, only heart rate control was attempted by titrating inhaled anesthetic concentration to a heart rate of 300-350 beats/min across the cohorts. Therefore, the potential for bias could not be fully eliminated13.
Using Doppler to evaluate transmitral blood flow has several advantages. First, it is a noninvasive technique that can be done in sedated animals at various timepoints without exposing the animals to surgical manipulation and its associated complications. It is suitable to assess velocities across valves and other cardiovascular structures at a precise location. The cost of the setup is only a fraction of the equipment and upkeep cost of a 2D echocardiography setup, and it can be done without core facility support or specialized training in advanced 2D echocardiography.
Limitations
The downsides of this technique are first and foremost that it is a non-image guided technique, and it may be difficult to align the ultrasound beam with blood flow across the mitral valve, with a risk of assessing tricuspid valve velocities instead of mitral valve flows21. However, the tricuspid valve annulus is larger than the mitral valve annulus; therefore, even if the pressure is equal across both valves, the flow pattern across the tricuspid valve would still be slightly less than across the mitral valve. Though with additional valve disease, a rare situation of these velocities being similar may occur, however, in the vast majority of models that focus on a single disease, it is unlikely to be identical. Also, this technique requires a dedicated operator who has practiced this on many animals to be facile with the fine motor adjustments required to obtain a reliable and valid signal. Lastly, this technique does not allow for visual inspection of other cardiac structures or function that might be relevant for a certain pathology. It is therefore most useful for investigators with training in cardiovascular science and ultrasound imaging, who would need to longitudinally follow diastolic heart function in a large number of animals, while minimizing cost. Utilizing this technique assumes an understanding of cardiac physiology and rodent anatomy.
Critical steps
To ensure good signal quality and valid measurements, proper animal positioning, probe selection and positioning, and selection of the optimal signal sampling range are required. Among the commonly encountered challenges is the difficulty in obtaining a Doppler signal. Here, it is recommended that the probe frequency be matched to the animal's body, and the ultrasound beam aligned with the blood flow to the extent possible. This is best accomplished with the ultrasound probe tip placed subcostally, just beside the xyphoid, to the left of the midline and pointed towards the left ear of the animal (see step 4.2). Signal quality is considered adequate when the spectrogram is crisp and in line with the expected flow velocities (Figure 1A). This can be improved by fine adjustment of the range of signal sampling (see step 4.4). An adequate diastolic flow pattern below the baseline warrants the setting of the "Direction" switch to "TWD" (toward). Blunted, low biphasic velocities above the baseline raise suspicion for detecting trans-tricuspidal flow (Figure 1B). Here, a slight leftward movement of the probe tip is recommended, without changing its angle in relation to the skin (see step 4.2).
An important research application of this technique is the noninvasive longitudinal assessment of diastolic function in rodent preclinical models of PH and HFpEF to identify a cohesive time course in this disease, in particular in females pre-, peri-, and post-menopause.
In summary, this low-cost technique identifies relevant markers of diastolic dysfunction associated with PH-LHD, and reproducibly detects its progression over time, when 2D echocardiography is not preferred or available. The primary limitation of this method is related to the absence of visualization of cardiac structures and the potential for inaccuracies due to the angle dependence of the Doppler signal, as well as the animal's body.
The authors have no financial relationships with any manufacturers of the materials used.
This work was supported by the Career Development Award of the American Heart Association to MB (24CDA1267633), an NHLBI R56 award to JS (1R56HL169285), and an NHLBI award to LS (R01HL14811201).
| Name | Company | Catalog Number | Comments |
|---|---|---|---|
| Anesthesia setup | VetEquip Inc. | 901806 | |
| Computer | Dell | 7420 Plus | |
| Cotton swab | Puritan Medical Products Company | 806-WC | |
| Depilator cream | Church and Dwight | 300725 | |
| Doppler probe 10 and 20MHz | INDUS Instruments | NA | |
| Doppler setup | INDUS Instruments | NA | |
| ECG pad | INDUS Instruments | NA | |
| Electrode gel | Parker Laboratories | 45993 | |
| Isoflurane vaporizer | VetEquip Inc. | 911103 | |
| Ultrasound gel | Parker Laboratories | 45659 | |
| Vented anesthetizing induction chamber | VetEquip Inc. | 942102 |