Professor Rodney Scott

Background/Objectives: Copper (Cu) and zinc (Zn) are essential trace elements, and an imbalance in their levels may influence the progression of cancer. The role of Cu and Zn levels in blood and serum, as well as 10-year survival rates in kidney cancer patients, remains unclear. Our objective was to determine the association between these micronutrients and mortality of kidney cancer patients. In this prospective study, we examined 284 consecutive, unselected kidney cancer patients and assessed their 10-year survival in relation to Cu and Zn levels. Methods: Micronutrient levels were measured using an inductively coupled plasma mass spectrometer. Each patient was categorized into one of four groups based on the distribution of Cu and Zn levels, ranked in increasing order. The multivariable models included factors such as age at diagnosis, gender, smoking history, type of surgery, and histopathological results. Results: We observed a significantly higher risk of all-cause mortality in patients with the highest blood or serum copper levels compared to those with the lower levels (blood: HR = 4.89; p < 0.001; serum: HR = 3.75; p < 0.001). With regard to zinc, we found a trend where lower blood or serum zinc levels (I quartile) were associated with higher mortality. Additionally, we identified a significant correlation between the Zn/Cu ratio and mortality. Conclusions: Patients in the lowest Zn/Cu ratio quartile had elevated hazard ratios compared to those in the higher quartile with HRs of 3.05 (p < 0.002) in blood and 5.72 (p < 0.001) in serum. To our knowledge, this study is the first to investigate the relationship between blood and serum levels of copper and zinc and kidney cancer survival.

Background: Implementing pharmacogenomic-guided management in cancer patients equitably and effectively in a large population presents challenges. DPYD genotyping determines clinically significant variants of patients at increased risk of developing grade3¿5 fluoropyrimidine (FP) toxicity. FP chemotherapies are prescribed for ~16,000 Australians with a 10%¿40% grade3¿4 toxicity incidence and 1% mortality. Variant carriers can have FP dosing adjusted to improve treatment tolerance without compromising anticancer effect. This strategy has not been formally adopted within Australia, despite widespread international standardisation. Aim: This pilot study determined genotyping turnaround-times (TAT) for 4 DPYD variants (c.1905+1G>A, c.1679T>G, c.2846A>T and c.1236G>A/Haplotype B3) in Australian patients. Secondary objectives were identification of FP toxicities of DPYD variant carriers, and analysis of healthcare stakeholder perspectives, including enablers/barriers to implementation. Methods: Genotyping was determined by Real-Time Polymerase Chain Reaction. Qualitative data were determined through semi-structured questionnaire. Results: 104 patients recruited over 24 months had a mean TAT of 7.2 days, 5.2 business days (range 1¿30). Grade3¿4 toxicity occurred in 9/16 DPYD variant carriers, including 2 ICU admissions and 1 death. Themes from 30 questionnaire respondents suggest that clinical environment and resources were fundamental barriers, and motivation to improve patient care was the predominant enabler of change. Conclusion: DPYD genotyping is feasible for improving precision-oncology for patients requiring FP chemotherapies. A TAT of 7 days is acceptable by both stakeholder respondents and national oncology clinician groups. This pilot study, although small, informs a large national project evaluating prospective DPYD genotyping and its impact on FP tolerability, patient safety and cost-effectiveness in Australia.

The COVID-19 pandemic, spanning from early 2020 to late 2022, posed unprecedented challenges for global public health. However, it also spurred innovative approaches to pandemic management, notably the development of pathogen detection in wastewater. It was successfully demonstrated that wastewater analysis can not only reflect ongoing COVID-19 infections but also serve as an early indicator of disease prevalence within communities. Recognizing the value of longitudinal analyses of various pathogens, we identified the need for wastewater biobanking. This practice allows for the retrospective analysis of samples, offering critical public health insights at the population level. Moreover, the potential to transport and store biobanked samples at ambient temperature or in a dry state could greatly enhance the utility of this technology, especially in resource-limited settings such as low- and middle-income countries. This article also addresses the ethical considerations and public health implications of wastewater-based epidemiology. While this approach holds significant potential beyond pathogen detection, it is essential to evaluate the benefits and potential risks carefully.

Background: The Swedish Family-Cancer Database (FCD) is the largest source of data on familial cancer in the world, including practically complete family structures and individual cancer diagnoses from the high-quality cancer registry. We present a novel application of FCD by analyzing age-specific familial risks and interpreting them through likely causes, such as germline pathogenic variants and/or environmental exposures. Main body: The basic assumption for this approach is that a discrete familial clustering in a narrow age-interval is not random but may provide causal clues. For this analysis we selected reasonably common cancers to meaningfully scrutinize familial risk through adulthood in which cancers are diagnosed, that included colorectal (CRC) and endometrial cancers, prostate and kidney cancers and breast and lung cancers. The interpretation is based on the literature. The highest familial relative risks for CRC and endometrial cancers were found at ages 40¿44¿years, matching the peak impact of mismatch repair gene mutations. However endometrial cancer showed also a small early onset component which could not be explained. Age-related familial risks for breast, prostate and kidney cancers also matched data from large-scale sequencing; these included the early onset component in kidney cancer which was likely due to VHL mutations. Age distribution of familial lung cancer was unique in showing a wide peak extending from middle to old ages, which would be consistent with a combination of direct genetic effects and indirect influence on inheritance of smoking dependence. Conclusions: The present review of age-specific familial risks and age-of-onset data from the literature may allow an interpretation that the familial and germline landscapes are reasonably harmonious for relatively early onset cancers but at higher ages no discrete peaks can be found which may implicate attenuated impact of high-risk genes and polygenic influence.

