The Complexity of Tobacco Smoke-Induced Mutagenesis in Head and Neck Cancer.

Tobacco smoke, alone or combined with alcohol, is the predominant cause of head and neck cancer (HNC). Here, we further explore how tobacco exposure contributes to cancer development by mutational signature analysis of 265 whole-genome sequenced HNC from eight countries. Six tobacco-associated mutational signatures were detected, including some not previously reported. Differences in HNC incidence between countries corresponded with differences in mutation burdens of tobacco-associated signatures, consistent with the dominant role of tobacco in HNC causation. Differences were found in the burden of tobacco-associated signatures between anatomical subsites, suggesting that tissue-specific factors modulate mutagenesis. We identified an association between tobacco smoking and three additional alcohol-related signatures indicating synergism between the two exposures. Tobacco smoking was associated with differences in the mutational spectra and repertoire of driver mutations in cancer genes, and in patterns of copy number change. Together, the results demonstrate the multiple pathways by which tobacco smoke can influence the evolution of cancer cell clones.

Identifying proteomic risk factors for cancer using prospective and exome analyses of 1463 circulating proteins and risk of 19 cancers in the UK Biobank.

The availability of protein measurements and whole exome sequence data in the UK Biobank enables investigation of potential observational and genetic protein-cancer risk associations. We investigated associations of 1463 plasma proteins with incidence of 19 cancers and 9 cancer subsites in UK Biobank participants (average 12 years follow-up). Emerging protein-cancer associations were further explored using two genetic approaches, cis-pQTL and exome-wide protein genetic scores (exGS). We identify 618 protein-cancer associations, of which 107 persist for cases diagnosed more than seven years after blood draw, 29 of 618 were associated in genetic analyses, and four had support from long time-to-diagnosis ( > 7 years) and both cis-pQTL and exGS analyses: CD74 and TNFRSF1B with NHL, ADAM8 with leukemia, and SFTPA2 with lung cancer. We present multiple blood protein-cancer risk associations, including many detectable more than seven years before cancer diagnosis and that had concordant evidence from genetic analyses, suggesting a possible role in cancer development.

Genetic variants for head size share genes and pathways with cancer.

The size of the human head is highly heritable, but genetic drivers of its variation within the general population remain unmapped. We perform a genome-wide association study on head size (N = 80,890) and identify 67 genetic loci, of which 50 are novel. Neuroimaging studies show that 17 variants affect specific brain areas, but most have widespread effects. Gene set enrichment is observed for various cancers and the p53, Wnt, and ErbB signaling pathways. Genes harboring lead variants are enriched for macrocephaly syndrome genes (37-fold) and high-fidelity cancer genes (9-fold), which is not seen for human height variants. Head size variants are also near genes preferentially expressed in intermediate progenitor cells, neural cells linked to evolutionary brain expansion. Our results indicate that genes regulating early brain and cranial growth incline to neoplasia later in life, irrespective of height. This warrants investigation of clinical implications of the link between head size and cancer.

Alteration of DNA Methylation and Epigenetic Scores Associated With Features of Schizophrenia and Common Variant Genetic Risk.

BACKGROUND: Unpacking molecular perturbations associated with features of schizophrenia is a critical step toward understanding phenotypic heterogeneity in this disorder. Recent epigenome-wide association studies have uncovered pervasive dysregulation of DNA methylation in schizophrenia; however, clinical features of the disorder that account for a large proportion of phenotypic variability are relatively underexplored. METHODS: We comprehensively analyzed patterns of DNA methylation in a cohort of 381 individuals with schizophrenia from the deeply phenotyped Australian Schizophrenia Research Bank. Epigenetic changes were investigated in association with cognitive status, age of onset, treatment resistance, Global Assessment of Functioning scores, and common variant polygenic risk scores for schizophrenia. We subsequently explored alterations within genes previously associated with psychiatric illness, phenome-wide epigenetic covariance, and epigenetic scores. RESULTS: Epigenome-wide association studies of the 5 primary traits identified 662 suggestively significant (p < 6.72 × 10-5) differentially methylated probes, with a further 432 revealed after controlling for schizophrenia polygenic risk on the remaining 4 traits. Interestingly, we uncovered many probes within genes associated with a variety of psychiatric conditions as well as significant epigenetic covariance with phenotypes and exposures including acute myocardial infarction, C-reactive protein, and lung cancer. Epigenetic scores for treatment-resistant schizophrenia strikingly exhibited association with clozapine administration, while epigenetic proxies of plasma protein expression, such as CCL17, MMP10, and PRG2, were associated with several features of schizophrenia. CONCLUSIONS: Our findings collectively provide novel evidence suggesting that several features of schizophrenia are associated with alteration of DNA methylation, which may contribute to interindividual phenotypic variation in affected individuals.

