Skip to main content

ORIGINAL RESEARCH article

Front. Genet., 02 October 2020
Sec. Behavioral and Psychiatric Genetics

Schizophrenia Polygenic Risk and Brain Structural Changes in Methamphetamine-Associated Psychosis in a South African Population

  • 1SA MRC Unit on Risk & Resilience in Mental Disorders, Department of Psychiatry and Neuroscience Institute, University of Cape Town, Cape Town, South Africa
  • 2Department of Child and Adolescent Psychiatry and Psychotherapy, TU Dresden, Dresden, Germany
  • 3The Stanley Center for Psychiatric Research, Broad Institute, Cambridge, MA, United States

Background: The genetic architecture of psychotic disorders is complex, with hundreds of genetic risk loci contributing to a polygenic model of disease. Overlap in the genetics of psychotic disorders and brain measures has been found in European populations, but has not been explored in populations of African ancestry. The aim of this study was to determine whether a relationship exists between a schizophrenia-derived PRS and (i) methamphetamine associated psychosis (MAP), and (ii) brain structural measures, in a South African population.

Methods: The study sample consisted of three participant groups: 31 individuals with MAP, 48 with apsychotic methamphetamine dependence, and 49 healthy controls. Using PRSice, PRS was generated for each of the participants with GWAS summary statistics from the Psychiatric Genomics Consortium Schizophrenia working group (PGC-SCZ2) as the discovery dataset. Regression analyses were performed to determine associations of PRS, with diagnosis, whole brain, and regional gray and white matter measures.

Results: Schizophrenia-derived PRS did not significantly predict MAP diagnosis. After correction for multiple testing, no significant associations were found between PRS and brain measures across all groups.

Discussion: The lack of significant associations here may indicate that the study is underpowered, that brain volumes in MAP are due to factors other than polygenic risk for schizophrenia, or that PRS derived from a largely European discovery set has limited utility in individuals of African ancestry. Larger studies, that include diverse populations, and more nuanced brain measures, may help elucidate the relationship between schizophrenia-PRS, brain structural changes, and psychosis.

Conclusion: This research presents the first PRS study to investigate shared genetic effects across psychotic disorders and brain structural measures in an African population. Ancestrally comparable discovery datasets may be useful for future African genetic research.

Introduction

Methamphetamine-associated psychosis (MAP) is the development of psychosis during, or soon after, intoxication or withdrawal from methamphetamine (MA). MA substance use disorder is a substantial local and global public health burden (Shin et al., 2017) with up to 40% of those who use MA going on to develop psychosis (Glasner-Edwards and Mooney, 2014). Genetic susceptibility has been recognized as a significant risk factor for the development of psychotic disorders such as MAP and schizophrenia (Chen et al., 2003; Grant et al., 2012). However, the genetic architecture of psychotic disorders is complex, having hundreds of risk loci contributing to the polygenic model of disease (Lvovs et al., 2012).

Polygenic risk scoring (PRS) is a method used to elucidate the polygenic nature of complex disorders by measuring the common variant contribution to the phenotype of interest. PRS has become an established method to determine genetic risk not only within, but also across disorders that share similar phenotypes (International Schizophrenia Consortium, 2009). With the use of PRS, evidence has emerged that there is shared genetic risk across psychotic diagnoses, including schizophrenia and bipolar disorder (International Schizophrenia Consortium, 2009; Hamshere et al., 2011; Tesli et al., 2014; Harrisberger et al., 2016; Vassos et al., 2017; Jonas et al., 2019).

Brain structural measures demonstrate heritability and are altered in psychotic disorders (Keshavan et al., 2007; Grant et al., 2012; Miller and Rockstroh, 2016; Jia and Ck, 2018). More recently, significant associations between schizophrenia-derived PRS and variation in brain structural measures, including global cortical thickness and hippocampal volume have been demonstrated in healthy individuals (Jalbrzikowski et al., 2019; Neilson et al., 2019).

To the knowledge of the authors, there have not been any studies using PRS to investigate shared genetic effects across psychotic disorders and brain structural measures in African populations. The majority of neuropsychiatric research, including PRS studies, has investigated European populations (Duncan et al., 2019). Using genetic and imaging data, the aim of this study was to determine whether a relationship exists between a schizophrenia-derived PRS and (i) MAP diagnosis, and (ii) brain structural measures, in a South African population.

Materials and Methods

Study Participants

Data for this study were from the case-control study, Neural correlates of deficits in affect regulation in methamphetamine dependence with and without a history of psychosis (Uhlmann, 2015). Ethical approval was obtained from the University of Cape Town Human Research Ethics Committee (684/2017). This study comprised individuals with a diagnosis of MAP (n = 31), methamphetamine dependence without psychosis (MD) (n = 48), and healthy controls (HC) (n = 49). MAP and MD diagnoses were made using the Structured Clinical Interview for DSM-IV Axis I disorders (First et al., 2016). Participants were excluded if they had a history of head trauma, other psychiatric or neurological disease, additional substance dependence (other than nicotine), or a seropositive test for HIV. The participants were of different ancestry, including African, European and mixed ancestry. Participants were matched for age and gender.

