Rare variants in the endocytic pathway are associated with Alzheimer's disease, its related phenotypes, and functional consequences

Integrating rare pathogenic variant prioritization with gene-based association analysis to identify novel genes and relevant multimodal traits for Alzheimer's disease

INTRODUCTION

Increasing evidence has highlighted rare variants in Alzheimer's disease (AD). However, insufficient sample sizes, especially in underrepresented ethnic groups, hinder their investigation. Additionally, their impact on endophenotypes remains largely unexplored.

METHODS

We prioritized rare likely-deleterious variants based on whole-genome sequencing data from a Chinese AD cohort (n = 988). Gene-based optimal sequence kernel association tests were conducted between AD cases and normal controls to identify AD-related genes. Network clustering, endophenotype association, and cellular experiments were conducted to evaluate their functional consequences.

RESULTS

We identified 11 novel AD candidate genes, which captured AD-related pathways and enhanced AD risk prediction performance. Key genes (RABEP1, VIPR1, RPL3L, and CABIN1) were linked to cognitive decline and brain atrophy. Experiments showed RABEP1 p.R845W inducing endocytosis dysregulation and exacerbating toxic amyloid β accumulation, underscoring its therapeutic potential.

DISCUSSION

Our findings highlighted the contributions of rare variants to AD and provided novel insights into AD therapeutics.

HIGHLIGHTS

Identified 11 novel AD candidate genes in a Chinese AD cohort. Correlated candidate genes with AD-related cognitive and brain imaging traits. Indicated RABEP1 p.R845W as a critical AD contributor in the endocytic pathway.

ASACO: Automatic and Serial Analysis of CO-expression to discover gene modifiers with potential use in drug repurposing

Massive gene expression analyses are widely used to find differentially expressed genes under specific conditions. The results of these experiments are often available in public databases that are undergoing a growth similar to that of molecular sequence databases in the past. This now allows novel secondary computational tools to emerge that use such information to gain new knowledge. If several genes have a similar expression profile across heterogeneous transcriptomics experiments, they could be functionally related. These associations are usually useful for the annotation of uncharacterized genes. In addition, the search for genes with opposite expression profiles is useful for finding negative regulators and proposing inhibitory compounds in drug repurposing projects. Here we present a new web application, Automatic and Serial Analysis of CO-expression (ASACO), which has the potential to discover positive and negative correlator genes to a given query gene, based on thousands of public transcriptomics experiments. In addition, examples of use are presented, comparing with previous contrasted knowledge. The results obtained propose ASACO as a useful tool to improve knowledge about genes associated with human diseases and noncoding genes. ASACO is available at http://www.bioinfocabd.upo.es/asaco/.

Biostatistical Aspects of Whole Genome Sequencing Studies: Preprocessing and Quality Control

Rapid advances in high-throughput DNA sequencing technologies have enabled large-scale whole genome sequencing (WGS) studies. Before performing association analysis between phenotypes and genotypes, preprocessing and quality control (QC) of the raw sequence data need to be performed. Because many biostatisticians have not been working with WGS data so far, we first sketch Illumina's short-read sequencing technology. Second, we explain the general preprocessing pipeline for WGS studies. Third, we provide an overview of important QC metrics, which are applied to WGS data: on the raw data, after mapping and alignment, after variant calling, and after multisample variant calling. Fourth, we illustrate the QC with the data from the GENEtic SequencIng Study Hamburg-Davos (GENESIS-HD), a study involving more than 9000 human whole genomes. All samples were sequenced on an Illumina NovaSeq 6000 with an average coverage of 35× using a PCR-free protocol. For QC, one genome in a bottle (GIAB) trio was sequenced in four replicates, and one GIAB sample was successfully sequenced 70 times in different runs. Fifth, we provide empirical data on the compression of raw data using the DRAGEN original read archive (ORA). The most important quality metrics in the application were genetic similarity, sample cross-contamination, deviations from the expected Het/Hom ratio, relatedness, and coverage. The compression ratio of the raw files using DRAGEN ORA was 5.6:1, and compression time was linear by genome coverage. In summary, the preprocessing, joint calling, and QC of large WGS studies are feasible within a reasonable time, and efficient QC procedures are readily available.