PURPOSE It is uncertain whether, and to what extent, hormonal contraceptives increase breast cancer (BC) risk for germline BRCA1 or BRCA2 mutation carriers.METHODS Using pooled observational data from four prospective cohort studies, associations between hormonal contraceptive use and BC risk for unaffected female BRCA1 and BRCA2 mutation carriers were assessed using Cox regression.RESULTSOf 3,882 BRCA1 and 1,509 BRCA2 mutation carriers, 53% and 71%, respectively, had ever used hormonal contraceptives for at least 1 year (median cumulative duration of use, 4.8 and 5.7 years, respectively). Overall, 488 BRCA1 and 191 BRCA2 mutation carriers developed BC during median follow-up of 5.9 and 5.6 years, respectively. Although for BRCA1 mutation carriers, neither current nor past use of hormonal contraceptives for at least 1 year was statistically significantly associated with BC risk (hazard ratio [HR], 1.40 [95% CI, 0.94 to 2.08], P =.10 for current use; 1.16 [0.80 to 1.69], P =.4, 1.40 [0.99 to 1.97], P =.05, and 1.27 [0.98 to 1.63], P =.07 for past use 1-5, 6-10, and >10 years before, respectively), ever use was associated with increased risk (HR, 1.29 [95% CI, 1.04 to 1.60], P =.02). Furthermore, BC risk increased with longer cumulative duration of use, with an estimated proportional increase in risk of 3% (1%-5%, P =.002) for each additional year of use. For BRCA2 mutation carriers, there was no evidence that current or ever use was associated with increased BC risk (HR, 0.70 [95% CI, 0.33 to 1.47], P =.3 and 1.07 [0.73 to 1.57], P =.7, respectively).CONCLUSION Hormonal contraceptives were associated with increased BC risk for BRCA1 mutation carriers, especially if used for longer durations. Decisions about their use in women with BRCA1 mutations should carefully weigh the risks and benefits for each individual.

Structural neuroimaging studies have identified a combination of shared and disorder-specific patterns of gray matter (GM) deficits across psychiatric disorders. Pooling large data allows for examination of a possible common neuroanatomical basis that may identify a certain vulnerability for mental illness. Large-scale collaborative research is already facilitated by data repositories, institutionally supported databases, and data archives. However, these data-sharing methodologies can suffer from significant barriers. Federated approaches augment these approaches by enabling access or more sophisticated, shareable and scaled-up analyses of large-scale data. We examined GM alterations using Collaborative Informatics and Neuroimaging Suite Toolkit for Anonymous Computation, an open-source, decentralized analysis application. Through federated analysis of eight sites, we identified significant overlap in the GM patterns (n = 4,102) of individuals with schizophrenia, major depressive disorder, and autism spectrum disorder. These results show cortical and subcortical regions that may indicate a shared vulnerability to psychiatric disorders.

Stratified medicine holds great promise to tailor treatment to the needs of individual patients. While genetics holds great potential to aid¿patient stratification, it remains a major challenge to operationalize complex genetic risk factor profiles to deconstruct clinical heterogeneity. Contemporary approaches to this problem rely on polygenic risk scores (PRS), which provide only limited clinical utility and lack a clear biological foundation. To overcome these limitations, we develop the CASTom-iGEx approach to stratify individuals based on the aggregated impact of their genetic risk factor profiles on tissue specific gene expression levels. The paradigmatic application of this approach to coronary artery disease or schizophrenia patient cohorts identified diverse strata or biotypes. These biotypes are characterized by distinct endophenotype profiles as well as clinical parameters and are fundamentally distinct from PRS based groupings. In stark contrast to the latter, the CASTom-iGEx strategy discovers biologically meaningful and clinically actionable patient subgroups, where complex genetic liabilities are not randomly distributed across individuals but rather converge onto distinct disease relevant biological processes. These results support the notion of different patient biotypes characterized by partially distinct pathomechanisms. Thus, the universally applicable approach presented here has the potential to constitute an important component of future personalized medicine paradigms.

Rare coding variants that substantially affect function provide insights into the biology of a gene1¿3. However, ascertaining the frequency of such variants requires large sample sizes4¿8. Here we present a catalogue of human protein-coding variation, derived from exome sequencing of 983,578 individuals across diverse populations. In total, 23% of the Regeneron Genetics Center Million Exome (RGC-ME) data come from individuals of African, East Asian, Indigenous American, Middle Eastern and South Asian ancestry. The catalogue includes more than 10.4 million missense and 1.1 million predicted loss-of-function (pLOF) variants. We identify individuals with rare biallelic pLOF variants in 4,848 genes, 1,751 of which have not been previously reported. From precise quantitative estimates of selection against heterozygous loss of function (LOF), we identify 3,988 LOF-intolerant genes, including 86 that were previously assessed as tolerant and 1,153 that lack established disease annotation. We also define regions of missense depletion at high resolution. Notably, 1,482 genes have regions that are depleted of missense variants despite being tolerant of pLOF variants. Finally, we estimate that 3% of individuals have a clinically actionable genetic variant, and that 11,773 variants reported in ClinVar with unknown significance are likely to be deleterious cryptic splice sites. To facilitate variant interpretation and genetics-informed precision medicine, we make this resource of coding variation from the RGC-ME dataset publicly accessible through a variant allele frequency browser.

Introduction: Blood pressure (BP) is a highly heritable trait with over 2000 underlying genomic loci identified to date. Although the kidney plays a key role, little is known about specific cell types involved in the genetic regulation of BP. Methods: Here, we applied stratified linkage disequilibrium score (LDSC) regression to connect BP genome-wide association studies (GWAS) results to specific cell types of the mature human kidney. We used the largest single-stage BP genome-wide analysis to date, including up to 1,028,980 adults of European ancestry, and single-cell transcriptomic data from 14 mature human kidneys, with mean age of 41 years. Results: Our analyses prioritized myofibroblasts and endothelial cells, among the total of 33 annotated cell type, as specifically involved in BP regulation (P < 0.05/33, i.e., 0.001515). Enrichment of heritability for systolic BP (SBP) was observed in myofibroblast cells in mature human kidney cortex, and enrichment of heritability for diastolic BP (DBP) was observed in descending vasa recta and peritubular capillary endothelial cells as well as stromal myofibroblast cells. The new finding of myofibroblast, the significant cell type for both BP traits, was consistent in 8 replication efforts using 7 sets of independent data, including in human fetal kidney, in East-Asian (EAS) ancestry, using mouse single-cell RNA sequencing (scRNA-seq) data, and when using another prioritization method. Conclusion: Our findings provide a solid basis for follow-up studies to further identify genes and mechanisms in myofibroblast cells that underlie the regulation of BP.