Beyond the Global Brain Differences: Intraindividual Variability Differences in 1q21.1 Distal and 15q11.2 BP1-BP2 Deletion Carriers.

BACKGROUND: Carriers of the 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants exhibit regional and global brain differences compared with noncarriers. However, interpreting regional differences is challenging if a global difference drives the regional brain differences. Intraindividual variability measures can be used to test for regional differences beyond global differences in brain structure. METHODS: Magnetic resonance imaging data were used to obtain regional brain values for 1q21.1 distal deletion (n = 30) and duplication (n = 27) and 15q11.2 BP1-BP2 deletion (n = 170) and duplication (n = 243) carriers and matched noncarriers (n = 2350). Regional intra-deviation scores, i.e., the standardized difference between an individual's regional difference and global difference, were used to test for regional differences that diverge from the global difference. RESULTS: For the 1q21.1 distal deletion carriers, cortical surface area for regions in the medial visual cortex, posterior cingulate, and temporal pole differed less and regions in the prefrontal and superior temporal cortex differed more than the global difference in cortical surface area. For the 15q11.2 BP1-BP2 deletion carriers, cortical thickness in regions in the medial visual cortex, auditory cortex, and temporal pole differed less and the prefrontal and somatosensory cortex differed more than the global difference in cortical thickness. CONCLUSIONS: We find evidence for regional effects beyond differences in global brain measures in 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants. The results provide new insight into brain profiling of the 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants, with the potential to increase understanding of the mechanisms involved in altered neurodevelopment.

Identifying proteomic risk factors for overall, aggressive and early onset prostate cancer using Mendelian randomization and tumor spatial transcriptomics.

Background: Understanding the role of circulating proteins in prostate cancer risk can reveal key biological pathways and identify novel targets for cancer prevention. Methods: We investigated the association of 2,002 genetically predicted circulating protein levels with risk of prostate cancer overall, and of aggressive and early onset disease, using cis-pQTL Mendelian randomization (MR) and colocalization. Findings for proteins with support from both MR, after correction for multiple-testing, and colocalization were replicated using two independent cancer GWAS, one of European and one of African ancestry. Proteins with evidence of prostate-specific tissue expression were additionally investigated using spatial transcriptomic data in prostate tumor tissue to assess their role in tumor aggressiveness. Finally, we mapped risk proteins to drug and ongoing clinical trials targets. Results: We identified 20 proteins genetically linked to prostate cancer risk (14 for overall [8 specific], 7 for aggressive [3 specific], and 8 for early onset disease [2 specific]), of which a majority were novel and replicated. Among these were proteins associated with aggressive disease, such as PPA2 [Odds Ratio (OR) per 1 SD increment = 2.13, 95% CI: 1.54-2.93], PYY [OR = 1.87, 95% CI: 1.43-2.44] and PRSS3 [OR = 0.80, 95% CI: 0.73-0.89], and those associated with early onset disease, including EHPB1 [OR = 2.89, 95% CI: 1.99-4.21], POGLUT3 [OR = 0.76, 95% CI: 0.67-0.86] and TPM3 [OR = 0.47, 95% CI: 0.34-0.64]. We confirm an inverse association of MSMB with prostate cancer overall [OR = 0.81, 95% CI: 0.80-0.82], and also find an inverse association with both aggressive [OR = 0.84, 95% CI: 0.82-0.86] and early onset disease [OR = 0.71, 95% CI: 0.68-0.74]. Using spatial transcriptomics data, we identified MSMB as the genome-wide top-most predictive gene to distinguish benign regions from high grade cancer regions that had five-fold lower MSMB expression. Additionally, ten proteins that were associated with prostate cancer risk mapped to existing therapeutic interventions. Conclusion: Our findings emphasize the importance of proteomics for improving our understanding of prostate cancer etiology and of opportunities for novel therapeutic interventions. Additionally, we demonstrate the added benefit of in-depth functional analyses to triangulate the role of risk proteins in the clinical aggressiveness of prostate tumors. Using these integrated methods, we identify a subset of risk proteins associated with aggressive and early onset disease as priorities for investigation for the future prevention and treatment of prostate cancer.

