A computational study of red blood cell deformability effect on hemodynamic alteration in capillary vessel networks

The effect of the endothelial surface layer on cell-cell interactions in microvessel bifurcations

Red blood cells (RBCs) carry oxygen and make up 40-45% of blood by volume in large vessels down to 10% or less in smaller capillaries. Because of their finite size and large volume fraction, they are heterogeneously distributed throughout the body. This is partially because RBCs are distributed or partitioned nonuniformly at diverging vessel bifurcations where blood flows from one vessel into two. Despite its increased recognition as an important player in the microvasculature, few studies have explored how the endothelial surface layer (ESL; a vessel wall coating) may affect partitioning and RBC dynamics at diverging vessel bifurcations. Here, we use a mathematical and computational model to consider how altering ESL properties, as can occur in pathological scenarios, change RBC partitioning, deformation, and penetration of the ESL. The two-dimensional finite element model considers pairs of cells, represented by interconnected viscoelastic elements, passing through an ESL-lined diverging vessel bifurcation. The properties of the ESL include the hydraulic resistivity and an osmotic pressure difference modeling how easily fluid flows through the ESL and how easily the ESL is structurally compressed, respectively. We find that cell-cell interaction leads to more uniform partitioning and greatly enhances the effects of ESL properties, especially for deformation and penetration. This includes the trend that increased hydraulic resistivity leads to more uniform partitioning, increased deformation, and decreased penetration. It also includes the trend that decreased osmotic pressure increases penetration.

Impact of elevated IOP on lamina cribrosa oxygenation; A combined experimental-computational study on monkeys

Purpose

Our goal is to evaluate how lamina cribrosa (LC) oxygenation is affected by the tissue distortions resulting from elevated IOP.

Design

Experimental study on monkeys.

Subjects

Four healthy monkey eyes with OCT scans with IOP of 10 to 50 mmHg, and then with histological sections of LC.

Methods

Since in-vivo LC oxygenation measurement is not yet possible, we used 3D eye-specific numerical models of the LC vasculature which we subjected to experimentally-derived tissue deformations. We reconstructed 3D models of the LC vessel networks of 4 healthy monkey eyes from histological sections. We also obtained in-vivo IOP-induced tissue deformations from a healthy monkey using OCT images and digital volume correlation analysis techniques. The extent that LC vessels distort under a given OCT-derived tissue strain remains unknown. We biomechanics-based mapping techniques: cross-sectional and isotropic. The hemodynamics and oxygenations of the four vessel networks were simulated for deformations at several IOPs up to 60mmHg. The results were used to determine the effects of IOP on LC oxygen supply, assorting the extent of tissue mild and severe hypoxia.

Main Outcome Measures

IOP-induced deformation, vasculature structure, blood supply, and oxygen supply for LC region.

Result

IOP-induced deformations reduced LC oxygenation significantly. More than 20% of LC tissue suffered from mild hypoxia when IOP reached 30 mmHg. Extreme IOP(>50mmHg) led to large severe hypoxia regions (>30%) in the isotropic mapping cases.

Conclusion

Our models predicted that moderately elevated IOP can lead to mild hypoxia in a substantial part of the LC, which, if sustained chronically, may contribute to neural tissue damage. For extreme IOP elevations, severe hypoxia was predicted, which would potentially cause more immediate damage. Our findings suggest that despite the remarkable LC vascular robustness, IOP-induced distortions can potentially contribute to glaucomatous neuropathy.

In vivo real-time positron emission particle tracking (PEPT) and single particle PET

