Combining microfluidic spleen-like filtering unit with machine learning algorithms to characterize rare hereditary hemolytic anemia

doi:10.21203/rs.3.rs-227247/v1

Download PDF

Research Article

Combining microfluidic spleen-like filtering unit with machine learning algorithms to characterize rare hereditary hemolytic anemia

https://doi.org/10.21203/rs.3.rs-227247/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Combining microfluidics technology with machine learning represents an innovative approach to conduct massive quantitative cell behavior study and implement smart decision-making systems in support of clinical diagnostics.

The spleen plays a key-role in rare hereditary hemolytic anemia (RHHA), being the organ responsible for the premature removal of defective red blood cells (RBCs). The goal is to adapt the physiological spleen filtering strategy for in vitro study and monitoring of blood diseases through RBCs shape analysis. Then, a microfluidic spleen-like filtering unit and video data analysis are combined for the characterization of RBCs in RHHA. A filtering unit is designed to measure deformability by maintaining fixed the RBC orientation to study its capacity of restoring the original shape after crossing microconstrictions. Two cooperative learning approaches are used for the analysis: the majority voting scheme, in which the most voted label for all the cell images is the class assigned to the entire video; and the maximum sum of scores to decide the maximally scored class to assign.

The proposed platform shows the capability to discriminate healthy controls and patients with an average efficiency of 91%, but also to distinguish between RHHA subtypes, with an efficiency of 82%.

Biotechnology and Bioengineering

Electronic Materials and Devices

Biomedical Engineering

RHHA

RBCs

microfluidics technology

machine learning

IES

The spleen is the organ specialized in filtering the blood, removing old and defective red blood cells (RBCs). This filtering function is specifically performed by the red pulp, where the blood is forced to flow through Inter Endothelial Slits (IES), whose diameter can reach 1 µm.^1,2

Human RBCs are 120-day lifespan cells, known for containing hemoglobin (Hb), necessary to accomplish their main function as oxygen transporters. They have a diameter and thickness of 8 and 2–3 µm respectively, presenting a unique flexibility characteristic needed to rapidly deform themselves and squeeze into the blood flow; particularly into the capillaries which diameter can reach half of the size of an RBC.^3,4

Defects affecting RBCs can lead to rare hereditary hemolytic anemia (RHHA), a group of heterogeneous inherited disorders characterized by premature removal of RBCs (hemolysis) in the spleen, as they are not able to pass through the IES. Consequently, Hb level concentration decreases, sometimes insufficient to fulfill physiology needs.⁵ RHHA RBCs can show changes in volume, plasticity, shape and deformability leading to compromised mechanical behavior.

Depending on the affected RBC’s component, RHHAs are classified as: i) hemoglobinopathies, which include sickle cell disease (SCD) and thalassemia syndromes (THAL), ii) membranopathies, as hereditary spherocytosis (HS), and iii) enzymopathies. SCD is characterized by the presence of abnormal HbS as a consequence of a single mutation in the gene encoding the Hb subunit β (HBB). Under deoxygenated conditions, HbS polymerizes leading to RBCs sickling. HbS polymerization changes shape and physical properties of RBC interfering with their flexibility and rheological properties. Repeated episodes of HbS polymerization and RBC sickling in conditions of low pO2 and unsickling in conditions of high pO2 can lead to severe alterations in the membrane structure and function and eventually result in the formation of an irreversibly sickle cell. Deformed sickle RBCs can occlude blood flow in the microvascular circulation producing vascular damage, organ infarcts, painful episodes and other symptoms.^6,7 Concerning HS, it is caused by mutations in genes encoding RBC membrane proteins, causing defects in the vertical interactions and leading to spherically-shaped RBCs with decreased deformability. Finally, Thalassemia is characterized by genetic abnormalities responsible for the absence or the decrease of the synthesis of one or more globin chains, responsible for a malfunctioning and for the lost of deformability of the RBCs. Besides the heterogeneity among patients in RHHAs pathophysiology, therapies have been traditionally non-specific, limited to symptomatic control of the anemia. Conventional techniques for monitoring patients are based on the quantitative or qualitative measurement of the different blood components, while other technologies, as ectacytometry or flow-cytometry, are not always implemented in the lab routine of patients monitoring, probably due to their high cost and not easy access.⁸ Therefore, there is a need for finding more accessible and easy techniques for RHHA studies. Indeed, the ability to monitor the disease’s phenotype would open to the possibility of the development of new treatment strategies and could have potential in prediction of disease complications and patient’s response to therapy.

