Development of circulating microRNA-based biomarkers for medical decision-making: a friendly reminder of what should NOT be done.

Circulating cell-free microRNAs (miRNAs) represent a major reservoir for biomarker discovery. Unfortunately, their implementation in clinical practice is limited due to a profound lack of reproducibility. The great technical variability linked to major pre-analytical and analytical caveats makes the interpretation of circulating cell-free miRNA data challenging and leads to inconsistent findings. Additional efforts directed to standardization are fundamental. Several well-established protocols are currently used by independent groups worldwide. Nonetheless, there are some specific aspects in specimen collection and processing, sample handling, miRNA quantification, and data analysis that should be considered to ensure reproducibility of results. Here, we have addressed this challenge using an alternative approach. We have highlighted and discussed common pitfalls that negatively impact the robustness of circulating miRNA quantification and their application for clinical decision-making. Furthermore, we provide a checklist usable by investigators to facilitate and ensure the control of the whole miRNA quantification and analytical process. We expect that these recommendations improve the reproducibility of findings, and ultimately, facilitate the incorporation of circulating miRNA profiles into clinical practice as the next generation of disease biomarkers.

Circulating miR-133a-3p defines a low-risk subphenotype in patients with heart failure and central sleep apnea: a decision tree machine learning approach.

BACKGROUND: Patients with heart failure with reduced ejection fraction (HFrEF) and central sleep apnea (CSA) are at a very high risk of fatal outcomes. OBJECTIVE: To test whether the circulating miRNome provides additional information for risk stratification on top of clinical predictors in patients with HFrEF and CSA. METHODS: The study included patients with HFrEF and CSA from the SERVE-HF trial. A three-step protocol was applied: microRNA (miRNA) screening (n = 20), technical validation (n = 60), and biological validation (n = 587). The primary outcome was either death from any cause, lifesaving cardiovascular intervention, or unplanned hospitalization for worsening of heart failure, whatever occurred first. MiRNA quantification was performed in plasma samples using miRNA sequencing and RT-qPCR. RESULTS: Circulating miR-133a-3p levels were inversely associated with the primary study outcome. Nonetheless, miR-133a-3p did not improve a previously established clinical prognostic model in terms of discrimination or reclassification. A customized regression tree model constructed using the Classification and Regression Tree (CART) algorithm identified eight patient subphenotypes with specific risk patterns based on clinical and molecular characteristics. MiR-133a-3p entered the regression tree defining the group at the lowest risk; patients with log(NT-proBNP) ≤ 6 pg/mL (miR-133a-3p levels above 1.5 arbitrary units). The overall predictive capacity of suffering the event was highly stable over the follow-up (from 0.735 to 0.767). CONCLUSIONS: The combination of clinical information, circulating miRNAs, and decision tree learning allows the identification of specific risk subphenotypes in patients with HFrEF and CSA.

Altered insulin secretion dynamics relate to oxidative stress and inflammasome activation in children with obesity and insulin resistance.

BACKGROUND: Insulin resistance (IR) is considered the main driver of obesity related metabolic complications, and is related to oxidative stress and inflammation, which in turn promote each other. There is currently no specific definition of IR in children, rather, that for adult population is used by pediatric endocrinologists instead. Altered insulin secretion dynamics are associated with worse metabolic profiles and type 2 diabetes mellitus development, thus we aimed to test whether insulin response relates to oxidative stress and inflammation in children. METHODS: We conducted a case-control study, including 132 children classified as follows: 33 children without obesity (Lean); 42 with obesity but no IR according to the American Diabetes Association criteria for adults (OBIR-); 25 with obesity and IR and an early insulin response to an oral glucose tolerance test (OGTT) (EP-OBIR +); 32 with obesity, IR, and a late insulin peak (LP-OBIR +); and studied variables associated with lipid and carbohydrate metabolism, oxidative stress, inflammation and inflammasome activation. RESULTS: The measured parameters of children with obesity, IR, and an early insulin response were similar to those of children with obesity but without IR. It was late responders who presented an impaired antioxidant system and elevated oxidative damage in erythrocytes and plasma, and inflammasome activation at their white blood cells, despite lower classical inflammation markers. Increased uric acid levels seems to be one of the underlying mechanisms for inflammasome activation. CONCLUSIONS: It is insulin response to an OGTT that identifies children with obesity suffering oxidative stress and inflammasome activation more specifically. Uric acid could be mediating this pathological inflammatory response by activating NLRP3 in peripheral blood mononuclear cells.