Common genetic variations in telomere length genes and lung cancer: a Mendelian randomisation study and its novel application in lung tumour transcriptome.

Background: Genome-wide association studies (GWASs) have identified genetic susceptibility variants for both leukocyte telomere length (LTL) and lung cancer susceptibility. Our study aims to explore the shared genetic basis between these traits and investigate their impact on somatic environment of lung tumours. Methods: We performed genetic correlation, Mendelian randomisation (MR), and colocalisation analyses using the largest available GWASs summary statistics of LTL (N=464,716) and lung cancer (N=29,239 cases and 56,450 controls). Principal components analysis based on RNA-sequencing data was used to summarise gene expression profile in lung adenocarcinoma cases from TCGA (N=343). Results: Although there was no genome-wide genetic correlation between LTL and lung cancer risk, longer LTL conferred an increased risk of lung cancer regardless of smoking status in the MR analyses, particularly for lung adenocarcinoma. Of the 144 LTL genetic instruments, 12 colocalised with lung adenocarcinoma risk and revealed novel susceptibility loci, including MPHOSPH6, PRPF6, and POLI. The polygenic risk score for LTL was associated with a specific gene expression profile (PC2) in lung adenocarcinoma tumours. The aspect of PC2 associated with longer LTL was also associated with being female, never smokers, and earlier tumour stages. PC2 was strongly associated with cell proliferation score and genomic features related to genome stability, including copy number changes and telomerase activity. Conclusions: This study identified an association between longer genetically predicted LTL and lung cancer and sheds light on the potential molecular mechanisms related to LTL in lung adenocarcinomas. Funding: Institut National du Cancer (GeniLuc2017-1-TABAC-03-CIRC-1-TABAC17-022), INTEGRAL/NIH (5U19CA203654-03), CRUK (C18281/A29019), and Agence Nationale pour la Recherche (ANR-10-INBS-09).

Genetics-informed precision treatment formulation in schizophrenia and bipolar disorder.

Genetically informed drug development and repurposing is an attractive prospect for improving patient outcomes in psychiatry; however, the effectiveness of these endeavors is confounded by heterogeneity. We propose an approach that links interventions implicated by disorder-associated genetic risk, at the population level, to a framework that can target these compounds to individuals. Specifically, results from genome-wide association studies are integrated with expression data to prioritize individual "directional anchor" genes for which the predicted risk-increasing direction of expression could be counteracted by an existing drug. While these compounds represent plausible therapeutic candidates, they are not likely to be equally efficacious for all individuals. To account for this heterogeneity, we constructed polygenic scores restricted to variants annotated to the network of genes that interact with each directional anchor gene. These metrics, which we call a pharmagenic enrichment score (PES), identify individuals with a higher burden of genetic risk, localized in biological processes related to the candidate drug target, to inform precision drug repurposing. We used this approach to investigate schizophrenia and bipolar disorder and reveal several compounds targeting specific directional anchor genes that could be plausibly repurposed. These genetic risk scores, mapped to the networks associated with target genes, revealed biological insights that cannot be observed in undifferentiated genome-wide polygenic risk score (PRS). For example, an enrichment of these partitioned scores in schizophrenia cases with otherwise low PRS. In summary, genetic risk could be used more specifically to direct drug repurposing candidates that target particular genes implicated in psychiatric and other complex disorders.