Positron emission particle tracking (PEPT) enables 3D localization and tracking of single positron-emitting radiolabelled particles with high spatiotemporal resolution. The translation of PEPT to the biomedical imaging field has been limited due to the lack of methods to radiolabel biocompatible particles with sufficient specific activity and protocols to isolate a single particle in the sub-micrometre size range, below the threshold for capillary embolization. Here we report two key developments: the synthesis and 68Ga-radiolabelling of homogeneous silica particles of 950 nm diameter with unprecedented specific activities (2.1 ± 1.4 kBq per particle), and the isolation and manipulation of a single particle. We have combined these developments to perform in vivo PEPT and dynamic positron emission tomography (PET) imaging of a single radiolabelled sub-micrometre size particle using a pre-clinical positron emission tomography/computed tomography scanner. This work opens possibilities for quantitative assessment of haemodynamics in vivo in real time, at the whole-body level using minimal amounts of injected radioactive dose and material.

The interaction between particles and vascular endothelium in blood flow

Particle-based drug delivery systems have shown promising application potential to treat human diseases; however, an incomplete understanding of their interactions with vascular endothelium in blood flow prevents their inclusion into mainstream clinical applications. The flow performance of nano/micro-sized particles in the blood are disturbed by many external/internal factors, including blood constituents, particle properties, and endothelium bioactivities, affecting the fate of particles in vivo and therapeutic effects for diseases. This review highlights how the blood constituents, hemodynamic environment and particle properties influence the interactions and particle activities in vivo. Moreover, we briefly summarized the structure and functions of endothelium and simulated devices for studying particle performance under blood flow conditions. Finally, based on particle-endothelium interactions, we propose future opportunities for novel therapeutic strategies and provide solutions to challenges in particle delivery systems for accelerating their clinical translation. This review helps provoke an increasing in-depth understanding of particle-endothelium interactions and inspires more strategies that may benefit the development of particle medicine.

Predicting capillary vessel network hemodynamics in silico by machine learning

Blood velocity and red blood cell (RBC) distribution profiles in a capillary vessel cross-section in the microcirculation are generally complex and do not follow Poiseuille's parabolic or uniform pattern. Existing imaging techniques used to map large microvascular networks in vivo do not allow a direct measurement of full 3D velocity and RBC concentration profiles, although such information is needed for accurate evaluation of the physiological variables, such as the wall shear stress (WSS) and near-wall cell-free layer (CFL), that play critical roles in blood flow regulation, disease progression, angiogenesis, and hemostasis. Theoretical network flow models, often used for hemodynamic predictions in experimentally acquired images of the microvascular network, cannot provide the full 3D profiles either. In contrast, such information can be readily obtained from high-fidelity computational models that treat blood as a suspension of deformable RBCs. These models, however, are computationally expensive and not feasible for extension to the microvascular network at large spatial scales up to an organ level. To overcome such limitations, here we present machine learning (ML) models that bypass such expensive computations but provide highly accurate and full 3D profiles of the blood velocity, RBC concentration, WSS, and CFL in every vessel in the microvascular network. The ML models, which are based on artificial neural networks and convolution-based U-net models, predict hemodynamic quantities that compare very well against the true data but reduce the prediction time by several orders. This study therefore paves the way for ML to make detailed and accurate hemodynamic predictions in spatially large microvascular networks at an organ-scale.

Photoplethysmography for Assessing Microcirculation in Hypertensive Patients After Taking Antihypertensive Drugs: A Review

High blood pressure (BP) is a common disease and is associated with many chronic diseases. Measuring BP is essential for the treatment and management of many diseases, and therefore there is a growing need for a non-invasive, sleeveless and continuous BP monitoring device. With the development of technology, pulse waveform analysis using photoplethysmography (PPG) has become more feasible for evaluating BP. This study aimed to evaluate the changes of vascular elasticity and blood volume over time by using the characteristic parameters extracted by PPG. We reviewed the latest progress and literature on the observation of capillary network characteristics in hypertensive and non-hypertensive patients by PPG, the influence of different drugs on microcirculation characteristics in hypertensive patients with PPG, and further explored the key relationship between microcirculation and hypertension. We found that the PPG waveform produced by the fingertips of hypertensive patients is very different from that of healthy people, and the PPG waveform changes significantly during diastolic period after antihypertensive treatment. With the rapid development of medical technology, people can get more intuitive microcirculation image data, which provides beneficial help for the comprehensive understanding of hypertension.

