Non-coding genetic variants outside of protein-coding genome regions play an important role in genetic and epigenetic regulation. It has become increasingly important to understand their roles, as non-coding variants often make up the majority of top findings of genome-wide association studies (GWAS). In addition, the growing popularity of disease-specific whole-genome sequencing (WGS) efforts expands the library of and offers unique opportunities for investigating both common and rare non-coding variants, which are typically not detected in more limited GWAS approaches. However, the sheer size and breadth of WGS data introduce additional challenges to predicting functional impacts in terms of data analysis and interpretation. This review focuses on the recent approaches developed for efficient, at-scale annotation and prioritization of non-coding variants uncovered in WGS analyses. In particular, we review the latest scalable annotation tools, databases and functional genomic resources for interpreting the variant findings from WGS based on both experimental data and in silico predictive annotations. We also review machine learning-based predictive models for variant scoring and prioritization. We conclude with a discussion of future research directions which will enhance the data and tools necessary for the effective functional analyses of variants identified by WGS to improve our understanding of disease etiology.

Cumulative evidence has revealed the association between mitochondrial dysfunction and Alzheimer’s disease (AD). Because the number of mitochondrial genes is very limited, the mitochondrial pathogenesis of AD must involve certain nuclear genes. In this study, we employed systems genetic methods to identify mitochondrion-associated nuclear genes that may participate in the pathogenesis of AD. First, we performed a mitochondrial genome-wide association study (MiWAS, n = 809) to identify mitochondrial single-nucleotide polymorphisms (MT-SNPs) associated with AD. Then, epistasis analysis was performed to examine interacting SNPs between the mitochondrial and nuclear genomes. Weighted co-expression network analysis (WGCNA) was applied to transcriptomic data from the same sample (n = 743) to identify AD-related gene modules, which were further enriched by mitochondrion-associated genes. Using hub genes derived from these modules, random forest models were constructed to predict AD risk in four independent datasets (n = 743, n = 542, n = 161, and n = 540). In total, 9 potentially significant MT-SNPs and 14,340 nominally significant MT-nuclear interactive SNPs were identified for AD, which were validated by functional analysis. A total of 6 mitochondrion-related modules involved in AD pathogenesis were found by WGCNA, from which 91 hub genes were screened and used to build AD risk prediction models. For the four independent datasets, these models perform better than those derived from AD genes identified by genome-wide association studies (GWASs) or differential expression analysis (DeLong’s test, p < 0.05). Overall, through systems genetics analyses, mitochondrion-associated SNPs/genes with potential roles in AD pathogenesis were identified and preliminarily validated, illustrating the power of mitochondrial genetics in AD pathogenesis elucidation and risk prediction.

Background: The 17q21.31 region with various structural forms characterized by the H1/H2 haplotypes and three large copy number variations (CNVs) represents the strongest risk locus in progressive supranuclear palsy (PSP).Objective: To investigate the association between CNVs and structural forms on 17q.21.31 with the risk of PSP.Methods: Utilizing whole genome sequencing data from 1684 PSP cases and 2392 controls, the three large CNVs (α, β, and γ) and structural forms within 17q21.31 were identified and analyzed for their association with PSP.Results: We found that the copy number of γ was associated with increased PSP risk (odds ratio [OR] = 1.10, P = 0.0018). From H1β1γ1 (OR = 1.21) and H1β2γ1 (OR = 1.24) to H1β1γ4 (OR = 1.57), structural forms of H1 with additional copies of γ displayed a higher risk for PSP. The frequency of the risk sub‐haplotype H1c rises from 1% in individuals with two γ copies to 88% in those with eight copies. Additionally, γ duplication up‐regulates expression of ARL17B, LRRC37A/LRRC37A2, and NSFP1, while down‐regulating KANSL1. Single‐nucleus RNA‐seq of the dorsolateral prefrontal cortex analysis reveals γ duplication primarily up‐regulates LRRC37A/LRRC37A2 in neuronal cells.Conclusions: The copy number of γ is associated with the risk of PSP after adjusting for H1/H2, indicating that the complex structure at 17q21.31 is an important consideration when evaluating the genetic risk of PSP. © 2025 The Author(s). Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Genetic risk scores (GRS) are crucial tools for estimating an individual's genetic liability to various traits and diseases, computed as a weighted sum of trait‐associated allele counts. Traditionally, GRS models assume additive, linear effects of risk variants. However, complex traits often involve nonadditive interactions, such as epistasis, which are not captured by these conventional methods. In this study, we investigate the use of random forest (RF) models as a model‐free approach for constructing GRS, leveraging RF's capacity to capture complex, nonlinear interactions among genetic variants. Specifically, we introduce two new RF‐based GRS strategies to boost RF performance and to incorporate base data information if available, including (1) ctRF, which optimizes linkage disequilibrium (LD) clumping and p ‐value thresholds within RF; and (2) wRF, which adjusts the chance of SNP inclusion in tree nodes based on their association strength. Through simulation studies and real data applications of Alzheimer's disease, body mass index, and atopy, we find that ctRF consistently outperforms other RF‐based methods and classical additive models when traits exhibit complex genetic architectures. Additionally, incorporating informative base data into RF‐GRS construction can enhance predictive accuracy. Our findings suggest that RF‐based GRS can effectively capture intricate genetic interactions, and offer a robust alternative to traditional GRS methods, especially for complex traits with nonlinear genetic effects.