Micro and nano technologies have gained an interesting role in this scenario, offering simple, low cost and rapid platform able to mimic the microvasculature properties reflecting cellular/tissue level response. Furthermore, microfluidics devices constitute also a way to facilitate and speed the efficacy validation for monitoring.⁸ Microfluidics and Lab-on-a-Chip (LoC) technology constitute a valid choice for the study of red blood cells characteristics and to mimic spleen filtering behavior, using precise manipulation of fluid dynamics.⁹ Evaluation of the spleen function based on the characterization of RBC mechanical properties, and their observation, may play a key role in the prognosis and morbidity of RHHA patients. Several microfluidics devices have been developed for the study of RBCs deformability, considered as a valid mechanical parameter for disease diagnosis^10–17. Faustino et al developed a device for the estimation of single-cell deformability through a hyperbolic converging micro channels to measure deformation and cell motion from healthy and end-stage kidney disease patients.¹¹ Lizarralde et al. fabricated a microfluidic device to mimic the mechanical stress on flowing sickle RBCs, being able to evaluate their resistance to lysis.¹² This device was also adapted to use in a bioimpendance-based approach to evaluate RBCs elasticity by electrically measuring their transit time in pathological and healthy RBCs.¹⁷ Mehri et al created a device combined with optical viscometer techniques for sorting RBCs to simultaneously measure viscosity with shear rates and aggregates sizes, and capable of maintaining RBCs integrity during cold storage and so to improve the efficiency of transfusions.¹⁴

Some works are focused on the study of RBC morphology, by applying deep learning and shape feature extraction.^18–21 Devices designed for in flux measurement, usually combined microfluidic techniques with sophisticated equipment such as laser diffraction viscometer, laser scanner electron microscope, real-time deformability cytometry, sensors or supernatants measurements with the scope to evaluate RBC deformability properties on healthy donors or in presence of blood disorders.^12,22 Most of them considered a limited number of donors and the manual user-dependent analysis by means of commercial tools (e.g., ImageJ, Adobe Photoshop).^12,19,22,23

A new frontier for microfluidic devices appears when they meet machine learning algorithms for image analysis; aiming not only to simulate organ functionalities but also to conduct massive quantitative cell behavior analysis and implement smart decision-making systems in support to clinical diagnostic. The use of non-invasive acquisition techniques such as Time-Lapse Microscop (TLM) allows label-free analysis (with no fluorescence). Then, one of the straightforward applications of such an approach is to adapt the physiological spleen filtering strategy for in vitro study and monitoring of blood diseases through Red Blood Cell shape analysis.

We propose a microfluidic platform that mimics the flow conditions of the spleen in the IES region, combined with image analysis based on deep learning algorithms for the massive study of RBCs mechanical properties with the aim to distinguish between specific types of RHHA, using for our study human blood samples from healthy donors and from patients (SCD, THAL, and HS).

Specifically, we assumed that RBCs have the capacity of restoring their original shape after crossing IES which could be used as a parameter to measure deformability, going far beyond the analysis of old and/or defective RBCs. Furthermore, the exploitation of the deep learning algorithms adds the chance to conduct user-independent and automatic analysis. The automatic localization of cells in the image, performed by the design of image analysis modules, compensates the label-free acquisition condition and prevents laborious manual cropping of each cell. TLM permits to implement a continuous monitoring using low-cost equipment for the acquisition, on the contrary of some of the previously mentioned costly devices. At present, to author’s knowledge, it does not exist a platform that combines in an automatic way microfluidic device with non-invasive Time-Lapse Microscopy (TLM), video and smart data analysis of flowing cells for RBC cells.

Figure 1 shows the study’s workflow followed to perform these experiments: samples collection in the hospital, sample preparation, video recording under the microscope and video analysis.

Samples collection

Human blood samples from healthy donors, SCD, HS, and THAL patients were collected from University Hospital Vall d’Hebron for our study. (Figure 1A).

Experimental process

RBCs separation and experimental solution

RBCs were isolated from the blood samples and diluted before perfusing them through the microfluidic device. (Figure 1B)

Considering that shear stress and changes in physiological solution can cause osmotic stress on RBC, rendering echinocytes or blocking visibility when recording videos;^3,24 we prepared different RBCs solutions and dilutions. We found that the optimal working solution for our experiments was physiological serum containing 1% BSA, 0, 25% of 0,25M EDTA and 15% of glycerol. (see supplemental Figure S1A, Supplementary material).

Figure 2 (A. Spleen filtering unit on a chip designed in the study. It consists of a main channel branched until forming eight parallel microchannels. Each microchannels contain a row of filtering funnel-shaped micro-constriction to mimic the IES section of the spleen. B. zoom out of filtering funnel-shaped constriction. C. Representation of a healthy RBCs and a RHHA RBC passing through the microconstricion. A healthy RBC deforms its shape and recovers it soon after passed the slits. On the contrary, in a RHHA patient the RBC capacity of returning to the original shape is compromise)

Spleen filtering unit on a chip and RBCs perfusion