A blood microRNA classifier for the prediction of ICU mortality in COVID-19 patients: a multicenter validation study.

BACKGROUND: The identification of critically ill COVID-19 patients at risk of fatal outcomes remains a challenge. Here, we first validated candidate microRNAs (miRNAs) as biomarkers for clinical decision-making in critically ill patients. Second, we constructed a blood miRNA classifier for the early prediction of adverse outcomes in the ICU. METHODS: This was a multicenter, observational and retrospective/prospective study including 503 critically ill patients admitted to the ICU from 19 hospitals. qPCR assays were performed in plasma samples collected within the first 48 h upon admission. A 16-miRNA panel was designed based on recently published data from our group. RESULTS: Nine miRNAs were validated as biomarkers of all-cause in-ICU mortality in the independent cohort of critically ill patients (FDR < 0.05). Cox regression analysis revealed that low expression levels of eight miRNAs were associated with a higher risk of death (HR from 1.56 to 2.61). LASSO regression for variable selection was used to construct a miRNA classifier. A 4-blood miRNA signature composed of miR-16-5p, miR-192-5p, miR-323a-3p and miR-451a predicts the risk of all-cause in-ICU mortality (HR 2.5). Kaplan‒Meier analysis confirmed these findings. The miRNA signature provides a significant increase in the prognostic capacity of conventional scores, APACHE-II (C-index 0.71, DeLong test p-value 0.055) and SOFA (C-index 0.67, DeLong test p-value 0.001), and a risk model based on clinical predictors (C-index 0.74, DeLong test-p-value 0.035). For 28-day and 90-day mortality, the classifier also improved the prognostic value of APACHE-II, SOFA and the clinical model. The association between the classifier and mortality persisted even after multivariable adjustment. The functional analysis reported biological pathways involved in SARS-CoV infection and inflammatory, fibrotic and transcriptional pathways. CONCLUSIONS: A blood miRNA classifier improves the early prediction of fatal outcomes in critically ill COVID-19 patients.

Childhood obesity, metabolic syndrome, and oxidative stress: microRNAs go on stage.

The incidence of childhood obesity and metabolic syndrome has grown notably in the last years, becoming major public health burdens in developed countries. Nowadays, oxidative stress is well-recognized to be closely associated with the onset and progression of several obesity-related complications within the framework of a complex crosstalk involving other intertwined pathogenic events, such as inflammation, insulin disturbances, and dyslipidemia. Thus, understanding the molecular basis behind these oxidative dysregulations could provide new approaches for the diagnosis, prevention, and treatment of childhood obesity and associated disorders. In this respect, the transcriptomic characterization of miRNAs bares great potential because of their involvement in post-transcriptional modulation of genetic expression. Herein, we provide a comprehensive literature revision gathering state-of-the-art research into the association between childhood obesity, metabolic syndrome, and miRNAs. We put special emphasis on the potential role of miRNAs in modulating obesity-related pathogenic events, with particular focus on oxidative stress.

Chromatin regulation by Histone H4 acetylation at Lysine 16 during cell death and differentiation in the myeloid compartment.

