microRNA and the Post-Transcriptional Response to Oxidative Stress during Neuronal Differentiation: Implications for Neurodevelopmental and Psychiatric Disorders.

Oxidative stress is one of the most important environmental exposures associated with psychiatric disorders, but the underlying molecular mechanisms remain to be elucidated. In a previous study, we observed a substantial alteration of the gene expression landscape in neuron-like cells that were differentiated from SH-SY5Y cells after or during exposure to oxidative stress, with a subset of dysregulated genes being enriched for neurodevelopmental processes. To further explore the regulatory mechanisms that might account for such profound perturbations, we have now applied small RNA-sequencing to investigate changes in the expression of miRNAs. These molecules are known to play crucial roles in brain development and response to stress through their capacity to suppress gene expression and influence complex biological networks. Through these analyses, we observed more than a hundred differentially expressed miRNAs, including 80 previously reported to be dysregulated in psychiatric disorders. The seven most influential miRNAs associated with pre-treatment exposure, including miR-138-5p, miR-96-5p, miR-34c-5p, miR-1287-5p, miR-497-5p, miR-195-5p, and miR-16-5p, supported by at least 10 negatively correlated mRNA connections, formed hubs in the interaction network with 134 genes enriched with neurobiological function, whereas in the co-treatment condition, miRNA-mRNA interaction pairs were enriched in cardiovascular and immunity-related disease ontologies. Interestingly, 12 differentially expressed miRNAs originated from the DLK1-DIO3 location, which encodes a schizophrenia-associated miRNA signature. Collectively, our findings suggest that early exposure to oxidative stress, before and during prenatal neuronal differentiation, might increase the risk of mental illnesses in adulthood by disturbing the expression of miRNAs that regulate neurodevelopmentally significant genes and networks.

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.

miRNA cargo in circulating vesicles from neurons is altered in individuals with schizophrenia and associated with severe disease.

While RNA expression appears to be altered in several brain disorders, the constraints of postmortem analysis make it impractical for well-powered population studies and biomarker development. Given that the unique molecular composition of neurons are reflected in their extracellular vesicles (EVs), we hypothesized that the fractionation of neuron derived EVs provides an opportunity to specifically profile their encapsulated contents noninvasively from blood. To investigate this hypothesis, we determined miRNA expression in microtubule associated protein 1B (MAP1B)-enriched serum EVs derived from neurons from a large cohort of individuals with schizophrenia and nonpsychiatric comparison participants. We observed dysregulation of miRNA in schizophrenia subjects, in particular those with treatment-resistance and severe cognitive deficits. These data support the hypothesis that schizophrenia is associated with alterations in posttranscriptional regulation of synaptic gene expression and provides an example of the potential utility of tissue-specific EV analysis in brain disorders.

MicroRNA binding site variation is enriched in psychiatric disorders.

Psychiatric disorders have a polygenic architecture, often associated with dozens or hundreds of independent genomic loci. Most associated loci impact noncoding regions of the genome, suggesting that the majority of disease heritability originates from the disruption of regulatory sequences. While most research has focused on variants that modify regulatory DNA elements, those affecting cis-acting RNA sequences, such as miRNA binding sites, are also likely to have a significant impact. We intersected genome-wide association study (GWAS) summary statistics with the dbMTS database of predictions for miRNA binding site variants (MBSVs). We compared the distributions of MBSV association statistics to non-MBSVs within brain-expressed 3'UTR regions. We aggregated GWAS p values at the gene, pathway, and miRNA family levels to investigate cellular functions and miRNA families strongly associated with each trait. We performed these analyses in several psychiatric disorders as well as nonpsychiatric traits for comparison. We observed significant enrichment of MBSVs in schizophrenia, depression, bipolar disorder, and anorexia nervosa, particularly in genes targeted by several miRNA families, including miR-335-5p, miR-21-5p/590-5p, miR-361-5p, and miR-557, and a nominally significant association between miR-323b-3p MBSVs and schizophrenia risk. We identified evidence for the association between MBSVs in synaptic gene sets in schizophrenia and bipolar disorder. We also observed a significant association of MBSVs in other complex traits including type 2 diabetes. These observations support the role of miRNA in the pathophysiology of psychiatric disorders and suggest that MBSVs are an important class of regulatory variants that have functional implications for many disorders, as well as other complex human traits.

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.

The genetic architecture of pneumonia susceptibility implicates mucin biology and a relationship with psychiatric illness.