Optical inspection for experiments were performed using an optical Zeiss microscope and videos were recorded using a mono camera coupled to microscope. A Precision Pressure Control System was used to regulate the flow pressure in the microfluidic device. Videos were recorded for analysis while RBCs were perfused through the chip (Figure 1C).

A microfluidic device was designed and fabricated with the aim of mimicking the filtering function of the red pulp’s spleen. It consisted of a main channel branched until forming eight parallel microchannels. Each microchannels contained a row of filtering funnel-shaped micro-constriction to mimic the IES section of the spleen. Due to their funnel shape, the distance between two slits varied from 1.5 μm to 6.8 μm. Of particular importance, the 1,5 μm narrowest distance defined to ensured RBC deformation when crossing a slit. Also, along each canal, there was a matrix of pillars (Figure 2A-B). This matrix was designed to mimic the reticular mesh of the spleen. Total length for channel was 9.5 mm and height was 4.5 μm. Note that the height of the device was selected to limit RBC movement and maintain it in a planar orientation for a better visualization of its shape under the microscope when deforming itself while flowing through the micro-constrictions.

Video analysis

Video analysis consisted of three steps: cell localization and automatic Region Of Interest (ROI) extraction, feature selection and data refining, and finally classification step. (Figure 1D)

Dataset collection and analysis of images from human RBCs

79 different videos were recorded and analyzed from the 32 subjects included in this study. Cells were automatically localized in each video frame, and a ROI around each cell was extracted and stored as a single data sample. The cell deforms its shape and recovers it soon after passed the slits in the case of a healthy RBC. On the contrary, when the RBC was from a RHHA patient, the RBC capacity of returning to the original shape is compromised (Figure 2C). According to this assumption, we focused the analysis on the ROIs collected after the barrier.

Table 1 shows the number of ROIs for each anemia condition (Control vs RHHA) after the slit barrier.

Sample	Number of VIDEO	Number of ROIs after the barrier
Total	79	3442
Control	30	1259
RHHA	49	2183
SCD	11	876
THAL	10	406
HS	28	901

ROI (Region Of Interest); RHHA (rare hereditary hemolytic anemia); SCD (sickle cell disease);

THAL (thalassemia); HS (hereditary spherocytosis)

Table 1. (Number of videos and ROIs analyzed in the study. The numbers in the table represent the number of videos recorded during the experiments and the number of ROIs extracted from the corresponding videos. Numbers are listed considering the total of the individuals included in the study, then divided in categories.)

Each extracted ROI was processed through a Deep Learning architecture (Figure 1D) and coded into a list of numerical descriptors with the aim of assessing the lack of deformability. Due to the fact that we used a pretrained DL network, and no further fine tuning procedure is performed over the dataset acquired in this study, this procedure is called “transfer learning”.

We implemented Leave-One-Experiment-Out (LOEO) cross-validation procedure to assess the performance of the proposed methodology. Such a strategy was preferred to demonstrate the robustness of the approach, avoiding depending on different acquisition conditions (i.e. image illumination, channel pressure, RBC positioning of the chip). In this strategy, the ROIs extracted in a given video are left out for testing and the remaining ones acquired in the other videos are used for training the model. The procedure is then repeated exhaustively over all the 79 videos considered in the analysis.

We performed two different studies for the validation procedure: first, we calculated the accuracy of recognition of healthy (label 0) vs unhealthy (label 1) videos. Then, we also considered a more challenging scenario in which we tried to recognize individual kinds of anemia, by discriminating healthy (label 0) vs SCD (label 1), THAL (label 2), and HS (label 3) categories.

Figure 3 (Confusion matrices A. Confusion matrices reporting the results of the majority voting and the unhealthy percentage limit criteria for the two-class problem. 30% and 40% were the two percentage limit values considered. B. Confusion matrices reporting the results of the majority voting and of the maximum trustiness criteria for the four-class problem)

Two-class problem: healthy vs unhealthy

In this first study, results were initially collected in terms of classification accuracy at the single cell level and then in terms of the final label assigned to the entire experiment. Two distinct cooperative strategies are used for the task: majority voting and unhealthy percentagelimit criteria.

The rationale was that even in the case of RBCs in RHHA not all RBCs show the same loss of deformability. Anyway, by applying a cooperative strategy procedure over the labels assigned to the cells in a given experiment, the approach allows understanding a global RBCs phenotype rather than an individual RBC behavior.

In addition, we also argue that RBC deformation persistence is exposed to temporal variation thus leading to the fact that in some frames the same cell appears as a “normal” cell while in some other frames it appears with a persistent loss of elasticity. Not less relevant is the fact that each frame represents the 2D view of a quasi-3D scene in which cells move in a 3D space and are visualized over a 2D domain with focus z plan automatically set by microscope. The projection errors may also contribute to the visual lack of persistence of the cell deformation.

