Improved human disease candidate gene prioritization using mouse phenotype

IL36G-producing neutrophil-like monocytes promote cachexia in cancer

Most patients with advanced cancer develop cachexia, a multifactorial syndrome characterized by progressive skeletal muscle wasting. Despite its catastrophic impact on survival, the critical mediators responsible for cancer cachexia development remain poorly defined. Here, we show that a distinct subset of neutrophil-like monocytes, which we term cachexia-inducible monocytes (CiMs), emerges in the advanced cancer milieu and promotes skeletal muscle loss. Unbiased transcriptome analysis reveals that interleukin 36 gamma (IL36G)-producing CD38+ CiMs are induced in chronic monocytic blood cancer characterized by prominent cachexia. Notably, the emergence of CiMs and the activation of CiM-related gene signatures in monocytes are confirmed in various advanced solid cancers. Stimuli of toll-like receptor 4 signaling are responsible for the induction of CiMs. Genetic inhibition of IL36G-mediated signaling attenuates skeletal muscle loss and rescues cachexia phenotypes in advanced cancer models. These findings indicate that the IL36G-producing subset of neutrophil-like monocytes could be a potential therapeutic target in cancer cachexia.

Jumonji histone demethylases are therapeutic targets in small cell lung cancer

Small cell lung cancer (SCLC) is a recalcitrant cancer of neuroendocrine (NE) origin. Changes in therapeutic approaches against SCLC have been lacking over the decades. Here, we use preclinical models to identify a new therapeutic vulnerability in SCLC consisting of the targetable Jumonji lysine demethylase (KDM) family. We show that Jumonji demethylase inhibitors block malignant growth and that etoposide-resistant SCLC cell lines are particularly sensitive to Jumonji inhibition. Mechanistically, small molecule-mediated inhibition of Jumonji KDMs activates endoplasmic reticulum (ER) stress genes, upregulates ER stress signaling, and triggers apoptotic cell death. Furthermore, Jumonji inhibitors decrease protein levels of SCLC NE markers INSM1 and Secretogranin-3 and of driver transcription factors ASCL1 and NEUROD1. Genetic knockdown of KDM4A, a Jumonji demethylase highly expressed in SCLC and a known regulator of ER stress genes, induces ER stress response genes, decreases INSM1, Secretogranin-3, and NEUROD1 and inhibits proliferation of SCLC in vitro and in vivo. Lastly, we demonstrate that two different small molecule Jumonji KDM inhibitors (pan-inhibitor JIB-04 and KDM4 inhibitor SD70) block the growth of SCLC tumor xenografts in vivo. Our study highlights the translational potential of Jumonji KDM inhibitors against SCLC, a clinically feasible approach in light of recently opened clinical trials evaluating this drug class, and establishes KDM4A as a relevant target across SCLC subtypes.

Missense Mutations in Myc Box I Influence Nucleocytoplasmic Transport to Promote Leukemogenesis

PURPOSE

Somatic missense mutations in the phosphodegron domain of the MYC gene (MYC Box I or MBI) are detected in the dominant clones of a subset of patients with acute myeloid leukemia (AML), but the mechanisms by which they contribute to AML are unknown.

EXPERIMENTAL DESIGN

To investigate the effects of MBI MYC mutations on hematopoietic cells, we employed a multi-omic approach to systematically compare the cellular and molecular consequences of expressing oncogenic doses of wild type, threonine-58 and proline-59 mutant MYC proteins in hematopoietic cells, and we developed a knockin mouse harboring the germline MBI mutation p.T58N in the Myc gene.

RESULTS

Both wild-type and MBI mutant MYC proteins promote self-renewal programs and expand highly selected subpopulations of progenitor cells in the bone marrow. Compared with their wild-type counterparts, mutant cells display decreased cell death and accelerated leukemogenesis in vivo, changes that are recapitulated in the transcriptomes of human AML-bearing MYC mutations. The mutant phenotypes feature decreased stability and translation of mRNAs encoding proapoptotic and immune-regulatory genes, increased translation of RNA binding proteins and nuclear export machinery, and distinct nucleocytoplasmic RNA profiles. MBI MYC mutant proteins also show a higher propensity to aggregate in perinuclear regions and cytoplasm. Like the overexpression model, heterozygous p.T58N knockin mice displayed similar changes in subcellular MYC localization, progenitor expansion, transcriptional signatures, and develop hematopoietic tumors.

CONCLUSIONS

This study uncovers that MBI MYC mutations alter RNA nucleocytoplasmic transport mechanisms to contribute to the development of hematopoietic malignancies.

Multi-omics Analysis of Umbilical Cord Hematopoietic Stem Cells from a Multi-ethnic Cohort of Hawaii Reveals the Transgenerational Effect of Maternal Pre-Pregnancy Obesity

Background

Maternal obesity is a health concern that may predispose newborns to a high risk of medical problems later in life. To understand the transgenerational effect of maternal obesity, we conducted a multi-omics study, using DNA methylation and gene expression in the CD34+/CD38-/Lin- umbilical cord blood hematopoietic stem cells (uHSCs) and metabolomics of the cord blood, all from a multi-ethnic cohort (n=72) from Kapiolani Medical Center for Women and Children in Honolulu, Hawaii (collected between 2016 and 2018).

Results

Differential methylation (DM) analysis unveiled a global hypermethylation pattern in the maternal pre-pregnancy obese group (BH adjusted p<0.05), after adjusting for major clinical confounders. Comprehensive functional analysis showed hypermethylation in promoters of genes involved in cell cycle, protein synthesis, immune signaling, and lipid metabolism. Utilizing Shannon entropy on uHSCs methylation, we discerned notably higher quiescence of uHSCs impacted by maternal obesity. Additionally, the integration of multi-omics data-including methylation, gene expression, and metabolomics-provided further evidence of dysfunctions in adipogenesis, erythropoietin production, cell differentiation, and DNA repair, aligning with the findings at the epigenetic level. Furthermore, the CpG sites associated with maternal obesity from these pathways also predicted highly accurately (average AUC = 0.8687) between cancer vs. normal tissues in 14 cancer types in The Cancer Genome Atlas (TCGA).

Conclusions

This study revealed the significant correlation between pre-pregnancy maternal obesity and multi-omics level molecular changes in the uHSCs of offspring, particularly in DNA methylation. Moreover, these maternal obesity epigenetic markers in uHSCs may predispose offspring to higher cancer risks.

MECP2 directly interacts with RNA polymerase II to modulate transcription in human neurons

Mutations in the methyl-DNA-binding protein MECP2 cause the neurodevelopmental disorder Rett syndrome (RTT). How MECP2 contributes to transcriptional regulation in normal and disease states is unresolved; it has been reported to be an activator and a repressor. We describe here the first integrated CUT&Tag, transcriptome, and proteome analyses using human neurons with wild-type (WT) and mutant MECP2 molecules. MECP2 occupies CpG-rich promoter-proximal regions in over four thousand genes in human neurons, including a plethora of autism risk genes, together with RNA polymerase II (RNA Pol II). MECP2 directly interacts with RNA Pol II, and genes occupied by both proteins showed reduced expression in neurons with MECP2 patient mutations. We conclude that MECP2 acts as a positive cofactor for RNA Pol II gene expression at many neuronal genes that harbor CpG islands in promoter-proximal regions and that RTT is due, in part, to the loss of gene activity of these genes in neurons.

KMT2D regulates activation, localization, and integrin expression by T-cells

Individuals with Kabuki syndrome present with immunodeficiency; however, how pathogenic variants in the gene encoding the histone-modifying enzyme lysine methyltransferase 2D (KMT2D) lead to immune alterations remain poorly understood. Following up on our prior report of KMT2D-altered integrin expression in B-cells, we performed targeted analyses of KMT2D's influence on integrin expression in T-cells throughout development (thymocytes through peripheral T-cells) in murine cells with constitutive- and conditional-targeted Kmt2d deletion. Using high-throughput RNA-sequencing and flow cytometry, we reveal decreased expression (both at the transcriptional and translational levels) of a cluster of leukocyte-specific integrins, which perturb aspects of T-cell activation, maturation, adhesion/localization, and effector function. H3K4me3 ChIP-PCR suggests that these evolutionary similar integrins are under direct control of KMT2D. KMT2D loss also alters multiple downstream programming/signaling pathways, including integrin-based localization, which can influence T-cell populations. We further demonstrated that KMT2D deficiency is associated with the accumulation of murine CD8+ single-positive (SP) thymocytes and shifts in both human and murine peripheral T-cell populations, including the reduction of the CD4+ recent thymic emigrant (RTE) population. Together, these data show that the targeted loss of Kmt2d in the T-cell lineage recapitulates several distinct features of Kabuki syndrome-associated immune deficiency and implicates epigenetic mechanisms in the regulation of integrin signaling.

Characterization of Alternative Splicing in High-Risk Wilms' Tumors