Genotyping and Quality Control

Blood or saliva samples (using the Oragene DNA OG-500 kits) were obtained from each of the participants with the appropriate informed consent. DNA was extracted from blood using the salting out method (Miller et al., 1988) and where saliva samples were collected, DNA was extracted using the manufacturer’s guidelines1. Samples were genotyped using the Illumina® Infinium PsychArray at the Broad Institute (Cambridge, MA, United States). Genotyping data was available for 588,454 variants. Using Plink v1.9 (Chang et al., 2015; Purcell and Chang, 2015), the following quality control (QC) steps were performed: removal of duplicate samples, relatedness check where individuals with pi-hat >0.2 were removed, Hardy-Weinberg Equilibrium (HWE) checks (p < 0.001), variants with a minor allele frequency (MAF) <0.05 were removed; variants with >10% missing genotype rate were removed.

Imputation and Post-imputation Quality Control

The Michigan Imputation Server (U.S. National Institutes of Health, 2020) was used to impute the genotype data. For this, the 1000 Genomes (1000G) phase3 v5 ref panel (Auton et al., 2015) was used and the rsq filter was set at 0.3 (estimate of the squared correlation between imputed and true genotypes). Phasing was performed using Eagle v2.4. Post-imputation QC consisted of Hardy-Weinberg Equilibrium (HWE) checks (p < 1e-6), variants with a minor allele frequency (MAF) <0.05 were removed; and variants with >10% missing genotype rate were removed. After QC, 128 individuals and 8,249,215 variants remained for downstream analysis.

Structural Brain Imaging

Structural MRI images were acquired using a 3T Siemens Magnetom Allegra at the Cape Universities Brain Imaging Centre. A radiologist, blinded to diagnosis, examined each scan for structural abnormalities. MRI scans were analyzed using the FreeSurfer software package v5.32 and images were quality controlled following the ENIGMA protocol3. Brain measures and regions of interest which have been previously found to be structurally or genetically associated with psychotic disorders, including MAP and schizophrenia, were chosen for analysis: (1) total brain volume; (2) total white matter (WM) volume; (3) cortical thickness; (4) hippocampal volume; (5) inferior temporal gyrus WM volume; and (6) superior temporal gyrus WM volume.

Principal Component Analysis (PCA)

To account for population stratification, principal components (PC) were calculated from a set of independent SNPs using flashPCA (Abraham and Inouye, 2014), with the 1000G dataset as a reference panel (Auton et al., 2015). The first two PCs were plotted using R (package ggplot2) (R Core Team, 2016; Wickham, 2016; Supplementary Figure 1). Prior to calculating the PCs, pairwise LD pruning was conducted with a window size of 1,000 variants and an r2 threshold of 0.05 using Plink v2.0 (Chang et al., 2015; Purcell and Chang, 2015).

PRS

PRS combines the effects of risk variants at specified p-value thresholds from a “discovery” GWAS, into a single risk score. The number of alleles an individual in the “target” dataset possesses for a particular variant, is weighted by the effect size of that variant in the “discovery” dataset (International Schizophrenia Consortium, 2009). To assess whether the aggregate scores reflect risk of disease, a mean risk score in target cases is compared to controls (International Schizophrenia Consortium, 2009). The discovery dataset was the Psychiatric Genetic Consortium Schizophrenia (PGC-SCZ2) GWAS summary statistics, comprising 102,636 SNPs4 (Ripke et al., 2014; Psychiatric Genomics Consortium, 2016). This large database is made up of 49 ancestry matched, non-overlapping case-control samples (46 of European and three of east Asian ancestry) with 34,241 cases and 45,604 controls; and 3 family based samples of European ancestry (1,235 parent affected-offspring trios). PRS were calculated in the target sample (n = 128) using PRSice-2 (Choi and O’reilly, 2019) at multiple p-value thresholds (PT) (0.001, 0.05, 0.1, 0.2, 0.3, 0.4, 0.5, 1). As a default in PRSice, LD pruning was performed whereby variants were pruned if they were nearby (within 250 kb) and in LD (r2 > 0.1) with the leading variant (lowest p-value) in a given region. Distribution of explained variance (R2) and optimal p-value threshold (pT) were calculated. To generate the best-fit PRS, MAP diagnosis was used as the primary outcome measure where the first five PCs, age and sex were included as covariates. The proportion of variance explained by PRS was estimated as the difference in Nagelkerke’s R2 between the full model (including PRS plus covariates) and the null model (only covariates).

Statistical Analysis