The majority voting assigns to the video the most voted class among those assigned to the ROIs extracted. In the unhealthy percentage limit, the system assigns the unhealthy label (i.e., label=1) if the percentage of ROIs assigned to an unhealthy label exceeds a predetermined limit value (e.g., 30%, 40%, etc.) over the entire video.

Figure 3A reports the accuracy results for the majority voting and the unhealthy percentage limit criteria for the two-class problem, considering two percentage limit values of 30% and 40%. Note that, using the 30% limit value we totally recognize RHHA subjects with no false-negative values.

Four-class problem: healthy vs SCD, THAL, and HS

In the four-class problem, in addition to the majority voting cooperative strategy, we applied a maximum trustiness criterion. Given that the classification model provides a label and a score associated, in this second approach, the video is assigned to a certain class if the sum of the scores assigned to the ROIs of the same video to that class (out of the four considered) is the highest one. Figure 3B illustrates the confusion matrices of the majority voting and of the maximum trustiness criteria.

To fully understand the role of the scores and evidence the potential of the method, we also show a sketch of the scores assigned to each of the 79 videos and related ground truth label in Figure 4.

The height of the vertical bars indicates the score values assigned to each video (video index on the x-axis). Colors indicate the category assigned according to the legend to the top-right corner. The black solid stair line indicates the expected category for each video as indicated by the right y-axis labels.

Note that, for controls, not only the scores were generally higher (the range of the score is [0,1]) but also there are not any healthy videos incorrectly assigned to a different category. Regarding the three anemia conditions, the values of the scores are smaller indicating the critical task to solve, but also in this case, there are a very few errors of classification, mostly due to the misclassification between THAL and HS samples (Figure 4).

We present a methodology that combines the use of a microfluidic spleen filtering unit with machine learning data analysis from video recorded images for evaluating the premature lack of plasticity of RBCs in patients with RHHA. A better characterization of RBCs in RHHA through machine learning algorithms will enable the stratification of patients based on severity and/or response to treatment. Thus, facilitating personalized medicine and development of new treatment strategies. Beyond standard decision system implementation (e.g., single cell classification model) the considered strategies use two additional cooperative learning approaches. In the first one, i.e., the majority voting scheme, the most voted label over all the cell images from the same video is finally assigned to the entire video. In a second strategy, the maximum sum of scores assigned by the classifier to each class is finally used to decide the maximally scored class to assign. These strategies allow us overturning the problem of so called false instances (diseased cells detected as healthy cells and viceversa) leading to focus on unhealthy patient more than to single unhealthy cells. This approach starts from two main assumptions: the former is that not all cells belonging to a patient are necessarily unhealthy; the latter is that healthy cells can be appeared unhealthy due to the 2D projection errors of a 3D scene.

Previous work from our group described a functional spleen on a chip unit, incorporating the fast and slow spleen’s compartments.¹⁵ Mimicking the spleen function, only 10% of RBCs circulates through the slow compartment which involves the filtering function, as also it is estimated to occur in the spleen. We designed a filtering unit to process a higher number of RBCs mimicking the conditions experienced by human RBCs while passing through the microcapillaries of the red pulp in the spleen, which microcirculatory structures have a dimension from 1 to 3 µm.

The RBC can show two kinds of motions in a shear flow when suspended in a blood plasma-like medium. The RBC acts as a rigid cell with a flipping motion under flow in low viscosity condition. Increasing the shear rate, the RBC acquires a rolling behavior on movement, rotating into an orbit with the axis of symmetry perpendicular to the shear plane.²⁵ Then, the height of the designed microchannels aims to maintain these RBCs in a planar orientation to ensure that we can observe the deformation of their larger side when passing through the slits. Slits present a distance of 3 um in the narrowest part, essential to achieve our goal of mimicking the microcirculatory behave when the RBC deforms itself while flowing into the device.

To assess our spleen filtering unit on a chip as a valid device for the characterization of the RBCs, analysis was made using control and RHHA samples. RHHA subgroups were composed by SCD, THAL and HS, conditions in which, due to their structural defects, RBCs lose the ability to deform themselves and flew through the capillaries in the spleen. As also described in other studies, the loss of deformability can be mostly observed after the RBC passing through a microconstriction. ² Therefore, our hypothesis for analysis was that RBCs from patients and RBCs from controls should have different performances when passing through the slits. The variability of the RHHA disease and related morphological manifestation in RBCs shape are also accounted by considering a high number of human samples (79 videos from 32 different subjects) and related cell images (more than 600K).

The results show the capability of our system not only to distinguish between healthy controls and patients with an average efficiency of 91%, but also to distinguish between RHHA subtypes, with an efficiency of 82%, proving the possibility of using our platform for the characterization of the RBCs in RHHA disorders. With the maximum trustiness criterion, the number of false positives vanishes, and all the healthy subjects are perfectly recognized. As expected, the more challenging scenario exhibits lower performance with respect to the two-class problem, but it is possible to differentiate among the different RHHA.