Co-observation of a gene variant with a pathogenic variant in another gene that explains the disease presentation has been designated as evidence against pathogenicity for commonly used variant classification guidelines. Multiple variant curation expert panels have specified, from consensus opinion, that this evidence type is not applicable for the classification of breast cancer predisposition gene variants. Statistical analysis of sequence data for 55,815 individuals diagnosed with breast cancer from the BRIDGES sequencing project was undertaken to formally assess the utility of co-observation data for germline variant classification. Our analysis included expected loss-of-function variants in 11 breast cancer predisposition genes and pathogenic missense variants in BRCA1, BRCA2, and TP53. We assessed whether co-observation of pathogenic variants in two different genes occurred more or less often than expected under the assumption of independence. Co-observation of pathogenic variants in each of BRCA1, BRCA2, and PALB2 with the remaining genes was less frequent than expected. This evidence for depletion remained after adjustment for age at diagnosis, study design (familial versus population-based), and country. Co-observation of a variant of uncertain significance in BRCA1, BRCA2, or PALB2 with a pathogenic variant in another breast cancer gene equated to supporting evidence against pathogenicity following criterion strength assignment based on the likelihood ratio and showed utility in reclassification of missense BRCA1 and BRCA2 variants identified in BRIDGES. Our approach has applicability for assessing the value of co-observation as a predictor of variant pathogenicity in other clinical contexts, including for gene-specific guidelines developed by ClinGen Variant Curation Expert Panels.

Colorectal cancer (CRC) is the third most common cancer worldwide. Screening programs allow early diagnosis and have improved the clinical management of this disease. Aberrant DNA methylation is increasingly being explored as potential biomarkers for many types of cancers. In this study we investigate the methylation of ten target genes in 105 CRC and paired normal adjacent colonic tissue samples using a MethylLight droplet digital PCR (ML-ddPCR) assay. Receiver operator characteristic (ROC) curves were used to determine the diagnostic performance of all target genes individually and in combination. All 515 different combinations of genes showed significantly higher levels of methylation in CRC tissue. The combination of multiple target genes into a single test generally resulted in greater diagnostic accuracy when compared to single target genes. Our data confirms that ML-ddPCR is able to reliably detect significant differences in DNA methylation between CRC tissue and normal adjacent colonic tissue in a specific selection of target genes.

Evidence from literature, including the BRIDGES study, indicates that germline protein truncating variants (PTVs) in FANCM confer moderately increased risk of ER-negative and triple-negative breast cancer (TNBC), especially for women with a family history of the disease. Association between FANCM missense variants (MVs) and breast cancer risk has been postulated. In this study, we further used the BRIDGES study to test 689 FANCM MVs for association with breast cancer risk, overall and in ER-negative and TNBC subtypes, in 39,885 cases (7566 selected for family history) and 35,271 controls of European ancestry. Sixteen common MVs were tested individually; the remaining rare 673 MVs were tested by burden analyses considering their position and pathogenicity score. We also conducted a meta-analysis of our results and those from published studies. We did not find evidence for association for any of the 16 variants individually tested. The rare MVs were significantly associated with increased risk of ER-negative breast cancer by burden analysis comparing familial cases to controls (OR = 1.48; 95% CI 1.07¿2.04; P = 0.017). Higher ORs were found for the subgroup of MVs located in functional domains or predicted to be pathogenic. The meta-analysis indicated that FANCM MVs overall are associated with breast cancer risk (OR = 1.22; 95% CI 1.08¿1.38; P = 0.002). Our results support the definition from previous analyses of FANCM as a moderate-risk breast cancer gene and provide evidence that FANCM MVs could be low/moderate risk factors for ER-negative and TNBC subtypes. Further genetic and functional analyses are necessary to clarify better the increased risks due to FANCM MVs.

Background: Breast cancer (BC) patients with a germline CHEK2 c.1100delC variant have an increased risk of contralateral BC (CBC) and worse BC-specific survival (BCSS) compared to non-carriers. Aim: To assessed the associations of CHEK2 c.1100delC, radiotherapy, and systemic treatment with CBC risk and BCSS. Methods: Analyses were based on 82,701 women diagnosed with a first primary invasive BC including 963 CHEK2 c.1100delC carriers; median follow-up was 9.1 years. Differential associations with treatment by CHEK2 c.1100delC status were tested by including interaction terms in a multivariable Cox regression model. A multi-state model was used for further insight into the relation between CHEK2 c.1100delC status, treatment, CBC risk and death. Results: There was no evidence for differential associations of therapy with CBC risk by CHEK2 c.1100delC status. The strongest association with reduced CBC risk was observed for the combination of chemotherapy and endocrine therapy [HR (95% CI): 0.66 (0.55¿0.78)]. No association was observed with radiotherapy. Results from the multi-state model showed shorter BCSS for CHEK2 c.1100delC carriers versus non-carriers also after accounting for CBC occurrence [HR (95% CI): 1.30 (1.09¿1.56)]. Conclusion: Systemic therapy was associated with reduced CBC risk irrespective of CHEK2 c.1100delC status. Moreover, CHEK2 c.1100delC carriers had shorter BCSS, which appears not to be fully explained by their CBC risk.

Evidence linking coding germline variants in breast cancer (BC)-susceptibility genes other than BRCA1, BRCA2, and CHEK2 with contralateral breast cancer (CBC) risk and breast cancer-specific survival (BCSS) is scarce. The aim of this study was to assess the association of protein-truncating variants (PTVs) and rare missense variants (MSVs) in nine known (ATM, BARD1, BRCA1, BRCA2, CHEK2, PALB2, RAD51C, RAD51D, and TP53) and 25 suspected BC-susceptibility genes with CBC risk and BCSS. Hazard ratios (HRs) and 95% confidence intervals (CIs) were estimated with Cox regression models. Analyses included 34,401 women of European ancestry diagnosed with BC, including 676 CBCs and 3,449 BC deaths; the median follow-up was 10.9 years. Subtype analyses were based on estrogen receptor (ER) status of the first BC. Combined PTVs and pathogenic/likely pathogenic MSVs in BRCA1, BRCA2, and TP53 and PTVs in CHEK2 and PALB2 were associated with increased CBC risk [HRs (95% CIs): 2.88 (1.70¿4.87), 2.31 (1.39¿3.85), 8.29 (2.53¿27.21), 2.25 (1.55¿3.27), and 2.67 (1.33¿5.35), respectively]. The strongest evidence of association with BCSS was for PTVs and pathogenic/likely pathogenic MSVs in BRCA2 (ER-positive BC) and TP53 and PTVs in CHEK2 [HRs (95% CIs): 1.53 (1.13¿2.07), 2.08 (0.95¿4.57), and 1.39 (1.13¿1.72), respectively, after adjusting for tumor characteristics and treatment]. HRs were essentially unchanged when censoring for CBC, suggesting that these associations are not completely explained by increased CBC risk, tumor characteristics, or treatment. There was limited evidence of associations of PTVs and/or rare MSVs with CBC risk or BCSS for the 25 suspected BC genes. The CBC findings are relevant to treatment decisions, follow-up, and screening after BC diagnosis.