The Shapiro-Wilk test was used to determine whether the PRS were normally distributed and Levene’s test was used to determine whether there was homogeneity of variance across groups (Shapiro and Wilk, 1965; Schultz, 1985). A t-test was performed to test for difference in mean PRS between the MAP and the combined MD and HC groups (apsychotic group). A Welch one-way test was then performed to test for difference in mean PRS between all three groups, MAP, MD, and HC (Welch, 1947). Area under the receiver-operator characteristic curve (AUROC) was calculated to evaluate the ability of a schizophrenia-derived PRS to accurately predict MAP diagnosis in this population.

Linear regression was used to determine the association between best-fit PRS and brain structural measures, with the first five PCs, sex, age, and intracranial volume included as covariates. Participants with missing brain measures were removed. Brain measures were log transformed to improve the model fit and reduce the residual standard error. Significance was set at a Bonferroni-corrected p-value threshold of <0.005 (0.05/10 brain measures). Logistic regression was done to test association between MAP diagnosis and brain regional measures. All analyses were performed in the statistical environment, R (R Core Team, 2016).

Post-hoc Power Calculation

The “avengeme” R package (Dudbridge, 2013) was used to calculate the sample size necessary to achieve 80% power. With a disease prevalence of 1%, 40,755 overlapping independent SNPs between the discovery and target datasets, and n = 128 for the target dataset, a total of 98,885 samples are required in the discovery dataset to attain this level of power. However, to achieve 80% power with the specified target sample size, a p-value threshold of 0.99, for selecting markers into the PRS, is required.

Results

Sample Demographics

The demographics of the study sample are listed in Table 1. The majority of the sample was male (n = 100) with a mean age of 27 (6) years. The ancestries of the target sample were 97 (75.8%) mixed ancestry, 21 (16.4%) African, and 10 (7.8%) European.

TABLE 1
www.frontiersin.org

Table 1. Participant demographics.

PRS Nagelkerke R2 and AUC

The best fit PRS at a PT = 0.0099 (explaining 4.2% of the variance in MAP, p = 0.05) was used for downstream analysis (Supplementary Figure 2). PRS showed normal distribution across groups. There was no significant difference in mean PRS between MAP and AP participant groups [t = −1.5086, df = 126, p = 0.1339, 95% CI (−2.30 × 10–4, 3.10 × 10–5)]. There was no significant difference in means when groups were further split into MAP, MD and HC [F(2, 70) = 2.49, p = 0.09]. Using the PRS as the predictor, and MAP participant group as the outcome, the corrected partial area under the curve was 53.4% (95% CI: 49.7–61.6%).

Brain Measures

A total of 18 samples were removed due to missing brain measures, leaving 27 MAP patients and 83 AP patients in the regression analyses. After correction for multiple testing, there were no statistically significant association between PRS and any of the log transformed brain measures across all groups. These results are summarized in Table 2. The strongest associations were observed in left inferior temporal WM volume (β = −9.82 × 102, p = 0.02) and left superior temporal WM volume (β = −1.12 × 103, p = 0.009). No significant associations were found between any of the brain measures and MAP diagnosis. The differences in brain measures between groups was not the focus of this research, but has been published elsewhere (Uhlmann et al., 2016).

TABLE 2
www.frontiersin.org

Table 2. Linear regression results for brain structural measures and polygenic risk score.

Discussion

In this study, we sought to determine whether PRS derived from a well-powered schizophrenia GWAS, comprising mainly European populations, was able to predict MAP and brain volume and thickness in an ancestrally diverse South African target sample. After correction for multiple testing, we did not identify a significant association between schizophrenia-derived PRS and MAP diagnosis, or any of the brain structural measures. This is in contrast to a study conducted in an Asian population group which showed that a large number of “risk” alleles for MAP were over-represented in individuals with schizophrenia. However, it is worth noting that the overlap from this previous study was only able to explain 0.7% of the variance in schizophrenia liability (Ikeda et al., 2013).

There are a number of different explanations for the null findings obtained here. In particular, polygenic risk scores for schizophrenia may not be associated with measures of brain volume in healthy individuals or in those with psychotic disorders. This explanation is consistent with findings from previous research in populations of largely European ancestry (Reus et al., 2017; Harrisberger et al., 2018; Lancaster et al., 2018; Simões et al., 2020). For example, no associations between PRS for schizophrenia and for bipolar disorder with either subcortical volume or WM microstructure, were found in the United Kingdom Biobank (Reus et al., 2017). Similarly, in healthy subjects with higher genetic risk for schizophrenia, based on loci found to be associated with schizophrenia (Psychiatric Genomics Consortium, 2016), no consistent associated brain volume changes were observed (Van Der Auwera et al., 2017). Furthermore, a systematic review established that schizophrenia-derived PRSs were not significantly associated with brain structural changes in five out of the seven studies included (Van Der Merwe et al., 2018). As the knowledge and evidence of imaging genetics increases, more nuanced structural and functional brain measures are being studied. For example, measures of brain connectivity and of task-dependent recruitment of multiple brain regions may be associated with genetic variation in psychosis (Ranlund et al., 2017; Dezhina et al., 2019; Cao et al., 2020). Added to this, advanced methodological approaches integrating environmental exposures, gene-gene interactions, and epigenetics from a variety of ancestral cohorts, are needed to fully appreciate missing heritability (Mufford et al., 2017; Alnæs et al., 2019).

