Gene-based Tests Research Articles

COX15 deficiency causes oocyte ferroptosis

Mitochondria play diverse roles in mammalian physiology. The architecture, activity, and physiological functions of mitochondria in oocytes are largely different from those in somatic cells, but the mitochondrial proteins related to oocyte quality and reproductive longevity remain largely unknown. Here, using whole-exome sequencing data from 1,024 women (characterized by oocyte maturation arrest and degenerated or morphologically abnormal oocytes) and 2,868 healthy controls, we performed a population and gene-based burden test for mitochondrial genes and identified a candidate gene, cytochrome c oxidase assembly protein 15 (COX15). We report that biallelic COX15 pathogenic variants cause human oocyte ferroptosis and female infertility in a recessive inheritance pattern. COX15 variants impaired mitochondrial respiration in Saccharomyces cerevisiae and led to reduced protein levels in HeLa cells. Oocyte-specific deletion of Cox15 led to impaired Fe2+ and reactive oxygen species homeostasis that caused mitochondrial dysfunction and ultimately sensitized oocytes to ferroptosis. In addition, ferrostatin-1 (an inhibitor of ferroptosis) could rescue the oocyte ferroptosis phenotype in vitro and ex vivo. Our findings not only provide a genetic diagnostic marker for oocyte development defects but also expand the spectrum of mitochondrial disorders to female infertility and contribute to unique insights into the role of ferroptosis in human oocyte defects.

Gene-based burden tests of rare germline variants identify six cancer susceptibility genes.

Discovery of cancer risk variants in the sequence of the germline genome can shed light on carcinogenesis. Here we describe gene burden association analyses, aggregating rare missense and loss of function variants, at 22 cancer sites, including 130,991 cancer cases and 733,486 controls from Iceland, Norway and the United Kingdom. We identified four genes associated with increased cancer risk; the pro-apoptotic BIK for prostate cancer, the autophagy involved ATG12 for colorectal cancer, TG for thyroid cancer and CMTR2 for both lung cancer and cutaneous melanoma. Further, we found genes with rare variants that associate with decreased risk of cancer; AURKB for any cancer, irrespective of site, and PPP1R15A for breast cancer, suggesting that inhibition of PPP1R15A may be a preventive strategy for breast cancer. Our findings pinpoint several new cancer risk genes and emphasize autophagy, apoptosis and cell stress response as a focus point for developing new therapeutics.

6846 Rare Variants in SEMA3A in Individuals with Constitutional Delay of Puberty (CDP)

