GediNET for discovering gene associations across diseases using knowledge based machine learning approach.

The most common approaches to discovering genes associated with specific diseases are based on machine learning and use a variety of feature selection techniques to identify significant genes that can serve as biomarkers for a given disease. More recently, the integration in this process of prior knowledge-based approaches has shown significant promise in the discovery of new biomarkers with potential translational applications. In this study, we developed a novel approach, GediNET, that integrates prior biological knowledge to gene Groups that are shown to be associated with a specific disease such as a cancer. The novelty of GediNET is that it then also allows the discovery of significant associations between that specific disease and other diseases. The initial step in this process involves the identification of gene Groups. The Groups are then subjected to a Scoring component to identify the top performing classification Groups. The top-ranked gene Groups are then used to train a Machine Learning Model. The process of Grouping, Scoring and Modelling (G-S-M) is used by GediNET to identify other diseases that are similarly associated with this signature. GediNET identifies these relationships through Disease-Disease Association (DDA) based machine learning. DDA explores novel associations between diseases and identifies relationships which could be used to further improve approaches to diagnosis, prognosis, and treatment. The GediNET KNIME workflow can be downloaded from: https://github.com/malikyousef/GediNET.git or https://kni.me/w/3kH1SQV_mMUsMTS .

Telomere dysfunction implicates POT1 in patients with idiopathic pulmonary fibrosis.

Exonic sequencing identified a family with idiopathic pulmonary fibrosis (IPF) containing a previously unreported heterozygous mutation in POT1 p.(L259S). The family displays short telomeres and genetic anticipation. We found that POT1(L259S) is defective in binding the telomeric overhang, nuclear accumulation, negative regulation of telomerase, and lagging strand maintenance. Patient cells containing the mutation display telomere loss, lagging strand defects, telomere-induced DNA damage, and premature senescence with G1 arrest. Our data suggest POT1(L259S) is a pathogenic driver of IPF and provide insights into gene therapy options.

Functional Annotation of MicroRNAs Using Existing Resources.

MicroRNAs (miRNAs) are endogenous small noncoding RNAs that are involved in most biological signaling pathways, including the cell cycle, apoptosis, proliferation, immune response, metabolism as well as in biological processes including organ development and in human diseases like cancers. During the past two decades, high-throughput transcriptomic profiling using next generation sequencing and microarrays have been extensively utilized to identify differentially expressed miRNAs across different conditions and diseases. A natural extension of miRNA identification is to the process of functionally annotating known or predicted gene targets of those miRNAs and, by inference, revealing their potential influences on diverse biological pathways and functions. In this chapter, we provide a stepwise guideline on how to perform functional enrichment analyses on miRNAs of interest using publicly available resources such as miRWalk2.0.

Lack of Atorvastatin Effect on Monocyte Gene Expression and Inflammatory Markers in HIV-1-infected ART-suppressed Individuals at Risk of non-AIDS Comorbidities.

BACKGROUND: Many people living with HIV have persistent monocyte activation despite viral suppression by antiretroviral therapy (ART), which contributes to non-AIDS complications including neurocognitive and other disorders. Statins have immunomodulatory properties that might be beneficial by reducing monocyte activation. METHODS: We previously characterized monocyte gene expression and inflammatory markers in 11 HIV-positive individuals on long-term ART (HIV/ART) at risk for non-AIDS complications because of low nadir CD4+ counts (median 129 cells/uL) and elevated hsCRP. Here, these individuals participated in a double-blind, randomized, placebo-controlled crossover study of 12 weeks of atorvastatin treatment. Monocyte surface markers were assessed by flow cytometry, plasma mediators by ELISA and Luminex, and monocyte gene expression by microarray analysis. RESULTS: Among primary outcome measures, 12 weeks of atorvastatin treatment led to an unexpected increase in CCR2+ monocytes (P=0.04), but did not affect CD16+ or CD163+ monocytes, nor levels in plasma of CCL2/MCP-1 or sCD14. Among secondary outcomes, atorvastatin treatment was associated with decreased plasma hsCRP (P=0.035) and IL-2R (P=0.012). Treatment was also associated with increased total CD14+ monocytes (P=0.015), and increased plasma CXCL9 (P=0.003) and IL-12 (P<0.001). Comparable results were seen in a subgroup that had inflammatory marker elevations at baseline. Atorvastatin treatment did not significantly alter monocyte gene expression or normalize aberrant baseline transcriptional patterns. CONCLUSIONS: In this study of aviremic HIV+ individuals at high risk of non-AIDS events, 12 weeks of atorvastatin did not normalize monocyte gene expression patterns nor lead to significant changes in monocyte surface markers or plasma mediators linked to non-AIDS comorbidities.

