Unlocking biomarker discovery: Large scale application of aptamer proteomic technology for early detection of lung cancer

Lung cancer is the leading cause of cancer deaths, because ~84% of cases are diagnosed at an advanced stage. Worldwide in 2008, ~1.5 million people were diagnosed and ~1.3 million died – a survival rate unchanged since 1960. However, patients diagnosed at an early stage and have surgery experience an 86% overall 5-year survival. New diagnostics are therefore needed to identify lung cancer at this stage. Here we present the first large scale clinical use of aptamers to discover blood protein biomarkers in disease with our breakthrough proteomic technology. This multi-center case-control study was conducted in archived samples from 1,326 subjects from four independent studies of non-small cell lung cancer (NSCLC) in long-term tobacco-exposed populations. We measured >800 proteins in 15uL of serum, identified 44 candidate biomarkers, and developed a 12-protein panel that distinguished NSCLC from controls with 91% sensitivity and 84% specificity in a training set and 89% sensitivity and 83% specificity in a blinded, independent verification set. Performance was similar for early and late stage NSCLC. This is a significant advance in proteomics in an area of high clinical need

Mechanisms of sodium-glucose cotransporter-2 inhibition: insights from large-scale proteomics

OBJECTIVE To assess the effects of empagliflozin, a selective sodium–glucose cotransporter 2 (SGLT2) inhibitor, on broad biological systems through proteomics. RESEARCH DESIGN AND METHODS Aptamer-based proteomics was used to quantify 3,713 proteins in 144 paired plasma samples obtained from 72 participants across the spectrum of glucose tolerance before and after 4 weeks of empagliflozin 25 mg/day. The biology of the plasma proteins significantly changed by empagliflozin (at false discovery rate–corrected P < 0.05) was discerned through Ingenuity Pathway Analysis. RESULTS Empagliflozin significantly affected levels of 43 proteins, 6 related to cardiomyocyte function (fatty acid–binding protein 3 and 4 [FABPA], neurotrophic receptor tyrosine kinase, renin, thrombospondin 4, and leptin receptor), 5 to iron handling (ferritin heavy chain 1, transferrin receptor protein 1, neogenin, growth differentiation factor 2 [GDF2], and β2-microglobulin), and 1 to sphingosine/ceramide metabolism (neutral ceramidase), a known pathway of cardiovascular disease. Among the protein changes achieving the strongest statistical significance, insulin-like binding factor protein-1 (IGFBP-1), transgelin-2, FABPA, GDF15, and sulphydryl oxidase 2 precursor were increased, while ferritin, thrombospondin 3, and Rearranged during Transfection (RET) were decreased by empagliflozin administration. CONCLUSIONS SGLT2 inhibition is associated, directly or indirectly, with multiple biological effects, including changes in markers of cardiomyocyte contraction/relaxation, iron handling, and other metabolic and renal targets. The most significant differences were detected in protein species (GDF15, ferritin, IGFBP-1, and FABP) potentially related to the clinical and metabolic changes that were actually measured in the same patients. These novel results may inform further studies using targeted proteomics and a prospective design

HFrEF subphenotypes based on 4210 repeatedly measured circulating proteins are driven by different biological mechanisms

Background: HFrEF is a heterogenous condition with high mortality. We used serial assessments of 4210 circulating proteins to identify distinct novel protein-based HFrEF subphenotypes and to investigate underlying dynamic biological mechanisms. Herewith we aimed to gain pathophysiological insights and fuel opportunities for personalised treatment. Methods: In 382 patients, we performed trimonthly blood sampling during a median follow-up of 2.1 [IQR:1.1–2.6] years. We selected all baseline samples and two samples closest to the primary endpoint (PEP; composite of cardiovascular mortality, HF hospitalization, LVAD implantation, and heart transplantation) or censoring, and applied an aptamer-based multiplex proteomic approach. Using unsupervised machine learning methods, we derived clusters from 4210 repeatedly measured proteomic biomarkers. Sets of proteins that drove cluster allocation were analysed via an enrichment analysis. Differences in clinical characteristics and PEP occurrence were evaluated. Findings: We identified four subphenotypes with different protein profiles, prognosis and clinical characteristics, including age (median [IQR] for subphenotypes 1–4, respectively:70 [64, 76], 68 [60, 79], 57 [47, 65], 59 [56, 66]years), EF (30 [26, 36], 26 [20, 38], 26 [22, 32], 33 [28, 37]%), and chronic renal failure (45%, 65%, 36%, 37%). Subphenotype allocation was driven by subsets of proteins associated with various biological functions, such as oxidative stress, inflammation and extracellular matrix organisation. Clinical characteristics of the subphenotypes were aligned with these associations. Subphenotypes 2 and 3 had the worst prognosis compared to subphenotype 1 (adjHR (95%CI):3.43 (1.76–6.69), and 2.88 (1.37–6.03), respectively). Interpretation: Four circulating-protein based subphenotypes are present in HFrEF, which are driven by varying combinations of protein subsets, and have different clinical characteristics and prognosis. Clinical Trial Registration: ClinicalTrials.gov Identifier: NCT01851538 https://clinicaltrials.gov/ct2/show/NCT01851538. Funding: EU/ EFPIA IMI2JU BigData@Heart grant n° 116074, Jaap Schouten Foundation and Noordwest Academie.</p