Abstract Disclosure: W. He: None. E.C. Schafer: None. R. Brauner: None. A. Delaney: None. L. Dunkel: None. R. Grinspon: None. J.E. Hall: None. J.N. Hirschhorn: None. S. Howard: None. A. Latronico: None. A.A. Jorge: None. K. McElreavey: None. V. Mericq: None. P.M. Merino: None. M.R. Palmert: None. R.A. Rey: None. R.C. Rezende: None. S.B. Seminara: None. N. Shaw: None. T. Delayed Puberty Genetics Consortium: None. Y. Chan: None. Background: Constitutional delay of puberty (CDP), also known as self-limited delayed puberty, refers to late onset of puberty with relatively normal progression afterward. In contrast, idiopathic hypogonadotropic hypogonadism (IHH) results in a long-lasting defect in reproductive endocrine activity. Both CDP and IHH are highly heritable and are sometimes observed in different family members in the same pedigree. Previous work by our group and others has suggested genetic overlap between the two conditions. In this study, we aimed to identify specific IHH-related genes that are also genetic causes of CDP. Methods: The Delayed Puberty Consortium assembled exome data for 233 individuals of diverse ancestries with CDP, with onset of puberty between 12 and 18 y in girls and between 13.5 and 18 y in boys. For control data, we used exome sequencing data from gnomAD v2.1.1 (n=125,748). After data harmonization and quality control, we used a modification of the “TRAPD” method to perform gene-based burden testing for 63 genes associated with IHH. We analyzed rare variants (minor allele frequency < 0.01) with high impact (premature stop, frameshift, start-lost, stop-lost) in all individuals across ancestries (n=233), as allele frequencies for high-impact variants are unlikely to vary widely across populations. We also analyzed rare variants with moderate impact (missense, inframe insertions/deletions) predicted to be deleterious by the Primate AI algorithm (score >0.6) in non-Finnish European (NFE) individuals with CDP (n=99), then meta-analyzed the results for high- and moderate-impact variants. We considered p-values < 7.9 x 10-4 (0.05 with Bonferroni adjustment for 63 genes) significant. Results: Of the 63 genes analyzed, SEMA3A (which encodes semaphorin-3A) exhibited significant enrichment for moderate-impact variants in the NFE CDP cohort compared to controls (variant frequency in CDP vs control: 0.064 vs 0.0076, p = 1.3 x 10-4) and in the meta-analysis of both high- and moderate-impact variants (p = 3.5 x 10-6). Novel SEMA3A variants identified exclusively in CDP individuals but not in controls were a frameshift variant, p.F546fsX2, and a missense variant, I252R. Also, three missense variants in SEMA3A were observed at a higher frequency in NFE CDP individuals compared with controls: p.R66T, p.V461A, and p.T717I. No other IHH genes showed significant enrichment of rare variants in the CDP cohort. Conclusion: Rare variants in SEMA3A contribute to both CDP and IHH. While functional studies are required to assess the effects of these specific variants on protein function, previous in vitro studies of similar SEMA3A variants demonstrated defects in secretion or signaling activity. Subsequent studies employing individual-level ancestry-matched controls and larger CDP cohorts are planned to identify additional genes that contribute to CDP. Presentation: 6/1/2024

12330 Rare Variants in SEMA3A in Individuals with Constitutional Delay of Puberty (CDP)

Abstract Disclosure: W. He: None. E.C. Schafer: None. R. Brauner: None. A. Delaney: None. L. Dunkel: None. R. Grinspon: None. J.E. Hall: None. J.N. Hirschhorn: None. S. Howard: None. A. Latronico: None. A.A. Jorge: None. K. McElreavey: None. V. Mericq: None. P.M. Merino: None. M.R. Palmert: None. R.A. Rey: None. R.C. Rezende: None. S.B. Seminara: None. N. Shaw: None. T. Delayed Puberty Genetics Consortium: None. Y. Chan: None. Background: Constitutional delay of puberty (CDP), also known as self-limited delayed puberty, refers to late onset of puberty with relatively normal progression afterward. In contrast, idiopathic hypogonadotropic hypogonadism (IHH) results in a long-lasting defect in reproductive endocrine activity. Both CDP and IHH are highly heritable and are sometimes observed in different family members in the same pedigree. Previous work by our group and others has suggested genetic overlap between the two conditions. In this study, we aimed to identify specific IHH-related genes that are also genetic causes of CDP. Methods: The Delayed Puberty Consortium assembled exome data for 233 individuals of diverse ancestries with CDP, with onset of puberty between 12 and 18 y in girls and between 13.5 and 18 y in boys. For control data, we used exome sequencing data from gnomAD v2.1.1 (n=125,748). After data harmonization and quality control, we used a modification of the “TRAPD” method to perform gene-based burden testing for 63 genes associated with IHH. We analyzed rare variants (minor allele frequency < 0.01) with high impact (premature stop, frameshift, start-lost, stop-lost) in all individuals across ancestries (n=233), as allele frequencies for high-impact variants are unlikely to vary widely across populations. We also analyzed rare variants with moderate impact (missense, inframe insertions/deletions) predicted to be deleterious by the Primate AI algorithm (score >0.6) in non-Finnish European (NFE) individuals with CDP (n=99), then meta-analyzed the results for high- and moderate-impact variants. We considered p-values < 7.9 x 10-4 (0.05 with Bonferroni adjustment for 63 genes) significant. Results: Of the 63 genes analyzed, SEMA3A (which encodes semaphorin-3A) exhibited significant enrichment for moderate-impact variants in the NFE CDP cohort compared to controls (variant frequency in CDP vs control: 0.064 vs 0.0076, p = 1.3 x 10-4) and in the meta-analysis of both high- and moderate-impact variants (p = 3.5 x 10-6). Novel SEMA3A variants identified exclusively in CDP individuals but not in controls were a frameshift variant, p.F546fsX2, and a missense variant, I252R. Also, three missense variants in SEMA3A were observed at a higher frequency in NFE CDP individuals compared with controls: p.R66T, p.V461A, and p.T717I. No other IHH genes showed significant enrichment of rare variants in the CDP cohort. Conclusion: Rare variants in SEMA3A contribute to both CDP and IHH. While functional studies are required to assess the effects of these specific variants on protein function, previous in vitro studies of similar SEMA3A variants demonstrated defects in secretion or signaling activity. Subsequent studies employing individual-level ancestry-matched controls and larger CDP cohorts are planned to identify additional genes that contribute to CDP. Presentation: 6/1/2024

