exRNA Atlas Analysis Reveals Distinct Extracellular RNA Cargo Types and Their Carriers Present across Human Biofluids.

To develop a map of cell-cell communication mediated by extracellular RNA (exRNA), the NIH Extracellular RNA Communication Consortium created the exRNA Atlas resource (https://exrna-atlas.org). The Atlas version 4P1 hosts 5,309 exRNA-seq and exRNA qPCR profiles from 19 studies and a suite of analysis and visualization tools. To analyze variation between profiles, we apply computational deconvolution. The analysis leads to a model with six exRNA cargo types (CT1, CT2, CT3A, CT3B, CT3C, CT4), each detectable in multiple biofluids (serum, plasma, CSF, saliva, urine). Five of the cargo types associate with known vesicular and non-vesicular (lipoprotein and ribonucleoprotein) exRNA carriers. To validate utility of this model, we re-analyze an exercise response study by deconvolution to identify physiologically relevant response pathways that were not detected previously. To enable wide application of this model, as part of the exRNA Atlas resource, we provide tools for deconvolution and analysis of user-provided case-control studies.

Tracking Distinct RNA Populations Using Efficient and Reversible Covalent Chemistry.

We describe a chemical method to label and purify 4-thiouridine (s(4)U)-containing RNA. We demonstrate that methanethiosulfonate (MTS) reagents form disulfide bonds with s(4)U more efficiently than the commonly used HPDP-biotin, leading to higher yields and less biased enrichment. This increase in efficiency allowed us to use s(4)U labeling to study global microRNA (miRNA) turnover in proliferating cultured human cells without perturbing global miRNA levels or the miRNA processing machinery. This improved chemistry will enhance methods that depend on tracking different populations of RNA, such as 4-thiouridine tagging to study tissue-specific transcription and dynamic transcriptome analysis (DTA) to study RNA turnover.

The Dynamic Platelet Transcriptome in Obesity and Weight Loss.

OBJECTIVE: Adiposity is associated with oxidative stress, inflammation, and glucose intolerance. Previous data suggest that platelet gene expression is associated with key cardiometabolic phenotypes, including body mass index but stable in healthy individuals over time. However, modulation of gene expression in platelets in response to metabolic shifts (eg, weight reduction) is unknown and may be important to defining mechanism. Approach and Results: Platelet RNA sequencing and aggregation were performed from 21 individuals with massive weight loss (>45 kg) following bariatric surgery. Based on RNA sequencing data, we measured the expression of 67 genes from isolated platelet RNA using high-throughput quantitative reverse transcription quantitative PCR in 1864 FHS (Framingham Heart Study) participants. Many transcripts not previously studied in platelets were differentially expressed with bariatric surgical weight loss, appeared specific to platelets (eg, not differentially expressed in leukocytes), and were enriched for a nonalcoholic fatty liver disease pathway. Platelet aggregation studies did not detect alteration in platelet function after significant weight loss. Linear regression models demonstrated several platelet genes modestly associated with cross-sectional cardiometabolic phenotypes, including body mass index. There were no associations between studied transcripts and incident diabetes or cardiovascular end points. CONCLUSIONS: In summary, while there is no change in platelet aggregation function after significant weight loss, the human platelet experiences a dramatic transcriptional shift that implicates pathways potentially relevant to improved cardiometabolic risk postweight loss (eg, nonalcoholic fatty liver disease). Further studies are needed to determine the mechanistic importance of these observations.

Evaluation of serological lateral flow assays for severe acute respiratory syndrome coronavirus-2.

BACKGROUND: COVID-19 has resulted in significant morbidity and mortality worldwide. Lateral flow assays can detect anti-Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) antibodies to monitor transmission. However, standardized evaluation of their accuracy and tools to aid in interpreting results are needed. METHODS: We evaluated 20 IgG and IgM assays selected from available tests in April 2020. We evaluated the assays' performance using 56 pre-pandemic negative and 56 SARS-CoV-2-positive plasma samples, collected 10-40 days after symptom onset, confirmed by a molecular test and analyzed by an ultra-sensitive immunoassay. Finally, we developed a user-friendly web app to extrapolate the positive predictive values based on their accuracy and local prevalence. RESULTS: Combined IgG + IgM sensitivities ranged from 33.9 to 94.6%, while combined specificities ranged from 92.6 to 100%. The highest sensitivities were detected in Lumiquick for IgG (98.2%), BioHit for both IgM (96.4%), and combined IgG + IgM sensitivity (94.6%). Furthermore, 11 LFAs and 8 LFAs showed perfect specificity for IgG and IgM, respectively, with 15 LFAs showing perfect combined IgG + IgM specificity. Lumiquick had the lowest estimated limit-of-detection (LOD) (0.1 µg/mL), followed by a similar LOD of 1.5 µg/mL for CareHealth, Cellex, KHB, and Vivachek. CONCLUSION: We provide a public resource of the accuracy of select lateral flow assays with potential for home testing. The cost-effectiveness, scalable manufacturing process, and suitability for self-testing makes LFAs an attractive option for monitoring disease prevalence and assessing vaccine responsiveness. Our web tool provides an easy-to-use interface to demonstrate the impact of prevalence and test accuracy on the positive predictive values.

