Glucocorticoid receptor regulates protein chaperone, circadian clock and affective disorder genes in the zebrafish brain.

Glucocorticoid resistance is commonly observed in depression, and has been linked to reduced expression and/or function of the glucocorticoid receptor (NR3C1 in human, hereafter referred to as GR). Previous studies have shown that GR-mutant zebrafish exhibit behavioural abnormalities that are indicative of an affective disorder, suggesting that GR plays a role in brain function. We compared the brain methylomes and brain transcriptomes of adult wild-type and GR-mutant zebrafish, and identified 249 differentially methylated regions (DMRs) that are regulated by GR. These include a cluster of CpG sites within the first intron of fkbp5, the gene encoding the glucocorticoid-inducible heat shock protein co-chaperone Fkbp5. RNA-sequencing analysis revealed that genes associated with chaperone-mediated protein folding, the regulation of circadian rhythm and the regulation of metabolism are particularly sensitive to loss of GR function. In addition, we identified subsets of genes exhibiting GR-regulated transcription that are known to regulate behaviour, and are linked to unipolar depression and anxiety. Taken together, our results identify key biological processes and novel molecular mechanisms through which the GR is likely to mediate responses to stress in the adult zebrafish brain, and they provide further support for the zebrafish GR mutant as a model for the study of affective disorders.

Plasma Cell-Free DNA Predicts Survival and Maps Specific Sources of Injury in Pulmonary Arterial Hypertension.

BACKGROUND: Cell-free DNA (cfDNA) is a noninvasive marker of cellular injury. Its significance in pulmonary arterial hypertension (PAH) is unknown. METHODS: Plasma cfDNA was measured in 2 PAH cohorts (A, n=48; B, n=161) and controls (n=48). Data were collected for REVEAL 2.0 (Registry to Evaluate Early and Long-Term PAH Disease Management) scores and outcome determinations. Patients were divided into the following REVEAL risk groups: low (≤6), medium (7-8), and high (≥9). Total cfDNA concentrations were compared among controls and PAH risk groups by 1-way analysis of variance. Log-rank tests compared survival between cfDNA tertiles and REVEAL risk groups. Areas under the receiver operating characteristic curve were estimated from logistic regression models. A sample subset from cohort B (n=96) and controls (n=16) underwent bisulfite sequencing followed by a deconvolution algorithm to map cell-specific cfDNA methylation patterns, with concentrations compared using t tests. RESULTS: In cohort A, median (interquartile range) age was 62 years (47-71), with 75% female, and median (interquartile range) REVEAL 2.0 was 6 (4-9). In cohort B, median (interquartile range) age was 59 years (49-71), with 69% female, and median (interquartile range) REVEAL 2.0 was 7 (6-9). In both cohorts, cfDNA concentrations differed among patients with PAH of varying REVEAL risk and controls (analysis of variance P≤0.002) and were greater in the high-risk compared with the low-risk category (P≤0.002). In cohort B, death or lung transplant occurred in 14 of 54, 23 of 53, and 35 of 54 patients in the lowest, middle, and highest cfDNA tertiles, respectively. cfDNA levels stratified as tertiles (log-rank: P=0.0001) and REVEAL risk groups (log-rank: P<0.0001) each predicted transplant-free survival. The addition of cfDNA to REVEAL improved discrimination (area under the receiver operating characteristic curve, 0.72-0.78; P=0.02). Compared with controls, methylation analysis in patients with PAH revealed increased cfDNA originating from erythrocyte progenitors, neutrophils, monocytes, adipocytes, natural killer cells, vascular endothelium, and cardiac myocytes (Bonferroni adjusted P<0.05). cfDNA concentrations derived from erythrocyte progenitor cells, cardiac myocytes, and vascular endothelium were greater in patients with PAH with high-risk versus low-risk REVEAL scores (P≤0.02). CONCLUSIONS: Circulating cfDNA is elevated in patients with PAH, correlates with disease severity, and predicts worse survival. Results from cfDNA methylation analyses in patients with PAH are consistent with prevailing paradigms of disease pathogenesis.