Joint testing of rare variant burden scores using non-negative least squares

Gene-based burden tests are a popular and powerful approach for analysis of exome-wide association studies. These approaches combine sets of variants within a gene into a single burden score that is then tested for association. Typically, a range of burden scores are calculated and tested across a range of annotation classes and frequency bins. Correlation between these tests can complicate the multiple testing correction and hamper interpretation of the results. We introduce a method called the sparse burden association test (SBAT) that tests the joint set of burden scores under the assumption that causal burden scores act in the same effect direction. The method simultaneously assesses the significance of the model fit and selects the set of burden scores that best explain the association at the same time. Using simulated data, we show that the method is well calibrated and highlight scenarios where the test outperforms existing gene-based tests. We apply the method to 73 quantitative traits from the UK Biobank, showing that SBAT is a valuable additional gene-based test when combined with other existing approaches. This test is implemented in the REGENIE software.

Genome scan reveals several loci associated with torus palatinus.

Torus Palatinus (TP) is a common trait with an unclear aetiology. Although prior studies suggest a hereditary component, the genetic factors that influence TP risk remain unknown. The purpose of this study is to identify genetic variants associated with TP. We assessed the TP status of 829 individuals from various ancestral backgrounds using 3D palate scans. We then carried out a genome-wide association study (GWAS) to identify common variants associated with TP. We also performed gene-based tests across the exome to investigate the role of low-frequency coding variants. Our GWAS did not identify any genome-wide significant signals but identified suggestive associations including hits on chromosomes 2, 5 and 17 with p-values less than 5 × 10-6. Candidate genes at these suggestive loci have been implicated in normal-range craniofacial features, syndromes with facial and oral malformations, and bone density. We did not find evidence that low-frequency coding variants influence TP risk. In addition, we failed to replicate associations identified in prior genetic studies of TP. These findings suggest that multiple genes likely influence the development of TP. Independent replication will be required to confirm our suggestive associations.

Lack of a genetic risk continuum between pubertal timing in the general population and idiopathic hypogonadotropic hypogonadism.