Histone H4 acetylation at Lysine 16 (H4K16ac) is a key epigenetic mark involved in gene regulation, DNA repair and chromatin remodeling, and though it is known to be essential for embryonic development, its role during adult life is still poorly understood. Here we show that this lysine is massively hyperacetylated in peripheral neutrophils. Genome-wide mapping of H4K16ac in terminally differentiated blood cells, along with functional experiments, supported a role for this histone post-translational modification in the regulation of cell differentiation and apoptosis in the hematopoietic system. Furthermore, in neutrophils, H4K16ac was enriched at specific DNA repeats. These DNA regions presented an accessible chromatin conformation and were associated with the cleavage sites that generate the 50 kb DNA fragments during the first stages of programmed cell death. Our results thus suggest that H4K16ac plays a dual role in myeloid cells as it not only regulates differentiation and apoptosis, but it also exhibits a non-canonical structural role in poising chromatin for cleavage at an early stage of neutrophil cell death.

Loss of 5hmC identifies a new type of aberrant DNA hypermethylation in glioma.

Aberrant DNA hypermethylation is a hallmark of cancer although the underlying molecular mechanisms are still poorly understood. To study the possible role of 5-hydroxymethylcytosine (5hmC) in this process we analyzed the global and locus-specific genome-wide levels of 5hmC and 5-methylcytosine (5mC) in human primary samples from 12 non-tumoral brains and 53 gliomas. We found that the levels of 5hmC identified in non-tumoral samples were significantly reduced in gliomas. Strikingly, hypo-hydroxymethylation at 4627 (9.3%) CpG sites was associated with aberrant DNA hypermethylation and was strongly enriched in CpG island shores. The DNA regions containing these CpG sites were enriched in H3K4me2 and presented a different genuine chromatin signature to that characteristic of the genes classically aberrantly hypermethylated in cancer. As this 5mC gain is inversely correlated with loss of 5hmC and has not been identified with classical sodium bisulfite-based technologies, we conclude that our data identifies a novel 5hmC-dependent type of aberrant DNA hypermethylation in glioma.

Age-associated hydroxymethylation in human bone-marrow mesenchymal stem cells.

BACKGROUND: Age-associated changes in genomic DNA methylation have been primarily attributed to 5-methylcytosine (5mC). However, the recent discovery of 5-hydroxymethylcytosine (5hmC) suggests that this epigenetic mark might also play a role in the process. METHODS: Here, we analyzed the genome-wide profile of 5hmc in mesenchymal stem cells (MSCs) obtained from bone-marrow donors, aged 2-89 years. RESULTS: We identified 10,685 frequently hydroxymethylated CpG sites in MSCs that were, as in other cell types, significantly associated with low density CpG regions, introns, the histone posttranslational modification H3k4me1 and enhancers. Study of the age-associated changes to 5hmC identified 785 hyper- and 846 hypo-hydroxymethylated CpG sites in the MSCs obtained from older individuals. CONCLUSIONS: DNA hyper-hydroxymethylation in the advanced-age group was associated with loss of 5mC, which suggests that, at specific CpG sites, this epigenetic modification might play a role in DNA methylation changes during lifetime. Since bone-marrow MSCs have many clinical applications, and the fact that the epigenomic alterations in this cell type associated with aging identified in this study could have associated functional effects, the age of donors should be taken into account in clinical settings.

Exploring the translational landscape of the long noncoding RNA transcriptome in acute respiratory distress syndrome: it is a long way to the top.

Acute respiratory distress syndrome (ARDS) poses a significant and widespread public health challenge. Extensive research conducted in recent decades has considerably improved our understanding of the disease pathophysiology. Nevertheless, ARDS continues to rank among the leading causes of mortality in intensive care units and its management remains a formidable task, primarily due to its remarkable heterogeneity. As a consequence, the syndrome is underdiagnosed, prognostication has important gaps and selection of the appropriate therapeutic approach is laborious. In recent years, the noncoding transcriptome has emerged as a new area of attention for researchers interested in biomarker development. Numerous studies have confirmed the potential of long noncoding RNAs (lncRNAs), transcripts with little or no coding information, as noninvasive tools for diagnosis, prognosis and prediction of the therapeutic response across a broad spectrum of ailments, including respiratory conditions. This article aims to provide a comprehensive overview of lncRNAs with specific emphasis on their role as biomarkers. We review current knowledge on the circulating lncRNAs as potential markers that can be used to enhance decision making in ARDS management. Additionally, we address the primary limitations and outline the steps that will be essential for integration of the use of lncRNAs in clinical laboratories. Our ultimate objective is to provide a framework for the implementation of lncRNAs in the management of ARDS.