To detect novel endometrial cancer risk variants, we leveraged information from endometrial cancer risk factors in a multi-trait GWAS analysis. We first assessed causal relationships between established and suspected endometrial cancer risk factors, and endometrial cancer using Mendelian randomization. Following multivariable analysis, five independent risk factors (waist circumference, testosterone levels, sex hormone binding globulin levels, age at menarche, and age at natural menopause) were included in a multi-trait Bayesian GWAS analysis. We identified three potentially novel loci that associate with endometrial cancer risk, one of which (7q22.1) replicated in an independent endometrial cancer GWAS dataset and was genome-wide significant in a meta-analysis. This locus may affect endometrial cancer risk through altered testosterone levels. Consistent with this, we observed colocalization between the signals for endometrial cancer risk and expression of CYP3A7, a gene involved in testosterone metabolism. Thus, our findings suggest opportunities for hormone therapy to prevent or treat endometrial cancer.

Genetics have nominated many schizophrenia risk genes and identified convergent signals between schizophrenia and neurodevelopmental disorders. However, functional interpretation of the nominated genes in the relevant brain cell types is often lacking. We executed interaction proteomics for six schizophrenia risk genes that have also been implicated in neurodevelopment in human induced cortical neurons. The resulting protein network is enriched for common variant risk of schizophrenia in Europeans and East Asians, is down-regulated in layer 5/6 cortical neurons of individuals affected by schizophrenia, and can complement fine-mapping and eQTL data to prioritize additional genes in GWAS loci. A sub-network centered on HCN1 is enriched for common variant risk and contains proteins (HCN4 and AKAP11) enriched for rare protein-truncating mutations in individuals with schizophrenia and bipolar disorder. Our findings showcase brain cell-type-specific interactomes as an organizing framework to facilitate interpretation of genetic and transcriptomic data in schizophrenia and its related disorders.

Background: Low-frequency variants play an important role in breast cancer (BC) susceptibility. Gene-based methods can increase power by combining multiple variants in the same gene and help identify target genes. Methods: We evaluated the potential of gene-based aggregation in the Breast Cancer Association Consortium cohorts including 83,471 cases and 59,199 controls. Low-frequency variants were aggregated for individual genes' coding and regulatory regions. Association results in European ancestry samples were compared to single-marker association results in the same cohort. Gene-based associations were also combined in meta-analysis across individuals with European, Asian, African, and Latin American and Hispanic ancestry. Results: In European ancestry samples, 14 genes were significantly associated (q < 0.05) with BC. Of those, two genes, FMNL3 (P = 6.11 × 10-6) and AC058822.1 (P = 1.47 × 10-4), represent new associations. High FMNL3 expression has previously been linked to poor prognosis in several other cancers. Meta-analysis of samples with diverse ancestry discovered further associations including established candidate genes ESR1 and CBLB. Furthermore, literature review and database query found further support for a biologically plausible link with cancer for genes CBLB, FMNL3, FGFR2, LSP1, MAP3K1, and SRGAP2C. Conclusions: Using extended gene-based aggregation tests including coding and regulatory variation, we report identification of plausible target genes for previously identified single-marker associations with BC as well as the discovery of novel genes implicated in BC development. Including multi ancestral cohorts in this study enabled the identification of otherwise missed disease associations as ESR1 (P = 1.31 × 10-5), demonstrating the importance of diversifying study cohorts.

Background: Pregnancy in women with multiple sclerosis (wwMS) is associated with a reduction of long-term disability progression. The mechanism that drives this effect is unknown, but converging evidence suggests a role for epigenetic mechanisms altering immune and/or central nervous system function. In this study, we aimed to identify whole blood and immune cell-specific DNA methylation patterns associated with parity in relapse-onset MS. Results: We investigated the association between whole blood and immune cell-type-specific genome-wide methylation patterns and parity in 192 women with relapse-onset MS, matched for age and disease severity. The median time from last pregnancy to blood collection was 16.7¿years (range = 1.5¿44.4¿years). We identified 2965 differentially methylated positions in whole blood, 68.5% of which were hypermethylated in parous women; together with two differentially methylated regions on Chromosomes 17 and 19 which mapped to TMC8 and ZNF577, respectively. Our findings validated 22 DMPs and 366 differentially methylated genes from existing literature on epigenetic changes associated with parity in wwMS. Differentially methylated genes in whole blood were enriched in neuronal structure and growth-related pathways. Immune cell-type-specific analysis using cell-type proportion estimates from statistical deconvolution of whole blood revealed further differential methylation in T cells specifically (four in CD4+ and eight in CD8+ T cells). We further identified reduced methylation age acceleration in parous women, demonstrating slower biological aging compared to nulligravida women. Conclusion: Differential methylation at genes related to neural plasticity offers a potential molecular mechanism driving the long-term effect of pregnancy on MS outcomes. Our results point to a potential 'CNS signature' of methylation in peripheral immune cells, as previously described in relation to MS progression, induced by parity. As the first epigenome-wide association study of parity in wwMS reported, validation studies are needed to confirm our findings.

Schizophrenia (SCZ) is a chronic mental illness and among the most debilitating conditions encountered in medical practice. A recent landmark SCZ study of the protein-coding regions of the genome identified a causal role for ten genes and a concentration of rare variant signals in evolutionarily constrained genes1. This recent study¿and most other large-scale human genetics studies¿was mainly composed of individuals of European (EUR) ancestry, and the generalizability of the findings in non-EUR populations remains unclear. To address this gap, we designed a custom sequencing panel of 161 genes selected based on the current knowledge of SCZ genetics and sequenced a new cohort of 11,580 SCZ cases and 10,555 controls of diverse ancestries. Replicating earlier work, we found that cases carried a significantly higher burden of rare protein-truncating variants (PTVs) among evolutionarily constrained genes (odds ratio = 1.48; P = 5.4 × 10-6). In meta-analyses with existing datasets totaling up to 35,828 cases and 107,877 controls, this excess burden was largely consistent across five ancestral populations. Two genes (SRRM2 and AKAP11) were newly implicated as SCZ risk genes, and one gene (PCLO) was identified as shared by individuals with SCZ and those with autism. Overall, our results lend robust support to the rare allelic spectrum of the genetic architecture of SCZ being conserved across diverse human populations.