Pubertal timing is a highly heritable trait in the general population. Recently, a large-scale exome-wide association study has implicated rare variants in six genes (KDM4C, MC3R, MKRN3, PDE10A, TACR3, and ZNF483) as genetic determinants of pubertal timing within the general population. Two of the genes (TACR3, MKRN3) are already implicated in extreme disorders of pubertal timing. This observation suggests that there may be a pervasive "genetic risk continuum" wherein genes that govern pubertal timing in the general population, by extension, may also be causal for rare Mendelian disorders of pubertal timing. Hence, we hypothesized that the four novel genes linked to pubertal timing in the population will also contribute to idiopathic hypogonadotropic hypogonadism (IHH), a genetic disorder characterized by absent puberty. Exome sequencing data from 1322 unrelated IHH probands were reviewed for rare sequence variants (RSVs) (minor allele frequency bins: <1%; <0.1%; <0.01%) in the six genes linked to puberty in the general population. A gene-based rare variant association testing (RVAT) was performed between the IHH cohort and a reference public genomic sequences repository-the Genome Aggregation Database (gnomAD). As expected, RVAT analysis showed that RSVs in TACR3, a known IHH gene, were significantly enriched in the IHH cohort compared to gnomAD cohort across all three MAF bins. However, RVAT analysis of the remaining five genes failed to show any RSV enrichment in the IHH cohort across all MAF bins. Our findings argue strongly against a pervasive genetic risk continuum between pubertal timing in the general population and extreme pubertal phenotypes. The biologic basis of such distinct genetic architectures' merits further evaluation.

Association analysis of mitochondrial DNA heteroplasmic variants: Methods and application

We rigorously assessed a comprehensive association testing framework for heteroplasmy, employing both simulated and real-world data. This framework employed a variant allele fraction (VAF) threshold and harnessed multiple gene-based tests for robust identification and association testing of heteroplasmy. Our simulation studies demonstrated that gene-based tests maintained an appropriate type I error rate at α = 0.001. Notably, when 5 % or more heteroplasmic variants within a target region were linked to an outcome, burden-extension tests (including the adaptive burden test, variable threshold burden test, and z-score weighting burden test) outperformed the sequence kernel association test (SKAT) and the original burden test. Applying this framework, we conducted association analyses on whole-blood derived heteroplasmy in 17,507 individuals of African and European ancestries (31 % of African Ancestry, mean age of 62, with 58 % women) with whole genome sequencing data. We performed both cohort- and ancestry-specific association analyses, followed by meta-analysis on both pooled samples and within each ancestry group. Our results suggest that mtDNA-encoded genes/regions are likely to exhibit varying rates in somatic aging, with the notably strong associations observed between heteroplasmy in the RNR1 and RNR2 genes (p < 0.001) and advance aging by the Original Burden test. In contrast, SKAT identified significant associations (p < 0.001) between diabetes and the aggregated effects of heteroplasmy in several protein-coding genes. Further research is warranted to validate these findings. In summary, our proposed statistical framework represents a valuable tool for facilitating association testing of heteroplasmy with disease traits in large human populations.

A New Method for Conditional Gene-Based Analysis Effectively Accounts for the Regional Polygenic Background.

Gene-based association analysis is a powerful tool for identifying genes that explain trait variability. An essential step of this analysis is a conditional analysis. It aims to eliminate the influence of SNPs outside the gene, which are in linkage disequilibrium with intragenic SNPs. The popular conditional analysis method, GCTA-COJO, accounts for the influence of several top independently associated SNPs outside the gene, correcting the z statistics for intragenic SNPs. We suggest a new TauCOR method for conditional gene-based analysis using summary statistics. This method accounts the influence of the full regional polygenic background, correcting the genotype correlations between intragenic SNPs. As a result, the distribution of z statistics for intragenic SNPs becomes conditionally independent of distribution for extragenic SNPs. TauCOR is compatible with any gene-based association test. TauCOR was tested on summary statistics simulated under different scenarios and on real summary statistics for a 'gold standard' gene list from the Open Targets Genetics project. TauCOR proved to be effective in all modelling scenarios and on real data. The TauCOR's strategy showed comparable sensitivity and higher specificity and accuracy than GCTA-COJO on both simulated and real data. The method can be successfully used to improve the effectiveness of gene-based association analyses.

Pheno- and Genotypic Epidemiological Characterization of Vancomycin-Resistant Enterococcus faecium Isolates from Intensive Care Unit Patients in Central Türkiye.

