Persistent complement dysregulation with signs of thromboinflammation in active Long Covid

Long Covid is a debilitating condition of unknown etiology. We performed multimodal proteomics analyses of blood serum from COVID-19 patients followed up to 12 months after confirmed severe acute respiratory syndrome coronavirus 2 infection. Analysis of >6500 proteins in 268 longitudinal samples revealed dysregulated activation of the complement system, an innate immune protection and homeostasis mechanism, in individuals experiencing Long Covid. Thus, active Long Covid was characterized by terminal complement system dysregulation and ongoing activation of the alternative and classical complement pathways, the latter associated with increased antibody titers against several herpesviruses possibly stimulating this pathway. Moreover, markers of hemolysis, tissue injury, platelet activation, and monocyte-platelet aggregates were increased in Long Covid. Machine learning confirmed complement and thromboinflammatory proteins as top biomarkers, warranting diagnostic and therapeutic interrogation of these systems.

Dual-specificity protein phosphatase 6 (DUSP6) overexpression reduces amyloid load and improves memory deficits in male 5xFAD mice

Background

Dual specificity protein phosphatase 6 (DUSP6) was recently identified as a key hub gene in a causal network that regulates late-onset Alzheimer's disease. Importantly, decreased DUSP6 levels are correlated with an increased clinical dementia rating in human subjects, and DUSP6 levels are additionally decreased in the 5xFAD amyloidopathy mouse model.

Methods

AAV5-DUSP6 or AAV5-GFP (control) were stereotactically injected into the dorsal hippocampus (dHc) of female and male 5xFAD or wild type mice to overexpress DUSP6 or GFP. Spatial learning memory of these mice was assessed in the Barnes maze, after which hippocampal tissues were isolated for downstream analysis.

Results

Barnes maze testing indicated that DUSP6 overexpression in the dHc of 5xFAD mice improved memory deficits and was associated with reduced amyloid plaque load, Aß 1-40 and Aß 1-42 levels, and amyloid precursor protein processing enzyme BACE1, in male but not in female mice. Microglial activation and microgliosis, which are increased in 5xFAD mice, were significantly reduced by dHc DUSP6 overexpression in both males and females. Transcriptomic profiling of female 5xFAD hippocampus revealed upregulated expression of genes involved in inflammatory and extracellular signal-regulated kinase (ERK) pathways, while dHc DUSP6 overexpression in female 5xFAD mice downregulated a subset of genes in these pathways. A limited number of differentially expressed genes (DEGs) (FDR<0.05) were identified in male mice; gene ontology analysis of DEGs (p<0.05) identified a greater number of synaptic pathways that were regulated by DUSP6 overexpression in male compared to female 5xFAD. Notably, the msh homeobox 3 gene, Msx3 , previously shown to regulate microglial M1/M2 polarization and reduce neuroinflammation, was one of the most robustly upregulated genes in female and male wild type and 5xFAD mice overexpressing DUSP6.

Conclusions

In summary, our data indicate that DUSP6 overexpression in dHc reduced amyloid deposition and memory deficits in male but not female 5xFAD mice, whereas reduced neuroinflammation and microglial activation were observed in both males and females. The sex-dependent regulation of synaptic pathways by DUSP6 overexpression, however, correlated with the improvement of spatial memory deficits in male but not female 5xFAD.

A study of gene expression in the living human brain

A goal of medical research is to determine the molecular basis of human brain health and illness. One way to achieve this goal is through observational studies of gene expression in human brain tissue. Due to the unavailability of brain tissue from living people, most such studies are performed using tissue from postmortem brain donors. An assumption underlying this practice is that gene expression in the postmortem human brain is an accurate representation of gene expression in the living human brain. Here, this assumption - which, until now, had not been adequately tested - is tested by comparing human prefrontal cortex gene expression between 275 living samples and 243 postmortem samples. Expression levels differed significantly for nearly 80% of genes, and a systematic examination of alternative explanations for this observation determined that these differences are not a consequence of cell type composition, RNA quality, postmortem interval, age, medication, morbidity, symptom severity, tissue pathology, sample handling, batch effects, or computational methods utilized. Analyses integrating the data generated for this study with data from earlier landmark studies that used tissue from postmortem brain donors showed that postmortem brain gene expression signatures of neurological and mental illnesses, as well as of normal traits such as aging, may not be accurate representations of these gene expression signatures in the living brain. By using tissue from large cohorts living people, future observational studies of human brain biology have the potential to (1) determine the medical research questions that can be addressed using postmortem tissue as a proxy for living tissue and (2) expand the scope of medical research to include questions about the molecular basis of human brain health and illness that can only be addressed in living people (e.g., "What happens at the molecular level in the brain as a person experiences an emotion?").

Circulating proteins to predict COVID-19 severity

Predicting COVID-19 severity is difficult, and the biological pathways involved are not fully understood. To approach this problem, we measured 4701 circulating human protein abundances in two independent cohorts totaling 986 individuals. We then trained prediction models including protein abundances and clinical risk factors to predict COVID-19 severity in 417 subjects and tested these models in a separate cohort of 569 individuals. For severe COVID-19, a baseline model including age and sex provided an area under the receiver operator curve (AUC) of 65% in the test cohort. Selecting 92 proteins from the 4701 unique protein abundances improved the AUC to 88% in the training cohort, which remained relatively stable in the testing cohort at 86%, suggesting good generalizability. Proteins selected from different COVID-19 severity were enriched for cytokine and cytokine receptors, but more than half of the enriched pathways were not immune-related. Taken together, these findings suggest that circulating proteins measured at early stages of disease progression are reasonably accurate predictors of COVID-19 severity. Further research is needed to understand how to incorporate protein measurement into clinical care.

Molecular states during acute COVID-19 reveal distinct etiologies of long-term sequelae

Post-acute sequelae of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection are debilitating, clinically heterogeneous and of unknown molecular etiology. A transcriptome-wide investigation was performed in 165 acutely infected hospitalized individuals who were followed clinically into the post-acute period. Distinct gene expression signatures of post-acute sequelae were already present in whole blood during acute infection, with innate and adaptive immune cells implicated in different symptoms. Two clusters of sequelae exhibited divergent plasma-cell-associated gene expression patterns. In one cluster, sequelae associated with higher expression of immunoglobulin-related genes in an anti-spike antibody titer-dependent manner. In the other, sequelae associated independently of these titers with lower expression of immunoglobulin-related genes, indicating lower non-specific antibody production in individuals with these sequelae. This relationship between lower total immunoglobulins and sequelae was validated in an external cohort. Altogether, multiple etiologies of post-acute sequelae were already detectable during SARS-CoV-2 infection, directly linking these sequelae with the acute host response to the virus and providing early insights into their development.