The significant heterogeneity of Wilms' tumors between different patients is thought to arise from genetic and epigenetic distortions that occur during various stages of fetal kidney development in a way that is poorly understood. To address this, we characterized the heterogeneity of alternative mRNA splicing in Wilms' tumors using a publicly available RNAseq dataset of high-risk Wilms' tumors and normal kidney samples. Through Pareto task inference and cell deconvolution, we found that the tumors and normal kidney samples are organized according to progressive stages of kidney development within a triangle-shaped region in latent space, whose vertices, or "archetypes", resemble the cap mesenchyme, the nephrogenic stroma, and epithelial tubular structures of the fetal kidney. We identified a set of genes that are alternatively spliced between tumors located in different regions of latent space and found that many of these genes are associated with the epithelial-to-mesenchymal transition (EMT) and muscle development. Using motif enrichment analysis, we identified putative splicing regulators, some of which are associated with kidney development. Our findings provide new insights into the etiology of Wilms' tumors and suggest that specific splicing mechanisms in early stages of development may contribute to tumor development in different patients.

Epicardial deletion of Sox9 leads to myxomatous valve degeneration and identifies Cd109 as a novel gene associated with valve development

Epicardial-derived cells (EPDCs) are involved in the regulation of myocardial growth and coronary vascularization and are critically important for proper development of the atrioventricular (AV) valves. SOX9 is a transcription factor expressed in a variety of epithelial and mesenchymal cells in the developing heart, including EPDCs. To determine the role of SOX9 in epicardial development, an epicardial-specific Sox9 knockout mouse model was generated. Deleting Sox9 from the epicardial cell lineage impairs the ability of EPDCs to invade both the ventricular myocardium and the developing AV valves. After birth, the mitral valves of these mice become myxomatous with associated abnormalities in extracellular matrix organization. This phenotype is reminiscent of that seen in humans with myxomatous mitral valve disease (MVD). An RNA-seq analysis was conducted in an effort to identify genes associated with this myxomatous degeneration. From this experiment, Cd109 was identified as a gene associated with myxomatous valve pathogenesis in this model. Cd109 has never been described in the context of heart development or valve disease. This study highlights the importance of SOX9 in the regulation of epicardial cell invasion-emphasizing the importance of EPDCs in regulating AV valve development and homeostasis-and reports a novel expression profile of Cd109, a gene with previously unknown relevance in heart development.

RAAS-deficient organoids indicate delayed angiogenesis as a possible cause for autosomal recessive renal tubular dysgenesis

Autosomal Recessive Renal Tubular Dysgenesis (AR-RTD) is a fatal genetic disorder characterized by complete absence or severe depletion of proximal tubules (PT) in patients harboring pathogenic variants in genes involved in the Renin-Angiotensin-Aldosterone System. To uncover the pathomechanism of AR-RTD, differentiation of ACE-/- and AGTR1-/- induced pluripotent stem cells (iPSCs) and AR-RTD patient-derived iPSCs into kidney organoids is leveraged. Comprehensive marker analyses show that both mutant and control organoids generate indistinguishable PT in vitro under normoxic (21% O2) or hypoxic (2% O2) conditions. Fully differentiated (d24) AGTR1-/- and control organoids transplanted under the kidney capsule of immunodeficient mice engraft and mature well, as do renal vesicle stage (d14) control organoids. By contrast, d14 AGTR1-/- organoids fail to engraft due to insufficient pro-angiogenic VEGF-A expression. Notably, growth under hypoxic conditions induces VEGF-A expression and rescues engraftment of AGTR1-/- organoids at d14, as does ectopic expression of VEGF-A. We propose that PT dysgenesis in AR-RTD is primarily a non-autonomous consequence of delayed angiogenesis, starving PT at a critical time in their development.

Missense mutations in Myc Box I influence MYC cellular localization, mRNA partitioning and turnover to promote leukemogenesis

Somatic missense mutations in the phosphodegron domain of the MYC gene ( M YC Box I) are detected in the dominant clones of a subset of acute myeloid leukemia (AML) patients, but the mechanisms by which they contribute to AML are unknown. To unveil unique proprieties of MBI MYC mutant proteins, we systematically compared the cellular and molecular consequences of expressing similar oncogenic levels of wild type and MBI mutant MYC. We found that MBI MYC mutants can accelerate leukemia by driving unique transcriptional signatures in highly selected, myeloid progenitor subpopulations. Although these mutations increase MYC stability, they overall dampen MYC chromatin localization and lead to a cytoplasmic accumulation of the mutant proteins. This phenotype is coupled with increased translation of RNA binding proteins and nuclear export machinery, which results in altered RNA partitioning and accelerated decay of select transcripts encoding proapoptotic and proinflammatory genes. Heterozygous knockin mice harboring the germline MBI mutation Myc p.T73N exhibit cytoplasmic MYC localization, myeloid progenitors' expansion with similar transcriptional signatures to the overexpression model, and eventually develop hematological malignancies. This study uncovers that MBI MYC mutations alter MYC localization and disrupt mRNA subcellular distribution and turnover of select transcripts to accelerate tumor initiation and growth.

Complex urban atmosphere alters alveolar stem cells niche properties and drives lung fibrosis

