17q21.31 sub-haplotypes underlying H1-associated risk for Parkinson's disease are associated with LRRC37A/2 expression in astrocytes

Induction of acrosome reaction by 4-Br-A23187 alters the glycoproteomic profile of boar spermatozoa

Protein glycosylation is a post-translational modification involved in wide range of biological processes. In mammalian spermatozoa this modification has been identified in numerous proteins, and membrane glycoproteins are involved in the fertilization process. The objective of the present study was to identify changes in protein glycosylation after acrosome reaction (AR) induction using the 4-Br-A23187 ionophore. Our results showed that treatment with 10 μM of 4-Br-A23187 for 20 min significantly increased the percentage of live acrosome-reacted spermatozoa compared to the control (69.8 ± 0.8 vs. 6.4 ± 0.5; mean % ± SEM, respectively). Also, we observed an increase in 32 kDa tyrosine-phosphorylated protein (p32) and a decrease in serine/threonine phosphorylation of the protein kinase A substrates (phospho-PKA-substrates) after ionophore treatment. Furthermore, changes in glycosylated proteins following AR induction were analyzed using different HRP-conjugated lectins (GNA, DSA, and SNA), revealing changes in mannose and sialic acid residues. Proteomic analysis of isolated proteins using GNA lectin revealed that 50 proteins exhibited significantly different abundance (q-value < 0.01). Subsequent analysis using Uniprot database identified 39 downregulated and 11 upregulated proteins in the presence of 4-Br-A23187. Notably, six of these proteins were classified as transmembrane proteins, namely LRRC37A/B like protein 1 C-terminal domain-containing protein, Membrane metalloendopeptidase like 1, VWFA domain-containing protein, Syndecan, Membrane spanning 4-domains A14 and Serine protease 54. This study shows a novel protocol to induce acrosome reaction in boar spermatozoa and identifies new transmembrane proteins containing mannose residues. Further work is needed to elucidate the role of these proteins in sperm-oocyte fusion.

Bioinformatics analysis of an immunotherapy responsiveness-related gene signature in predicting lung adenocarcinoma prognosis

Background

Immune therapy has become first-line treatment option for patients with lung cancer, but some patients respond poorly to immune therapy, especially among patients with lung adenocarcinoma (LUAD). Novel tools are needed to screen potential responders to immune therapy in LUAD patients, to better predict the prognosis and guide clinical decision-making. Although many efforts have been made to predict the responsiveness of LUAD patients, the results were limited. During the era of immunotherapy, this study attempts to construct a novel prognostic model for LUAD by utilizing differentially expressed genes (DEGs) among patients with differential immune therapy responses.

Methods

Transcriptome data of 598 patients with LUAD were downloaded from The Cancer Genome Atlas (TCGA) database, which included 539 tumor samples and 59 normal control samples, with a mean follow-up time of 29.69 months (63.1% of patients remained alive by the end of follow-up). Other data sources including three datasets from the Gene Expression Omnibus (GEO) database were analyzed, and the DEGs between immunotherapy responders and nonresponders were identified and screened. Univariate Cox regression analysis was applied with the TCGA cohort as the training set and GSE72094 cohort as the validation set, and least absolute shrinkage and selection operator (LASSO) Cox regression were applied in the prognostic-related genes which fulfilled the filter criteria to establish a prognostic formula, which was then tested with time-dependent receiver operating characteristic (ROC) analysis. Enriched pathways of the prognostic-related genes were analyzed with Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses, and tumor immune microenvironment (TIME), tumor mutational burden, and drug sensitivity tests were completed with appropriate packages in R (The R Foundation of Statistical Computing). Finally, a nomogram incorporating the prognostic formula was established.

Results

A total of 1,636 DEGs were identified, 1,163 prognostic-related DEGs were extracted, and 34 DEGs were selected and incorporated into the immunotherapy responsiveness-related risk score (IRRS) formula. The IRRS formula had good performance in predicting the overall prognoses in patients with LUAD and had excellent performance in prognosis prediction in all LUAD subgroups. Moreover, the IRRS formula could predict anticancer drug sensitivity and immunotherapy responsiveness in patients with LUAD. Mechanistically, immune microenvironments varied profoundly between the two IRRS groups; the most significantly varied pathway between the high-IRRS and low-IRRS groups was ribonucleoprotein complex biogenesis, which correlated closely with the TP53 and TTN mutation burdens. In addition, we established a nomogram incorporating the IRRS, age, sex, clinical stage, T-stage, N-stage, and M-stage as predictors that could predict the prognoses of 1-year, 3-year, and 5-year survival in patients with LUAD, with an area under curve (AUC) of 0.718, 0.702, and 0.68, respectively.

