Whole genome sequencing study of identical twins discordant for psychosis.

Monozygotic (MZ) twins are often thought to have identical genomes, but recent work has shown that early post-zygotic events can result in a spectrum of DNA variants that are different between MZ twins. Such variants may explain phenotypic discordance and contribute to disease etiology. Here we performed whole genome sequencing in 17 pairs of MZ twins discordant for a psychotic disorder (schizophrenia, schizoaffective disorder or bipolar disorder). We examined various classes of rare variants that are discordant within a twin pair. We identified four genes harboring rare, predicted deleterious missense variants that were private to an affected individual in the cohort. Variants in FOXN1 and FLOT2 would have been categorized as damaging from recent schizophrenia and bipolar exome sequencing studies. Additionally, we identified four rare genic copy number variants (CNVs) private to an affected sample, two of which overlapped genes that have shown evidence for association with schizophrenia or bipolar disorder. One such CNV was a 3q29 duplication previously implicated in autism and developmental delay. We have performed the largest MZ twin study for discordant psychotic phenotypes to date. These findings warrant further investigation using other analytical approaches.

Investigating copy number variants in schizophrenia pedigrees using a new consensus pipeline called PECAN.

Copy number variants (CNVs) have been implicated in many human diseases, including psychiatric disorders. Whole genome sequencing offers advantages in CNV calling compared to previous array-based methods. Here we present a robust and transparent CNV calling pipeline, PECAN (PEdigree Copy number vAriaNt calling), for short-read, whole genome sequencing data, comprised of a novel combination of four calling methods and structural variant genotyping. This method is scalable and can incorporate pedigree information to retain lower-confidence CNVs that would otherwise be discarded. We have robustly benchmarked PECAN using gold-standard CNV calls for two well-established evaluation samples, NA12878 and HG002, showing that PECAN performs with high precision and recall on both datasets, outperforming another pedigree-based CNV calling pipeline. As part of this work, we provide a list of high-confidence gold standard CNVs for the NA12878 reference sample, curated from multiple studies. We applied PECAN to a collection of pedigrees multiply affected with schizophrenia and identified a rare deletion that perfectly co-segregates with schizophrenia in one of the pedigrees. The CNV overlaps the gene PITRM1, which has been implicated in a complex phenotype including ataxia, developmental delay, and schizophrenia-like episodes in affected adults.

The impact of perinatal maternal stress on the maternal and infant gut and human milk microbiomes: A scoping review protocol.

OBJECTIVE: The objective of this scoping review is to review the research evidence regarding the impact of perinatal maternal stress on the maternal and infant gut and human milk microbiomes. INTRODUCTION: Perinatal stress which refers to psychological stress experienced by individuals during pregnancy and the postpartum period is emerging as a public health concern. Early exposure of infants to perinatal maternal stress can potentially lead to metabolic, immune, and neurobehavioral disorders that extend into adulthood. The role of the gut and human milk microbiome in the microbiome-gut-brain axis as a mechanism of stress transfer has been previously reported. A transfer of colonised aberrant microbiota from mother to infant is proposed to predispose the infant to a pro- inflammatory microbiome with dysregulated metabolic process thereby initiating early risk of chronic diseases. The interplay of perinatal maternal stress and its relationship to the maternal and infant gut and human milk microbiome requires further systematic examination in the literature. INCLUSION CRITERIA: This scoping review is an exploratory mapping review which will focus on the population of mothers and infants with the exploration of the key concepts of maternal stress and its impact on the maternal and infant gut and human milk microbiome in the context of the perinatal period. It will focus on the pregnancy and the post-natal period up to 6 months with infants who are exclusively breastfed. METHODS: This study will be guided by the Joanna Briggs Institute's (JBI) methodology for scoping reviews along with use of the Prisma Scr reporting guideline. A comprehensive search will be conducted using the following databases, CINAHL Complete; MEDLINE; PsycINFO, Web of Science and Scopus. A search strategy with pre-defined inclusion and exclusion criteria will be used to retrieve peer reviewed data published in English from 2014 to present. Screening will involve a three-step process with screening tool checklists. Results will be presented in tabular and narrative summaries, covering thematic concepts and their relationships. This protocol is registered with Open Science Framework DOI 10.17605/OSF.IO/5SRMV.

Ultrarare Missense Variants Implicated in Utah Pedigrees Multiply Affected With Schizophrenia.