Dual-Specificity Protein Phosphatase 4 (DUSP4) Overexpression Improves Learning Behavior Selectively in Female 5xFAD Mice, and Reduces β-Amyloid Load in Males and Females

Recent multiscale network analyses of banked brains from subjects who died of late-onset sporadic Alzheimer's disease converged on VGF (non-acronymic) as a key hub or driver. Within this computational VGF network, we identified the dual-specificity protein phosphatase 4 (DUSP4) [also known as mitogen-activated protein kinase (MAPK) phosphatase 2] as an important node. Importantly, DUSP4 gene expression, like that of VGF, is downregulated in postmortem Alzheimer's disease (AD) brains. We investigated the roles that this VGF/DUSP4 network plays in the development of learning behavior impairment and neuropathology in the 5xFAD amyloidopathy mouse model. We found reductions in DUSP4 expression in the hippocampi of male AD subjects, correlating with increased CDR scores, and in 4-month-old female and 12-18-month-old male 5xFAD hippocampi. Adeno-associated virus (AAV5)-mediated overexpression of DUSP4 in 5xFAD mouse dorsal hippocampi (dHc) rescued impaired Barnes maze performance in females but not in males, while amyloid loads were reduced in both females and males. Bulk RNA sequencing of the dHc from 5-month-old mice overexpressing DUSP4, and Ingenuity Pathway and Enrichr analyses of differentially expressed genes (DEGs), revealed that DUSP4 reduced gene expression in female 5xFAD mice in neuroinflammatory, interferon-gamma (IFNγ), programmed cell death protein-ligand 1/programmed cell death protein 1 (PD-L1/PD-1), and extracellular signal-regulated kinase (ERK)/MAPK pathways, via which DUSP4 may modulate AD phenotype with gender-specificity.

Exome-wide association study to identify rare variants influencing COVID-19 outcomes: Results from the Host Genetics Initiative

Host genetics is a key determinant of COVID-19 outcomes. Previously, the COVID-19 Host Genetics Initiative genome-wide association study used common variants to identify multiple loci associated with COVID-19 outcomes. However, variants with the largest impact on COVID-19 outcomes are expected to be rare in the population. Hence, studying rare variants may provide additional insights into disease susceptibility and pathogenesis, thereby informing therapeutics development. Here, we combined whole-exome and whole-genome sequencing from 21 cohorts across 12 countries and performed rare variant exome-wide burden analyses for COVID-19 outcomes. In an analysis of 5,085 severe disease cases and 571,737 controls, we observed that carrying a rare deleterious variant in the SARS-CoV-2 sensor toll-like receptor TLR7 (on chromosome X) was associated with a 5.3-fold increase in severe disease (95% CI: 2.75-10.05, p = 5.41x10-7). This association was consistent across sexes. These results further support TLR7 as a genetic determinant of severe disease and suggest that larger studies on rare variants influencing COVID-19 outcomes could provide additional insights.

A shift in lung macrophage composition is associated with COVID-19 severity and recovery

Although it has been more than 2 years since the start of the coronavirus disease 2019 (COVID-19) pandemic, COVID-19 continues to be a worldwide health crisis. Despite the development of preventive vaccines, therapies to treat COVID-19 and other inflammatory diseases remain a major unmet need in medicine. Our study sought to identify drivers of disease severity and mortality to develop tailored immunotherapy strategies to halt disease progression. We assembled the Mount Sinai COVID-19 Biobank, which was composed of almost 600 hospitalized patients followed longitudinally through the peak of the pandemic in 2020. Moderate disease and survival were associated with a stronger antigen presentation and effector T cell signature. In contrast, severe disease and death were associated with an altered antigen presentation signature, increased numbers of inflammatory immature myeloid cells, and extrafollicular activated B cells that have been previously associated with autoantibody formation. In severely ill patients with COVID-19, lung tissue-resident alveolar macrophages not only were drastically depleted but also had an altered antigen presentation signature, which coincided with an influx of inflammatory monocytes and monocyte-derived macrophages. In addition, we found that the size of the alveolar macrophage pool correlated with patient outcome and that alveolar macrophage numbers and functionality were restored to homeostasis in patients who recovered from COVID-19. These data suggest that local and systemic myeloid cell dysregulation are drivers of COVID-19 severity and modulation of alveolar macrophage numbers and activity in the lung may be a viable therapeutic strategy for the treatment of critical inflammatory lung diseases.

A translational genomics approach identifies IL10RB as the top candidate gene target for COVID-19 susceptibility

Recent efforts have identified genetic loci that are associated with coronavirus disease 2019 (COVID-19) infection rates and disease outcome severity. Translating these genetic findings into druggable genes that reduce COVID-19 host susceptibility is a critical next step. Using a translational genomics approach that integrates COVID-19 genetic susceptibility variants, multi-tissue genetically regulated gene expression (GReX), and perturbagen signatures, we identified IL10RB as the top candidate gene target for COVID-19 host susceptibility. In a series of validation steps, we show that predicted GReX upregulation of IL10RB and higher IL10RB expression in COVID-19 patient blood is associated with worse COVID-19 outcomes and that in vitro IL10RB overexpression is associated with increased viral load and activation of disease-relevant molecular pathways.

Shift of lung macrophage composition is associated with COVID-19 disease severity and recovery

