Shannon K. McDonnell

Summary: Long-read nanopore sequencing enables simultaneous detection of germline variation and native DNA base modifications on individual DNA molecules, providing a unique opportunity to investigate allele-specific epigenetic regulation. Here, we performed whole-genome nanopore sequencing on normal and tumor prostate tissues to characterize differential methylation, methylation entropy, and allele-specific methylation (ASM) associated with noncoding genetic variants. Genome-wide analysis identified extensive cancer-associated differentially methylated regions (DMRs), with hypermethylated DMRs significantly enriched near transcription start sites and transcriptional regulatory regions. Integration with transcriptomic datasets revealed strong inverse relationships between promoter methylation and gene expression, while 5-hydroxymethylcytosine (5hmC) levels positively correlated with transcriptional activity across gene bodies. Using fragment-level methylation patterns enabled by long-read sequencing, we further quantified methylation entropy incorporating both 5mCG and 5hmCG states. Cancer-hypermethylated DMRs exhibited markedly reduced entropy, consistent with clonal fixation of methylation states during tumor progression. Entropy profiling across chromatin annotations demonstrated maximal epigenetic heterogeneity at partially modified enhancer-associated regions. To investigate cis-regulatory genetic effects, we developed a simple ASM framework (nanoASM) that can partition sequencing reads by allelic state and identifies allele-specific DMRs directly from long-read data. Compared with conventional population-level mQTL analysis, ASM demonstrated substantially improved statistical efficiency by leveraging within-individual contrasts and reducing sample-level heterogeneity. Although germline single nucleotide polymorphisms (SNPs) were largely shared between normal and tumor tissues, ASM patterns differed substantially, with tumor-associated ASM regions displaying significantly larger genomic span and stronger allelic methylation differences. Comparative analysis with TCGA prostate mQTL and GTEx prostate eQTL datasets demonstrated substantial concordance between ASM directionality and downstream transcriptional effects, particularly for variants located within DMRs and near transcription start sites. At the IRX4 prostate cancer risk locus, ASM identified an androgen-responsive regulatory domain overlapping AR ChIP-seq and H3K27ac peaks, nominating rs6885084 as a candidate functional variant. At the PSCA locus, ASM anchored by rs4736369 was associated with allele-specific methylation, chromatin activation, transcript abundance, and isoform usage. Together, these findings establish nanopore-based ASM analysis as a powerful approach for resolving functional noncoding variants and their regulatory domains they control in prostate cancer.

Key Points: Circulating proteins independently associated with kidney stone diagnosis are related to kidney stone matrix. Lower circulating uromodulin and higher scavenger receptor cysteine-rich domain-containing group B protein levels were associated with higher kidney stone risk, whereas kidney stone presence elevates plasma matrix metalloproteinase 7 level. The protective effect of uromodulin against kidney stones was independent of kidney function. Background: Kidney stones are increasingly recognized as a systemic disorder with a high global prevalence. However, large proteomics studies are lacking. Methods: An individual-level proteomics study with rigorous adjustments was performed on 35,331 UK Biobank participants to uncover the independent associations between 2922 circulating proteins and prevalent kidney stone disease. Mendelian randomization analysis was used to assess causal relationships. Findings were validated using genomic data from the Mayo Clinic Biobank ( N =43,744), concentration-response analysis, transcriptomics analyses, and additional genome-wide association studies. Results: Nine plasma proteins were independently associated with a kidney stone diagnosis, including reduced uromodulin (UMOD; β , −0.10; 95% confidence interval [CI], −0.15 to −0.05) and elevated scavenger receptor cysteine-rich domain-containing group B protein (SSC4D; β , 0.28; 95% CI, 0.17 to 0.38) and were enriched in extracellular matrix pathways. Mendelian randomization analysis revealed that the presence of kidney stone contributed to elevated levels of matrix metalloproteinase 7. Conversely, lower plasma UMOD (odds ratio [OR], 0.93; 95% CI, 0.90 to 0.97) and higher plasma SSC4D (OR, 1.10; 95% CI, 1.02 to 1.18) were associated with kidney stone risk. These associations were consistently replicated in the Mayo Clinic Biobank dataset (UMOD; OR, 0.92; 95% CI, 0.86 to 0.98; SSC4D; OR, 1.13; 95% CI, 1.01 to 1.27) and further validated by a concentration-response analysis. Single-nucleus RNA sequencing and quantitative trait loci analyses confirmed consistent associations between thick ascending limb UMOD expression and stone former status and with blood and urine UMOD concentrations. Genome-wide association study analysis, adjusted for eGFR, suggested that the protective role of UMOD against kidney stones was independent of kidney function. Conclusions: This study highlights significant associations between concentrations of specific blood proteins and a history of kidney stones. Several implicated proteins are related to kidney stone matrix, with UMOD independently associated with lower risk and SSC4D with higher risk of kidney stones.