Partitioning of dense RBC suspensions in single microfluidic bifurcations: role of cell deformability and bifurcation angle

Red blood cells (RBCs) are a key determinant of human physiology and their behaviour becomes extremely heterogeneous as they navigate in narrow, bifurcating vessels in the microvasculature, affecting local haemodynamics. This is due to partitioning in bifurcations which is dependent on the biomechanical properties of RBCs, especially deformability. We examine the effect of deformability on the haematocrit distributions of dense RBC suspensions flowing in a single, asymmetric Y-shaped bifurcation, experimentally. Human RBC suspensions (healthy and artificially hardened) at 20% haematocrit (Ht) were perfused through the microchannels at different flow ratios between the outlet branches, and negligible inertia, and imaged to infer cell distributions. Notable differences in the shape of the haematocrit distributions were observed between healthy and hardened RBCs near the bifurcation apex. These lead to more asymmetric distributions for healthy RBCs in the daughter and outlet branches with cells accumulating near the inner channel walls, exhibiting distinct hematocrit peaks which are sharper for healthy RBCs. Although the hematocrit distributions differed locally, similar partitioning characteristics were observed for both suspensions. Comparisons with RBC distributions measured in a T-shaped bifurcation showed that the bifurcation angle affects the haematocrit characteristics of the healthy RBCs and not the hardened ones. The extent of RBC partitioning was found similar in both geometries and suspensions. The study highlights the differences between local and global characteristics which impact RBC distribution in more complex, multi-bifurcation networks.

AAPH-induced oxidative damage reduced anion exchanger 1 (SLC4A1/AE1) activity in human red blood cells: protective effect of an anthocyanin-rich extract

Introduction: During their lifespan in the bloodstream, red blood cells (RBCs) are exposed to multiple stressors, including increased oxidative stress, which can affect their morphology and function, thereby contributing to disease. Aim: This investigation aimed to explore the cellular and molecular mechanisms related to oxidative stress underlying anion exchanger 1 activity (band 3, SLC4A1/AE1) in human RBCs. To achieve this aim, the relationship between RBC morphology and functional and metabolic activity has been explored. Moreover, the potential protective effect of an anthocyanin-enriched fraction extracted from Callistemon citrinus flowers was studied. Methods: Cellular morphology, parameters of oxidative stress, as well as the anion exchange capability of band 3 have been analyzed in RBCs treated for 1 h with 50 mM of the pro-oxidant 2,2'-azobis (2-methylpropionamide)-dihydrochloride (AAPH). Before or after the oxidative insult, subsets of cells were exposed to 0.01 μg/mL of an anthocyanin-enriched fraction for 1 h. Results: Exposure to AAPH caused oxidative stress, exhaustion of reduced glutathione, and over-activation of the endogenous antioxidant machinery, resulting in morphological alterations of RBCs, specifically the formation of acanthocytes, increased lipid peroxidation and oxidation of proteins, as well as abnormal distribution and hyper-phosphorylation of band 3. Expected, oxidative stress was also associated with a decreased band 3 ion transport activity and an increase of oxidized haemoglobin, which led to abnormal clustering of band 3. Exposure of cells to the anthocyanin-enriched fraction prior to, but not after, oxidative stress efficiently counteracted oxidative stress-related alterations. Importantly, protection of band3 function from oxidative stress could only be achieved in intact cells and not in RBC ghosts. Conclusion: These findings contribute a) to clarify oxidative stress-related physiological and biochemical alterations in human RBCs, b) propose anthocyanins as natural antioxidants to neutralize oxidative stress-related modifications, and 3) suggest that cell integrity, and therefore a cytosolic component, is required to reverse oxidative stress-related pathophysiological derangements in human mature RBCs.