Though it has been 2 years since the start of the Coronavirus Disease 19 (COVID-19) pandemic, COVID-19 continues to be a worldwide health crisis. Despite the development of preventive vaccines, very little progress has been made to identify curative therapies to treat COVID-19 and other inflammatory diseases which remain a major unmet need in medicine. Our study sought to identify drivers of disease severity and death to develop tailored immunotherapy strategies to halt disease progression. Here we assembled the Mount Sinai COVID-19 Biobank which was comprised of ~600 hospitalized patients followed longitudinally during the peak of the pandemic. Moderate disease and survival were associated with a stronger antigen (Ag) presentation and effector T cell signature, while severe disease and death were associated with an altered Ag presentation signature, increased numbers of circulating inflammatory, immature myeloid cells, and extrafollicular activated B cells associated with autoantibody formation. Strikingly, we found that in severe COVID-19 patients, lung tissue resident alveolar macrophages (AM) were not only severely depleted, but also had an altered Ag presentation signature, and were replaced by inflammatory monocytes and monocyte-derived macrophages (MoMΦ). Notably, the size of the AM pool correlated with recovery or death, while AM loss and functionality were restored in patients that recovered. These data therefore suggest that local and systemic myeloid cell dysregulation is a driver of COVID-19 severity and that modulation of AM numbers and functionality in the lung may be a viable therapeutic strategy for the treatment of critical lung inflammatory illnesses.

Acute COVID-19 gene-expression profiles show multiple etiologies of long-term sequelae

Two years into the SARS-CoV-2 pandemic, the post-acute sequelae of infection are compounding the global health crisis. Often debilitating, these sequelae are clinically heterogeneous and of unknown molecular etiology. Here, a transcriptome-wide investigation of this new condition was performed in a large cohort of acutely infected patients followed clinically into the post-acute period. Gene expression signatures of post-acute sequelae were already present in whole blood during the acute phase of infection, with both innate and adaptive immune cells involved. Plasma cells stood out as driving at least two distinct clusters of sequelae, one largely dependent on circulating antibodies against the SARS-CoV-2 spike protein and the other antibody-independent. Altogether, multiple etiologies of post-acute sequelae were found concomitant with SARS-CoV-2 infection, directly linking the emergence of these sequelae with the host response to the virus.

Prognostic value of polygenic risk scores for adults with psychosis

Polygenic risk scores (PRS) summarize genetic liability to a disease at the individual level, and the aim is to use them as biomarkers of disease and poor outcomes in real-world clinical practice. To date, few studies have assessed the prognostic value of PRS relative to standards of care. Schizophrenia (SCZ), the archetypal psychotic illness, is an ideal test case for this because the predictive power of the SCZ PRS exceeds that of most other common diseases. Here, we analyzed clinical and genetic data from two multi-ethnic cohorts totaling 8,541 adults with SCZ and related psychotic disorders, to assess whether the SCZ PRS improves the prediction of poor outcomes relative to clinical features captured in a standard psychiatric interview. For all outcomes investigated, the SCZ PRS did not improve the performance of predictive models, an observation that was generally robust to divergent case ascertainment strategies and the ancestral background of the study participants.

Downregulation of exhausted cytotoxic T cells in gene expression networks of multisystem inflammatory syndrome in children

Multisystem inflammatory syndrome in children (MIS-C) presents with fever, inflammation and pathology of multiple organs in individuals under 21 years of age in the weeks following severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. Although an autoimmune pathogenesis has been proposed, the genes, pathways and cell types causal to this new disease remain unknown. Here we perform RNA sequencing of blood from patients with MIS-C and controls to find disease-associated genes clustered in a co-expression module annotated to CD56dimCD57+ natural killer (NK) cells and exhausted CD8+ T cells. A similar transcriptome signature is replicated in an independent cohort of Kawasaki disease (KD), the related condition after which MIS-C was initially named. Probing a probabilistic causal network previously constructed from over 1,000 blood transcriptomes both validates the structure of this module and reveals nine key regulators, including TBX21, a central coordinator of exhausted CD8+ T cell differentiation. Together, this unbiased, transcriptome-wide survey implicates downregulation of NK cells and cytotoxic T cell exhaustion in the pathogenesis of MIS-C.

IL10RB as a key regulator of COVID-19 host susceptibility and severity

BACKGROUND

Recent efforts have identified genetic loci that are associated with coronavirus disease 2019 (COVID-19) infection rates and disease outcome severity. Translating these genetic findings into druggable genes and readily available compounds that reduce COVID-19 host susceptibility is a critical next step.

METHODS

We integrate COVID-19 genetic susceptibility variants, multi-tissue genetically regulated gene expression (GReX) and perturbargen signatures to identify candidate genes and compounds that reverse the predicted gene expression dysregulation associated with COVID-19 susceptibility. The top candidate gene is validated by testing both its GReX and observed blood transcriptome association with COVID-19 severity, as well as by in vitro perturbation to quantify effects on viral load and molecular pathway dysregulation. We validate the in silico drug repositioning analysis by examining whether the top candidate compounds decrease COVID-19 incidence based on epidemiological evidence.

RESULTS

We identify IL10RB as the top key regulator of COVID-19 host susceptibility. Predicted GReX up-regulation of IL10RB and higher IL10RB expression in COVID-19 patient blood is associated with worse COVID-19 outcomes. In vitro IL10RB overexpression is associated with increased viral load and activation of immune-related molecular pathways. Azathioprine and retinol are prioritized as candidate compounds to reduce the likelihood of testing positive for COVID-19.

CONCLUSIONS

We establish an integrative data-driven approach for gene target prioritization. We identify and validate IL10RB as a suitable molecular target for modulation of COVID-19 host susceptibility. Finally, we provide evidence for a few readily available medications that would warrant further investigation as drug repositioning candidates.

Toward a fine-scale population health monitoring system

Understanding population health disparities is an essential component of equitable precision health efforts. Epidemiology research often relies on definitions of race and ethnicity, but these population labels may not adequately capture disease burdens and environmental factors impacting specific sub-populations. Here, we propose a framework for repurposing data from electronic health records (EHRs) in concert with genomic data to explore the demographic ties that can impact disease burdens. Using data from a diverse biobank in New York City, we identified 17 communities sharing recent genetic ancestry. We observed 1,177 health outcomes that were statistically associated with a specific group and demonstrated significant differences in the segregation of genetic variants contributing to Mendelian diseases. We also demonstrated that fine-scale population structure can impact the prediction of complex disease risk within groups. This work reinforces the utility of linking genomic data to EHRs and provides a framework toward fine-scale monitoring of population health.

A Molecular network approach reveals shared cellular and molecular signatures between chronic fatigue syndrome and other fatiguing illnesses