microRNA Expression Profile in Obesity-Induced Kidney Disease Driven by High-Fat Diet in Mice.

Obesity is one of the main causes of chronic kidney disease; however, the precise molecular mechanisms leading to the onset of kidney injury and dysfunction in obesity-associated nephropathy remain unclear. The present study aimed to unveil the kidney microRNA (miRNA) expression profile in a model of obesity-induced kidney disease in C57BL/6J mice using next-generation sequencing (NGS) analysis. High-fat diet (HFD)-induced obesity led to notable structural alterations in tubular and glomerular regions of the kidney, increased renal expression of proinflammatory and profibrotic genes, as well as an elevated renal expression of genes involved in cellular lipid metabolism. The miRNA sequencing analysis identified a set of nine miRNAs differentially expressed in the kidney upon HFD feeding, with miR-5099, miR-551b-3p, miR-223-3p, miR-146a-3p and miR-21a-3p showing the most significant differential expression between standard diet (STD) and HFD mice. A validation analysis showed that the expression levels of miR-5099, miR-551b-3p and miR-146a-3p were consistent with NGS results, while Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment analyses revealed that these three validated miRNAs modulated target genes involved in metabolic and adipocytokine pathways, fatty acid and lipid metabolism, and inflammatory, senescence and profibrotic pathways. Our results suggest that differentially expressed miRNAs play pivotal roles in the intricate pathophysiology of obesity-associated kidney disease and could potentially create novel treatment strategies to counteract the deleterious effects of obesity on kidney function.

Addressing the unsolved challenges in microRNA-based biomarker development: Suitable endogenous reference microRNAs for SARS-CoV-2 infection severity.

Circulating cell-free microRNAs (miRNAs) are promising biomarkers for medical decision-making. Suitable endogenous controls are essential to ensure reproducibility. We aimed to identify and validate endogenous reference miRNAs for qPCR data normalization in samples from SARS-CoV-2-infected hospitalized patients. We used plasma samples (n = 170) from COVID-19 patients collected at hospital admission (COVID-Ponent project, www.clinicaltrials.gov/NCT04824677). First, 179 miRNAs were profiled using RT-qPCR. After stability assessment, candidates were validated using the same methodology. miRNA stability was analyzed using the geNorm, NormFinder and BestKeeper algorithms. Stability was further evaluated using an RNA-seq dataset derived from COVID-19 hospitalized patients, along with plasma samples from patients with critical COVID-19 profiled using RT-qPCR. In the screening phase, after strict control of expression levels, stability assessment selected eleven candidates (miR-17-5p, miR-20a-5p, miR-30e-5p, miR-106a-5p, miR-151a-5p, miR-185-5p, miR-191-5p, miR-423-3p, miR-425-5p, miR-484 and miR-625-5p). In the validation phase, all algorithms identified miR-106a-5p and miR-484 as top endogenous controls. No association was observed between these miRNAs and clinical or sociodemographic characteristics. Both miRNAs were stably detected and showed low variability in the additional analyses. In conclusion, a 2-miRNA panel composed of miR-106a-5p and miR-484 constitutes a first-line normalizer for miRNA-based biomarker development using qPCR in hospitalized patients infected with SARS-CoV-2.

MicroRNA-guided drug discovery for mitigating persistent pulmonary complications in critical COVID-19 survivors: A longitudinal pilot study.

