the plasma samples included in this study were collected from patients suspected of having lung cancer due to clinical symptoms .
patients were examined from may 2011 until january 2014 at the department of pulmonary medicine , aarhus university hospital , denmark .
one hundred and nine patients with nsclc adenocarcinoma stage iiia  iv were consecutively collected .
a control group was made of 110 patients who were diagnosed to be cancer free based on computed axial tomography and , if relevant , negative biopsies .
are frequently observed amongst patients with lung cancer , acute or chronic infection or inflammation of any kind were not considered reason for exclusion . by doing so
, a control group resembling the clinical situation in a diagnostic outpatient department could be made .
microarray printing was performed on a spotbot extreme protein edition microarray printer with a 946mp4 pin ( arrayit corporation , ca , usa ) . as a positive control 100 g / ml of biotinylated human igg
antibodies were printed on epoxy - coated slides ( 75.6 mm25.0 mm ; schott nexterion , de ) , which were left to dry at room temperature overnight prior to further analysis .
a total of 37 anti - human antibodies were used and listed in the following with the corresponding clone , if available : cd146 ( p1h12 ) , flotillin-1 , hb - egf ( 4g10 ) , her3 ( 2f9 ) , her4 ( h4.77.16 ) , hsp90 ( igf1 ) , n - cadherin ( 8c11 ) , p53 ( pab240 ) , tag72 ( 0.n.561 ) , and tsg101 ( 5b7 ) from abcam ( ma , usa ) ; sftpd ( vif11 ) from acris antibodies gmbh ( de ) ; egfr and egfrviii from antibodies-online.com ( de ) ; cd63 ( mem-259 ) from biolegend ( ca , usa ) ; her2 ( 29d8 ) from cell signalling ( ma , usa ) ; cd9 and cd81 from lifespan biosciences ( wa , usa ) ; spa ( 6f10 ) from novus biologicals ( co , usa ) ; cd13 ( 498001 ) , cd14 ( 50040 ) , cd142 ( 323514 ) , cd151 ( 210127 ) , cd206 , cea ( 487609 ) , muc1 ( 604804 ) , tnf ri , tnf rii and tspan8 ( 458811 ) from r&d systems ( mn , usa ) ; areg ( s-13 ) , epcam ( 0.n.277 ) , mucin16 ( x306 ) , ny - eso-1 ( e978 ) and plap ( 8b6 ) from santa cruz biotechnologies ( tx , usa ) ; cd171 from sigma - aldrich ( mo , usa ) ; c - met ( 016 ) and pd - l1 from sino biological inc .
( china ) ; cd163 ( mac2 - 158 ) from trillium diagnostics , llc ( me , usa ) .
all antibodies were diluted with pbs with 5% glycerol and printed in triplicates at 75200 g / ml .
( 15 ) with modifications . in short , the microarray slides were initially blocked ( 50 mm ethanolamine , 100 mm tris , 0.1% sds , ph 9.0 ) prior to incubation with 10 l plasma sample diluted ( 1:10 ) in wash - buffer ( 0.05% tween20 in pbs ) .
the incubation was performed in multi - well hybridization cassettes ( arrayit corporation ) at rt for 2 hours followed by overnight incubation at 4c . following a short wash ,
the slides were incubated with a cocktail of biotinylated detection antibodies ( anti - human - cd9 , -cd63 and -cd81 , lifespan biosciences , wa , usa ) diluted 1:1,500 in wash - buffer . after a wash ,
a subsequent 30-minute incubation step with cy5-labelled streptavidin ( life technologies ) diluted 1:1,500 in wash - buffer was carried out for detection . prior to scanning
, the slides were washed in wash - buffer , then in ultrapure / deionized water and finally dried using a microarray high - speed centrifuge ( arrayit corporation ) . scanning and spot detection was performed as previously described ( 15 ) .
graphs and statistics were carried out using graphpad prism ( version 6.04 , graphpad software , inc .
, san diego , ca , usa ) , excel ( version 2013 , microsoft , redmond , wa , usa ) and roccet ( version 03/17/2014 , tmic , ab , canada ) . heat maps and hierarchical cluster analysis
were produced using genesis ( version 1.7.6 , igb tu graz , graz , austria ) .
random forests was used for classification and monte - carlo cross validation ( mccv ) using balanced subsampling was employed to compute roc curves and to calculate confidence intervals of their area under the curve ( auc ) . in each mccv , two thirds ( 2/3 ) of the samples were used to evaluate the feature importance .
the top important features were then used to build a classification model , which was validated on the 1/3 of the samples that were left out .
five hundred iterations were used to calculate performance of the model . prior to analysis missing value estimation was performed by replacing missing values by  min .  a top 20 of ratios ( based on p - values ) was computed and included into the analysis prior to data transformation and normalization by log2 .
the plasma samples included in this study were collected from patients suspected of having lung cancer due to clinical symptoms .
patients were examined from may 2011 until january 2014 at the department of pulmonary medicine , aarhus university hospital , denmark .
one hundred and nine patients with nsclc adenocarcinoma stage iiia  iv were consecutively collected .
a control group was made of 110 patients who were diagnosed to be cancer free based on computed axial tomography and , if relevant , negative biopsies .
are frequently observed amongst patients with lung cancer , acute or chronic infection or inflammation of any kind were not considered reason for exclusion . by doing so
, a control group resembling the clinical situation in a diagnostic outpatient department could be made .
microarray printing was performed on a spotbot extreme protein edition microarray printer with a 946mp4 pin ( arrayit corporation , ca , usa ) . as a positive control 100 g / ml of biotinylated human igg
antibodies were printed on epoxy - coated slides ( 75.6 mm25.0 mm ; schott nexterion , de ) , which were left to dry at room temperature overnight prior to further analysis .
a total of 37 anti - human antibodies were used and listed in the following with the corresponding clone , if available : cd146 ( p1h12 ) , flotillin-1 , hb - egf ( 4g10 ) , her3 ( 2f9 ) , her4 ( h4.77.16 ) , hsp90 ( igf1 ) , n - cadherin ( 8c11 ) , p53 ( pab240 ) , tag72 ( 0.n.561 ) , and tsg101 ( 5b7 ) from abcam ( ma , usa ) ; sftpd ( vif11 ) from acris antibodies gmbh ( de ) ; egfr and egfrviii from antibodies-online.com ( de ) ; cd63 ( mem-259 ) from biolegend ( ca , usa ) ; her2 ( 29d8 ) from cell signalling ( ma , usa ) ; cd9 and cd81 from lifespan biosciences ( wa , usa ) ; spa ( 6f10 ) from novus biologicals ( co , usa ) ; cd13 ( 498001 ) , cd14 ( 50040 ) , cd142 ( 323514 ) , cd151 ( 210127 ) , cd206 , cea ( 487609 ) , muc1 ( 604804 ) , tnf ri , tnf rii and tspan8 ( 458811 ) from r&d systems ( mn , usa ) ; areg ( s-13 ) , epcam ( 0.n.277 ) , mucin16 ( x306 ) , ny - eso-1 ( e978 ) and plap ( 8b6 ) from santa cruz biotechnologies ( tx , usa ) ; cd171 from sigma - aldrich ( mo , usa ) ; c - met ( 016 ) and pd - l1 from sino biological inc . ( china ) ; cd163 ( mac2 - 158 ) from trillium diagnostics , llc ( me , usa ) .
all antibodies were diluted with pbs with 5% glycerol and printed in triplicates at 75200 g / ml .
the microarray slides were initially blocked ( 50 mm ethanolamine , 100 mm tris , 0.1% sds , ph 9.0 ) prior to incubation with 10 l plasma sample diluted ( 1:10 ) in wash - buffer ( 0.05% tween20 in pbs ) .
the incubation was performed in multi - well hybridization cassettes ( arrayit corporation ) at rt for 2 hours followed by overnight incubation at 4c . following a short wash ,
the slides were incubated with a cocktail of biotinylated detection antibodies ( anti - human - cd9 , -cd63 and -cd81 , lifespan biosciences , wa , usa ) diluted 1:1,500 in wash - buffer . after a wash ,
a subsequent 30-minute incubation step with cy5-labelled streptavidin ( life technologies ) diluted 1:1,500 in wash - buffer was carried out for detection . prior to scanning ,
the slides were washed in wash - buffer , then in ultrapure / deionized water and finally dried using a microarray high - speed centrifuge ( arrayit corporation ) . scanning and spot detection was performed as previously described ( 15 ) .
microarray printing was performed on a spotbot extreme protein edition microarray printer with a 946mp4 pin ( arrayit corporation , ca , usa ) . as a positive control 100 g / ml of biotinylated human igg
antibodies were printed on epoxy - coated slides ( 75.6 mm25.0 mm ; schott nexterion , de ) , which were left to dry at room temperature overnight prior to further analysis .
a total of 37 anti - human antibodies were used and listed in the following with the corresponding clone , if available : cd146 ( p1h12 ) , flotillin-1 , hb - egf ( 4g10 ) , her3 ( 2f9 ) , her4 ( h4.77.16 ) , hsp90 ( igf1 ) , n - cadherin ( 8c11 ) , p53 ( pab240 ) , tag72 ( 0.n.561 ) , and tsg101 ( 5b7 ) from abcam ( ma , usa ) ; sftpd ( vif11 ) from acris antibodies gmbh ( de ) ; egfr and egfrviii from antibodies-online.com ( de ) ; cd63 ( mem-259 ) from biolegend ( ca , usa ) ; her2 ( 29d8 ) from cell signalling ( ma , usa ) ; cd9 and cd81 from lifespan biosciences ( wa , usa ) ; spa ( 6f10 ) from novus biologicals ( co , usa ) ; cd13 ( 498001 ) , cd14 ( 50040 ) , cd142 ( 323514 ) , cd151 ( 210127 ) , cd206 , cea ( 487609 ) , muc1 ( 604804 ) , tnf ri , tnf rii and tspan8 ( 458811 ) from r&d systems ( mn , usa ) ; areg ( s-13 ) , epcam ( 0.n.277 ) , mucin16 ( x306 ) , ny - eso-1 ( e978 ) and plap ( 8b6 ) from santa cruz biotechnologies ( tx , usa ) ; cd171 from sigma - aldrich ( mo , usa ) ; c - met ( 016 ) and pd - l1 from sino biological inc .
( china ) ; cd163 ( mac2 - 158 ) from trillium diagnostics , llc ( me , usa ) .
all antibodies were diluted with pbs with 5% glycerol and printed in triplicates at 75200 g / ml .
the ev array was prepared as described by jrgensen et al . ( 15 ) with modifications . in short ,
the microarray slides were initially blocked ( 50 mm ethanolamine , 100 mm tris , 0.1% sds , ph 9.0 ) prior to incubation with 10 l plasma sample diluted ( 1:10 ) in wash - buffer ( 0.05% tween20 in pbs ) .
the incubation was performed in multi - well hybridization cassettes ( arrayit corporation ) at rt for 2 hours followed by overnight incubation at 4c . following a short wash ,
the slides were incubated with a cocktail of biotinylated detection antibodies ( anti - human - cd9 , -cd63 and -cd81 , lifespan biosciences , wa , usa ) diluted 1:1,500 in wash - buffer . after a wash ,
a subsequent 30-minute incubation step with cy5-labelled streptavidin ( life technologies ) diluted 1:1,500 in wash - buffer was carried out for detection . prior to scanning ,
the slides were washed in wash - buffer , then in ultrapure / deionized water and finally dried using a microarray high - speed centrifuge ( arrayit corporation ) . scanning and spot detection was performed as previously described ( 15 ) .
graphs and statistics were carried out using graphpad prism ( version 6.04 , graphpad software , inc . ,
san diego , ca , usa ) , excel ( version 2013 , microsoft , redmond , wa , usa ) and roccet ( version 03/17/2014 , tmic , ab , canada ) .
heat maps and hierarchical cluster analysis were produced using genesis ( version 1.7.6 , igb tu graz , graz , austria ) .
random forests was used for classification and monte - carlo cross validation ( mccv ) using balanced subsampling was employed to compute roc curves and to calculate confidence intervals of their area under the curve ( auc ) . in each mccv , two thirds ( 2/3 ) of the samples were used to evaluate the feature importance .
the top important features were then used to build a classification model , which was validated on the 1/3 of the samples that were left out .
five hundred iterations were used to calculate performance of the model . prior to analysis missing value estimation was performed by replacing missing values by  min .  a top 20 of ratios ( based on p - values ) was computed and included into the analysis prior to data transformation and normalization by log2 .
the ev array analysis was performed on 109 plasma samples from patients diagnosed with stage iiia  iv nsclc adenocarcinoma and on 110 plasma samples from patients without a cancer diagnosis ( table i ) .
even though only 10 l unpurified plasma sample was loaded for each patient , sufficient signal for further analysis was obtained from all samples .
selected baseline patient characteristics for the two study groups thirty - seven different antibodies were used for the analysis of plasma exosomes .
in addition , antibodies targeting markers involved in nsclc pathology were optimised for the ev array analysis . for this purpose ,
members of the egf system and selected markers of inflammation ( cd163 , cd206 ) ( 16 , 17 ) were selected ( table ii ) .
description of the markers selected for the ev array together with the outcome ( p - value ) of the non - parametric t - test comparing the log2 transformed data from the control group and the cancer group ( supplementary fig .
1a and b ) unless otherwise stated the marker level was higher in the control group . elevated in cancer patients . the univariate analysis of each marker using log2-transformed signal - values revealed generally a significantly higher marker - expression in control subjects compared to lung cancer subjects ( table ii and supplementary fig . 1a and b ) .
hierarchical cluster analysis for each of the groups depicted a high variety of expression among the individuals , which was also seen in a previously study of healthy donors ( 15 ) .
the hierarchical cluster analysis also showed co - variance among certain markers , for example cd9 , cd81 and cd151 ( marked with a box ) and egfrviii , flotilin-1 , her4 , n - cadherin and cd163 ( marked with a dashed box ) ( fig .
the control group had a significantly higher ( p<0.05 ) cd9 , cd81 and cd151 expression when looking at log2-transformed data ( fig .
for the co - varying group consisting of egfrviii , flotilin-1 , her4 , n - cadherin and cd163 , all markers were higher expressed in the control group , as well as when looking at log2-transformed data ( fig .
2b ) . as the cancerous patients showed a tendency to have a lower amount of positive antigens over all ( fig .
3a ) , a normalisation ( calculation of the percentage ) of the signals for each antigen to the total signal for each individual patient was performed ( fig .
3a ) . taken together with the fact that the ev array signal detection relies on cd9- , cd63- and/or cd81-antibody binding , the relative expression of these markers was examined further ( pie charts in fig .
3b ) , the cd9 , cd63 and cd81 percentage was significantly higher ( p<0.05 ) in the cancer group .
two groups of markers show co - variance both in the control group and in the cancer group ( marked with boxes ) .
b ) heat map illustration of all markers in the cancer group . ev array signal intensities for selected antigens .
a ) the ev array signal intensities for the exosomal markers cd9 , cd63 and cd81 displayed in box plots .
the co - variation of the signal intensities across the patient samples can be seen in fig . 1 and supplementary fig
b ) box plot of a group of antigens ( flotilin-1 , her4 , egfrviii , n - cadherin and cd163 ) showing a high degree of co - variation ( see fig .
* p<0.05 ; * * * p<0.001 ; * * * * p<0.0001 ; ns = not significant .
a ) the signals for all analytes were summed for each individual patient and plotted ; controls indicated with green and cancer with red . for each individual patient
the pie charts illustrate an example of the normalised data for a patient in each group with a total amount of signal of ~40 .
b ) and c ) box plot of the relative expression of markers from figure 2a and b in percentage ( in relation to the total sum of exosomal signal ) . * p < 0.05 ; * * p<0.01 ; * * * * p<0.0001 ; ns = not significant .
normalised data for the co - varying group egfrviii , flotilin-1 , her4 , n - cadherin and cd163 was also investigated . for these markers
none of the individual markers alone gained an auc indicating a diagnostic potential ( fig .
however , the diagnostic potential of a panel of markers was evaluated using multivariate analysis . for the multivariate analysis
the univariate analysis of the normalised data for cd9 , cd63 and cd81 ( fig .
therefore , these normalised values were accepted in the multivariate analysis together with the individual relations between the analysed markers . 

figure 4 shows the roc curves for the 3 , 5 , 10 and 30 best rated marker combinations .
the 30-marker model has the largest auc of 0.83 , ci : 0.770.90 , a sensitivity of 0.75 and a specificity of 0.76 .
was observed when expanding the number of markers included beyond 30 ( data not shown ) .
the area under curve ( auc ) for top 3- , 5- , 10- , and 30-marker panels are given together with the 95% confidence interval .
the importance of each marker or their individual relations for the model was assessed ( fig .
the most important marker was found to be cd81 followed by the relation between tag72 and the normalised values of cd63 .
the top 10 
 ranking of the importance did not change between the models including 3 , 5 , 10 or 30 markers .
the mean average importance to the classification model using the 30-marker panel illustrated in fig . 4 for each of the analysed exosomal antigens , the normalised values ( indicated by  *  ) and their internal relations ( indicated by  /  ) .
the top 10 ranking did not change between the models including 3- , 5- , 10- or 30-markers and the markers included in each model are visualised by the coloured lines .
the ev array analysis was performed on 109 plasma samples from patients diagnosed with stage iiia  iv nsclc adenocarcinoma and on 110 plasma samples from patients without a cancer diagnosis ( table i ) .
even though only 10 l unpurified plasma sample was loaded for each patient , sufficient signal for further analysis was obtained from all samples .
selected baseline patient characteristics for the two study groups thirty - seven different antibodies were used for the analysis of plasma exosomes .
in addition , antibodies targeting markers involved in nsclc pathology were optimised for the ev array analysis . for this purpose ,
members of the egf system and selected markers of inflammation ( cd163 , cd206 ) ( 16 , 17 ) were selected ( table ii ) .
description of the markers selected for the ev array together with the outcome ( p - value ) of the non - parametric t - test comparing the log2 transformed data from the control group and the cancer group ( supplementary fig .
1a and b ) unless otherwise stated the marker level was higher in the control group . elevated in cancer patients .
the univariate analysis of each marker using log2-transformed signal - values revealed generally a significantly higher marker - expression in control subjects compared to lung cancer subjects ( table ii and supplementary fig .
hierarchical cluster analysis for each of the groups depicted a high variety of expression among the individuals , which was also seen in a previously study of healthy donors ( 15 ) .
the hierarchical cluster analysis also showed co - variance among certain markers , for example cd9 , cd81 and cd151 ( marked with a box ) and egfrviii , flotilin-1 , her4 , n - cadherin and cd163 ( marked with a dashed box ) ( fig .
the control group had a significantly higher ( p<0.05 ) cd9 , cd81 and cd151 expression when looking at log2-transformed data ( fig .
for the co - varying group consisting of egfrviii , flotilin-1 , her4 , n - cadherin and cd163 , all markers were higher expressed in the control group , as well as when looking at log2-transformed data ( fig .
2b ) . as the cancerous patients showed a tendency to have a lower amount of positive antigens over all ( fig .
3a ) , a normalisation ( calculation of the percentage ) of the signals for each antigen to the total signal for each individual patient was performed ( fig .
3a ) . taken together with the fact that the ev array signal detection relies on cd9- , cd63- and/or cd81-antibody binding , the relative expression of these markers was examined further ( pie charts in fig .
3b ) , the cd9 , cd63 and cd81 percentage was significantly higher ( p<0.05 ) in the cancer group .
two groups of markers show co - variance both in the control group and in the cancer group ( marked with boxes ) .
b ) heat map illustration of all markers in the cancer group . ev array signal intensities for selected antigens .
a ) the ev array signal intensities for the exosomal markers cd9 , cd63 and cd81 displayed in box plots .
the co - variation of the signal intensities across the patient samples can be seen in fig . 1 and supplementary fig
b ) box plot of a group of antigens ( flotilin-1 , her4 , egfrviii , n - cadherin and cd163 ) showing a high degree of co - variation ( see fig .
* p<0.05 ; * * * p<0.001 ; * * * * p<0.0001 ; ns = not significant .
a ) the signals for all analytes were summed for each individual patient and plotted ; controls indicated with green and cancer with red . for each individual patient
the pie charts illustrate an example of the normalised data for a patient in each group with a total amount of signal of ~40 .
b ) and c ) box plot of the relative expression of markers from figure 2a and b in percentage ( in relation to the total sum of exosomal signal ) . * p < 0.05 ; * * p<0.01 ; * * * * p<0.0001 ; ns = not significant .
normalised data for the co - varying group egfrviii , flotilin-1 , her4 , n - cadherin and cd163 was also investigated . for these markers
none of the individual markers alone gained an auc indicating a diagnostic potential ( fig .
however , the diagnostic potential of a panel of markers was evaluated using multivariate analysis . for the multivariate analysis
the univariate analysis of the normalised data for cd9 , cd63 and cd81 ( fig .
therefore , these normalised values were accepted in the multivariate analysis together with the individual relations between the analysed markers . 

figure 4 shows the roc curves for the 3 , 5 , 10 and 30 best rated marker combinations .
the 30-marker model has the largest auc of 0.83 , ci : 0.770.90 , a sensitivity of 0.75 and a specificity of 0.76 .
was observed when expanding the number of markers included beyond 30 ( data not shown ) .
the area under curve ( auc ) for top 3- , 5- , 10- , and 30-marker panels are given together with the 95% confidence interval .
the importance of each marker or their individual relations for the model was assessed ( fig .
the most important marker was found to be cd81 followed by the relation between tag72 and the normalised values of cd63 .
the top 10 
 ranking of the importance did not change between the models including 3 , 5 , 10 or 30 markers .
the mean average importance to the classification model using the 30-marker panel illustrated in fig . 4 for each of the analysed exosomal antigens , the normalised values ( indicated by  *  ) and their internal relations ( indicated by  /  )
. the top 10 ranking did not change between the models including 3- , 5- , 10- or 30-markers and the markers included in each model are visualised by the coloured lines .
so far , only a limited number of studies have examined the diagnostic potential in relation to nsclc . until now
, the studies have been limited to investigating exosomes captured by antibodies towards a single antigen , for example , epcam or cd9 ( 9 , 13 ) . in spite of these constraints ,
promising results were obtained through these studies and a subset of exosomal mirnas and the exosomal - bound cd91-protein was revealed as promising biomarkers .
the present study examined plasma exosomes from 219 patients , initially suspected of having lung cancer , from a multi - marker aspect .
investigation of 37 exosome capturing antibodies and the use of a cocktail of 3 different detecting antibodies against known exosome markers gave comprehensive information paving the way for multivariate analysis .
our multivariate ev array analysis results in a 30-marker panel classifying 75% of the patients correctly .
this result is superior to exosomal cd91-expression alone ( 61% ) as recently reported by ueda et al .
our data suggests that the ev array analysis could be a valuable complementary tool in diagnosing nsclc .
the univariate t - test analysis revealed higher exosomal marker abundance on the control patient exosomes compared to the cancer patient exosomes .
this is in concordance with the recent study of serum samples from nsclc patients and healthy blood donors performed by ueda et al .
. the decreased overall signal could be due to less compact exosomes in lung cancer patients , for example , fewer proteins per exosome .
presumably , this would result in fewer lung cancer exosomes being captured by the antibodies .
previously , differences in the protein load of the exosomes of melanoma patients have also been shown by peinado et al . to have an influence on the survival rate ( 61 ) .
another possibility is that exosomes from nsclc patients express another subset of markers than the ones targeted in this analysis or that the exosomes carrying the cancerous markers are cleared by an activated immune system in the cancerous patients ( 62 ) .
additional markers could be explored to optimise the ev array and other detection antibodies could be evaluated for further investigation of this phenomenon .
cd9 and cd81 remain well - described exosomal markers and may be expected to co - vary .
interestingly , cd63 , another well - known marker , does not follow this pattern , which has also been described previously ( 15 , 63 ) .
the difference in expression of cd63 in relation to cd9 and cd81 highlights the advantage of detecting the signal with a cocktail of cd9 , cd63 and cd81 antibodies in order to achieve sufficient information .
egfrviii together with flotilin-1 , her4 , n - cadherin and cd163 constituted an additional co - varying group ( fig .
nevertheless , egfr truncation is not normally observed in nsclc and neither in normal tissue ( 64 , 65 ) .
an obvious first thought to explain the level of egfrviii signal in both control and cancer patients could point to an unspecific antibody binding to normal egfr .
however , since the patterns of these 2 markers ( egfr and egfrviii ) vary from each other this explanation does not seem to be the case .
this phenomenon should be explored further in order to assess the biological bearing of this truncated receptor . as seen from the ranking of the 30 markers ,
their contribution differs in importance to the model . increasing the average importance with new and more efficient markers
furthermore , as no increase in diagnostic potential was seen when including more markers , it is essential to substitute weak markers with better candidate markers in order to increase the average importance and thereby provide an even more valuable tool . in conclusion ,
the ev array constitutes a broad spectra tool offering exhaustive information on exosomal content with a minimal of sample requirement , which can be optimized and adjusted to fit individual sets of samples .
the authors have not received any funding or benefits from industry or elsewhere to conduct this study .
 backgroundlung cancer is one of the leading causes of cancer - related death . at the time of diagnosis , 
 more than half of the patients will have disseminated disease and , yet , diagnosing can be challenging . 
 new methods are desired to improve the diagnostic work - up . 
 exosomes are cell - derived vesicles displaying various proteins on their membrane surfaces . in addition , they are readily available in blood samples where they constitute potential biomarkers of human diseases , such as cancer . here 
 , we examine the potential of distinguishing non - small cell lung carcinoma ( nsclc ) patients from control subjects based on the differential display of exosomal protein markers.methodsplasma was isolated from 109 nsclc patients with advanced stage ( iiia  iv ) disease and 110 matched control subjects initially suspected of having cancer , but diagnosed to be cancer free . the extracellular vesicle array ( ev array ) 
 was used to phenotype exosomes directly from the plasma samples . 
 the array contained 37 antibodies targeting lung cancer - related proteins and was used to capture exosomes , which were visualised with a cocktail of biotin - conjugated cd9 , cd63 and cd81 antibodies.resultsthe ev array analysis was capable of detecting and phenotyping exosomes in all samples from only 10 l of unpurified plasma . 
 multivariate analysis using the random forests method produced a combined 30-marker model separating the two patient groups with an area under the curve of 0.83 , ci : 0.770.90 . 
 the 30-marker model has a sensitivity of 0.75 and a specificity of 0.76 , and it classifies patients with 75.3% accuracy.conclusionthe ev array technique is a simple , minimal - invasive tool with potential to identify lung cancer patients . 
Material and methods
Patient samples
Phenotyping of extracellular vesicles using the EV Array
Production of microarrays
Antibodies for production of the EV Array
Catching and visualization of exosomes
Data analysis
Multivariant ROC analysis
Results
Exosomal markers can be detected in small amounts of plasma
Univariate analysis of exosomal markers
Multivariate exosomal marker analysis as a diagnostic tool in NSCLC
Discussion
Supplementary Material
Conflict of interest and funding