The molecular mechanisms of chronic fatigue syndrome (CFS, or Myalgic encephalomyelitis), a disease defined by extreme, long-term fatigue, remain largely uncharacterized, and presently no molecular diagnostic test and no specific treatments exist to diagnose and treat CFS patients. While CFS has historically had an estimated prevalence of 0.1-0.5% [1], concerns of a “long hauler” version of Coronavirus disease 2019 (COVID-19) that symptomatically overlaps CFS to a significant degree(Supplemental Table-1)and appears to occur in 10% of COVID-19 patients[2], has raised concerns of a larger spike in CFS [3]. Here, we established molecular signatures of CFS and a corresponding network-based disease context from RNA-sequencing data generated on whole blood and FACs sorted specific peripheral blood mononuclear cells (PBMCs) isolated from CFS cases and non-CFS controls. The immune cell type specific molecular signatures of CFS we identified, overlapped molecular signatures from other fatiguing illnesses, demonstrating a common molecular etiology. Further, after constructing a probabilistic causal model of the CFS gene expression data, we identified master regulator genes modulating network states associated with CFS, suggesting potential therapeutic targets for CFS.

An integrative multiomic network model links lipid metabolism to glucose regulation in coronary artery disease

Elevated plasma cholesterol and type 2 diabetes (T2D) are associated with coronary artery disease (CAD). Individuals treated with cholesterol-lowering statins have increased T2D risk, while individuals with hypercholesterolemia have reduced T2D risk. We explore the relationship between lipid and glucose control by constructing network models from the STARNET study with sequencing data from seven cardiometabolic tissues obtained from CAD patients during coronary artery by-pass grafting surgery. By integrating gene expression, genotype, metabolomic, and clinical data, we identify a glucose and lipid determining (GLD) regulatory network showing inverse relationships with lipid and glucose traits. Master regulators of the GLD network also impact lipid and glucose levels in inverse directions. Experimental inhibition of one of the GLD network master regulators, lanosterol synthase (LSS), in mice confirms the inverse relationships to glucose and lipid levels as predicted by our model and provides mechanistic insights.

Transformative Network Modeling of Multi-omics Data Reveals Detailed Circuits, Key Regulators, and Potential Therapeutics for Alzheimer's Disease

To identify the molecular mechanisms and novel therapeutic targets of late-onset Alzheimer's Disease (LOAD), we performed an integrative network analysis of multi-omics profiling of four cortical areas across 364 donors with varying cognitive and neuropathological phenotypes. Our analyses revealed thousands of molecular changes and uncovered neuronal gene subnetworks as the most dysregulated in LOAD. ATP6V1A was identified as a key regulator of a top-ranked neuronal subnetwork, and its role in disease-related processes was evaluated through CRISPR-based manipulation in human induced pluripotent stem cell-derived neurons and RNAi-based knockdown in Drosophila models. Neuronal impairment and neurodegeneration caused by ATP6V1A deficit were improved by a repositioned compound, NCH-51. This study provides not only a global landscape but also detailed signaling circuits of complex molecular interactions in key brain regions affected by LOAD, and the resulting network models will serve as a blueprint for developing next-generation therapeutic agents against LOAD.

Predictive network modeling in human induced pluripotent stem cells identifies key driver genes for insulin responsiveness

Insulin resistance (IR) precedes the development of type 2 diabetes (T2D) and increases cardiovascular disease risk. Although genome wide association studies (GWAS) have uncovered new loci associated with T2D, their contribution to explain the mechanisms leading to decreased insulin sensitivity has been very limited. Thus, new approaches are necessary to explore the genetic architecture of insulin resistance. To that end, we generated an iPSC library across the spectrum of insulin sensitivity in humans. RNA-seq based analysis of 310 induced pluripotent stem cell (iPSC) clones derived from 100 individuals allowed us to identify differentially expressed genes between insulin resistant and sensitive iPSC lines. Analysis of the co-expression architecture uncovered several insulin sensitivity-relevant gene sub-networks, and predictive network modeling identified a set of key driver genes that regulate these co-expression modules. Functional validation in human adipocytes and skeletal muscle cells (SKMCs) confirmed the relevance of the key driver candidate genes for insulin responsiveness.

Insulin resistance is characterized by a defective response (“resistance”) to normal insulin concentrations to uptake the glucose present in the blood, and is the underlying condition that leads to type 2 diabetes (T2D) and increases the risk of cardiovascular disease. It is estimated that 25–33% of the US population are insulin resistant enough to be at risk of serious clinical consequences. For more than a decade, large population studies have tried to discover the genes that participate in the development of insulin resistance, but without much success. It is now increasingly clear that the complex genetic nature of insulin resistance requires novel approaches centered in patient specific cellular models. To fill this gap, we have generated an induced pluripotent stem cell (iPSC) library from individuals with accurate measurements of insulin sensitivity, and performed gene expression and key driver analyses. Our work demonstrates that iPSCs can be used as a revolutionary technology to model insulin resistance and to discover key genetic drivers. Moreover, they can develop our basic knowledge of the disease, and are ultimately expected to increase the therapeutic targets to treat insulin resistance and type 2 diabetes.

Machine Learning to Predict Mortality and Critical Events in a Cohort of Patients With COVID-19 in New York City: Model Development and Validation

BACKGROUND

COVID-19 has infected millions of people worldwide and is responsible for several hundred thousand fatalities. The COVID-19 pandemic has necessitated thoughtful resource allocation and early identification of high-risk patients. However, effective methods to meet these needs are lacking.

OBJECTIVE

The aims of this study were to analyze the electronic health records (EHRs) of patients who tested positive for COVID-19 and were admitted to hospitals in the Mount Sinai Health System in New York City; to develop machine learning models for making predictions about the hospital course of the patients over clinically meaningful time horizons based on patient characteristics at admission; and to assess the performance of these models at multiple hospitals and time points.

METHODS

We used Extreme Gradient Boosting (XGBoost) and baseline comparator models to predict in-hospital mortality and critical events at time windows of 3, 5, 7, and 10 days from admission. Our study population included harmonized EHR data from five hospitals in New York City for 4098 COVID-19-positive patients admitted from March 15 to May 22, 2020. The models were first trained on patients from a single hospital (n=1514) before or on May 1, externally validated on patients from four other hospitals (n=2201) before or on May 1, and prospectively validated on all patients after May 1 (n=383). Finally, we established model interpretability to identify and rank variables that drive model predictions.

RESULTS