Conclusions

The model we established in the present study could predict the prognosis of LUAD patients, help to identify patients with good responses to anticancer drugs and immunotherapy, and serve as a valuable tool to guide clinical decision-making.

Unraveling the complex role of MAPT-containing H1 and H2 haplotypes in neurodegenerative diseases

A ~ 1 Mb inversion polymorphism exists within the 17q21.31 locus of the human genome as direct (H1) and inverted (H2) haplotype clades. This inversion region demonstrates high linkage disequilibrium, but the frequency of each haplotype differs across ancestries. While the H1 haplotype exists in all populations and shows a normal pattern of genetic variability and recombination, the H2 haplotype is enriched in European ancestry populations, is less frequent in African ancestry populations, and nearly absent in East Asian ancestry populations. H1 is a known risk factor for several neurodegenerative diseases, and has been associated with many other traits, suggesting its importance in cellular phenotypes of the brain and entire body. Conversely, H2 is protective for these diseases, but is associated with predisposition to recurrent microdeletion syndromes and neurodevelopmental disorders such as autism. Many single nucleotide variants and copy number variants define H1/H2 haplotypes and sub-haplotypes, but identifying the causal variant(s) for specific diseases and phenotypes is complex due to the extended linkage equilibrium. In this review, we assess the current knowledge of this inversion region regarding genomic structure, gene expression, cellular phenotypes, and disease association. We discuss recent discoveries and challenges, evaluate gaps in knowledge, and highlight the importance of understanding the effect of the 17q21.31 haplotypes to promote advances in precision medicine and drug discovery for several diseases.

Critical reasoning on the co-expression module QTL in the dorsolateral prefrontal cortex

Expression quantitative trait locus (eQTL) analysis is a popular method of gaining insight into the function of regulatory variation. While cis-eQTL resources have been instrumental in linking genome-wide association study variants to gene function, complex trait heritability may be additionally mediated by other forms of gene regulation. Toward this end, novel eQTL methods leverage gene co-expression (module-QTL) to investigate joint regulation of gene modules by single genetic variants. Here we broadly define a "module-QTL" as the association of a genetic variant with a summary measure of gene co-expression. This approach aims to reduce the multiple testing burden of a trans-eQTL search through the consolidation of gene-based testing and provide biological context to eQTLs shared between genes. In this article we provide an in-depth examination of the co-expression module eQTL (module-QTL) through literature review, theoretical investigation, and real-data application of the module-QTL to three large prefrontal cortex genotype-RNA sequencing datasets. We find module-QTLs in our study that are disease associated and reproducible are not additionally informative beyond cis- or trans-eQTLs for module genes. Through comparison to prior studies, we highlight promises and limitations of the module-QTL across study designs and provide recommendations for further investigation of the module-QTL framework.

Deciphering the role of a SINE-VNTR-Alu retrotransposon polymorphism as a biomarker of Parkinson's disease progression

SINE-VNTR-Alu (SVA) retrotransposons are transposable elements which represent a source of genetic variation. We previously demonstrated that the presence/absence of a human-specific SVA, termed SVA_67, correlated with the progression of Parkinson's disease (PD). In the present study, we demonstrate that SVA_67 acts as expression quantitative trait loci, thereby exhibiting a strong regulatory effect across the genome using whole genome and transcriptomic data from the Parkinson's progression markers initiative cohort. We further show that SVA_67 is polymorphic for its variable number tandem repeat domain which correlates with both regulatory properties in a luciferase reporter gene assay in vitro and differential expression of multiple genes in vivo. Additionally, this variation's utility as a biomarker is reflected in a correlation with a number of PD progression markers. These experiments highlight the plethora of transcriptomic and phenotypic changes associated with SVA_67 polymorphism which should be considered when investigating the missing heritability of neurodegenerative diseases.

Validation of Enhancer Regions in Primary Human Neural Progenitor Cells using Capture STARR-seq

Genome-wide association studies (GWAS) and expression analyses implicate noncoding regulatory regions as harboring risk factors for psychiatric disease, but functional characterization of these regions remains limited. We performed capture STARR-sequencing of over 78,000 candidate regions to identify active enhancers in primary human neural progenitor cells (phNPCs). We selected candidate regions by integrating data from NPCs, prefrontal cortex, developmental timepoints, and GWAS. Over 8,000 regions demonstrated enhancer activity in the phNPCs, and we linked these regions to over 2,200 predicted target genes. These genes are involved in neuronal and psychiatric disease-associated pathways, including dopaminergic synapse, axon guidance, and schizophrenia. We functionally validated a subset of these enhancers using mutation STARR-sequencing and CRISPR deletions, demonstrating the effects of genetic variation on enhancer activity and enhancer deletion on gene expression. Overall, we identified thousands of highly active enhancers and functionally validated a subset of these enhancers, improving our understanding of regulatory networks underlying brain function and disease.

