Jonathan T.
Butement
*ab,
Paul M.
Holloway‡
bc,
Joshua A.
Welsh
b,
Judith A.
Holloway
bd,
Nicola A.
Englyst
b,
Peter
Horak
a,
Jonathan
West
bcd and
James S.
Wilkinson
ad
aOptoelectronics Research Centre, University of Southampton, UK. E-mail: jb3006@soton.ac.uk
bFaculty of Medicine, University of Southampton, UK
cCentre for Hybrid Biodevices, University of Southampton, UK
dInstitute for Life Sciences, University of Southampton, UK
First published on 2nd March 2020
Size measurement of extracellular vesicles is hampered by the high cost and measurement uncertainty of conventional flow cytometers which is mainly due to the use of non-specialised free space optics. Integrated cytometry, where the optics and fluidics are embedded in a monolithic chip shows promise for the production of low cost, micro-flow cytometers dedicated for extracellular vesicle (EV) analysis with improved size measurement accuracy and precision. This research demonstrates a unique integrated cytometer for sub-micron particle size measurement using multi-angle scattering analysis. A combination of three technologies is used: (i) Dean-based hydrodynamic focussing to deliver a tight sample core stream to the analysis region, (ii) integrated waveguides with multimode interference devices to focus a narrow excitation beam onto the sample stream, and (iii) an angular array of collection waveguides to measure particle scattering distribution and calculate diameter. Low index 200 nm liposomes could be detected and polystyrene size standards as small as 400 nm diameter could be measured with an uncertainty of ±21 nm (1/2 IQR) demonstrating a first step on the path to high performance integrated cytometry of EVs.
The small size of EVs presents a challenge for sensitive, precise and accurate size measurement. A variety of approaches are currently in use, each with various advantages and limitations which are discussed in detail in recent reviews.6,7 Two widely used techniques are nanoparticle tracking analysis (NTA) and conventional flow cytometry. NTA involves imaging EVs in suspension under dark-field illumination and digitally tracking the Brownian motion of individual particles to calculate size using the Stokes–Einstein equation. The technique can detect particles down to 50 nm in diameter and can be combined with fluorescence to indicate EV phenotype.8 The technique has a relatively low throughput however at around 16 particles per second9 compared to flow cytometry and concentration determination requires careful calibration and can be inaccurate with broad EV size distributions.10 Conventional flow cytometry has a higher limit of detection than NTA at 200 nm but is proving to be highly practical for EV analysis as it can measure the size and concentration of polydisperse samples at kHz rates11,12 along with detecting multiple molecular markers by antibody-based fluorescent labelling. The internal apparatus of a conventional flow cytometer has remained largely unchanged since its original inception and consists of a free space optical system of lenses, mirrors and filters aligned with a fluidic flow cell as well as photodetectors and data acquisition electronics. EVs are delivered as a tight core surrounded by sheath flow into one or more laser beams which excite scattering and fluorescence from antibody labels.13 The scattering signal is used to measure particle size and the fluorescence signal from antibody labels is used to detect phenotype.
The use of conventional cytometers for EV analysis has two primary drawbacks however. Firstly, the machines have substantial costs due to their complexity which limits their widespread use. Secondly and significantly, sizing accuracy and precision is limited due to the effects of Mie resonances which are prominent at the EV size range. To measure particle size on a conventional cytometer the magnitude of the optical power scattered into the side scatter channel (SSC) is typically used as a metric. The arbitrary SSC values are obtained for a range of diameters of size standards of known refractive index (usually polystyrene or silica) and these measurements are used to normalise the arbitrary scattered power calculated from a Mie scattering model onto the same scale. Literature refractive index values for EVs can then be used with the normalised Mie scattering model to calculate EV diameter from SSC values from samples of unknown diameter.12,14 Sizing inaccuracy and imprecision arise as the relationship between particle diameter and SSC value is non-monotonic, exhibiting plateaus and local minima at specific diameter ranges leading to ambiguity in sizing in these regions. An accurate model also requires knowledge of the system collection geometry which is often proprietary and can vary with optical alignment.
Multi-angle scattering approaches have been developed to improve sizing performance. These use the relative angular scattering distribution produced by a particle to calculate diameter as opposed to the magnitude of scattered power. Here the relative scattering distribution of a particle is obtained and an analytical model is fitted to calculate particle diameter. High accuracy and precision can be achieved because there are large form changes in scattering distribution for small incremental changes in particle diameter. There is also no need to take calibration measurements from size standards to normalise a scattering model. High performance multi-angle scattering approaches have been developed using cytometers with free space optics15–18 but have not been demonstrated in an integrated cytometer. Such optofluidic technology which combines microfluidics and integrated optical waveguides19–22 offers an ideal platform with which to develop cost effective and high performance flow cytometers dedicated for EV analysis using multi-angle scattering. For example, as will be demonstrated in this paper, an angular array of collection waveguides23 can be used to measure the light scattered at multiple discrete angles providing a relative scattering distribution for model fitting to calculate diameter. Further advantages are conferred by the intrinsic and permanent alignment of the monolithic platform, which effectively removes signal drift and background noise that would be caused by the vibration of independent optical and fluidic components. In addition, chips can be mass produced with conventional microfabrication techniques greatly lowering cost per device. In this contribution, we present the design, fabrication and characterisation of an integrated cytometer for measuring particle diameter in the extracellular vesicle size range using multi-angle scattering analysis.
Carefully designed on chip lenses were required to shape the excitation beam for intersection with the particles flowing along the centre of the fluidic channel, in order to optimise the performance of the integrated cytometer. The excitation beam dimensions orthogonal to the particle flow, here referred to as beam height, ideally had to exceed the height of the particle flow stream in order to ensure all flowing particles are illuminated. For this reason a beam launched into a microfluidic channel that diverged due to diffraction orthogonally to the waveguide plane was used. This also avoided the need for complex multilayer optical structures for vertical beam shaping. The vertical divergence of the beam crossing the microfluidic channel is shown in Fig. 2(c) and was estimated with an analytical approximation for a diverging Gaussian beam with the initial dimensions of the waveguide fundamental mode.27 The theoretical beam height of ≈3.5 μm in the centre of the channel dictates the largest detectable particle and is well outside the size range of EVs.
The excitation beam dimensions parallel to the particle flow and waveguide plane, here referred to as beam width, was required to be narrow in order to reduce the possibility of signals resulting from multiple particles simultaneously occurring within the beam14 and to minimise the amount of background light leaking into the collection waveguides. Previous work has demonstrated that multimode interference devices (MMIs) are well suited for use as on chip lenses25 and suggested their use in flow cytometry. We have implemented an MMI in our integrated cytometer design: the concept is shown schematically in Fig. 1(a). MMIs consist of a multimode waveguide with a narrow single mode input. Multiple modes are excited in the multimode waveguide and interfere with each other to produce a self-imaging phenomenon whereby the narrow single mode input is reimaged at a defined distance, known as the imaging length, from the start of the MMI.28 The MMI can be positioned relative to a microfluidic channel to focus light in the horizontal plane mid-way across the microfluidic channel.
A numerical simulation package (Beamprop, RSoft Synopsys) which uses the beam propagation method (BPM) was used to calculate imaging length for a range of MMI widths and compared to an analytical approximation as detailed in Hunt et al.25 The MMIs were designed to produce a single focal spot located mid-height and at the lateral midpoint of the fluidic channel. A restricted self-imaging configuration was chosen where the input light field from the single mode waveguide is symmetrical about the longitudinal axis of the multimode waveguide and excites only symmetrical modes. This allows self-imaging over a shorter MMI length and therefore a smaller device footprint. As the imaging length of an MMI is proportional to the square of its width, a narrower MMI was desirable to again offer a smaller device footprint and reduce the effect of fabrication tolerances on focal position. The numerical and analytical calculations presented in Fig. 1(b) show that an MMI width of 20.0 μm produced an imaging length at 1266 μm for a wavelength of 532 nm, enabling compact chip designs.25,28
The effects of the fluidic system on the performance of the waveguide circuit were considered. The inclusion of strip waveguide layers on a chip creates ridges in the top surface of the chip, which would make sealing a flat enclosing manifold against the top surface difficult. In order to remove these ridges and create a flat surface for bonding to a PDMS (polydimethylsiloxane) manifold, a fluidic border region was included in the optical layer design in a method similar to that used by Friis et al.29 This is shown in the schematic of the final design in Fig. 2(b), between the fluidic channel and waveguides. At the optical analysis region the fluidic border consists of a 10 μm deep border surrounding the fluidic channel. Light leaving the MMI must propagate through the fluidic border region which, having no lateral optical bounds, acts as a slab waveguide. A fluidic border region of too great a depth will therefore alter the imaging length of an MMI and so must be kept to the minimum depth required for effective fluidic sealing. Based on experimental experience a border depth of 10 μm was found to be acceptable for sealing against a PDMS manifold, and BPM simulations showed that a 10 μm deep fluidic border did not significantly affect the imaging length of a 20.0 μm wide MMI.
The final device design is shown as an overview in Fig. 2(a). Two MMIs are positioned to provide incident beams at 45° and 90° to the particle flow as shown in Fig. 2(b). The placement angles of the input MMIs are adjusted using Snell's law to account for refraction at the microfluidic channel walls. In this work only the 90° MMI was used for scattering excitation experiments. Importantly, having the additional optional input channel at 45° to particle flow allows the possibility of simultaneous fluorescence excitation for phenotyping EVs labelled with fluorescent antibodies. A wide multimode waveguide without an MMI is positioned to collect fluorescently emitted light at 90° to the fluorescence excitation beam for this future inclusion.
To collect the light scattered from particles at discrete angles in order to measure scattering distribution, 13 waveguides, 12 μm wide, are positioned at 12° intervals over a range of 0–72° relative to the excitation beam to create an angular array, as shown in Fig. 2(a) and in expanded detail in Fig. 2(b). The 12° angular interval of the waveguide placement was designed to be close to the waveguide full acceptance angle of 16.1° providing a balance between collecting sufficient optical power for a strong signal and providing 6 collection channels, an adequate number for least squares fitting of the theoretical scattering distribution to calculate particle diameter. Again, the placement angle of the collection waveguides is adjusted to account for refraction at the microfluidic channel wall. The output facets of the scatter collection waveguides at the chip edge are evenly spaced with a pitch of 1 mm to allow butt coupling to a fibre array or directly to a photomultiplier array of matching pitch.
The microchannel height in the integrated cytometer was constrained to 20 μm, as this was the maximum channel height which could be economically fabricated using the processes available in an academic setting. Channels of greater depth would require deposition of a thicker cladding and deeper etching greatly increasing fabrication time and cost. The microchannel width at the optical analysis region was chosen to be 20 μm in order to avoid disruption to MMI focussing. These dimensions substantially exceed the channel width:particle diameter maxima of 7, thereby avoiding inertial focussing of particles towards the channel walls, that would otherwise vertically migrate half of the particulate sample outside the excitation beam.31,32
A solution of 10 μM fluorescein sodium salt in phosphate buffered saline was used to image the sample stream in the Dean focussing element. Flow performance was assessed with standard fluorescence microscopy to straightforwardly infer the vertical compression by measuring the lateral stretching of the sample stream as a percentage of the channel width. Confocal microscopy was used to measure the lateral and vertical dimensions of the sample stream in the optical detection region of the microchannel. The sidewalls of the channels were visualised by staining the device for 30 minutes with 100 μg mL−1 of the co-polymer poly-L-lysine grafted to poly(ethylene glycol) and labelled with tetramethylrhodamine (PLL-g-PEG TRITC, SuSoS AG) immediately after plasma bonding. A Leica SP8 confocal microscope with a 20× objective and pinhole set to 1 airy unit was used to obtain confocal Z stack images with a 0.5 μm step. Image deconvolution was performed using the Huygens Essential software with a theoretical point spread function.
For the scattering measurements of size standard particles, the waveguides set at 12–72° to the excitation beam axes were used as shown in Fig. 3(c). The 0° waveguide functioned as a beam dump. The mirrored set of scattered collection waveguides at equivalent angles were not used, due to a limited number of DAQ card input channels. The scatter collection array outputs were butt-coupled into an array of multimode fibres (GR-S50/125-20P, Nufern, USA) of matching pitch mounted on a v-groove block under 5 axis positional control. The v-groove block had a total of 16 fibres available but only 6 fibres were used by offsetting the v-groove block from the chip. This was to prevent collection of light from the beam dump waveguide at 0° which could damage the PMT and could be avoided in future with adjustment of the fibre array configuration. In this case, the output end of the multimode fibre array was also mounted on a v-groove block of identical layout to the input end. This output v-groove block was butt-coupled to a photomultiplier tube linear array module (H11459-20, Hamamatsu, Japan) with each fibre aligned with an individual photocathode element on the PMT array.
To characterise the shape of the beam focussed in the fluidic channel by the MMI, the beam was imaged by pumping a 0.4 μM solution of the fluorescent dye R-phycoerythrin through the channel at 0.1 m s−1 and imaging from above as shown in the inset fluorescence micrograph in Fig. 3(c) and false colour images in Fig. 4.
Refractive index values for the liposome components at a wavelength of 532 nm were obtained from the literature and were 1.44 for the DPLC outer layer33 and 1.3367 for the PBS core.34 The refractive index of an entire liposome was estimated using the volumetric ratio of its components assuming a membrane thickness of 4 nm (ref. 33) and was calculated to be 1.3486. For acquisition of the scattering signals from a particle suspension, 3D fluidic focussing conditions were established at an average flow velocity of 0.75 m s−1 in the Dean curve channel, accelerating after lateral focussing to an average flow rate of 2.25 m s−1 within the analysis microchannel. The flow rate was held for 1 minute to ensure flow stabilisation before data acquisition commenced. Light was coupled into the excitation waveguide in the TM polarisation. Data from the six channels of the PMT were acquired on a 14 bit PCIE DAQ card (TPCE-2016-8S DAQ Module, Elsys, Switzerland). The data channels recorded a voltage proportional to the optical power incident on each of six individual PMT elements coupled to the scatter collection array waveguide outputs.
The DAQ card was set to AC coupling to remove constant background. A particle signal from the 12° scatter collection waveguide channel was used to trigger data acquisition of signal pulses on all channels. To set the trigger threshold, a background measurement of the 12° PMT voltage, 8 MS (Mega Samples) in length, was acquired at 20 MHz with fluid flowing under 3D focussing conditions without particles in order to acquire the noise on the detector when no particles are passing which was represented as the standard deviation of the detector voltage during this acquisition period. The trigger threshold voltage was then set at three standard deviations of this background signal above zero with a hysteresis of one standard deviation on a positive slope. A voltage above this threshold on the 12° channel was deemed to be due to particle-induced scattering and was used to trigger simultaneous acquisition of the pulse on all scatter collection waveguide channels at 20 MHz with a 1 KS acquisition length. A minimum of 1500 pulses were recorded for suspensions of each particle size.
The results from the experimental matrix to determine the optimum Dean flow geometry are presented in Fig. 5. Increasing the arc angle from 90° to 180° doubles the transport time and focussing duration to further the lateral stretch of the sample stream across the channel.35 Although the magnitude of the Dean secondary flows scales with (H/2R)0.5, larger bend radii also have extended travel and therefore focussing times. In our experiments a 50 μm radius of curvature produced marginally improved focussing than the 25 μm radius of curvature and was chosen with the further benefit of more available area to create a pressure-tolerant bond at this critical feature. The 10 μm-wide channel produced slightly greater proportional lateral stretching, albeit at the risk of increased local pressure and problems of clogging during prolonged operation. Based on the above observations the optimum geometry for reliable Dean focussing to a tight sample core involved a radius of curvature of 50 μm, a 180° arc angle and a channel width of 20 μm, and was implemented in the final integrated cytometer design.
The vertical compression of the sample stream in this optimised curved channel geometry was measured with confocal microscopy. It was found that increasing the average fluid flow velocity from 0.33 m s−1 (De = 2.94) to 0.66 m s−1 (De = 5.88), compressed the sample stream from a height of 8 μm to 6 μm. In addition, it was found that using a sample flow to sheath flow rate ratio of 1:7 was too unstable for long-term use as sheath fluid would occasionally back-flow into the sample channel. Reducing the sample to sheath flow ratio to 1:3 prevented these back-flows and still maintained a tight sample core with a height of 6 μm and a width of 4 μm which is shown in the confocal micrograph in Fig. 5(c) with an overlaid beam profile approximated using the equation for a propagating Gaussian beam launched from the waveguide.27 Here the sample stream passes below the excitation beam and is also slightly taller than the beam height of 3.9 μm and so some sample is expected to pass outside of the beam. Whilst this was expected to increase signal variability in absolute scattered power between individual particles, the effect on the shape of the relative scattering distribution was expected to be minimal and so would not critically effect the calculation of particle diameter. The flow focussing was sufficiently compact to ensure at least half the particles encountered the beam which would provide a sufficient number of measurements from a sample. Faster flow rates would be expected to further vertically compress the sample stream but were beyond the pressure tolerance limits of the PDMS microfluidic circuit used in the matrix study. The silica channel of the integrated cytometer, which was sealed with a flat PDMS manifold, showed a greater pressure tolerance and so here the average flow velocity in the Dean focussing channel could be increased to 0.75 m s−1 to achieve even greater vertical compression. This average flow velocity in the Dean focussing channel was used for the subsequent scattering measurements of size standard particles.
Fig. 6(a) shows an example of synchronised voltage pulses from the PMT array recorded from scattering from one individual 700 nm diameter particle into 6 collection waveguides set from 12–72° at 12° increments. ESI† Fig. S1 shows example voltage pulses for the other nominal particle diameters measured. In all cases as the collection angle increases, both the peak amplitude of the pulse and the pulse area decrease as less light is scattered into collection waveguides at greater angles to the excitation beam, as expected from scattering in the Mie regime. The pulse shapes exhibit multiple sub-peaks, which we propose are caused by optical interference between a fraction of the excitation beam light directly entering the collection waveguides and the phase-shifted light scattered by the particle. Specifically, if Pp(t) is the optical power scattered from the particle into one of the waveguides as a function of time and Pb is the power of background scattering e.g. from the fluid and from optical interfaces, then the total detected power is
(1) |
E(t) = ∫Pp(t)dt + ∫Pbdt = ∫Pp(t)dt + const | (2) |
The scattering distributions of thousands of flowing particles for each nominal diameter of size standard were measured using the device. For each individual particle, the pulse energies for each angular scattering channel were converted to a relative scale with the 12° channel corresponding to unity. Fig. 6(b) visualises the particle scattering measurements as histograms of relative pulse energy for each angular collection channel for a sample of 700 nm nominal diameter particles. The robust coefficient of variation (rCV = 100 × 0.75 × sample interquartile range/median relative pulse energy)13 is displayed inset for each angular channel as a metric for variation about the median pulse energy. Pulse energy histograms for other nominal diameters are shown in ESI† Fig. S2 and sample statistics are shown in ESI† Table S1. For all particle samples measurement variation increases with angle because less light is scattered into the wider-angle channels whilst random noise remains constant.
Data were collected in the same way with monodisperse samples ranging in diameter from 400 nm to 1600 nm. Particle size was calculated from each individual particle scattering distribution by fitting a model using MATLAB® (MathWorks, US) based on Mie scattering theory.2,12,14 This model predicts the angular scattering of light from a spherical particle, such as EVs36 and the relative optical energy collected by each collection waveguide in the integrated cytometer. An optimal fit was achieved for each particle by minimizing the weighted sum of squared residuals
(3) |
Fig. 7(a–d) represents the scattering distributions recorded for each sample of size standards as median relative pulse energies. To visualise the goodness of fit of the Mie scattering model, the scattering distribution of the median particle size calculated by the model is overlaid on the experimental data. The distribution of the calculated particle diameters for each size standard sample are shown in histograms in Fig. 7(e–h) with the sizing statistics for the sample inset. Particles as small as 400 nm nominal diameter could be detected. Sizing accuracy is comparable for the 400 nm, 700 nm and 1000 nm nominal diameter samples to within 34 nm of true size whilst the 1600 nm nominal diameter particles are within 116 nm. Sizing precision for all samples is between 3–7% rCV and is within the same order as the manufacturer stated CV of 1.5%. Absolute precision in diameter measurement, represented by the IQR, is highest for the 700 nm diameter sample because the waveguide configuration is most sensitive to scatter distribution form changes in this region. The scattering distribution for a ≈700 nm diameter particle is shown in Fig. 2(b), and exhibits a single local minimum in relative scattered energy around 46° which falls within the collection aperture of the 48° collection waveguide. Thus there is a large change in scattering distribution form for small changes in particle size whilst noise remains constant. Conversely, accuracy and precision are lowest for the 1600 nm nominal diameter particles, which are larger than typical EVs, because the number of local minima in the scattering distribution for particles at this size range approaches the number of waveguides and so local minima cannot be resolved. The change in scattering distribution form with diameter is therefore very small for particles ≈1600 nm in diameter and so there is an increased effect of random noise on the diameter calculation. Whilst this device was designed for sub-micron particles, future designs could tailor the number of waveguides or the excitation wavelength to a different size range of interest in order to ensure waveguide spacing is most sensitive to Mie resonance positioning. Overall the sensitivity accuracy and precision of the integrated cytometer for sizing polystyrene particles is comparable to the high performance achieved by free-space, multi-angle flow cytometers but crucially it has been demonstrated in a fully integrated chip.15
The integrated cytometer was also used to analyse samples of 200 nm diameter liposomes (n ≈ 1.35) in order to test the device performance with particles that more closely match the scattering response of EVs (n ≈ 1.43), which were not available. Fig. 8(a) presents an example detector pulse set recorded for an individual liposome showing that the 12° degree channel displays a clear pulse while the wider angle channels pulses are hidden within the background noise and so the full scattering distribution is not measurable. Without the full scattering distribution the least squares fitting process shown in ESI† Fig. S6 could not determine an accurate liposome diameter. The additional sensing challenge posed by EVs and liposomes is presented in Fig. 8(b) which compares the model-calculated scattering distributions for liposomes, EVs and polystyrene particles on an absolute scale and shows that EVs and liposomes scatter 1–2 orders of magnitude less light than polystyrene particles and thus require further improvements to the system in order to be measured. While these results do not show size determination they do demonstrate detection of a very low index particle in an integrated format and show promise that the integrated cytometer can be optimised for EV measurements.
The size measurement performance for low index particles can be improved by increasing the device signal to noise ratio. The scattering signal magnitude can be increased by using a greater excitation power and by individually increasing the gain on each PMT anode in the array. Background noise can be reduced by eliminating vibration from the system by gluing all optical connections and by using a laser source with improved power stability. Greater fluidic confinement of flowing particles in a tighter core would also reduce any effects of positional variation on the scattering distribution shape. Further improvements can also be made to the fitting process to improve size determination accuracy and precision. Characterising the individual optical losses in each optical channel including the propagation and coupling losses in the waveguides and fibre optics would reduce the effect of any systematic error on model fitting. In addition using a greater number of collection channels would improve fitting performance. Providing sufficient light is still collected, this could be achieved by using a new design with more waveguides with a smaller collection angle increment or by using the 45° excitation MMI in the existing design to collect light over 5–137°.
The data used in this publication can be accessed at DOI: https://doi.org/10.5258/SOTON/D1173.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI:10.1039/c9lc01182j |
‡ Now with: Investigative Medicine Division, Radcliffe Department of Medicine, University of Oxford, UK. |
This journal is © The Royal Society of Chemistry 2020 |