Clinical presentation of COVID-19 - a model derived by a machine learning algorithm.

COVID-19 pandemic has flooded all triage stations, making it difficult to carefully select those most likely infected. Data on total patients tested, infected, and hospitalized is fragmentary making it difficult to easily select those most likely to be infected. The Israeli Ministry of Health made public its registry of immediate clinical data and the respective status of infected/not infected for all viral DNA tests performed up to Apr. 18th, 2020 including almost 120,000 tests. We used a machine-learning algorithm to find out which immediate clinical elements mattered the most in identifying the true status of the tested persons including age or gender matter, to enable future better allocation of surveillance policy for those belonging to high-risk groups. In addition to the analyses applied on the first batch of the available data (Apr. 11th), we further tested the algorithm on the independent second batch (Apr. 12th to 18th). Fever, cough and headache were the most diagnostic, differing in degree of importance in different subgroups. Higher percentage of men were found positive (9.3 vs. 7.3%), but gender did not matter for the clinical presentation. The prediction power of the model was high, with accuracy of 0.84 and area under the curve 0.92. We provide a hand-held short checklist with verbal description of importance for the leading symptoms, which should expedite the triage and enable proper selection of people for further follow-up.

A multi-omics approach to Epstein-Barr virus immortalization of B-cells reveals EBNA1 chromatin pioneering activities targeting nucleotide metabolism.

Epstein-Barr virus (EBV) immortalizes resting B-lymphocytes through a highly orchestrated reprogramming of host chromatin structure, transcription and metabolism. Here, we use a multi-omics-based approach to investigate these underlying mechanisms. ATAC-seq analysis of cellular chromatin showed that EBV alters over a third of accessible chromatin during the infection time course, with many of these sites overlapping transcription factors such as PU.1, Interferon Regulatory Factors (IRFs), and CTCF. Integration of RNA-seq analysis identified a complex transcriptional response and associations with EBV nuclear antigens (EBNAs). Focusing on EBNA1 revealed enhancer-binding activity at gene targets involved in nucleotide metabolism, supported by metabolomic analysis which indicated that adenosine and purine metabolism are significantly altered by EBV immortalization. We further validated that adenosine deaminase (ADA) is a direct and critical target of the EBV-directed immortalization process. These findings reveal that purine metabolism and ADA may be useful therapeutic targets for EBV-driven lymphoid cancers.

Methyltransferase inhibitors restore SATB1 protective activity against cutaneous T cell lymphoma in mice.

Cutaneous T cell lymphoma (CTCL) has a poorly understood etiology and no known cure. Using conditional knockout mice, we found that ablation of the genomic organizer special AT-rich sequence-binding protein 1 (Satb1) caused malignant transformation of mature, skin-homing, Notch-activated CD4+ and CD8+ T cells into progressively fatal lymphoma. Mechanistically, Satb1 restrained Stat5 phosphorylation and the expression of skin-homing chemokine receptors in mature T cells. Notably, methyltransferase-dependent epigenetic repression of SATB1 was universally found in human Sézary syndrome, but not in other peripheral T cell malignancies. H3K27 and H3K9 trimethylation occluded the SATB1 promoter in Sézary cells, while inhibition of SUV39H1/2 methyltransferases (unlike EZH2 inhibition) restored protective SATB1 expression and selectively abrogated the growth of primary Sézary cells more effectively than romidepsin. Therefore, inhibition of methyltransferases that silence SATB1 could address an unmet need for patients with mycosis fungoides/Sézary syndrome, a set of incurable diseases.