An international report on bacterial communities in esophageal squamous cell carcinoma.

The incidence of esophageal squamous cell carcinoma (ESCC) is disproportionately high in the eastern corridor of Africa and parts of Asia. Emerging research has identified a potential association between poor oral health and ESCC. One possible link between poor oral health and ESCC involves the alteration of the microbiome. We performed an integrated analysis of four independent sequencing efforts of ESCC tumors from patients from high- and low-incidence regions of the world. Using whole genome sequencing (WGS) and RNA sequencing (RNAseq) of ESCC tumors from 61 patients in Tanzania, we identified a community of bacteria, including members of the genera Fusobacterium, Selenomonas, Prevotella, Streptococcus, Porphyromonas, Veillonella and Campylobacter, present at high abundance in ESCC tumors. We then characterized the microbiome of 238 ESCC tumor specimens collected in two additional independent sequencing efforts consisting of patients from other high-ESCC incidence regions (Tanzania, Malawi, Kenya, Iran, China). This analysis revealed similar ESCC-associated bacterial communities in these cancers. Because these genera are traditionally considered members of the oral microbiota, we next explored whether there was a relationship between the synchronous saliva and tumor microbiomes of ESCC patients in Tanzania. Comparative analyses revealed that paired saliva and tumor microbiomes were significantly similar with a specific enrichment of Fusobacterium and Prevotella in the tumor microbiome. Together, these data indicate that cancer-associated oral bacteria are associated with ESCC tumors at the time of diagnosis and support a model in which oral bacteria are present in high abundance in both saliva and tumors of some ESCC patients.

Association of germline TYK2 variation with lung cancer and non-Hodgkin lymphoma risk.

Deucravacitinib, a novel, selective inhibitor of TYK2 is currently under review at the FDA and EMA for treatment of moderate-to-severe plaque psoriasis. It is unclear whether recent safety concerns (ie, elevated rates of lung cancer and lymphoma) related to similar medications (ie, other JAK inhibitors) are shared with this novel TYK2 inhibitor. We used a partial loss-of-function variant in TYK2 (rs34536443), previously shown to protect against psoriasis and other autoimmune diseases, to evaluate the potential effect of therapeutic TYK2 inhibition on risk of lung cancer and non-Hodgkin lymphoma. Summary genetic association data on lung cancer risk were obtained from a GWAS meta-analysis of 29 266 cases and 56 450 controls in the Integrative Analysis of Lung Cancer Risk and Aetiology (INTEGRAL) consortium. Summary genetic association data on non-Hodgkin lymphoma risk were obtained from a GWAS meta-analysis of 8489 cases and 374 506 controls in the UK Biobank and InterLymph consortium. In the primary analysis, each copy of the minor allele of rs34536443, representing partial TYK2 inhibition, was associated with an increased risk of lung cancer (OR 1.15, 95% CI 1.09-1.23, P = 2.29 × 10-6 ) and non-Hodgkin lymphoma (OR 1.18, 95% CI 1.05-1.33, P = 5.25 × 10-3 ). Our analyses using an established partial loss-of-function mutation to mimic TYK2 inhibition provide genetic evidence that therapeutic TYK2 inhibition may increase risk of lung cancer and non-Hodgkin lymphoma. These findings, consistent with recent reports from postmarketing trials of similar JAK inhibitors, could have important implications for future safety assessment of deucravacitinib and other TYK2 inhibitors in development.

Genetic Analysis of Lung Cancer and the Germline Impact on Somatic Mutation Burden.