While germline copy-number variants (CNVs) contribute to schizophrenia (SCZ) risk, the contribution of somatic CNVs (sCNVs)¿present in some but not all cells¿remains unknown. We identified sCNVs using blood-derived genotype arrays from 12,834 SCZ cases and 11,648 controls, filtering sCNVs at loci recurrently mutated in clonal blood disorders. Likely early-developmental sCNVs were more common in cases (0.91%) than controls (0.51%, p = 2.68e-4), with recurrent somatic deletions of exons 1¿5 of the NRXN1 gene in five SCZ cases. Hi-C maps revealed ectopic, allele-specific loops forming between a potential cryptic promoter and non-coding cis-regulatory elements upon 5' deletions in NRXN1. We also observed recurrent intragenic deletions of ABCB11, encoding a transporter implicated in anti-psychotic response, in five treatment-resistant SCZ cases and showed that ABCB11 is specifically enriched in neurons forming mesocortical and mesolimbic dopaminergic projections. Our results indicate potential roles of sCNVs in SCZ risk.

In vivo experimental analysis of human brain tissue poses substantial challenges and ethical concerns. To address this problem, we developed a computational method called the Brain Gene Expression and Network-Imputation Engine (BrainGENIE) that leverages peripheral-blood transcriptomes to predict brain tissue-specific gene-expression levels. Paired blood¿brain transcriptomic data collected by the Genotype-Tissue Expression (GTEx) Project was used to train BrainGENIE models to predict gene-expression levels in ten distinct brain regions using whole-blood gene-expression profiles. The performance of BrainGENIE was compared to PrediXcan, a popular method for imputing gene expression levels from genotypes. BrainGENIE significantly predicted brain tissue-specific expression levels for 2947¿11,816 genes (false-discovery rate-adjusted p < 0.05), including many transcripts that cannot be predicted significantly by a transcriptome-imputation method such as PrediXcan. BrainGENIE recapitulated measured diagnosis-related gene-expression changes in the brain for autism, bipolar disorder, and schizophrenia better than direct correlations from blood and predictions from PrediXcan. We developed a convenient software toolset for deploying BrainGENIE, and provide recommendations for how best to implement models. BrainGENIE complements and, in some ways, outperforms existing transcriptome-imputation tools, providing biologically meaningful predictions and opening new research avenues.

Purpose: Genetic variants causing aberrant premessenger RNA splicing are increasingly being recognized as causal variants in genetic disorders. In this study, we devise standardized practices for polymerase chain reaction (PCR)-based RNA diagnostics using clinically accessible specimens (blood, fibroblasts, urothelia, biopsy). Methods: A total of 74 families with diverse monogenic conditions (31% prenatal-congenital onset, 47% early childhood, and 22% teenage-adult onset) were triaged into PCR-based RNA testing, with comparative RNA sequencing for 19 cases. Results: Informative RNA assay data were obtained for 96% of cases, enabling variant reclassification for 75% variants that can be used for genetic counseling (71%), to inform clinical care (32%) and prenatal counseling (41%). Variant-associated mis-splicing was highly reproducible for 28 cases with samples from =2 affected individuals or heterozygotes and 10 cases with =2 biospecimens. PCR amplicons encompassing another segregated heterozygous variant was vital for clinical interpretation of 22 of 79 variants to phase RNA splicing events and discern complete from partial mis-splicing. Conclusion: RNA diagnostics enabled provision of a genetic diagnosis for 64% of recruited cases. PCR-based RNA diagnostics has capacity to analyze 81.3% of clinically significant genes, with long amplicons providing an advantage over RNA sequencing to phase RNA splicing events. The Australasian Consortium for RNA Diagnostics (SpliceACORD) provide clinically-endorsed, standardized protocols and recommendations for interpreting RNA assay data.

We conduct a genome-wide association study (GWAS) of educational attainment (EA) in a sample of ~3 million individuals and identify 3,952 approximately uncorrelated genome-wide-significant single-nucleotide polymorphisms (SNPs). A genome-wide polygenic predictor, or polygenic index (PGI), explains 12¿16% of EA variance and contributes to risk prediction for ten diseases. Direct effects (i.e., controlling for parental PGIs) explain roughly half the PGI's magnitude of association with EA and other phenotypes. The correlation between mate-pair PGIs is far too large to be consistent with phenotypic assortment alone, implying additional assortment on PGI-associated factors. In an additional GWAS of dominance deviations from the additive model, we identify no genome-wide-significant SNPs, and a separate X-chromosome additive GWAS identifies 57.

Background: Genome-wide association studies (GWAS) have identified multiple common breast cancer susceptibility variants. Many of these variants have differential associations by estrogen receptor (ER) status, but how these variants relate with other tumor features and intrinsic molecular subtypes is unclear. Methods: Among 106,571 invasive breast cancer cases and 95,762 controls of European ancestry with data on 173 breast cancer variants identified in previous GWAS, we used novel two-stage polytomous logistic regression models to evaluate variants in relation to multiple tumor features (ER, progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2) and grade) adjusting for each other, and to intrinsic-like subtypes. Results: Eighty-five of 173 variants were associated with at least one tumor feature (false discovery rate < 5%), most commonly ER and grade, followed by PR and HER2. Models for intrinsic-like subtypes found nearly all of these variants (83 of 85) associated at p < 0.05 with risk for at least one luminal-like subtype, and approximately half (41 of 85) of the variants were associated with risk of at least one non-luminal subtype, including 32 variants associated with triple-negative (TN) disease. Ten variants were associated with risk of all subtypes in different magnitude. Five variants were associated with risk of luminal A-like and TN subtypes in opposite directions. Conclusion: This report demonstrates a high level of complexity in the etiology heterogeneity of breast cancer susceptibility variants and can inform investigations of subtype-specific risk prediction.