Marginated aberrant red blood cells induce pathologic vascular stress fluctuations in a computational model of hematologic disorders

Red blood cell (RBC) disorders such as sickle cell disease affect billions worldwide. While much attention focuses on altered properties of aberrant RBCs and corresponding hemodynamic changes, RBC disorders are also associated with vascular dysfunction, whose origin remains unclear and which provoke severe consequences including stroke. Little research has explored whether biophysical alterations of RBCs affect vascular function. We use a detailed computational model of blood that enables characterization of cell distributions and vascular stresses in blood disorders and compare simulation results with experimental observations. Aberrant RBCs, with their smaller size and higher stiffness, concentrate near vessel walls (marginate) because of contrasts in physical properties relative to normal cells. In a curved channel exemplifying the geometric complexity of the microcirculation, these cells distribute heterogeneously, indicating the importance of geometry. Marginated cells generate large transient stress fluctuations on vessel walls, indicating a mechanism for the observed vascular inflammation.

Cell-free layer development and spatial organization of healthy and rigid red blood cells in a microfluidic bifurcation

Bifurcations and branches in the microcirculation dramatically affect blood flow as they determine the spatiotemporal organization of red blood cells (RBCs). Such changes in vessel geometries can further influence the formation of a cell-free layer (CFL) close to the vessel walls. Biophysical cell properties, such as their deformability, which is impaired in various diseases, are often thought to impact blood flow and affect the distribution of flowing RBCs. This study investigates the flow behavior of healthy and artificially hardened RBCs in a bifurcating microfluidic T-junction. We determine the RBC distribution across the channel width at multiple positions before and after the bifurcation. Thus, we reveal distinct focusing profiles in the feeding mother channel for rigid and healthy RBCs that dramatically impact the cell organization in the successive daughter channels. Moreover, we experimentally show how the characteristic asymmetric CFLs in the daughter vessels develop along their flow direction. Complimentary numerical simulations indicate that the buildup of the CFL is faster for healthy than for rigid RBCs. Our results provide fundamental knowledge to understand the partitioning of rigid RBC as a model of cells with pathologically impaired deformability in complex in vitro networks.

Stretch-activated ion channel TMEM63B associates with developmental and epileptic encephalopathies and progressive neurodegeneration

By converting physical forces into electrical signals or triggering intracellular cascades, stretch-activated ion channels allow the cell to respond to osmotic and mechanical stress. Knowledge of the pathophysiological mechanisms underlying associations of stretch-activated ion channels with human disease is limited. Here, we describe 17 unrelated individuals with severe early-onset developmental and epileptic encephalopathy (DEE), intellectual disability, and severe motor and cortical visual impairment associated with progressive neurodegenerative brain changes carrying ten distinct heterozygous variants of TMEM63B, encoding for a highly conserved stretch-activated ion channel. The variants occurred de novo in 16/17 individuals for whom parental DNA was available and either missense, including the recurrent p.Val44Met in 7/17 individuals, or in-frame, all affecting conserved residues located in transmembrane regions of the protein. In 12 individuals, hematological abnormalities co-occurred, such as macrocytosis and hemolysis, requiring blood transfusions in some. We modeled six variants (p.Val44Met, p.Arg433His, p.Thr481Asn, p.Gly580Ser, p.Arg660Thr, and p.Phe697Leu), each affecting a distinct transmembrane domain of the channel, in transfected Neuro2a cells and demonstrated inward leak cation currents across the mutated channel even in isotonic conditions, while the response to hypo-osmotic challenge was impaired, as were the Ca2+ transients generated under hypo-osmotic stimulation. Ectopic expression of the p.Val44Met and p.Gly580Cys variants in Drosophila resulted in early death. TMEM63B-associated DEE represents a recognizable clinicopathological entity in which altered cation conductivity results in a severe neurological phenotype with progressive brain damage and early-onset epilepsy associated with hematological abnormalities in most individuals.

