Exploring susceptibility and therapeutic targets for kidney stones through proteome-wide Mendelian randomization
Given the high recurrence rate of kidney stones, surgical lithotripsy and stone removal are not the ultimate treatments for kidney stones. There's an urgent need to explore the genetic mechanisms behind the susceptibility to kidney stones and to identify potential targets for prevention, to reduce the renal damage caused by recurrent stone formation. In this study, we screened 4548 circulating proteins using proteome-wide Mendelian Randomization (MR) to find proteins with a causal relationship to kidney stone risk. Additionally, proteome-wide association study (PWAS) and colocalization analysis were used to validate and prioritize candidate proteins. Moreover, downstream analyses including single-cell analysis, enrichment analysis, protein-protein interaction (PPI), and druggability analysis were conducted on the proteins causally related to kidney stones, to further explore the genetic mechanisms of susceptibility and the potential of proteins as drug targets. Ultimately, 22 target proteins associated with the risk of kidney stones were identified. Six plasma proteins (COLGALT1, CLMP, LECT1, ITIH1, CDHR3, CPLX2) were negatively correlated with kidney stone risk, while the genetic overexpression of 16 target proteins (GJA1, STOM, IRF9, F9, TMPRSS11D, ADH1B, SPINK13, CRYBB2, TNS2, DOCK9, OXSM, MST1, IL2, LMAN2, ITIH3, KLRF1) increased the risk of kidney stones. Based on the PWAS and colocalization analysis results, the 22 target proteins were classified into 3 tiers: IL2, CPLX2, and LMAN2 as tier 1 proteins with the most compelling evidence, MST1, ITIH1, and ITIH3 as tier 2 proteins, and the rest as tier 3 proteins. Enrichment analysis and PPI showed that target proteins mainly affect the occurrence of kidney stones through leukocyte activation and cell junction assembly. Druggability analysis suggested that IL2, MST1, and ITIH1 have potential as drug targets, and potential drugs were evaluated through molecular docking. In summary, this study employed multiple analytical methods to screen plasma proteins related to susceptibility to kidney stones, providing new insights into the genetic mechanisms of kidney stones and potential targets for treatment and prevention.
4-Phenylbutyric acid mitigates ER stress-induced neurodegeneration in the spinal cords of a GM2 gangliosidosis mouse model
Sandhoff disease (SD), a fatal and rare lysosomal storage disorder (LSD), is caused by a deficiency of the enzyme β-hexosaminidase B and leads to severe accumulation of GM2 gangliosides in lysosomes, primarily within the central nervous system (CNS). This accumulation results in severe neurological impairment, lower motor neuron disease, and death. Currently, there are no effective therapies available for SD. Here, we explored the role of endoplasmic reticulum (ER) stress in the spinal cord during disease progression in an established mouse model of SD and revealed the beneficial outcome of off-label treatment with the FDA-approved drug, 4-phenylbutyric acid (4-PBA). We analyzed the expression and localization of ER stress and cellular apoptosis markers, which revealed significant upregulation of these factors within motor neurons. Additionally, we observed a > 50% reduction in neuronal numbers throughout all spinal cord regions. Our studies also tested the impact of the chemical chaperone 4-PBA on ER stress in mice, and following administration, we observed significant improvements in motor neuromuscular function and life span throughout disease progression. 4-PBA treatment significantly reduced apoptosis in spinal cord neurons and increased the number of choline acetyltransferase (ChAT)-positive neurons, with little effect on astrogliosis or sensory interneurons. Overall, this study provides strong evidence for the role of chronic ER stress in the pathophysiology of SD and highlights 4-PBA as a promising therapeutic treatment for SD and potentially other related LSDs.
Analysis of genotype effects and inter-individual variability in iPSC-derived trisomy 21 neural progenitor cells
Trisomy of human chromosome 21 (T21) gives rise to Down syndrome (DS), the most frequent live-born autosomal aneuploidy. T21 triggers genome-wide transcriptomic alterations that result in multiple atypical phenotypes with highly variable penetrance and expressivity in individuals with DS. Many of these phenotypes, including atypical neurodevelopment, emerge prenatally. To enable in vitro analyses of the cellular and molecular mechanisms leading to the neurological alterations associated with T21, we created and characterized a panel of genomically diverse T21 and euploid induced pluripotent stem cells (iPSCs). We subsequently differentiated these iPSCs to generate a panel of neural progenitor cells (NPCs). Alongside characterizing genotype effects from T21, we found that T21 NPCs showed inter-individual variability in growth rates, oxidative stress, senescence characteristics, and gene and protein expression. Pathway enrichment analyses of T21 NPCs identified vesicular transport, DNA repair, and cellular response to stress pathways. These results demonstrate T21-associated variability at the cellular level and suggest that cell lines from individuals with DS should not solely be analyzed as a homogenous population. Examining large cohorts of genetically diverse samples may more fully reveal the effects of aneuploidy on transcriptomic and phenotypic characteristics in T21 cell types. A panel of genomically diverse T21 and euploid induced pluripotent stem cells (iPSCs) were created and subsequently differentiated into neural progenitor cells (NPCs). T21 NPCs showed reduced growth, increased oxidative stress, and inter-individual variability in gene and protein expression. This inter-individual variability suggests that studies with large cohorts of genetically diverse T21 samples may more fully reveal the effects of aneuploidy.
Loss of paired immunoglobin-like type 2 receptor B gene associated with age-related macular degeneration impairs photoreceptor function in mouse retina
Genome-wide association studies have uncovered mostly non-coding variants at over 60 genetic loci linked to susceptibility for age-related macular degeneration (AMD). To ascertain the causal gene at the PILRB/PILRA locus, we used a CRISPR strategy to produce germline deletions in the mouse paired immunoglobin-like type 2 receptor (Pilr) genes that encode highly related activating (PILRB) and inhibitory (PILRA) receptors. We show that a combined loss of Pilrb1 and Pilrb2, but not Pilra, leads to an early but relatively stationary defect as the electroretinography (ERG) amplitudes of Pilrb1/2-/- mice exhibit a marked reduction as early as postnatal day 15 and do not show additional significant decrease at 3 and 12-months. No alterations are evident in Müller glia, microglia, bipolar, amacrine and horizontal cells based on immunohistochemistry using cell-type specific markers. PILRB immunostaining is specifically detected at the proximal part of photoreceptor outer segment. Reduced expression of select calcium-regulated phototransduction and synapse-associated proteins, including GCAP1 and 2, PDE6b, AIPL1, PSD95, and CTBP1 indicates dysregulation of calcium homeostasis as a possible mechanism of retinal phenotype in Pilrb1/2-/- mice. Our studies suggest a novel function of PILRB in retinal photoreceptors and an association of PILRB, but not PILRA, with AMD pathogenesis.
Developing and validating a comprehensive polygenic risk score to enhance keratoconus risk prediction
This study aimed to develop and validate a comprehensive polygenic risk score (PRS) for keratoconus, enhancing the predictive accuracy for identifying individuals at increased risk, which is crucial for preventing keratoconus-associated visual impairment such as post-Laser-assisted in situ keratomileusis (LASIK) ectasia.
De novo missense variants in the PP2A regulatory subunit PPP2R2B in a neurodevelopmental syndrome: potential links to mitochondrial dynamics and spinocerebellar ataxias
The heterotrimeric protein phosphatase 2A (PP2A) complex catalyzes about half of Ser/Thr dephosphorylations in eukaryotic cells. A CAG repeat expansion in the neuron-specific protein PP2A regulatory subunit PPP2R2B gene causes spinocerebellar ataxia type 12 (SCA12). We established five monoallelic missense variants in PPP2R2B (four confirmed as de novo) as a cause of intellectual disability with developmental delay (R149P, T246K, N310K, E37K, I427T). In addition to moderate to severe intellectual disability and developmental delay, affected individuals presented with seizures, microcephaly, aggression, hypotonia, as well as broad-based or stiff gait. We used biochemical and cellular assays, including a novel luciferase complementation assay to interrogate PP2A holoenzyme assembly and activity, as well as deregulated mitochondrial dynamics as possible pathogenic mechanisms. Cell-based assays documented impaired ability of PPP2R2B missense variants to incorporate into the PP2A holoenzyme, localize to mitochondria, induce fission of neuronal mitochondria, and dephosphorylate the mitochondrial fission enzyme dynamin-related protein 1. AlphaMissense-based pathogenicity prediction suggested that an additional seven unreported missense variants may be pathogenic. In conclusion, our studies identify loss-of-function at the PPP2R2B locus as the basis for syndromic intellectual disability with developmental delay. They also extend PPP2R2B-related pathologies from neurodegenerative (SCA12) to neurodevelopmental disorders and suggests that altered mitochondrial dynamics may contribute to mechanisms.
Stem cell models of TAFAZZIN deficiency reveal novel tissue-specific pathologies in Barth syndrome
Barth syndrome (BTHS) is a rare mitochondrial disease caused by pathogenic variants in the gene TAFAZZIN, which leads to abnormal cardiolipin (CL) metabolism on the inner mitochondrial membrane. Although TAFAZZIN is ubiquitously expressed, BTHS involves a complex combination of tissue specific phenotypes including cardiomyopathy, neutropenia, skeletal myopathy, and growth delays, with a relatively minimal neurological burden. To understand both the developmental and functional effects of TAZ-deficiency in different tissues, we generated isogenic TAZ knockout (TAZ-KO) and WT cardiomyocytes (CMs) and neural progenitor cells (NPCs) from CRISPR-edited induced pluripotent stem cells (iPSCs). In TAZ-KO CMs we discovered evidence of dysregulated mitophagy including dysmorphic mitochondria and mitochondrial cristae, differential expression of key autophagy-associated genes, and an inability of TAZ-deficient CMs to properly initiate stress-induced mitophagy. In TAZ-deficient NPCs we identified novel phenotypes including a reduction in CIV abundance and CIV activity in the CIII2&CIV2 intermediate complex. Interestingly, while CL acyl chain manipulation was unable to alter mitophagy defects in TAZ-KO CMs, we found that linoleic acid or oleic acid supplementation was able to partially restore CIV abundance in TAZ-deficient NPCs. Taken together, our results have implications for understanding the tissue-specific pathology of BTHS and potential for tissue-specific therapeutic targeting. Moreover, our results highlight an emerging role for mitophagy in the cardiac pathophysiology of BTHS and reveal a potential neuron-specific bioenergetic phenotype.
Alternative polyadenylation shapes the molecular and clinical features of lung adenocarcinoma
Alternative polyadenylation (APA) is a major mechanism of post-transcriptional regulation that affects mRNA stability, localization and translation efficiency. Previous pan-cancer studies have revealed that APA is frequently disrupted in cancer and is associated with patient outcomes. Yet, little is known about cancer type-specific APA alterations. Here, we integrated RNA-sequencing data from a Korean cohort (GEO: GSE40419) and The Cancer Genome Atlas (TCGA) to comprehensively analyze APA alterations in lung adenocarcinomas (LUADs). Comparing expression levels of core genes involved in polyadenylation, we find that overall, the set of 28 of 31 genes are upregulated, with CSTF2 particularly upregulated. We observed broad and recurrent APA changes in LUAD growth-promoting genes. In addition, we find enrichment of APA events in genes associated with known LUAD pathways and an increased heterogeneity in polyadenylation (polyA) site usage of proliferation-associated genes. Upon further investigation, we report smoking-specific APA changes are also highly relevant to LUAD development. Overall, our in-depth analysis reveals APA as an important driver for the molecular and clinical features of lung adenocarcinoma.
Correction to: Expanded CAG/CTG repeats resist gene silencing mediated by targeted epigenome editing
Uplift of genetic diagnosis of rare respiratory disease using airway epithelium transcriptome analysis
Rare genetic respiratory disease has an incidence rate of more than 1:2500 live births in Northern Europe and carries significant disease burden. Early diagnosis improves outcomes, but many individuals remain without a confident genetic diagnosis. Improved and expanded molecular testing methods are required to improve genetic diagnosis rates and thereby improve clinical outcomes. Using primary ciliary dyskinesia (PCD) as an exemplar rare genetic respiratory disease, we developed a standardized method to identify pathogenic variants using whole transcriptome RNA-sequencing (RNA-seq) of nasal epithelial cells cultured at air-liquid interface (ALI). The method was optimized using cells from healthy volunteers, and people with rhino-pulmonary disease but no diagnostic indication of PCD. We validated the method using nasal epithelial cells from PCD patients with known genetic cause. We then assessed the ability of RNA-seq to identify pathogenic variants and the disease mechanism in PCD likely patients but in whom DNA genetic testing was inconclusive. The majority of 49 targeted PCD genes were optimally identified in RNA-seq data from nasal epithelial cells grown for 21 days at ALI culture. Four PCD-likely patients without a previous genetic diagnosis received a confirmed genetic diagnosis from the findings of the RNA-seq data. We demonstrate the clinical potential of RNA-seq of nasal epithelial cells to identify variants in individuals with genetically unsolved PCD. This uplifted genetic diagnosis should improve genetic counselling, enables family cascade screening, opens the door to potential personalised treatment and care approaches. This methodology could be implemented in other rare lung diseases such as cystic fibrosis.
The regulation of miRNAs using curcumin and other polyphenols during the prevention and treatment of Alzheimer's disease
Alzheimer's disease (AD), a prevalent neurodegenerative disorder, predominantly affects individuals over the age of 65 and poses significant challenges in terms of effective management and treatment. The disease's pathogenesis involves complex molecular pathways including misfolded proteins accumulation, neuroinflammation, and synaptic dysfunction. Recent insights have highlighted the role of microRNAs (miRNAs) as critical regulators within these pathways, where they influence gene expression and contribute to the pathophysiological landscape of AD. Notably, emerging research has demonstrated that polyphenols, including curcumin, might modulate miRNA activity, thus offering a novel approach to mitigate AD symptoms and progression. This review explores the potential mechanisms through which polyphenols regulate miRNA expression and activity, specifically focusing on autophagy enhancement and inflammation reduction in the context of AD. We provide a detailed examination of key studies linking miRNA dysregulation to AD pathogenesis and discuss how polyphenols might correct these aberrations. The findings presented here underscore the therapeutic potential of polyphenols in AD treatment via miRNA modulation, pointing to new directions in disease management strategies and highlighting the need for targeted research into miRNA-based interventions.
ZBTB24 is a conserved multifaceted transcription factor at genes and centromeres that governs the DNA methylation state and expression of satellite repeats
Since its discovery as a causative gene of the Immunodeficiency with Centromeric instability and Facial anomalies syndrome, ZBTB24 has emerged as a key player in DNA methylation, immunity and development. By extensively analyzing ZBTB24 genomic functions in ICF-relevant mouse and human cellular models, we document here its multiple facets as a transcription factor, with key roles in immune response-related genes expression and also in early embryonic development. Using a constitutive Zbtb24 ICF-like mutant and an auxin-inducible degron system in mouse embryonic stem cells, we showed that ZBTB24 is recruited to centromeric satellite DNA where it is required to establish and maintain the correct DNA methylation patterns through the recruitment of DNMT3B. The ability of ZBTB24 to occupy centromeric satellite DNA is conserved in human cells. Together, our results unveiled an essential and underappreciated role for ZBTB24 at mouse and human centromeric satellite repeat arrays by controlling their DNA methylation and transcription status.
Brain multi-omic Mendelian randomisation to identify novel drug targets for gliomagenesis
Genetic variants associated with molecular traits that are also associated with liability to glioma can provide causal evidence for the identification and prioritisation of drug targets.
SMN depletion impairs skeletal muscle formation and maturation in a mouse model of SMA
Spinal muscular atrophy (SMA) is characterized by low levels of the ubiquitously expressed Survival Motor Neuron (SMN) protein, leading to progressive muscle weakness and atrophy. Skeletal muscle satellite cells play a crucial role in muscle fiber maintenance, repair, and remodelling. While the effects of SMN depletion in muscle are well documented, its precise role in satellite cell function remains largely unclear. Using the Smn2B/- mouse model, we investigated SMN-depleted satellite cell biology through single fiber culture studies. Myofibers from Smn2B/- mice were smaller in size, shorter in length, had reduced myonuclear domain size, and reduced sub-synaptic myonuclear clusters-all suggesting impaired muscle function and integrity. These changes were accompanied by a reduction in the number of myonuclei in myofibers from Smn2B/- mice across all disease stages examined. Although the number of satellite cells in myofibers was significantly reduced, those remaining retained their capacity for myogenic activation and proliferation. These findings support the idea that a dysregulated myogenic process could be occurring as early in muscle stem cells during muscle formation and maturation in SMA. Targeting those pathways could offer additional options for combinatorial therapies for SMA.
Ligand distances as key predictors of pathogenicity and function in NMDA receptors
Genetic variants in the genes GRIN1, GRIN2A, GRIN2B, and GRIN2D, which encode subunits of the N-methyl-D-aspartate receptor (NMDAR), have been associated with severe and heterogeneous neurologic and neurodevelopmental disorders, including early onset epilepsy, developmental and epileptic encephalopathy, intellectual disability, and autism spectrum disorders. Missense variants in these genes can result in gain or loss of the NMDAR function, requiring opposite therapeutic treatments. Computational methods that predict pathogenicity and molecular functional effects of missense variants are therefore crucial for therapeutic applications. We assembled 223 missense variants from patients, 631 control variants from the general population, and 160 missense variants characterized by electrophysiological readouts that show whether they can enhance or reduce the function of the receptor. This includes new functional data from 33 variants reported here, for the first time. By mapping these variants onto the NMDAR protein structures, we found that pathogenic/benign variants and variants that increase/decrease the channel function were distributed unevenly on the protein structure, with spatial proximity to ligands bound to the agonist and antagonist binding sites being a key predictive feature for both variant pathogenicity and molecular functional consequences. Leveraging distances from ligands, we developed two machine-learning based predictors for NMDA variants: a pathogenicity predictor which outperforms currently available predictors and the first molecular function (increase/decrease) predictor. Our findings can have direct application to patient care by improving diagnostic yield for genetic neurodevelopmental disorders and by guiding personalized treatment informed by the knowledge of the molecular disease mechanism.
Characterisation of LGP2 complex multitranscript system in humans: role in the innate immune response and evolution from non-human primates
Retinoic acid inducible gene I (RIG-I)-like receptors (RLRs), including RIG-I, MDA5 and LGP2, recognize viral RNA to mount an antiviral interferon (IFN) response RLRs share three different protein domains: C-terminal domain, DExD/H box RNA helicase domain, and an N-terminal domain with two tandem repeats (CARDs). LGP2 lacks tandem CARD and is not able to induce an IFN response. However, LGP2 positively enhances MDA5 and negatively regulates RIG-I signaling. In this study, we determined the LGP2 alternative transcripts in humans to further comprehend the mechanism of its regulation, their evolutionary origin, and the isoforms functionallity. The results showed new eight alternative transcripts in the samples tested. The presence of these transcripts demonstrated that the main mechanisms for the regulation of LGP2 expression are both by insertion of introns and by the loss of exons. The phylogenetic analysis of the comparison between sequences from exon 1 to exon 3 of humans and those previously described in non-human primates showed three well-differentiated groups (lineages) originating from gorillas, suggesting that the transspecies evolution has been maintained for 10 million years. The corresponding protein models (isoforms) were also established, obtaining four isoforms: one complete and three others lacking the C-terminal domain or this domain and the partial or total He2 Helicase domain, which would compromise the functionality of LGP2. In conclusion, this is the first study that elucidate the large genomic organization and complex transcriptional regulation of human LGP2, its pattern of sequence generation, and a mode of evolutionary inheritance across species.
Genetic analysis implicates ERAP1 and HLA as risk factors for severe Puumala virus infection
Puumala virus (PUUV) infections can cause severe illnesses such as Hemorrhagic Fever with Renal Syndrome in humans. However, human genetic risk factors contributing to disease severity are still poorly understood. Our goal was to elucidate genetic factors contributing to PUUV infections and understand the biological mechanisms underlying individual vulnerability to PUUV infections. Leveraging data from the FinnGen study, we conducted a genome-wide association study on severe Hemorrhagic Fever with Renal Syndrome caused by PUUV with 2227 cases. We identified associations at the Human Leukocyte Antigen (HLA) locus and ERAP1 with severe PUUV infection. HLA molecules are canonical mediators for immune recognition and response. ERAP1 facilitates immune system recognition and activation by cleaving viral proteins into smaller peptides which are presented to the immune system via HLA class I molecules. Notably, we identified that the lead variant (rs26653, OR = 0.84, P = 2.9 × 10-8) in the ERAP1 gene was a missense variant changing amino acid arginine to proline. From the HLA region, we showed independent and significant associations with both HLA class I and II genes. Furthermore, we showed independent associations with four HLA alleles with severe PUUV infection using conditional HLA fine mapping. The strongest association was found with the HLA-C*07:01 allele (OR = 1.54, P = 4.0 × 10-24) followed by signals at HLA-B*13:02, HLA-DRB1*01:01, and HLA-DRB1*11:01 alleles (P < 5 × 10-8). Our findings suggest an association of viral peptide processing with ERAP1 and antigen presentation through HLA alleles that may contribute to the development of severe PUUV disease.
The genetic architecture of hip shape and its role in the development of hip osteoarthritis and fracture
Hip shape is thought to be an important causal risk factor for hip osteoarthritis and fracture. We aimed to identify genetic determinants of hip shape and use these to assess causal relationships with hip osteoarthritis.
A novel GFAP frameshift variant identified in a family with optico-retinal dysplasia and vision impairment
Gain-of-function variants in GFAP leads to protein aggregation and is the cause of the severe neurodegenerative disorder Alexander Disease (AxD), while loss of GFAP function has been considered benign. Here, we investigated a six-generation family, where multiple individuals presented with gliosis of the optic nerve head and visual impairment. Whole genome sequencing (WGS) revealed a frameshift variant in GFAP (c.928dup, p.(Met310Asnfs*113)) segregating with disease. Analysis of human embryonic tissues revealed strong expression of GFAP in retinal neural progenitors. A zebrafish model verified that c.928dup does not result in extensive GFAP protein aggregation and zebrafish gfap loss-of-function mutants showed vision impairment and retinal dysplasia, characterized by a significant loss of Müller glia cells and photoreceptor cells. Our findings show how different mutational mechanisms can cause diverging phenotypes and reveal a novel function of GFAP in vertebrate eye development.
Correction to: Reduced levels of MRE11 cause disease phenotypes distinct from ataxia telangiectasia-like disorder
Increased intrinsic membrane excitability is associated with olivary hypertrophy in spinocerebellar ataxia type 1
One of the characteristic regions of brainstem degeneration across multiple spinocerebellar ataxias (SCAs) is the inferior olive (IO), a medullary nucleus that plays a key role in motor learning. The vulnerability of IO neurons remains a poorly-understood area of SCA pathology. In this work, we address this by evaluating IO disease in SCA1, a prototypic inherited olivopontocerebellar atrophy, using the genetically-precise SCA1 knock-in (SCA1-KI) mouse. We find that these mice exhibit olivary hypertrophy, a phenotype reminiscent of a degenerative disorder known as hypertrophic olivary degeneration (HOD). Similar to early stages of HOD, SCA1-KI IO neurons display early dendritic lengthening and later somatic expansion without frank cell loss. Though HOD is known to be caused by brainstem lesions that disrupt IO inhibitory innervation, we observe no loss of inhibitory terminals in the SCA1-KI IO. Additionally, we find that a separate mouse model of SCA1 in which mutant ATXN1 is expressed solely in cerebellar Purkinje cells shows no evidence of olivary hypertrophy. Patch-clamp recordings from brainstem slices indicate that SCA1-KI IO neurons are hyperexcitable, generating spike trains in response to membrane depolarization. Transcriptome analysis further reveals reduced medullary expression of ion channels responsible for IO neuron spike afterhyperpolarization (AHP)-a result that appears to have a functional consequence, as SCA1-KI IO neuron spikes exhibit a diminished AHP. These findings suggest that expression of mutant ATXN1 in IO neurons results in an HOD-like olivary hypertrophy, in association with increased intrinsic membrane excitability and ion channel transcriptional dysregulation.