Caution is also needed when interpreting studies with differing ancestries in the discovery and target datasets. Correlations between true (discovery population) and inferred (target population) risk are highest in the population from which summary statistics are derived (Martin et al., 2017; Mostafavi et al., 2019). Therefore, scores are fundamentally less informative in populations more diverged from the discovery GWAS study cohorts (Scutari et al., 2016; Martin et al., 2019). The majority of neuropsychiatric genetic studies have been undertaken in high-income settings, thus GWAS summary statistics are limited to mostly individuals of European and, more recently, Asian ancestry (Martin et al., 2017; Duncan et al., 2019). A review of PRS performance in diverse human populations revealed that, with current available genetic discovery datasets, people of African descent have the lowest polygenic score performance, compared to all other populations tested, including Latino, Middle Eastern, and East and South Asian. Indeed, predictive performance of European ancestry-derived PRS in populations of African descent is only 42% of that of matched European ancestry samples- almost halving the potential of studies such as our own (Martin et al., 2017; Vassos et al., 2017; Duncan et al., 2019).

This points to what is needed next: targeted large-scale genetic investigation of schizophrenia and other psychiatric disorders in African populations. Such research might ensure that future developments, and the clinical utility, of PRS will be equally applicable to health care users of African descent and limit exacerbating already existing health disparities (Martin et al., 2019). Projects under way, such as the Neuropsychiatric Genetics of African Populations-Psychosis (NeuroGAP-Psychosis), are expanding knowledge of the genetic and environmental risk architecture of neuropsychiatric disorders in African populations of South Africa, Ethiopia, Kenya and Uganda (Stevenson et al., 2019). This will improve the availability of ancestrally comparable discovery datasets and meaningful results for African genetic research. In further support of this call to action, it has been recognized that African genomes contain more diversity than any other population group, having uniquely shorter haplotype blocks with more variants per individual (Genomes Project Consortium, 2015). As a result, studies of African populations could uncover additional pathogenic variants and identify novel disease-associated loci (Dalvie et al., 2015). Therefore, future African neuropsychiatric research may yield genomic insights into the risk, resilience (Wojcik et al., 2019), and treatment of psychiatric disorders, advancing precision medicine across global populations (Dalvie et al., 2015).

This study has several limitations that deserve emphasis. Firstly, the target group was relatively small, and comprised participants of different ancestral groups. Although the use of PCAs as covariates in the PRS is a robust method for correcting for population stratification (Wu et al., 2011), it cannot resolve all bias due to differences in the population sample (Price et al., 2006). Ideally, the study sample should be stratified into ancestrally homogenous groups and analyzed separately (Duncan et al., 2019). Secondly, psychiatric diagnosis was made on the basis of a single diagnostic interview, which may lead to misclassification (Schijven et al., 2020). Future research on MAP may benefit from longitudinal study designs which ensure longitudinal expert assessment and diagnosis. Lastly, the brain structural changes observed in previous psychiatric research have been subtle, are highly heterogeneous, and are influenced by multiple potential confounding factors (Brent et al., 2013; Weinberger and Radulescu, 2015). Such confounding effects include: demographics, frequency of MA use, long term neuroleptic treatment, participant movement, hydration and stress during scan, which were not accounted for in this study (Yudofsky and Hales, 2004; Streitbürger et al., 2012; Weinberger and Radulescu, 2015; Mufford et al., 2017; Arunogiri et al., 2018).

Conclusion

This research is the first to use PRS to investigate shared genetic effects across psychotic disorders and brain structural measures in an African population. Ancestrally comparable datasets and more nuanced structural and functional brain measures may be useful in further elucidating the genetic risks for psychotic disorders including MAP.

Data Availability Statement

The datasets generated for this article are not publicly available The primary data was collected using a patient consent form that did not include consent to deposit genotyped data into a public repository. Requests to access the datasets should be directed to AU, dWhsbWFubi5hYUBnbWFpbC5jb20=.

Ethics Statement

This study was reviewed and approved by the University of Cape Town Human Research Ethics Committee (approved 684/2017). The participants provided written informed consent.

Author Contributions

All authors contributed to the article and approved the submitted version. RP: study design, statistical and genetic analysis, draft of manuscript, and revisions. DS: study design, review of manuscript, and editing. AU: collection of primary data and review of manuscript. CM: review of manuscript. SD: study design, statistical and genetic analysis, review and revisions of manuscript.

Funding

This study was funded by SA MRC Unit on Risk & Resilience in Mental Disorders.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

Computations were performed using facilities provided by the University of Cape Town’s ICTS High Performance Computing team: hpc.uct.ac.za.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2020.01018/full#supplementary-material