Background: Recent work from the Schizophrenia Exome Sequencing Meta-analysis (SCHEMA) consortium showed significant enrichment of ultrarare variants in schizophrenia cases. Family-based studies offer a unique opportunity to evaluate rare variants because risk in multiplex pedigrees is more likely to be influenced by the same collection of variants than an unrelated cohort. Methods: Here, we examine whole genome sequencing data from 35 individuals across 6 pedigrees multiply affected by schizophrenia. We applied a rigorous filtering pipeline to search for classes of protein-coding variants that cosegregated with disease status, and we examined these for evidence of enrichment in the SCHEMA dataset. Additionally, we applied a family-based consensus approach to call copy number variants and screen against a list of schizophrenia-associated risk variants. Results: We identified deleterious missense variants in 3 genes (ATP2B2, SLC25A28, and GSK3A) that cosegregated with disease in 3 of the pedigrees. In the SCHEMA, the gene ATP2B2 shows highly suggestive evidence for deleterious missense variants in schizophrenia cases (p = .000072). ATP2B2 is involved in intracellular calcium homeostasis, expressed in multiple brain tissue types, and predicted to be intolerant to loss-of-function and missense variants. Conclusions: We have identified genes that are likely to increase schizophrenia risk in 3 of the 6 pedigrees examined, the strongest evidence being for a gene involved in calcium homeostasis. Further work is required to examine other classes of variants that may be contributing to disease burden.

Outcomes of Tofacitinib Use in an Irish Pediatric Cohort.

Background: Pediatric ulcerative colitis (UC) is typically more extensive and severe at diagnosis compared with adult disease. Tofacitinib, a Janus kinase inhibitor, has been used since 2018 to induce and maintain remission in UC. There are limited pediatric data regarding its use, either as a monotherapy or in combination with other treatments. Objectives: To determine the real-world experience and outcomes of tofacitinib therapy in the Irish national cohort with pediatric UC. Methods: A retrospective study of tofacitinib outcomes was undertaken at Ireland's single national center for pediatric inflammatory bowel disease. All patients commenced on tofacitinib since its availability in 2019 were included. Baseline and follow-up clinical characteristics, phenotype, Pediatric Ulcerative Colitis Activity Index (PUCAI) scores, and treatments before and after tofacitinib commenced were recorded. The primary outcome was remission by 8 weeks, with other clinical outcomes being recorded to maximal available follow-up. Results: Between November 1, 2019 and June 30, 2022, 15 children (M:F 1:2) were prescribed tofacitinib, 5 as monotherapy. Thirteen had baseline pancolitis at diagnosis and all patients had prior infliximab exposure. The mean time from diagnosis to starting tofacitinib was 381 days (±SD 265). Dual therapy included 5 with infliximab, 4 with vedolizumab, and 1 with adalimumab. The average length of treatment on tofacitinib was 232 days (±SD 170) with 2 patients transitioning to adult services while in remission on tofacitinib therapy. The mean PUCAI score was 48.7 (±SD 14.1) pre-tofacitinib, 16.7 (±SD 15.6) at week 8, and 22.5 (±SD 29.6) by week 16, with a significant reduction in PUCAI by week 16 (P = 0.0004). Eight patients (3 monotherapy) achieved clinical remission, with 4 of the 5 dual therapy patients on infliximab. There were no significant outcome differences between those on mono- or dual therapy. Three patients with combined vedolizumab therapy did not achieve remission, 2 of whom required colectomy by week 24. There were no malignancies, 1 patient developed shingles and another developed herpangina post-tofacitinib. Failure to achieve clinical remission by week 16 was seen in all children who progressed to colectomy (n = 4). Conclusion: Combining tofacitinib with other biologics is effective in select children with refractory UC. Early responders were more likely to achieve a sustained response at week 16. Failure to achieve remission by week 16 of tofacitinib therapy was strongly associated with progression to colectomy.

Investigating the dark-side of the genome: a barrier to human disease variant discovery?

The human genome contains regions that cannot be adequately assembled or aligned using next generation short-read sequencing technologies. More than 2500 genes are known contain such 'dark' regions. In this study, we investigate the negative consequences of dark regions on gene discovery across a range of disease and study types, showing that dark regions are likely preventing researchers from identifying genetic variants relevant to human disease.

Evidence for parent-of-origin effects in autism spectrum disorder: a narrative review.