Background: The variation in multiple sclerosis (MS) disease severity is incompletely explained by genetics, suggesting genetic and environmental interactions are involved. Moreover, the lack of prognostic biomarkers makes it difficult for clinicians to optimise care. DNA methylation is one epigenetic mechanism by which gene¿environment interactions can be assessed. Here, we aimed to identify DNA methylation patterns associated with mild and severe relapse-onset MS (RMS) and to test the utility of methylation as a predictive biomarker. Methods: We conducted an epigenome-wide association study between 235 females with mild (n = 119) or severe (n = 116) with RMS. Methylation was measured with the Illumina methylationEPIC array and analysed using logistic regression. To generate hypotheses about the functional consequence of differential methylation, we conducted gene set enrichment analysis using ToppGene. We compared the accuracy of three machine learning models in classifying disease severity: (1) clinical data available at baseline (age at onset and first symptoms) built using elastic net (EN) regression, (2) methylation data using EN regression and (3) a weighted methylation risk score of differentially methylated positions (DMPs) from the main analysis using logistic regression. We used a conservative 70:30 test:train split for classification modelling. A false discovery rate threshold of 0.05 was used to assess statistical significance. Results: Females with mild or severe RMS had 1472 DMPs in whole blood (839 hypermethylated, 633 hypomethylated in the severe group). Differential methylation was enriched in genes related to neuronal cellular compartments and processes, and B-cell receptor signalling. Whole-blood methylation levels at 1708 correlated CpG sites classified disease severity more accurately (machine learning model 2, AUC = 0.91) than clinical data (model 1, AUC = 0.74) or the wMRS (model 3, AUC = 0.77). Of the 1708 selected CpGs, 100 overlapped with DMPs from the main analysis at the gene level. These overlapping genes were enriched in neuron projection and dendrite extension, lending support to our finding that neuronal processes, rather than immune processes, are implicated in disease severity. Conclusion: RMS disease severity is associated with whole-blood methylation at genes related to neuronal structure and function. Moreover, correlated whole-blood methylation patterns can assign disease severity in females with RMS more accurately than clinical data available at diagnosis.

Schizophrenia has a heritability of 60¿80%1, much of which is attributable to common risk alleles. Here, in a two-stage genome-wide association study of up to 76,755 individuals with schizophrenia and 243,649 control individuals, we report common variant associations at 287 distinct genomic loci. Associations were concentrated in genes that are expressed in excitatory and inhibitory neurons of the central nervous system, but not in other tissues or cell types. Using fine-mapping and functional genomic data, we identify 120 genes (106 protein-coding) that are likely to underpin associations at some of these loci, including 16 genes with credible causal non-synonymous or untranslated region variation. We also implicate fundamental processes related to neuronal function, including synaptic organization, differentiation and transmission. Fine-mapped candidates were enriched for genes associated with rare disruptive coding variants in people with schizophrenia, including the glutamate receptor subunit GRIN2A and transcription factor SP4, and were also enriched for genes implicated by such variants in neurodevelopmental disorders. We identify biological processes relevant to schizophrenia pathophysiology; show convergence of common and rare variant associations in schizophrenia and neurodevelopmental disorders; and provide a resource of prioritized genes and variants to advance mechanistic studies.

Purpose: The association between the NOD2 c.3020insC allele and CDKN2A missense variant c.442G>A (p.P.A148T) and survival of patients with bladder or kidney cancer remains controversial. Materials and Methods: We compared the allele frequencies of NOD2 c.3020insC and CDKN2A p.A148T allele in 706 patients with bladder cancer, 410 cases with kidney cancer against two control groups. The Cox proportional hazards model was used to determine whether there were any survival differences between carriers of the NOD2 c.3020insC or the CDKN2A p.A148T variant. Results: Among the three patient subgroups: patients under 60 years of age, non-smokers and a third with histological features of low grade noninvasive papillary bladder cancer, we observed that the c.3020insC allele had a nominal statistically significant effect on survival. We also observed that the NOD2 c.3020insC variant was more frequent in patients with bladder cancer aged between 51 and 60 years. There was some nominal evidence that the CDKN2A p.A148T polymorphism reduced survival in the subgroup of bladder cancer patients under 60 years of age. We observed that in kidney cancer patients, the incidence of the NOD2 variant appeared to be lower in the group aged between 60 and 70 years, however, this was not statistically significant. In addition, in patients with histological features of grade III chromophobic kidney cancer, the c.3020insC allele also appeared to be overrepresented but this too was not statistically significant. Conclusion: We have shown that the NOD2 c.3020insC allele and the CDKN2A p.A148T polymorphism does not play a role in the survival of patients with bladder cancer.

Conventional genome-wide association studies (GWASs) of complex traits, such as Multiple Sclerosis (MS), are reliant on per-SNP p-values and are therefore heavily burdened by multiple testing correction. Thus, in order to detect more subtle alterations, ever increasing sample sizes are required, while ignoring potentially valuable information that is readily available in existing datasets. To overcome this, we used penalised regression incorporating elastic net with a stability selection method by iterative subsampling to detect the potential interaction of loci with MS risk. Through re-analysis of the ANZgene dataset (1617 cases and 1988 controls) and an IMSGC dataset as a replication cohort (1313 cases and 1458 controls), we identified new association signals for MS predisposition, including SNPs above and below conventional significance thresholds while targeting two natural killer receptor loci and the well-established HLA loci. For example, rs2844482 (98.1% iterations), otherwise ignored by conventional statistics (p = 0.673) in the same dataset, was independently strongly associated with MS in another GWAS that required more than 40 times the number of cases (~45 K). Further comparison of our hits to those present in a large-scale meta-analysis, confirmed that the majority of SNPs identified by the elastic net model reached conventional statistical GWAS thresholds (p < 5 × 10-8 ) in this much larger dataset. Moreover, we found that gene variants involved in oxidative stress, in addition to innate immunity, were associated with MS. Overall, this study highlights the benefit of using more advanced statistical methods to (re-)analyse subtle genetic variation among loci that have a biological basis for their contribution to disease risk.

