STIGMA: Single-cell tissue-specific gene prioritization using machine learning

Clinical exome and genome sequencing have revolutionized the understanding of human disease genetics. Yet many genes remain functionally uncharacterized, complicating the establishment of causal disease links for genetic variants. While several scoring methods have been devised to prioritize these candidate genes, these methods fall short of capturing the expression heterogeneity across cell subpopulations within tissues. Here, we introduce single-cell tissue-specific gene prioritization using machine learning (STIGMA), an approach that leverages single-cell RNA-seq (scRNA-seq) data to prioritize candidate genes associated with rare congenital diseases. STIGMA prioritizes genes by learning the temporal dynamics of gene expression across cell types during healthy organogenesis. To assess the efficacy of our framework, we applied STIGMA to mouse limb and human fetal heart scRNA-seq datasets. In a cohort of individuals with congenital limb malformation, STIGMA prioritized 469 variants in 345 genes, with UBA2 as a notable example. For congenital heart defects, we detected 34 genes harboring nonsynonymous de novo variants (nsDNVs) in two or more individuals from a set of 7,958 individuals, including the ortholog of Prdm1, which is associated with hypoplastic left ventricle and hypoplastic aortic arch. Overall, our findings demonstrate that STIGMA effectively prioritizes tissue-specific candidate genes by utilizing single-cell transcriptome data. The ability to capture the heterogeneity of gene expression across cell populations makes STIGMA a powerful tool for the discovery of disease-associated genes and facilitates the identification of causal variants underlying human genetic disorders.

Evaluation of long acting GLP1R/GCGR agonist in a DIO and biopsy-confirmed mouse model of NASH suggest a beneficial role of GLP-1/glucagon agonism in NASH patients

OBJECTIVE

The metabolic benefits of GLP-1 receptor (GLP-1R) agonists on glycemic and weight control are well established as therapy for type 2 diabetes and obesity. Glucagon's ability to increase energy expenditure is well described, and the combination of these mechanisms-of-actions has the potential to further lower hepatic steatosis in metabolic disorders and could therefore be attractive for the treatment for non-alcoholic steatohepatitis (NASH). Here, we have investigated the effects of a dual GLP-1/glucagon receptor agonist NN1177 on hepatic steatosis, fibrosis, and inflammation in a preclinical mouse model of NASH. Having observed strong effects on body weight loss in a pilot study with NN1177, we hypothesized that direct engagement of the hepatic glucagon receptor (GCGR) would result in a superior effect on steatosis and other liver related parameters as compared to the GLP-1R agonist semaglutide at equal body weight.

METHODS

Male C57Bl/6 mice were fed a diet high in trans-fat, fructose, and cholesterol (Diet-Induced Obese (DIO)-NASH) for 36 weeks. Following randomization based on the degree of fibrosis at baseline, mice were treated once daily with subcutaneous administration of a vehicle or three different doses of NN1177 or semaglutide for 8 weeks. Hepatic steatosis, inflammation and fibrosis were assessed by immunohistochemistry and morphometric analyses. Plasma levels of lipids and liver enzymes were determined, and hepatic gene expression was analyzed by RNA sequencing.

RESULTS

NN1177 dose-dependently reduced body weight up to 22% compared to vehicle treatment. Plasma levels of ALT, a measure of liver injury, were reduced in all treatment groups with body weight loss. The dual agonist reduced hepatic steatosis to a greater extent than semaglutide at equal body weight loss, as demonstrated by three independent methods. Both the co-agonist and semaglutide significantly decreased histological markers of inflammation such as CD11b and Galectin-3, in addition to markers of hepatic stellate activation (αSMA) and fibrosis (Collagen I). Interestingly, the maximal beneficial effects on above mentioned clinically relevant endpoints of NN1177 treatment on hepatic health appear to be achieved with the middle dose tested. Administering the highest dose resulted in a further reduction of liver fat and accompanied by a massive induction in genes involved in oxidative phosphorylation and resulted in exaggerated body weight loss and a downregulation of a module of co-expressed genes involved in steroid hormone biology, bile secretion, and retinol and linoleic acid metabolism that are also downregulated due to NASH itself.

CONCLUSIONS

These results indicate that, in a setting of overnutrition, the liver health benefits of activating the fasting-related metabolic pathways controlled by the glucagon receptor displays a bell-shaped curve. This observation is of interest to the scientific community, due to the high number of ongoing clinical trials attempting to leverage the positive effects of glucagon biology to improve metabolic health.

Enhancer hijacking at the ARHGAP36 locus is associated with connective tissue to bone transformation

