Abstract

The jugular venous pulse (JVP) provides valuable information about cardiac haemodynamics and filling pressures and is an indirect estimate of the central venous pressure (CVP). Recently it has been proven that JVP can be obtained by measuring the cross-sectional area (CSA) of the IJV on each sonogram of an ultrasound B-mode sonogram sequence. It has also been proven that during its pulsation the IJV is distended and hence that the pressure gradient drives the IJV haemodynamics. If this is true, then it will imply the following: (i) the blood velocity in the IJV is a periodic function of the time with period equal to the cardiac period and (ii) the instantaneous blood velocity is given by a time function that can be derived from a flow-dynamics theory that uses the instantaneous pressure gradient as a parameter. The aim of the present study is to confirm the hypothesis that JVP regulates the IJV blood flow and that pressure waves are transmitted from the heart toward the brain through the IJV wall.

1. Background

The evaluation of the jugular venous pulse (JVP), defined as the movement of expansion of the jugular vein due to changes in pressure in the right atrium, provides valuable information about cardiac haemodynamics and filling pressures [1], characteristic wave patterns pathognomic of cardiac diseases [2], and an indirect estimate of the central venous pressure (CVP). The JVP evaluation can be useful in the diagnosis and/or prognosis of many heart diseases [3]. Such a pulse consists of three positive waves (, , and ) and two descents, defined, respectively, as and . Wave corresponds to the atrial contraction and is synchronized with the wave of the electrocardiogram (ECG). Descent corresponds to the lowering of the atrioventricular septum, interrupted by a small positive wave in relation to the closure of the tricuspid valve; the third wave corresponds to the cardiac systole and is followed by the descent which corresponds to the opening of the tricuspid valve.

In a recent paper [4] we proved that JVP can be obtained by a simple ultrasound (US) B-mode investigation that consists in measuring the cross-sectional area (CSA) of the IJV on each sonogram of video clip acquired in the transversal plane. In that paper, we acquired the time-dependent CSA datasets of three healthy subjects and then calculated the autocorrelation function of the datasets to show that they were periodic and finally we showed that their wave form presented the same , , and waves as the JVP. Our study gave a quantification of the IJV CSA variation during the cardiac cycle. We also have seen that the IJV perimeter, measured on each sonogram of the video clip, was correlating with the IJV CSA (). On this point, we believe that it is desirable to have a confirmation of this finding also using a different imaging modality. However, a direct explanation of this finding is that, when in supine position, the pulsation of the IJV is a distension of its wall. This result is very important because it means that there is a time varying transmural pressure whose effect is clearly visible and cannot be neglected; in fact, large changes in transmural pressure are required to induce CSA deformation accompanied by a stretching of the wall [5, 6]. If this is true then the following points are also true: (i) the blood velocity in the IJV is a periodic function of the time with period equal to the cardiac period and (ii) the instantaneous blood velocity is given by a time function that can be derived from a flow-dynamics theory that uses the instantaneous pressure gradient as a parameter, for example, the linear Womersley solution of the Navier-Stokes equations [7].

It is worth noting that, in an elastic medium as is the IJV wall, an impulse that gives rise to a wall distension will be propagated following the wave equation of d’Alembert. The time periodic variation of the CSA and hence of pressure measured for the IJV is hence related to the propagation of pressure waves propagating from the heart toward the brain. Thus, in supine position, when the IJV is distended, the blood flow is governed by an oscillating pressure gradient [7] and pressure waves generated by the cardiac revolution are then transmitted from the heart toward the brain (see Figure 1).

The aim of the present study is to verify such a claim by comparing the blood velocity assessed by the current US Doppler method with the blood velocity calculated using the linear Womersley equations. The correspondence between them will confirm the hypothesis that JVP regulates the IJV blood flow and that theoretically pressure waves are transmitted from the heart toward the brain through the IJV wall.

2. Methods

2.1. Subjects Scanning and Protocol