BACKGROUND AND PURPOSE: The post-acute sequelae of SARS-CoV-2 infection pose a significant global challenge, with nearly 50% of critical COVID-19 survivors manifesting persistent lung abnormalities. The lack of understanding about the molecular mechanisms and effective treatments hampers their management. Here, we employed microRNA (miRNA) profiling to decipher the systemic molecular underpinnings of the persistent pulmonary complications. EXPERIMENTAL APPROACH: We conducted a longitudinal investigation including 119 critical COVID-19 survivors. A comprehensive pulmonary evaluation was performed in the short-term (median = 94.0 days after hospital discharge) and long-term (median = 358 days after hospital discharge). Plasma miRNAs were quantified at the short-term evaluation using the gold-standard technique, RT-qPCR. The analyses combined machine learning feature selection techniques with bioinformatic investigations. Two additional datasets were incorporated for validation. KEY RESULTS: In the short-term, 84% of the survivors exhibited impaired lung diffusion (DLCO  < 80% of predicted). One year post-discharge, 54.4% of this patient subgroup still presented abnormal DLCO . Four feature selection methods identified two specific miRNAs, miR-9-5p and miR-486-5p, linked to persistent lung dysfunction. The downstream experimentally validated targetome included 1473 genes, with heterogeneous enriched pathways associated with inflammation, angiogenesis and cell senescence. Validation studies using RNA-sequencing and proteomic datasets emphasized the pivotal roles of cell migration and tissue repair in persistent lung dysfunction. The repositioning potential of the miRNA targets was limited. CONCLUSION AND IMPLICATIONS: Our study reveals early mechanistic pathways contributing to persistent lung dysfunction in critical COVID-19 survivors, offering a promising approach for the development of targeted disease-modifying agents.

Association between Acute Respiratory Distress Syndrome Due to COVID-19 and Long-Term Sleep and Circadian Sleep-Wake Disorders.

Current studies agree on the impact of sleep and circadian rest-activity rhythm alterations in acute respiratory distress syndrome (ARDS) survivors. However, research on the duration of this impact is scarce. In this study, we evaluate the impact of ARDS on the sleep and circadian rest-activity rhythm of COVID-19 survivors twelve months after hospital discharge. This is a prospective study including COVID-19 survivors with and without ARDS during hospitalization. Data was collected four and twelve months after hospital discharge. The interventions included one-week wrist actigraphy and a home sleep apnea test (HSAT), and evaluations were conducted according to the Pittsburgh sleep quality index (PSQI), Epworth sleepiness scale (ESS), and insomnia severity index (ISI). Fifty-two patients were evaluated (ARDS = 31 and non-ARDS = 21); they had a median age of 49.0 [39.0;57.2] years and 53.8% were male. After twelve months, 91.3% presented poor sleep quality, 58.7% presented insomnia, 50% presented daytime somnolence, and 37% presented comorbid insomnia and obstructive sleep apnea (COMISA). No significant improvement was observed in relation to sleep or the circadian rest-activity rhythm between four and twelve months. A tendency of poor sleep quality, insomnia, daytime somnolence, and COMISA was observed. Finally, there was no significant impact on the circadian rest-activity rhythm between four and twelve months or between the groups.

Circulating lipoprotein-carried miRNome analysis reveals novel VLDL-enriched microRNAs that strongly correlate with the HDL-microRNA profile.