Aptamer-based multiplexed proteomic technology for biomarker discovery

Interrogation of the human proteome in a highly multiplexed and efficient manner remains a coveted and challenging goal in biology. We present a new aptamer-based proteomic technology for biomarker discovery capable of simultaneously measuring thousands of proteins from small sample volumes (15 [mu]L of serum or plasma). Our current assay allows us to measure ~800 proteins with very low limits of detection (1 pM average), 7 logs of overall dynamic range, and 5% average coefficient of variation. This technology is enabled by a new generation of aptamers that contain chemically modified nucleotides, which greatly expand the physicochemical diversity of the large randomized nucleic acid libraries from which the aptamers are selected. Proteins in complex matrices such as plasma are measured with a process that transforms a signature of protein concentrations into a corresponding DNA aptamer concentration signature, which is then quantified with a DNA microarray. In essence, our assay takes advantage of the dual nature of aptamers as both folded binding entities with defined shapes and unique sequences recognizable by specific hybridization probes. To demonstrate the utility of our proteomics biomarker discovery technology, we applied it to a clinical study of chronic kidney disease (CKD). We identified two well known CKD biomarkers as well as an additional 58 potential CKD biomarkers. These results demonstrate the potential utility of our technology to discover unique protein signatures characteristic of various disease states. More generally, we describe a versatile and powerful tool that allows large-scale comparison of proteome profiles among discrete populations. This unbiased and highly multiplexed search engine will enable the discovery of novel biomarkers in a manner that is unencumbered by our incomplete knowledge of biology, thereby helping to advance the next generation of evidence-based medicine

Biomarkers for staging fibrosis and non-alcoholic steatohepatitis in non-alcoholic fatty liver disease (the LITMUS project): a comparative diagnostic accuracy study.

BACKGROUND The reference standard for detecting non-alcoholic steatohepatitis (NASH) and staging fibrosis-liver biopsy-is invasive and resource intensive. Non-invasive biomarkers are urgently needed, but few studies have compared these biomarkers in a single cohort. As part of the Liver Investigation: Testing Marker Utility in Steatohepatitis (LITMUS) project, we aimed to evaluate the diagnostic accuracy of 17 biomarkers and multimarker scores in detecting NASH and clinically significant fibrosis in patients with non-alcoholic fatty liver disease (NAFLD) and identify their optimal cutoffs as screening tests in clinical trial recruitment. METHODS This was a comparative diagnostic accuracy study in people with biopsy-confirmed NAFLD from 13 countries across Europe, recruited between Jan 6, 2010, and Dec 29, 2017, from the LITMUS metacohort of the prospective European NAFLD Registry. Adults (aged ≥18 years) with paired liver biopsy and serum samples were eligible; those with excessive alcohol consumption or evidence of other chronic liver diseases were excluded. The diagnostic accuracy of the biomarkers was expressed as the area under the receiver operating characteristic curve (AUC) with liver histology as the reference standard and compared with the Fibrosis-4 index for liver fibrosis (FIB-4) in the same subgroup. Target conditions were the presence of NASH with clinically significant fibrosis (ie, at-risk NASH; NAFLD Activity Score ≥4 and F≥2) or the presence of advanced fibrosis (F≥3), analysed in all participants with complete data. We identified thres holds for each biomarker for reducing the number of biopsy-based screen failures when recruiting people with both NASH and clinically significant fibrosis for future trials. FINDINGS Of 1430 participants with NAFLD in the LITMUS metacohort with serum samples, 966 (403 women and 563 men) were included after all exclusion criteria had been applied. 335 (35%) of 966 participants had biopsy-confirmed NASH and clinically significant fibrosis and 271 (28%) had advanced fibrosis. For people with NASH and clinically significant fibrosis, no single biomarker or multimarker score significantly reached the predefined AUC 0·80 acceptability threshold (AUCs ranging from 0·61 [95% CI 0·54-0·67] for FibroScan controlled attenuation parameter to 0·81 [0·75-0·86] for SomaSignal), with accuracy mostly similar to FIB-4. Regarding detection of advanced fibrosis, SomaSignal (AUC 0·90 [95% CI 0·86-0·94]), ADAPT (0·85 [0·81-0·89]), and FibroScan liver stiffness measurement (0·83 [0·80-0·86]) reached acceptable accuracy. With 11 of 17 markers, histological screen failure rates could be reduced to 33% in trials if only people who were marker positive had a biopsy for evaluating eligibility. The best screening performance for NASH and clinically significant fibrosis was observed for SomaSignal (number needed to test [NNT] to find one true positive was four [95% CI 4-5]), then ADAPT (six [5-7]), MACK-3 (seven [6-8]), and PRO-C3 (nine [7-11]). INTERPRETATION None of the single markers or multimarker scores achieved the predefined acceptable AUC for replacing biopsy in detecting people with both NASH and clinically significant fibrosis. However, several biomarkers could be applied in a prescreening strategy in clinical trial recruitment. The performance of promising markers will be further evaluated in the ongoing prospective LITMUS study cohort. FUNDING The Innovative Medicines Initiative 2 Joint Undertaking