We have chosen two subjects, labeled as , (a -year-old female and a -year-old male, resp., with no history of cardiovascular, hepatic, gastrointestinal, renal, and cerebral diseases. ECD screening for CCSVI was completely negative [810]), out of our database of US scan of the neck by using My-Lab 70 x-vision system (ESAOTE, Genoa, Italy) with a linear array probe 7.5–11 MHz and Vivid-q ultrasound system (GE Medical Systems Ultrasound, Horten, Norway) equipped with a linear probe (L11-812 MHz). They were chosen for this study because they were representative of a small and normal IJV CSA, respectively (with respect to the average diameter value of about 1.3 cm reported in [11]).

The assessment of the jugular CSA was performed using a B-mode scan in the transverse plane of the right IJV at c5/c6 level. Such region corresponds to the segment close to the junction of the IJV with the subclavian vein. For each subject we recorded first a transversal video clip of the CSA (i.e., a sonogram sequence) and immediately after a velocity spectral Doppler trace, both for at least four cardiac cycles. The trace was obtained selecting the “mean velocity mode,” which automatically produces a weighted average of the velocity over the whole US reflected spectrum. This results in a green trace (see Figure 3) that corresponds to the blood velocity averaged across the sample volume. The latter is selected to roughly correspond to the IJV lumen.

This study was conducted in accordance with the Ethical Standards of the Committee on Human Experimentation of the Azienda Ospedaliera Universitaria di Ferrara. All the volunteers signed an informed consent form.

2.2. Pulsed Flow-Dynamics Assumptions

Following the original paper of Womersley [7], we denote as the instantaneous blood velocity in direction and averaged over the IJV CSA. The velocity is calculated, at any time , using the Womersley equations, for which the time-dependent pressure gradient is required. The internal pressure was here assumed to be linearly related to the [6, 12, 13] as follows:where is a convenient additive constant and is the IJV compliance for unit of length defined as Given a pressure wave velocity , the wave equation states that

Substituting from (1) into (3), the pressure gradient is obtained. The minus sign on the right side of (3) represents a wave moving toward the positive direction of the -axis while the plus sign represents a wave moving toward the negative . In order to simplify the interpretation of the results, in this paper we chose to assume the positive direction of the from the head toward the heart (see Figure 1). The pressure waves are generated from the heart pulsation and propagate toward the jugular vein in the negative direction; for this reason they are represented by positive sign in (3). In this work we neglect the effects of reflected waves propagating from the brain toward the heart.

Moreover, the Womersley equations only give the instantaneous oscillating part of the blood velocity. The net instantaneous value of the blood velocity is given by summing the steady velocity with the oscillating velocity. The steady component of the velocity can be obtained by the Poiseuille law once the pressure gradient between the head and the heart is known. In this study we are not interested in a full assessment of the mean velocity so only the oscillating velocity is considered and analyzed.

Having assumed a positive from the head to the heart, the net blood velocity is positive (see Figure 1).

2.3. IJV Blood Velocity Assessment: The Oscillating Component

For each subject, the CSA dataset, representing the instantaneous value of the CSA during a cardiac cycle, was derived from the acquired video clip using the semiautomatic algorithm described in [4]. The period of the IJV pulsation was obtained by its discrete Fourier transform calculated using the Grace software [14]. Then we write so that runs from to during one pulse period. The mean CSA, denoted by , over a period is calculated. The Fourier coefficients of pressure gradient up to the tenth harmonic are calculated. These coefficients are then converted to modulus and phase . The blood velocity is then calculated from Equation in [7] by substituting the pressure gradient expressed as a Fourier series. The blood velocity is given by the summation of ten terms as in Equation in [7]:where we have substituted the term in [7] by , which represents the time-dependent vein CSA. The parameter is the blood viscosity (0.04 Pa), that runs from 1 to 10 is the order of the Fourier coefficient, and and are derived from the Bessel function and tabulated in [7]. Finally, is the well-known Womersley number given bywhere is the blood density (1.05 g/mL) and .

2.4. Spectral Doppler Averaged Velocity