Neuronal MAPT expression is mediated by long-range interactions with cis-regulatory elements

Tauopathies are a group of neurodegenerative diseases defined by abnormal aggregates of tau, a microtubule-associated protein encoded by MAPT. MAPT expression is near absent in neural progenitor cells (NPCs) and increases during differentiation. This temporally dynamic expression pattern suggests that MAPT expression could be controlled by transcription factors and cis-regulatory elements specific to differentiated cell types. Given the relevance of MAPT expression to neurodegeneration pathogenesis, identification of such elements is relevant to understanding disease risk and pathogenesis. Here, we performed chromatin conformation assays (HiC & Capture-C), single-nucleus multiomics (RNA-seq+ATAC-seq), bulk ATAC-seq, and ChIP-seq for H3K27ac and CTCF in NPCs and differentiated neurons to nominate candidate cis-regulatory elements (cCREs). We assayed these cCREs using luciferase assays and CRISPR interference (CRISPRi) experiments to measure their effects on MAPT expression. Finally, we integrated cCRE annotations into an analysis of genetic variation in neurodegeneration-affected individuals and control subjects. We identified both proximal and distal regulatory elements for MAPT and confirmed the regulatory function for several regions, including three regions centromeric to MAPT beyond the H1/H2 haplotype inversion breakpoint. We also found that rare and predicted damaging genetic variation in nominated CREs was nominally depleted in dementia-affected individuals relative to control subjects, consistent with the hypothesis that variants that disrupt MAPT enhancer activity, and thereby reduced MAPT expression, may be protective against neurodegenerative disease. Overall, this study provides compelling evidence for pursuing detailed knowledge of CREs for genes of interest to permit better understanding of disease risk.

Broken silence: 22,841 predicted deleterious synonymous variants identified in the human exome through computational analysis

Synonymous single nucleotide variants (sSNVs) do not alter the primary structure of a protein, thus it was previously accepted that they were neutral. Recently, several studies demonstrated their significance to a range of diseases. Still, variant prioritization strategies lack focus on sSNVs. Here, we identified 22,841 deleterious synonymous variants in 125,748 human exomes using two in silico predictors (SilVA and CADD). While 98.2% of synonymous variants are classified as neutral, 1.8% are predicted to be deleterious, yielding an average of 9.82 neutral and 0.18 deleterious sSNVs per exome. Further investigation of prediction features via Heterogeneous Ensemble Feature Selection revealed that impact on amino acid sequence and conservation carry the most weight for a deleterious prediction. Thirty nine detrimental sSNVs are not rare and are located on disease associated genes. Ten distinct putatively non-deleterious sSNVs are likely to be under positive selection in the North-Western European and East Asian populations. Taken together our analysis gives voice to the so-called silent mutations as we propose a robust framework for evaluating the deleteriousness of sSNVs in variant prioritization studies.

MAPT Locus in Parkinson's Disease Patients of Ashkenazi Origin: A Stratified Analysis

Introduction: MAPT locus is associated with Parkinson's disease (PD), which is located within a large inversion region of high linkage disequilibrium (LD). We aimed to determine whether the H2-haplotype protective effect and its effect size depends on the GBA1 or LRRK2 risk allele carrier status, and to further characterize genetic alterations that might contribute to its effect. Methods: LD analysis was performed using whole-genome sequencing data of 202 unrelated Ashkenazi Jewish (AJ) PDs. A haplotype-divergent variant was genotyped in a cohort of 1200 consecutively recruited AJ-PDs. The odd ratios were calculated using AJ-non-neuro cases from the gnomAD database as the controls in an un-stratified and a stratified manner according to the mutation carrier status, and the effect on the Age at Motor Symptom Onset (AMSO) was examined. Expression and splicing quantitative trait locus (eQTL and sQTL) analyses were carried out using brain tissues from a database. Results: The H2 haplotype exhibited significant association with PD protection, with a similar effect size in GBA1 carriers, LRRK2-G2019S carriers, and non-carriers (OR = 0.77, 0.69, and 0.82, respectively), and there was no effect on AMSO. The LD interval was narrowed to approximately 1.2 Mb. The H2 haplotype carried potential variants in candidate genes (MAPT and SPPL2C); structural deletions and segmental duplication (KANSL1); and variants affecting gene expression and intron excision ratio in brain tissues (LRRC37A/2). Conclusions: Our results demonstrate that H2 is associated with PD and its protective effect is not influenced by the GBA1/LRRK2 risk allele carrier status. This effect may be genetically complex, resulting from different levels of variations such as missense mutations in relevant genes, structural variations, epigenetic modifications, and RNA expression changes, which may operate independently or in synergy.