Vancomycin-resistant Enterococcus faecium (VRE) has been detected in Türkiye. Only limited information is available on its dissemination in the central regions of the country. This study describes the first epidemiological characterization of VRE clinical isolates detected in patients in a hospital in the province of Aksaray. In this one-year study conducted between 2021 and 2022, stool samples from intensive care unit patients were screened for VRE using the phenotypic E-test method, and the antibiotic sensitivity test was analyzed by using the VITEK® 2 system. A molecular assay for confirmation of species level was carried out by 16S rRNA gene-based sequencing and testing for antibiotic resistance (vanA or vanB) and virulence factor-encoding genes (esp, asa1, and hyl). Further, genotypic characterization was determined by macro-restriction fragment pattern analysis (MRFPA) of genomic DNA digested with SmaI restriction enzyme. Of the total 350 Enterococcus positive patients from different hospital intensive care units, 22 (6.3%) were positive for VRE using the phenotypic E-test method. All isolates showed resistance to ampicillin, ciprofloxacin, vancomycin, and teicoplanin and positive amplification for the vanA gene. However, none of the isolates was positive for the vanB gene. The most prevalent virulence gene was esp. The results indicate that the isolates are persistent in the hospital environment and subsequently transmitted to hospitalized patients, thus representing challenges to an outbreak and infection control. These study results would also help formulate more effective strategies to reduce the transmission and propagation of VRE contamination in various hospital settings.

Rare coding variant analysis for human diseases across biobanks and ancestries.

Large-scale sequencing has enabled unparalleled opportunities to investigate the role of rare coding variation in human phenotypic variability. Here, we present a pan-ancestry analysis of sequencing data from three large biobanks, including the All of Us research program. Using mixed-effects models, we performed gene-based rare variant testing for 601 diseases across 748,879 individuals, including 155,236 with ancestry dissimilar to European. We identified 363 significant associations, which highlighted core genes for the human disease phenome and identified potential novel associations, including UBR3 for cardiometabolic disease and YLPM1 for psychiatric disease. Pan-ancestry burden testing represented an inclusive and useful approach for discovery in diverse datasets, although we also highlight the importance of ancestry-specific sensitivity analyses in this setting. Finally, we found that effect sizes for rare protein-disrupting variants were concordant between samples similar to European ancestry and other genetic ancestries (βDeming = 0.7-1.0). Our results have implications for multi-ancestry and cross-biobank approaches in sequencing association studies for human disease.

The performance of AlphaMissense to identify genes influencing disease

A novel algorithm, AlphaMissense, has been shown to have an improved ability to predict the pathogenicity of rare missense genetic variants. However, it is not known whether AlphaMissense improves the ability of gene-based testing to identify disease-influencing genes. Using whole-exome sequencing data from the UK Biobank, we compared gene-based association analysis strategies including sets of deleterious variants: predicted loss-of-function (pLoF) variants only, pLoF plus AlphaMissense pathogenic variants, pLoF with missense variants predicted to be deleterious by any of five commonly utilized annotation methods (Missense (1/5)) or only variants predicted to be deleterious by all five methods (Missense (5/5)). We measured performance to identify 519 previously identified positive control genes, which can lead to Mendelian diseases, or are the targets of successfully developed medicines. These strategies identified 850k pLoF variants and 5 million deleterious missense variants, including 22k likely pathogenic missense variants identified exclusively by AlphaMissense. The gene-based association tests found 608 significant gene associations (at P<1.25x10-7) across 24 common traits and diseases. Compared to pLoFs plus Missense (5/5), tests using pLoFs and AlphaMissense variants found slightly more significant gene-disease and gene-trait associations, albeit with a marginally lower proportion of positive control genes. Nevertheless, their overall performance was similar. Merging AlphaMissense with Missense (5/5), whether through their intersection or union, did not yield any further enhancement in performance. In summary, employing AlphaMissense to select deleterious variants for gene-based testing did not improve the ability to identify genes that are known to influence disease.

A genome-wide association study identifies novel loci of vertigo in an Asian population-based cohort