BACKGROUND: Germline genetic variation contributes to lung cancer (LC) susceptibility. Previous genome-wide association studies (GWAS) have implicated susceptibility loci involved in smoking behaviors and DNA repair genes, but further work is required to identify susceptibility variants. METHODS: To identify LC susceptibility loci, a family history-based genome-wide association by proxy (GWAx) of LC (48 843 European proxy LC patients, 195 387 controls) was combined with a previous LC GWAS (29 266 patients, 56 450 controls) by meta-analysis. Colocalization was used to explore candidate genes and overlap with existing traits at discovered susceptibility loci. Polygenic risk scores (PRS) were tested within an independent validation cohort (1 666 LC patients vs 6 664 controls) using variants selected from the LC susceptibility loci and a novel selection approach using published GWAS summary statistics. Finally, the effects of the LC PRS on somatic mutational burden were explored in patients whose tumor resections have been profiled by exome (n = 685) and genome sequencing (n = 61). Statistical tests were 2-sided. RESULTS: The GWAx-GWAS meta-analysis identified 8 novel LC loci. Colocalization implicated DNA repair genes (CHEK1), metabolic genes (CYP1A1), and smoking propensity genes (CHRNA4 and CHRNB2). PRS analysis demonstrated that these variants, as well as subgenome-wide significant variants related to expression quantitative trait loci and/or smoking propensity, assisted in LC genetic risk prediction (odds ratio = 1.37, 95% confidence interval = 1.29 to 1.45; P < .001). Patients with higher genetic PRS loads of smoking-related variants tended to have higher mutation burdens in their lung tumors. CONCLUSIONS: This study has expanded the number of LC susceptibility loci and provided insights into the molecular mechanisms by which these susceptibility variants contribute to LC development.

Genetic estimates of correlation and causality between blood-based biomarkers and psychiatric disorders.

There is a long-standing interest in exploring the relationship between blood-based biomarkers and psychiatric disorders, despite their causal role being difficult to resolve in observational studies. In this study, we leverage genome-wide association study data for a large panel of heritable serum biochemical traits to refine our understanding of causal effect in biochemical-psychiatric trait pairings. We observed widespread positive and negative genetic correlation between psychiatric disorders and biochemical traits. Causal inference was then implemented to distinguish causation from correlation, with strong evidence that C-reactive protein (CRP) exerts a causal effect on psychiatric disorders. Notably, CRP demonstrated both protective and risk-increasing effects on different disorders. Multivariable models that conditioned CRP effects on interleukin-6 signaling and body mass index supported that the CRP-schizophrenia relationship was not driven by these factors. Collectively, these data suggest that there are shared pathways that influence both biochemical traits and psychiatric illness.

Cell type-specific manifestations of cortical thickness heterogeneity in schizophrenia.

Brain morphology differs markedly between individuals with schizophrenia, but the cellular and genetic basis of this heterogeneity is poorly understood. Here, we sought to determine whether cortical thickness (CTh) heterogeneity in schizophrenia relates to interregional variation in distinct neural cell types, as inferred from established gene expression data and person-specific genomic variation. This study comprised 1849 participants in total, including a discovery (140 cases and 1267 controls) and a validation cohort (335 cases and 185 controls). To characterize CTh heterogeneity, normative ranges were established for 34 cortical regions and the extent of deviation from these ranges was measured for each individual with schizophrenia. CTh deviations were explained by interregional gene expression levels of five out of seven neural cell types examined: (1) astrocytes; (2) endothelial cells; (3) oligodendrocyte progenitor cells (OPCs); (4) excitatory neurons; and (5) inhibitory neurons. Regional alignment between CTh alterations with cell type transcriptional maps distinguished broad patient subtypes, which were validated against genomic data drawn from the same individuals. In a predominantly neuronal/endothelial subtype (22% of patients), CTh deviations covaried with polygenic risk for schizophrenia (sczPRS) calculated specifically from genes marking neuronal and endothelial cells (r = -0.40, p = 0.010). Whereas, in a predominantly glia/OPC subtype (43% of patients), CTh deviations covaried with sczPRS calculated from glia and OPC-linked genes (r = -0.30, p = 0.028). This multi-scale analysis of genomic, transcriptomic, and brain phenotypic data may indicate that CTh heterogeneity in schizophrenia relates to inter-individual variation in cell-type specific functions. Decomposing heterogeneity in relation to cortical cell types enables prioritization of schizophrenia subsets for future disease modeling efforts.

The MIR137 VNTR rs58335419 Is Associated With Cognitive Impairment in Schizophrenia and Altered Cortical Morphology.