Lipoproteins have been described as microRNAs (miRNAs) carriers. Unfortunately, the bibliography on this topic is scarce and shows a high variability between independent investigations. In addition, the miRNA profiles of the LDL and VLDL fractions have not been completely elucidated. Here, we profiled the human circulating lipoprotein-carried miRNome. Lipoprotein fractions (VLDL, LDL and HDL) were isolated from the serum of healthy subjects by ultracentrifugation and purified by size-exclusion chromatography. A panel of 179 miRNAs commonly expressed in circulation was evaluated in the lipoprotein fractions using quantitative real-time PCR (qPCR) assays. A total of 14, 4 and 24 miRNAs were stably detected in the VLDL, LDL and HDL fractions, respectively. VLDL- and HDL-miRNA signatures were highly correlated (rho 0.814), and miR-16-5p, miR-142-3p, miR-223-3p and miR-451a were among the top 5 expressed miRNAs in both fractions. miR-125a-5p, miR-335-3p and miR-1260a, were detected in all lipoprotein fractions. miR-107 and miR-221-3p were uniquely detected in the VLDL fraction. HDL showed the larger number of specifically detected miRNAs (n = 13). Enrichment in specific miRNA families and genomic clusters was observed for HDL-miRNAs. Two sequence motifs were also detected for this group of miRNAs. Functional enrichment analysis including the miRNA signatures from each lipoprotein fraction suggested a potential role in mechanistic pathways previously associated with cardiovascular disease: fibrosis, senescence, inflammation, immune response, angiogenesis, and cardiomyopathy. Collectively, our results not only support the role of lipoproteins as circulating miRNA carriers but also describe for the first time the role of VLDL as a miRNA transporter.

Circulating miRNAs signature on breast cancer: the MCC-Spain project.

PURPOSE: To build models combining circulating microRNAs (miRNAs) able to identify women with breast cancer as well as different types of breast cancer, when comparing with controls without breast cancer. METHOD: miRNAs analysis was performed in two phases: screening phase, with a total n = 40 (10 controls and 30 BC cases) analyzed by Next Generation Sequencing, and validation phase, which included 131 controls and 269 cases. For this second phase, the miRNAs were selected combining the screening phase results and a revision of the literature. They were quantified using RT-PCR. Models were built using logistic regression with LASSO penalization. RESULTS: The model for all cases included seven miRNAs (miR-423-3p, miR-139-5p, miR-324-5p, miR-1299, miR-101-3p, miR-186-5p and miR-29a-3p); which had an area under the ROC curve of 0.73. The model for cases diagnosed via screening only took in one miRNA (miR-101-3p); the area under the ROC curve was 0.63. The model for disease-free cases in the follow-up had five miRNAs (miR-101-3p, miR-186-5p, miR-423-3p, miR-142-3p and miR-1299) and the area under the ROC curve was 0.73. Finally, the model for cases with active disease in the follow-up contained six miRNAs (miR-101-3p, miR-423-3p, miR-139-5p, miR-1307-3p, miR-331-3p and miR-21-3p) and its area under the ROC curve was 0.82. CONCLUSION: We present four models involving eleven miRNAs to differentiate healthy controls from different types of BC cases. Our models scarcely overlap with those previously reported.

Study protocol for the epigenetic characterization of angor pectoris according to the affected coronary compartment: Global and comprehensive assessment of the relationship between invasive coronary physiology and microRNAs.

BACKGROUND: MicroRNAs (miRNAs) are noncoding RNAs involved in post-transcriptional genetic regulation with a proposed role in intercellular communication. miRNAs are considered promising biomarkers in ischemic heart disease. Invasive physiological evaluation allows a precise assessment of each affected coronary compartment. Although some studies have associated the expression of circulating miRNAs with invasive physiological indexes, their global relationship with coronary compartments has not been assessed. Here, we will evaluate circulating miRNAs profiles according to the coronary pattern of the vascular compartment affectation. STUDY AND DESIGN: This is an investigator-initiated, multicentre, descriptive study to be conducted at three centres in Spain (NCT05374694). The study will include one hundred consecutive patients older than 18 years with chest pain of presumed coronary cause undergoing invasive physiological evaluation, including fractional flow reserve (FFR) and index of microvascular resistance (IMR). Patients will be initially classified into four groups, according to FFR and IMR: macrovascular and microvascular affectation (FFR≤0.80 / IMR≥25), isolated macrovascular affectation (FFR≤0.80 / IMR<25), isolated microvascular affectation (FFR>0.80 / IMR ≥25) and normal coronary indexes (FFR>0.80 / IMR<25). Patients with isolated microvascular affectation or normal indexes will also undergo the acetylcholine test and may be reclassified as a fifth group in the presence of spasm. A panel of miRNAs previously associated with molecular mechanisms linked to chronic coronary syndrome will be analysed using RT-qPCR. CONCLUSIONS: The results of this study will identify miRNA profiles associated with patterns of coronary affectation and will contribute to a better understanding of the mechanistic pathways of coronary pathology.