Footnotes

  1. ^ https://www.dnagenotek.com/us/pdf/PD-PR-006.pdf
  2. ^ http://surfer.nmr.mgh.harvard.edu/
  3. ^ http://enigma.ini.usc.edu/protocols/imaging-protocols
  4. ^ http://pgc.unc.edu

References

Abraham, G., and Inouye, M. (2014). Fast principal component analysis of large-scale genome-wide data. PLoS One 9:e93766. doi: 10.1371/journal.pone.0093766

PubMed Abstract | CrossRef Full Text | Google Scholar

Alnæs, D., Kaufmann, T., Van Der Meer, D., Córdova-Palomera, A., Rokicki, J., Moberget, T., et al. (2019). Brain heterogeneity in schizophrenia and its association with polygenic risk. JAMA Psychiatry 76, 739–748. doi: 10.1001/jamapsychiatry.2019.0257

PubMed Abstract | CrossRef Full Text | Google Scholar

Arunogiri, S., Foulds, J. A., Mcketin, R., and Lubman, D. I. (2018). A systematic review of risk factors for methamphetamine-associated psychosis. Aust. N. Z. J. Psychiatry 52, 514–529. doi: 10.1177/0004867417748750

PubMed Abstract | CrossRef Full Text | Google Scholar

Auton, A., Abecasis, G. R., Altshuler, D. M., Durbin, R. M., Abecasis, G. R., Bentley, D. R., et al. (2015). A global reference for human genetic variation. Nature 526, 68–74. doi: 10.1038/nature15393

PubMed Abstract | CrossRef Full Text | Google Scholar

Brent, B. K., Thermenos, H. W., Keshavan, M. S., and Seidman, L. J. (2013). Gray matter alterations in schizophrenia high-risk youth and early-onset schizophrenia: a review of structural MRI findings. Child Adolesc. Psychiatr. Clin. 22, 689–714. doi: 10.1016/j.chc.2013.06.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Cao, H., Zhou, H., and Cannon, T. D. (2020). Functional connectome-wide associations of schizophrenia polygenic risk. Mol. Psychiatry 25, 1–9. doi: 10.1038/s41380-020-0699-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Chang, C. C., Chow, C. C., Tellier, L. C., Vattikuti, S., Purcell, S. M., and Lee, J. J. (2015). Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience 4:7. doi: 10.1186/s13742-015-0047-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, C.-K., Lin, S.-K., Sham, P., Ball, D., Loh, E.-W., Hsiao, C.-C., et al. (2003). Pre-morbid characteristics and co-morbidity of methamphetamine users with and without psychosis. Psychol. Med. 33, 1407–1414. doi: 10.1017/s0033291703008353

PubMed Abstract | CrossRef Full Text | Google Scholar

Choi, S. W., and O’reilly, P. F. (2019). PRSice-2: polygenic risk score software for biobank-scale data. Gigascience 8:giz082. doi: 10.1093/gigascience/giz082

PubMed Abstract | CrossRef Full Text | Google Scholar

Dalvie, S., Koen, N., Duncan, L., Abbo, C., Akena, D., Atwoli, L., et al. (2015). Large scale genetic research on neuropsychiatric disorders in African populations is needed. EBioMedicine 2, 1259–1261. doi: 10.1016/j.ebiom.2015.10.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Dezhina, Z., Ranlund, S., Kyriakopoulos, M., Williams, S. C. R., and Dima, D. (2019). A systematic review of associations between functional MRI activity and polygenic risk for schizophrenia and bipolar disorder. Brain Imaging Behav. 13, 862–877. doi: 10.1007/s11682-018-9879-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Dudbridge, F. (2013). Power and predictive accuracy of polygenic risk scores. PLoS Genet. 9:e1003348. doi: 10.1371/journal.pgen.1003348

PubMed Abstract | CrossRef Full Text | Google Scholar

Duncan, L., Shen, H., Gelaye, B., Meijsen, J., Ressler, K., Feldman, M., et al. (2019). Analysis of polygenic risk score usage and performance in diverse human populations. Nat. Commun. 10:3328. doi: 10.1038/s41467-019-11112-0

PubMed Abstract | CrossRef Full Text | Google Scholar

First, M. B., Williams, J. B., Karg, R. S., and Spitzer, R. L. (2016). Scid-5-Cv: Structured Clinical Interview for DSM-5 Disorders, Clinician Version. Washington, DC: American Psychiatric Association.

Google Scholar

Genomes Project Consortium (2015). A global reference for human genetic variation. Nature 526, 68–74. doi: 10.1038/nature15393

PubMed Abstract | CrossRef Full Text | Google Scholar

Glasner-Edwards, S., and Mooney, L. J. (2014). Methamphetamine psychosis: epidemiology and management. CNS Drugs 28, 1115–1126. doi: 10.1007/s40263-014-0209-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Grant, K. M., Levan, T. D., Wells, S. M., Li, M., Stoltenberg, S. F., Gendelman, H. E., et al. (2012). Methamphetamine-associated psychosis. J. Neuroimmune Pharmacol. 7, 113–139. doi: 10.1007/s11481-011-9288-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Hamshere, M. L., O’donovan, M. C., Jones, I. R., Jones, L., Kirov, G., Green, E. K., et al. (2011). Polygenic dissection of the bipolar phenotype. Br. J. Psychiatry 198, 284–288. doi: 10.1192/bjp.bp.110.087866