Frontotemporal lobar degeneration with neuronal inclusions of the TAR DNA-binding protein 43 (FTLD-TDP) is a fatal neurodegenerative disorder with only a limited number of risk loci identified. We report our comprehensive genome-wide association study as part of the International FTLD-TDP Whole-Genome Sequencing Consortium, including 985 patients and 3,153 controls compiled from 26 institutions/brain banks in North America, Europe and Australia, and meta-analysis with the Dementia-seq cohort. We confirm UNC13A as the strongest overall FTLD-TDP risk factor and identify TNIP1 as a novel FTLD-TDP risk factor. In subgroup analyzes, we further identify genome-wide significant loci specific to each of the three main FTLD-TDP pathological subtypes (A, B and C), as well as enrichment of risk loci in distinct tissues, brain regions, and neuronal subtypes, suggesting distinct disease aetiologies in each of the subtypes. Rare variant analysis confirmed TBK1 and identified C3AR1, SMG8, VIPR1, RBPJL, L3MBTL1 and ANO9, as novel subtype-specific FTLD-TDP risk genes, further highlighting the role of innate and adaptive immunity and notch signaling pathway in FTLD-TDP, with potential diagnostic and novel therapeutic implications.

Key Points: Loss-of-function ALG8 and ALG9 variants were enriched in polycystic kidney/liver groups and International Classification of Diseases–coded cystic individuals in population cohorts. The ALG8 and ALG9 kidney phenotypes were usually mild to moderate, and lower eGFR or kidney failure was rare. ALG8 pathogenic variants sometimes resulted in severe polycystic liver disease. Background: Autosomal dominant polycystic kidney disease (ADPKD) is a common, inherited nephropathy often resulting in kidney failure. It is genetically heterogeneous; along with the major genes, PKD1 and PKD2 , at least eight others have been suggested. ALG8 pathogenic variants have been associated with autosomal dominant polycystic liver disease and implicated in ADPKD, while ALG9 has been suggested as an ADPKD gene, but details of the phenotypes and penetrance are unclear. Methods: We screened >3900 families with cystic kidneys and/or livers using global approaches to detect ALG8 or ALG9 pathogenic variants. In addition, population cohorts with sequence data (Genomics England 100K Genomics Project, UK Biobank, and Mayo Clinic Biobank [MCBB]) were screened for ALG8 / ALG9 pathogenic variants. Results: Multicenter screening of individuals with polycystic kidney and/or liver disease identified 51 (1.3%) ALG8 (7 multiplex) and 23 (0.6%) ALG9 (5 multiplex) families—frequencies that were approximately 10× and approximately 24× greater than nonpolycystic kidney disease controls. Analysis of individuals with polycystic kidney disease phenotypes in 100K Genomics Project, UK Biobank, and MCBB identified nine ALG8 (0.39%) and nine ALG9 (0.39%) families, an enriched frequency over controls. Two individuals had PKD1 and ALG8 pathogenic changes. Eighty-nine percent of individuals with ALG8 mutations with imaging in the entire MCBB had kidney cysts (50%, >10 cysts), with greater median kidney and liver cyst numbers than controls. For ALG9, 78% had kidney cysts (27%, >10 cysts). Individuals with ALG8 mutations typically had mild cystic kidneys with limited enlargement. Liver cysts were common (71%), with enlarged livers (>2L) found in 11 of 62 patients, although surgical intervention was rare. The ALG9 kidney phenotype was also of mild cystic kidneys, but enlarged livers were rare; for both genes, CKD or kidney failure were rare. Conclusions: ALG8 and ALG9 are defined as cystic kidney/liver genes but with limited penetrance for lower eGFR.