Genetic variants of PKLR are associated with acute pain in sickle cell disease.

Acute pain, the most prominent complication of sickle cell disease (SCD), results from vaso-occlusion triggered by sickling of deoxygenated red blood cells (RBCs). Concentration of 2,3-diphosphoglycerate (2,3-DPG) in RBCs promotes deoxygenation by preferentially binding to the low-affinity T conformation of HbS. 2,3-DPG is an intermediate substrate in the glycolytic pathway in which pyruvate kinase (gene PKLR, protein PKR) is a rate-limiting enzyme; variants in PKLR may affect PKR activity, 2,3-DPG levels in RBCs, RBC sickling, and acute pain episodes (APEs). We performed a candidate gene association study using 2 cohorts: 242 adult SCD-HbSS patients and 977 children with SCD-HbSS or SCD-HbSß0 thalassemia. Seven of 47 PKLR variants evaluated in the adult cohort were associated with hospitalization: intron 4, rs2071053; intron 2, rs8177970, rs116244351, rs114455416, rs12741350, rs3020781, and rs8177964. All 7 variants showed consistent effect directions in both cohorts and remained significant in weighted Fisher's meta-analyses of the adult and pediatric cohorts using P < .0071 as threshold to correct for multiple testing. Allele-specific expression analyses in an independent cohort of 52 SCD adults showed that the intronic variants are likely to influence APE by affecting expression of PKLR, although the causal variant and mechanism are not defined.

Amygdala and anterior cingulate transcriptomes from individuals with bipolar disorder reveal downregulated neuroimmune and synaptic pathways.

Recent genetic studies have identified variants associated with bipolar disorder (BD), but it remains unclear how brain gene expression is altered in BD and how genetic risk for BD may contribute to these alterations. Here, we obtained transcriptomes from subgenual anterior cingulate cortex and amygdala samples from post-mortem brains of individuals with BD and neurotypical controls, including 511 total samples from 295 unique donors. We examined differential gene expression between cases and controls and the transcriptional effects of BD-associated genetic variants. We found two coexpressed modules that were associated with transcriptional changes in BD: one enriched for immune and inflammatory genes and the other with genes related to the postsynaptic membrane. Over 50% of BD genome-wide significant loci contained significant expression quantitative trait loci (QTL) (eQTL), and these data converged on several individual genes, including SCN2A and GRIN2A. Thus, these data implicate specific genes and pathways that may contribute to the pathology of BP.

Integrated cell-free DNA and cytokine analysis uncovers distinct tissue injury and immune response patterns in solid organ transplant recipients with COVID-19.

COVID-19 pathogenesis is associated with an exuberant inflammatory response. However, the tissue injury pattern and immune response in solid-organ transplant recipients (SOTRs) taking immunosuppressive therapy have not been well characterized. Here, we perform both cfDNA and cytokine profiling on plasma samples to map tissue damage, including allograft injury and delineate underlying immunopathology. We identified injuries from multiple-tissue types, including hematopoietic cells, vascular endothelium, hepatocyte, adipocyte, pancreas, kidney, heart, and lung in SOTRs with COVID-19 that correlates with disease severity. SOTRs with COVID-19 have higher plasma levels of cytokines such as IFN-λ1, IFN-γ, IL-15, IL-18 IL-1RA, IL-6, MCP-2, and TNF-α as compared to healthy controls, and the levels of GM-CSF, IL-15, IL-6, IL-8, and IL-10 were associated with disease severity in SOTRs. Strikingly, IFN-λ and IP-10 were markedly increased in SOTRs compared to immunocompetent patients with COVID-19. Correlation analyses showed a strong association between monocyte-derived cfDNA and inflammatory cytokines/chemokines in SOTRs with COVID-19. Moreover, compared to other respiratory viral infections, COVID-19 induced pronounced injury in hematopoietic, vascular endothelial and endocrine tissues. Allograft injury, measured as donor-derived cfDNA was elevated in SOTRs with COVID-19, including allografts distant from the primary site of infection. Allograft injury correlated with inflammatory cytokines and cfDNA from immune cells. Furthermore, longitudinal analysis identified a gradual decrease of cfDNA and inflammatory cytokine levels in patients with a favorable outcome. Our findings highlight distinct tissue injury and cytokine features in SOTRs with COVID-19 that correlate with disease severity.

