key: cord-0013023-7r2rk9ri authors: Lin, Jia-Yun; Zhang, Chi-Hao; Zheng, Lei; Song, Chen-Lu; Deng, Wen-Sheng; Zhu, Yi-Ming; Zheng, Li; Wu, Li-Zhong; Sun, Long-Ci; Luo, Meng title: Assessment of a biofluid mechanics-based model for calculating portal pressure in canines date: 2020-08-26 journal: BMC Vet Res DOI: 10.1186/s12917-020-02478-1 sha: 6505e9a1e922fc446c033b265d80579ab1af931c doc_id: 13023 cord_uid: 7r2rk9ri BACKGROUND: Portal hypertension is a severe complication caused by various chronic liver diseases. The standard methods for detecting portal hypertension (hepatic venous pressure gradient and free portal pressure) are available in only a few hospitals due to their technical difficulty and invasiveness; thus, non-invasive measuring methods are needed. This study aimed to establish and assess a novel model to calculate free portal pressure based on biofluid mechanics. RESULT: Comparison of each dog’s virtual and actual free portal pressure showed that a biofluid mechanics-based model could accurately predict free portal pressure (mean difference: -0.220, 95% CI: − 0.738 to 0.298; upper limit of agreement: 2.24, 95% CI: 1.34 to 3.14; lower limit of agreement: -2.68, 95% CI: − 3.58 to − 1.78; intraclass correlation coefficient: 0.98, 95% CI: 0.96 to 0.99; concordance correlation coefficient: 0.97, 95% CI: 0.93 to 0.99) and had a high AUC (0.984, 95% CI: 0.834 to 1.000), sensitivity (92.3, 95% CI: 64.0 to 99.8), specificity (91.7, 95% CI: 61.5 to 99.8), positive likelihood ratio (11.1, 95% CI: 1.7 to 72.8), and low negative likelihood ratio (0.08, 95% CI: 0.01 to 0.6) for detecting portal hypertension. CONCLUSIONS: Our study suggests that the biofluid mechanics-based model was able to accurately predict free portal pressure and detect portal hypertension in canines. With further research and validation, this model might be applicable for calculating human portal pressure, detecting portal hypertensive patients, and evaluating disease progression and treatment efficacy. Therefore, an accurate and non-invasive portal pressure measuring method is needed and would be useful in the diagnosis and evaluation of portal hypertension. Biofluid mechanics is the study of biological flow mechanisms and the inter-relationships with physiological and pathological processes using the fundamental principles of fluid mechanics [12] . Using biofluid mechanics, cardiologists succeeded in calculating fractional flow reserve, the standard assessment of the haemodynamics of coronary stenoses [13, 14] . This motivated us to establish a noninvasive method for accurately calculating FPP. The aim of this study was to establish and assess a biofluid mechanics-based model for predicting FPP and detecting portal hypertension in canines. To validate canine portal hypertensive models, we focused on canine liver fibrogenesis ( Fig. 1 ), liver function, FPP (Table 1) , as well as blood vessel diameters ( Table 2 ) and blood flow velocity (Table 3 ) of their portal venous system. We found that canines from CCl 4 -treated groups developed liver fibrogenesis, had lower blood flow velocity, wider blood vessel diameters, and higher liver function and FPP values as compared with those from a control group. These results proved the success of our canine portal hypertensive model. We imported canine abdominal CT images into the IQQA-Liver system to produce a three-dimensional portal venous model, in which the portal vein and its main branches, including the left and right portal vein, the splenic vein, the superior mesenteric vein, and the inferior mesenteric vein, were visualized precisely. We used FLU-ENT software to divide the model (one case is shown in Fig. 2a ), mesh the model surfaces into triangular surface grids and create the body meshes accordingly (one case is shown in Fig. 2b ). Blood density and overall viscosity were used as the properties of the blood. The "pressure outlet boundary conditions" module was used for the portal vein and the "velocity inlet boundary conditions" module was used for each inlet and outlet branch. The blood flow velocity at the boundaries of each branch was calculated (formula 1-2). FLUENT software was used to solve the Navier-Stokes equations (formula 3-8), simulate the blood flow within the portal venous system (one case is shown in Fig. 3) , and obtain the virtual free portal pressure (vFPP) (one case is shown in Fig. 4) . To assess the numeric correlation between vFPP and FPP, we performed Bland and Altman's limits of agreement analysis (Fig. 5a) , the intraclass correlation coefficient (Fig. 5b) , and Lin's concordance correlation coefficient ( Fig. 5c ) between FPP and vFPP. For Bland and Altman's limits of agreement analysis, the mean of the difference was − 0.220 (95% CI: − 0.738 to 0.298); the upper limit of agreement was 2.24 (95% CI: 1.34 to 3.14); the lower limit of agreement was − 2.68 (95% CI: − 3.58 to − 1.78). The intraclass correlation coefficient was 0.98 (95% CI: 0.96 to 0.99, P < 0.0001). The concordance correlation coefficient was 0.97 (95% CI: 0.93 to 0.99). These results showed that the vFPP model provided a good prediction of FPP. To evaluate the performance of vFPP in the diagnostics of portal hypertension, we carried out receiver operating characteristic (ROC) curve analysis (Fig. 5d ). The area under the ROC curve (AUC) was 0.984 (95% CI: 0.834 to 1.000, P < 0.0001). We selected vFPP = 12 mmHg as the criterion value. The sensitivity was 92.3 (95% CI: 64.0 to 99.8); the specificity was 91.7 (95% CI: 61.5 to 99.8); the positive likelihood ratio was 11.1 (95% CI: 1.7 to 72.8); the negative likelihood ratio was 0.08 (95% CI: 0.01 to 0.6). All canines survived and had no important adverse events during the experiment. Portal hypertension is a life-threatening disease. Portal pressure measurement is of great importance because it is the standard method for detecting portal hypertension. Moreover, portal pressure strongly correlates with severe complications, such as hepatocellular carcinomas, gastrooesophageal variceal haemorrhaging, hepatic encephalopathy, and ascites [15] . In this study, we performed canine experiments to evaluate a biofluid mechanics-based model for calculating vFPP and detecting portal hypertension. Both FPP and HVPG are the standard methods for detecting portal hypertension [1, 2] , however, we used only FPP as a reference because canines' FPP can be easily measured, and balloon-tipped catheters for canine HVPG measurements were unavailable due to the narrowness of the canine vessels. The IQQA system, the FLUENT software, and the Navier-Stokes equations have been widely used to reconstruct blood vessels, simulate blood flow, and precisely calculate the haemodynamics [16] [17] [18] [19] [20] [21] [22] [23] [24] . In the portal venous system, due to the relatively tiny size of blood cells (compared with the size of the vessels) and the steadily fast blood flow, the blood within could be modelled as an incompressible Newtonian fluid; therefore, it is appropriate to apply the IQQA system, the FLUENT software and the Navier-Stokes equations to our study [25, 26] . Based on previous studies, we built the vFPP calculation model using IQQA and FLUENT software to reconstruct canine portal venous system, simulate the blood flow, and calculate the vFPP by the Navier-Stokes equations. After assessing the model by Bland and Altman's limits of agreement analysis, the intraclass correlation coefficient, Lin's concordance correlation coefficient and ROC curve analysis, we showed that this model was able to predict FPP and diagnose portal hypertension accurately. Iranmanesh and Liu also established models for detecting portal hypertensive patients, which showed good diagnostic performance with high sensitivity and specificity [27, 28] . Although their works were able to detect patients suffering from clinically significant portal hypertension (HVPG greater than 10 mmHg) well, their models were not suitable for mild or moderate portal There are some limitations in this study. Firstly, although we demonstrated a high correlation between vFPP and FPP, this result was acquired from canines, not humans. We chose beagles as experimental animals because they have a moderate body size and a digestive system similar to humans. It is also easy to perform examinations on them. In addition, the methods and results of this study might be applied to humans more appropriately. After this study, we applied this model to several portal hypertensive patients who underwent portosystemic shunts or splenectomy with periesophagogastric devascularisation and found similar results (unpublished observations). We are currently carrying out a prospective, randomised, non-controlled, multicentre trial (trial registration number: NCT03470389) to further validate this model in humans [29] . Secondly, the intra-abdominal pressure might have changed after The non-invasive and biofluid mechanics-based model could accurately predict FPP and had high sensitivity and specificity for detecting portal hypertension in canines. With further research and validation, this model might be applicable for calculating human FPP, detecting patients with portal hypertension, and evaluating disease progressions and treatment efficacies. The canines used in this study were purchased from the laboratory animal department of Shanghai Jiagan Biotechnology Co., Ltd. as experimental animals. The study protocol was reviewed and approved by the Animal Care and Use Committee of Shanghai Ninth People's Hospital, School of Medicine, Shanghai Jiao Tong University. All procedures were conducted according to the Animal Experimentation Guidelines of Shanghai Jiao Tong University. The study was performed on 25 adult male beagles (10.5 to 12.5 kg), which were caged with constant temperature (25°C), humidity (60 ± 10%), and circadianrhythmic lighting in the laboratory animal department of Shanghai Jiagan Biotechnology Co., Ltd. To reflect different stages of cirrhosis and portal hypertension, canine portal hypertensive models were induced by continuous subcutaneous injection of CCl 4 . The 25 dogs were divided randomly into five 5-member groups: Control, CCl 4 (1 month), CCl 4 (2 months), CCl 4 (3 months), and CCl 4 (4 months). Randomisation was based on a computergenerated random digits table. Each dog's group was blind to the researchers responsible for histological staining, laboratory tests, Doppler ultrasound, CT, FPP measurement, and vFPP computation in order to prevent biases. The 20 dogs in the CCl 4 -treated groups began receiving CCl 4 administration 1, 2, 3 or 4 months before the end of the study, so that all dogs reached the end of the study simultaneously. CCl 4 was dissolved in olive oil (60% CCl 4 and 40% olive oil) and injected subcutaneously in the Fig. 3 The blood flow simulation of the portal venous system Fig. 4 The blood pressure (Pa) of the portal venous system dorsal area of the canines. This administration was repeated every 10 days at a dose of 1.0 to 1.3 ml/kg. Histological staining, laboratory tests, Doppler ultrasound, CT, FPP measurement, and vFPP computation were performed at the end of the study. All canines survived and continued to live in the laboratory animal department of Shanghai Jiagan Biotechnology Co., Ltd. after the study. Hepatic tissue was taken and fixed in 10% formalin and embedded in paraffin. The paraffin-embedded tissue was sectioned at 5 μm and then placed on slides, deparaffinized in xylene, hydrated in decreasing concentrations of ethanol, and washed in water. After hematoxylin-eosin staining, Masson's trichrome staining, and Sirius red staining, the sections were examined under a microscope. Each dog's peripheral blood samples were taken from the small saphenous vein of the hind leg for blood viscosity and liver function tests. The blood density test was performed by weighing 1 millilitre of blood using an electronic balance. The blood density measurements were repeated at least three times and then averaged. Each canine underwent an abdominal Doppler ultrasound scan after an overnight fast and a venous injection of pentobarbital sodium (30 mg/kg) in accordance with previously published methods [30] . Each canine was in dorsal recumbency throughout the scan. The ultrasound specialists measured the inner diameters and the maximum blood flow velocity of the portal vein and its main branches, including the right branch of the portal vein, the left branch of the portal vein, the portal vein, the splenic vein, the superior mesenteric vein, and the inferior mesenteric vein. The 3-5 MHz Doppler ultrasound probes were used and the insonation angles were between 45°and 65°. Each measurement was repeated twice; both intra-observer variability and inter-observer variability were less than 10%. Each dog was fixed in the supine position and underwent an abdominal contrast-enhanced CT after an overnight fast and a venous injection of pentobarbital sodium (30 mg/kg) in accordance with previously published methods [31, 32] . A non-ionic iodinated contrast agent (600 mg of iodine per kilogram of body weight, 300 mg of iodine per ml, 5 ml per second) was injected. The arterial phase imaging began 10 s after the beginning of the intravenous contrast injection, and portal phase imaging began 30 to 40 s after the beginning of the intravenous contrast injection. Each canine's FPP was measured after general anaesthesia. An abdominal midline incision was made, exposing the right gastroepiploic vein, and a pressure sensor-connected catheter was inserted into the portal vein through the right gastroepiploic vein. The FPP was recorded by a physiological signal acquisition system. The right atrium pressure was defined as the zero reference point. The simulation model of the canine portal venous system, which was created from canine CT images by the IQQA-Liver system version 2.0 (EDDA Technology, Inc., USA), was imported into the Fluent software version 6.3 (ANSYS, Inc., USA). The model was divided into different parts, with each inlet and outlet branch identified as a separate part. The model surface was then meshed into 0.2-1.0 mm triangular surface grids and the body meshes were created accordingly. The laminar viscous model was used. The material type was set to fluid, blood density and overall viscosity were used as the properties of the fluid. The pressure outlet boundary conditions module was used for the portal vein, with the following parameters: backflow reference frame: absolute; gauge pressure: 0; backflow direction specification method: normal to boundary; radial equilibrium pressure distribution: disabled; average pressure specification: disabled; target mass flow rate: disabled. The velocity inlet boundary conditions module was used for each inlet and outlet branch, with the following parameters: velocity specification method: magnitude, normal to boundary; reference frame: absolute; supersonic/initial gauge pressure: 0. The velocity magnitude value was set to the blood flow velocity at the boundaries of each branch, which was calculated according to the inner diameter, blood flow velocity and direction measured by Doppler ultrasound, the inner diameter at the boundaries obtained from the simulation model, and the principle of mass conservation. The equations are as follows: Q: volume flow rate; A: cross-sectional area; v: velocity; d b : inner diameter at the boundaries; v b : velocity at the boundaries; d us : inner diameter measured by Doppler ultrasound; v us : velocity measured by Doppler ultrasound. According to Doppler ultrasound images (one case is shown in Fig. 6 ), the mean blood flow velocity approximately equals 0.7 times the maximum blood flow velocity. Since the blood within the portal venous system can be assumed to be an incompressible Newtonian fluid, blood flow was modelled by the Navier-Stokes equations as follows: The mass conservation equation: ρ: density; t: time; S m : mass added to the continuous phase. For axisymmetric geometries, the mass conservation equation can be written as follows: x: axial coordinate; r: radial coordinate; v x : axial velocity; v r : radial velocity. The momentum conservation equation: p: static pressure; τ : stress tensor; p g ! : gravitational body force; F ! : external body force. τ, the stress tensor, can be written as follows: μ: molecular viscosity; I: unit tensor; T: temperature. The energy conservation equation: E: total energy of fluid; k eff : effective conductive coefficient; h: enthalpy; J: diffusion flux; S h : volumetric heat sources. E, the total energy of fluid, can be written as follows: The FLUENT software was used to solve the Navier-Stokes equations, simulate the blood flow within the portal venous system, and obtain the pressure on each volume grid. The vFPP equals the pressure at the centre of the portal vein. Continuous variables were checked for normal distribution and summarized by either mean ± standard deviation or median and inter-quartile range, as appropriate. Comparison of continuous variables were performed using Student's t-test or analysis of variance for normally distributed variables and the Mann-Whitney U test or the Kruskal-Wallis test for non-normally distributed variables as appropriate. The numeric correlation between FPP and vFPP was analysed by Bland and Altman's limits of agreement analysis [33] . Bias was defined as the mean of the difference between FPP and vFPP. Upper and lower limits of agreement were defined as mean difference ± 1.96 standard deviation of the difference. The numeric correlation between FPP and vFPP was also analysed by the intraclass correlation coefficient and Lin's concordance correlation coefficient. The diagnostic accuracy of vFPP for detecting portal hypertension (FPP greater than 12 mmHg) was assessed by ROC curve analysis, the sensitivity, the specificity, and the likelihood ratio. All tests of significance were at the 5% significance level. Analyses were performed using SPSS version 24.0 (IBM, USA) and MedCalc Statistical Software version 18.11 (MedCalc Software bvba, Belgium). The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request. The canines used in this study were purchased from the laboratory animal department of Shanghai Jiagan Biotechnology Co., Ltd. as experimental animals. The study protocol was reviewed and approved by the Animal Care and Use Committee of Shanghai Ninth People's Hospital, School of Medicine, Shanghai Jiao Tong University. All procedures were conducted according to the Animal Experimentation Guidelines of Shanghai Jiao Tong University. Fig. 6 The Doppler ultrasound image of the portal venous system Expanding consensus in portal hypertension: report of the Baveno VI consensus workshop: stratifying risk and individualizing care for portal hypertension Liver cirrhosis Science to practice: noninvasive assessment of portal hypertension--can US aid in the prediction of portal pressure and monitoring of therapy? The clinical use of HVPG measurements in chronic liver disease Invasive and noninvasive methods to diagnose portal hypertension and esophageal varices Duplex Doppler ultrasound examination of the portal venous system: an emerging novel technique for the estimation of portal vein pressure EUS for portal hypertension: a comprehensive and critical appraisal of clinical and experimental indications Emerging non-invasive approaches for diagnosis and monitoring of portal hypertension A novel noninvasive assessment of hepatic venous pressure gradient and portal pressure computed from computed tomography angiography Virtual portal pressure from anatomic CT angiography Virtual portal pressure gradient from anatomic CT angiography Biofluid mechanics: innovations and challenges Computational fluid dynamics applied to cardiac computed tomography for noninvasive quantification of fractional flow reserve: scientific basis Diagnosis of ischemia-causing coronary stenoses by noninvasive fractional flow reserve computed from coronary computed tomographic angiograms. Results from the prospective multicenter DISCOVER-FLOW (diagnosis of ischemia-causing Stenoses obtained via noninvasive fractional flow reserve) study guidelines on the management of variceal haemorrhage in cirrhotic patients Three-dimensional quantitative evaluation of the segmental functional Reserve in the Cirrhotic Liver Using Multi-Modality Imaging Clinical implications of biliary confluence pattern for bismuth-Corlette type IV Hilar Cholangiocarcinoma applied to Hemihepatectomy Application of a three-dimensional reconstruction technique in liver autotransplantation for end-stage hepatic alveolar Echinococcosis Application of 3D reconstruction for surgical treatment of hepatic alveolar echinococcosis The impact of MRI-based inflow for the hemodynamic evaluation of aortic coarctation Is MRI-based CFD able to improve clinical treatment of coarctations of aorta? Real-time intravascular shear stress in the rabbit abdominal aorta Functional assessment of cerebral artery stenosis: a pilot study based on computational fluid dynamics Blood flow in the cerebral venous system: modeling and simulation Considerations of blood properties, outlet boundary conditions and energy loss approaches in computational fluid dynamics modeling Regularization of blood motion fields by modified Navier-stokes equations Accurate computed tomography-based portal pressure assessment in patients with hepatocellular carcinoma Development and validation of a radiomics signature for clinically significant portal hypertension in cirrhosis (CHESS1701): a prospective multicenter study Establishment and assessment of the hepatic venous pressure gradient using biofluid mechanics (HVPGBFM): protocol for a prospective, randomised, non-controlled, multicentre study Normal duplex Doppler waveforms of major abdominal blood vessels in dogs: a review Computed tomography in small animals--basic principles and state of the art applications Dual-phase CT angiography of the normal canine portal and hepatic vasculature Statistical methods for assessing agreement between two methods of clinical measurement Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations The authors thank Xiao-Yu Yang for his assistance in biofluid mechanics.Authors' contributions ML and L-CS conceived the study. J-YL and LZ1 were responsible for FPP measurement. C-HZ was responsible for vFPP computation. C-LS was responsible for statistical analysis. W-SD was responsible for histological staining. Y-MZ was responsible for laboratory test. LZ2 was responsible for Doppler ultrasound. L-ZW was responsible for CT. J-YL, C-HZ, LZ1 and C-LS drafted the protocol manuscript. All authors contributed to the design of the study and critically revised the successive versions of the manuscript and approved the final version. Not applicable. The authors declare that they have no competing interests.