Autism spectrum disorder (ASD) is a heterogeneous group of early-onset neurodevelopmental disorders known to be highly heritable with a complex genetic architecture. Abnormal brain developmental trajectories that impact synaptic functioning, excitation-inhibition balance and brain connectivity are now understood to play a central role in ASD. Ongoing efforts to identify the genetic underpinnings still prove challenging, in part due to phenotypic and genetic heterogeneity.This review focuses on parent-of-origin effects (POEs), where the phenotypic effect of an allele depends on its parental origin. POEs include genomic imprinting, transgenerational effects, mitochondrial DNA, sex chromosomes and mutational transmission bias. The motivation for investigating these mechanisms in ASD has been driven by their known impacts on early brain development and brain functioning, in particular for the most well-documented POE, genomic imprinting. Moreover, imprinting is implicated in syndromes such as Angelman and Prader-Willi, which frequently share comorbid symptoms with ASD. In addition to other regions in the genome, this comprehensive review highlights the 15q11-q13 and 7q chromosomal regions as well as the mitochondrial DNA as harbouring the majority of currently identified POEs in ASD.

Identity-by-descent analysis of a large Tourette's syndrome pedigree from Costa Rica implicates genes involved in neuronal development and signal transduction.

Tourette Syndrome (TS) is a heritable, early-onset neuropsychiatric disorder that typically begins in early childhood. Identifying rare genetic variants that make a significant contribution to risk in affected families may provide important insights into the molecular aetiology of this complex and heterogeneous syndrome. Here we present a whole-genome sequencing (WGS) analysis from the 11-generation pedigree (>500 individuals) of a densely affected Costa Rican family which shares ancestry from six founder pairs. By conducting an identity-by-descent (IBD) analysis using WGS data from 19 individuals from the extended pedigree we have identified putative risk haplotypes that were not seen in controls, and can be linked with four of the six founder pairs. Rare coding and non-coding variants present on the haplotypes and only seen in haplotype carriers show an enrichment in pathways such as regulation of locomotion and signal transduction, suggesting common mechanisms by which the haplotype-specific variants may be contributing to TS-risk in this pedigree. In particular we have identified a rare deleterious missense variation in RAPGEF1 on a chromosome 9 haplotype and two ultra-rare deleterious intronic variants in ERBB4 and IKZF2 on the same chromosome 2 haplotype. All three genes play a role in neurodevelopment. This study, using WGS data in a pedigree-based approach, shows the importance of investigating both coding and non-coding variants to identify genes that may contribute to disease risk. Together, the genes and variants identified on the IBD haplotypes represent biologically relevant targets for investigation in other pedigree and population-based TS data.

A critical review of life cycle assessment studies of woody biomass conversion to sugars.

Woody biomass could potentially become a viable raw material for the future sustainable chemical industry. For this, a suitable regulatory framework must exist, that would create favourable economic conditions for wood biorefineries. Such policies must be developed on the basis of scientific evidence-in this case, data supporting the environmental advantages of the bio-based feedstocks to the chemical industry. The most suitable methodology for comprehensive evaluation of environmental performance of technologies is life cycle assessment (LCA). In this review, the available LCA studies of woody biomass fractionation and conversion to bulk chemical feedstocks are critically evaluated. It has been revealed that the majority of the openly available studies do not contain transparent inventory data and, therefore, cannot be verified or re-used; studies containing inventory data are reported in this review. The lack of inventory data also prevents comparison between studies of the same processes performed with different evaluation methods or using different system boundaries. Recommendations are proposed on how to overcome issues of commercial data sensitivity by using black-box modelling when reporting environmental information. From several comparable LCA studies, it has been concluded that today the most environmentally favourable technology for wood biomass fractionation is organosolv. This article is part of the theme issue 'Bio-derived and bioinspired sustainable advanced materials for emerging technologies (part 1)'.

Converting single nucleotide variants between genome builds: from cautionary tale to solution.