Comparative analysis of the transcriptome across distant species.

The transcriptome is the readout of the genome. Identifying common features in it across distant species can reveal fundamental principles. To this end, the ENCODE and modENCODE consortia have generated large amounts of matched RNA-sequencing data for human, worm and fly. Uniform processing and comprehensive annotation of these data allow comparison across metazoan phyla, extending beyond earlier within-phylum transcriptome comparisons and revealing ancient, conserved features. Specifically, we discover co-expression modules shared across animals, many of which are enriched in developmental genes. Moreover, we use expression patterns to align the stages in worm and fly development and find a novel pairing between worm embryo and fly pupae, in addition to the embryo-to-embryo and larvae-to-larvae pairings. Furthermore, we find that the extent of non-canonical, non-coding transcription is similar in each organism, per base pair. Finally, we find in all three organisms that the gene-expression levels, both coding and non-coding, can be quantitatively predicted from chromatin features at the promoter using a 'universal model' based on a single set of organism-independent parameters.

Isoform-Level Interpretation of High-Throughput Proteomics Data Enabled by Deep Integration with RNA-seq.

Cellular control of gene expression is a complex process that is subject to multiple levels of regulation, but ultimately it is the protein produced that determines the biosynthetic state of the cell. One way that a cell can regulate the protein output from each gene is by expressing alternate isoforms with distinct amino acid sequences. These isoforms may exhibit differences in localization and binding interactions that can have profound functional implications. High-throughput liquid chromatography tandem mass spectrometry proteomics (LC-MS/MS) relies on enzymatic digestion and has lower coverage and sensitivity than transcriptomic profiling methods such as RNA-seq. Digestion results in predictable fragmentation of a protein, which can limit the generation of peptides capable of distinguishing between isoforms. Here we exploit transcript-level expression from RNA-seq to set prior likelihoods and enable protein isoform abundances to be directly estimated from LC-MS/MS, an approach derived from the principle that most genes appear to be expressed as a single dominant isoform in a given cell type or tissue. Through this deep integration of RNA-seq and LC-MS/MS data from the same sample, we show that a principal isoform can be identified in >80% of gene products in homogeneous HEK293 cell culture and >70% of proteins detected in complex human brain tissue. We demonstrate that the incorporation of translatome data from ribosome profiling further refines this process. Defining isoforms in experiments with matched RNA-seq/translatome and proteomic data increases the functional relevance of such data sets and will further broaden our understanding of multilevel control of gene expression.

Development of a highly automated and multiplexed targeted proteome pipeline and assay for 112 rat brain synaptic proteins.