Key Points: Common kidney stones are unlikely to be an independent and direct cause of CKD in the general population. CKD may protect against kidney stones because of changes in key urinary factors critical for stone formation. Background: Kidney stones and CKD are common disorders with a substantial interaction. Although observational studies have suggested a potential for enhanced CKD risk after prior kidney stones, the exact relationship remains ambiguous. Methods: Shared comorbidities between two diseases were identified using unbiased screening. Genome-wide association study summary statistics were obtained from the UK Biobank (UKBB), FinnGen, and CKDGen, followed by genetic association analyses across various traits. Bidirectional Mendelian randomization (MR) analyses were performed to define causal links, complemented by multivariable MR that included the shared comorbidities including hypertension, diabetes, and obesity. Observational analyses were undertaken using cohorts from the Mayo Clinic and a UKBB subset. Results: Despite identifying a total of 123 conditions as shared comorbidities, there was no significant genetic correlation between kidney stones and CKD. Unadjusted MR analysis revealed no significant association between kidney stones and CKD risk (UKBB [exposure]/FinnGen [outcome]: odds ratio [OR]=0.97, 95% confidence interval [CI], 0.88 to 1.06; FinnGen/UKBB: OR=1.17, 95% CI, 0.98 to 1.39). Kidney stones did significantly associate with a higher urinary albumin-creatinine ratio (β=0.014, 95% CI, 0.002 to –0.025), but this association disappeared in the multivariable MR model (β=0.009, 95% CI, −0.003 to 0.020). Furthermore, in a cross-sectional analysis limited to the UKBB cohort, a robust regression model did not detect an independent association between kidney stones and urinary albumin-creatinine ratio (β=0.16, 95% CI, −0.04 to 0.35) or eGFR (β=0.10, 95% CI, −0.07 to 0.28). Conversely, CKD associated with a diminished risk of kidney stones in multivariable MR models (UKBB/FinnGen: OR=0.77, 95% CI, 0.69 to 0.87; FinnGen/UKBB: OR=0.73, 95% CI, 0.66 to 0.81). Furthermore, in the Mayo Clinic cohort with available urinary biochemistries, lower eGFR was associated with lower urinary calcium excretion and urinary calcium oxalate/phosphate supersaturation. Conclusions: In this study, kidney stones were not independently associated with CKD. Conversely, CKD was associated with a lower risk of calcium kidney stones likely via changes in key urinary traits, including lower calcium excretion.

We aimed to determine the effect of sample size on performance of polygenic hazard score (PHS) models in predicting the age at onset of prostate cancer. Age and genotypes were obtained for 40,861 men from the PRACTICAL consortium. The dataset included 201,590 SNPs per subject, and was split into training (34,444 samples) and testing (6,417 samples) sets. Two PHS model-building strategies were investigated. Established-SNP model considered 65 SNPs that had been associated with prostate cancer in the literature. A stepwise SNP selection was used to develop Discovery-SNP models. The performance of each PHS model was calculated for random sizes of the training set (1 to 30 thousand). The performance of a representative Established-SNP model was estimated for random sizes of the testing set (0.5 to 6 thousand). Mean HR _98/50 (hazard ratio of top 2% to the average in the test set) of the Established-SNP model increased from 1.73[95%CI: 1.69-1.77] to 2.41[2.40-2.43] when the number of training samples was increased from 1 to 30 thousand. The corresponding HR _98/50 of the Discovery-SNP model increased from 1.05[0.93-1.18] to 2.19[2.16-2.23]. HR _98/50 of a representative Established-SNP model using testing set sample sizes of 0.6 and 6 thousand observations were 1.78[1.70-1.85] and 1.73[1.71-1.76], respectively. We estimate that a study population of 20 to 30 thousand men is required to develop Discovery-SNP PHS models for prostate cancer. The required sample size could be reduced to 10 thousand samples, if a set of SNPs associated with the disease has already been established. Author summary: Polygenic hazard scores represent a recent advancement in polygenic prediction to model the age of onset of various diseases, such as Alzheimer’s disease or prostate cancer. These scores accumulate small effect sizes from several tens of genetic variants and can be used to establish an individual’s risk of experiencing an event relative to a control population across time. The largest barrier to the development of polygenic hazard scores is the large number of study subjects needed to develop the underlying models. We sought to understand the effect of varying the total number of samples on the performance of a polygenic hazard score in the context of prostate cancer. We found that the performance of the score did not appreciably change beyond 20 to 30 thousand observations when developing the model from scratch. However, when the discovery of the genetic variants can be borrowed from those already identified in the literature to be associated with the disease, the required number of samples is reduced to 10 thousand with no appreciable detriment in performance. We hope that these results can guide the design of future studies of polygenic scores in other diseases and demonstrate the importance of genome-wide association studies.