Heterotopic ossification is a disorder caused by abnormal mineralization of soft tissues in which signaling pathways such as BMP, TGFβ and WNT are known key players in driving ectopic bone formation. Identifying novel genes and pathways related to the mineralization process are important steps for future gene therapy in bone disorders. In this study, we detect an inter-chromosomal insertional duplication in a female proband disrupting a topologically associating domain and causing an ultra-rare progressive form of heterotopic ossification. This structural variant lead to enhancer hijacking and misexpression of ARHGAP36 in fibroblasts, validated here by orthogonal in vitro studies. In addition, ARHGAP36 overexpression inhibits TGFβ, and activates hedgehog signaling and genes/proteins related to extracellular matrix production. Our work on the genetic cause of this heterotopic ossification case has revealed that ARHGAP36 plays a role in bone formation and metabolism, outlining first details of this gene contributing to bone-formation and -disease.

Biallelic variants in ADAMTS15 cause a novel form of distal arthrogryposis

PURPOSE

We aimed to identify the underlying genetic cause for a novel form of distal arthrogryposis.

METHODS

Rare variant family-based genomics, exome sequencing, and disease-specific panel sequencing were used to detect ADAMTS15 variants in affected individuals. Adamts15 expression was analyzed at the single-cell level during murine embryogenesis. Expression patterns were characterized using in situ hybridization and RNAscope.

RESULTS

We identified homozygous rare variant alleles of ADAMTS15 in 5 affected individuals from 4 unrelated consanguineous families presenting with congenital flexion contractures of the interphalangeal joints and hypoplastic or absent palmar creases. Radiographic investigations showed physiological interphalangeal joint morphology. Additional features included knee, Achilles tendon, and toe contractures, spinal stiffness, scoliosis, and orthodontic abnormalities. Analysis of mouse whole-embryo single-cell sequencing data revealed a tightly regulated Adamts15 expression in the limb mesenchyme between embryonic stages E11.5 and E15.0. A perimuscular and peritendinous expression was evident in in situ hybridization in the developing mouse limb. In accordance, RNAscope analysis detected a significant coexpression with Osr1, but not with markers for skeletal muscle or joint formation.

CONCLUSION

In aggregate, our findings provide evidence that rare biallelic recessive trait variants in ADAMTS15 cause a novel autosomal recessive connective tissue disorder, resulting in a distal arthrogryposis syndrome.

Repression and 3D-restructuring resolves regulatory conflicts in evolutionarily rearranged genomes

Regulatory landscapes drive complex developmental gene expression, but it remains unclear how their integrity is maintained when incorporating novel genes and functions during evolution. Here, we investigated how a placental mammal-specific gene, Zfp42, emerged in an ancient vertebrate topologically associated domain (TAD) without adopting or disrupting the conserved expression of its gene, Fat1. In ESCs, physical TAD partitioning separates Zfp42 and Fat1 with distinct local enhancers that drive their independent expression. This separation is driven by chromatin activity and not CTCF/cohesin. In contrast, in embryonic limbs, inactive Zfp42 shares Fat1's intact TAD without responding to active Fat1 enhancers. However, neither Fat1 enhancer-incompatibility nor nuclear envelope-attachment account for Zfp42's unresponsiveness. Rather, Zfp42's promoter is rendered inert to enhancers by context-dependent DNA methylation. Thus, diverse mechanisms enabled the integration of independent Zfp42 regulation in the Fat1 locus. Critically, such regulatory complexity appears common in evolution as, genome wide, most TADs contain multiple independently expressed genes.

Single-cell sequencing of human midbrain reveals glial activation and a Parkinson-specific neuronal state