We present a comprehensive workflow for large scale (>1000 transitions/run) label-free LC-MRM proteome assays. Innovations include automated MRM transition selection, intelligent retention time scheduling that improves S/N by twofold, and automatic peak modeling. Improvements to data analysis include a novel Q/C metric, normalized group area ratio, MLR normalization, weighted regression analysis, and data dissemination through the Yale protein expression database. As a proof of principle we developed a robust 90 min LC-MRM assay for mouse/rat postsynaptic density fractions which resulted in the routine quantification of 337 peptides from 112 proteins based on 15 observations per protein. Parallel analyses with stable isotope dilution peptide standards (SIS), demonstrate very high correlation in retention time (1.0) and protein fold change (0.94) between the label-free and SIS analyses. Overall, our method achieved a technical CV of 11.4% with >97.5% of the 1697 transitions being quantified without user intervention, resulting in a highly efficient, robust, and single injection LC-MRM assay.

Data-driven identification of predictive risk biomarkers for subgroups of osteoarthritis using interpretable machine learning.

Osteoarthritis (OA) is increasing in prevalence and has a severe impact on patients' lives. However, our understanding of biomarkers driving OA risk remains limited. We developed a model predicting the five-year risk of OA diagnosis, integrating retrospective clinical, lifestyle and biomarker data from the UK Biobank (19,120 patients with OA, ROC-AUC: 0.72, 95%CI (0.71-0.73)). Higher age, BMI and prescription of non-steroidal anti-inflammatory drugs contributed most to increased OA risk prediction ahead of diagnosis. We identified 14 subgroups of OA risk profiles. These subgroups were validated in an independent set of patients evaluating the 11-year OA risk, with 88% of patients being uniquely assigned to one of the 14 subgroups. Individual OA risk profiles were characterised by personalised biomarkers. Omics integration demonstrated the predictive importance of key OA genes and pathways (e.g., GDF5 and TGF-ß signalling) and OA-specific biomarkers (e.g., CRTAC1 and COL9A1). In summary, this work identifies opportunities for personalised OA prevention and insights into its underlying pathogenesis.

Longitudinal analysis of serum-derived extracellular vesicle RNA to monitor dacomitinib treatment response in EGFR-amplified recurrent glioblastoma patients.

Background: Glioblastoma (GBM) is a highly aggressive and invasive brain tumor associated with high patient mortality. A large fraction of GBM tumors have been identified as epidermal growth factor receptor (EGFR) amplified and ~50% also are EGFRvIII mutant positive. In a previously reported multicenter phase II study, we have described the response of recurrent GBM (rGBM) patients to dacomitinib, an EGFR tyrosine kinase inhibitor (TKI). As a continuation of that report, we leverage the tumor cargo-encapsulating extracellular vesicles (EVs) and explore their genetic composition as carriers of tumor biomarker. Methods: Serum samples were longitudinally collected from EGFR-amplified rGBM patients who clinically benefitted from dacomitinib therapy (responders) and those who did not (nonresponders), as well as from a healthy cohort of individuals. The serum EV transcriptome was evaluated to map the RNA biotype distribution and distinguish GBM disease. Results: Using long RNA sequencing, we show enriched detection of over 10 000 coding RNAs from serum EVs. The EV transcriptome yielded a unique signature that facilitates differentiation of GBM patients from healthy donors. Further analysis revealed genetic enrichment that enables stratification of responders from nonresponders prior to dacomitinib treatment as well as following administration. Conclusion: This study demonstrates that genetic composition analysis of serum EVs may aid in therapeutic stratification to identify patients with dacomitinib-responsive GBM.

Hepatocyte mARC1 promotes fatty liver disease.