Red Blood Cell Partitioning Using a Microfluidic Channel with Ladder Structure

This study investigated the partitioning characteristics of red blood cells (RBCs) within capillaries, with a specific focus on ladder structures observed near the end of the capillaries. In vitro experiments were conducted using microfluidic channels with a ladder structure model comprising six bifurcating channels that exhibited an anti-parallel flow configuration. The effects of various factors, such as the parent channel width, distance between branches, and hematocrit, on RBC partitioning in bifurcating channels were evaluated. A decrease in the parent channel width resulted in an increase in the heterogeneity in the hematocrit distribution and a bias in the fractional RBC flux. Additionally, variations in the distance between branches affected the RBC distribution, with smaller distances resulting in greater heterogeneity. The bias of the RBC distribution in the microchannel cross section had a major effect on the RBC partitioning characteristics. The influence of hematocrit variations on the RBC distribution was also investigated, with lower hematocrit values leading to a more pronounced bias in the RBC distribution. Overall, this study provides valuable insights into RBC distribution characteristics in capillary networks, contributing to our understanding of the physiological mechanisms of RBC phase separation in the microcirculatory system. These findings have implications for predicting oxygen heterogeneity in tissues and could aid in the study of diseases associated with impaired microcirculation.

Red blood cell lingering modulates hematocrit distribution in the microcirculation

The distribution of red blood cells (RBCs) in the microcirculation determines the oxygen delivery and solute transport to tissues. This process relies on the partitioning of RBCs at successive bifurcations throughout the microvascular network, and it has been known since the last century that RBCs partition disproportionately to the fractional blood flow rate, therefore leading to heterogeneity of the hematocrit (i.e., volume fraction of RBCs in blood) in microvessels. Usually, downstream of a microvascular bifurcation, the vessel branch with a higher fraction of blood flow receives an even higher fraction of RBC flux. However, both temporal and time-average deviations from this phase-separation law have been observed in recent studies. Here, we quantify how the microscopic behavior of RBC lingering (i.e., RBCs temporarily residing near the bifurcation apex with diminished velocity) influences their partitioning, through combined in vivo experiments and in silico simulations. We developed an approach to quantify the cell lingering at highly confined capillary-level bifurcations and demonstrate that it correlates with deviations of the phase-separation process from established empirical predictions by Pries et al. Furthermore, we shed light on how the bifurcation geometry and cell membrane rigidity can affect the lingering behavior of RBCs; e.g., rigid cells tend to linger less than softer ones. Taken together, RBC lingering is an important mechanism that should be considered when studying how abnormal RBC rigidity in diseases such as malaria and sickle-cell disease could hinder the microcirculatory blood flow or how the vascular networks are altered under pathological conditions (e.g., thrombosis, tumors, aneurysm).

Red blood cell dynamics in extravascular biological tissues modelled as canonical disordered porous media

The dynamics of blood flow in the smallest vessels and passages of the human body, where the cellular character of blood becomes prominent, plays a dominant role in the transport and exchange of solutes. Recent studies have revealed that the microhaemodynamics of a vascular network is underpinned by its interconnected structure, and certain structural alterations such as capillary dilation and blockage can substantially change blood flow patterns. However, for extravascular media with disordered microstructure (e.g. the porous intervillous space in the placenta), it remains unclear how the medium's structure affects the haemodynamics. Here, we simulate cellular blood flow in simple models of canonical porous media representative of extravascular biological tissue, with corroborative microfluidic experiments performed for validation purposes. For the media considered here, we observe three main effects: first, the relative apparent viscosity of blood increases with the structural disorder of the medium; second, the presence of red blood cells (RBCs) dynamically alters the flow distribution in the medium; third, symmetry breaking introduced by moderate structural disorder can promote more homogeneous distribution of RBCs. Our findings contribute to a better understanding of the cell-scale haemodynamics that mediates the relationship linking the function of certain biological tissues to their microstructure.