Early Myeloid Derived Suppressor Cells (eMDSCs) Are Associated With High Donor Myeloid Chimerism Following Haploidentical HSCT for Sickle Cell Disease.

Haploidentical hematopoietic stem cell transplantation (haplo-HSCT) is a widely available curative option for patients with sickle cell disease (SCD). Our original non-myeloablative haplo-HSCT trial employing post-transplant (PT) cyclophosphamide had a low incidence of GVHD but had high rejection rates. Here, we aimed to evaluate immune reconstitution following haplo-HSCT and identify cytokines and cells associated with graft rejection/engraftment. 50 cytokines and 10 immune cell subsets were screened using multiplex-ELISA and flow cytometry, respectively, at baseline and PT-Days 30, 60, 100, and 180. We observed the most significant differences in cytokine levels between the engrafted and rejected groups at PT-Day 60, corresponding with clinical findings of secondary graft rejection. Of the 44 cytokines evaluated, plasma concentrations of 19 cytokines were different between the two groups at PT-Day 60. Factor analysis suggested two independent factors. The first factor (IL-17A, IL-10, IL-7, G-CSF, IL-2, MIP-1a, VEGF, and TGFb1 contributed significantly) was strongly associated with engraftment with OR = 2.7 (95%CI of 1.4 to 5.4), whereas the second factor (GROa and IL-18 contributed significantly) was not significantly associated with engraftment. Sufficient donor myeloid chimerism (DMC) is critical for the success of HSCT; here, we evaluated immune cells among high (H) DMC (DMC≥20%) and low (L) DMC (DMC<20%) groups along with engrafted and rejected groups. We found that early myeloid-derived suppressor cell (eMDSC) frequencies were elevated in engrafted patients and patients with HDMC at PT-Day 30 (P< 0.04 & P< 0.003, respectively). 9 of 20 patients were evaluated for the source of eMDSCs. The HDMC group had high mixed chimeric eMDSCs as compared to the LDMC group (P< 0.00001). We found a positive correlation between the frequencies of eMDSCs and Tregs at PT-Day 100 (r=0.72, P <0.0007); eMDSCs at BSL and Tregs at PT-Day 100 (r=0.63, P <0.004). Of 10 immune regulatory cells and 50 cytokines, we observed mixed chimeric eMDSCs and IL-17A, IL-10, IL-7, G-CSF, IL-2, MIP-1a, VEGF, TGFb1 as potential hits which could serve as prognostic markers in predicting allograft outcome towards engraftment following haploidentical HSCT employing post-transplant cyclophosphamide. The current findings need to be replicated and further explored in a larger cohort.

NOTCH-mediated ex vivo expansion of human hematopoietic stem and progenitor cells by culture under hypoxia.

Activation of NOTCH signaling in human hematopoietic stem/progenitor cells (HSPCs) by treatment with an engineered Delta-like ligand (DELTA1ext-IgG [DXI]) has enabled ex vivo expansion of short-term HSPCs, but the effect on long-term repopulating hematopoietic stem cells (LTR-HSCs) remains uncertain. Here, we demonstrate that ex vivo culture of human adult HSPCs with DXI under low oxygen tension limits ER stress in LTR-HSCs and lineage-committed progenitors compared with normoxic cultures. A distinct HSC gene signature was upregulated in cells cultured with DXI in hypoxia and, after 21 days of culture, the frequency of LTR-HSCs increased 4.9-fold relative to uncultured cells and 4.2-fold compared with the normoxia + DXI group. NOTCH and hypoxia pathways intersected to maintain undifferentiated phenotypes in cultured HSPCs. Our work underscores the importance of mitigating ER stress perturbations to preserve functional LTR-HSCs in extended cultures and offers a clinically feasible platform for the expansion of human HSPCs.