Background & Aims: Non-alcoholic fatty liver disease (NAFLD) has a prevalence of ∼25% worldwide, with significant public health consequences yet few effective treatments. Human genetics can help elucidate novel biology and identify targets for new therapeutics. Genetic variants in mitochondrial amidoxime-reducing component 1 (MTARC1) have been associated with NAFLD and liver-related mortality; however, its pathophysiological role and the cell type(s) mediating these effects remain unclear. We aimed to investigate how MTARC1 exerts its effects on NAFLD by integrating human genetics with in vitro and in vivo studies of mARC1 knockdown. Methods: Analyses including multi-trait colocalisation and Mendelian randomisation were used to assess the genetic associations of MTARC1. In addition, we established an in vitro long-term primary human hepatocyte model with metabolic readouts and used the Gubra Amylin NASH (GAN)-diet non-alcoholic steatohepatitis mouse model treated with hepatocyte-specific N-acetylgalactosamine (GalNAc)-siRNA to understand the in vivo impacts of MTARC1. Results: We showed that genetic variants within the MTARC1 locus are associated with liver enzymes, liver fat, plasma lipids, and body composition, and these associations are attributable to the same causal variant (p.A165T, rs2642438 G>A), suggesting a shared mechanism. We demonstrated that increased MTARC1 mRNA had an adverse effect on these traits using Mendelian randomisation, implying therapeutic inhibition of mARC1 could be beneficial. In vitro mARC1 knockdown decreased lipid accumulation and increased triglyceride secretion, and in vivo GalNAc-siRNA-mediated knockdown of mARC1 lowered hepatic but increased plasma triglycerides. We found alterations in pathways regulating lipid metabolism and decreased secretion of 3-hydroxybutyrate upon mARC1 knockdown in vitro and in vivo. Conclusions: Collectively, our findings from human genetics, and in vitro and in vivo hepatocyte-specific mARC1 knockdown support the potential efficacy of hepatocyte-specific targeting of mARC1 for treatment of NAFLD. Impact and implications: We report that genetically predicted increases in MTARC1 mRNA associate with poor liver health. Furthermore, knockdown of mARC1 reduces hepatic steatosis in primary human hepatocytes and a murine NASH model. Together, these findings further underscore the therapeutic potential of targeting hepatocyte MTARC1 for NAFLD.

Maternal separation with early weaning: a rodent model providing novel insights into neglect associated developmental deficits.

Child neglect is the most prevalent form of child maltreatment in the United States, and poses a serious public health concern. Children who survive such episodes go on to experience long-lasting psychological and behavioral problems, including higher rates of post-traumatic stress disorder symptoms, depression, alcohol and drug abuse, attention-deficit/hyperactivity disorder, and cognitive deficits. To date, most research into the causes of these life-long problems has focused on well-established targets such as stress responsive systems, including the hypothalamus-pituitary-adrenal axis. Using the maternal separation and early weaning model, we have attempted to provide comprehensive molecular profiling of a model of early-life neglect in an organism amenable to genomic manipulation: the mouse. In this article, we report new findings generated with this model using chromatin immunoprecipitation sequencing, diffuse tensor magnetic resonance imaging, and behavioral analyses. We also review the validity of the maternal separation and early weaning model, which reflects behavioral deficits observed in neglected humans including hyperactivity, anxiety, and attentional deficits. Finally, we summarize the molecular characterization of these animals, including RNA profiling and label-free proteomics, which highlight protein translation and myelination as novel pathways of interest.

Differential global gene expression in cystic fibrosis nasal and bronchial epithelium.

Respiratory epithelium is the target of therapies, such as gene therapy, for cystic fibrosis (CF) lung disease. To determine the usefulness of the nasal epithelium as a pre-screen for lung-directed therapies, we profiled gene expression in CF and non-CF nasal and bronchial epithelium samples using Illumina HumanRef-8 Expression BeadChips. 863 genes were differentially expressed between CF and non-CF bronchial epithelium but only 15 were differentially expressed between CF and non-CF nasal epithelium (≥1.5-fold, P≤0.05). The most enriched pathway in CF bronchial epithelium was inflammatory response, whereas in CF nasal epithelium it was amino acid metabolism. We also compared nasal and bronchial epithelium in each group and identified differential expression of cellular movement genes in CF patients and cellular growth genes in non-CF subjects. We conclude that CF and non-CF nasal and bronchial epithelium are transcriptionally distinct and CF nasal epithelium is not a good surrogate for the lung respiratory epithelium.

Proteomic characterization of post-mortem human brain tissue following ultracentrifugation-based subcellular fractionation.