The contributing genetic factors of vertigo remain poorly characterized, particularly in individuals of non-European ancestries. Here we show the genetic landscape of vertigo in an Asian population-based cohort. In a two-stage genome-wide association study (Ncase = 6199; Ncontrol = 54,587), we identify vertigo-associated genomic loci in DROSHA and ZNF91/LINC01224, with the latter replicating the findings in European ancestries. Gene-based association testing corroborates these findings. Interestingly, both genes are enriched in cerebellum, a key structure receiving sensory input from the vestibular system. Subjects carrying risk alleles from lead SNPs of DROSHA and ZNF91 incur a 1.74-fold risk of vertigo than those without. Moreover, composite clinical-polygenic risk scores allow differentiation between patients and controls, yielding an area under receiver operating characteristic curve of 0.69. This study identified novel genomic loci for vertigo in an Asian population-based cohort, which may help identifying high risk subjects and provide mechanistic insight in understanding the pathogenesis of vertigo.

Using gene and gene-set association tests to identify lethal prostate cancer genes.

Recent advances in the detection and treatment of prostate cancer (PCa) have reduced morbidity and mortality from this common cancer. Despite these improvements, PCa remains the second leading cause of cancer death in men in the United States. Further understanding of the genetic underpinnings of lethal PCa is required to drive risk detection and prevention and ultimately reduce mortality. We therefore set out to identify germline variants associated with cases of lethal prostate cancer (LPCa). Using a two-stage study design, we compared whole-exome sequencing data of 550 LPCa patients to 488 healthy male controls. Men were classified as having LPCa based on medical record review. Candidate genes were identified using gene- and gene-set-based rare truncating variant association tests. Case-control burden testing through Firth's penalized logistic regression and case-gnomAD allelic burden testing through a one-sided mid-p Fisher's exact test were conducted. Each gene's p-values from these tests were combined into an omnibus p-value for candidate gene selection. In the subsequent validation stage, genes were assessed using the UK Biobank and Firth's penalized logistic regression for each ancestry, combined through meta-analysis. Gene-based rare variant association tests identified 12 genes nominally associated with LPCa. Rare-variant association tests identified a gene set with a significantly higher burden of truncating germline mutations in LPCa patients than controls. Combining gene- and gene-set test results, four nominally significant genes (PPP1R3A, TG, PPFIBP2, and BTN3A3) were selected as candidates. Subsequent validation using the UK Biobank found that PPP1R3A was significantly associated with LPCa risk (odds ratio 2.34, CI 1.20-4.59). Specifically, pGln662ArgfsTer7 was identified as the predominant variant in PPP1R3A among LPCa patients in our dataset. Both individual gene and gene-set analyses identified candidates associated with LPCa. The novel association of PPP1R3A and LPCa risk merits further investigation.

Loss-of-function variants in ITSN1 confer high risk of Parkinson’s disease

Parkinson’s disease (PD) is a debilitating neurodegenerative disorder and its rising global incidence highlights the need for the identification of modifiable risk factors. In a gene-based burden test of rare variants (8647 PD cases and 777,693 controls) we discovered a novel association between loss-of-function variants in ITSN1 and PD. This association was further supported with burden data from the Neurodegenerative Disease Knowledge Portal and the Accelerating Medicines Partnership Parkinson’s Disease Knowledge Platform. Our findings show that Rho GTPases and disruptions in synaptic vesicle transport may be involved in the pathogenesis of PD, pointing to the possibility of novel therapeutic approaches.

A genome-wide investigation into the underlying genetic architecture of personality traits and overlap with psychopathology.