Methylomic and transcriptomic predictors of one-month exposure to cortisol in healthy individuals.

Allostatic load (AL) refers to the cumulative "wear and tear" on an organism throughout its lifetime. One of the primary contributing factors to AL is prolonged exposure to stress or its primary catabolic agent cortisol. Chronic exposure to stress or cortisol is associated with numerous diseases, including cardiovascular disease, metabolic disorders, and psychiatric disorders. Therefore, a molecular marker capable of integrating a past history of cortisol exposure would be of great utility for assessing disease risk. To this end, we recruited 87 healthy males and females of European ancestry between 18 and 60 years old, extracted genomic DNA and RNA from leukocytes, and implemented a gene-centric DNA enrichment method coupled with bisulfite sequencing and RNA-Seq of total RNA for the determination of genome-wide methylation and gene transcription, respectively. Sequencing data were analyzed against awakening and bedtime cortisol data to identify differentially methylated regions (DMRs) and CpGs (DMCs) and differentially expressed genes (DEGs). Six candidate DMCs (punadjusted < 0.005) and nine DEGs (punadjusted < 0.0005) were used to construct a prediction model that could capture past 30+ days of both bedtime and awakening cortisol levels. Utilizing a cross-validation approach, we obtained a regression coefficient of R2 = 0.308 for predicting continuous awakening cortisol and an area under the curve (AUC) = 0.753 for dichotomous (high vs. low tertile) awakening cortisol, and R2 = 0.224 and AUC = 0.723 for continuous and dichotomous bedtime cortisol levels, respectively. To our knowledge, the current study represents the first attempt to identify genome-wide predictors of cortisol exposure that utilizes both methylation and transcription targets. The utility of our approach needs to be replicated in an independent cohort of samples for which similar cortisol metrics are available.

Network Analysis and Transcriptome Profiling Identify Autophagic and Mitochondrial Dysfunctions in SARS-CoV-2 Infection.

Analyzing host cells' transcriptional response to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection will help delineate biological processes underlying viral pathogenesis. First, analysis of expression profiles of lung cell lines A549 and Calu3 revealed upregulation of antiviral interferon signaling genes in response to all three SARS-CoV-2, MERS-CoV, or influenza A virus (IAV) infections. However, perturbations in expression of genes involved in inflammatory, mitochondrial, and autophagy processes were specifically observed in SARS-CoV-2-infected cells. Next, a validation study in infected human nasopharyngeal samples also revealed perturbations in autophagy and mitochondrial processes. Specifically, mTOR expression, mitochondrial ribosomal, mitochondrial complex I, lysosome acidification, and mitochondrial fission promoting genes were concurrently downregulated in both infected cell lines and human samples. SARS-CoV-2 infection impeded autophagic flux either by upregulating GSK3B in lung cell lines or by downregulating autophagy genes, SNAP29, and lysosome acidification genes in human samples, contributing to increased viral replication. Therefore, drugs targeting lysosome acidification or autophagic flux could be tested as intervention strategies. Finally, age-stratified SARS-CoV-2-positive human data revealed impaired upregulation of chemokines, interferon-stimulated genes, and tripartite motif genes that are critical for antiviral signaling. Together, this analysis has revealed specific aspects of autophagic and mitochondrial function that are uniquely perturbed in SARS-CoV-2-infected host cells.

Fasting-induced FOXO4 blunts human CD4+ T helper cell responsiveness.