Red blood cells tracking and cell-free layer formation in a microchannel with hyperbolic contraction: A CFD model validation

BACKGROUND AND OBJECTIVE

In recent years, progress in microfabrication technologies has attracted the attention of researchers across disciplines. Microfluidic devices have the potential to be developed into powerful tools that can elucidate the biophysical behavior of blood flow in microvessels. Such devices can also be used to separate the suspended physiological fluid from whole in vitro blood, which includes cells. Therefore, it is essential to acquire a detailed description of the complex interaction between erythrocytes (red blood cells; RBCs) and plasma. RBCs tend to undergo axial migration caused by occurrence of the Fåhræus-Lindqvist effect. These dynamics result in a cell-free layer (CFL), or a low volume fraction of cells, near the vessel wall. The aim of the paper is to develop a numerical model capable of reproducing the behavior of multiphase flow in a microchannel obtained under laboratory conditions and to compare two multiphase modelling techniques Euler-Euler and Euler-Lagrange.

METHODS

In this work, we employed a numerical Computational Fluid Dynamics (CFD) model of the blood flow within microchannels with two hyperbolic contraction shapes. The simulation was used to reproduce the blood flow behavior in a microchannel under laboratory conditions, where the CFL formation is visible downstream of the hyperbolic contraction. The multiphase numerical model was developed using Euler-Euler and hybrid Euler-Lagrange approaches. The hybrid CFD simulation of the RBC transport model was performed using a Discrete Phase Model. Blood was assumed to be a nonhomogeneous mixture of two components: dextran, whose properties are consistent with plasma, and RBCs, at a hematocrit of 5% (percent by volume of RBCs).

RESULTS

The results show a 5 μm thick CFL in a microchannel with a broader contraction and a 35 μm thick CFL in a microchannel with a narrower contraction. The RBC volume fraction in the CFL is less than 2%, compared to 7-8% in the core flow. The results are consistent for both multiphase simulation techniques used. The simulation results were then validated against the experimentally-measured CFL in each of the studied microchannel geometries.

CONCLUSIONS

Reasonable agreement between experiments and simulations was achieved. A validated model such as the one tested in this study can expedite the microchannel design process by minimizing the need to prefabricate prototypes and test them under laboratory conditions.

Circulating cell clusters aggravate the hemorheological abnormalities in COVID-19

Microthrombi and circulating cell clusters are common microscopic findings in patients with coronavirus disease 2019 (COVID-19) at different stages in the disease course, implying that they may function as the primary drivers in disease progression. Inspired by a recent flow imaging cytometry study of the blood samples from patients with COVID-19, we perform computational simulations to investigate the dynamics of different types of circulating cell clusters, namely white blood cell (WBC) clusters, platelet clusters, and red blood cell clusters, over a range of shear flows and quantify their impact on the viscosity of the blood. Our simulation results indicate that the increased level of fibrinogen in patients with COVID-19 can promote the formation of red blood cell clusters at relatively low shear rates, thereby elevating the blood viscosity, a mechanism that also leads to an increase in viscosity in other blood diseases, such as sickle cell disease and type 2 diabetes mellitus. We further discover that the presence of WBC clusters could also aggravate the abnormalities of local blood rheology. In particular, the extent of elevation of the local blood viscosity is enlarged as the size of the WBC clusters grows. On the other hand, the impact of platelet clusters on the local rheology is found to be negligible, which is likely due to the smaller size of the platelets. The difference in the impact of WBC and platelet clusters on local hemorheology provides a compelling explanation for the clinical finding that the number of WBC clusters is significantly correlated with thrombotic events in COVID-19 whereas platelet clusters are not. Overall, our study demonstrates that our computational models based on dissipative particle dynamics can serve as a powerful tool to conduct quantitative investigation of the mechanism causing the pathological alterations of hemorheology and explore their connections to the clinical manifestations in COVID-19.