PubMed Abstract | CrossRef Full Text | Google Scholar

Harrisberger, F., Smieskova, R., Egli, T., Simon, A. E., Riecher-Rössler, A., Fusar-Poli, P., et al. (2018). Impact on the onset of psychosis of a polygenic schizophrenia-related risk score and changes in white matter volume. Cell. Physiol. Biochem. 48, 1201–1214. doi: 10.1159/000491986

PubMed Abstract | CrossRef Full Text | Google Scholar

Harrisberger, F., Smieskova, R., Vogler, C., Egli, T., Schmidt, A., Lenz, C., et al. (2016). Impact of polygenic schizophrenia-related risk and hippocampal volumes on the onset of psychosis. Transl. Psychiatry 6:e868. doi: 10.1038/tp.2016.143

PubMed Abstract | CrossRef Full Text | Google Scholar

Ikeda, M., Oahisa, Y., Aleksic, B., Won, M., Kondo, N., Naruse, N., et al. (2013). Evidence for shared genetic risk between methamphetamine-induced psychosis and schizophrenia. Neuropsychopharmacology 38, 1864–1870. doi: 10.1038/npp.2013.94

PubMed Abstract | CrossRef Full Text | Google Scholar

International Schizophrenia Consortium (2009). Common polygenic variation contributes to risk of schizophrenia and bipolar disorder. Nature 460, 748–752. doi: 10.1038/nature08185

PubMed Abstract | CrossRef Full Text | Google Scholar

Jalbrzikowski, M., Klei, L., and Devlin, B. (2019). T20. Polygenic risk for schizophrenia is associated with hippocampal volume in typically developing youth. Schizophr. Bull. 45:S211. doi: 10.1093/schbul/sbz019.300

PubMed Abstract | CrossRef Full Text | Google Scholar

Jia, H., and Ck, C. R. (2018). Progress in the studies on endophenotypes of schizophrenia. Chin. Sci. Bull. 63, 127–135. doi: 10.1360/n972015-01383

CrossRef Full Text | Google Scholar

Jonas, K. G., Lencz, T., Li, K., Malhotra, A. K., Perlman, G., Fochtmann, L. J., et al. (2019). Schizophrenia polygenic risk score and 20-year course of illness in psychotic disorders. Transl. Psychiatry 9:300. doi: 10.1038/s41398-019-0612-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Keshavan, M. S., Prasad, K. M., and Pearlson, G. (2007). Are brain structural abnormalities useful as endophenotypes in schizophrenia? Int. Rev. Psychiatry 19, 397–406. doi: 10.1080/09540260701486233

PubMed Abstract | CrossRef Full Text | Google Scholar

Lancaster, T. M., Dimitriadis, S. L., Tansey, K. E., Perry, G., Ihssen, N., Jones, D. K., et al. (2018). Structural and functional neuroimaging of polygenic risk for schizophrenia: a recall-by-genotype–based approach. Schizophr. Bull. 45, 405–414. doi: 10.1093/schbul/sby037

PubMed Abstract | CrossRef Full Text | Google Scholar

Lvovs, D., Favorova, O., and Favorov, A. (2012). A polygenic approach to the study of polygenic diseases. Acta Naturae 4, 59–71. doi: 10.32607/20758251-2012-4-3-59-71

CrossRef Full Text | Google Scholar

Martin, A. R., Gignoux, C. R., Walters, R. K., Wojcik, G. L., Neale, B. M., Gravel, S., et al. (2017). Human demographic history impacts genetic risk prediction across diverse populations. Am. J. Hum. Genet. 100, 635–649. doi: 10.1016/j.ajhg.2017.03.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Martin, A. R., Kanai, M., Kamatani, Y., Okada, Y., Neale, B. M., and Daly, M. J. (2019). Clinical use of current polygenic risk scores may exacerbate health disparities. Nat. Genet. 51, 584–591. doi: 10.1038/s41588-019-0379-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Miller, G. A., and Rockstroh, B. S. (2016). “Chapter 2 - progress and prospects for endophenotypes for schizophrenia in the time of genomics, epigenetics, oscillatory brain dynamics, and the research domain criteria,” in The Neurobiology of Schizophrenia, eds T. Abel, and T. Nickl-Jockschat (San Diego, CA: Academic Press). doi: 10.1016/B978-0-12-801829-3.00010-0

CrossRef Full Text | Google Scholar

Miller, S., Dykes, D., and Polesky, H. (1988). A simple salting out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res. 16:1215. doi: 10.1093/nar/16.3.1215

PubMed Abstract | CrossRef Full Text | Google Scholar