Sleep and breast and prostate cancer risk in the MCC-Spain study.

Breast and prostate cancers have been associated with circadian disruption. Some previous studies examined associations of sleep duration and breast or prostate cancer risk though findings remain inconsistent. This study examines associations of a range of detailed sleep characteristics and breast and prostate cancer risk in a large-scale population-based case-control study, MCC-Spain. A total of 1738 incident breast cancer cases, 1112 prostate cancer cases and frequency matched controls (n = 1910, and 1493 respectively) were recruited. Detailed data on habitual sleep duration, quality, timing, and daytime napping ("siesta") were collected at recruitment. Additional data on sleep habits during both the previous year and at age 40 years were also subsequently captured. Adjusted odds ratios (ORs) and 95% confidence intervals (CI) were estimated. There were no associations of habitual sleep duration (h), timing of sleep, or any or specific sleep problems, and either breast and prostate cancer risk. There was a significant positive association of ever taking habitual siestas at recruitment and breast cancer risk (OR = 1.22, 95% CI 1.06-1.42), which strengthened with increased frequency or duration. There were also significant positive associations observed for both breast and prostate cancer, among those reporting recent sleep problems, but not sleep problems at age 40 years, in a subsequent circadian questionnaire. Adverse associations with siesta and disturbed sleep during the previous year likely reflect symptoms of developing/diagnosed cancer and comorbidities. Overall, there was no clear association between various sleep characteristics and breast or prostate cancer risk observed.

Multiplex protein profiling of bronchial aspirates reveals disease-, mortality- and respiratory sequelae-associated signatures in critically ill patients with ARDS secondary to SARS-CoV-2 infection.

Introduction: Bronchial aspirates (BAS) obtained during invasive mechanical ventilation (IMV) constitutes a useful tool for molecular phenotyping and decision making. Aim: To identify the proteomic determinants associated with disease pathogenesis, all-cause mortality and respiratory sequelae in BAS samples from critically ill patients with SARS-CoV-2-induced ARDS. Methods: Multicenter study including 74 critically ill patients with COVID-19 and non-COVID-19 ARDS. BAS were obtained by bronchoaspiration after IMV initiation. Three hundred sixty-four proteins were quantified using proximity extension assay (PEA) technology. Random forest models were used to assess predictor importance. Results: After adjusting for confounding factors, CST5, NADK, SRPK2 and TGF-α were differentially detected in COVID-19 and non-COVID-19 patients. In random forest models for COVID-19, CST5, DPP7, NADK, KYAT1 and TYMP showed the highest variable importance. In COVID-19 patients, reduced levels of ENTPD2 and PTN were observed in nonsurvivors of ICU stay, even after adjustment. AGR2, NQO2, IL-1α, OSM and TRAIL showed the strongest associations with in-ICU mortality and were used to construct a protein-based prediction model. Kaplan-Meier curves revealed a clear separation in mortality risk between subgroups of PTN, ENTPD2 and the prediction model. Cox regression models supported these findings. In survivors, the levels of FCRL1, NTF4 and THOP1 in BAS samples obtained during the ICU stay correlated with lung function (i.e., DLCO levels) 3 months after hospital discharge. Similarly, Flt3L and THOP1 levels were correlated with radiological features (i.e., TSS). These proteins are expressed in immune and nonimmune lung cells. Poor host response to viral infectivity and an inappropriate reparative mechanism seem to be linked with the pathogenesis of the disease and fatal outcomes, respectively. Conclusion: BAS proteomics identified novel factors associated with the pathology of SARS-CoV-2-induced ARDS and its adverse outcomes. BAS-based protein testing emerges as a novel tool for risk assessment in the ICU.