Modeling a 3-D multiscale blood-flow and heat-transfer framework for realistic vascular systems

Modeling of biological domains and simulation of biophysical processes occurring in them can help inform medical procedures. However, when considering complex domains such as large regions of the human body, the complexities of blood vessel branching and variation of blood vessel dimensions present a major modeling challenge. Here, we present a Voxelized Multi-Physics Simulation (VoM-PhyS) framework to simulate coupled heat transfer and fluid flow using a multi-scale voxel mesh on a biological domain obtained. In this framework, flow in larger blood vessels is modeled using the Hagen–Poiseuille equation for a one-dimensional flow coupled with a three-dimensional two-compartment porous media model for capillary circulation in tissue. The Dirac distribution function is used as Sphere of Influence (SoI) parameter to couple the one-dimensional and three-dimensional flow. This blood flow system is coupled with a heat transfer solver to provide a complete thermo-physiological simulation. The framework is demonstrated on a frog tongue and further analysis is conducted to study the effect of convective heat exchange between blood vessels and tissue, and the effect of SoI on simulation results.

Application of machine learning in predicting blood flow and red cell distribution in capillary vessel networks

Capillary blood vessels in the body partake in the exchange of gas and nutrients with tissues. They are interconnected via multiple vascular junctions forming the microvascular network. Distributions of blood flow and red cells (RBCs) in such networks are spatially uneven and vary in time. Since they dictate the pathophysiology of tissues, their knowledge is important. Theoretical models used to obtain flow and RBC distribution in large networks have limitations as they treat each vessel as a one-dimensional segment and do not explicitly consider cell-cell and cell-vessel interactions. High-fidelity computational models that accurately model each individual RBC are computationally too expensive to predict haemodynamics in large vascular networks and over a long time. Here we investigate the applicability of machine learning (ML) techniques to predict blood flow and RBC distributions in physiologically realistic vascular networks. We acquire data from high-fidelity simulations of deformable RBC suspension flowing in the networks. With the flow and haematocrit specified at an inlet of vasculature, the ML models predict the time-averaged flow rate and RBC distributions in the entire network, time-dependent flow rate and haematocrit in each vessel and vascular bifurcation in isolation over a long time, and finally, simultaneous spatially and temporally evolving quantities through the vessel hierarchy in the networks.

Physical Properties of Blood and their Relationship to Clinical Conditions

It has been long known that blood health heavily influences optimal physiological function. Abnormalities affecting the physical properties of blood have been implicated in the pathogenesis of various disorders, although the exact mechanistic links between hemorheology and clinical disease manifestations remain poorly understood. Often overlooked in current medical practice, perhaps due to the promises offered in the molecular and genetic era, the physical properties of blood which remain a valuable and definitive indicator of circulatory health and disease. Bridging this gap, the current manuscript provides an introduction to hemorheology. It reviews the properties that dictate bulk and microcirculatory flow by systematically dissecting the biomechanics that determine the non-Newtonian behavior of blood. Specifically, the impact of hematocrit, the mechanical properties and tendency of red blood cells to aggregate, and various plasma factors on blood viscosity will be examined. Subsequently, the manner in which the physical properties of blood influence hemodynamics in health and disease is discussed. Special attention is given to disorders such as sickle cell disease, emphasizing the clinical impact of severely abnormal blood rheology. This review expands into concepts that are highly topical; the relation between mechanical stress and intracellular homeostasis is examined through a contemporary cell-signaling lens. Indeed, accumulating evidence demonstrates that nitric oxide is not only transported by erythrocytes, but is locally produced by mechanically-sensitive enzymes, which appears to have intracellular and potentially extracellular effects. Finally, given the importance of shear forces in the developing field of mechanical circulatory support, we review the role of blood rheology in temporary and durable mechanical circulatory support devices, an increasingly utilized method of life support. This review thus provides a comprehensive overview for interested trainees, scientists, and clinicians.