Interindividual variation in human cortical cell type abundance and expression

Single-cell transcriptomic studies have identified a conserved set of neocortical cell types from small postmortem cohorts. We extended these efforts by assessing cell type variation across 75 adult individuals undergoing epilepsy and tumor surgeries. Nearly all nuclei map to one of 125 robust cell types identified in the middle temporal gyrus. However, we found interindividual variance in abundances and gene expression signatures, particularly in deep-layer glutamatergic neurons and microglia. A minority of donor variance is explainable by age, sex, ancestry, disease state, and cell state. Genomic variation was associated with expression of 150 to 250 genes for most cell types. This characterization of cellular variation provides a baseline for cell typing in health and disease.

Shared Genetics and Comorbid Genes of Amyotrophic Lateral Sclerosis and Parkinson's Disease

BACKGROUND

Comorbidity exists between amyotrophic lateral sclerosis (ALS) and Parkinson's disease (PD), but the role of genetic factors is unclear.

OBJECTIVE

We aim to investigate genetic correlation, causal relationship, and comorbid genes between ALS and PD.

METHODS

Leveraging the largest genome-wide association study data (ALS: 27,205 cases, 110,881 controls; PDG: 33,674 cases, 449,056 controls), we used linkage disequilibrium score regression and Mendelian randomization analysis for genetic correlation and causal inference. We performed genome-wide cross-trait analysis via Multi-Trait Analysis of Genome-Wide Association Studies and Cross-Phenotype Association to identify specific single-nucleotide polymorphisms, followed by functional mapping and annotation. Integrating expression quantitative trait loci data from 13 brain regions, we conducted a transcriptome-wide association study via functional summary-based imputation and joint-tissue imputation to explore comorbid genes, followed by pathway enrichment analysis.

RESULTS

We found that PD positively correlates with ALS (rg  = 0.144, P = 0.026) and confers a causal effect (odds ratio = 1.09, 95% confidence interval: 1.03-1.15, P = 3.00 × 10-3 ). We identified nine single-nucleotide polymorphisms (eight new), associating with three risk loci (chromosomes 4, 10, and 17) and seven genes (TMEM175, MAPT, NSF, LRRC37A2, ARHGAP27, GAK, and FGFRL1). In transcriptome-wide association study analysis, we showed six previously unreported pleiotropic genes (KANSL1, ARL17B, EFNA1, WNT3, ERCC8, and ADAM15), and we found these candidate genes are mainly enriched in negative regulation of neuron projection development (GO:0010977).

CONCLUSIONS

Our work demonstrates shared genetic architecture between ALS and PD, reports new pleiotropic genes, and sheds light on the comorbid mechanism. © 2023 International Parkinson and Movement Disorder Society.

CRISPR deletion of a SINE-VNTR-Alu (SVA_67) retrotransposon demonstrates its ability to differentially modulate gene expression at the MAPT locus

Background

SINE-VNTR-Alu (SVA) retrotransposons are hominid-specific elements which have been shown to play important roles in processes such as chromatin structure remodelling and regulation of gene expression demonstrating that these repetitive elements exert regulatory functions. We have previously shown that the presence or absence of a specific SVA element, termed SVA_67, was associated with differential expression of several genes at the MAPT locus, a locus associated with Parkinson's Disease (PD) and frontotemporal dementia. However, we were not able to demonstrate that causation of differential gene expression was directed by the SVA due to lack of functional validation.

Methods

We performed CRISPR to delete SVA_67 in the HEK293 cell line. Quantification of target gene expression was performed using qPCR to assess the effects on expression in response to the deletion of SVA_67. Differences between CRISPR edit and control cell lines were analysed using two-tailed t-test with a minimum 95% confidence interval to determine statistical significance.

Results

In this study, we provide data highlighting the SVA-specific effect on differential gene expression. We demonstrate that the hemizygous deletion of the endogenous SVA_67 in CRISPR edited cell lines was associated with differential expression of several genes at the MAPT locus associated with neurodegenerative diseases including KANSL1, MAPT and LRRC37A.

Discussion