Pneumonia remains one of the leading causes of death worldwide. In this study, we use genome-wide meta-analysis of lifetime pneumonia diagnosis (N = 391,044) to identify four association signals outside of the previously implicated major histocompatibility complex region. Integrative analyses and finemapping of these signals support clinically tractable targets, including the mucin MUC5AC and tumour necrosis factor receptor superfamily member TNFRSF1A. Moreover, we demonstrate widespread evidence of genetic overlap with pneumonia susceptibility across the human phenome, including particularly significant correlations with psychiatric phenotypes that remain significant after testing differing phenotype definitions for pneumonia or genetically conditioning on smoking behaviour. Finally, we show how polygenic risk could be utilised for precision treatment formulation or drug repurposing through pneumonia risk scores constructed using variants mapped to pathways with known drug targets. In summary, we provide insights into the genetic architecture of pneumonia susceptibility and genetics informed targets for drug development or repositioning.

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.

Genetic association and causal inference converge on hyperglycaemia as a modifiable factor to improve lung function.

Measures of lung function are heritable, and thus, we sought to utilise genetics to propose drug-repurposing candidates that could improve respiratory outcomes. Lung function measures were found to be genetically correlated with seven druggable biochemical traits, with further evidence of a causal relationship between increased fasting glucose and diminished lung function. Moreover, we developed polygenic scores for lung function specifically within pathways with known drug targets and investigated their relationship with pulmonary phenotypes and gene expression in independent cohorts to prioritise individuals who may benefit from particular drug-repurposing opportunities. A transcriptome-wide association study (TWAS) of lung function was then performed which identified several drug-gene interactions with predicted lung function increasing modes of action. Drugs that regulate blood glucose were uncovered through both polygenic scoring and TWAS methodologies. In summary, we provided genetic justification for a number of novel drug-repurposing opportunities that could improve lung function.

Chronic respiratory disorders like asthma affect around 600 million people worldwide. Although these illnesses are widespread, they can have several different underlying causes, making them difficult to treat. Drugs that work well on one type of respiratory disorder may be completely ineffective on another. Understanding the biological and environmental factors that cause these illnesses will allow them to be treated more effectively by tailoring therapies to each patient. Reduced lung function is a factor in respiratory disorders and it can have many genetic causes. Studying the genes of patients with reduced lung function can reveal the genes involved, some of which may already be targets of existing drugs for other illnesses. So, could a patient's genetics be used to repurpose existing drugs to treat their respiratory disorders? Reay et al. combined three methods to link genetics and biological processes to the causes of reduced lung function. The results reveal several factors that could lead to new treatments. In one example, reduced lung function showed a link to genes associated with high blood sugar. As such, treatments used in diabetes might help improve lung function in some patients. Reay et al. also developed a scoring system that could predict the efficacy of a treatment based on a patient's genetics. The study suggests that COVID-19 infection could be affected by blood sugar levels too. Chronic respiratory disorders are a critical issue worldwide and have proven difficult to treat, but these results suggest a way to identify new therapies and target them to the right patients. The findings also support a connection between lung function and blood sugar levels. This implies that perhaps existing diabetes treatments ­ including diet and lifestyle changes aimed at reducing or limiting blood sugar ­ could be repurposed to treat respiratory disorders in some patients. The next step will be to perform clinical trials to test whether these therapies are in fact effective.

Investigation of glycaemic traits in psychiatric disorders using Mendelian randomisation revealed a causal relationship with anorexia nervosa.

Data from observational studies have suggested an involvement of abnormal glycaemic regulation in the pathophysiology of psychiatric illness. This may be an attractive target for clinical intervention as glycaemia can be modulated by both lifestyle factors and pharmacological agents. However, observational studies are inherently confounded, and therefore, causal relationships cannot be reliably established. We employed genetic variants rigorously associated with three glycaemic traits (fasting glucose, fasting insulin, and glycated haemoglobin) as instrumental variables in a two-sample Mendelian randomisation analysis to investigate the causal effect of these measures on the risk for eight psychiatric disorders. A significant protective effect of a natural log transformed pmol/L increase in fasting insulin levels was observed for anorexia nervosa after the application of multiple testing correction (OR = 0.48 [95% CI: 0.33-0.71]-inverse-variance weighted estimate). There was no consistently strong evidence for a causal effect of glycaemic factors on the other seven psychiatric disorders considered. The relationship between fasting insulin and anorexia nervosa was supported by a suite of sensitivity analyses, with no statistical evidence of instrument heterogeneity or horizontal pleiotropy. Further investigation is required to explore the relationship between insulin levels and anorexia.

Oxidative Stress Impact on the Transcriptome of Differentiating Neuroblastoma Cells: Implication for Psychiatric Disorders.