Upon cross-validation, the XGBoost classifier outperformed baseline models, with an area under the receiver operating characteristic curve (AUC-ROC) for mortality of 0.89 at 3 days, 0.85 at 5 and 7 days, and 0.84 at 10 days. XGBoost also performed well for critical event prediction, with an AUC-ROC of 0.80 at 3 days, 0.79 at 5 days, 0.80 at 7 days, and 0.81 at 10 days. In external validation, XGBoost achieved an AUC-ROC of 0.88 at 3 days, 0.86 at 5 days, 0.86 at 7 days, and 0.84 at 10 days for mortality prediction. Similarly, the unimputed XGBoost model achieved an AUC-ROC of 0.78 at 3 days, 0.79 at 5 days, 0.80 at 7 days, and 0.81 at 10 days. Trends in performance on prospective validation sets were similar. At 7 days, acute kidney injury on admission, elevated LDH, tachypnea, and hyperglycemia were the strongest drivers of critical event prediction, while higher age, anion gap, and C-reactive protein were the strongest drivers of mortality prediction.

CONCLUSIONS

We externally and prospectively trained and validated machine learning models for mortality and critical events for patients with COVID-19 at different time horizons. These models identified at-risk patients and uncovered underlying relationships that predicted outcomes.

An inflammatory cytokine signature predicts COVID-19 severity and survival

Several studies have revealed that the hyper-inflammatory response induced by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a major cause of disease severity and death. However, predictive biomarkers of pathogenic inflammation to help guide targetable immune pathways are critically lacking. We implemented a rapid multiplex cytokine assay to measure serum interleukin (IL)-6, IL-8, tumor necrosis factor (TNF)-α and IL-1β in hospitalized patients with coronavirus disease 2019 (COVID-19) upon admission to the Mount Sinai Health System in New York. Patients (n = 1,484) were followed up to 41 d after admission (median, 8 d), and clinical information, laboratory test results and patient outcomes were collected. We found that high serum IL-6, IL-8 and TNF-α levels at the time of hospitalization were strong and independent predictors of patient survival (P < 0.0001, P = 0.0205 and P = 0.0140, respectively). Notably, when adjusting for disease severity, common laboratory inflammation markers, hypoxia and other vitals, demographics, and a range of comorbidities, IL-6 and TNF-α serum levels remained independent and significant predictors of disease severity and death. These findings were validated in a second cohort of patients (n = 231). We propose that serum IL-6 and TNF-α levels should be considered in the management and treatment of patients with COVID-19 to stratify prospective clinical trials, guide resource allocation and inform therapeutic options.

An inflammatory cytokine signature predicts COVID-19 severity and survival

Several studies have revealed that the hyper-inflammatory response induced by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a major cause of disease severity and death. However, predictive biomarkers of pathogenic inflammation to help guide targetable immune pathways are critically lacking. We implemented a rapid multiplex cytokine assay to measure serum interleukin (IL)-6, IL-8, tumor necrosis factor (TNF)-α and IL-1β in hospitalized patients with coronavirus disease 2019 (COVID-19) upon admission to the Mount Sinai Health System in New York. Patients (n = 1,484) were followed up to 41 d after admission (median, 8 d), and clinical information, laboratory test results and patient outcomes were collected. We found that high serum IL-6, IL-8 and TNF-α levels at the time of hospitalization were strong and independent predictors of patient survival (P < 0.0001, P = 0.0205 and P = 0.0140, respectively). Notably, when adjusting for disease severity, common laboratory inflammation markers, hypoxia and other vitals, demographics, and a range of comorbidities, IL-6 and TNF-α serum levels remained independent and significant predictors of disease severity and death. These findings were validated in a second cohort of patients (n = 231). We propose that serum IL-6 and TNF-α levels should be considered in the management and treatment of patients with COVID-19 to stratify prospective clinical trials, guide resource allocation and inform therapeutic options.

Cytotoxic lymphocytes are dysregulated in multisystem inflammatory syndrome in children

Multisystem inflammatory syndrome in children (MIS-C) presents with fever, inflammation and multiple organ involvement in individuals under 21 years following severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. To identify genes, pathways and cell types driving MIS-C, we sequenced the blood transcriptomes of MIS-C cases, pediatric cases of coronavirus disease 2019, and healthy controls. We define a MIS-C transcriptional signature partially shared with the transcriptional response to SARS-CoV-2 infection and with the signature of Kawasaki disease, a clinically similar condition. By projecting the MIS-C signature onto a co-expression network, we identified disease gene modules and found genes downregulated in MIS-C clustered in a module enriched for the transcriptional signatures of exhausted CD8 + T-cells and CD56 dim CD57 + NK cells. Bayesian network analyses revealed nine key regulators of this module, including TBX21 , a central coordinator of exhausted CD8 + T-cell differentiation. Together, these findings suggest dysregulated cytotoxic lymphocyte response to SARS-Cov-2 infection in MIS-C.

Multiscale causal networks identify VGF as a key regulator of Alzheimer's disease

Though discovered over 100 years ago, the molecular foundation of sporadic Alzheimer's disease (AD) remains elusive. To better characterize the complex nature of AD, we constructed multiscale causal networks on a large human AD multi-omics dataset, integrating clinical features of AD, DNA variation, and gene- and protein-expression. These probabilistic causal models enabled detection, prioritization and replication of high-confidence master regulators of AD-associated networks, including the top predicted regulator, VGF. Overexpression of neuropeptide precursor VGF in 5xFAD mice partially rescued beta-amyloid-mediated memory impairment and neuropathology. Molecular validation of network predictions downstream of VGF was also achieved in this AD model, with significant enrichment for homologous genes identified as differentially expressed in 5xFAD brains overexpressing VGF. Our findings support a causal role for VGF in protecting against AD pathogenesis and progression.

An inflammatory cytokine signature helps predict COVID-19 severity and death

