Background. Pharmacogenetic (PGx) testing can guide drug prescribing but remains limited by the genomic assay used. Genotyping arrays are widely implemented yet limited to predefined variants, whereas low-pass whole-genome sequencing (LP-WGS) is not constrained by fixed probe design and may provide broader PGx variant availability after imputation. Methods. We compared Illumina Global Screening Array (GSA) v3 with ~1x LP-WGS for PGx profiling in 500 hospital biobank participants with electronic health record evidence of exposure to pharmacogenetically actionable drugs and reported adverse drug reactions. Concordance was evaluated genome-wide, at 20 actionable pharmacogenes for PharmCAT-derived star alleles and metabolizer phenotypes, and for HLA alleles. Results. Genome-wide concordance between imputed array and LP-WGS data was high (median 99.63%; interquartile range, 99.59%-99.64%). For pharmacogenetically relevant variants, LP-WGS captured a larger fraction, particularly rare alleles absent from the array data, whilst maintaining high concordance at shared sites. Predicted phenotype concordance exceeded 98% for most genes, although gene-specific differences in phenotype classification were observed. LP-WGS reduced missing phenotype assignments for selected loci, particularly CYP2C19 and NAT2, by improving resolution of star-allele structure. However, in structurally complex or incompletely characterized genes such as CYP2C9 and CYP2D6, broader variant recovery increased indeterminate classifications rather than consistently improving clinical interpretability. For HLA loci, concordance varied by imputation strategy, with SNP2HLA performing marginally better utilizing the GSA array compared to the LP-WGS approach. Conclusions. Overall, LP-WGS provides broader variant coverage and improved resolution for selected pharmacogenes but did not resolve all clinically important loci. These findings support further evaluation of LP-WGS as a scalable PGx screening approach, especially where long-term genomic data reuse is a priority.
F. Hodel, C. Thorball, D. Haefliger et al.· medRxiv· 0 citations
Mutational biases can influence genome composition, but their contribution to protein evolution remains difficult to quantify. Here we utilize a nearly neutral framework that translates nucleotide mutational spectra into expected amino acid substitution patterns and equilibrium amino acid compositions. Using SARS-CoV-2 as a model system, we show that the viral mutational spectrum explains more than 50% of the variation in observed single-nucleotide amino acid substitutions and predicts the overall direction of proteome-wide amino acid composition change during the COVID-19 pandemic. The predictive power of the model varies with selection regime: effectively neutral and weakly deleterious substitutions conform most closely to the mutational expectation, whereas strongly constrained sites and mutational hotspots show larger deviations. This indicates that departures from the nearly neutral baseline provide a quantitative proxy for purifying and positive selection. Extending the analysis across 34 RNA virus species, we find that positive-sense, negative-sense and double-stranded RNA viruses differ systematically in their mutational spectra, and that these differences are associated with predictable shifts in proteome composition. The same relationship is detectable in RNA-dependent RNA polymerase sequences from more than 77,000 viral species. These results indicate that taxon-specific mutational bias contributes persistently to protein evolution across evolutionary scales.
B. Efimenko, Alexander Voronka, Victoriya Skripskaya et al.· bioRxiv· 1 citation