Proteomic characterization of human brain tissue is increasingly utilized to identify potential novel biomarkers and drug targets for a variety of neurological diseases. In whole-tissue studies, results may be driven by changes in the proportion of the largest and most abundant organelles or tissue cell-type composition. Spatial proteomics approaches enhance our knowledge of disease mechanisms and changing signalling pathways at the subcellular level by taking into account the importance of cellular localization, which critically influences protein function. Density gradient-based ultracentrifugation methods allow for subcellular fractionation and have been utilized in cell lines, mouse and human brain tissue to quantify thousands of proteins in specific enriched organelles such as the pre- and post-synapse. Serial ultracentrifugation methods allow for the analysis of multiple cellular organelles from the same biological sample, and to our knowledge have not been previously applied to frozen post-mortem human brain tissue. The use of frozen human tissue for tissue fractionation faces two major challenges, the post-mortem interval, during which proteins may leach from their usual location into the cytosol, and freezing, which results in membrane breakdown. Despite these challenges, in this proof-of-concept study, we show that the majority of proteins segregate reproducibly into crude density-based centrifugation fractions, that the fractions are enriched for the appropriate organellar markers and that significant differences in protein localization can be observed between tissue from individuals with Alzheimer's disease and control individuals.

Technical Performance Evaluation of Olink Proximity Extension Assay for Blood-Based Biomarker Discovery in Longitudinal Studies of Alzheimer's Disease.

The core Alzheimer's disease (AD) cerebrospinal fluid (CSF) biomarkers; amyloid-ß (Aß), total tau (t-tau), and phosphorylated tau (p-tau181), are strong indicators of the presence of AD pathology, but do not correlate well with disease progression, and can be difficult to implement in longitudinal studies where repeat biofluid sampling is required. As a result, blood-based biomarkers are increasingly being sought as alternatives. In this study, we aimed to evaluate a promising blood biomarker discovery technology, Olink Proximity Extension Assays for technical reproducibility characteristics in order to highlight the advantages and disadvantages of using this technology in biomarker discovery in AD. We evaluated the performance of five Olink Proteomic multiplex proximity extension assays (PEA) in plasma samples. Three technical control samples included on each plate allowed calculation of technical variability. Biotemporal stability was measured in three sequential annual samples from 54 individuals with and without AD. Coefficients of variation (CVs), analysis of variance (ANOVA), and variance component analyses were used to quantify technical and individual variation over time. We show that overall, Olink assays are technically robust, with the largest experimental variation stemming from biological differences between individuals for most analytes. As a powerful illustration of one of the potential pitfalls of using a multi-plexed technology for discovery, we performed power calculations using the baseline samples to demonstrate the size of study required to overcome the need for multiple test correction with this technology. We show that the power of moderate effect size proteins was strongly reduced, and as a result investigators should strongly consider pooling resources to perform larger studies using this multiplexed technique where possible.

Exosome based analysis for Space Associated Neuro-Ocular Syndrome and health risks in space exploration.

Molecular profiling to characterize the effects of environmental exposures is important from the human health and performance as well as the occupational medicine perspective in space exploration. We have developed a novel exosome-based platform that allows profiling of biological processes in the body from a variety of body fluids. The technology is suitable for diagnostic applications as well as studying the pathophysiology of the Space Associated Neuro-Ocular Syndrome in astronauts and monitoring patients with chronically impaired cerebrospinal fluid drainage or elevated intracranial pressure. In this proof-of-concept, we demonstrate that: (a) exosomes from different biofluids contain a specific population of RNA transcripts; (b) urine collection hardware aboard the ISS is compatible with exosome gene expression technology; (c) cDNA libraries from exosomal RNA can be stored in dry form and at room temperature, representing an interesting option for the creation of longitudinal molecular catalogs that can be stored as a repository for retrospective analysis.

Plasma Proteomics of COVID-19-Associated Cardiovascular Complications: Implications for Pathophysiology and Therapeutics.

To gain insights into the mechanisms driving cardiovascular complications in COVID-19, we performed a case-control plasma proteomics study in COVID-19 patients. Our results identify the senescence-associated secretory phenotype, a marker of biological aging, as the dominant process associated with disease severity and cardiac involvement. FSTL3, an indicator of senescence-promoting Activin/TGFß signaling, and ADAMTS13, the von Willebrand Factor-cleaving protease whose loss-of-function causes microvascular thrombosis, were among the proteins most strongly associated with myocardial stress and injury. Findings were validated in a larger COVID-19 patient cohort and the hamster COVID-19 model, providing new insights into the pathophysiology of COVID-19 cardiovascular complications with therapeutic implications.