Genome-wide association studies (GWAS) of schizophrenia have strongly implicated a risk locus in close proximity to the gene for miR-137. While there are candidate single-nucleotide polymorphisms (SNPs) with functional implications for the microRNA's expression encompassed by the common haplotype tagged by rs1625579, there are likely to be others, such as the variable number tandem repeat (VNTR) variant rs58335419, that have no proxy on the SNP genotyping platforms used in GWAS to date. Using whole-genome sequencing data from schizophrenia patients (n = 299) and healthy controls (n = 131), we observed that the MIR137 4-repeats VNTR (VNTR4) variant was enriched in a cognitive deficit subtype of schizophrenia and associated with altered brain morphology, including thicker left inferior temporal gyrus and deeper right postcentral sulcus. These findings suggest that the MIR137 VNTR4 may impact neuroanatomical development that may, in turn, influence the expression of more severe cognitive symptoms in patients with schizophrenia.

The genetic architecture of the human cerebral cortex.

The cerebral cortex underlies our complex cognitive capabilities, yet little is known about the specific genetic loci that influence human cortical structure. To identify genetic variants that affect cortical structure, we conducted a genome-wide association meta-analysis of brain magnetic resonance imaging data from 51,665 individuals. We analyzed the surface area and average thickness of the whole cortex and 34 regions with known functional specializations. We identified 199 significant loci and found significant enrichment for loci influencing total surface area within regulatory elements that are active during prenatal cortical development, supporting the radial unit hypothesis. Loci that affect regional surface area cluster near genes in Wnt signaling pathways, which influence progenitor expansion and areal identity. Variation in cortical structure is genetically correlated with cognitive function, Parkinson's disease, insomnia, depression, neuroticism, and attention deficit hyperactivity disorder.

Pharmacological enrichment of polygenic risk for precision medicine in complex disorders.

Individuals with complex disorders typically have a heritable burden of common variation that can be expressed as a polygenic risk score (PRS). While PRS has some predictive utility, it lacks the molecular specificity to be directly informative for clinical interventions. We therefore sought to develop a framework to quantify an individual's common variant enrichment in clinically actionable systems responsive to existing drugs. This was achieved with a metric designated the pharmagenic enrichment score (PES), which we demonstrate for individual SNP profiles in a cohort of cases with schizophrenia. A large proportion of these had elevated PES in one or more of eight clinically actionable gene-sets enriched with schizophrenia associated common variation. Notable candidates targeting these pathways included vitamins, antioxidants, insulin modulating agents, and cholinergic drugs. Interestingly, elevated PES was also observed in individuals with otherwise low common variant burden. The biological saliency of PES profiles were observed directly through their impact on gene expression in a subset of the cohort with matched transcriptomic data, supporting our assertion that this gene-set orientated approach could integrate an individual's common variant risk to inform personalised interventions, including drug repositioning, for complex disorders such as schizophrenia.

The Medical Genome Reference Bank contains whole genome and phenotype data of 2570 healthy elderly.

Population health research is increasingly focused on the genetic determinants of healthy ageing, but there is no public resource of whole genome sequences and phenotype data from healthy elderly individuals. Here we describe the first release of the Medical Genome Reference Bank (MGRB), comprising whole genome sequence and phenotype of 2570 elderly Australians depleted for cancer, cardiovascular disease, and dementia. We analyse the MGRB for single-nucleotide, indel and structural variation in the nuclear and mitochondrial genomes. MGRB individuals have fewer disease-associated common and rare germline variants, relative to both cancer cases and the gnomAD and UK Biobank cohorts, consistent with risk depletion. Age-related somatic changes are correlated with grip strength in men, suggesting blood-derived whole genomes may also provide a biologic measure of age-related functional deterioration. The MGRB provides a broadly applicable reference cohort for clinical genetics and genomic association studies, and for understanding the genetics of healthy ageing.

Polygenic disruption of retinoid signalling in schizophrenia and a severe cognitive deficit subtype.