PAF-R on activated T cells: Role in the IL-23/Th17 pathway and relevance to multiple sclerosis.

IL-23 is a potent stimulus for Th17 cells. These cells have a distinct developmental pathway from Th1 cells induced by IL-12 and are implicated in autoimmune and inflammatory disorders including multiple sclerosis (MS). TGF-ß, IL-6, and IL-1, the transcriptional regulator RORγt (RORC) and IL-23 are implicated in Th17 development and maintenance. In human polyclonally activated T cells, IL-23 enhances IL-17 production. The aims of our study were: 1). To validate microarray results showing preferential expression of platelet activating factor receptor (PAF-R) on IL-23 stimulated T cells. 2). To determine whether PAF-R on activated T cells is functional, whether it is co-regulated with Th17-associated molecules, and whether it is implicated in Th17 function. 3). To determine PAF-R expression in MS. We show that PAF-R is expressed on activated T cells, and is inducible by IL-23 and IL-17, which in turn are induced by PAF binding to PAF-R. PAF-R is co-expressed with IL-17 and regulated similarly with Th17 markers IL-17A, IL-17F, IL-22 and RORC. PAF-R is upregulated on PBMC and T cells of MS patients, and levels correlate with IL-17 and with MS disability scores. Our results show that PAF-R on T cells is associated with the Th17 phenotype and function. Clinical Implications Targeting PAF-R may interfere with Th17 function and offer therapeutic intervention in Th17-associated conditions, including MS.

Epigenetic specifications of host chromosome docking sites for latent Epstein-Barr virus.

Epstein-Barr virus (EBV) genomes persist in latently infected cells as extrachromosomal episomes that attach to host chromosomes through the tethering functions of EBNA1, a viral encoded sequence-specific DNA binding protein. Here we employ circular chromosome conformation capture (4C) analysis to identify genome-wide associations between EBV episomes and host chromosomes. We find that EBV episomes in Burkitt's lymphoma cells preferentially associate with cellular genomic sites containing EBNA1 binding sites enriched with B-cell factors EBF1 and RBP-jK, the repressive histone mark H3K9me3, and AT-rich flanking sequence. These attachment sites correspond to transcriptionally silenced genes with GO enrichment for neuronal function and protein kinase A pathways. Depletion of EBNA1 leads to a transcriptional de-repression of silenced genes and reduction in H3K9me3. EBV attachment sites in lymphoblastoid cells with different latency type show different correlations, suggesting that host chromosome attachment sites are functionally linked to latency type gene expression programs.

Recursive Cluster Elimination based Rank Function (SVM-RCE-R) implemented in KNIME.

In our earlier study, we proposed a novel feature selection approach, Recursive Cluster Elimination with Support Vector Machines (SVM-RCE) and implemented this approach in Matlab. Interest in this approach has grown over time and several researchers have incorporated SVM-RCE into their studies, resulting in a substantial number of scientific publications. This increased interest encouraged us to reconsider how feature selection, particularly in biological datasets, can benefit from considering the relationships of those genes in the selection process, this led to our development of SVM-RCE-R. The usefulness of SVM-RCE-R is further supported by development of maTE tool, which uses a similar approach to identify microRNA (miRNA) targets. We have now implemented the SVM-RCE-R algorithm in Knime in order to make it easier to apply and to make it more accessible to the biomedical community. The use of SVM-RCE-R in Knime is simple and intuitive, allowing researchers to immediately begin their data analysis without having to consult an information technology specialist. The input for the Knime tool is an EXCEL file (or text or CSV) with a simple structure and the output is also an EXCEL file. The Knime version also incorporates new features not available in the previous version. One of these features is a user-specific ranking function that enables the user to provide the weights of the accuracy, sensitivity, specificity, f-measure, area under curve and precision in the ranking function, allowing the user to select for greater sensitivity or greater specificity as needed. The results show that the ranking function has an impact on the performance of SVM-RCE-R. Some of the clusters that achieve high scores for a specified ranking can also have high scores in other metrics. This finding motivates future studies to suggest the optimal ranking function.

Involvement of condensin in cellular senescence through gene regulation and compartmental reorganization.