Below we describe the procedure adopted to obtain a dataset sampled from the US Doppler velocity trace. Common commercial US systems do not allow the Doppler dataset to be exported; for this reason, in this work, the mean velocity dataset is obtained by digitally processing the image shot of the Doppler trace. The mean blood velocity is represented as a line overlapping with the spectral Doppler trace (Figure 3). The line is composed of pixels. An in-house developed procedure identifies the pixel belonging to the line by its RGB values. The coordinates and of each pixel were automatically recorded by the procedure, where is an index going from 1 to . The function has been obtained from and values using the following procedure:where is the index for the pixels and is the distance in pixels between two time divisions (separated by 1 s) of the time axes while is the distance in pixels between the 0 and 100 cm/s divisions, measured along the -axis. Combining the two expressions in (6) produces

The sampled Doppler velocity dataset is then plotted together with the calculated Womersley averaged velocity (). The instant of time corresponding to the maximum value of velocity was detected for both the Doppler and and it was used as reference marker to overlap the two plots.

2.4.1. Detailed Calculation

The detailed algorithm of all the calculations needed to produce the following results is reported in the appendix.

3. Results

3.1. IJV Blood Velocity Assessment

Each result is produced for both subjects and and presented in the text following such order. Each figure starting from Figure 2 is divided into top and bottom for subjects and , respectively. The main numerical results are summarized in Table 1.

The US-JVP are plotted for both subjects in Figure 2; their repetition periods are 0.9 and 0.92 s for subjects and , respectively; the time averaged cross-sectional areas are 0.25 and 1.25 cm2. Both subjects show clearly visible and waves, while the one is barely detectable.

The time derivatives of the internal pressure and pressure gradient were obtained from the CSA dataset as described in the appendix.

The mean blood velocity was calculated following (4) and is plotted together with the CSA in Figure 4. It is worth noting that for both subjects the blood velocity changes inversely with respect to the IJV CSA; that is, when the IJV CSA goes toward its maximum the blood velocity goes toward its minimum and vice versa.

In Figure 4, is plotted together with the pressure gradient .

3.2. Spectral Doppler Averaged Velocity

Figure 3 depicts the spectral Doppler velocity traces. The time periods of such traces turned out to be 0.95 and 0.88 s for subjects and , respectively (see Table 1). Such values are very close to that of the US-JVP period; this is an expected result since the spectral Doppler trace was acquired as soon as the B-mode IJV investigation was completed so as to avoid changes in physiological conditions.

It is evident that the traces of both subjects and show a “two-wave” profile. Figure 5 depicts the blood velocity sampled from the spectral Doppler versus the calculated ones for subjects and , respectively. For both subjects, sampled and calculated velocity values show the same time dependence or, in other words, the same wave form. A reasonable agreement of their amplitudes was obtained by requiring a CVP up to 10 mmHg (see Appendix A.5 for details).

4. Discussion

4.1. IJV Blood Velocity Assessment

The main result of the present study is the well apparent relationship between the blood velocity in the IJV measured by the current Doppler system and the same parameter derived from its pulse. For both subjects the presence of the two positive peak waves, “” and “” waves, respectively, in the US-JVP (the CSA diagram), generates two descents in the assessed velocity time diagram within the same cardiac period. Although the relationship between the JVP and the variation of pressure in the heart has been known since 1902 [15], until now, this finding has not been taken into consideration for assessing the brain outflow. Currently, the Doppler flow quantification is provided by the product of the time average velocity of flow and the CSA [1619]. This quantitative approach currently provides an accepted objective criterion for assessing the flow rate. However, despite the pulsating variation during the cardiac cycle of the IJV CSA, the flow, according to this method, is calculated by using just a single CSA value. This raises the question as to whether overlooking the pulsatile variation of the IJV CSA may affect the US evaluation of cerebral venous return in the clinical setting. In addition, the results of the present study further corroborate the role of the mechanical propagation of the cardiac contraction in the upward direction. This wave propagation is not haemodynamically negligible because the present study clearly demonstrates how this phenomenon influences blood flow velocity. There are several studies which suggest that flow waves directed up towards the brain can contribute to neurodegenerative diseases by affecting cerebrospinal fluid absorption and/or brain perfusion, such as Alzheimer [8, 9, 20].