This data is consistent with our previous bioinformatic work of differential gene expression analysis using transcriptomic data from the Parkinson's Progression Markers Initiative (PPMI) cohort. As SVAs have regulatory influences on gene expression, and insertion polymorphisms contribute to interpersonal differences in expression patterns, these results highlight the potential contribution of these elements to complex diseases with potentially many genetic components, such as PD.

MAPT allele and haplotype frequencies in Nigerian Africans: Population distribution and association with Parkinson's disease risk and age at onset

INTRODUCTION

The association between MAPT and PD risk may be subject to ethnic variability even within populations of similar geographical origin. Data on MAPT haplotype frequencies, and its association with PD risk in black Africans are lacking. We aimed to determine the frequencies of MAPT haplotypes and their role as risk factors for PD and age at onset in Nigerians.

METHODS

The haplotype and genotype frequencies of MAPT rs1052553 were analysed in 907 individuals with PD and 1022 age-matched healthy controls from the Nigeria Parkinson's Disease Research network cohort. Clinical data related to PD included age at study, age at onset (AAO), and disease duration.

RESULTS

The frequency of the H1 haplotype was 98.7% in PD, and 99.1% in controls (p = 0.19). The H2 haplotype was present in - 1.3% of PD and 0.9% of controls (p = 0.24). The most frequent MAPT genotype was H1H1 (PD - 97.5%, controls - 98.2%). The H1 haplotype was not associated with PD risk after accounting for gender and AAO (Odds ratio for H1/H1 vs H1/H2 and H2/H2: 0.68 (95% CI:0.39-1.28); p = 0.23).

CONCLUSIONS

Our findings support previous studies that report a low frequency of the MAPT H2 haplotype in black ancestry Africans but document its occurrence in Nigerians. The MAPT H1 haplotype was not associated with an increased risk or age at onset of PD in this cohort.

MAPT allele and haplotype frequencies in Nigerian Africans: population distribution and association with Parkinson's disease risk and age at onset

BACKGROUND

The microtubule-associated protein tau ( MAPT ) gene is critical because of its putative role in the causal pathway of neurodegenerative diseases including Parkinson's disease (PD). However, there is a lack of clarity regarding the link between the main H1 haplotype and risk of PD. Inconsistencies in reported association may be driven by genetic variability in the populations studied to date. Data on MAPT haplotype frequencies in the general population and association studies exploring the role of MAPT haplotypes in conferring PD risk in black Africans are lacking.

OBJECTIVES

To determine the frequencies of MAPT haplotypes and explore the role of the H1 haplotype as a risk factor for PD risk and age at onset in Nigerian Africans.

METHODS

The haplotype and genotype frequencies of MAPT rs1052553 were analysed using PCR-based KASP™ in 907 individuals with PD and 1,022 age-matched neurologically normal controls from the Nigeria Parkinson's Disease Research (NPDR) network cohort. Clinical data related to PD included age at study, age at onset, and disease duration.

RESULTS

The frequency of the main MAPT H1 haplotype in this cohort was 98.7% in individuals with PD, and 99.1% in healthy controls (p=0.19). The H2 haplotype was present in 41/1929 (2.1%) of the cohort (PD - 1.3%; Controls - 0.9%; p=0.24). The most frequent MAPT genotype was H1H1 (PD - 97.5%, controls - 98.2%). The H1 haplotype was not associated with PD risk after accounting for gender and age at onset (Odds ratio for H1/H1 vs H1/H2 and H2/H2: 0.68 (95% CI:0.39-1.28); p=0.23).

CONCLUSIONS

Our findings support previous studies that report a low frequency of the MAPT H2 haplotype in black ancestry Africans, but document its occurrence in the Nigerian population (2.1%). In this cohort of black Africans with PD, the MAPT H1 haplotype was not associated with an increased risk or age at onset of PD.

Functional regulatory variants implicate distinct transcriptional networks in dementia

Predicting the function of noncoding variation is a major challenge in modern genetics. In this study, we used massively parallel reporter assays to screen 5706 variants identified from genome-wide association studies for both Alzheimer's disease (AD) and progressive supranuclear palsy (PSP), identifying 320 functional regulatory variants (frVars) across 27 loci, including the complex 17q21.31 region. We identified and validated multiple risk loci using CRISPR interference or excision, including complement 4 (C4A) and APOC1 in AD and PLEKHM1 and KANSL1 in PSP. Functional variants disrupt transcription factor binding sites converging on enhancers with cell type-specific activity in PSP and AD, implicating a neuronal SP1-driven regulatory network in PSP pathogenesis. These analyses suggest that noncoding genetic risk is driven by common genetic variants through their aggregate activity on specific transcriptional programs.