Concerning the presence of false negatives, some factors should be taken into account as possible cause. Among others, a mild phenotype of the disease or patients under transfusion therapy as RBCs from patients are mixed with RBCs from donors. Then, not all RBCs in a RHHA patients are “unhealthy”, and not all RBCs in a healthy control are “healthy”, the “healthy” and “unhealthy” status is primarily age dependent, meaning that old RBCs will be more likely to be labeled as “unhealthy” by the platform.

Thresholds can be established for labeling a sample as “healthy” or “unhealthy” within the majority voting approach. The establishment of a specific threshold would depend on the use of the platform for patient screening or for disease evolution monitoring. For example, using a majority voting with a 30% unhealthy percentage limit criteria the platform shows high sensitivity, recognizing all RHHA subjects with no false negative value and only 10,9 % of the healthy subjects were recognized as false positive (6 of 30 individuals). In this way, no patients will be excluded.

The ability of having results based on single ROIs (cells) represents an advantage of our platform, something not possible at the moment using other techniques in the market, where the results obtained corresponds to an average of sample’s cells. This competence can have an interesting impact on the development of new techniques for the understanding of the complexity of the disease or of the heterogeneity, among patients, in response to treatment. At the moment, such monitoring techniques based on assessment of single cells are not implemented, on the contrary of other techniques for RHHA diagnosis, among these EMA binding test and Ektacytometry for HS, and Hb’s fraction for hemoglobinopathies. Indeed, choosing the proper treatment is becoming more challenging due to the recent advances in the molecular basis underlying RHHAs which have led to many new drugs and gene therapy approaches that are still in clinical trials.

The percentage of abnormal RBCs in a patient sample could be used as a parameter of disease’s evolution or of response to treatment to be monitored over time.

We may be able to observe the ongoing variations in the percentage of cells with lost ability of changing shape and to correlate this information with the clinical picture of the patients, allowing their stratification according to disease severity or prediction of acute event in SCD.

The approach of combining microfluidic platform together with deep learning described in this study constitutes an interesting tool for the study of RBCs disorders.

A low performing camera can be used with the aim to acquire a large field of view, in which many RBCs can be visualized simultaneously in the same frame. The relatively low spatial resolution achieved on each cell is compensated by the high performing image analysis module and by the flow control system. Furthermore, the use of deep learning analysis helps to overcome the difficulty faced when working with high-speed movies, such as unavailable or unaffordable recording using a speed camera.

The results presented here indicates that our platform is a valid tool for discriminating RHHA patients from healthy donors, with an efficiency of 91%, and among specific disorders with an efficiency of 82%. Nevertheless, further studies to improve sensitivity and increase robustness of the results are needed to validate the clinical use of these microfluidic platforms.

The characterization of RBCs rheological behavior opens the possibility of application for patient's stratification according to disease severity and response to treatments and consequently it would open to the chance of personalized treatment.

This platform could be beneficious also in the biobank field for the optimization of protocol for blood storage or stratification of donors or for other disease models, such as cardiovascular diseases.

Experimental process

Samples collection

The study was conducted in accordance with the Declaration of Helsinki and was approved by the Vall d’Hebron University Hospital Ethics Committee (n.367 and n.464 CPMP/ICH/135/95). Informed consent for blood samples collection was obtained from participants.

Experiments were performed collecting on ethylenediaminetetraacetic acid (EDTA) peripheral blood samples from healthy donors, SCD patients, HS patients, Thalassemia. Samples were stored at +4°C and analyzed within 48 hours after their collection.

For the normal control, subjects with normal complete blood count (CBC) parameters for anemia were taken.

RBCs separation and experimental solution

Packed RBCs were isolated from whole blood. RBCs count from CBC parameters were considered to calculate the volume of whole blood needed in order to set the number of RBCs to a same value for each experiment. The correspondent amount of whole blood was centrifuged at 0.4 rcf for 5 minutes at +4°C. After removing buffy coat and plasma, packed RBCs were suspended and washed in physiological serum (0.9% NaCl).

After centrifugation, supernatant was removed, and packed RBCs re-suspended in 1 mL of physiological serum. Our working solution consisted of physiological serum containing 1% BSA, 0,25% of 0,25M EDTA and 15% of glycerol (Sigma Aldrich). Finally, working solution containing 1% of re-suspended packed RBCs was prepared.

Spleen filtering unit on a chip fabrication