Intermittent fasting blunts inflammation in asthma1 and rheumatoid arthritis2, suggesting that fasting may be exploited as an immune-modulatory intervention. However, the mechanisms underpinning the anti-inflammatory effects of fasting are poorly characterized3-5. Here, we show that fasting in humans is sufficient to blunt CD4+ T helper cell responsiveness. RNA sequencing and flow cytometry immunophenotyping of peripheral blood mononuclear cells from volunteers subjected to overnight or 24-h fasting and 3 h of refeeding suggest that fasting blunts CD4+ T helper cell activation and differentiation. Transcriptomic analysis reveals that longer fasting has a more robust effect on CD4+ T-cell biology. Through bioinformatics analyses, we identify the transcription factor FOXO4 and its canonical target FK506-binding protein 5 (FKBP5) as a potential fasting-responsive regulatory axis. Genetic gain- or loss-of-function of FOXO4 and FKBP5 is sufficient to modulate TH1 and TH17 cytokine production. Moreover, we find that fasting-induced or genetic overexpression of FOXO4 and FKBP5 is sufficient to downregulate mammalian target of rapamycin complex 1 signalling and suppress signal transducer and activator of transcription 1/3 activation. Our results identify FOXO4-FKBP5 as a new fasting-induced, signal transducer and activator of transcription-mediated regulatory pathway to blunt human CD4+ T helper cell responsiveness.

Circulating mitochondrial DNA is a proinflammatory DAMP in sickle cell disease.

The pathophysiology of sickle cell disease (SCD) is driven by chronic inflammation fueled by damage associated molecular patterns (DAMPs). We show that elevated cell-free DNA (cfDNA) in patients with SCD is not just a prognostic biomarker, it also contributes to the pathological inflammation. Within the elevated cfDNA, patients with SCD had a significantly higher ratio of cell-free mitochondrial DNA (cf-mtDNA)/cell-free nuclear DNA compared with healthy controls. Additionally, mitochondrial DNA in patient samples showed significantly disproportionately increased hypomethylation compared with healthy controls, and it was increased further in crises compared with steady-state. Using flow cytometry, structured illumination microscopy, and electron microscopy, we showed that circulating SCD red blood cells abnormally retained their mitochondria and, thus, are likely to be the source of the elevated cf-mtDNA in patients with SCD. Patient plasma containing high levels of cf-mtDNA triggered the formation of neutrophil extracellular traps (NETs) that was substantially reduced by inhibition of TANK-binding kinase 1, implicating activation of the cGAS-STING pathway. cf-mtDNA is an erythrocytic DAMP, highlighting an underappreciated role for mitochondria in sickle pathology. These trials were registered at www.clinicaltrials.gov as #NCT00081523, #NCT03049475, and #NCT00047996.

Cell-free DNA maps COVID-19 tissue injury and risk of death and can cause tissue injury.

INTRODUCTIONThe clinical course of coronavirus 2019 (COVID-19) is heterogeneous, ranging from mild to severe multiorgan failure and death. In this study, we analyzed cell-free DNA (cfDNA) as a biomarker of injury to define the sources of tissue injury that contribute to such different trajectories.METHODSWe conducted a multicenter prospective cohort study to enroll patients with COVID-19 and collect plasma samples. Plasma cfDNA was subject to bisulfite sequencing. A library of tissue-specific DNA methylation signatures was used to analyze sequence reads to quantitate cfDNA from different tissue types. We then determined the correlation of tissue-specific cfDNA measures to COVID-19 outcomes. Similar analyses were performed for healthy controls and a comparator group of patients with respiratory syncytial virus and influenza.RESULTSWe found markedly elevated levels and divergent tissue sources of cfDNA in COVID-19 patients compared with patients who had influenza and/or respiratory syncytial virus and with healthy controls. The major sources of cfDNA in COVID-19 were hematopoietic cells, vascular endothelium, hepatocytes, adipocytes, kidney, heart, and lung. cfDNA levels positively correlated with COVID-19 disease severity, C-reactive protein, and D-dimer. cfDNA profile at admission identified patients who subsequently required intensive care or died during hospitalization. Furthermore, the increased cfDNA in COVID-19 patients generated excessive mitochondrial ROS (mtROS) in renal tubular cells in a concentration-dependent manner. This mtROS production was inhibited by a TLR9-specific antagonist.CONCLUSIONcfDNA maps tissue injury that predicts COVID-19 outcomes and may mechanistically propagate COVID-19-induced tissue injury.FUNDINGIntramural Targeted Anti-COVID-19 grant, NIH.