There is growing evidence suggesting that urban pollution has adverse effects on lung health. However, how urban pollution affects alveolar mesenchymal and epithelial stem cell niches remains unknown. This study aimed to determine how complex representative urban atmospheres alter alveolar stem cell niche properties. Mice were placed in an innovative chamber realistically simulating the atmosphere of a megalopolis, or "clean air," for 7 days. Lungs were collected, and fibroblasts and epithelial cells (EpCAM+) were isolated. Proliferative capacities of fibroblasts were tested by population doubling levels (PDL), and microarray analyses were performed. Fibroblasts and EpCAM+ cells from exposed, nonexposed, or naive mice were cocultured in organoid assays to assess the stem cell properties. Collagen deposition (Sirius red), lipofibroblasts (ADRP, COL1A1), myofibroblasts (αSMA), alveolar type 2 cells (AT2, SFTPC+), and alveolar differentiation intermediate cell [ADI, keratin-8-positive (KRT8+)/claudin-4-positive (CLDN4+)] markers were quantified in the lungs. Fibroblasts obtained from mice exposed to urban atmosphere had lower PDL and survival and produced fewer and smaller organoids. Microarray analysis showed a decrease of adipogenesis and an increase of genes associated with fibrosis, suggesting a lipofibroblast to myofibroblast transition. Collagen deposition and myofibroblast number increased in the lungs of urban atmosphere-exposed mice. AT2 number was reduced and associated with an increase in ADI cells KRT8+/CLDN4+. Furthermore, EpCAM+ cells from exposed mice also produced fewer and smaller organoids. In conclusion, urban atmosphere alters alveolar mesenchymal stem cell niche properties by inducing a lipofibroblast to myofibroblast shift. It also results in alveolar epithelial dysfunction and a fibrotic-like phenotype.NEW & NOTEWORTHY Urban pollution is known to have major adverse effects on lung health. To assess the effect of pollution on alveolar regeneration, we exposed adult mice to a simulated high-pollution urban atmosphere, using an innovative CESAM simulation chamber (Multiphase Atmospheric Experimental Simulation Chamber, https://cesam.cnrs.fr/). We demonstrated that urban atmosphere alters alveolar mesenchymal stem cell niche properties by inducing a lipofibroblast to myofibroblast shift and induces alveolar epithelial dysfunction.

An extended transcription factor regulatory network controls hepatocyte identity

Cell identity is specified by a core transcriptional regulatory circuitry (CoRC), typically limited to a small set of interconnected cell-specific transcription factors (TFs). By mining global hepatic TF regulons, we reveal a more complex organization of the transcriptional regulatory network controlling hepatocyte identity. We show that tight functional interconnections controlling hepatocyte identity extend to non-cell-specific TFs beyond the CoRC, which we call hepatocyte identity (Hep-ID)CONNECT TFs. Besides controlling identity effector genes, Hep-IDCONNECT TFs also engage in reciprocal transcriptional regulation with TFs of the CoRC. In homeostatic basal conditions, this translates into Hep-IDCONNECT TFs being involved in fine tuning CoRC TF expression including their rhythmic expression patterns. Moreover, a role for Hep-IDCONNECT TFs in the control of hepatocyte identity is revealed in dedifferentiated hepatocytes where Hep-IDCONNECT TFs are able to reset CoRC TF expression. This is observed upon activation of NR1H3 or THRB in hepatocarcinoma or in hepatocytes subjected to inflammation-induced loss of identity. Our study establishes that hepatocyte identity is controlled by an extended array of TFs beyond the CoRC.

Associative gene networks reveal novel candidates important for ADHD and dyslexia comorbidity

BACKGROUND

Attention deficit hyperactivity disorder (ADHD) is commonly associated with developmental dyslexia (DD), which are both prevalent and complicated pediatric neurodevelopmental disorders that have a significant influence on children's learning and development. Clinically, the comorbidity incidence of DD and ADHD is between 25 and 48%. Children with DD and ADHD may have more severe cognitive deficiencies, a poorer level of schooling, and a higher risk of social and emotional management disorders. Furthermore, patients with this comorbidity are frequently treated for a single condition in clinical settings, and the therapeutic outcome is poor. The development of effective treatment approaches against these diseases is complicated by their comorbidity features. This is often a major problem in diagnosis and treatment. In this study, we developed bioinformatical methodology for the analysis of the comorbidity of these two diseases. As such, the search for candidate genes related to the comorbid conditions of ADHD and DD can help in elucidating the molecular mechanisms underlying the comorbid condition, and can also be useful for genotyping and identifying new drug targets.

RESULTS

Using the ANDSystem tool, the reconstruction and analysis of gene networks associated with ADHD and dyslexia was carried out. The gene network of ADHD included 599 genes/proteins and 148,978 interactions, while that of dyslexia included 167 genes/proteins and 27,083 interactions. When the ANDSystem and GeneCards data were combined, a total of 213 genes/proteins for ADHD and dyslexia were found. An approach for ranking genes implicated in the comorbid condition of the two diseases was proposed. The approach is based on ten criteria for ranking genes by their importance, including relevance scores of association between disease and genes, standard methods of gene prioritization, as well as original criteria that take into account the characteristics of an associative gene network and the presence of known polymorphisms in the analyzed genes. Among the top 20 genes with the highest priority DRD2, DRD4, CNTNAP2 and GRIN2B are mentioned in the literature as directly linked with the comorbidity of ADHD and dyslexia. According to the proposed approach, the genes OPRM1, CHRNA4 and SNCA had the highest priority in the development of comorbidity of these two diseases. Additionally, it was revealed that the most relevant genes are involved in biological processes related to signal transduction, positive regulation of transcription from RNA polymerase II promoters, chemical synaptic transmission, response to drugs, ion transmembrane transport, nervous system development, cell adhesion, and neuron migration.

CONCLUSIONS

The application of methods of reconstruction and analysis of gene networks is a powerful tool for studying the molecular mechanisms of comorbid conditions. The method put forth to rank genes by their importance for the comorbid condition of ADHD and dyslexia was employed to predict genes that play key roles in the development of the comorbid condition. The results can be utilized to plan experiments for the identification of novel candidate genes and search for novel pharmacological targets.

The A53T mutation in α-synuclein enhances pro-inflammatory activation in human microglia

Parkinson's disease (PD) is characterized by the aggregation of α-synuclein into Lewy bodies and Lewy neurites in the brain. Microglia-driven neuroinflammation may contribute to neuronal death in PD, however the exact role of microglia remains unclear and has been understudied. The A53T mutation in the gene coding for α-synuclein has been linked to early-onset PD, and exposure to A53T-mutant human α-synuclein increases the potential for inflammation of murine microglia. To date, its effect has not been studied in human microglia. Here, we used 2-dimensional cultures of human iPSC-derived microglia and transplantation of these cells into the mouse brain to assess the effects of the A53T mutation on human microglia. We found that A53T-mutant human microglia had an intrinsically increased propensity towards pro-inflammatory activation upon inflammatory stimulus. Additionally, A53T mutant microglia showed a strong decrease in catalase expression in non-inflammatory conditions, and increased oxidative stress. Our results indicate that A53T mutant human microglia display cell-autonomous phenotypes that may worsen neuronal damage in early-onset PD.

In Silico Analysis of the Effect of Hydrastis canadensis on Controlling Breast Cancer

Background and Objectives: Breast cancer is a significant type of cancer among women worldwide. Studies have reported the anti-carcinogenic activity of Hydrastis Canadensis (Goldenseal) in cancer cell lines. Hydrastis Canadensis could help eliminate toxic substances due to its anti-cancer, anti-inflammatory, and other properties. The design phase includes the identification of potential and effective molecules through modern computational techniques. Objective: This work aims to study Hydrastis Canadensis's effect in controlling hormone-independent breast cancer through in-silico analysis. Materials and Methods: The preliminary screening of reported phytochemicals includes biomolecular networking. Identifying functionally relevant phytochemicals and the respective target mutations/genes leads to selecting 3D proteins of the desired mutations being considered the target. Interaction studies have been conducted using docking. The kinetic and thermodynamic stability of complexes was studied through molecular dynamic simulation and MM-PBSA/GBSA analysis. Pharmacodynamic and pharmacokinetic features have been predicted. The mechanism-wise screening, functional enrichment, and interactional studies suggest that canadaline and Riboflavin effectively interact with the target proteins. Results: Hydrastis Canadensis has been identified as the effective formulation containing all these constituents. The phytoconstituents; Riboflavin and Canadensis showed good interaction with the targets of hormone-independent breast cancer. The complexes were found to be kinetically and thermodynamically stable. Conclusions: Hydrastis Canadensis has been identified as effective in controlling 'hormone-independent or basal-like breast cancer' followed by 'hormone-dependent breast cancer: Luminal A' and Luminal B.

Relationships between gene expression and behavior in mice in response to systemic modulation of the O-GlcNAcylation pathway

Enhancing protein O-GlcNAcylation by pharmacological inhibition of the enzyme O-GlcNAcase (OGA), which removes the O-GlcNAc modification from proteins, has been explored in mouse models of amyloid-beta and tau pathology. However, the O-GlcNAcylation-dependent link between gene expression and neurological behavior remains to be explored. Using chronic administration of Thiamet G (TG, an OGA inhibitor) in vivo, we used a protocol designed to relate behavior with the transcriptome and selected biochemical parameters from the cortex of individual animals. TG-treated mice showed improved working memory as measured using a Y-maze test. RNA sequencing analysis revealed 151 top differentially expressed genes with a Log2fold change >0.33 and adjusted p-value <0.05. Top TG-dependent upregulated genes were related to learning, cognition and behavior, while top downregulated genes were related to IL-17 signaling, inflammatory response and chemotaxis. Additional pathway analysis uncovered 3 pathways, involving gene expression including 14 cytochrome c oxidase subunits/regulatory components, chaperones or assembly factors, and 5 mTOR (mechanistic target of rapamycin) signaling factors. Multivariate Kendall correlation analyses of behavioral tests and the top TG-dependent differentially expressed genes revealed 91 statistically significant correlations in saline-treated mice and 70 statistically significant correlations in TG-treated mice. These analyses provide a network regulation landscape that is important in relating the transcriptome to behavior and the potential impact of the O-GlcNAC pathway.

Wide identification of chemical constituents in fermented licorice and explore its efficacy of anti-neurodegeneration by combining quasi-targeted metabolomics and in-depth bioinformatics

Licorice (Gan-Cao in Chinese) is one of the most famous herbal medicines around the world. The fermentation of probiotics and herbs can change the chemical constituents and significantly improve the efficacy. However, it is still unknown whether licorice fermented with probiotics would produce beneficial therapeutic effects. This study aimed to comprehensively analyze the chemical constituents in fermented licorice via quasi-targeted metabolomics, predict the potential efficacy of fermentation products via diverse bioinformatic methods, and further verify the efficacy of fermentation products through in vitro and in vivo experiments. As a result, 1,435 compounds were identified totally. Among them, 424 natural medicinal products were classified with potentially important bioactivities, including 11 anthocyanins, 10 chalcones and dihydrochalcones, 25 flavanones, 45 flavones and flavonols, 117 flavonoids, 34 isoflavonoids, 21 phenols and its derivatives, 20 phenylpropanoids and polyketides, 96 terpenoids and 25 coumarins and derivatives. Interestingly, bioinformatic prediction showed that the targets of some important compounds were related to neurodegeneration, oxidoreductase activity and response to stress. In vitro and in vivo tests further verified that fermented licorice had excellent effects of DPPH clearance, anti-oxidation, anti-neurodegeneration, and anti-stress. Thus, this study would provide a reference method for related research and the development of fermented licorice-related products.

Resources and tools for rare disease variant interpretation

Collectively, rare genetic disorders affect a substantial portion of the world's population. In most cases, those affected face difficulties in receiving a clinical diagnosis and genetic characterization. The understanding of the molecular mechanisms of these diseases and the development of therapeutic treatments for patients are also challenging. However, the application of recent advancements in genome sequencing/analysis technologies and computer-aided tools for predicting phenotype-genotype associations can bring significant benefits to this field. In this review, we highlight the most relevant online resources and computational tools for genome interpretation that can enhance the diagnosis, clinical management, and development of treatments for rare disorders. Our focus is on resources for interpreting single nucleotide variants. Additionally, we present use cases for interpreting genetic variants in clinical settings and review the limitations of these results and prediction tools. Finally, we have compiled a curated set of core resources and tools for analyzing rare disease genomes. Such resources and tools can be utilized to develop standardized protocols that will enhance the accuracy and effectiveness of rare disease diagnosis.

Identification of a heterogeneous and dynamic ciliome during embryonic development and cell differentiation

Primary cilia are nearly ubiquitous organelles that transduce molecular and mechanical signals. Although the basic structure of the cilium and the cadre of genes that contribute to ciliary formation and function (the ciliome) are believed to be evolutionarily conserved, the presentation of ciliopathies with narrow, tissue-specific phenotypes and distinct molecular readouts suggests that an unappreciated heterogeneity exists within this organelle. Here, we provide a searchable transcriptomic resource for a curated primary ciliome, detailing various subgroups of differentially expressed genes within the ciliome that display tissue and temporal specificity. Genes within the differentially expressed ciliome exhibited a lower level of functional constraint across species, suggesting organism and cell-specific function adaptation. The biological relevance of ciliary heterogeneity was functionally validated by using Cas9 gene-editing to disrupt ciliary genes that displayed dynamic gene expression profiles during osteogenic differentiation of multipotent neural crest cells. Collectively, this novel primary cilia-focused resource will allow researchers to explore longstanding questions related to how tissue and cell-type specific functions and ciliary heterogeneity may contribute to the range of phenotypes associated with ciliopathies.

Comparative transcriptomes reveal pro-survival and cytotoxic programs of mucosal-associated invariant T cells upon Bacillus Calmette-Guérin stimulation

Mucosal-associated invariant T (MAIT) cells are protective against tuberculous and non-tuberculous mycobacterial infections with poorly understood mechanisms. Despite an innate-like nature, MAIT cell responses remain heterogeneous in bacterial infections. To comprehensively characterize MAIT activation programs responding to different bacteria, we stimulated MAIT cells with E. coli to compare with Bacillus Calmette-Guérin (BCG), which remains the only licensed vaccine and a feasible tool for investigating anti-mycobacterial immunity in humans. Upon sequencing mRNA from the activated and inactivated CD8+ MAIT cells, results demonstrated the altered MAIT cell gene profiles by each bacterium with upregulated expression of activation markers, transcription factors, cytokines, and cytolytic mediators crucial in anti-mycobacterial responses. Compared with E. coli, BCG altered more MAIT cell genes to enhance cell survival and cytolysis. Flow cytometry analyses similarly displayed a more upregulated protein expression of B-cell lymphoma 2 and T-box transcription factor Eomesodermin in BCG compared to E.coli stimulations. Thus, the transcriptomic program and protein expression of MAIT cells together displayed enhanced pro-survival and cytotoxic programs in response to BCG stimulation, supporting BCG induces cell-mediated effector responses of MAIT cells to fight mycobacterial infections.

Identify potential driver genes for PAX-FOXO1 fusion-negative rhabdomyosarcoma through frequent gene co-expression network mining

Background

Rhabdomyosarcoma (RMS) is a soft tissue sarcoma usually originated from skeletal muscle. Currently, RMS classification based on PAX-FOXO1 fusion is widely adopted. However, compared to relatively clear understanding of the tumorigenesis in the fusion-positive RMS, little is known for that in fusion-negative RMS (FN-RMS).

Methods

We explored the molecular mechanisms and the driver genes of FN-RMS through frequent gene co-expression network mining (fGCN), differential copy number (CN) and differential expression analyses on multiple RMS transcriptomic datasets.

Results

We obtained 50 fGCN modules, among which five are differentially expressed between different fusion status. A closer look showed 23% of Module 2 genes are concentrated on several cytobands of chromosome 8. Upstream regulators such as MYC, YAP1, TWIST1 were identified for the fGCN modules. Using in a separate dataset we confirmed that, comparing to FP-RMS, 59 Module 2 genes show consistent CN amplification and mRNA overexpression, among which 28 are on the identified chr8 cytobands. Such CN amplification and nearby MYC (also resides on one of the above cytobands) and other upstream regulators (YAP1, TWIST1) may work together to drive FN-RMS tumorigenesis and progression. Up to 43.1% downstream targets of Yap1 and 45.8% of the targets of Myc are differentially expressed in FN-RMS vs. normal comparisons, which also confirmed the driving force of these regulators.

Discussion

We discovered that copy number amplification of specific cytobands on chr8 and the upstream regulators MYC, YAP1 and TWIST1 work together to affect the downstream gene co-expression and promote FN-RMS tumorigenesis and progression. Our findings provide new insights for FN-RMS tumorigenesis and offer promising targets for precision therapy. Experimental investigation about the functions of identified potential drivers in FN-RMS are in progress.

Characterization of Continuous Transcriptional Heterogeneity in High-Risk Blastemal-Type Wilms' Tumors Using Unsupervised Machine Learning

Wilms' tumors are pediatric malignancies that are thought to arise from faulty kidney development. They contain a wide range of poorly differentiated cell states resembling various distorted developmental stages of the fetal kidney, and as a result, differ between patients in a continuous manner that is not well understood. Here, we used three computational approaches to characterize this continuous heterogeneity in high-risk blastemal-type Wilms' tumors. Using Pareto task inference, we show that the tumors form a triangle-shaped continuum in latent space that is bounded by three tumor archetypes with "stromal", "blastemal", and "epithelial" characteristics, which resemble the un-induced mesenchyme, the cap mesenchyme, and early epithelial structures of the fetal kidney. By fitting a generative probabilistic "grade of membership" model, we show that each tumor can be represented as a unique mixture of three hidden "topics" with blastemal, stromal, and epithelial characteristics. Likewise, cellular deconvolution allows us to represent each tumor in the continuum as a unique combination of fetal kidney-like cell states. These results highlight the relationship between Wilms' tumors and kidney development, and we anticipate that they will pave the way for more quantitative strategies for tumor stratification and classification.

Ultra-Rare Variants Identify Biological Pathways and Candidate Genes in the Pathobiology of Non-Syndromic Cleft Palate Only

In recent decades, many efforts have been made to elucidate the genetic causes of non-syndromic cleft palate (nsCPO), a complex congenital disease caused by the interaction of several genetic and environmental factors. Since genome-wide association studies have evidenced a minor contribution of common polymorphisms in nsCPO inheritance, we used whole exome sequencing data to explore the role of ultra-rare variants in this study. In a cohort of 35 nsCPO cases and 38 controls, we performed a gene set enrichment analysis (GSEA) and a hypergeometric test for assessing significant overlap between genes implicated in nsCPO pathobiology and genes enriched in ultra-rare variants in our cohort. GSEA highlighted an enrichment of ultra-rare variants in genes principally belonging to cytoskeletal protein binding pathway (Probability Density Function corrected p-value = 1.57 × 10^-4); protein-containing complex binding pathway (p-value = 1.06 × 10^-2); cell adhesion molecule binding pathway (p-value = 1.24 × 10^-2); ECM-receptor interaction pathway (p-value = 1.69 × 10^-2); and in the Integrin signaling pathway (p-value = 1.28 × 10^-2). Two genes implicated in nsCPO pathobiology, namely COL2A1 and GLI3, ranked among the genes (n = 34) with nominal enrichment in the ultra-rare variant collapsing analysis (Fisher's exact test p-value < 0.05). These genes were also part of an independent list of genes highly relevant to nsCPO biology (n = 25). Significant overlap between the two sets of genes (hypergeometric test p-value = 5.86 × 10^-3) indicated that enriched genes are likely to be implicated in physiological palate development and/or the pathological processes of oral clefting. In conclusion, ultra-rare variants collectively impinge on biological pathways crucial to nsCPO pathobiology and point to candidate genes that may contribute to the individual risk of disease. Sequencing can be an effective approach to identify candidate genes and pathways for nsCPO.

The meningeal transcriptional response to traumatic brain injury and aging

Emerging evidence suggests that the meningeal compartment plays instrumental roles in various neurological disorders, however, we still lack fundamental knowledge about meningeal biology. Here, we utilized high-throughput RNA sequencing (RNA-seq) techniques to investigate the transcriptional response of the meninges to traumatic brain injury (TBI) and aging in the sub-acute and chronic time frames. Using single-cell RNA sequencing (scRNA-seq), we first explored how mild TBI affects the cellular and transcriptional landscape in the meninges in young mice at one-week post-injury. Then, using bulk RNA-seq, we assessed the differential long-term outcomes between young and aged mice following TBI. In our scRNA-seq studies, we highlight injury-related changes in differential gene expression seen in major meningeal cell populations including macrophages, fibroblasts, and adaptive immune cells. We found that TBI leads to an upregulation of type I interferon (IFN) signature genes in macrophages and a controlled upregulation of inflammatory-related genes in the fibroblast and adaptive immune cell populations. For reasons that remain poorly understood, even mild injuries in the elderly can lead to cognitive decline and devastating neuropathology. To better understand the differential outcomes between the young and the elderly following brain injury, we performed bulk RNA-seq on young and aged meninges 1.5 months after TBI. Notably, we found that aging alone induced upregulation of meningeal genes involved in antibody production by B cells and type I IFN signaling. Following injury, the meningeal transcriptome had largely returned to its pre-injury signature in young mice. In stark contrast, aged TBI mice still exhibited upregulation of immune-related genes and downregulation of genes involved in extracellular matrix remodeling. Overall, these findings illustrate the dynamic transcriptional response of the meninges to mild head trauma in youth and aging.

Fragile X Syndrome Patient-Derived Neurons Developing in the Mouse Brain Show FMR1-Dependent Phenotypes

BACKGROUND

Fragile X syndrome (FXS) is characterized by physical abnormalities, anxiety, intellectual disability, hyperactivity, autistic behaviors, and seizures. Abnormal neuronal development in FXS is poorly understood. Data on patients with FXS remain scarce, and FXS animal models have failed to yield successful therapies. In vitro models do not fully recapitulate the morphology and function of human neurons.

METHODS

To mimic human neuron development in vivo, we coinjected neural precursor cells derived from FXS patient-derived induced pluripotent stem cells and neural precursor cells derived from corrected isogenic control induced pluripotent stem cells into the brain of neonatal immune-deprived mice.

RESULTS

The transplanted cells populated the brain and a proportion differentiated into neurons and glial cells. Immunofluorescence and single and bulk RNA sequencing analyses showed accelerated maturation of FXS neurons after an initial delay. Additionally, we found increased percentages of Arc- and Egr-1-positive FXS neurons and wider dendritic protrusions of mature FXS striatal medium spiny neurons.

CONCLUSIONS

This transplantation approach provides new insights into the alterations of neuronal development in FXS by facilitating physiological development of cells in a 3-dimensional context.

Dual Oxidase 2 (DUOX2) as a Proteomic Biomarker for Predicting Treatment Response to Chemoradiation Therapy for Locally Advanced Rectal Cancer: Using High-Throughput Proteomic Analysis and Machine Learning Algorithm

High-throughput mass-spectrometry-based quantitative proteomic analysis was performed using formalin-fixed, paraffin-embedded (FFPE) biopsy samples obtained before treatment from 13 patients with locally advanced rectal cancer (LARC), who were treated with concurrent chemoradiation therapy (CCRT) followed by surgery. Patients were divided into complete responder (CR) and non-complete responder (nCR) groups. Immunohistochemical (IHC) staining of 79 independent FFPE tissue samples was performed to validate the predictive ability of proteomic biomarker candidates. A total of 3637 proteins were identified, and the expression of 498 proteins was confirmed at significantly different levels (differentially expressed proteins—DEPs) between two groups. In Gene Ontology enrichment analyses, DEPs enriched in biological processes in the CR group included proteins linked to cytoskeletal organization, immune response processes, and vesicle-associated protein transport processes, whereas DEPs in the nCR group were associated with biosynthesis, transcription, and translation processes. Dual oxidase 2 (DUOX2) was selected as the most predictive biomarker in machine learning algorithm analysis. Further IHC validation ultimately confirmed DUOX2 as a potential biomarker for predicting the response of nCR to CCRT. In conclusion, this study suggests that the treatment response to RT may be affected by the pre-treatment tumor microenvironment. DUOX2 is a potential biomarker for the early prediction of nCR after CCRT.

Aging exacerbates the brain inflammatory micro-environment contributing to α-synuclein pathology and functional deficits in a mouse model of DLB/PD

BACKGROUND

Although ɑ-synuclein (ɑ-syn) spreading in age-related neurodegenerative diseases such as Parkinson's disease (PD) and Dementia with Lewy bodies (DLB) has been extensively investigated, the role of aging in the manifestation of disease remains unclear.

METHODS

We explored the role of aging and inflammation in the pathogenesis of synucleinopathies in a mouse model of DLB/PD initiated by intrastriatal injection of ɑ-syn preformed fibrils (pff).

RESULTS

We found that aged mice showed more extensive accumulation of ɑ-syn in selected brain regions and behavioral deficits that were associated with greater infiltration of T cells and microgliosis. Microglial inflammatory gene expression induced by ɑ-syn-pff injection in young mice had hallmarks of aged microglia, indicating that enhanced age-associated pathologies may result from inflammatory synergy between aging and the effects of ɑ-syn aggregation. Based on the transcriptomics analysis projected from Ingenuity Pathway Analysis, we found a network that included colony stimulating factor 2 (CSF2), LPS related genes, TNFɑ and poly rl:rC-RNA as common regulators.

CONCLUSIONS

We propose that aging related inflammation (eg: CSF2) influences outcomes of pathological spreading of ɑ-syn and suggest that targeting neuro-immune responses might be important in developing treatments for DLB/PD.

IoMT-Based Mitochondrial and Multifactorial Genetic Inheritance Disorder Prediction Using Machine Learning

A genetic disorder is a serious disease that affects a large number of individuals around the world. There are various types of genetic illnesses, however, we focus on mitochondrial and multifactorial genetic disorders for prediction. Genetic illness is caused by a number of factors, including a defective maternal or paternal gene, excessive abortions, a lack of blood cells, and low white blood cell count. For premature or teenage life development, early detection of genetic diseases is crucial. Although it is difficult to forecast genetic disorders ahead of time, this prediction is very critical since a person's life progress depends on it. Machine learning algorithms are used to diagnose genetic disorders with high accuracy utilizing datasets collected and constructed from a large number of patient medical reports. A lot of studies have been conducted recently employing genome sequencing for illness detection, but fewer studies have been presented using patient medical history. The accuracy of existing studies that use a patient's history is restricted. The internet of medical things (IoMT) based proposed model for genetic disease prediction in this article uses two separate machine learning algorithms: support vector machine (SVM) and K-Nearest Neighbor (KNN). Experimental results show that SVM has outperformed the KNN and existing prediction methods in terms of accuracy. SVM achieved an accuracy of 94.99% and 86.6% for training and testing, respectively.

A SOX17-PDGFB signaling axis regulates aortic root development

Developmental etiologies causing complex congenital aortic root abnormalities are unknown. Here we show that deletion of Sox17 in aortic root endothelium in mice causes underdeveloped aortic root leading to a bicuspid aortic valve due to the absence of non-coronary leaflet and mispositioned left coronary ostium. The respective defects are associated with reduced proliferation of non-coronary leaflet mesenchyme and aortic root smooth muscle derived from the second heart field cardiomyocytes. Mechanistically, SOX17 occupies a Pdgfb transcriptional enhancer to promote its transcription and Sox17 deletion inhibits the endothelial Pdgfb transcription and PDGFB growth signaling to the non-coronary leaflet mesenchyme. Restoration of PDGFB in aortic root endothelium rescues the non-coronary leaflet and left coronary ostium defects in Sox17 nulls. These data support a SOX17-PDGFB axis underlying aortic root development that is critical for aortic valve and coronary ostium patterning, thereby informing a potential shared disease mechanism for concurrent anomalous aortic valve and coronary arteries.

Comparative Proteomic Analyses of Poorly Motile Swamp Buffalo Spermatozoa Reveal Low Energy Metabolism and Deficiencies in Motility-Related Proteins

The acquisition of mammalian sperm motility is a main indicator of epididymal sperm maturation and helps ensure fertilization. Poor sperm motility will prevent sperm cells from reaching the fertilization site, resulting in fertilization failure. To investigate the proteomic profiling of normal and poorly motile buffalo spermatozoa, a strategy applying liquid chromatography tandem mass spectrometry combined with tandem mass targeting was used. As a result, 145 differentially expressed proteins (DEPs) were identified in poorly motile spermatozoa (fold change > 1.5), including 52 upregulated and 93 downregulated proteins. The upregulated DEPs were mainly involved in morphogenesis and regulation of cell differentiation. The downregulated DEPs were involved with transport, oxidation-reduction, sperm motility, regulation of cAMP metabolism and regulation of DNA methylation. The mRNA and protein levels of PRM1 and AKAP3 were lower in poorly motile spermatozoa, while the expressions of SDC2, TEKT3 and IDH1 were not correlated with motility, indicating that their protein changes were affected by transcription or translation. Such changes in the expression of these proteins suggest that the formation of poorly motile buffalo spermatozoa reflects a low efficiency of energy metabolism, decreases in sperm protamine proteins, deficiencies in motility-related proteins, and variations in tail structural proteins. Such proteins could be biomarkers of poorly motile spermatozoa. These results illustrate some of the molecular mechanisms associated with poorly motile spermatozoa and provide clues for finding molecular markers of these pathways.

Association between resting-state functional brain connectivity and gene expression is altered in autism spectrum disorder

Gene expression covaries with brain activity as measured by resting state functional magnetic resonance imaging (MRI). However, it is unclear how genomic differences driven by disease state can affect this relationship. Here, we integrate from the ABIDE I and II imaging cohorts with datasets of gene expression in brains of neurotypical individuals and individuals with autism spectrum disorder (ASD) with regionally matched brain activity measurements from fMRI datasets. We identify genes linked with brain activity whose association is disrupted in ASD. We identified a subset of genes that showed a differential developmental trajectory in individuals with ASD compared with controls. These genes are enriched in voltage-gated ion channels and inhibitory neurons, pointing to excitation-inhibition imbalance in ASD. We further assessed differences at the regional level showing that the primary visual cortex is the most affected region in ASD. Our results link disrupted brain expression patterns of individuals with ASD to brain activity and show developmental, cell type, and regional enrichment of activity linked genes.

A Bayesian framework to integrate multi-level genome-scale data for Autism risk gene prioritization

Background

Autism spectrum disorder (ASD) is a group of complex neurodevelopment disorders with a strong genetic basis. Large scale sequencing studies have identified over one hundred ASD risk genes. Nevertheless, the vast majority of ASD risk genes remain to be discovered, as it is estimated that more than 1000 genes are likely to be involved in ASD risk. Prioritization of risk genes is an effective strategy to increase the power of identifying novel risk genes in genetics studies of ASD. As ASD risk genes are likely to exhibit distinct properties from multiple angles, we reason that integrating multiple levels of genomic data is a powerful approach to pinpoint genuine ASD risk genes.

Results

We present BNScore, a Bayesian model selection framework to probabilistically prioritize ASD risk genes through explicitly integrating evidence from sequencing-identified ASD genes, biological annotations, and gene functional network. We demonstrate the validity of our approach and its improved performance over existing methods by examining the resulting top candidate ASD risk genes against sets of high-confidence benchmark genes and large-scale ASD genome-wide association studies. We assess the tissue-, cell type- and development stage-specific expression properties of top prioritized genes, and find strong expression specificity in brain tissues, striatal medium spiny neurons, and fetal developmental stages.

Conclusions

In summary, we show that by integrating sequencing findings, functional annotation profiles, and gene-gene functional network, our proposed BNScore provides competitive performance compared to current state-of-the-art methods in prioritizing ASD genes. Our method offers a general and flexible strategy to risk gene prioritization that can potentially be applied to other complex traits as well.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04616-y.

Endosomal Sequestration of TLR4 Antibody Induces Myeloid-Derived Suppressor Cells and Reverses Acute Type 1 Diabetes

We previously showed that treating NOD mice with an agonistic monoclonal anti-TLR4/MD2 antibody (TLR4-Ab) reversed acute type 1 diabetes (T1D). Here, we show that TLR4-Ab reverses T1D by induction of myeloid-derived suppressor cells (MDSCs). Unbiased gene expression analysis after TLR4-Ab treatment demonstrated upregulation of genes associated with CD11b+Ly6G+ myeloid cells and downregulation of T-cell genes. Further RNA sequencing of purified, TLR4-Ab-treated CD11b+ cells showed significant upregulation of genes associated with bone marrow-derived CD11b+ cells and innate immune system genes. TLR4-Ab significantly increased percentages and numbers of CD11b+ cells. TLR4-Ab-induced CD11b+ cells, derived ex vivo from TLR4-Ab-treated mice, suppress T cells, and TLR4-Ab-conditioned bone marrow cells suppress acute T1D when transferred into acutely diabetic mice. Thus, the TLR4-Ab-induced CD11b+ cells, by the currently accepted definition, are MDSCs able to reverse T1D. To understand the TLR4-Ab mechanism, we compared TLR4-Ab with TLR4 agonist lipopolysaccharide (LPS), which cannot reverse T1D. TLR4-Ab remains sequestered at least 48 times longer than LPS within early endosomes, alters TLR4 signaling, and downregulates inflammatory genes and proteins, including nuclear factor-κB. TLR4-Ab in the endosome, therefore, induces a sustained, attenuated inflammatory response, providing an ideal "second signal" for the activation/maturation of MDSCs that can reverse acute T1D.

The genetic architecture of pediatric cardiomyopathy

To understand the genetic contribution to primary pediatric cardiomyopathy, we performed exome sequencing in a large cohort of 528 children with cardiomyopathy. Using clinical interpretation guidelines and targeting genes implicated in cardiomyopathy, we identified a genetic cause in 32% of affected individuals. Cardiomyopathy sub-phenotypes differed by ancestry, age at diagnosis, and family history. Infants < 1 year were less likely to have a molecular diagnosis (p < 0.001). Using a discovery set of 1,703 candidate genes and informatic tools, we identified rare and damaging variants in 56% of affected individuals. We see an excess burden of damaging variants in affected individuals as compared to two independent control sets, 1000 Genomes Project (p < 0.001) and SPARK parental controls (p < 1 × 10^-16). Cardiomyopathy variant burden remained enriched when stratified by ancestry, variant type, and sub-phenotype, emphasizing the importance of understanding the contribution of these factors to genetic architecture. Enrichment in this discovery candidate gene set suggests multigenic mechanisms underlie sub-phenotype-specific causes and presentations of cardiomyopathy. These results identify important information about the genetic architecture of pediatric cardiomyopathy and support recommendations for clinical genetic testing in children while illustrating differences in genetic architecture by age, ancestry, and sub-phenotype and providing rationale for larger studies to investigate multigenic contributions.

Effects of photodeoxygenation on cell biology using dibenzothiophene S-oxide derivatives as O(3P)-precursors

Photodeoxygenation of dibenzothiophene S-oxide and its derivatives have been used to generate atomic oxygen [O(³P)] to examine its effect on proteins, nucleic acids, and lipids. The unique reactivity and selectivity of O(³P) have shown distinct oxidation products and outcomes in biomolecules and cell-based studies. To understand the scope of its global impact on the cell, we treated MDA-MB-231 cells with 2,8-diacetoxymethyldibenzothiophene S-oxide and UV-A light to produce O(³P) without targeting a specific cell organelle. Cellular responses to O(³P)-release were analyzed using cell viability and cell cycle phase determination assays. Cell death was observed when cells were treated with higher concentrations of sulfoxides and UV-A light. However, significant differences in cell cycle phases due to UV-A irradiation of the sulfoxide were not observed. We further performed RNA-Seq analysis to study the underlying biological processes at play, and while UV-irradiation itself influenced gene expression, there were 9 upregulated and 8 downregulated genes that could be attributed to photodeoxygenation.

Quantitative phosphoproteomics reveals GSK3A substrate network is involved in the cryodamage of sperm motility

During sperm cryopreservation, the most significant phenotype of cryodamage is the decrease in sperm motility. Several proteomics studies have already been performed to search for key regulators at the protein level. However, sperm functions are known to be highly regulated by phosphorylation signaling. Here, we constructed a quantitative phosphoproteome to investigate the expression change of phosphorylated sites during sperm cryopreservation. A total of 3107 phosphorylated sites are identified and 848 of them are found to be significantly differentially expressed (DE). Bioinformatics analysis showed that the corresponding genes of these regulated sites are highly associated with sperm motility, providing a connection between the molecular basis and the phenotype of cryodamage. We then performed kinase enrichment analysis and successfully identified glycogen synthase kinase-3α (GSK3A) as the key kinase that may play an important role in the regulation of sperm motility. We further constructed a GSK3A centric network that could help us better understand the molecular mechanism of cryodamage in sperm motility. Finally, we also verified that GSK3A was abnormally activated during this process. The presented phosphoproteome and functional associations provide abundant research resources for us to learn the regulation of sperm functions, as well as to optimize the cryoprotectant for sperm cryopreservation.

Gene expression profiling utilizing extremely sensitive CDNA arrays and enrichment-based network study of major bone cancer genes

Background:

The gene interaction network is a set of genes interconnected by functional interactions among the genes. The gene interaction networks are studied to determine pathways and regulatory mechanisms in model organisms. In this research, the enrichment study of bone cancer-causing genes is undertaken to identify several hub genes associated to the development of bone cancer.

Materials and Methods:

Data on bone cancer is obtained from mutated gene samples; highly mutated genes are selected for the enrichment analysis. Due to certain interactions with each other the interaction network model for the hub genes is developed and simulations are produced to determine the levels of expression. For the array analyses, a total of 100 tumor specimens are collected. Cell cultures are prepared, RNA is extracted, cDNA arrays probes are generated, and the expressions analysis of Hub genes is determined.

Results:

Out of cDNA array findings, only 7 genes: CDKN2A, AKT1, NRAS, PIK3CA, RB1, BRAF, and TP53 are differentially expressed and shown as significant in the development of bone tumors, approximately 15 pathways have been identified, including pathways for non-small cell lung cancer, prostate cancer, pancreatic cancer, chronic myeloid leukemia, and glioma, consisting of all the identified 7 genes. After clinical validations of tumor samples, the IDH1 and TP53 gene revealed significant number of mutations similar to other genes. Specimens analysis showed that RB1, P53, and NRAS are amplified in brain tumor, while BRAF, CDKN2A, and AKT1 are amplified in sarcoma. Maximum deletion mutations of the PIK3CA gene are observed in leukemia. CDKN2A gene amplifications have been observed in virtually all tumor specimens.

Conclusion:

This study points to a recognizable evidence of novel superimposed pathways mechanisms strongly linked to cancer.

CD14 Expressing Precursors Give Rise to Highly Functional Conventional Dendritic Cells for Use as Dendritic Cell Vaccine

Simple Summary

Dendritic cells are attractive candidates for immunotherapy to prevent disease recurrence in cancer patients. Dendritic cells are a plastic population of antigen presenting cells and a variety of protocols have been described to generate dendritic cells from either monocytes or stem cells. To induce long lasting immunity, dendritic cells need to have a fully mature phenotype and activate naïve T-cells. Here, we describe a good manufacturer protocol to generate very potent conventional DC-like cells, derived from cord blood stem cells via a CD14+CD115+ precursor stage. They express high levels of CD1c and strongly activate both naïve as well as antigen-experienced T-cells. Implementation of this protocol in the clinic could advance the efficiency of dendritic cell based vaccines.

Abstract

Induction of long-lasting immunity by dendritic cells (DCs) makes them attractive candidates for anti-tumor vaccination. Although DC vaccinations are generally considered safe, clinical responses remain inconsistent in clinical trials. This initiated studies to identify subsets of DCs with superior capabilities to induce effective and memory anti-tumor responses. The use of primary DCs has been suggested to overcome the functional limitations of ex vivo monocyte-derived DCs (moDC). The ontogeny of primary DCs has recently been revised by the introduction of DC3, which phenotypically resembles conventional (c)DC2 as well as moDC. Previously, we developed a protocol to generate cDC2s from cord blood (CB)-derived stem cells via a CD115-expressing precursor. Here, we performed index sorting and single-cell RNA-sequencing to define the heterogeneity of in vitro developed DC precursors and identified CD14+CD115+ expressing cells that develop into CD1c++DCs and the remainder cells brought about CD123+DCs, as well as assessed their potency. The maturation status and T-cell activation potential were assessed using flow cytometry. CD123+DCs were specifically prone to take up antigens but only modestly activated T-cells. In contrast, CD1c++ are highly mature and specialized in both naïve as well as antigen-experienced T-cell activation. These findings show in vitro functional diversity between cord blood stem cell-derived CD123+DC and CD1c++DCs and may advance the efficiency of DC-based vaccines.

Identification of Novel Susceptible Genes of Gastric Cancer Based on Integrated Omics Data

Gastric cancer (GC) is one of the most common causes of cancer-related deaths in the world. This cancer has been regarded as a biological and genetically heterogeneous disease with a poorly understood carcinogenesis at the molecular level. Thousands of biomarkers and susceptible loci have been explored via experimental and computational methods, but their effects on disease outcome are still unknown. Genome-wide association studies (GWAS) have identified multiple susceptible loci for GC, but due to the linkage disequilibrium (LD), single-nucleotide polymorphisms (SNPs) may fall within the non-coding region and exert their biological function by modulating the gene expression level. In this study, we collected 1,091 cases and 410,350 controls from the GWAS catalog database. Integrating with gene expression level data obtained from stomach tissue, we conducted a machine learning-based method to predict GC-susceptible genes. As a result, we identified 787 novel susceptible genes related to GC, which will provide new insight into the genetic and biological basis for the mechanism and pathology of GC development.

Disease gene prediction with privileged information and heteroscedastic dropout

Motivation

Recently, machine learning models have achieved tremendous success in prioritizing candidate genes for genetic diseases. These models are able to accurately quantify the similarity among disease and genes based on the intuition that similar genes are more likely to be associated with similar diseases. However, the genetic features these methods rely on are often hard to collect due to high experimental cost and various other technical limitations. Existing solutions of this problem significantly increase the risk of overfitting and decrease the generalizability of the models.

Results

In this work, we propose a graph neural network (GNN) version of the Learning under Privileged Information paradigm to predict new disease gene associations. Unlike previous gene prioritization approaches, our model does not require the genetic features to be the same at training and test stages. If a genetic feature is hard to measure and therefore missing at the test stage, our model could still efficiently incorporate its information during the training process. To implement this, we develop a Heteroscedastic Gaussian Dropout algorithm, where the dropout probability of the GNN model is determined by another GNN model with a mirrored GNN architecture. To evaluate our method, we compared our method with four state-of-the-art methods on the Online Mendelian Inheritance in Man dataset to prioritize candidate disease genes. Extensive evaluations show that our model could improve the prediction accuracy when all the features are available compared to other methods. More importantly, our model could make very accurate predictions when >90% of the features are missing at the test stage.

Availability and implementation

Our method is realized with Python 3.7 and Pytorch 1.5.0 and method and data are freely available at: https://github.com/juanshu30/Disease-Gene-Prioritization-with-Privileged-Information-and-Heteroscedastic-Dropout.

ACE2 abrogates tumor resistance to VEGFR inhibitors suggesting angiotensin-(1-7) as a therapy for clear cell renal cell carcinoma

Angiotensin converting enzyme 2 (ACE2) is an enzyme that belongs to the renin-angiotensin system (RAS) and antagonizes the classical angiotensin (Ang) II/angiotensin II receptor type 1 (AT1) receptor pathway. Here, we report that higher ACE2 expression correlates with better overall survival in patients with clear cell renal cell carcinoma (ccRCC). Moreover, ACE2 has inhibitory effects on tumor proliferation in ccRCC in vitro and in preclinical animal models of ccRCC. We further show that Ang-(1-7), a heptapeptide generated by ACE2, is the likely mediator of this effect. Vascular endothelial growth factor receptor-tyrosine kinase inhibitor (VEGFR-TKI) treatment of ccRCC xenografts decreased ACE2 expression, and combination treatment with VEGFR-TKI and Ang-(1-7) generated additive suppression of tumor growth and improved survival outcomes. Last, the addition of Ang-(1-7) to programmed death-ligand 1 (PD-L1) pathway inhibitor and VEGFR-TKI showed further growth suppression in an immunocompetent RCC model. Together, these results suggest that targeting the ACE2/Ang-(1-7) axis is a promising therapeutic strategy against ccRCC.

Experimental and Computational Approaches to Improve Binding Affinity in Chemical Biology and Drug Discovery

Drug discovery is one of the most complicated processes and establishment of a single drug may require multidisciplinary attempts to design efficient and commercially viable drugs. The main purpose of drug design is to identify a chemical compound or inhibitor that can bind to an active site of a specific cavity on a target protein. The traditional drug design methods involved various experimental based approaches including random screening of chemicals found in nature or can be synthesized directly in chemical laboratories. Except for the long cycle design and time, high cost is also the major issue of concern. Modernized computer-based algorithm including structure-based drug design has accelerated the drug design and discovery process adequately. Surprisingly from the past decade remarkable progress has been made concerned with all area of drug design and discovery. CADD (Computer Aided Drug Designing) based tools shorten the conventional cycle size and also generate chemically more stable and worthy compounds and hence reduce the drug discovery cost. This special edition of editorial comprises the combination of seven research and review articles set emphasis especially on the computational approaches along with the experimental approaches using a chemical synthesizing for the binding affinity in chemical biology and discovery as a salient used in de-novo drug designing. This set of articles exfoliates the role that systems biology and the evaluation of ligand affinity in drug design and discovery for the future.

Activin-A Induces Early Differential Gene Expression Exclusively in Periodontal Ligament Fibroblasts from Fibrodysplasia Ossificans Progressiva Patients

Fibrodysplasia Ossificans Progressiva (FOP) is a rare genetic disease characterized by heterotopic ossification (HO). It is caused by mutations in the Activin receptor type 1 (ACVR1) gene, resulting in enhanced responsiveness to ligands, specifically to Activin-A. Though it has been shown that capturing Activin-A protects against heterotopic ossification in animal models, the exact underlying mechanisms at the gene expression level causing ACVR1 R206H-mediated ossifications and progression are thus far unknown. We investigated the early transcriptomic changes induced by Activin-A of healthy control and patient-derived periodontal ligament fibroblasts (PLF) isolated from extracted teeth by RNA sequencing analysis. To study early differences in response to Activin-A, periodontal ligament fibroblasts from six control teeth and from six FOP patient teeth were cultured for 24 h without and with 50 ng/mL Activin-A and analyzed with RNA sequencing. Pathway analysis on genes upregulated by Activin-A in FOP cells showed an association with pathways involved in, among others, Activin, TGFβ, and BMP signaling. Differential gene expression induced by Activin-A was exclusively seen in the FOP cells. Median centered supervised gene expression analysis showed distinct clusters of up- and downregulated genes in the FOP cultures after stimulation with Activin-A. The upregulated genes with high fold changes like SHOC2, TTC1, PAPSS2, DOCK7, and LOX are all associated with bone metabolism. Our open-ended approach to investigating the early effect of Activin-A on gene expression in control and FOP PLF shows that the molecule exclusively induces differential gene expression in FOP cells and not in control cells.

Alopecia-mental retardation syndrome: Molecular genetics of a rare neuro-dermal disorder

Alopecia-mental retardation syndrome (APMR) is a rare autosomal recessive neuro-dermal disorder. It is characterized by heterogeneous phenotypic features, that is, absence of hair on the scalp, eyelashes, and eyebrows and mild to severe intellectual disability. So far, approximately 14 families (i.e., Iranian, Pakistani, and Swiss) with APMR have been reported in the scientific literature. Its precise prevalence is still unknown, but according to a predictive estimate, it prevails with the ratio of 1 in 1,000,000 persons worldwide. Until now, only four loci (two characterized and two uncharacterized) have been reported to be involved in APMR. The pathogenic variants in alpha-2-HS-glycoprotein [AHSG; APMR1 (MIM#203650)] and lanosterol synthase [LSS; APMR4 (MIM#618840)] are the characterized genetic factors associated with APMR. Among them, AHSG was reported in a consanguineous Iranian family and LSS gene in a Swiss origin family, while the remaining two uncharacterized loci, that is, APMR2 and APMR3, are reported in the Pakistani population. The current mini-report discusses the molecular genetics and mutational spectrum of APMR syndrome, its differential diagnosis from related disorders, and prediction of plausible candidate genes in two uncharacterized loci.

AUTS2 isoforms control neuronal differentiation

Mutations in AUTS2 are associated with autism, intellectual disability, and microcephaly. AUTS2 is expressed in the brain and interacts with polycomb proteins, yet it is still unclear how mutations in AUTS2 lead to neurodevelopmental phenotypes. Here we report that when neuronal differentiation is initiated, there is a shift in expression from a long isoform to a short AUTS2 isoform. Yeast two-hybrid screen identified the splicing factor SF3B1 as an interactor of both isoforms, whereas the polycomb group proteins, PCGF3 and PCGF5, were found to interact exclusively with the long AUTS2 isoform. Reporter assays showed that the first exons of the long AUTS2 isoform function as a transcription repressor, but the part that consist of the short isoform acts as a transcriptional activator, both influenced by the cellular context. The expression levels of PCGF3 influenced the ability of the long AUTS2 isoform to activate or repress transcription. Mouse embryonic stem cells (mESCs) with heterozygote mutations in Auts2 had an increase in cell death during in vitro corticogenesis, which was significantly rescued by overexpressing the human AUTS2 transcripts. mESCs with a truncated AUTS2 protein (missing exons 12-20) showed premature neuronal differentiation, whereas cells overexpressing AUTS2, especially the long transcript, showed increase in expression of pluripotency markers and delayed differentiation. Taken together, our data suggest that the precise expression of AUTS2 isoforms is essential for regulating transcription and the timing of neuronal differentiation.

Combined the SMAC mimetic and BCL2 inhibitor sensitizes neoadjuvant chemotherapy by targeting necrosome complexes in tyrosine aminoacyl-tRNA synthase-positive breast cancer

Background

Chemotherapy is the standard treatment for breast cancer; however, the response to chemotherapy is disappointingly low. Here, we investigated the alternative therapeutic efficacy of novel combination treatment with necroptosis-inducing small molecules to overcome chemotherapeutic resistance in tyrosine aminoacyl-tRNA synthetase (YARS)-positive breast cancer.

Methods

Pre-chemotherapeutic needle biopsy of 143 invasive ductal carcinomas undergoing the same chemotherapeutic regimen was subjected to proteomic analysis. Four different machine learning algorithms were employed to determine signature protein combinations. Immunoreactive markers were selected using three common candidate proteins from the machine-learning algorithms and verified by immunohistochemistry using 123 cases of independent needle biopsy FFPE samples. The regulation of chemotherapeutic response and necroptotic cell death was assessed using lentiviral YARS overexpression and depletion 3D spheroid formation assay, viability assays, LDH release assay, flow cytometry analysis, and transmission electron microscopy. The ROS-induced metabolic dysregulation and phosphorylation of necrosome complex by YARS were assessed using oxygen consumption rate analysis, flow cytometry analysis, and 3D cell viability assay. The therapeutic roles of SMAC mimetics (LCL161) and a pan-BCL2 inhibitor (ABT-263) were determined by 3D cell viability assay and flow cytometry analysis. Additional biologic process and protein-protein interaction pathway analysis were performed using Gene Ontology annotation and Cytoscape databases.

Results

YARS was selected as a potential biomarker by proteomics-based machine-learning algorithms and was exclusively associated with good response to chemotherapy by subsequent immunohistochemical validation. In 3D spheroid models of breast cancer cell lines, YARS overexpression significantly improved chemotherapy response via phosphorylation of the necrosome complex. YARS-induced necroptosis sequentially mediated mitochondrial dysfunction through the overproduction of ROS in breast cancer cell lines. Combination treatment with necroptosis-inducing small molecules, including a SMAC mimetic (LCL161) and a pan-BCL2 inhibitor (ABT-263), showed therapeutic efficacy in YARS-overexpressing breast cancer cells.

Conclusions

Our results indicate that, before chemotherapy, an initial screening of YARS protein expression should be performed, and YARS-positive breast cancer patients might consider the combined treatment with LCL161 and ABT-263; this could be a novel stepwise clinical approach to apply new targeted therapy in breast cancer patients in the future.

Identification of candidate mediators of chemoresponse in breast cancer through therapy-driven selection of somatic variants

Purpose

More than a third of primary breast cancer patients are treated with cytotoxic chemotherapy, typically without guidance from predictive markers. Increased use of neoadjuvant chemotherapy provides opportunities for identification of molecules associated with treatment response, by comparing matched tumour samples before and after therapy. Our hypothesis was that somatic variants of increased prevalence after therapy promote resistance, while variants with reduced prevalence cause sensitivity.

Methods

We performed systematic analyses of matched pairs of cancer exomes from primary oestrogen receptor-positive/HER2-negative breast cancers (n = 6) treated with neoadjuvant epirubicin/cyclophosphamide. We identified candidate genes as mediators of chemotherapy response by consistent subclonal changes in somatic variant prevalence through therapy, predicted variant impact on gene function, and enrichment of specific functional pathways. Influence of candidate genes on breast cancer outcome was tested using publicly available breast cancer expression data (n = 1903).

Results

We identified 14 genes as the strongest candidate mediators of chemoresponse: TCHH, MUC17, ARAP2, FLG2, ABL1, CENPF, COL6A3, DMBT1, ITGA7, PLXNA1, S100PBP, SYNE1, ZFHX4, and CACNA1C. Genes contained somatic variants showing prevalence changes in up to 4 patients, with up to 3 being predicted as damaging. Genes coding for extra-cellular matrix components or related signalling pathways were significantly over-represented among variants showing prevalence changes. Expression of 5 genes (TCHH, ABL1, CENPF, S100PBP, and ZFHX4) was significantly associated with patient survival.

Conclusions

Genomic analysis of paired pre- and post-therapy samples resulting from neoadjuvant therapy provides a powerful method for identification of mediators of response. Genes we identified should be assessed as predictive markers or targets in chemo-sensitization.

Electronic supplementary material

The online version of this article (10.1007/s10549-020-05836-7) contains supplementary material, which is available to authorized users.

AP-1 activity induced by co-stimulation is required for chromatin opening during T cell activation

Activation of T cells is dependent on the organized and timely opening and closing of chromatin. Herein, we identify AP-1 as the transcription factor that directs most of this remodeling. Chromatin accessibility profiling showed quick opening of closed chromatin in naive T cells within 5 h of activation. These newly opened regions were strongly enriched for the AP-1 motif, and indeed, ChIP-seq demonstrated AP-1 binding at >70% of them. Broad inhibition of AP-1 activity prevented chromatin opening at AP-1 sites and reduced the expression of nearby genes. Similarly, induction of anergy in the absence of co-stimulation during activation was associated with reduced induction of AP-1 and a failure of proper chromatin remodeling. The translational relevance of these findings was highlighted by the substantial overlap of AP-1-dependent elements with risk loci for multiple immune diseases, including multiple sclerosis, inflammatory bowel disease, and allergic disease. Our findings define AP-1 as the key link between T cell activation and chromatin remodeling.

Genome-wide Enrichment of De Novo Coding Mutations in Orofacial Cleft Trios

Although de novo mutations (DNMs) are known to increase an individual's risk of congenital defects, DNMs have not been fully explored regarding orofacial clefts (OFCs), one of the most common human birth defects. Therefore, whole-genome sequencing of 756 child-parent trios of European, Colombian, and Taiwanese ancestry was performed to determine the contributions of coding DNMs to an individual's OFC risk. Overall, we identified a significant excess of loss-of-function DNMs in genes highly expressed in craniofacial tissues, as well as genes associated with known autosomal dominant OFC syndromes. This analysis also revealed roles for zinc-finger homeobox domain and SOX2-interacting genes in OFC etiology.

Machine learning-based approaches for disease gene prediction

Disease gene prediction is an essential issue in biomedical research. In the early days, annotation-based approaches were proposed for this problem. With the development of high-throughput technologies, interaction data between genes/proteins have grown quickly and covered almost genome and proteome; thus, network-based methods for the problem become prominent. In parallel, machine learning techniques, which formulate the problem as a classification, have also been proposed. Here, we firstly show a roadmap of the machine learning-based methods for the disease gene prediction. In the beginning, the problem was usually approached using a binary classification, where positive and negative training sample sets are comprised of disease genes and non-disease genes, respectively. The disease genes are ones known to be associated with diseases; meanwhile, non-disease genes were randomly selected from those not yet known to be associated with diseases. However, the later may contain unknown disease genes. To overcome this uncertainty of defining the non-disease genes, more realistic approaches have been proposed for the problem, such as unary and semi-supervised classification. Recently, more advanced methods, including ensemble learning, matrix factorization and deep learning, have been proposed for the problem. Secondly, 12 representative machine learning-based methods for the disease gene prediction were examined and compared in terms of prediction performance and running time. Finally, their advantages, disadvantages, interpretability and trust were also analyzed and discussed.