Protein Signature of Lung Cancer Tissues

Lung cancer remains the most common cause of cancer-related mortality. We applied a highly multiplexed proteomic technology (SOMAscan) to compare protein expression signatures of non small-cell lung cancer (NSCLC) tissues with healthy adjacent and distant tissues from surgical resections. In this first report of SOMAscan applied to tissues, we highlight 36 proteins that exhibit the largest expression differences between matched tumor and non-tumor tissues. The concentrations of twenty proteins increased and sixteen decreased in tumor tissue, thirteen of which are novel for NSCLC. NSCLC tissue biomarkers identified here overlap with a core set identified in a large serum-based NSCLC study with SOMAscan. We show that large-scale comparative analysis of protein expression can be used to develop novel histochemical probes. As expected, relative differences in protein expression are greater in tissues than in serum. The combined results from tissue and serum present the most extensive view to date of the complex changes in NSCLC protein expression and provide important implications for diagnosis and treatment

ProteinSeq: High-Performance Proteomic Analyses by Proximity Ligation and Next Generation Sequencing

Despite intense interest, methods that provide enhanced sensitivity and specificity in parallel measurements of candidate protein biomarkers in numerous samples have been lacking. We present herein a multiplex proximity ligation assay with readout via realtime PCR or DNA sequencing (ProteinSeq). We demonstrate improved sensitivity over conventional sandwich assays for simultaneous analysis of sets of 35 proteins in 5 µl of blood plasma. Importantly, we observe a minimal tendency to increased background with multiplexing, compared to a sandwich assay, suggesting that higher levels of multiplexing are possible. We used ProteinSeq to analyze proteins in plasma samples from cardiovascular disease (CVD) patient cohorts and matched controls. Three proteins, namely P-selectin, Cystatin-B and Kallikrein-6, were identified as putative diagnostic biomarkers for CVD. The latter two have not been previously reported in the literature and their potential roles must be validated in larger patient cohorts. We conclude that ProteinSeq is promising for screening large numbers of proteins and samples while the technology can provide a much-needed platform for validation of diagnostic markers in biobank samples and in clinical use

Global age-sex-specific fertility, mortality, healthy life expectancy (HALE), and population estimates in 204 countries and territories, 1950-2019 : a comprehensive demographic analysis for the Global Burden of Disease Study 2019