Prenatal environmental exposures that have been shown to induce oxidative stress (OS) during pregnancy, such as smoking and alcohol consumption, are risk factors for the onset of schizophrenia and other neurodevelopmental disorders (NDDs). While the OS role in the etiology of neurodegenerative diseases is well known, its contribution to the genomic dysregulation associated with psychiatric disorders is less well defined. In this study we used the SH-SY5Y cell line and applied RNA-sequencing to explore transcriptomic changes in response to OS before or during neural differentiation. We observed differential expression of many genes, most of which localised to the synapse and were involved in neuronal differentiation. These genes were enriched in schizophrenia-associated signalling pathways, including PI3K/Akt, axon guidance, and signalling by retinoic acid. Interestingly, circulatory system development was affected by both treatments, which is concordant with observations of increased prevalence of cardiovascular disease in patients with NDDs. We also observed a very significant increase in the expression of immunity-related genes, supporting current hypotheses of immune system involvement in psychiatric disorders. While further investigation of this influence in other cell and animal models is warranted, our data suggest that early life exposure to OS has a disruptive influence on neuronal gene expression that may perturb normal differentiation and neurodevelopment, thereby contributing towards overall risk for developing psychiatric diseases.

Characterising the Transcriptional and Translational Impact of the Schizophrenia-Associated miR-1271-5p in Neuronal Cells.

MicroRNA (miRNA) coordinate complex gene expression networks in cells that are vital to support highly specialised morphology and cytoarchitecture. Neurons express a rich array of miRNA, including many that are specific or enriched, which have important functions in this context and implications for neurological conditions. While the neurological function of a number of brain-derived miRNAs have been examined thoroughly, the mechanistic basis of many remain obscure. In this case, we investigated the transcriptome-wide impact of schizophrenia-associated miR-1271-5p in response to bidirectional modulation. Alteration of miR-1271-5p induced considerable changes to mRNA abundance and translation, which spanned a diverse range of cellular functions, including directly targeted genes strongly associated with cytoskeletal dynamics and cellular junctions. Mechanistic analyses additionally revealed that upregulation of miR-1271-5p predominantly repressed mRNAs through destabilisation, wherein 3'UTR and coding sequence binding sites exhibited similar efficacy. Knockdown, however, produced no discernible trend in target gene expression and strikingly resulted in increased expression of the highly conserved miR-96-5p, which shares an identical seed region with miR-1271-5p, suggesting the presence of feedback mechanisms that sense disruptions to miRNA levels. These findings indicate that, while bidirectional regulation of miR-1271-5p results in substantial remodeling of the neuronal transcriptome, these effects are not inverse in nature. In addition, we provide further support for the idea that destabilisation of mRNA is the predominant mechanism by which miRNAs regulate complementary mRNAs.

Alteration of miRNA-mRNA interactions in lymphocytes of individuals with schizophrenia.

The aetiology of schizophrenia is complex, heterogeneous, and involves interplay of many genetic and environmental influences. While significant progress has been made in the understanding the common heritable component, we are still grappling with the genomic encoding of environmental risk. One class of molecule that has tremendous potential is miRNA. These molecules are regulated by genetic and environmental factors associated with schizophrenia and have a very significant impact on temporospatial patterns of gene expression. To better understand the relationship between miRNA and gene expression in the disorder we analysed these molecules in RNA isolated from peripheral blood mononuclear cells (PBMCs) obtained from an Australian cohort of 36 individuals with schizophrenia and 15 healthy controls using next-generation RNA sequencing. Significant changes in both mRNA and miRNA expression profiles were observed implicating important interaction networks involved in immune activity and development. We also observed sexual dimorphism, particularly in relation to variation in mRNA, with males showing significantly more differentially expressed genes. Interestingly, while we explored expression in lymphocytes, the systems biology of miRNA-mRNA interactions was suggestive of significant pleiotropy with enrichment of networks related to neuronal activity.

Circular RNA biogenesis is decreased in postmortem cortical gray matter in schizophrenia and may alter the bioavailability of associated miRNA.

Circular RNAs (circRNAs) are a covalently closed subclass of non-coding RNA molecules formed by back splicing of linear precursor RNA. These molecules are relatively stable and particularly abundant in the mammalian brain and therefore may participate in neural development and function. With the emergence of circRNAs activity in gene regulation, these molecules have been implicated in several biological processes, including synaptic plasticity, and we therefore suspect they may have a role in neurobehavioral disorders. Here, we profile cortical circRNAs expression in 35 postmortem cortical gray matter (BA46) schizophrenia and a non-psychiatric comparison group, using circRNA enrichment sequencing. While more than 90,000 circRNAs species were identified in the dorsolateral prefrontal cortex (DLPFC), we observed lower complexity and substantial depletion in subjects with the disorder. Although circRNAs expression was independent of their host gene transcription, alternative splicing rates were lower in samples from cases compared to controls. Gene set analysis of differentially expressed circRNAs host genes revealed significant enrichment of neural functions and neurological disorders. Many of these depleted circRNAs are also predicted to sequester miRNAs that were shown previously to be increased in the disorder, potentially exacerbating the functional impact of their dysregulation through posttranscriptional gene silencing. While this is the first reported exploration of circRNAs in schizophrenia, there is significant potential for dysregulation more broadly in other major mental illnesses and behavioral disorders. Given their capacity for modulating miRNA function, circRNA may play a significant role in the pathophysiology of disease and even be targeted for therapeutic manipulation.