Senescence is induced by various stimuli such as oncogene expression and telomere shortening, referred to as oncogene-induced senescence (OIS) and replicative senescence (RS), respectively, and accompanied by global transcriptional alterations and 3D genome reorganization. Here, we demonstrate that the human condensin II complex participates in senescence via gene regulation and reorganization of euchromatic A and heterochromatic B compartments. Both OIS and RS are accompanied by A-to-B and B-to-A compartmental transitions, the latter of which occur more frequently and are undergone by 14% (430 Mb) of the human genome. Mechanistically, condensin is enriched in A compartments and implicated in B-to-A transitions. The full activation of senescence genes (SASP genes and p53 targets) requires condensin; its depletion impairs senescence markers. This study describes that condensin reinforces euchromatic A compartments and promotes B-to-A transitions, both of which are coupled to optimal expression of senescence genes, thereby allowing condensin to contribute to senescent processes.

NK Response Correlates with HIV Decrease in Pegylated IFN-α2a-Treated Antiretroviral Therapy-Suppressed Subjects.

We previously reported that pegylated IFN-α2a (Peg-IFN-α2a) added to antiretroviral therapy (ART)-suppressed, HIV-infected subjects resulted in plasma HIV control and integrated HIV DNA decrease. We now evaluated whether innate NK cell activity or PBMC transcriptional profiles were associated with decreases in HIV measures. Human peripheral blood was analyzed prior to Peg-IFN-α2a administration (ART, baseline), after 5 wk of ART+Peg-IFN-α2a, and after 12 wk of Peg-IFN-α2a monotherapy (primary endpoint). After 5 wk of ART+Peg-IFN-α2a, immune subset frequencies were preserved, and induction of IFN-stimulated genes was noted in all subjects except for a subset in which the lack of IFN-stimulated gene induction was associated with increased expression of microRNAs. Viral control during Peg-IFN-α2a monotherapy was associated with 1) higher levels of NK cell activity and IFN-γ-induced protein 10 (IP-10) on ART (preimmunotherapy) and 2) downmodulation of NK cell KIR2DL1 and KIR2DL2/DL3 expression, transcriptional enrichment of expression of genes associated with NK cells in HIV controller subjects, and higher ex vivo IFN-α-induced NK cytotoxicity after 5 wk of ART+Peg-IFN-α2a. Integrated HIV DNA decline after immunotherapy was also associated with gene expression patterns indicative of cell-mediated activation and NK cytotoxicity. Overall, an increase in innate activity and NK cell cytotoxicity were identified as correlates of Peg-IFN-α2a-mediated HIV control.

N6-Methylation of Adenosine of FZD10 mRNA Contributes to PARP Inhibitor Resistance.

Despite the high initial response rates to PARP inhibitors (PARPi) in BRCA-mutated epithelial ovarian cancers (EOC), PARPi resistance remains a major challenge. Chemical modifications of RNAs have emerged as a new layer of epigenetic gene regulation. N6-methyladenosine (m6A) is the most abundant chemical modification of mRNA, yet the role of m6A modification in PARPi resistance has not previously been explored. Here, we show that m6A modification of FZD10 mRNA contributes to PARPi resistance by upregulating the Wnt/ß-catenin pathway in BRCA-mutated EOC cells. Global m6A profile revealed a significant increase in m6A modification in FZD10 mRNA, which correlated with increased FZD10 mRNA stability and an upregulation of the Wnt/ß-catenin pathway. Depletion of FZD10 or inhibition of the Wnt/ß-catenin sensitizes resistant cells to PARPi. Mechanistically, downregulation of m6A demethylases FTO and ALKBH5 was sufficient to increase FZD10 mRNA m6A modification and reduce PARPi sensitivity, which correlated with an increase in homologous recombination activity. Moreover, combined inhibition of PARP and Wnt/ß-catenin showed synergistic suppression of PARPi-resistant cells in vitro and in vivo in a xenograft EOC mouse model. Overall, our results show that m6A contributes to PARPi resistance in BRCA-deficient EOC cells by upregulating the Wnt/ß-catenin pathway via stabilization of FZD10. They also suggest that inhibition of the Wnt/ß-catenin pathway represents a potential strategy to overcome PARPi resistance. SIGNIFICANCE: These findings elucidate a novel regulatory mechanism of PARPi resistance in EOC by showing that m6A modification of FZD10 mRNA contributes to PARPi resistance in BRCA-deficient EOC cells via upregulation of Wnt/ß-catenin pathway.