Prostate cancer is a polygenic disease with a large heritable component. A number of common, low-penetrance prostate cancer risk loci have been identified through GWAS. Here we apply the Bayesian multivariate variable selection algorithm JAM to fine-map 84 prostate cancer susceptibility loci, using summary data from a large European ancestry meta-analysis. We observe evidence for multiple independent signals at 12 regions and 99 risk signals overall. Only 15 original GWAS tag SNPs remain among the catalogue of candidate variants identified; the remainder are replaced by more likely candidates. Biological annotation of our credible set of variants indicates significant enrichment within promoter and enhancer elements, and transcription factor-binding sites, including AR, ERG and FOXA1. In 40 regions at least one variant is colocalised with an eQTL in prostate cancer tissue. The refined set of candidate variants substantially increase the proportion of familial relative risk explained by these known susceptibility regions, which highlights the importance of fine-mapping studies and has implications for clinical risk profiling.

Genome-wide association studies (GWAS) and fine-mapping efforts to date have identified more than 100 prostate cancer (PrCa)-susceptibility loci. We meta-analyzed genotype data from a custom high-density array of 46,939 PrCa cases and 27,910 controls of European ancestry with previously genotyped data of 32,255 PrCa cases and 33,202 controls of European ancestry. Our analysis identified 62 novel loci associated (P < 5.0 × 10⁻⁸) with PrCa and one locus significantly associated with early-onset PrCa (≤55 years). Our findings include missense variants rs1800057 (odds ratio (OR) = 1.16; P = 8.2 × 10⁻⁹; G>C, p.Pro1054Arg) in ATM and rs2066827 (OR = 1.06; P = 2.3 × 10⁻⁹; T>G, p.Val109Gly) in CDKN1B. The combination of all loci captured 28.4% of the PrCa familial relative risk, and a polygenic risk score conferred an elevated PrCa risk for men in the ninetieth to ninety-ninth percentiles (relative risk = 2.69; 95% confidence interval (CI): 2.55–2.82) and first percentile (relative risk = 5.71; 95% CI: 5.04–6.48) risk stratum compared with the population average. These findings improve risk prediction, enhance fine-mapping, and provide insight into the underlying biology of PrCa¹.

After decades of identifying risk factors using array-based genome-wide association studies (GWAS), genetic research of complex diseases has shifted to sequencing-based rare variants discovery. This requires large sample sizes for statistical power and has brought up questions about whether the current variant calling practices are adequate for large cohorts. It is well-known that there are discrepancies between variants called by different pipelines, and that using a single pipeline always misses true variants exclusively identifiable by other pipelines. Nonetheless, it is common practice today to call variants by one pipeline due to computational cost and assume that false negative calls are a small percent of total. We analyzed 10,000 exomes from the Alzheimer’s Disease Sequencing Project (ADSP) using multiple analytic pipelines consisting of different read aligners and variant calling strategies. We compared variants identified by using two aligners in 50,100, 200, 500, 1000, and 1952 samples; and compared variants identified by adding single-sample genotyping to the default multi-sample joint genotyping in 50,100, 500, 2000, 5000 and 10,000 samples. We found that using a single pipeline missed increasing numbers of high-quality variants correlated with sample sizes. By combining two read aligners and two variant calling strategies, we rescued 30% of pass-QC variants at sample size of 2000, and 56% at 10,000 samples. The rescued variants had higher proportions of low frequency (minor allele frequency [MAF] 1–5%) and rare (MAF < 1%) variants, which are the very type of variants of interest. In 660 Alzheimer’s disease cases with earlier onset ages of ≤65, 4 out of 13 (31%) previously-published rare pathogenic and protective mutations in APP, PSEN1, and PSEN2 genes were undetected by the default one-pipeline approach but recovered by the multi-pipeline approach. Identification of the complete variant set from sequencing data is the prerequisite of genetic association analyses. The current analytic practice of calling genetic variants from sequencing data using a single bioinformatics pipeline is no longer adequate with the increasingly large projects. The number and percentage of quality variants that passed quality filters but are missed by the one-pipeline approach rapidly increased with sample size.