The CAPP2 trial investigated the long-term effects of aspirin and resistant starch on cancer incidence in patients with Lynch syndrome (LS). Participants with LS were randomized double-blind to 30 g resistant starch (RS) daily or placebo for up to 4 years. We present long-term cancer outcomes based on the planned 10-year follow-up from recruitment, supplemented by National Cancer Registry data to 20 years in England, Wales, and Finland. Overall, 463 participants received RS and 455 participants received placebo. After up to 20 years follow-up, there was no difference in colorectal cancer incidence (n = 52 diagnosed with colorectal cancer among those randomized to RS against n = 53 on placebo) but fewer participants had non-colorectal LS cancers in those randomized to RS (n = 27) compared with placebo (n = 48); intention-to-treat (ITT) analysis [HR, 0.54; 95% confidence interval (CI), 0.33-0.86; P = 0.010]. In ITT analysis, allowing for multiple primary cancer diagnoses among participants by calculating incidence rate ratios (IRR) confirmed the protective effect of RS against non-colorectal cancer LS cancers (IRR, 0.52; 95% CI, 0.32-0.84; P = 0.0075). These effects are particularly pronounced for cancers of the upper GI tract; 5 diagnoses in those on RS versus 21 diagnoses on placebo. The reduction in non-colorectal cancer LS cancers was detectable in the first 10 years and continued in the next decade. For colorectal cancer, ITT analysis showed no effect of RS on colorectal cancer risk (HR, 0.92; 95% CI, 0.62-1.34; P = 0.63). There was no interaction between aspirin and RS treatments. In conclusion, 30 g daily RS appears to have a substantial protective effect against non-colorectal cancer cancers for patients with LS. Prevention Relevance: Regular bowel screening and aspirin reduce colorectal cancer among patients with LS but extracolonic cancers are difficult to detect and manage. This study suggests that RS reduces morbidity associated with extracolonic cancers.

Understanding the genomic basis of memory processes may help in combating neurodegenerative disorders. Hence, we examined the associations of common genetic variants with verbal short-term memory and verbal learning in adults without dementia or stroke (N = 53,637). We identified novel loci in the intronic region of CDH18, and at 13q21 and 3p21.1, as well as an expected signal in the APOE/APOC1/TOMM40 region. These results replicated in an independent sample. Functional and bioinformatic analyses supported many of these loci and further implicated POC1. We showed that polygenic score for verbal learning associated with brain activation in right parieto-occipital region during working memory task. Finally, we showed genetic correlations of these memory traits with several neurocognitive and health outcomes. Our findings suggest a role of several genomic loci in verbal memory processes.

Triple negative breast cancer (TNBC) is a highly aggressive subtype with a high rate of metastasis, early distant recurrence and resistance to therapy leading to worse survival than other breast cancer subtypes. There are no well-established biomarkers that can determine women who will do better and those who are likely to have poorer outcomes with TNBC, nor are there targeted therapies. Thus, the identification of prognostic and/or predictive biomarkers will enable tailored therapies based on their likelihood of disease outcomes and may prevent over- and under-diagnosis. Previous studies from our laboratory have identified four genes (ANP32E, DSC2, ANKRD30A and IL6ST/gp130) that are specific to TNBC and were associated with lymph node metastasis (LNmets), the earliest indicator of tumor progression via distal spread. This study aimed to validate these findings using absolute quantitation by digital droplet PCR (ddPCR) and to determine relationships with clinicopathological features and survival. Our analysis confirmed all four genes displayed significant expression differences between TNBC cases and non-TNBC cases. Moreover, low IL6ST expression was significantly associated with grade 3 disease, hormone receptor negativity and earlier age at diagnosis; low ANKRD30A expression was associated with tumor size; and high ANP32E expression was significantly associated with grade and the number of positive lymph nodes. Individually, three of the four genes were associated with relapse-free survival in TNBC and in combination, all four genes were significantly associated with TNBC survival, but not in hormone receptor-positive cases. Collectively our results suggest that the four genes may have utility in TNBC prognostication.

Germline copy number variants (CNVs) are pervasive in the human genome but potential disease associations with rare CNVs have not been comprehensively assessed in large datasets. We analysed rare CNVs in genes and non-coding regions for 86,788 breast cancer cases and 76,122 controls of European ancestry with genome-wide array data. Gene burden tests detected the strongest association for deletions in BRCA1 (P = 3.7E-18). Nine other genes were associated with a p-value < 0.01 including known susceptibility genes CHEK2 (P = 0.0008), ATM (P = 0.002) and BRCA2 (P = 0.008). Outside the known genes we detected associations with p-values < 0.001 for either overall or subtype-specific breast cancer at nine deletion regions and four duplication regions. Three of the deletion regions were in established common susceptibility loci. To the best of our knowledge, this is the first genome-wide analysis of rare CNVs in a large breast cancer case-control dataset. We detected associations with exonic deletions in established breast cancer susceptibility genes. We also detected suggestive associations with non-coding CNVs in known and novel loci with large effects sizes. Larger sample sizes will be required to reach robust levels of statistical significance.

Background: Epidemiological studies provide strong evidence for a role of endogenous sex hormones in the aetiology of breast cancer. The aim of this analysis was to identify genetic variants that are associated with urinary sex-hormone levels and breast cancer risk. Methods: We carried out a genome-wide association study of urinary oestrone-3-glucuronide and pregnanediol-3-glucuronide levels in 560 premenopausal women, with additional analysis of progesterone levels in 298 premenopausal women. To test for the association with breast cancer risk, we carried out follow-up genotyping in 90,916 cases and 89,893 controls from the Breast Cancer Association Consortium. All women were of European ancestry. Results: For pregnanediol-3-glucuronide, there were no genome-wide significant associations; for oestrone-3-glucuronide, we identified a single peak mapping to the CYP3A locus, annotated by rs45446698. The minor rs45446698-C allele was associated with lower oestrone-3-glucuronide (-49.2%, 95% CI -56.1% to -41.1%, P = 3.1 × 10¿18); in follow-up analyses, rs45446698-C was also associated with lower progesterone (-26.7%, 95% CI -39.4% to -11.6%, P = 0.001) and reduced risk of oestrogen and progesterone receptor-positive breast cancer (OR = 0.86, 95% CI 0.82¿0.91, P = 6.9 × 10¿8). Conclusions: The CYP3A7*1C allele is associated with reduced risk of hormone receptor-positive breast cancer possibly mediated via an effect on the metabolism of endogenous sex hormones in premenopausal women.