Next-generation sequencing studies are dependent on a high-quality reference genome for single nucleotide variant (SNV) calling. Although the two most recent builds of the human genome are widely used, position information is typically not directly comparable between them. Re-alignment gives the most accurate position information, but this procedure is often computationally expensive, and therefore, tools such as liftOver and CrossMap are used to convert data from one build to another. However, the positions of converted SNVs do not always match SNVs derived from aligned data, and in some instances, SNVs are known to change chromosome when converted. This is a significant problem when compiling sequencing resources or comparing results across studies. Here, we describe a novel algorithm to identify positions that are unstable when converting between human genome reference builds. These positions are detected independent of the conversion tools and are determined by the chain files, which provide a mapping of contiguous positions from one build to another. We also provide the list of unstable positions for converting between the two most commonly used builds GRCh37 and GRCh38. Pre-excluding SNVs at these positions, prior to conversion, results in SNVs that are stable to conversion. This simple procedure gives the same final list of stable SNVs as applying the algorithm and subsequently removing variants at unstable positions. This work highlights the care that must be taken when converting SNVs between genome builds and provides a simple method for ensuring higher confidence converted data. Unstable positions and algorithm code, available at https://github.com/cathaloruaidh/genomeBuildConversion.

Methyl-CpG-binding protein 2 mediates overlapping mechanisms across brain disorders.

MECP2 and its product, Methyl-CpG binding protein 2 (MeCP2), are mostly known for their association to Rett Syndrome (RTT), a rare neurodevelopmental disorder. Additional evidence suggests that MECP2 may underlie other neuropsychiatric and neurological conditions, and perhaps modulate common presentations and pathophysiology across disorders. To clarify the mechanisms of these interactions, we develop a method that uses the binding properties of MeCP2 to identify its targets, and in particular, the genes recognized by MeCP2 and associated to several neurological and neuropsychiatric disorders. Analysing mechanisms and pathways modulated by these genes, we find that they are involved in three main processes: neuronal transmission, immuno-reactivity, and development. Also, while the nervous system is the most relevant in the pathophysiology of the disorders, additional systems may contribute to MeCP2 action through its target genes. We tested our results with transcriptome analysis on Mecp2-null models and cells derived from a patient with RTT, confirming that the genes identified by our procedure are directly modulated by MeCP2. Thus, MeCP2 may modulate similar mechanisms in different pathologies, suggesting that treatments for one condition may be effective for related disorders.

Parental opinion of consent in neonatal research.

BACKGROUND: Neonatal research, particularly neonatal emergency research is a challenging area, notably in relation to obtaining valid prospective informed consent. The aim of this study is to determine parental perceptions of the consent process involved in performing research in newborn care, to explore methods used to obtain consent and their acceptability to parents. METHODS: A parental questionnaire was developed that examined attitudes towards research and hypothetical research studies, in which the acceptability of various methods of consent was examined (informed, waived, deferred). These research scenarios were of varying time sensitivity and perceived risk level. The study setting was an Irish maternity hospital. RESULTS: There were 600 responses to the questionnaire. In 93% of cases, parents felt that their involvement in the consent process was essential. In emergency situations, 52% felt full prospective informed consent was necessary; however, almost 28% of parents would feel pressure to consent. Most (75%) parents would prefer to be approached to discuss neonatal research studies antenatally, irrespective of study type and 40% of parents felt that neonates involved in research studies received overall better care. Acceptability of deferred consent was greater than waived, and was highest for the more emergency-based scenarios presented. DISCUSSION: Parents feel that they should play a central role in research involving their children. There were differences in the acceptability of various consent methods with strongest agreement for informed consent and lowest agreement for waived consent. Parents were more willing to accede to deferred consent in the cardiopulmonary resuscitation scenario study. These findings provide useful insights to consent strategies in future newborn research studies.

Independent Computerized Cognitive Remediation for Psychosis: An Investigation of Patient Experiences.

Cognitive remediation (CR) training improves cognition and functioning in patients with psychosis. To date, however, few studies have investigated CR from a subjective patient perspective. We recently conducted a randomized control trial demonstrating the effectiveness of a new, low therapist support, computer-based training program. This study aims to assess the service user experience of this program. Twenty CR completers with psychosis were interviewed using both rated and open-ended questions. Thematic analysis identified three broad themes: a) benefits of doing CR, b) costs of doing CR, and c) experience of doing CR. Positive experiences of therapy participation included improved cognition, improved positive self-regard, a development of life skills, and a transfer of benefits to everyday life. Negative experiences included therapy being difficult and tiring, leading to frustration and anxiety. These findings demonstrate the feasibility of low therapist support CR and point to variables that might impact therapy adherence.

DNA sequence-level analyses reveal potential phenotypic modifiers in a large family with psychiatric disorders.