Background: Accurate and up-to-date assessment of demographic metrics is crucial for understanding a wide range of social, economic, and public health issues that affect populations worldwide. The Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) 2019 produced updated and comprehensive demographic assessments of the key indicators of fertility, mortality, migration, and population for 204 countries and territories and selected subnational locations from 1950 to 2019. Methods: 8078 country-years of vital registration and sample registration data, 938 surveys, 349 censuses, and 238 other sources were identified and used to estimate age-specific fertility. Spatiotemporal Gaussian process regression (ST-GPR) was used to generate age-specific fertility rates for 5-year age groups between ages 15 and 49 years. With extensions to age groups 10–14 and 50–54 years, the total fertility rate (TFR) was then aggregated using the estimated age-specific fertility between ages 10 and 54 years. 7417 sources were used for under-5 mortality estimation and 7355 for adult mortality. ST-GPR was used to synthesise data sources after correction for known biases. Adult mortality was measured as the probability of death between ages 15 and 60 years based on vital registration, sample registration, and sibling histories, and was also estimated using ST-GPR. HIV-free life tables were then estimated using estimates of under-5 and adult mortality rates using a relational model life table system created for GBD, which closely tracks observed age-specific mortality rates from complete vital registration when available. Independent estimates of HIV-specific mortality generated by an epidemiological analysis of HIV prevalence surveys and antenatal clinic serosurveillance and other sources were incorporated into the estimates in countries with large epidemics. Annual and single-year age estimates of net migration and population for each country and territory were generated using a Bayesian hierarchical cohort component model that analysed estimated age-specific fertility and mortality rates along with 1250 censuses and 747 population registry years. We classified location-years into seven categories on the basis of the natural rate of increase in population (calculated by subtracting the crude death rate from the crude birth rate) and the net migration rate. We computed healthy life expectancy (HALE) using years lived with disability (YLDs) per capita, life tables, and standard demographic methods. Uncertainty was propagated throughout the demographic estimation process, including fertility, mortality, and population, with 1000 draw-level estimates produced for each metric. Findings: The global TFR decreased from 2·72 (95% uncertainty interval [UI] 2·66–2·79) in 2000 to 2·31 (2·17–2·46) in 2019. Global annual livebirths increased from 134·5 million (131·5–137·8) in 2000 to a peak of 139·6 million (133·0–146·9) in 2016. Global livebirths then declined to 135·3 million (127·2–144·1) in 2019. Of the 204 countries and territories included in this study, in 2019, 102 had a TFR lower than 2·1, which is considered a good approximation of replacement-level fertility. All countries in sub-Saharan Africa had TFRs above replacement level in 2019 and accounted for 27·1% (95% UI 26·4–27·8) of global livebirths. Global life expectancy at birth increased from 67·2 years (95% UI 66·8–67·6) in 2000 to 73·5 years (72·8–74·3) in 2019. The total number of deaths increased from 50·7 million (49·5–51·9) in 2000 to 56·5 million (53·7–59·2) in 2019. Under-5 deaths declined from 9·6 million (9·1–10·3) in 2000 to 5·0 million (4·3–6·0) in 2019. Global population increased by 25·7%, from 6·2 billion (6·0–6·3) in 2000 to 7·7 billion (7·5–8·0) in 2019. In 2019, 34 countries had negative natural rates of increase; in 17 of these, the population declined because immigration was not sufficient to counteract the negative rate of decline. Globally, HALE increased from 58·6 years (56·1–60·8) in 2000 to 63·5 years (60·8–66·1) in 2019. HALE increased in 202 of 204 countries and territories between 2000 and 2019

Global burden of 369 diseases and injuries in 204 countries and territories, 1990–2019: a systematic analysis for the Global Burden of Disease Study 2019