Identification of human long noncoding RNAs associated with nonalcoholic fatty liver disease and metabolic homeostasis.

A growing number of long noncoding RNAs (lncRNAs) have emerged as vital metabolic regulators. However, most human lncRNAs are nonconserved and highly tissue specific, vastly limiting our ability to identify human lncRNA metabolic regulators (hLMRs). In this study, we established a pipeline to identify putative hLMRs that are metabolically sensitive, disease relevant, and population applicable. We first progressively processed multilevel human transcriptome data to select liver lncRNAs that exhibit highly dynamic expression in the general population, show differential expression in a nonalcoholic fatty liver disease (NAFLD) population, and respond to dietary intervention in a small NAFLD cohort. We then experimentally demonstrated the responsiveness of selected hepatic lncRNAs to defined metabolic milieus in a liver-specific humanized mouse model. Furthermore, by extracting a concise list of protein-coding genes that are persistently correlated with lncRNAs in general and NAFLD populations, we predicted the specific function for each hLMR. Using gain- and loss-of-function approaches in humanized mice as well as ectopic expression in conventional mice, we validated the regulatory role of one nonconserved hLMR in cholesterol metabolism by coordinating with an RNA-binding protein, PTBP1, to modulate the transcription of cholesterol synthesis genes. Our work overcame the heterogeneity intrinsic to human data to enable the efficient identification and functional definition of disease-relevant human lncRNAs in metabolic homeostasis.

Bioinformatics Approaches for Functional Prediction of Long Noncoding RNAs.

There is accumulating evidence that long noncoding RNAs (lncRNAs) play crucial roles in biological processes and diseases. In recent years, computational models have been widely used to predict potential lncRNA-disease relations. In this chapter, we systematically describe various computational algorithms and prediction tools that have been developed to elucidate the roles of lncRNAs in diseases, coding potential/functional characterization, or ascertaining their involvement in critical biological processes as well as provide a comprehensive summary of these applications.

Genome-Wide Analysis of Off-Target CRISPR/Cas9 Activity in Single-Cell-Derived Human Hematopoietic Stem and Progenitor Cell Clones.

CRISPR/Cas9 (clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9)-mediated genome editing holds remarkable promise for the treatment of human genetic diseases. However, the possibility of off-target Cas9 activity remains a concern. To address this issue using clinically relevant target cells, we electroporated Cas9 ribonucleoprotein (RNP) complexes (independently targeted to two different genomic loci, the CXCR4 locus on chromosome 2 and the AAVS1 locus on chromosome 19) into human mobilized peripheral blood-derived hematopoietic stem and progenitor cells (HSPCs) and assessed the acquisition of somatic mutations in an unbiased, genome-wide manner via whole genome sequencing (WGS) of single-cell-derived HSPC clones. Bioinformatic analysis identified >20,000 total somatic variants (indels, single nucleotide variants, and structural variants) distributed among Cas9-treated and non-Cas9-treated control HSPC clones. Statistical analysis revealed no significant difference in the number of novel non-targeted indels among the samples. Moreover, data analysis showed no evidence of Cas9-mediated indel formation at 623 predicted off-target sites. The median number of novel single nucleotide variants was slightly elevated in Cas9 RNP-recipient sample groups compared to baseline, but did not reach statistical significance. Structural variants were rare and demonstrated no clear causal connection to Cas9-mediated gene editing procedures. We find that the collective somatic mutational burden observed within Cas9 RNP-edited human HSPC clones is indistinguishable from naturally occurring levels of background genetic heterogeneity.

Acetylation-mediated remodeling of the nucleolus regulates cellular acetyl-CoA responses.