Psychiatric disorders are a group of genetically related diseases with highly polygenic architectures. Genome-wide association analyses have made substantial progress towards understanding the genetic architecture of these disorders. More recently, exome- and whole-genome sequencing of cases and families have identified rare, high penetrant variants that provide direct functional insight. There remains, however, a gap in the heritability explained by these complementary approaches. To understand how multiple genetic variants combine to modify both severity and penetrance of a highly penetrant variant, we sequenced 48 whole genomes from a family with a high loading of psychiatric disorder linked to a balanced chromosomal translocation. The (1;11)(q42;q14.3) translocation directly disrupts three genes: DISC1, DISC2, DISC1FP and has been linked to multiple brain imaging and neurocognitive outcomes in the family. Using DNA sequence-level linkage analysis, functional annotation and population-based association, we identified common and rare variants in GRM5 (minor allele frequency (MAF) > 0.05), PDE4D (MAF > 0.2) and CNTN5 (MAF < 0.01) that may help explain the individual differences in phenotypic expression in the family. We suggest that whole-genome sequencing in large families will improve the understanding of the combined effects of the rare and common sequence variation underlying psychiatric phenotypes.

Altered DNA methylation associated with a translocation linked to major mental illness.

Recent work has highlighted a possible role for altered epigenetic modifications, including differential DNA methylation, in susceptibility to psychiatric illness. Here, we investigate blood-based DNA methylation in a large family where a balanced translocation between chromosomes 1 and 11 shows genome-wide significant linkage to psychiatric illness. Genome-wide DNA methylation was profiled in whole-blood-derived DNA from 41 individuals using the Infinium HumanMethylation450 BeadChip (Illumina Inc., San Diego, CA). We found significant differences in DNA methylation when translocation carriers (n = 17) were compared to related non-carriers (n = 24) at 13 loci. All but one of the 13 significant differentially methylated positions (DMPs) mapped to the regions surrounding the translocation breakpoints. Methylation levels of five DMPs were associated with genotype at SNPs in linkage disequilibrium with the translocation. Two of the five genes harbouring significant DMPs, DISC1 and DUSP10, have been previously shown to be differentially methylated in schizophrenia. Gene Ontology analysis revealed enrichment for terms relating to neuronal function and neurodevelopment among the genes harbouring the most significant DMPs. Differentially methylated region (DMR) analysis highlighted a number of genes from the MHC region, which has been implicated in psychiatric illness previously through genetic studies. We show that inheritance of a translocation linked to major mental illness is associated with differential DNA methylation at loci implicated in neuronal development/function and in psychiatric illness. As genomic rearrangements are over-represented in individuals with psychiatric illness, such analyses may be valuable more widely in the study of these conditions.

Identification of a novel locus on chromosome 2q13, which predisposes to clinical vertebral fractures independently of bone density.

OBJECTIVES: To identify genetic determinants of susceptibility to clinical vertebral fractures, which is an important complication of osteoporosis. METHODS: Here we conduct a genome-wide association study in 1553 postmenopausal women with clinical vertebral fractures and 4340 controls, with a two-stage replication involving 1028 cases and 3762 controls. Potentially causal variants were identified using expression quantitative trait loci (eQTL) data from transiliac bone biopsies and bioinformatic studies. RESULTS: A locus tagged by rs10190845 was identified on chromosome 2q13, which was significantly associated with clinical vertebral fracture (P=1.04×10-9) with a large effect size (OR 1.74, 95% CI 1.06 to 2.6). Bioinformatic analysis of this locus identified several potentially functional SNPs that are associated with expression of the positional candidate genes TTL (tubulin tyrosine ligase) and SLC20A1 (solute carrier family 20 member 1). Three other suggestive loci were identified on chromosomes 1p31, 11q12 and 15q11. All these loci were novel and had not previously been associated with bone mineral density or clinical fractures. CONCLUSION: We have identified a novel genetic variant that is associated with clinical vertebral fractures by mechanisms that are independent of BMD. Further studies are now in progress to validate this association and evaluate the underlying mechanism.

Genome-wide haplotype-based association analysis of major depressive disorder in Generation Scotland and UK Biobank.