Genome-wide transcriptional profiling of pulmonary functional sequelae in ARDS- secondary to SARS-CoV-2 infection.

BACKGROUND: Up to 80% of patients surviving acute respiratory distress syndrome (ARDS) secondary to SARS-CoV-2 infection present persistent anomalies in pulmonary function after hospital discharge. There is a limited understanding of the mechanistic pathways linked to post-acute pulmonary sequelae. AIM: To identify the molecular underpinnings associated with severe lung diffusion involvement in survivors of SARS-CoV-2-induced ARDS. METHODS: Survivors attended to a complete pulmonary evaluation 3 months after hospital discharge. RNA sequencing (RNA-seq) was performed using Illumina technology in whole-blood samples from 50 patients with moderate to severe diffusion impairment (DLCO<60%) and age- and sex-matched individuals with mild-normal lung function (DLCO≥60%). A transcriptomic signature for optimal classification was constructed using random forest. Transcriptomic data were analyzed for biological pathway enrichment, cellular deconvolution, cell/tissue-specific gene expression and candidate drugs. RESULTS: RNA-seq identified 1357 differentially expressed transcripts. A model composed of 14 mRNAs allowed the optimal discrimination of survivors with severe diffusion impairment (AUC=0.979). Hallmarks of lung sequelae involved cell death signaling, cytoskeleton reorganization, cell growth and differentiation and the immune response. Resting natural killer (NK) cells were the most important immune cell subtype for the prediction of severe diffusion impairment. Components of the signature correlated with neutrophil, lymphocyte and monocyte counts. A variable expression profile of the transcripts was observed in lung cell subtypes and bodily tissues. One upregulated gene, TUBB4A, constitutes a target for FDA-approved drugs. CONCLUSIONS: This work defines the transcriptional programme associated with post-acute pulmonary sequelae and provides novel insights for targeted interventions and biomarker development.

Toenail zinc as a biomarker: Relationship with sources of environmental exposure and with genetic variability in MCC-Spain study.

BACKGROUND: Toenails are commonly used as biomarkers of exposure to zinc (Zn), but there is scarce information about their relationship with sources of exposure to Zn. OBJECTIVES: To investigate the main determinants of toenail Zn, including selected sources of environmental exposure to Zn and individual genetic variability in Zn metabolism. METHODS: We determined toenail Zn by inductively coupled plasma mass spectrometry in 3,448 general population controls from the MultiCase-Control study MCC-Spain. We assessed dietary and supplement Zn intake using food frequency questionnaires, residential proximity to Zn-emitting industries and residential topsoil Zn levels through interpolation methods. We constructed a polygenic score of genetic variability based on 81 single nucleotide polymorphisms in genes involved in Zn metabolism. Geometric mean ratios of toenail Zn across categories of each determinant were estimated from multivariate linear regression models on log-transformed toenail Zn. RESULTS: Geometric mean toenail Zn was 104.1 µg/g in men and 100.3 µg/g in women. Geometric mean toenail Zn levels were 7 % lower (95 % confidence interval 1-13 %) in men older than 69 years and those in the upper tertile of fibre intake, and 9 % higher (3-16 %) in smoking men. Women residing within 3 km from Zn-emitting industries had 4 % higher geometric mean toenail Zn levels (0-9 %). Dietary Zn intake and polygenic score were unrelated to toenail Zn. Overall, the available determinants only explained 9.3 % of toenail Zn variability in men and 4.8 % in women. DISCUSSION: Sociodemographic factors, lifestyle, diet, and environmental exposure explained little of the individual variability of toenail Zn in the study population. The available genetic variants related to Zn metabolism were not associated with toenail Zn.