The metabolite acetyl-coenzyme A (acetyl-CoA) serves as an essential element for a wide range of cellular functions including adenosine triphosphate (ATP) production, lipid synthesis, and protein acetylation. Intracellular acetyl-CoA concentrations are associated with nutrient availability, but the mechanisms by which a cell responds to fluctuations in acetyl-CoA levels remain elusive. Here, we generate a cell system to selectively manipulate the nucleo-cytoplasmic levels of acetyl-CoA using clustered regularly interspaced short palindromic repeat (CRISPR)-mediated gene editing and acetate supplementation of the culture media. Using this system and quantitative omics analyses, we demonstrate that acetyl-CoA depletion alters the integrity of the nucleolus, impairing ribosomal RNA synthesis and evoking the ribosomal protein-dependent activation of p53. This nucleolar remodeling appears to be mediated through the class IIa histone deacetylases (HDACs). Our findings highlight acetylation-mediated control of the nucleolus as an important hub linking acetyl-CoA fluctuations to cellular stress responses.

lncRNAKB, a knowledgebase of tissue-specific functional annotation and trait association of long noncoding RNA.

Long non-coding RNA Knowledgebase (lncRNAKB) is an integrated resource for exploring lncRNA biology in the context of tissue-specificity and disease association. A systematic integration of annotations from six independent databases resulted in 77,199 human lncRNA (224,286 transcripts). The user-friendly knowledgebase covers a comprehensive breadth and depth of lncRNA annotation. lncRNAKB is a compendium of expression patterns, derived from analysis of RNA-seq data in thousands of samples across 31 solid human normal tissues (GTEx). Thousands of co-expression modules identified via network analysis and pathway enrichment to delineate lncRNA function are also accessible. Millions of expression quantitative trait loci (cis-eQTL) computed using whole genome sequence genotype data (GTEx) can be downloaded at lncRNAKB that also includes tissue-specificity, phylogenetic conservation and coding potential scores. Tissue-specific lncRNA-trait associations encompassing 323 GWAS (UK Biobank) are also provided. LncRNAKB is accessible at http://www.lncrnakb.org/ , and the data are freely available through Open Science Framework ( https://doi.org/10.17605/OSF.IO/RU4D2 ).

Network Analysis and Transcriptome Profiling Identify Autophagic and Mitochondrial Dysfunctions in SARS-CoV-2 Infection.

Analyzing host transcriptional changes in response to SARS-CoV-2 infection will help delineate biological processes underlying viral pathogenesis. Comparison of expression profiles of lung cell lines A549 (infected with either SARS-CoV-2 (with ACE2 expression)) or Influenza A virus (IAV)) and Calu3 (infected with SARS-CoV-2 or MERS-CoV) revealed upregulation of the antiviral interferon signaling in all three viral infections. However, perturbations in inflammatory, mitochondrial, and autophagy processes were specifically observed in SARS-CoV-2 infected cells. Validation of findings from cell line data revealed perturbations in autophagy and mitochondrial processes in the infected human nasopharyngeal samples. Specifically, downregulation of mTOR expression, mitochondrial ribosomal, mitochondrial complex I, and lysosome acidification genes were concurrently observed in both infected cell lines and human datasets. Furthermore, SARS-CoV-2 infection impedes autophagic flux by upregulating GSK3B in lung cell lines, or by downregulating autophagy genes, SNAP29 and lysosome acidification genes in human samples, contributing to increased viral replication. Therefore, drugs targeting lysosome acidification or autophagic flux could be tested as intervention strategies. Additionally, downregulation of MTFP1 (in cell lines) or SOCS6 (in human samples) results in hyperfused mitochondria and impede proper interferon response. Coexpression networks analysis identifies correlated clusters of genes annotated to inflammation and mitochondrial processes that are misregulated in SARS-CoV-2 infected cells. Finally, comparison of age stratified human gene expression data revealed impaired upregulation of chemokines, interferon stimulated and tripartite motif genes that are critical for antiviral signaling. Together, this analysis has revealed specific aspects of autophagic and mitochondrial function that are uniquely perturbed in SARS-CoV-2 infection.

ßT87Q-Globin Gene Therapy Reduces Sickle Hemoglobin Production, Allowing for Ex Vivo Anti-sickling Activity in Human Erythroid Cells.