Background: In an era of shifting global agendas and expanded emphasis on non-communicable diseases and injuries along with communicable diseases, sound evidence on trends by cause at the national level is essential. The Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) provides a systematic scientific assessment of published, publicly available, and contributed data on incidence, prevalence, and mortality for a mutually exclusive and collectively exhaustive list of diseases and injuries. Methods: GBD estimates incidence, prevalence, mortality, years of life lost (YLLs), years lived with disability (YLDs), and disability-adjusted life-years (DALYs) due to 369 diseases and injuries, for two sexes, and for 204 countries and territories. Input data were extracted from censuses, household surveys, civil registration and vital statistics, disease registries, health service use, air pollution monitors, satellite imaging, disease notifications, and other sources. Cause-specific death rates and cause fractions were calculated using the Cause of Death Ensemble model and spatiotemporal Gaussian process regression. Cause-specific deaths were adjusted to match the total all-cause deaths calculated as part of the GBD population, fertility, and mortality estimates. Deaths were multiplied by standard life expectancy at each age to calculate YLLs. A Bayesian meta-regression modelling tool, DisMod-MR 2.1, was used to ensure consistency between incidence, prevalence, remission, excess mortality, and cause-specific mortality for most causes. Prevalence estimates were multiplied by disability weights for mutually exclusive sequelae of diseases and injuries to calculate YLDs. We considered results in the context of the Socio-demographic Index (SDI), a composite indicator of income per capita, years of schooling, and fertility rate in females younger than 25 years. Uncertainty intervals (UIs) were generated for every metric using the 25th and 975th ordered 1000 draw values of the posterior distribution. Findings: Global health has steadily improved over the past 30 years as measured by age-standardised DALY rates. After taking into account population growth and ageing, the absolute number of DALYs has remained stable. Since 2010, the pace of decline in global age-standardised DALY rates has accelerated in age groups younger than 50 years compared with the 1990–2010 time period, with the greatest annualised rate of decline occurring in the 0–9-year age group. Six infectious diseases were among the top ten causes of DALYs in children younger than 10 years in 2019: lower respiratory infections (ranked second), diarrhoeal diseases (third), malaria (fifth), meningitis (sixth), whooping cough (ninth), and sexually transmitted infections (which, in this age group, is fully accounted for by congenital syphilis; ranked tenth). In adolescents aged 10–24 years, three injury causes were among the top causes of DALYs: road injuries (ranked first), self-harm (third), and interpersonal violence (fifth). Five of the causes that were in the top ten for ages 10–24 years were also in the top ten in the 25–49-year age group: road injuries (ranked first), HIV/AIDS (second), low back pain (fourth), headache disorders (fifth), and depressive disorders (sixth). In 2019, ischaemic heart disease and stroke were the top-ranked causes of DALYs in both the 50–74-year and 75-years-and-older age groups. Since 1990, there has been a marked shift towards a greater proportion of burden due to YLDs from non-communicable diseases and injuries. In 2019, there were 11 countries where non-communicable disease and injury YLDs constituted more than half of all disease burden. Decreases in age-standardised DALY rates have accelerated over the past decade in countries at the lower end of the SDI range, while improvements have started to stagnate or even reverse in countries with higher SDI. Interpretation: As disability becomes an increasingly large component of disease burden and a larger component of health expenditure, greater research and developm nt investment is needed to identify new, more effective intervention strategies. With a rapidly ageing global population, the demands on health services to deal with disabling outcomes, which increase with age, will require policy makers to anticipate these changes. The mix of universal and more geographically specific influences on health reinforces the need for regular reporting on population health in detail and by underlying cause to help decision makers to identify success stories of disease control to emulate, as well as opportunities to improve. Funding: Bill & Melinda Gates Foundation. © 2020 The Author(s). Published by Elsevier Ltd. This is an Open Access article under the CC BY 4.0 licens

Increased serum miR-193a-5p during non-alcoholic fatty liver disease progression: diagnostic and mechanistic relevance

Background & Aims: Serum microRNAs (miRNAs) levels are known to change in non-alcoholic fatty liver disease (NAFLD) and may serve as useful biomarkers. This study aimed to profile miRNAs comprehensively at all NAFLD stages.Methods: We profiled 2,083 serum miRNAs in a discovery cohort (183 NAFLD cases representing the complete NAFLD spectrum and 10 population controls). MiRNA libraries generated by HTG EdgeSeq were sequenced by Illumina NextSeq. Selected serum miRNAs were profiled in 372 additional NAFLD cases and 15 population controls by quantitative reverse transcriptase-polymerase chain reaction.Results: Levels of 275 miRNAs differed between cases and population controls. Fewer differences were seen within individual NAFLD stages but miR-193a-5p consistently the showed increased levels in all comparisons. Relative to NAFL/NASH with mild fibrosis (stage 0/1), three miRNAs (miR-193a-5p, miR-378d and miR378d) were increased in cases with NASH and clinically significant fibrosis (stage 2-4), seven (miR193a-5p, miR-378d, miR-378e, miR-320b, c, d & e) increased in cases with NAFLD Activity Score (NAS) 5-8 compared with lower NAS, and three (miR-193a-5p, miR-378d, miR-378e) increased but one (miR-19b-3p) decreased in steatosis, activity, and fibrosis "activity" (SAF-A) score 2-4 compared with lower SAF-A. The significant findings for miR-193a-5p were replicated in the additional NAFLD cohort. Studies in Hep G2 cells showed that following palmitic acid treatment, miR-193a-5p expression decreased significantly. Gene targets for miR-193a-5p were investigated in liver RNAseq data for a case subgroup (n=80); liver GPX8 levels correlated positively with serum miR-193a-5p. Conclusions: Serum miR-193a-5p levels correlate strongly with NAFLD activity grade and fibrosis stage. MiR-193a-5p may have a role in the hepatic response to oxidative stress and is a potential clinically tractable circulating biomarker for progressive NAFLD