Genome-wide association studies using genotype data have had limited success in the identification of variants associated with major depressive disorder (MDD). Haplotype data provide an alternative method for detecting associations between variants in weak linkage disequilibrium with genotyped variants and a given trait of interest. A genome-wide haplotype association study for MDD was undertaken utilising a family-based population cohort, Generation Scotland: Scottish Family Health Study (n = 18,773), as a discovery cohort with UK Biobank used as a population-based replication cohort (n = 25,035). Fine mapping of haplotype boundaries was used to account for overlapping haplotypes potentially tagging the same causal variant. Within the discovery cohort, two haplotypes exceeded genome-wide significance (P < 5 × 10-8) for an association with MDD. One of these haplotypes was nominally significant in the replication cohort (P < 0.05) and was located in 6q21, a region which has been previously associated with bipolar disorder, a psychiatric disorder that is phenotypically and genetically correlated with MDD. Several haplotypes with P < 10-7 in the discovery cohort were located within gene coding regions associated with diseases that are comorbid with MDD. Using such haplotypes to highlight regions for sequencing may lead to the identification of the underlying causal variants.

Balanced translocation linked to psychiatric disorder, glutamate, and cortical structure/function.

Rare genetic variants of large effect can help elucidate the pathophysiology of brain disorders. Here we expand the clinical and genetic analyses of a family with a (1;11)(q42;q14.3) translocation multiply affected by major psychiatric illness and test the effect of the translocation on the structure and function of prefrontal, and temporal brain regions. The translocation showed significant linkage (LOD score 6.1) with a clinical phenotype that included schizophrenia, schizoaffective disorder, bipolar disorder, and recurrent major depressive disorder. Translocation carriers showed reduced cortical thickness in the left temporal lobe, which correlated with general psychopathology and positive psychotic symptom severity. They showed reduced gyrification in prefrontal cortex, which correlated with general psychopathology severity. Translocation carriers also showed significantly increased activation in the caudate nucleus on increasing verbal working memory load, as well as statistically significant reductions in the right dorsolateral prefrontal cortex glutamate concentrations. These findings confirm that the t(1;11) translocation is associated with a significantly increased risk of major psychiatric disorder and suggest a general vulnerability to psychopathology through altered cortical structure and function, and decreased glutamate levels.

Preliminary assessment of pre-morbid DNA methylation in individuals at high genetic risk of mood disorders.

OBJECTIVES: Accumulating evidence implicates altered DNA methylation in psychiatric disorders, including bipolar disorder (BD) and major depressive disorder (MDD). It is not clear, however, whether these changes are causative or result from illness progression or treatment. To disentangle these possibilities we profiled genome-wide DNA methylation in well, unrelated individuals at high familial risk of mood disorder. DNA methylation was compared between individuals who subsequently developed BD or MDD [ill later (IL)] and those who remained well [well later (WL)]. METHODS: DNA methylation profiles were obtained from whole-blood samples from 22 IL and 23 WL individuals using the Infinium HumanMethylation450 BeadChip. Differential methylation was assessed on a single-locus and regional basis. Pathway analysis was performed to assess enrichment for particular biological processes amongst nominally significantly differentially methylated loci. RESULTS: Although no locus withstood correction for multiple testing, uncorrected P-values provided suggestive evidence for altered methylation at sites within genes previously implicated in neuropsychiatric conditions, such as Transcription Factor 4 (TCF4) and Interleukin 1 Receptor Accessory Protein-Like 1 ([IL1RAPL1]; P≤3.11×10(-5) ). Pathway analysis revealed significant enrichment for several neurologically relevant pathways and functions, including Nervous System Development and Function and Behavior; these findings withstood multiple testing correction (q≤0.05). Analysis of differentially methylated regions identified several within the major histocompatibility complex (P≤.000 479), a region previously implicated in schizophrenia and BD. CONCLUSIONS: Our data provide provisional evidence for the involvement of altered whole-blood DNA methylation in neurologically relevant genes in the aetiology of mood disorders. These findings are convergent with the findings of genome-wide association studies.

SuRFing the genomics wave: an R package for prioritising SNPs by functionality.

Identifying functional non-coding variants is one of the greatest unmet challenges in genetics. To help address this, we introduce an R package, SuRFR, which integrates functional annotation and prior biological knowledge to prioritise candidate functional variants. SuRFR is publicly available, modular, flexible, fast, and simple to use. We demonstrate that SuRFR performs with high sensitivity and specificity and provide a widely applicable and scalable benchmarking dataset for model training and validation. Website: http://www.cgem.ed.ac.uk/resources/