Evaluation of gene-based family-based methods to detect novel genes associated with familial late onset Alzheimer disease

AbstractGene-based tests to study the combined effect of rare variants towards a particular phenotype have been widely developed for case-control studies, but their evolution and adaptation for family-based studies, especially for complex incomplete families, has been slower. In this study, we have performed a practical examination of all the latest gene-based methods available for family-based study designs using both simulated and real datasets. We have examined the performance of several collapsing, variance-component and transmission disequilibrium tests across eight different software and twenty-two models utilizing a cohort of 285 families (N=1,235) with late-onset Alzheimer disease (LOAD). After a thorough examination of each of these tests, we propose a methodological approach to identify, with high confidence, genes associated with the studied phenotype with high confidence and we provide recommendations to select the best software and model for family-based gene-based analyses. Additionally, in our dataset, we identified PTK2B, a GWAS candidate gene for sporadic AD, along with six novel genes (CHRD, CLCN2, HDLBP, CPAMD8, NLRP9, MAS1L) as candidates genes for familial LOAD.</jats:p

An atlas of cortical circular RNA expression in Alzheimer disease brains demonstrates clinical and pathological associations.

Parietal cortex RNA-sequencing (RNA-seq) data were generated from individuals with and without Alzheimer disease (AD; ncontrol = 13; nAD = 83) from the Knight Alzheimer Disease Research Center (Knight ADRC). Using this and an independent (Mount Sinai Brain Bank (MSBB)) AD RNA-seq dataset, cortical circular RNA (circRNA) expression was quantified in the context of AD. Significant associations were identified between circRNA expression and AD diagnosis, clinical dementia severity and neuropathological severity. It was demonstrated that most circRNA-AD associations are independent of changes in cognate linear messenger RNA expression or estimated brain cell-type proportions. Evidence was provided for circRNA expression changes occurring early in presymptomatic AD and in autosomal dominant AD. It was also observed that AD-associated circRNAs co-expressed with known AD genes. Finally, potential microRNA-binding sites were identified in AD-associated circRNAs for miRNAs predicted to target AD genes. Together, these results highlight the importance of analyzing non-linear RNAs and support future studies exploring the potential roles of circRNAs in AD pathogenesis

Hypoxia Reduces Arylsulfatase B Activity and Silencing Arylsulfatase B Replicates and Mediates the Effects of Hypoxia

This report presents evidence of 1) a role for arylsulfatase B (ARSB; N-acetylgalactosamine-4-sulfatase) in mediating intracellular oxygen signaling; 2) replication between the effects of ARSB silencing and hypoxia on sulfated glycosaminoglycan content, cellular redox status, and expression of hypoxia-associated genes; and 3) a mechanism whereby changes in chondroitin-4-sulfation that follow either hypoxia or ARSB silencing can induce transcriptional changes through galectin-3. ARSB removes 4-sulfate groups from the non-reducing end of chondroitin-4-sulfate and dermatan sulfate and is required for their degradation. For activity, ARSB requires modification of a critical cysteine residue by the formylglycine generating enzyme and by molecular oxygen. When primary human bronchial and human colonic epithelial cells were exposed to 10% O2×1 h, ARSB activity declined by ∼41% and ∼30% from baseline, as nuclear hypoxia inducible factor (HIF)-1α increased by ∼53% and ∼37%. When ARSB was silenced, nuclear HIF-1α increased by ∼81% and ∼61% from baseline, and mRNA expression increased to 3.73 (±0.34) times baseline. Inversely, ARSB overexpression reduced nuclear HIF-1α by ∼37% and ∼54% from baseline in the epithelial cells. Hypoxia, like ARSB silencing, significantly increased the total cellular sulfated glycosaminoglycans and chondroitin-4-sulfate (C4S) content. Both hypoxia and ARSB silencing had similar effects on the cellular redox status and on mRNA expression of hypoxia-associated genes. Transcriptional effects of both ARSB silencing and hypoxia may be mediated by reduction in galectin-3 binding to more highly sulfated C4S, since the galectin-3 that co-immunoprecipitated with C4S declined and the nuclear galectin-3 increased following ARSB knockdown and hypoxia

Circular RNA detection identifies circPSEN1 alterations in brain specific to autosomal dominant Alzheimer's disease