The COVID-19 pandemic caused by infection with Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has led to more than 100,000 deaths in the United States. Several studies have revealed that the hyper-inflammatory response induced by SARS-CoV-2 is a major cause of disease severity and death in infected patients. However, predictive biomarkers of pathogenic inflammation to help guide targetable immune pathways are critically lacking. We implemented a rapid multiplex cytokine assay to measure serum IL-6, IL-8, TNF-α, and IL-1β in hospitalized COVID-19 patients upon admission to the Mount Sinai Health System in New York. Patients (n=1484) were followed up to 41 days (median 8 days) and clinical information, laboratory test results and patient outcomes were collected. In 244 patients, cytokine measurements were repeated over time, and effect of drugs could be assessed. Kaplan-Meier methods were used to compare survival by cytokine strata, followed by Cox regression models to evaluate the independent predictive value of baseline cytokines. We found that high serum IL-6, IL-8, and TNF-α levels at the time of hospitalization were strong and independent predictors of patient survival. Importantly, when adjusting for disease severity score, common laboratory inflammation markers, hypoxia and other vitals, demographics, and a range of comorbidities, IL-6 and TNF-α serum levels remained independent and significant predictors of disease severity and death. We propose that serum IL-6 and TNF-α levels should be considered in the management and treatment of COVID-19 patients to stratify prospective clinical trials, guide resource allocation and inform therapeutic options. We also propose that patients with high IL-6 and TNF-α levels should be assessed for combinatorial blockade of pathogenic inflammation in this disease.

Late-onset megaconial myopathy in mice lacking group I Paks

BACKGROUND

Group I Paks are serine/threonine kinases that function as major effectors of the small GTPases Rac1 and Cdc42, and they regulate cytoskeletal dynamics, cell polarity, and transcription. We previously demonstrated that Pak1 and Pak2 function redundantly to promote skeletal myoblast differentiation during postnatal development and regeneration in mice. However, the roles of Pak1 and Pak2 in adult muscle homeostasis are unknown. Choline kinase β (Chk β) is important for adult muscle homeostasis, as autosomal recessive mutations in CHKβ are associated with two human muscle diseases, megaconial congenital muscular dystrophy and proximal myopathy with focal depletion of mitochondria.

METHODS

We analyzed mice conditionally lacking Pak1 and Pak2 in the skeletal muscle lineage (double knockout (dKO) mice) over 1 year of age. Muscle integrity in dKO mice was assessed with histological stains, immunofluorescence, electron microscopy, and western blotting. Assays for mitochondrial respiratory complex function were performed, as was mass spectrometric quantification of products of choline kinase. Mice and cultured myoblasts deficient for choline kinase β (Chk β) were analyzed for Pak1/2 phosphorylation.

RESULTS

dKO mice developed an age-related myopathy. By 10 months of age, dKO mouse muscles displayed centrally-nucleated myofibers, fibrosis, and signs of degeneration. Disease severity occurred in a rostrocaudal gradient, hindlimbs more strongly affected than forelimbs. A distinctive feature of this myopathy was elongated and branched intermyofibrillar (megaconial) mitochondria, accompanied by focal mitochondrial depletion in the central region of the fiber. dKO muscles showed reduced mitochondrial respiratory complex I and II activity. These phenotypes resemble those of rmd mice, which lack Chkβ and are a model for human diseases associated with CHKβ deficiency. Pak1/2 and Chkβ activities were not interdependent in mouse skeletal muscle, suggesting a more complex relationship in regulation of mitochondria and muscle homeostasis.

CONCLUSIONS

Conditional loss of Pak1 and Pak2 in mice resulted in an age-dependent myopathy with similarity to mice and humans with CHKβ deficiency. Protein kinases are major regulators of most biological processes but few have been implicated in muscle maintenance or disease. Pak1/Pak2 dKO mice offer new insights into these processes.

The Mount Sinai cohort of large-scale genomic, transcriptomic and proteomic data in Alzheimer's disease

Alzheimer's disease (AD) affects half the US population over the age of 85 and is universally fatal following an average course of 10 years of progressive cognitive disability. Genetic and genome-wide association studies (GWAS) have identified about 33 risk factor genes for common, late-onset AD (LOAD), but these risk loci fail to account for the majority of affected cases and can neither provide clinically meaningful prediction of development of AD nor offer actionable mechanisms. This cohort study generated large-scale matched multi-Omics data in AD and control brains for exploring novel molecular underpinnings of AD. Specifically, we generated whole genome sequencing, whole exome sequencing, transcriptome sequencing and proteome profiling data from multiple regions of 364 postmortem control, mild cognitive impaired (MCI) and AD brains with rich clinical and pathophysiological data. All the data went through rigorous quality control. Both the raw and processed data are publicly available through the Synapse software platform.

The human brainome: network analysis identifies HSPA2 as a novel Alzheimer&rsquo;s disease target

Our hypothesis is that changes in gene and protein expression are crucial to the development of late-onset Alzheimer’s disease. Previously we examined how DNA alleles control downstream expression of RNA transcripts and how those relationships are changed in late-onset Alzheimer’s disease. We have now examined how proteins are incorporated into networks in two separate series and evaluated our outputs in two different cell lines. Our pipeline included the following steps: (i) predicting expression quantitative trait loci; (ii) determining differential expression; (iii) analysing networks of transcript and peptide relationships; and (iv) validating effects in two separate cell lines. We performed all our analysis in two separate brain series to validate effects. Our two series included 345 samples in the first set (177 controls, 168 cases; age range 65–105; 58% female; KRONOSII cohort) and 409 samples in the replicate set (153 controls, 141 cases, 115 mild cognitive impairment; age range 66–107; 63% female; RUSH cohort). Our top target is heat shock protein family A member 2 (HSPA2), which was identified as a key driver in our two datasets. HSPA2 was validated in two cell lines, with overexpression driving further elevation of amyloid-β40 and amyloid-β42 levels in APP mutant cells, as well as significant elevation of microtubule associated protein tau and phosphorylated-tau in a modified neuroglioma line. This work further demonstrates that studying changes in gene and protein expression is crucial to understanding late onset disease and further nominates HSPA2 as a specific key regulator of late-onset Alzheimer’s disease processes.10.1093/brain/awy215_video1awy215media15824729224001.

Multiscale Analysis of Independent Alzheimer's Cohorts Finds Disruption of Molecular, Genetic, and Clinical Networks by Human Herpesvirus

Investigators have long suspected that pathogenic microbes might contribute to the onset and progression of Alzheimer's disease (AD) although definitive evidence has not been presented. Whether such findings represent a causal contribution, or reflect opportunistic passengers of neurodegeneration, is also difficult to resolve. We constructed multiscale networks of the late-onset AD-associated virome, integrating genomic, transcriptomic, proteomic, and histopathological data across four brain regions from human post-mortem tissue. We observed increased human herpesvirus 6A (HHV-6A) and human herpesvirus 7 (HHV-7) from subjects with AD compared with controls. These results were replicated in two additional, independent and geographically dispersed cohorts. We observed regulatory relationships linking viral abundance and modulators of APP metabolism, including induction of APBB2, APPBP2, BIN1, BACE1, CLU, PICALM, and PSEN1 by HHV-6A. This study elucidates networks linking molecular, clinical, and neuropathological features with viral activity and is consistent with viral activity constituting a general feature of AD.