Lentiviral addition of ßT87Q-globin, a modified ß-globin with an anti-sickling mutation, is currently being used in gene therapy trials for sickle cell disease (SCD) and ß-thalassemia patients. ßT87Q-globin interferes with sickle hemoglobin (HbS) polymerization. Here, we generated the SCD mutation in an immortalized human erythroid cell line (HUDEP-2) to investigate the anti-sickling activity of ßT87Q-globin. Sickle HUDEP-2 (sHUDEP-2) cells produced robust HbS after differentiation and sickled under deoxygenated conditions, comparable with SCD CD34+ progeny. Lentiviral transduction provided 9.5-26.8 pg/cell ßT87Q-globin (R2 = 0.83) in a vector copy number (VCN)-dependent manner, resulting in a significant reduction of sickling ratios (R2 = 0.92). Interestingly, ßT87Q-globin transduction markedly reduced endogenous ßS-globin (R2 = 0.84) to an undetectable level (0.4-16.8 pg/cell) in sHUDEP-2 cells, as well as endogenous ß-globin in human CD34+ cell-derived erythroid cells. RNA sequencing (RNA-seq) analysis with ßT87Q-transduced sHUDEP-2 and human CD34+-derived cells revealed activation of inflammation- and proliferation-related programs, suggesting minimal changes in background gene expression except for ßT87Q-globin expression and endogenous ß/ßS-globin suppression. In summary, using sHUDEP-2 and CD34+-derived cells, we demonstrated that lentiviral addition of ßT87Q-globin strongly reduced endogenous ß-/ßS-globin expression, resulting in an anti-sickling effect. Our findings should be helpful to understand the anti-sickling effects of therapeutic genes in SCD gene therapy.

Apolipoprotein E Signals via TLR4 to Induce CXCL5 Secretion by Asthmatic Airway Epithelial Cells.

The primary function of APOE (apolipoprotein E) is to mediate the transport of cholesterol- and lipid-containing lipoprotein particles into cells by receptor-mediated endocytosis. APOE also has pro- and antiinflammatory effects, which are both context and concentration dependent. For example, Apoe-/- mice exhibit enhanced airway remodeling and hyperreactivity in experimental asthma, whereas increased APOE levels in lung epithelial lining fluid induce IL-1ß secretion from human asthmatic alveolar macrophages. However, APOE-mediated airway epithelial cell inflammatory responses and signaling pathways have not been defined. Here, RNA sequencing of human asthmatic bronchial brushing cells stimulated with APOE identified increased expression of mRNA transcripts encoding multiple proinflammatory genes, including CXCL5 (C-X-C motif chemokine ligand 5), an epithelial-derived chemokine that promotes neutrophil activation and chemotaxis. We subsequently characterized the APOE signaling pathway that induces CXCL5 secretion by human asthmatic small airway epithelial cells (SAECs). Neutralizing antibodies directed against TLR4 (Toll-like receptor 4), but not TLR2, attenuated APOE-mediated CXCL5 secretion by human asthmatic SAECs. Inhibition of TAK1 (transforming growth factor-ß-activated kinase 1), IκKß (inhibitor of nuclear factor κ B kinase subunit ß), TPL2 (tumor progression locus 2), and JNK (c-Jun N-terminal kinase), but not p38 MAPK (mitogen-activated protein kinase) or MEK1/2 (MAPK kinase 1/2), attenuated APOE-mediated CXCL5 secretion. The roles of TAK1, IκKß, TPL2, and JNK in APOE-mediated CXCL5 secretion were verified by RNA interference. Furthermore, RNA interference showed that after APOE stimulation, both NF-κB p65 and TPL2 were downstream of TAK1 and IκKß, whereas JNK was downstream of TPL2. In summary, elevated levels of APOE in the airway may activate a TLR4/TAK1/IκKß/NF-κB/TPL2/JNK signaling pathway that induces CXCL5 secretion by human asthmatic SAECs. These findings identify new roles for TLR4 and TPL2 in APOE-mediated proinflammatory responses in asthma.