Evidence for Persistent Monocyte and Immune Dysregulation After Prolonged Viral Suppression Despite Normalization of Monocyte Subsets, sCD14 and sCD163 in HIV-Infected Individuals.

BACKGROUND: People living with HIV on antiretroviral therapy (HIV/ART) experience excess non-AIDS comorbidities, and also remain at increased risk for certain infections and viral malignancies. Monocytes/macrophages are central to many of these comorbidities, and elevated plasma cytokines and immune activation during untreated infection are often incompletely reversed by ART and are also associated with comorbidities. METHODS: We investigated monocyte surface markers, gene expression, and plasma cytokines in 11 HIV-infected older individuals (median 53 years) who started therapy with low CD4 counts (median 129 cells/µl), with elevated hsCRP (≥ 2mg/L) despite long-term ART (median 7.4 years), along with matched controls. RESULTS: Frequency of monocyte subsets (based on CD14/CD16/CD163), were not different from controls, but surface expression of CD163 was increased (P = 0.021) while PD1 was decreased (P = 0.013) along with a trend for higher tissue factor (P = 0.096). As a group, HIV/ART participants had elevated plasma CCL2 (MCP-1; P = 0.0001), CXCL9 (MIG; P = 0.04), and sIL2R (P = 0.015), which were correlated, while sCD14 was not elevated. Principal component analysis of soluble markers revealed that 6/11 HIV/ART participants clustered with controls, while 5 formed a distinct group, driven by IL-10, CCL11, CXCL10, CCL2, CXCL9, and sIL2R. These individuals were significantly older than those clustering with controls. Transcriptomic analysis revealed multiple genes linked to immune functions including inflammation, immune cell development, and cell-cell signaling that were downregulated in HIV/ART monocytes and distinct from patterns in untreated subjects. CONCLUSIONS: Long-term ART-treated individuals normalize monocyte subsets but exhibit immune dysregulation involving both aberrant inflammation and monocyte dysfunction, as well as inter-individual heterogeneity, suggesting complex mechanisms linking monocytes and HIV/ART comorbidities.

A Gene Expression Classifier from Whole Blood Distinguishes Benign from Malignant Lung Nodules Detected by Low-Dose CT.

Low-dose CT (LDCT) is widely accepted as the preferred method for detecting pulmonary nodules. However, the determination of whether a nodule is benign or malignant involves either repeated scans or invasive procedures that sample the lung tissue. Noninvasive methods to assess these nodules are needed to reduce unnecessary invasive tests. In this study, we have developed a pulmonary nodule classifier (PNC) using RNA from whole blood collected in RNA-stabilizing PAXgene tubes that addresses this need. Samples were prospectively collected from high-risk and incidental subjects with a positive lung CT scan. A total of 821 samples from 5 clinical sites were analyzed. Malignant samples were predominantly stage 1 by pathologic diagnosis and 97% of the benign samples were confirmed by 4 years of follow-up. A panel of diagnostic biomarkers was selected from a subset of the samples assayed on Illumina microarrays that achieved a ROC-AUC of 0.847 on independent validation. The microarray data were then used to design a biomarker panel of 559 gene probes to be validated on the clinically tested NanoString nCounter platform. RNA from 583 patients was used to assess and refine the NanoString PNC (nPNC), which was then validated on 158 independent samples (ROC-AUC = 0.825). The nPNC outperformed three clinical algorithms in discriminating malignant from benign pulmonary nodules ranging from 6-20 mm using just 41 diagnostic biomarkers. Overall, this platform provides an accurate, noninvasive method for the diagnosis of pulmonary nodules in patients with non-small cell lung cancer. SIGNIFICANCE: These findings describe a minimally invasive and clinically practical pulmonary nodule classifier that has good diagnostic ability at distinguishing benign from malignant pulmonary nodules.