The device was fabricated using standard photolithography and softlithography techniques. A wafer was first cleaned in three consecutive solvent baths of acetone, isopropanol and water. Then, after dehydration at 150 °C a precise multistep procedure was performed in order to obtain the desire thickness for the microfluidic device. 1) A layer of SU8-2005 was spin coated above the silicon wafer through two steps (first 15 s at 500 with an acceleration of 100 rpm/s and second 30 s at 8000 rpm with an acceleration of 300 rpm/s). Then, wafer was soft baked through three steps of temperature (1 min at 65 °C, 2 min at 95°C and 1 min at 65°C) and exposed to UV light (8,03 mJ/cm²). Finally, it was post-baked through three steps of temperature (2 min at 65 °C, 3 min at 95°C and 1 min at 65° and developed using SU-8 developer for 30 s. Dimensions of the channels were verified by microscopy and profilometer measurements. The microfluidic device replicas were made of PDMS prepolymer mixture (Sylgard®184, Dow Corning), a silicone elastomer diluted 10:1 with a cross linker using microfabrication and molding techniques previously described.¹⁵ After polymerization, PDMS was peeled off and inlet and outlet holes made using a 1 μm puncher.

, PDMS and glass slides were exposed to O2 plasma, immediately pressed together and heated at 95 °C for 10 minutes to form a permanent bond.

RBCs perfusion

The experimental set-up consisted of a Precision Pressure Control System (P2CS, Biophysical Tools GMBH) to regulate the flow pressure in the microfluidic device. 1mm flexible plastic tubing (Tygon) at the inlet hole connected pump and microfluidic chip.

Experiments were performed at room temperature at a constant pressure of 100-150 mBar.

Optical inspection for experiments were performed using a microscope (Zeiss) and an axiocam 503 mono camera. While RBCs perfusion through a filtering unit on a chip, several videos were recorded for analysis.

Video preprocessing

Each video was cropped to extract a limited area of interest around the slit barrier. This procedure assisted in eliminating confounding structures that may alter the automatic cell localization result. This step was also performed to face the problem of cell velocity variability around the barrier due to unpredictable change in the flux. In fact, we observed that when an RBC cell stops in the barrier, fluid dynamic phenomena arise for which the other cells are slowed down. On the contrary, when a cell restarts opening a gap in the slit barrier it creates a suction effect that accelerates the arriving. Cell velocity variability makes the cell appearance changing over time in an unpredictable manner. As a consequence, in some cases, cells disappear or go out of focus.

Video analysis

Cell localization and automatic ROI extraction

RBCs are localized through the Circular Hough Transform, an algorithm specifically designed to localize almost circular object in a given scene. ²⁶ This algorithm has already demonstrated to be efficient in localizing cancer cells and immune cells in previous work.^27–29 A ROI of fixed size is extracted around each localized cell (Figure 1D), frame by frame, by cropping the frame in a square region around the center of the cell. In this way, we collected two sets of ROIs. ROIs whose cells appear before the barrier (ROI-before) and ROIs corresponding to cells after the barrier (ROI-after). Our rationale was that RBCs in patients with anemia exhibit loss of deformability properties. Due to this, in patients with anemia, RBCs crossing the barrier cannot rapidly reset their shape. The comparative analysis of ROI-after the slit barrier taken from control subjects and patients could evidence such assumption.

Deep learning analysis and feature extraction

Multivariate analysis is even more crucial for recognize patterns in a given dataset and to provide proof of concepts to biological phenomena. However, an even more critical step is the extraction of specific features to characterize object shape, morphology, and movement variation. The advent of Deep Learning (DL) CNN architecture has opened new possibilities for the image and video analysis community.³⁰ In fact, DL allows to analyze an image as it appears and to transform it into a vector of low-level features automatically extracted through the internal layers of the architecture. Such approach is called “transfer learning” since it allows to transfer the result of a learning procedure elsewhere performed on very different datasets of images (e.g., so-called pre-trained DL networks) into new images, by exploiting the ability of the network to extract the relevant information of a given image in the form of numerical descriptors.³¹ In this way, by applying the well-known AlexNET DL network, ³² we transformed each ROI extracted after the slit barrier into a vector of features for further pattern recognition analysis.

Feature selection

The features extracted from the internal layers of AlexNET DL network (here we used the pool5 layer that corresponded to the 3rd Max Pooling Layer), had to be preliminarily processed. In fact, the high number of information the network provided contained some redundant features that commonly are related to the background of the images (i.e., the uniform, gray-shaded areas in the ROI) or to image characteristics not related to the cell morphological changes. Such useless information should present a strong uniformity over the different ROIs in the same video sequence since they are not related to RBC shape. For this reason, a preliminary unsupervised automatic feature selection procedure was applied with the aim of eliminating features exhibiting a small variation over the entire dataset under a given limit value. The procedure reduced the number of features of a factor of 50.

RBC classification set-up for healthy vs unhealthy

Support Vector Machine (SVM) classifiers with linear kernel had been used for the scope of constructing a classification model able to discriminate experiments with RBCs from control subjects and RBCs of patients with RHHA.³³ A leave one-experiment-out (LOEO) was used as cross-validation of the method in order to obtain experiment-independent, i.e., more general results. Single cells were labelled as healthy (0) or unhealthy (1). However, the assumption was that not all RBCs cells exhibit an equal lack of deformability properties. To evidence this assumption, we also applied two kinds of cooperative learning approach: first, we applied the majority voting procedure over the labels assigned to the cells of the same experiment; second, we considered the opportunity to assign an unhealthy label to an experiment that exhibits a percentage of unhealthy cells larger than a given limit value (i.e., 30%, 40%, or larger).

RBC classification set-up for healthy vs SCD, THAL, and HS.

With the aim to test the platform to the recognition of three distinct kinds of anemia, SCD, THAL, and HS, from control subjects, we also developed a specific set-up for the four-classes problem associated. The selected classification model was a four-classes SVM trained through a one-vs-all strategy. At each cell was then assigned a normalized score to belong to each class. In the standard approach, the class with the highest score was then selected for that cell. With the aim to perform cooperative learning, we also applied majority voting over the entire experiment and also a maximum-trustiness criterion. In the latter rule, all the scores to belong to the same class over the entire experiments were summed up. The class with the highest sum of score was finally assigned to the experiment.

Acknowledgements

Authors want to acknowledge MicroFabSpace and Microscopy Characterization Facility, Unit 7 of ICTS “NANBIOSIS” from CIBER-BBN at IBEC and Dr. Ortega-Machuca and Dr. Ramón-Azcón for excellent technical assistance. This work was supported by Networking Biomedical Research Center (CIBER), Spain. CIBER is an initiative funded by the VI National R&D&i Plan 2008–2011, Iniciativa Ingenio 2010, Consolider Program, CIBER Actions, and the Instituto de Salud Carlos III (RD16/0006/0012), with the support of the European Regional Development Fund. This work was funded by the European Commission H2020-MSCA-ITN-2019, Grant Agreement N860436, “EVIDENCE” and by the CERCA Programme and by the Commission for Universities and Research of the Department of Innovation, Universities, and Enterprise of the Generalitat de Catalunya (2017 SGR 1079), it has been developed in the context of AdvanceCat with the support of ACCIÓ (Catalonia Trade and Investment; Generalitat de Catalunya) under the Catalonian ERDF operational program (European Regional Development Fund) 2014–2020.

Author contributions

M.J.L.M. and B.A.G. designed the microfluidic spleen-like filtering unit. V.R. and M.J.L.M. designed experimental protocol and performed the experiments. D.B. and M.M.M.P. contributed on samples collection. A.M., D.D.G. and E.M. designed deep learning algorithms and performed video and statistical analysis. V.R. M.J.L.M, A.M. and E.M. drafted the manuscript and prepared the figures. J.S., M.J.L.M., M.M.M.P., A.M and E.M. reviewed and edited the manuscript. All authors critically revised the final version of the manuscript.

Conflict of Interest

The authors report no conflict of interest. The authors alone are responsible for the content and writing of this article.

Data avaiability

The videos generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Lewis, S. M., Williams, A. & Eisenbarth, S. C. Structure and function of the immune system in the spleen. Sci. Immunol.4, (2019).
Pivkin, I. V. et al. Biomechanics of red blood cells in human spleen and consequences for physiology and disease. Proc. Natl. Acad. Sci. U. S. A.113, 7804–7809 (2016).
Huisjes, R. et al. Squeezing for life - Properties of red blood cell deformability. Front. Physiol.9, 1–22 (2018).
Narla, J. & Mohandas, N. Red cell membrane disorders. Int. J. Lab. Hematol.39, 47–52 (2017).
Alaarg, A., Schiffelers, R. M., Van Solinge, W. W. & Van Wijk, R. Red blood cell vesiculation in hereditary hemolytic anemia. Front. Physiol.4 DEC, 1–15 (2013).
Bunn HF. Pathogenesis and T Reatment of S Ickle C Ell D Isease. N. Engl. J. Med.337, 762–769 (1997).
Steinberg, M. H. & Sebastiani, P. Critical Review Genetic modifiers of sickle cell disease. 795–803 (2012) doi:10.1002/ajh.23232.
Alapan, Y. et al. Emerging point-of-care technologies for sickle cell disease screening and monitoring. Expert Rev. Med. Devices13, 1073–1093 (2016).
Sackmann, E. K., Fulton, A. L. & Beebe, D. J. The present and future role of microfluidics in biomedical research. Nature507, 181–189 (2014).
Mao, X. & Huang, T. J. Exploiting mechanical biomarkers in microfluidics. Lab Chip12, 4006–4009 (2012).
Faustino, V. et al. A microfluidic deformability assessment of pathological red blood cells flowing in a hyperbolic converging microchannel. Micromachines10, (2019).
Alejandra, M. et al. Lab on a Chip mechanical stress on erythrocytes in sickle cell disease †. 2975–2984 (2018) doi:10.1039/c8lc00637g.
Cluitmans, J. C. A. et al. Alterations in Red Blood Cell Deformability during Storage: A Microfluidic Approach. Biomed Res. Int.2014, (2014).
Mehri, R., Mavriplis, C. & Fenech, M. Red blood cell aggregates and their effect on non-Newtonian blood viscosity at low hematocrit in a two-fluid low shear rate microfluidic system. PLoS One13, (2018).
Rigat-Brugarolas, L. G. et al. A functional microengineered model of the human splenon-on-a-chip. Lab Chip14, 1715–1724 (2014).
Picot, J. et al. A biomimetic microfluidic chip to study the circulation and mechanical retention of red blood cells in the spleen. Am. J. Hematol.90, 339–345 (2015).
Xu, T. et al. Characterization of red blood cell microcirculatory parameters using a bioimpedance microfluidic device. Sci. Rep.10, 1–10 (2020).
Aliyu, H. A., Sudirman, R., Abdul Razak, M. A. & Abd Wahab, M. A. Red Blood Cell Classification: Deep Learning Architecture Versus Support Vector Machine. in 2018 2nd International Conference on BioSignal Analysis, Processing and Systems (ICBAPS) 142–147 (2018). doi:10.1109/ICBAPS.2018.8527398.
Alapan, Y. et al. Sickle cell disease biochip: a functional red blood cell adhesion assay for monitoring sickle cell disease. Transl. Res.173, 74-91.e8 (2016).
Lange, J. R. et al. Microconstriction Arrays for High-Throughput Quantitative Measurements of Cell Mechanical Properties. Biophys. J.109, 26–34 (2015).
Eulenberg, P. et al. Reconstructing cell cycle and disease progression using deep learning. Nat. Commun.8, 1–6 (2017).
Aurich, K. et al. Label-free on chip quality assessment of cellular blood products using real-time deformability cytometry. Lab Chip 2306–2316 (2020) doi:10.1039/d0lc00258e.
Gambhire, P. et al. High Aspect Ratio Sub-Micrometer Channels Using Wet Etching: Application to the Dynamics of Red Blood Cell Transiting through Biomimetic Splenic Slits. Small13, 1–11 (2017).
Baba, T. et al. Ultrastructural study of echinocytes induced by poly (ethylene glycol)-cholesterol. Histochem. Cell Biol.122, 587–592 (2004).
Viallat, A. & Abkarian, M. Red blood cell : from its mechanics to its motion in shear flow. 237–243 (2014) doi:10.1111/ijlh.12233.
Gonzalez, R. C., Woods, R. E. & Eddins, S. L. Digital Image Processing Using MATLAB. (Pearson Education, 2004).
Di Giuseppe, D. et al. Learning Cancer-Related Drug Efficacy Exploiting Consensus in Coordinated Motility within Cell Clusters. IEEE Trans. Biomed. Eng.66, 2882–2888 (2019).
Nguyen, M. et al. Dissecting Effects of Anti-cancer Drugs and Cancer-Associated Fibroblasts by On-Chip Reconstitution of Immunocompetent Tumor Microenvironments. Cell Rep.25, 3884-3893.e3 (2018).
Comes, M. C. et al. A camera sensors-based system to study drug effects on in vitro motility: The case of PC-3 prostate cancer cells. Sensors (Switzerland)20, (2020).
Krizhevsky, A., Sutskever, I. & Hinton, G. E. ImageNet classification with deep convolutional neural networks. Commun. ACM60, 84–90 (2017).
Shin, H. C. et al. Deep Convolutional Neural Networks for Computer-Aided Detection: CNN Architectures, Dataset Characteristics and Transfer Learning. IEEE Trans. Med. Imaging35, 1285–1298 (2016).
Alom, M. Z. et al. The History Began from AlexNet: A Comprehensive Survey on Deep Learning Approaches. (2018).
Davies, E. R. Machine Vision: Theory, Algorithms, Practicalities. (Elsevier, 2004).

No competing interests reported.

Supplementarymaterial.pdf

Download PDF

Editorial decision: Major revision
27 Apr, 2021
Reviews received at journal
25 Mar, 2021
Reviewers agreed at journal
10 Mar, 2021
Reviewers invited by journal
01 Mar, 2021
Editor assigned by journal
23 Feb, 2021
Editor invited by journal
11 Feb, 2021
Submission checks completed at journal
11 Feb, 2021
First submitted to journal
09 Feb, 2021

You are reading this latest preprint version

Combining microfluidic spleen-like filtering unit with machine learning algorithms to characterize rare hereditary hemolytic anemia

Status:

Version 1

Abstract

Figures

Introduction

Results

Discussion

Conclusions

Methods

Declarations

References

Additional Declarations

Supplementary Files

Status:

Version 1