Mostafavi, H., Harpak, A., Conley, D., Pritchard, J. K., and Przeworski, M. (2019). Variable prediction accuracy of polygenic scores within an ancestry group. bioRxiv:629949. doi: 10.1101/629949

CrossRef Full Text | Google Scholar

Mufford, M. S., Stein, D. J., Dalvie, S., Groenewold, N. A., Thompson, P. M., and Jahanshad, N. (2017). Neuroimaging genomics in psychiatry—a translational approach. Genome Med. 9:102. doi: 10.1186/s13073-017-0496-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Neilson, E., Shen, X., Cox, S. R., Clarke, T.-K., Wigmore, E. M., Gibson, J., et al. (2019). Impact of polygenic risk for schizophrenia on cortical structure in UK biobank. Biol. Psychiatry 86, 536–544. doi: 10.1016/j.biopsych.2019.04.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Price, A. L., Patterson, N. J., Plenge, R. M., Weinblatt, M. E., Shadick, N. A., and Reich, D. (2006). Principal components analysis corrects for stratification in genome-wide association studies. Nat. Genet. 38, 904–909. doi: 10.1038/ng1847

PubMed Abstract | CrossRef Full Text | Google Scholar

Psychiatric Genomics Consortium (2016). University of North Carolina at Chapel Hill School of Medicine. Chapel Hill: UNC School of Medicine.

Google Scholar

Purcell, S., and Chang, C. (2015). PLINK 2.0 [Online]. Available online at: www.cog-genomics.org/plink/2.0/ (accessed August, 2020).

Google Scholar

R Core Team (2016). A Language and Environment for Statistical Computing. Version 4.0.2. Vienna, Austria: R Foundation for Statistical Computing. Available online at: https://www.R-project.org/ (accessed August, 2020).

Google Scholar

Ranlund, S., Calafato, S., Thygesen, J. H., Lin, K., Cahn, W., Crespo-Facorro, B., et al. (2017). A polygenic risk score analysis of psychosis endophenotypes across brain functional, structural, and cognitive domains. Am. J. Med. Genet. Part B Neuropsychiatr. Genet. 177, 21–34. doi: 10.1002/ajmg.b.32581

PubMed Abstract | CrossRef Full Text | Google Scholar

Reus, L. M., Shen, X., Gibson, J., Wigmore, E., Ligthart, L., Adams, M. J., et al. (2017). Association of polygenic risk for major psychiatric illness with subcortical volumes and white matter integrity in UK Biobank. Sci. Rep. 7:42140. doi: 10.1038/srep42140

PubMed Abstract | CrossRef Full Text | Google Scholar

Ripke, S., Neale, B. M., Corvin, A., Walters, J. T., Farh, K.-H., Holmans, P. A., et al. (2014). Biological insights from 108 schizophrenia-associated genetic loci. Nature 511, 421–427. doi: 10.1038/nature13595

PubMed Abstract | CrossRef Full Text | Google Scholar

Schijven, D., Veldink, J. H., and Luykx, J. J. (2020). Genetic cross-disorder analysis in psychiatry: from methodology to clinical utility. Br. J. Psychiatry 216, 246–249. doi: 10.1192/bjp.2019.72

PubMed Abstract | CrossRef Full Text | Google Scholar

Schultz, B. B. (1985). Levene’s test for relative variation. Syst. Zool. 34, 449–456. doi: 10.2307/2413207

CrossRef Full Text | Google Scholar

Scutari, M., Mackay, I., and Balding, D. (2016). Using genetic distance to infer the accuracy of genomic prediction. PLoS Genet. 12:e1006288. doi: 10.1371/journal.pgen.1006288

PubMed Abstract | CrossRef Full Text | Google Scholar

Shapiro, S. S., and Wilk, M. B. (1965). An analysis of variance test for normality (complete samples). Biometrika 52, 591–611. doi: 10.1093/biomet/52.3-4.591

CrossRef Full Text | Google Scholar

Shin, E.-J., Dang, D.-K., Tran, H.-Q., Jeong, J. H., Nah, S.-Y., Jang, C.-G., et al. (2017). Current understanding of methamphetamine-associated dopaminergic neurodegeneration and psychotoxic behaviors. Arch. Pharm. Res. 40, 403–428. doi: 10.1007/s12272-017-0897-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Simões, B., Vassos, E., Shergill, S., Mcdonald, C., Toulopoulou, T., Kalidindi, S., et al. (2020). Schizophrenia polygenic risk score influence on white matter microstructure. J. Psychiatr. Res. 121, 62–67. doi: 10.1016/j.jpsychires.2019.11.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Stevenson, A., Akena, D., Stroud, R. E., Atwoli, L., Campbell, M. M., Chibnik, L. B., et al. (2019). Neuropsychiatric genetics of African populations-psychosis (neurogap-psychosis): a case-control study protocol and GWAS in Ethiopia, Kenya, South Africa and Uganda. BMJ Open 9:e025469. doi: 10.1136/bmjopen-2018-025469