Personality is influenced by both genetic and environmental factors and is associated with other psychiatric traits such as anxiety and depression. The 'big five' personality traits, which include neuroticism, extraversion, agreeableness, conscientiousness and openness, are a widely accepted and influential framework for understanding and describing human personality. Of the big five personality traits, neuroticism has most often been the focus of genetic studies and is linked to various mental illnesses, including depression, anxiety and schizophrenia. Our knowledge of the genetic architecture of the other four personality traits is more limited. Here, utilizing the Million Veteran Program cohort, we conducted a genome-wide association study in individuals of European and African ancestry. Adding other published data, we performed genome-wide association study meta-analysis for each of the five personality traits with sample sizes ranging from 237,390 to 682,688. We identified 208, 14, 3, 2 and 7 independent genome-wide significant loci associated with neuroticism, extraversion, agreeableness, conscientiousness and openness, respectively. These findings represent 62 novel loci for neuroticism, as well as the first genome-wide significant loci discovered for agreeableness. Gene-based association testing revealed 254 genes showing significant association with at least one of the five personality traits. Transcriptome-wide and proteome-wide analysis identified altered expression of genes and proteins such as CRHR1, SLC12A5, MAPT and STX4. Pathway enrichment and drug perturbation analyses identified complex biology underlying human personality traits. We also studied the inter-relationship of personality traits with 1,437 other traits in a phenome-wide genetic correlation analysis, identifying new associations. Mendelian randomization showed positive bidirectional effects between neuroticism and depression and anxiety, while a negative bidirectional effect was observed for agreeableness and these psychiatric traits. This study improves our comprehensive understanding of the genetic architecture underlying personality traits and their relationship to other complex human traits.

Rare Sequence Variation Underlying Suspected Familial Cerebral Small-Vessel Disease.

Cerebral small-vessel disease (cSVD) is the leading monogenic cause of stroke. Despite genetic screening in routine diagnosis, many cases remain without a known causative variant. Using a cohort with suspected familial cSVD and whole-genome sequencing, we screened for variants in genes associated with monogenic cSVD and searched for novel variants associated with the disease. Rare variants were identified in whole-genome sequencing data from the NBR (National Institute for Health Research BioResource Rare Disease) study. Pathogenic variants in known monogenic cSVD genes were identified. Gene-based burden tests and family analysis were performed to identify novel variants associated with familial cSVD. A total of 257 suspected cSVD cases (mean ± SD age, 56.2 ± 16.1 years), and 13 086 controls with other nonstroke diseases (5874 [44.9%] men) were studied. A total of 8.9% of the cases carried a variant in known cSVD genes. Excluding these known causes, 23.6% of unrelated subjects with cSVD carried predicted deleterious variants in the Genomics England gene panel, but no association was found with cSVD in burden tests. We identified potential associations with cSVD in noncoding genes, including RP4-568F9.3 (adjusted P = 7.1 × 10-25), RP3-466I7.1 (adjusted P = 8.9 × 10-16), and ZNF209P (adjusted P = 1.0 × 10-15), and matrisomal genes (adjusted P = 5.1 × 10-6), including FAM20C, INHA, LAMC1, and VWA5B2. Predicted deleterious variants in known cSVD genes were present in 23.6% of unrelated cases with cSVD, but none of the genes were associated with the disease. Rare variants in noncoding and matrisomal genes could potentially contribute to cSVD development. These genes could play a role in tissue development and brain endothelial cell function. However, further studies are needed to confirm their pathophysiological roles.

Genome wide association study of clinical duration and age at onset of sporadic CJD.