Idiopathic Parkinson's disease is characterized by a progressive loss of dopaminergic neurons, but the exact disease aetiology remains largely unknown. To date, Parkinson's disease research has mainly focused on nigral dopaminergic neurons, although recent studies suggest disease-related changes also in non-neuronal cells and in midbrain regions beyond the substantia nigra. While there is some evidence for glial involvement in Parkinson's disease, the molecular mechanisms remain poorly understood. The aim of this study was to characterize the contribution of all cell types of the midbrain to Parkinson's disease pathology by single-nuclei RNA sequencing and to assess the cell type-specific risk for Parkinson's disease using the latest genome-wide association study. We profiled >41 000 single-nuclei transcriptomes of post-mortem midbrain from six idiopathic Parkinson's disease patients and five age-/sex-matched controls. To validate our findings in a spatial context, we utilized immunolabelling of the same tissues. Moreover, we analysed Parkinson's disease-associated risk enrichment in genes with cell type-specific expression patterns. We discovered a neuronal cell cluster characterized by CADPS2 overexpression and low TH levels, which was exclusively present in idiopathic Parkinson's disease midbrains. Validation analyses in laser-microdissected neurons suggest that this cluster represents dysfunctional dopaminergic neurons. With regard to glial cells, we observed an increase in nigral microglia in Parkinson's disease patients. Moreover, nigral idiopathic Parkinson's disease microglia were more amoeboid, indicating an activated state. We also discovered a reduction in idiopathic Parkinson's disease oligodendrocyte numbers with the remaining cells being characterized by a stress-induced upregulation of S100B. Parkinson's disease risk variants were associated with glia- and neuron-specific gene expression patterns in idiopathic Parkinson's disease cases. Furthermore, astrocytes and microglia presented idiopathic Parkinson's disease-specific cell proliferation and dysregulation of genes related to unfolded protein response and cytokine signalling. While reactive patient astrocytes showed CD44 overexpression, idiopathic Parkinson's disease microglia revealed a pro-inflammatory trajectory characterized by elevated levels of IL1B, GPNMB and HSP90AA1. Taken together, we generated the first single-nuclei RNA sequencing dataset from the idiopathic Parkinson's disease midbrain, which highlights a disease-specific neuronal cell cluster as well as 'pan-glial' activation as a central mechanism in the pathology of the movement disorder. This finding warrants further research into inflammatory signalling and immunomodulatory treatments in Parkinson's disease.

Parkin Deficiency Impairs Mitochondrial DNA Dynamics and Propagates Inflammation

BACKGROUND

Mutations in the E3 ubiquitin ligase parkin cause autosomal recessive Parkinson's disease (PD). Together with PTEN-induced kinase 1 (PINK1), parkin regulates the clearance of dysfunctional mitochondria. New mitochondria are generated through an interplay of nuclear- and mitochondrial-encoded proteins, and recent studies suggest that parkin influences this process at both levels. In addition, parkin was shown to prevent mitochondrial membrane permeability, impeding mitochondrial DNA (mtDNA) escape and subsequent neuroinflammation. However, parkin's regulatory roles independent of mitophagy are not well described in patient-derived neurons.

OBJECTIVES

We sought to investigate parkin's role in preventing neuronal mtDNA dyshomeostasis, release, and glial activation at the endogenous level.

METHODS

We generated induced pluripotent stem cell (iPSC)-derived midbrain neurons from PD patients with parkin (PRKN) mutations and healthy controls. Live-cell imaging, proteomic, mtDNA integrity, and gene expression analyses were employed to investigate mitochondrial biogenesis and genome maintenance. To assess neuroinflammation, we performed single-nuclei RNA sequencing in postmortem tissue and quantified interleukin expression in mtDNA/lipopolysaccharides (LPS)-treated iPSC-derived neuron-microglia co-cultures.

RESULTS

Neurons from patients with PRKN mutations revealed deficits in the mitochondrial biogenesis pathway, resulting in mtDNA dyshomeostasis. Moreover, the energy sensor sirtuin 1, which controls mitochondrial biogenesis and clearance, was downregulated in parkin-deficient cells. Linking mtDNA disintegration to neuroinflammation, in postmortem midbrain with PRKN mutations, we confirmed mtDNA dyshomeostasis and detected an upregulation of microglia overexpressing proinflammatory cytokines. Finally, parkin-deficient neuron-microglia co-cultures elicited an enhanced immune response when exposed to mtDNA/LPS.

CONCLUSIONS

Our findings suggest that parkin coregulates mitophagy, mitochondrial biogenesis, and mtDNA maintenance pathways, thereby protecting midbrain neurons from neuroinflammation and degeneration. © 2022 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Systems Immunology of Diabetes-Tuberculosis Comorbidity Reveals Signatures of Disease Complications

Comorbid diabetes mellitus (DM) increases tuberculosis (TB) risk and adverse outcomes but the pathological interactions between DM and TB remain incompletely understood. We performed an integrative analysis of whole blood gene expression and plasma analytes, comparing South Indian TB patients with and without DM to diabetic and non-diabetic controls without TB. Luminex assay of plasma cytokines and growth factors delineated a distinct biosignature in comorbid TBDM in this cohort. Transcriptional profiling revealed elements in common with published TB signatures from cohorts that excluded DM. Neutrophil count correlated with the molecular degree of perturbation, especially in TBDM patients. Body mass index and HDL cholesterol were negatively correlated with molecular degree of perturbation. Diabetic complication pathways including several pathways linked to epigenetic reprogramming were activated in TBDM above levels observed with DM alone. Our data provide a rationale for trials of host-directed therapies in TBDM, targeting neutrophilic inflammation and diabetic complication pathways to address the greater morbidity and mortality associated with this increasingly prevalent dual burden of communicable and non-communicable diseases.