PubMed Abstract | CrossRef Full Text | Google Scholar

Streitbürger, D.-P., Möller, H. E., Tittgemeyer, M., Hund-Georgiadis, M., Schroeter, M. L., and Mueller, K. (2012). Investigating structural brain changes of dehydration using voxel-based morphometry. PLoS One 7:e44195. doi: 10.1371/journal.pone.0044195

PubMed Abstract | CrossRef Full Text | Google Scholar

Tesli, M., Espeseth, T., Bettella, F., Mattingsdal, M., Aas, M., Melle, I., et al. (2014). Polygenic risk score and the psychosis continuum model. Acta Psychiatr. Scand. 130, 311–317. doi: 10.1111/acps.12307

PubMed Abstract | CrossRef Full Text | Google Scholar

U.S. National Institutes of Health (2020). Michigan Imputation Server. Bethesda, MD: U.S. National Institutes of Health.

Google Scholar

Uhlmann, A. (2015). Neural Correlates of Deficits in Affect Regulation in Methamphetamine Dependence with and without a History of Psychosis. Doctoral dissertation, University of Cape Town, Cape Town.

Google Scholar

Uhlmann, A., Fouche, J. P., Koen, N., Meintjes, E. M., Wilson, D., and Stein, D. J. (2016). Fronto-temporal alterations and affect regulation in methamphetamine dependence with and without a history of psychosis. Psychiatry Res. Neuroimaging 248, 30–38. doi: 10.1016/j.pscychresns.2016.01.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Der Auwera, S., Wittfeld, K., Shumskaya, E., Bralten, J., Zwiers, M. P., Onnink, A. M. H., et al. (2017). Predicting brain structure in population-based samples with biologically informed genetic scores for schizophrenia. Am. J. Med. Genet. Part B Neuropsychiatr. Genet. 174, 324–332. doi: 10.1002/ajmg.b.32519

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Der Merwe, C., Passchier, R., Mufford, M., Ramesar, R., Dalvie, S., and Stein, D. (2018). Polygenic risk for schizophrenia and associated brain structural changes: a systematic review. Compr. Psychiatry 88, 77–82. doi: 10.1016/j.comppsych.2018.11.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Vassos, E., Di Forti, M., Coleman, J., Iyegbe, C., Prata, D., Euesden, J., et al. (2017). An examination of polygenic score risk prediction in individuals with first-episode psychosis. Biol. Psychiatry 81, 470–477. doi: 10.1016/j.biopsych.2016.06.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Weinberger, D. R., and Radulescu, E. (2015). Finding the elusive psychiatric “lesion” with 21st-century neuroanatomy: a note of caution. Am. J. Psychiatry 173, 27–33. doi: 10.1176/appi.ajp.2015.15060753

PubMed Abstract | CrossRef Full Text | Google Scholar

Welch, B. L. (1947). The generalization ofstudent’s’ problem when several different population variances are involved. Biometrika 34, 28–35. doi: 10.2307/2332510

CrossRef Full Text | Google Scholar

Wickham, H. (2016). ggplot2: Elegant Graphics for Data Analysis. Berlin: Springer. doi: 10.1007/978-3-319-24277-4

CrossRef Full Text | Google Scholar

Wojcik, G. L., Graff, M., Nishimura, K. K., Tao, R., Haessler, J., Gignoux, C. R., et al. (2019). Genetic analyses of diverse populations improves discovery for complex traits. Nature 570, 514–518. doi: 10.1038/s41586-019-1310-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, C., Dewan, A., Hoh, J., and Wang, Z. (2011). A comparison of association methods correcting for population stratification in case–control studies. Ann. Hum. Genet. 75, 418–427. doi: 10.1111/j.1469-1809.2010.00639.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Yudofsky, S. C., and Hales, R. E. (2004). Essentials of Neuropsychiatry and Clinical Neurosciences. Washington, DC: American Psychiatric Publishing, Inc.

Google Scholar

Keywords: polygenic risk, methamphetamine associated psychosis, Africa, brain measures, schizophrenia

Citation: Passchier RV, Stein DJ, Uhlmann A, van der Merwe C and Dalvie S (2020) Schizophrenia Polygenic Risk and Brain Structural Changes in Methamphetamine-Associated Psychosis in a South African Population. Front. Genet. 11:1018. doi: 10.3389/fgene.2020.01018

Received: 28 November 2019; Accepted: 10 August 2020;
Published: 02 October 2020.

Edited by:

Francis J. McMahon, National Institutes of Health (NIH), United States

Reviewed by:

Donald Lyall, University of Glasgow, United Kingdom
Gabriel R. Fries, The University of Texas Health Science Center at Houston, United States

Copyright © 2020 Passchier, Stein, Uhlmann, van der Merwe and Dalvie. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Ruth V. Passchier, cnV0aHBhc3NjaGllckBnbWFpbC5jb20=

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.