Human prion diseases are rare, transmissible and often rapidly progressive dementias. The most common type, sporadic Creutzfeldt-Jakob disease (sCJD), is highly variable in clinical duration and age at onset. Genetic determinants of late onset or slower progression might suggest new targets for research and therapeutics. We assembled and array genotyped sCJD cases diagnosed in life or at autopsy. Clinical duration (median:4, interquartile range (IQR):2.5-9 (months)) was available in 3,773 and age at onset (median:67, IQR:61-73 (years)) in 3,767 cases. Phenotypes were successfully transformed to approximate normal distributions allowing genome-wide analysis without statistical inflation. 53 SNPs achieved genome-wide significance for the clinical duration phenotype; all of which were located at chromosome 20 (top SNP rs1799990, pvalue = 3.45x10-36, beta = 0.34 for an additive model; rs1799990, pvalue = 9.92x10-67, beta = 0.84 for a heterozygous model). Fine mapping, conditional and expression analysis suggests that the well-known non-synonymous variant at codon 129 is the obvious outstanding genome-wide determinant of clinical duration. Pathway analysis and suggestive loci are described. No genome-wide significant SNP determinants of age at onset were found, but the HS6ST3 gene was significant (pvalue = 1.93 x 10-6) in a gene-based test. We found no evidence of genome-wide genetic correlation between case-control (disease risk factors) and case-only (determinants of phenotypes) studies. Relative to other common genetic variants, PRNP codon 129 is by far the outstanding modifier of CJD survival suggesting only modest or rare variant effects at other genetic loci.

Gene-vegetarianism interactions in calcium, estimated glomerular filtration rate, and testosterone identified in genome-wide analysis across 30 biomarkers.

We examined the associations of vegetarianism with metabolic biomarkers using traditional and genetic epidemiology. First, we addressed inconsistencies in self-reported vegetarianism among UK Biobank participants by utilizing data from two dietary surveys to find a cohort of strict European vegetarians (N = 2,312). Vegetarians were matched 1:4 with nonvegetarians for non-genetic association analyses, revealing significant effects of vegetarianism in 15 of 30 biomarkers. Cholesterol measures plus vitamin D were significantly lower in vegetarians, while triglycerides were higher. A genome-wide association study revealed no genome-wide significant (GWS; 5×10-8) associations with vegetarian behavior. We performed genome-wide gene-vegetarianism interaction analyses for the biomarkers, and detected a GWS interaction impacting calcium at rs72952628 (P = 4.47×10-8). rs72952628 is in MMAA, a B12 metabolic pathway gene; B12 has major deficiency potential in vegetarians. Gene-based interaction tests revealed two significant genes, RNF168 in testosterone (P = 1.45×10-6) and DOCK4 in estimated glomerular filtration rate (eGFR) (P = 6.76×10-7), which have previously been associated with testicular and renal traits, respectively. These nutrigenetic findings indicate genotype can modify the associations between vegetarianism and health outcomes.

Association between rare pathogenic variants in established cancer risk genes and the diagnosis of single and multiple common cancers: A UK Biobank study.

10505 Background: Much of our understanding of cancer risk attributable to rare pathogenic variants (RPV) is derived from family-based and case-only cohort studies, which may be limited by ascertainment bias. We sought to identify associations between RPVs and cancer diagnoses in a large population-based study. Methods: We conducted a case-control study using the UK Biobank. We performed gene-based aggregate testing to examine the relationship between RPVs (defined as pathogenic or likely pathogenic variants) in 96 genes implicated in cancer risk and the diagnosis of 11 common cancers using the combined optimal unified sequence kernel association test (SKAT-O). Odds ratios (OR) and 95% confidence intervals (CI) were then calculated using Firth logistic regression. Results were adjusted for age, sex and genetic ancestry. Bonferroni-adjusted p-values were used to determine statistical significance. Variant pathogenicity was determined by American College of Medical Genetics criteria. Results: We identified 24 genes in which RPVs were statistically significantly associated with at least one of 11 common cancers among 183,626 individuals. The presence of an RPV in one of these 24 genes was associated with increased odds of one cancer (OR 1.8; 95% CI: 1.7-2.0), and multiple cancers (OR 2.5; 95% CI: 2.2-2.9). Conclusions: This large population-based study identified 24 established cancer risk genes in which RPVs were associated with common cancers. We identified both known and novel associations between cancer and RPVs. Some associations—i.e. renal/pancreas and MEN1—may be due to ascertainment and/or misclassification bias. We also show that individuals with RPVs in cancer risk genes not only have higher odds of one cancer diagnosis, but also multiple cancers.[Table: see text]