Our data clearly demonstrates how the assumption of the Poiseuille theorem in the case of pulsatile flow may expose the flow assessment to error. We focused our study on US but of course the MRV methodology used for IJV flow assessment could also be affected by the same incorrect assumption. We know that the distal segment of the IJV is characterized by significant CSA variation over time, an effect linked with the transmission of pressure waves from the right atrium. This phenomenon, very well known as JVP, is strongly connected with important regulators of the cardiovascular function: degree of filling, heart pulsation, and capacity of emptying of the IJV.

Moreover, the positive waves “,” “,” and “” clearly show a brain directed propagation of pressure waves with significance never investigated, as yet. We know that propagation is greater when the degree of filling of the venous system is higher and/or when the heart function is compromised. Acute symptoms in the latter rapidly appear in the pulmonary function but we have never investigated if this may have significance for cognitive or neurological decline already described in this particular category of patients. We speculate that cardiac, carotid, and jugular signals should be synchronized and used to assess the circulatory axis between the heart and brain, with significant changes in current clinical practice, in neuroimaging as well.

The results reported here can be reproduced using MRV because CSA and blood velocity can be measured at each instant of time as we did.

4.2. Pulsed Flow-Dynamics Assumptions

The Womersley approach presented in [7] derives the blood velocity from the pressure gradient assuming a rigid tube as a model. Such a model accounts for the pressure gradient time dependence but neglects the movement of the tube wall in the longitudinal and lateral direction. In later papers, Womersley made corrections accounting for such variability [13, 21] resulting in a new equation having the same form as the former (see Equation in [13]) and producing a flow about 10% higher in magnitude. In this paper we choose to follow the simpler approach presented in [7] for three reasons: the cited paper is widely known and widely adopted, the presented analytical model can be easily implemented making our results easily reproducible, and such a model results in an easily treatable linear equation and its results also hold for pulsing flow in an elastic tube [13, 21].

However, we encourage the adoption of more sophisticated models when the blood flow has to be calculated for clinical reasons (see, e.g., [12]).

The assumption of simple linear or quadratic relationships between CSA and pressure seems to be correct in the present case. In fact, considering the well-known “tube-law” [6], the CSA/perimeter correlation means that IJV pulsation lies in the curve region (i) where the transmural pressure is greater than the buckling pressure (see Figure in [5]); hence a linear relationship between pressure and CSA is straightforward. On the other hand, the elliptical shape of the IJV during its pulsation could suggest that the transmural pressure is below the buckling value; hence region (ii) could also be affected by the pulsation. Of course this finding needs to be further investigated.

4.3. Implication for Pressure Wave

The oscillating nature of the IJV is potentially relevant in the understanding of the relationship between brain drainage and several neurological disorders [2229]. In particular, the effect of pulsatile flow has been recently hypothesized to be associated with retrograde hypertension transmitted from the IJV which, presumably, underlies transient global amnesia [25] and perhaps other neurological disorders, as recently hypothesized [3032].

Appendix

A. Velocity Calculation Detailed Procedure

A.1. IJV Elasticity Model

In this paragraph we present an expression for the vein volumetric compliance , that is, the ratio of volume variation with respect to pressure variation , and consequently for the compliance per unit of length . We assumethat implies the idea that the vein may change in radius but not in length. From the definition of volumetric compliance it follows thatThe coefficient of functional elasticity is defined asthen, using the definitions given in (A.1), we haveFrom the two definitions of given in (A.2) and (A.4) it follows thatFor cylindrical geometry, the relationship between the Young Modulus and can be assumed to be [33]where is the vein wall thickness. The pressure wave velocity , that is, the velocity of propagation of the JVP in the jugular vein, is calculated following the Moens-Korteweg equation:which for the assumed relationship between and is dependent only on the parameters and :The central venous pressure (CVP), measured close to the IJV, is known to vary up to 5 mmHg during the cardiac cycle. Since during a cardiac cycle is given bythenEquation (A.10) can also be used with (A.5) to determine which is given as

A.2. Parameter Data