Comparative analyses of population-scale phenomic data in electronic medical records reveal race-specific disease networks

Motivation: Underrepresentation of racial groups represents an important challenge and major gap in phenomics research. Most of the current human phenomics research is based primarily on European populations; hence it is an important challenge to expand it to consider other population groups. One approach is to utilize data from EMR databases that contain patient data from diverse demographics and ancestries. The implications of this racial underrepresentation of data can be profound regarding effects on the healthcare delivery and actionability. To the best of our knowledge, our work is the first attempt to perform comparative, population-scale analyses of disease networks across three different populations, namely Caucasian (EA), African American (AA) and Hispanic/Latino (HL).

Results: We compared susceptibility profiles and temporal connectivity patterns for 1988 diseases and 37 282 disease pairs represented in a clinical population of 1 025 573 patients. Accordingly, we revealed appreciable differences in disease susceptibility, temporal patterns, network structure and underlying disease connections between EA, AA and HL populations. We found 2158 significantly comorbid diseases for the EA cohort, 3265 for AA and 672 for HL. We further outlined key disease pair associations unique to each population as well as categorical enrichments of these pairs. Finally, we identified 51 key ‘hub’ diseases that are the focal points in the race-centric networks and of particular clinical importance. Incorporating race-specific disease comorbidity patterns will produce a more accurate and complete picture of the disease landscape overall and could support more precise understanding of disease relationships and patient management towards improved clinical outcomes.

Contacts: rong.chen@mssm.edu or joel.dudley@mssm.edu

Supplementary information:Supplementary data are available at Bioinformatics online.

Meditation and vacation effects have an impact on disease-associated molecular phenotypes

Meditation is becoming increasingly practiced, especially for stress-related medical conditions. Meditation may improve cellular health; however, studies have not separated out effects of meditation from vacation-like effects in a residential randomized controlled trial. We recruited healthy women non-meditators to live at a resort for 6 days and randomized to either meditation retreat or relaxing on-site, with both groups compared with 'regular meditators' already enrolled in the retreat. Blood drawn at baseline and post intervention was assessed for transcriptome-wide expression patterns and aging-related biomarkers. Highly significant gene expression changes were detected across all groups (the 'vacation effect') that could accurately predict (96% accuracy) between baseline and post-intervention states and were characterized by improved regulation of stress response, immune function and amyloid beta (Aβ) metabolism. Although a smaller set of genes was affected, regular meditators showed post-intervention differences in a gene network characterized by lower regulation of protein synthesis and viral genome activity. Changes in well-being were assessed post intervention relative to baseline, as well as 1 and 10 months later. All groups showed equivalently large immediate post-intervention improvements in well-being, but novice meditators showed greater maintenance of lower distress over time compared with those in the vacation arm. Regular meditators showed a trend toward increased telomerase activity compared with randomized women, who showed increased plasma Aβ42/Aβ40 ratios and tumor necrosis factor alpha (TNF-α) levels. This highly controlled residential study showed large salutary changes in gene expression networks due to the vacation effect, common to all groups. For those already trained in the practice of meditation, a retreat appears to provide additional benefits to cellular health beyond the vacation effect.

A new molecular signature method for prediction of driver cancer pathways from transcriptional data

Assigning cancer patients to the most effective treatments requires an understanding of the molecular basis of their disease. While DNA-based molecular profiling approaches have flourished over the past several years to transform our understanding of driver pathways across a broad range of tumors, a systematic characterization of key driver pathways based on RNA data has not been undertaken. Here we introduce a new approach for predicting the status of driver cancer pathways based on signature functions derived from RNA sequencing data. To identify the driver cancer pathways of interest, we mined DNA variant data from TCGA and nominated driver alterations in seven major cancer pathways in breast, ovarian and colon cancer tumors. The activation status of these driver pathways were then characterized using RNA sequencing data by constructing classification signature functions in training datasets and then testing the accuracy of the signatures in test datasets. The signature functions differentiate well tumors with nominated pathway activation from tumors with no signs of activation: average AUC equals to 0.83. Our results confirm that driver genomic alterations are distinctively displayed at the transcriptional level and that the transcriptional signatures can generally provide an alternative to DNA sequencing methods in detecting specific driver pathways.

Detecting epigenetic motifs in low coverage and metagenomics settings

Background

It has recently become possible to rapidly and accurately detect epigenetic signatures in bacterial genomes using third generation sequencing data. Monitoring the speed at which a single polymerase inserts a base in the read strand enables one to infer whether a modification is present at that specific site on the template strand. These sites can be challenging to detect in the absence of high coverage and reliable reference genomes.

Methods

Here we provide a new method for detecting epigenetic motifs in bacteria on datasets with low-coverage, with incomplete references, and with mixed samples (i.e. metagenomic data). Our approach treats motif inference as a kmer comparison problem. First, genomes (or contigs) are deconstructed into kmers. Then, native genome-wide distributions of interpulse durations (IPDs) for kmers are compared with corresponding whole genome amplified (WGA, modification free) IPD distributions using log likelihood ratios. Finally, kmers are ranked and greedily selected by iteratively correcting for sequences within a particular kmer's neighborhood.

Conclusions

Our method can detect multiple types of modifications, even at very low-coverage and in the presence of mixed genomes. Additionally, we are able to predict modified motifs when genomes with "neighbor" modified motifs exist within the sample. Lastly, we show that these motifs can provide an alternative source of information by which to cluster metagenomics contigs and that iterative refinement on these clustered contigs can further improve both sensitivity and specificity of motif detection.

Availability

https://github.com/alibashir/EMMCKmer

A 600 kb deletion syndrome at 16p11.2 leads to energy imbalance and neuropsychiatric disorders

Background

The recurrent ∼600 kb 16p11.2 BP4-BP5 deletion is among the most frequent known genetic aetiologies of autism spectrum disorder (ASD) and related neurodevelopmental disorders.

Objective