Retinoid metabolites of vitamin A are intrinsically linked to neural development, connectivity and plasticity, and have been implicated in the pathophysiology of schizophrenia. We hypothesised that a greater burden of common and rare genomic variation in genes involved with retinoid biogenesis and signalling could be associated with schizophrenia and its cognitive symptoms. Common variants associated with schizophrenia in the largest genome-wide association study were aggregated in retinoid genes and used to formulate a polygenic risk score (PRSRet) for each participant in the Australian Schizophrenia Research Bank. In support of our hypothesis, we found PRSRet to be significantly associated with the disorder. Cases with severe cognitive deficits, while not further differentiated by PRSRet, were enriched with rare variation in the retinoic acid receptor beta gene RARB, detected through whole-genome sequencing. RARB rare variant burden was also associated with reduced cerebellar volume in the cases with marked cognitive deficit, and with covariation in grey matter throughout the brain. An excess of rare variation was further observed in schizophrenia in retinoic acid response elements proximal to target genes, which we show are differentially expressed in the disorder in two RNA sequencing datasets. Our results suggest that genomic variation may disrupt retinoid signalling in schizophrenia, with particular significance for cases with severe cognitive impairment.

Wnt receptor gene FZD1 was associated with schizophrenia in genome-wide SNP analysis of the Australian Schizophrenia Research Bank cohort.

OBJECTIVES: Large-scale genetic analysis of common variation in schizophrenia has been a powerful approach to understanding this complex but highly heritable psychotic disorder. To further investigate loci, genes and pathways associated more specifically in the well-characterized Australian Schizophrenia Research Bank cohort, we applied genome-wide single-nucleotide polymorphism analysis in these three annotation categories. METHODS: We performed a case-control genome-wide association study in 429 schizophrenia samples and 255 controls. Post-genome-wide association study analyses were then integrated with genomic annotations to explore the enrichment of variation at the gene and pathway level. We also examine candidate single-nucleotide polymorphisms with potential function within expression quantitative trait loci and investigate overall enrichment of variation within tissue-specific functional regulatory domains of the genome. RESULTS: The strongest finding (p = 2.01 × 10-6, odds ratio = 1.82, 95% confidence interval = [1.42, 2.33]) in genome-wide association study was with rs10252923 at 7q21.13, downstream of FZD1 (frizzled class receptor 1). While this did not stand alone after correction, the involvement of FZD1 was supported by gene-based analysis, which exceeded the threshold for genome-wide significance (p = 2.78 × 10-6). CONCLUSION: The identification of FZD1, as an independent association signal at the gene level, supports the hypothesis that the Wnt signalling pathway is altered in the pathogenesis of schizophrenia and may be an important target for therapeutic development.

The maternal immune activation model uncovers a role for the Arx gene in GABAergic dysfunction in schizophrenia.

A hallmark feature of schizophrenia is altered high frequency neural oscillations, including reduced auditory-evoked gamma oscillatory power, which is underpinned by parvalbumin (PV) interneuron dysfunction. Maternal immune activation (MIA) in rodents models an environmental risk factor for schizophrenia and recapitulates these PV interneuron changes. This study sought to link reduced PV expression in the MIA model with alterations to auditory-evoked gamma oscillations and transcript expression. We further aligned transcriptional findings from the animal model with human genome sequencing data. We show that MIA, induced by the viral mimetic, poly-I:C in C57Bl/6 mice, caused in adult offspring reduced auditory-evoked gamma and theta oscillatory power paralleled by reduced PV protein levels. We then showed the Arx gene, critical to healthy neurodevelopment of PV interneurons, is reduced in the forebrain of MIA exposed mice. Finally, in a whole-genome sequenced patient cohort, we identified a novel missense mutation of ARX in a patient with schizophrenia and in the Psychiatric Genomics Consortium 2 cohort, a nominal association of proximal ARX SNPs with the disorder. This suggests MIA, as a risk factor for schizophrenia, may be influencing Arx expression to induce the GABAergic dysfunction seen in schizophrenia and that the ARX gene may play a role in the prenatal origins of schizophrenia pathophysiology.