The following parameter value has been assumed [7, 33]:(i)Blood viscosity  Pa × s.(ii)Blood density  kg/m3.(iii)Vein wall thickness mm.

A.3. Input Data

From the US transversal scan of the IJV the following data represents the input for the velocity calculation procedure:(i)The US-JVP, consisting in a dataset with running from 1 to , where is the number of sonograms in the US video clip with each element of the dataset representing the IJV area measurement on the th sonogram, in pixels.(ii)The US video clip duration () in seconds.(iii)The US sonogram pixel dimension () in cm.

From the longitudinal US Doppler trace, consider the following:(i)The space averaged velocity time diagram, consisting in a dataset in pixels.(ii)The distance () between two divisions in the time axis of the Doppler trace.

A.4. Cross-Sectional Area

(1) dataset discrete Fourier transform (DFT) is operated (e.g., using Grace software).(2)The sonogram repetition period () of the dataset is obtained by the maximum of the DFT produced at .(3)The sonogram time interval () is calculated as .(4)The dataset time repetition period is then calculated as .(5)The dataset is multiplied by the pixel area () in order for each element of the dataset to represent the IJV area measurement on the th sonogram, in cm2.(6)The differential CSA dataset is obtained as .(7)The discrete partial time derivative dataset is obtained dividing the by .(8)Average CSA is calculated as .(9)The average vein diameter is calculated as .

A.5. Pressure

The following steps focus on the pressure dataset calculation:(1)The compliance for unit of length is calculated as , where is a convenient parameter representing the maximum value reached by the CVP during the cardiac cycle. A typical value is 5 mmHg.(2)The pulse wave velocity (, lower case) is calculated as .(3)Time derivative pressure dataset, in unit of mmHg, is given by .(4)It is then converted into Pa by multiplying the dataset by 133.322 Pa/mmHg.(5)The pressure gradient along the -axis direction is calculated as .(6)The parameter is defined as .(7)The pressure gradient is then defined as .

A.6. Pressure Gradient Fourier Series

The pressure gradient is then decomposed in its Fourier terms as

The series coefficients are calculated as follows:(1).(2).(3).(4).(5).

A.7. Blood Velocity Fourier Series

The blood velocity expression is obtained by the pressure gradient expression as(1)The parameter for each harmonic is calculated as .(2)The coefficients and are tabulated in [7] for each value from 0 to 10. For the asymptotic expression was used [7].

Abbreviations

Acronyms
CCSVI:Chronic venous cerebrospinal insufficiency
CSA:Cross-sectional area
CVP:Central venous pressure
DFT:Discrete Fourier transform
DZP:Pressure gradient
ECD:Echo color Doppler
ECG:Electrocardiogram
IJV:Internal jugular vein
JVP:Jugular venous pulse
MRV:Magnetic resonance venography
:Distance in pixels between two time divisions
:Distance in pixels between the 0 and 100 cm/s divisions
:Pearson coefficient
RGB:Red Green Blue
:Sonogram repetition period
US:Ultrasound.
Symbols
:Womersley number
:Pressure waves velocity
:Compliance
:Compliance for unit of length
:Sonogram time interval
:Coefficient of functional elasticity
:Vein wall thickness
:Blood viscosity
:Young Modulus
:Pixel dimension
:Pressure
:Blood density
:Cardiac period
:Volume
:Blood flow velocity averaged over the IJV CSA.

Conflict of Interests

The authors declare that they have no conflict of interests.

Authors’ Contribution

Francesco Sisini conceived the study, developed the algorithm for the data analysis, analyzed the US studies, performed the data analysis, and wrote the paper from the scratch. Eleuterio Toro wrote the paper and supervised all the mathematical aspects of the research. Mauro Gambaccini wrote the paper and critically revised it. Paolo Zamboni analyzed the jugular diagram, wrote the paper, and supervised all the medical aspects of the research.

Acknowledgments

This study was partially supported by the Italian Ministry of Education, University and Research (MIUR Programme PRIN 2010-2011), Grant no. 2010XE5L2R, and thanks are due to Jennifer Cowd, Xpress Translations Ltd. ([email protected]) for the proofreading of the paper.