To define the medical, neuropsychological, and behavioural phenotypes in carriers of this deletion.

Methods

We collected clinical data on 285 deletion carriers and performed detailed evaluations on 72 carriers and 68 intrafamilial non-carrier controls.

Results

When compared to intrafamilial controls, full scale intelligence quotient (FSIQ) is two standard deviations lower in carriers, and there is no difference between carriers referred for neurodevelopmental disorders and carriers identified through cascade family testing. Verbal IQ (mean 74) is lower than non-verbal IQ (mean 83) and a majority of carriers require speech therapy. Over 80% of individuals exhibit psychiatric disorders including ASD, which is present in 15% of the paediatric carriers. Increase in head circumference (HC) during infancy is similar to the HC and brain growth patterns observed in idiopathic ASD. Obesity, a major comorbidity present in 50% of the carriers by the age of 7 years, does not correlate with FSIQ or any behavioural trait. Seizures are present in 24% of carriers and occur independently of other symptoms. Malformations are infrequently found, confirming only a few of the previously reported associations.

Conclusions

The 16p11.2 deletion impacts in a quantitative and independent manner FSIQ, behaviour and body mass index, possibly through direct influences on neural circuitry. Although non-specific, these features are clinically significant and reproducible. Lastly, this study demonstrates the necessity of studying large patient cohorts ascertained through multiple methods to characterise the clinical consequences of rare variants involved in common diseases.

Genetic evidence for a role of adiponutrin in the metabolism of apolipoprotein B-containing lipoproteins

Adiponutrin (PNPLA3) is a predominantly liver-expressed transmembrane protein with phospholipase activity that is regulated by fasting and feeding. Recent genome-wide association studies identified PNPLA3 to be associated with hepatic fat content and liver function, thus pointing to a possible involvement in the hepatic lipoprotein metabolism. The aim of this study was to examine the association between two common variants in the adiponutrin gene and parameters of lipoprotein metabolism in 23,274 participants from eight independent West-Eurasian study populations including six population-based studies [Bruneck (n = 800), KORA S3/F3 (n = 1644), KORA S4/F4 (n = 1814), CoLaus (n = 5435), SHIP (n = 4012), Rotterdam (n = 5967)], the SAPHIR Study as a healthy working population (n = 1738) and the Utah Obesity Case-Control Study including a group of 1037 severely obese individuals (average BMI 46 kg/m2) and 827 controls from the same geographical region of Utah. We observed a strong additive association of a common non-synonymous variant within adiponutrin (rs738409) with age-, gender-, and alanine-aminotransferase-adjusted lipoprotein concentrations: each copy of the minor allele decreased levels of total cholesterol on average by 2.43 mg/dl (P = 8.87 x 10(-7)), non-HDL cholesterol levels by 2.35 mg/dl (P = 2.27 x 10(-6)) and LDL cholesterol levels by 1.48 mg/dl (P = 7.99 x 10(-4)). These associations remained significant after correction for multiple testing. We did not observe clear evidence for associations with HDL cholesterol or triglyceride concentrations. In conclusion, our study suggests that adiponutrin is involved in the metabolism of apoB-containing lipoproteins.

Population-based genome-wide association studies reveal six loci influencing plasma levels of liver enzymes

Plasma liver-enzyme tests are widely used in the clinic for the diagnosis of liver diseases and for monitoring the response to drug treatment. There is considerable evidence that human genetic variation influences plasma levels of liver enzymes. However, such genetic variation has not been systematically assessed. In the present study, we performed a genome-wide association study of plasma liver-enzyme levels in three populations (total n = 7715) with replication in three additional cohorts (total n = 4704). We identified two loci influencing plasma levels of alanine-aminotransferase (ALT) (CPN1-ERLIN1-CHUK on chromosome 10 and PNPLA3-SAMM50 on chromosome 22), one locus influencing gamma-glutamyl transferase (GGT) levels (HNF1A on chromosome 12), and three loci for alkaline phosphatase (ALP) levels (ALPL on chromosome 1, GPLD1 on chromosome 6, and JMJD1C-REEP3 on chromosome 10). In addition, we confirmed the associations between the GGT1 locus and GGT levels and between the ABO locus and ALP levels. None of the ALP-associated SNPs were associated with other liver tests, suggesting intestine and/or bone specificity. The mechanisms underlying the associations may involve cis- or trans-transcriptional effects (some of the identified variants were associated with mRNA transcription in human liver or lymphoblastoid cells), dysfunction of the encoded proteins (caused by missense variations at the functional domains), or other unknown pathways. These findings may help in the interpretation of liver-enzyme tests and provide candidate genes for liver diseases of viral, metabolic, autoimmune, or toxic origin. The specific associations with ALP levels may point to genes for bone or intestinal diseases.

Population-based genome-wide association studies reveal six loci influencing plasma levels of liver enzymes

Plasma liver-enzyme tests are widely used in the clinic for the diagnosis of liver diseases and for monitoring the response to drug treatment. There is considerable evidence that human genetic variation influences plasma levels of liver enzymes. However, such genetic variation has not been systematically assessed. In the present study, we performed a genome-wide association study of plasma liver-enzyme levels in three populations (total n = 7715) with replication in three additional cohorts (total n = 4704). We identified two loci influencing plasma levels of alanine-aminotransferase (ALT) (CPN1-ERLIN1-CHUK on chromosome 10 and PNPLA3-SAMM50 on chromosome 22), one locus influencing gamma-glutamyl transferase (GGT) levels (HNF1A on chromosome 12), and three loci for alkaline phosphatase (ALP) levels (ALPL on chromosome 1, GPLD1 on chromosome 6, and JMJD1C-REEP3 on chromosome 10). In addition, we confirmed the associations between the GGT1 locus and GGT levels and between the ABO locus and ALP levels. None of the ALP-associated SNPs were associated with other liver tests, suggesting intestine and/or bone specificity. The mechanisms underlying the associations may involve cis- or trans-transcriptional effects (some of the identified variants were associated with mRNA transcription in human liver or lymphoblastoid cells), dysfunction of the encoded proteins (caused by missense variations at the functional domains), or other unknown pathways. These findings may help in the interpretation of liver-enzyme tests and provide candidate genes for liver diseases of viral, metabolic, autoimmune, or toxic origin. The specific associations with ALP levels may point to genes for bone or intestinal diseases.