Background: Autosomal-dominant Alzheimer's disease (ADAD) is caused by pathogenic mutations in APP, PSEN1, and PSEN2, which usually lead to an early age at onset (< 65). Circular RNAs are a family of non-coding RNAs highly expressed in the nervous system and especially in synapses. We aimed to investigate differences in brain gene expression of linear and circular transcripts from the three ADAD genes in controls, sporadic AD, and ADAD. Methods: We obtained and sequenced RNA from brain cortex using standard protocols. Linear counts were obtained using the TOPMed pipeline; circular counts, using python package DCC. After stringent quality control (QC), we obtained the counts for PSEN1, PSEN2 and APP genes. Only circPSEN1 passed QC. We used DESeq2 to compare the counts across groups, correcting for biological and technical variables. Finally, we performed in-silico functional analyses using the Circular RNA interactome website and DIANA mirPath software. Results: Our results show significant differences in gene counts of circPSEN1 in ADAD individuals, when compared to sporadic AD and controls (ADAD = 21, AD = 253, Controls = 23-ADADvsCO: log2FC = 0.794, p = 1.63 × 10-04, ADADvsAD: log2FC = 0.602, p = 8.22 × 10-04). The high gene counts are contributed by two circPSEN1 species (hsa_circ_0008521 and hsa_circ_0003848). No significant differences were observed in linear PSEN1 gene expression between cases and controls, indicating that this finding is specific to the circular forms. In addition, the high circPSEN1 levels do not seem to be specific to PSEN1 mutation carriers; the counts are also elevated in APP and PSEN2 mutation carriers. In-silico functional analyses suggest that circPSEN1 is involved in several pathways such as axon guidance (p = 3.39 × 10-07), hippo signaling pathway (p = 7.38 × 10-07), lysine degradation (p = 2.48 × 10-05) or Wnt signaling pathway (p = 5.58 × 10-04) among other KEGG pathways. Additionally, circPSEN1 counts were able to discriminate ADAD from sporadic AD and controls with an AUC above 0.70. Conclusions: Our findings show the differential expression of circPSEN1 is increased in ADAD. Given the biological function previously ascribed to circular RNAs and the results of our in-silico analyses, we hypothesize that this finding might be related to neuroinflammatory events that lead or that are caused by the accumulation of amyloid-beta

RA-MAP, molecular immunological landscapes in early rheumatoid arthritis and healthy vaccine recipients

Rheumatoid arthritis (RA) is a chronic inflammatory disorder with poorly defined aetiology characterised by synovial inflammation with variable disease severity and drug responsiveness. To investigate the peripheral blood immune cell landscape of early, drug naive RA, we performed comprehensive clinical and molecular profiling of 267 RA patients and 52 healthy vaccine recipients for up to 18 months to establish a high quality sample biobank including plasma, serum, peripheral blood cells, urine, genomic DNA, RNA from whole blood, lymphocyte and monocyte subsets. We have performed extensive multi-omic immune phenotyping, including genomic, metabolomic, proteomic, transcriptomic and autoantibody profiling. We anticipate that these detailed clinical and molecular data will serve as a fundamental resource offering insights into immune-mediated disease pathogenesis, progression and therapeutic response, ultimately contributing to the development and application of targeted therapies for RA.</p

Evaluation of Gene-Based Family-Based Methods to Detect Novel Genes Associated With Familial Late Onset Alzheimer Disease

Gene-based tests to study the combined effect of rare variants on a particular phenotype have been widely developed for case-control studies, but their evolution and adaptation for family-based studies, especially studies of complex incomplete families, has been slower. In this study, we have performed a practical examination of all the latest gene-based methods available for family-based study designs using both simulated and real datasets. We examined the performance of several collapsing, variance-component, and transmission disequilibrium tests across eight different software packages and 22 models utilizing a cohort of 285 families (N = 1,235) with late-onset Alzheimer disease (LOAD). After a thorough examination of each of these tests, we propose a methodological approach to identify, with high confidence, genes associated with the tested phenotype and we provide recommendations to select the best software and model for family-based gene-based analyses. Additionally, in our dataset, we identified PTK2B, a GWAS candidate gene for sporadic AD, along with six novel genes (CHRD, CLCN2, HDLBP, CPAMD8, NLRP9, and MAS1L) as candidate genes for familial LOAD

Analysis of neurodegenerative Mendelian genes in clinically diagnosed Alzheimer Disease.

Alzheimer disease (AD), Frontotemporal lobar degeneration (FTD), Amyotrophic lateral sclerosis (ALS) and Parkinson disease (PD) have a certain degree of clinical, pathological and molecular overlap. Previous studies indicate that causative mutations in AD and FTD/ALS genes can be found in clinical familial AD. We examined the presence of causative and low frequency coding variants in the AD, FTD, ALS and PD Mendelian genes, in over 450 families with clinical history of AD and over 11,710 sporadic cases and cognitive normal participants from North America. Known pathogenic mutations were found in 1.05% of the sporadic cases, in 0.69% of the cognitively normal participants and in 4.22% of the families. A trend towards enrichment, albeit non-significant, was observed for most AD, FTD and PD genes. Only PSEN1 and PINK1 showed consistent association with AD cases when we used ExAC as the control population. These results suggest that current study designs may contain heterogeneity and contamination of the control population, and that current statistical methods for the discovery of novel genes with real pathogenic variants in complex late onset diseases may be inadequate or underpowered to identify genes carrying pathogenic mutations

Table2.DOCX

<p>Gene-based tests to study the combined effect of rare variants on a particular phenotype have been widely developed for case-control studies, but their evolution and adaptation for family-based studies, especially studies of complex incomplete families, has been slower. In this study, we have performed a practical examination of all the latest gene-based methods available for family-based study designs using both simulated and real datasets. We examined the performance of several collapsing, variance-component, and transmission disequilibrium tests across eight different software packages and 22 models utilizing a cohort of 285 families (N = 1,235) with late-onset Alzheimer disease (LOAD). After a thorough examination of each of these tests, we propose a methodological approach to identify, with high confidence, genes associated with the tested phenotype and we provide recommendations to select the best software and model for family-based gene-based analyses. Additionally, in our dataset, we identified PTK2B, a GWAS candidate gene for sporadic AD, along with six novel genes (CHRD, CLCN2, HDLBP, CPAMD8, NLRP9, and MAS1L) as candidate genes for familial LOAD.</p