Phase 2 of extracellular RNA communication consortium charts next-generation approaches for extracellular RNA research.

The extracellular RNA communication consortium (ERCC) is an NIH-funded program aiming to promote the development of new technologies, resources, and knowledge about exRNAs and their carriers. After Phase 1 (2013-2018), Phase 2 of the program (ERCC2, 2019-2023) aims to fill critical gaps in knowledge and technology to enable rigorous and reproducible methods for separation and characterization of both bulk populations of exRNA carriers and single EVs. ERCC2 investigators are also developing new bioinformatic pipelines to promote data integration through the exRNA atlas database. ERCC2 has established several Working Groups (Resource Sharing, Reagent Development, Data Analysis and Coordination, Technology Development, nomenclature, and Scientific Outreach) to promote collaboration between ERCC2 members and the broader scientific community. We expect that ERCC2's current and future achievements will significantly improve our understanding of exRNA biology and the development of accurate and efficient exRNA-based diagnostic, prognostic, and theranostic biomarker assays.

Relative impact of key sources of systematic noise in Affymetrix and Illumina gene-expression microarray experiments.

BACKGROUND: Systematic processing noise, which includes batch effects, is very common in microarray experiments but is often ignored despite its potential to confound or compromise experimental results. Compromised results are most likely when re-analysing or integrating datasets from public repositories due to the different conditions under which each dataset is generated. To better understand the relative noise-contributions of various factors in experimental-design, we assessed several Illumina and Affymetrix datasets for technical variation between replicate hybridisations of Universal Human Reference (UHRR) and individual or pooled breast-tumour RNA. RESULTS: A varying degree of systematic noise was observed in each of the datasets, however in all cases the relative amount of variation between standard control RNA replicates was found to be greatest at earlier points in the sample-preparation workflow. For example, 40.6% of the total variation in reported expressions were attributed to replicate extractions, compared to 13.9% due to amplification/labelling and 10.8% between replicate hybridisations. Deliberate probe-wise batch-correction methods were effective in reducing the magnitude of this variation, although the level of improvement was dependent on the sources of noise included in the model. Systematic noise introduced at the chip, run, and experiment levels of a combined Illumina dataset were found to be highly dependent upon the experimental design. Both UHRR and pools of RNA, which were derived from the samples of interest, modelled technical variation well although the pools were significantly better correlated (4% average improvement) and better emulated the effects of systematic noise, over all probes, than the UHRRs. The effect of this noise was not uniform over all probes, with low GC-content probes found to be more vulnerable to batch variation than probes with a higher GC-content. CONCLUSIONS: The magnitude of systematic processing noise in a microarray experiment is variable across probes and experiments, however it is generally the case that procedures earlier in the sample-preparation workflow are liable to introduce the most noise. Careful experimental design is important to protect against noise, detailed meta-data should always be provided, and diagnostic procedures should be routinely performed prior to downstream analyses for the detection of bias in microarray studies.

Statistical aspects of quantitative real-time PCR experiment design.

Experiments using quantitative real-time PCR to test hypotheses are limited by technical and biological variability; we seek to minimise sources of confounding variability through optimum use of biological and technical replicates. The quality of an experiment design is commonly assessed by calculating its prospective power. Such calculations rely on knowledge of the expected variances of the measurements of each group of samples and the magnitude of the treatment effect; the estimation of which is often uninformed and unreliable. Here we introduce a method that exploits a small pilot study to estimate the biological and technical variances in order to improve the design of a subsequent large experiment. We measure the variance contributions at several 'levels' of the experiment design and provide a means of using this information to predict both the total variance and the prospective power of the assay. A validation of the method is provided through a variance analysis of representative genes in several bovine tissue-types. We also discuss the effect of normalisation to a reference gene in terms of the measured variance components of the gene of interest. Finally, we describe a software implementation of these methods, powerNest, that gives the user the opportunity to input data from a pilot study and interactively modify the design of the assay. The software automatically calculates expected variances, statistical power, and optimal design of the larger experiment. powerNest enables the researcher to minimise the total confounding variance and maximise prospective power for a specified maximum cost for the large study.