Background: We assessed whether blood cadmium levels were associated with incident lung cancer and could be used in the context of a screening program for early-stage lung cancer. Material and methods: We measured blood cadmium levels among 205 lung cancer patients and 205 matched controls. Cases and controls were matched for sex, age and smoking history (total pack-years, years since cessation for former smokers). Results: The odds ratio for those in the highest quartile of cadmium level (versus lowest) was four-fold (OR = 4.41, 95 % CI:2.01¿9.67, p < 0.01). The association was present in former smokers (OR = 16.8, 95 % CI:3.96-71.2, p < 0.01), but not in current smokers (OR = 1.23, 95 % CI: 0.34¿4.38) or in never smokers (OR not defined). Among former smokers, the association was present in both early- and late-stage lung cancer. Conclusion: Blood cadmium levels may be a marker to help with the early detection of lung cancer among former smokers.

Patients with metastatic castration-resistant prostate cancer (mCRPC) have an average survival of only 13 months. Identification of novel predictive and actionable biomarkers in the homologous recombination repair (HRR) pathway in up to a quarter of patients with mCRPC has led to the approval of targeted therapies like poly-ADP ribose polymerase inhibitors (PARPi), with the potential to improve survival outcomes. The approval of PARPi has led to guideline bodies such as the National Comprehensive Cancer Network (NCCN) to actively recommend germline and or somatic HRR gene panel testing to identify patients who will benefit from PARPi. However, there are several challenges as genetic testing is still at an early stage especially in low- and middle-income countries, with cost and availability being major impediments. In addition, there are issues such as choice of optimal tissue for genetic testing, archival, storage, retrieval of tissue blocks, interpretation and classification of variants in the HRR pathway, and the need for pretest and post-test genetic counseling. This review provides insights into the HRR gene mutations prevalent in mCRPC and the challenges for a more widespread gene testing to identify actionable germline pathogenic variants and somatic mutations in the HRR pathway, and proposes a clinical algorithm to enhance the efficiency of the gene testing process.

Background and objective: Chronic obstructive pulmonary disease (COPD) is the third leading cause of illness and death worldwide. Current treatments aim to control symptoms with none able to reverse disease or stop its progression. We explored the major molecular changes in COPD pathogenesis. Methods: We employed quantitative label-based proteomics to map the changes in the lung tissue proteome of cigarette smoke-induced experimental COPD that is induced over 8 weeks and progresses over 12 weeks. Results: Quantification of 7324 proteins enabled the tracking of changes to the proteome. Alterations in protein expression profiles occurred in the induction phase, with 18 and 16 protein changes at 4- and 6-week time points, compared to age-matched controls, respectively. Strikingly, 269 proteins had altered expression after 8 weeks when the hallmark pathological features of human COPD emerge, but this dropped to 27 changes at 12 weeks with disease progression. Differentially expressed proteins were validated using other mouse and human COPD bronchial biopsy samples. Major changes in RNA biosynthesis (heterogeneous nuclear ribonucleoproteins C1/C2 [HNRNPC] and RNA-binding protein Musashi homologue 2 [MSI2]) and modulators of inflammatory responses (S100A1) were notable. Mitochondrial dysfunction and changes in oxidative stress proteins also occurred. Conclusion: We provide a detailed proteomic profile, identifying proteins associated with the pathogenesis and disease progression of COPD establishing a platform to develop effective new treatment strategies.

We evaluated the joint associations between a new 313-variant PRS (PRS313) and questionnaire-based breast cancer risk factors for women of European ancestry, using 72 284 cases and 80 354 controls from the Breast Cancer Association Consortium. Interactions were evaluated using standard logistic regression and a newly developed case-only method for breast cancer risk overall and by estrogen receptor status. After accounting for multiple testing, we did not find evidence that per-standard deviation PRS313 odds ratio differed across strata defined by individual risk factors. Goodness-of-fit tests did not reject the assumption of a multiplicative model between PRS313 and each risk factor. Variation in projected absolute lifetime risk of breast cancer associated with classical risk factors was greater for women with higher genetic risk (PRS313 and family history) and, on average, 17.5% higher in the highest vs lowest deciles of genetic risk. These findings have implications for risk prevention for women at increased risk of breast cancer.

Background: Despite a modest association between tobacco smoking and breast cancer risk reported by recent epidemiological studies, it is still equivocal whether smoking is causally related to breast cancer risk. Methods: We applied Mendelian randomisation (MR) to evaluate a potential causal effect of cigarette smoking on breast cancer risk. Both individual-level data as well as summary statistics for 164 single-nucleotide polymorphisms (SNPs) reported in genome-wide association studies of lifetime smoking index (LSI) or cigarette per day (CPD) were used to obtain MR effect estimates. Data from 108,420 invasive breast cancer cases and 87,681 controls were used for the LSI analysis and for the CPD analysis conducted among ever-smokers from 26,147 cancer cases and 26,072 controls. Sensitivity analyses were conducted to address pleiotropy. Results: Genetically predicted LSI was associated with increased breast cancer risk (OR 1.18 per SD, 95% CI: 1.07¿1.30, P = 0.11 × 10¿2), but there was no evidence of association for genetically predicted CPD (OR 1.02, 95% CI: 0.78¿1.19, P = 0.85). The sensitivity analyses yielded similar results and showed no strong evidence of pleiotropic effect. Conclusion: Our MR study provides supportive evidence for a potential causal association with breast cancer risk for lifetime smoking exposure but not cigarettes per day among smokers.

Individuals with Lynch syndrome (LS), one of the most common inherited cancer syndromes, are at increased risk of developing malignancies, in particular colorectal cancer (CRC). Regular colonoscopy with polypectomy is recommended to reduce CRC risk in LS individuals. However, recent independent studies demonstrated that a substantial proportion of LS individuals develop CRC despite regular colonoscopy. The reasons for this surprising observation confirmed by large prospective studies are a matter of debate. In this review, we collect existing evidence from clinical, epidemiological and molecular studies and interpret them with regard to the origins and progression of LS-associated CRC. Alongside with hypotheses addressing colonoscopy quality and pace of progression from adenoma to cancer, we discuss the role of alternative precursors and immune system in LS-associated CRC. We also identify gaps in current knowledge and make suggestions for future studies aiming at improved CRC prevention for LS individuals.