Patent Application
20180114315
The present invention is related to the use of digital holographic microscopy (DHM) or differential digital holographic microscopy (DDHM) for the quantitative analysis of platelets or platelets aggregates.
The platelet spreading and retraction play a pivotal role in the platelet plugging and the thrombus formation. In routine laboratory, platelet function tests include exhaustive information about the role of the different receptors present at the platelet surface without information on the 3D structure of platelet aggregates.
Blood leakage at the site of endothelium failure is counteracted by a two-step process, the haemostasis: the primary haemostasis step involves thrombocytes, the platelets coming together to form a plug, while in the secondary haemostasis step, the coagulation factors of the bloodstream form a fibrin meshwork. Platelets are 2-3 μm wide of cytoplasm with biconvex disc shape exclusively found in the blood of mammals.
They form a plug in a stepwise mechanism. First, circulating platelets are recruited at the site of injury by the exposure of deep tissue structures. The adhesion of these platelets to the surface triggers the recruitment and aggregation of novel platelets. Upon platelet adhesion, surface receptors are activated.
Thanks to these receptors, platelets are bridged together through the interaction with factors such as fibrinogen or Von Willebrand factor (vWF) that present multiple binding sites. Activation also triggers morphological changes of the platelets. Upon adhesion on the wound, the round shaped platelets of the bloodstream flatten on the surface and develop filopodia which intertwine and tighten the platelet aggregate.
Thus for more than 4 decades, factors inducing platelet aggregation seemed straightforward, requiring a stimulus, a soluble protein (fibrinogen), and a membrane-bound platelet receptor (integrin αIIbβ3 or GPIIb-IIIa), leading to a simple unified model of platelet aggregation as described by K S. Sakariassen, E. Fressinaud, J P. Girma, D. Meyer and Baumgartner in: “Role of platelet membrane glycoproteins and von Willebrand factor in adhesion of platelets to subendothelium and collagen”. Ann N Y Acad Sci. 516:52-65 (1987). However, recent technical advances allowing real time analysis of platelet aggregation in-vitro and in animal models demonstrated much more complex dynamical processes than previously expected (Jakson et Al. “Dynamics of platelet thrombus formation.” J Thromb Haemost. 7 Suppl 1:17-20 (2009)).
In particular, the mechanisms by which hemodynamic conditions lead to platelets adhesion and aggregation are still incompletely understood. Actual results suggest that platelet tethering requires different receptor/ligand pairs at low (up to 1000 s−1: fibrinogen/integrin αIIbβ3) and high (up to 10000 s−1: vWF/GPIb glycoprotein bonds) shear rates of the bloodstream (SP. Jackson. The growing complexity of platelet aggregation. Blood 109:5087-5095 (2007)).
In clinical practice, platelet function tests include exhaustive information about the role of the different receptors present at the surface of platelets. These tests are mainly based on turbidimetric optical detection, multiple electrode aggregometry and flow cytometry. With these tests, it is impossible to analyze an important process after platelets adhesion, the spreading and the retraction. And yet, the spreading and the retraction play a pivotal role in the platelet adhesion and the thrombus formation.
The present invention aims to provide a method for quantitative 3D morphology of platelet aggregation.
The present invention is related to a method for the quantitative analysis of cells or cells aggregates comprising the steps of:
Preferred embodiments of the present invention disclose at least one or a suitable combination of the following features:
The present invention discloses a method able to characterize the 3D platelets or aggregates shapes by using the quantitative phase contrast imaging provided by DHM or DDHM. This original method will be of a great interest in the study of platelets physiology in clinical practice and in the development of new drugs.
DHM suitable for the present invention are for example disclosed in patent documents EP1399730, EP1631788, and EP2357539.
Blood sampling was approved by the CHU Charleroi hospital ethics committee (Comite'd'Ethique I.S.P.PC: OM008). The studies conform to the principles outlined in the Declaration of Helsinki. Venous blood was drawn from healthy donors (from the Centre Hospitalier Universitaire de Charleroi, Belgium) into tubes with 3.2% sodium citrate solution, pH 7.4. A cone and plate device (Impact-R, Diamed©) was used (
The platelet aggregates formation was induced by exposing 130 μL of whole blood in a well to laminar flow using the disposable Teflon conical rotors. After washing, images on a circumferential plane from the wells were captured by the image analyzer on impact-R, which quantifies the platelet aggregates formed, in 2D, on the surface. The results were expressed as number of aggregates detected and the average area of the aggregates, according to the shear forces. Two sets of test were performed, with the blood samples exposed to a shear rate of 100 s−1 during 20 sec and 300 sec.
The
This section describes the off-axis DHM with a source of partially spatial coherence light to record the holographic information. The configuration is shown in
A coherent source (a mono-mode laser diode, λ=532 nm) is made partially spatial coherent by focusing the beam, by the lens ML1, close to the rotating plane of the ground glass (GG). The lens L1 collimates the beam that is divided by a beam splitter BS1. The object beam reflected by BS1 illuminates the sample in the analysis chamber in transmission. The plane, on which the platelet aggregates are sticking, is imaged by the couple of lenses ML3-L2 on the CCD camera sensor. The reference beam transmitted by the beam splitter BS1 has a similar optical path, excepted that there is no analysis chamber with a sample. The two beams are interfering on the CCD sensor.
The reference beam is slanted on the sensor with respect to the object beam in such a way that a grating-like thin interference pattern is recorded. This off-axis configuration enables to implement the Fourier method to compute the complex amplitude of the object beam for every recorded frame. Thanks to the partially coherent illumination, raw holograms directly displayed on the PC screen are with low noise and fully meaningful for the operator as with an usual microscope (
The microscope lenses ML2 and ML3 are Leica 40×, NA 0.6. The camera is a JAI, with a CCD providing holograms of 1024×1024 pixels, with a pixel size of 7.4 μm×7.4 μm. The field of view is 185 μm×185 μm.
As observed in
In digital holographic microscopy, optic elements can introduce minor distortions in the background phase. This is particularly pertinent when the main parameter analyzed is the phase information that allows to quantify the morphological shapes of the platelet aggregates. For that purpose, it is necessary to implement a phase background subtraction and permanent defects elimination. However, to do that, it is necessary to insert the sample in which there are already the aggregates, making it difficult to perform the background phase subtraction on the basis of a single hologram. In order to overcome this issue, we first recorded a sequence of N holograms (in the reported experiment, N=29) of the sample moved with lateral translations. This sequence of holograms allows us to implement corrections of the defects in the intensities and the phase maps in a self-consistent way. The defects correction uses the procedures that are described as follow.
The S complex amplitude gk(s,t), where (s,t) are the discrete spatial variables, with k=0, . . . , N−1, s,t=0, . . . , n−1 and n is the pixel number by side, are extracted from a set of recorded hologram hk(s,t). The original holograms have a size of 1024×1024 pixels giving rise to complex amplitudes of the same size.
A first step of correction on the intensity field is performed. For that purpose, the averaged intensity ia(s,t) of the intensities ik(s,t)=|gk(s,t)|2 is computed. The corrected intensities are computed thanks to:
i
ck(st)=ik(s,t)/ia(s,t) (1)
For the correction of the phase maps φk(s,t) associated to gk(s,t), the averaged phase map φa(s,t) is computed and subtracted to every phase map φk(s,t) to obtain the corrected phase maps φck(s,t) according to:
φck(s,t)=mod2π{φk(s,t)−φa(s,t)} (2)
The phase background is after set, in average, to a fixed phase value. In our case, for the phase values ranging on 255 levels (1 byte), we set the background to the level 30 to avoid phase jumps. The process is efficient and is illustrated by the
On each recorded hologram, the physical heights hk(s,t) are determined by:
Where n2 is the platelet refractive index (n2=1.399) and n1 is the air refractive index (n1=1).
In the aim to determine the phase error, the regions covered by the aggregates were detected to evaluate the phase fluctuations outside those regions. Those ones are then used to establish the accuracy of the phase measurements. The detection processing that we used is identical to the one described by C. Yourassowsky and F. Dubois in “High throughput holographic imaging-in-flow for the analysis of a wide plankton size range” Opt. Express 22, 6661-6673 (2014). It is based on a high-pass filtering of the gck(s,t). Indeed, when there is locally no aggregate in some areas of the field of view, gck(s,t) is almost constant in this region and the high-pass filtering process gives complex amplitude with very low module values that are eliminated by a simple threshold operation. On the contrary, the high-pass filtering enhances the presence of an aggregate by local strong complex amplitude.
To avoid border effects by the high-pass filter, it is constituted by an inverse Gaussian filter H(u,v) defined by:
H(u,v)=(1−exp{−(u2+V2)/2σ2}) (4)
Where (u, v) are the discrete spatial frequencies (u,v=−n/2, . . . ,n/2−1), and σ is the width of the high-pass filter. In practice, σ=10 gives good results. After the inverse Fourier transformation in the filtering process, a threshold is applied. The intensity image is computed and converted into 255 gray levels in such a way that aggregates give rise to bright regions, even with saturation, and that the background regions give intensity of few grey levels (typically less than 10). The threshold level we applied is 40. It results binary images constituted by a dark background in which there are unconnected bright regions. A surface analysis is performed in order to eliminate the smaller areas (<4pixels). A morphological dilatation is performed with a 10-pixels diameter structural element in order to guarantee that the bright areas partly cover already the background region.
As we want to assess the fluctuation of the phase background, we invert the contrast to obtain a mask on which the background fluctuations are computed. An example of result is shown in
The full set of phase images is used to assess the fluctuations of the phase background by computing the standard deviation that is given by StDev=3.2 grey level. Considering this value as the typical error on the phase, we obtain, tanks to Eq. (3), that the error Δh on the height of the aggregates is 17.05 nm. By using classical statistical tools, it results that the error on the volume of an aggregate ΔV can be expressed by:
ΔV=s√{square root over (NΔh)} (5)
Where s is the pixel area and N the number of pixel covered by the aggregate. Thanks to Eq. (5), assessment of the error on the volume computation is given here below.
The volume of the aggregate a in the hologram k is obtained by computing in the corresponding phase image:
Where s is the area of one pixel, n is the number of pixels by phase image side and wka (l,m) is a region of interest function that is equal to 1 when the pixel belongs to the aggregate a, and which is zero elsewhere. wka(l,m) can be achieved by performing:
However, with this method it is not excluded to have overlaps between some neighbor aggregates that could influence the assessment of the aggregate volume. For that reason, we decided to select regions of interest that are not too close to each other to avoid the overlaps effect. This was performed to have a statistical relevance. On the basis of 29 original holograms corresponding to shear rate exposition times of 20 s and 300 s, we processed a total of 340 platelet aggregates. To observe the platelet aggregates spreading process, we exposed whole blood to a constant shear rate (100 s−1) during 20 sec and 300 sec. As first indications, for those expositions, the averaged platelet volumes are, respectively, 4.50 μm3 and 3.74 μm3 with the average errors of 0.015 μm3 in both cases. More significant information were extracted as described below. This average error of less than 0.5% is surprisingly low.
On the
DHM allows to extract aggregates maximal highs, surfaces and volumes. On the
Platelet receptors and cytoplasmic molecules, such as calpain-1 and talin involved in cascades after adhesion and activation have been extensively studied. However, the effect of these cascades on plug morphology have been commonly analyzed through the transmission microscopy alone or combined with fluorescence but in 2D. The DHM is very convenient to use and it could give more data on the role of molecules involved in the spreading. By using a flow chamber directly on the DHM, it is possible to study dynamically the formation of platelets aggregates but in 3D.
COPD (chronic obstructive pulmonary disease) is a major cause of death worldwide with estimates projecting it as the third cause of death in 2020. One of the important characteristics of this disease is the presence of both the lung and systemic chronic low-grade inflammation. In addition, observational studies show a strong statistical association between COPD and cardiovascular disease.
Systemic inflammation has been identified as a causal factor for atherosclerosis although all the mechanisms involved are not all known. The systemic low-grade inflammation, hypoxia and oxidative stress are factors that may explain the increased cardiovascular risk and mortality in COPD. During exacerbations the systemic inflammation increases, hypoxia is more marked and it has been shown with autopsies of COPD patients died during an exacerbation, that the main causes of death were due either to a decompensated heart failure (37%) or pulmonary embolism (21%). Donaldson et al suggested in 2013 that the risk of cardiovascular events was increased in these patients during an exacerbation. They report a relative risk of myocardial infarction of 2.3 (1.1 to 4.7) to five days after the onset of exacerbation. In recent years, the potential role of platelet activation in the genesis of atherothrombotic events has been demonstrated in smokers and patients suffering from COPD. During this year a study in Pneumology has been performed to study in patients (11 COPD exacerbation, 16 stable COPD and 13 controls) platelet aggregates in 3D.
The present example shows the interest of the quantitative data that can be obtained by the method of the invention in helping diagnosis of COPD complications.
The diagnosis of COPD is based on the GOLD classification (Table 1) obtained from pulmonary function tests (PFT). The pulmonary function tests are those already previously performed in stable condition or performed at discharge.
The patients aged over 45 years, with moderate to severe COPD (Stage II-III-IV) with or without hypoxemia, in stable condition or exacerbation were included in the study.
Exacerbations were defined based on the following symptoms: increase of the dyspnea, cough and/or sputum. Exacerbations were classified when infectious pathogen was found in sputum or when the patient had received antibiotics. (Based on Anthonisen criteria, an inflammatory syndrome or a radiological image compatible with pneumonia)
COPD patients were matched for age and sex with non-smoking healthy sublects.
Patients with other respiratory diseases (such as pulmonary fibrosis) or unstabilized heart disease, or a cancer or cirrhosis, were excluded. Patients taking antiplatelet agents such as clopidogrel or anti vitamin K (acenocoumarol), were also excluded.
Among 23 patients hospitalized with exacerbation of COPD, 11 patients (10 superinfected and 1 not superinfected) could be included in the study. 12 patients were excluded (4 for associated cardiac failure, pulmonary fibrosis associated to 1, 2 for failure samples, 2 for refusal, 1 Hepatitis C and 2 for taking Sintrom)
Among 22 COPD patients selected by consultation, 16 in stable condition patients were included in the study. 6 patients could not be included (2 for refusal, 2 for severe heart failure associated exacerbation and 2 for COPD exacerbation at the time of consultation).
For each patient, a careful history was made and the computer records were consulted in search of the history of home treatment, to stop a smoking, long-term oxygen therapy at home, presence an associated emphysema.
Thirteen control subjects of exacerbation in patients were taken.
Patients in exacerbation and in stable condition were subjected to arterial blood samples (to assess the degree of hypoxemia) and venous blood samples for the study of hematological parameters (number of leukocytes, platelets and red blood cells, CRP, fibrinogen). The control patients have been venous samples. Blood samples were taken within 48 hours of admission for patients in exacerbation.
The
It is the first time that an alteration of the of the platelet aggregates morphology is observed in COPD patients. It is well known that platelets play a key role in the occurrence of cardiovascular events. Based on clinical studies it would be possible to show that the alteration of the morphology of platelet aggregates could be used in the prediction of a thromboembolic event. This can be applied to other pathologies such as diabetes, hypertension, autoimmune disease, sepsis, renal failure, etc. Thus, all diseases well known to increase cardiovascular events.
The DHM allows to study the morphological dynamic of the platelets adhesion, aggregation and spreading in in-vitro models. This original method is of a great interest in the study of platelets physiology, physiopathology in clinical practice and in the development of new drugs. It is the first time that platelets aggregates are analyzed by DHM.
| Number | Date | Country | Kind |
|---|---|---|---|
| 15164996.9 | Apr 2015 | EP | regional |
| Filing Document | Filing Date | Country | Kind |
|---|---|---|---|
| PCT/EP2016/059166 | 4/25/2016 | WO | 00 |