The effect of timing of influenza vaccination and sample collection on antibody titers and responses in the aged.

Antibody responses, B cell subset distribution in blood and the blood transcriptome were analyzed in younger and aged human subjects before and after vaccination with the inactivated influenza vaccine. In the aged, but not the younger, individuals we saw a clear difference in antibody titers including those at baseline depending on the time of vaccination and sample collection. Differences in baseline titers in aged individuals treated in the morning or afternoon in turn affected responsiveness to the vaccine. In both younger and aged individuals, the time of sample collection also affected relative numbers of some of the B cell subsets in blood. A global gene expression analysis with whole blood samples from the aged showed small but statistically significant differences depending on the time of sample collection. Our data do not indicate that timing of vaccination affects immune responsiveness of the aged, but rather shows that in clinical influenza vaccine trials timing of collection of samples can have a major and potentially misleading influence on study outcome. In future vaccine trials, timing of vaccination and sample collection should be recorded carefully to allow for its use as a study covariant.

Structural and functional analysis of the human POT1-TPP1 telomeric complex.

POT1 and TPP1 are part of the shelterin complex and are essential for telomere length regulation and maintenance. Naturally occurring mutations of the telomeric POT1-TPP1 complex are implicated in familial glioma, melanoma and chronic lymphocytic leukaemia. Here we report the atomic structure of the interacting portion of the human telomeric POT1-TPP1 complex and suggest how several of these mutations contribute to malignant cancer. The POT1 C-terminus (POT1C) forms a bilobal structure consisting of an OB-fold and a holiday junction resolvase domain. TPP1 consists of several loops and helices involved in extensive interactions with POT1C. Biochemical data shows that several of the cancer-associated mutations, partially disrupt the POT1-TPP1 complex, which affects its ability to bind telomeric DNA efficiently. A defective POT1-TPP1 complex leads to longer and fragile telomeres, which in turn promotes genomic instability and cancer.

ADAR1 controls apoptosis of stressed cells by inhibiting Staufen1-mediated mRNA decay.

Both p150 and p110 isoforms of ADAR1 convert adenosine to inosine in double-stranded RNA (dsRNA). ADAR1p150 suppresses the dsRNA-sensing mechanism that activates MDA5-MAVS-IFN signaling in the cytoplasm. In contrast, the biological function of the ADAR1p110 isoform, which is usually located in the nucleus, is largely unknown. Here, we show that stress-activated phosphorylation of ADAR1p110 by MKK6-p38-MSK MAP kinases promotes its binding to Exportin-5 and its export from the nucleus. After translocating to the cytoplasm, ADAR1p110 suppresses apoptosis in stressed cells by protecting many antiapoptotic gene transcripts that contain 3'-untranslated-region dsRNA structures primarily comprising inverted Alu repeats. ADAR1p110 competitively inhibits binding of Staufen1 to the 3'-untranslated-region dsRNAs and antagonizes Staufen1-mediated mRNA decay. Our study reveals a new stress-response mechanism in which human ADAR1p110 and Staufen1 regulate surveillance of a set of mRNAs required for survival of stressed cells.

Age-related changes in the gene expression profile of antigen-specific mouse CD8+ T cells can be partially reversed by blockade of the BTLA/CD160 pathways during vaccination.

We analyzed gene expression profiles of young and aged mouse CD8+ T cells specific for the nucleoprotein (NP) of influenza A/PR8/34 virus. CD8+ T cells were stimulated either by the NP antigen expressed in its native form or fused into the herpes virus (HSV)-1 glycoprotein D (gD) protein, which blocks signaling through the immunoinhibitory B and T lymphocyte attenuator (BTLA) and CD160 pathways. We show that NP-specific CD8+ T cells from aged mice exhibit numerous differences in gene expression compared to NP-specific CD8+ T cells from young mice, including a significant reduction of expression in genes involved in T cell receptor (TcR) and CD28 signaling. We also show that these changes can be reversed in a sub-population (~50%) of the aged mice by a BTLA/CD160 checkpoint blockade. These results suggest that BTLA/CD160 checkpoint blockade has potential value as a vaccine additive to induce better CD8+ T cell responses in the aged.