Follow-up of a major psychosis linkage site in 13q13-q14 reveals significant association in both case-control and family samples

Heterogeneity in the longitudinal courses of global functioning in children at familial risk of major psychiatric disorders: Association with trauma and familial characteristics

OBJECTIVES

The extent to which heterogeneity in childhood risk trajectories may underlie later heterogeneity in schizophrenia (SZ), bipolar disorder (BP), and major depressive disorder (MDD) remains a chief question. Answers may optimally be found by studying the longitudinal trajectories of children born to an affected parent. We aimed to differentiate trajectories of global functioning and their sensitive periods from the age of 6 to 17 years in children at familial risk (FHRs).

METHODS

First, a latent class mixed model analysis (LCMM) was applied to yearly ratings of the Children's Global Assessment Scale (CGAS) from the age of 6 to 17 years in 170 FHRs born to a parent affected by DSM-IV SZ (N = 37), BP (N = 82) or MDD (N = 51). Then, we compared the obtained Classes or trajectories of FHRs in terms of sex, parental diagnosis, IQ, child clinical status, childhood trauma, polygenic risk score (PRS), and outcome in transition to illness.

RESULTS

The LCMM on yearly CGAS trajectories identified a 4-class solution showing markedly different childhood and adolescence dynamic courses and temporal vulnerability windows marked by a functioning decline and a degree of specificity in parental diagnosis. Moreover, IQ, trauma exposure, PRS level, and timing of later transition to illness differentiated the trajectories. Almost half (46%) of the FHRs exhibited a good and stable global functioning trajectory.

CONCLUSIONS

FHRs of major psychiatric disorders show heterogeneous functional decline during development associated with parental diagnosis, polygenic risk loading, and childhood trauma.

Proof-of-concept study of a multi-gene risk score in adolescent bipolar disorder

BACKGROUND

Few studies have examined multiple genetic variants concurrently for the purpose of classifying bipolar disorder (BD); the literature among youth is particularly sparse. We selected 35 genetic variants, previously implicated in BD or associated characteristics, from which to identify the most robustly predictive group of genes.

METHODS

215 Caucasian adolescents (114 BD and 101 healthy controls (HC), ages 13-20 years) were included. Psychiatric diagnoses were determined based on semi-structured diagnostic interviews. Genomic DNA was extracted from saliva for genotyping. Two models were used to calculate a multi-gene risk score (MGRS). Model 1 used forward and backward regressions, and model 2 used a PLINK generated method.

RESULTS

In model 1, GPX3 rs3792797 was significant in the forward regression, DRD4 exonIII was significant in the backward regression; IL1β rs16944 and DISC1 rs821577 were significant in both the forward and backward regressions. These variants are involved in dopamine neurotransmission; inflammation and oxidative stress; and neuronal development. Model 1 MGRS did not significantly discriminate between BD and HC. In model 2, ZNF804A rs1344706 was significantly associated with BD; however, this association did not predict diagnosis when entered into the weighted model.

LIMITATIONS

This study was limited by the number of genetic variants examined and the modest sample size.

CONCLUSIONS

Whereas regression approaches identified four genetic variants that significantly discriminated between BD and HC, those same variants no longer discriminated between BD and HC when computed as a MGRS. Future larger studies are needed evaluating intermediate phenotypes such as neuroimaging and blood-based biomarkers.

A multimodal attempt to follow-up linkage regions using RNA expression, SNPs and CpG methylation in schizophrenia and bipolar disorder kindreds

The complexity of schizophrenia (SZ) and bipolar disorder (BD) has slowed down progress in understanding their genetic roots. Alternative genomic approaches are needed to bypass these difficulties. We attempted a multimodal approach to follow-up on reported linkage findings in SZ and BD from the Eastern Quebec kindreds in chromosomes 3q21, 4p34, 6p22, 8p21, 8p11, 13q11-q14, 15q13, 16p12, and 18q21. First, in 498 subjects, we measured RNA expression (47 K Illumina chips) in SZ and BD patients that we compared with their non-affected relatives (NARs) to identify, for each chromosomal region, genes showing the most significant differences in expression. Second, we performed SNP genotyping (700 K Illumina chips) and cis-eQTN analysis. Third, we measured DNA methylation on genes with RNA expression differences or eQTNs. We found a significant overexpression of the gene ITGB5 at 3q25 in SZ and BD after multiple testing p value adjustment. SPCS3 gene at 4q34, and FZD3 gene at 8p21, contained significant eQTNs after multiple testing corrections, while ITGB5 provided suggestive results. Methylation in associated genes did not explain the expression differences between patients and NARs. Our multimodal approach involving RNA expression, dense SNP genotyping and eQTN analyses, restricted to chromosomal regions having shown linkage, lowered the multiple testing burden and allowed for a deeper examination of candidate genes in SZ or BD.

Polygenic risk scores distinguish patients from non-affected adult relatives and from normal controls in schizophrenia and bipolar disorder multi-affected kindreds

Recent studies have used results on SNP association with schizophrenia (SZ) and bipolar disorder (BD) to create polygenic risk scores (PRS) discriminating non-familial unrelated patients from controls. Little is known about the role of PRS in densely affected multigenerational families. We tested PRS differences between affected SZ and BD family members from their non-affected adult relatives (NAARs) in Eastern Quebec Kindreds and from controls. We examined 1227 subjects: from 17 SZ and BD kindreds, we studied 153 patients (57 SZ, 13 schizoaffective, and 83 BD) and 180 NAARs, and 894 unrelated controls from the Eastern Quebec population. PRS were derived from published case-control association studies of SZ and BD. We also constructed a combined SZ and BD PRS by using SNPs from both SZ and BD PRS. SZ patients had higher SZ PRS than controls (p = 0.0039, R²  = 0.027) and BD patients had higher BD PRS than controls (p = 0.013, R²  = 0.027). Differences between affected subjects and NAARs and controls were significant with both SZ and BD PRS. Moreover, a combined SZ-BD PRS was also significantly associated with SZ and BD when compared to NAARs (p = 0.0019, R²  = 0.010) and controls (p = 0.0025, R²  = 0.028), revealing a SZ-BD commonality effect in PRS at the diagnosis level. The SZ and the BD PRS, however, showed a degree of specificity regarding thought disorder symptoms. Overall, our report would confirm the usefulness of PRS in capturing the contribution of common genetic variants to the risk of SZ and BD in densely affected families.

The longevity gene Klotho is differentially associated with cognition in subtypes of schizophrenia

Cognitive impairment is a core feature of schizophrenia and impacts negatively the functioning of affected individuals. Cognitive decline correlates with aging, and is the primary cause of loss of independence and reduced quality of life. The klotho gene is a key modulator of aging, with expression deficiency resulting in premature aging, while overexpression extends lifespan and enhances cognition. A haplotype and functional human variant of the gene, KL-VS, increases expression and promotes longevity. KL-VS heterozygosity is associated with enhanced cognition and a larger volume of the right dorsolateral prefrontal cortex, a region involved in planning and decision-making, which is especially susceptible to shrinkage with age. We examined the effect of KL-VS heterozygosity on cognition in 497 schizophrenia patients and 316 healthy controls from the Western Australian Family Study of Schizophrenia (WAFSS) who had been comprehensively characterised by neurocognitive tests and classified into cognitively deficient (CD) and cognitively "spared" (CS) clusters. An older, cognitively normal population sample from the Health in Men Study (HIMS) was included to allow assessment of heterozygosity and memory in aged individuals. We show that heterozygosity is associated with better learning and memory in the younger WAFSS healthy controls but not in the aging HIMS sample. However, in schizophrenia patients, KL-VS has a selective effect on memory, with heterozygotes in CD and CS clusters performing worse than non-carriers. This effect was significant and more severe in the CD cluster, reinforcing the utility of subtyping patients into CD and CS clusters that may differ in their genetic underpinnings.

The interaction of GSK3B and FXR1 genotypes may influence the mania and depression dimensions in mood disorders

BACKGROUND

Previous evidence in healthy subjects suggested that functional polymorphisms GSK3B rs12630592 and FXR1 rs496250 interact in regulating mood and emotional processing. We attempted to replicate this interaction primarily on manic and depressive dimensions in mood disorder patients, and secondarily on schizophrenia patients, diagnosis itself and age of onset.

METHODS

Symptom dimensions were derived from the Comprehensive Assessment of Symptoms and History 82 items rated lifetime in acute episodes and stabilized interepisode intervals in 384 patients from the Schizophrenia and Bipolar Disorder Eastern Quebec Kindred Study. Linear mixed effect models of symptom dimensions included rs12630592-rs496250 main and interaction fixed effects (obtained from TaqMan genotypes), and a polygenic random effect. The distribution of lifetime best-estimate DSM-IV diagnosis of 855 kindred members was studied versus genotype under a polytomous logistic model.

RESULTS

In mood disorder patients, the level of mania (in both acute and stabilized periods) and depression in stabilized periods was positively associated with GSK3B rs12630592 T only in FXR1 rs496250 A-allele carriers (Bonferroni-corrected interaction p=0.024, 0.052 and 0.017 respectively). The two polymorphisms explained 11% of mania variance and 5% of interepisode depression variance. The association was observed neither in schizophrenia patients nor with the psychotic dimension in mood disorder patients. Interaction with the diagnosis distribution (p=0.03) was driven by the decreasing prevalence of recurrent major depression with rs12630592 T also only in carriers of rs496250 A.

LIMITATIONS

Sample size was limited, but power was sufficient to detect the tested interaction effect in this replication sample.

CONCLUSIONS

We replicate in affective patients an interaction between the FXR1 rs496250 and GSK3B rs12630592 polymorphisms in regulating mood dimensions.

Polyunphased: an extension to polytomous outcomes of the Unphased package for family-based genetic association analysis

Polytomous phenotypes arise when a disease has multiple subtypes or when two dichotomous phenotypes are analyzed simultaneously. Few software programs offer the option to analyze such phenotypes in family studies, and none implements conditional polytomous logistic regression for within-family analysis robust to population stratification. We introduce Polyunphased, an extension to polytomous phenotypes of the Unphased package, a flexible software tool for genetic association analysis in nuclear families. Like Unphased, Polyunphased is written in C++ and runs from the command line or from a Java graphical user interface. Most Unphased options remain available in Polyunphased, including those handling missing parental genotypes while preserving robustness to population stratification, and the modelling options. Simulation studies confirmed the expected statistical behaviour of the maximum likelihood estimates of the association parameters of the conditional logistic regression model when the corresponding association parameters in the parental term of the likelihood function are set to 0, but revealed convergence problems when estimating these parental association parameters separately. The former approach is thus recommended with polytomous phenotypes.

Biological aspects and candidate biomarkers for psychotic bipolar disorder: A systematic review

AIM

We carried out a systematic review of the available literature about potential biomarkers of psychotic bipolar disorder (BD-P), a specific subset presenting worse outcome and greater risk of relapse than non-psychotic bipolar disorder (BD-NP).

METHODS

We searched the main psychiatric databases (PubMed, ISI Web of Knowledge, PsychInfo). Only original articles with the main topic of BD-P compared to schizophrenia/BD-NP/healthy controls (HC) written in English from 1994 to 2015 were included.

RESULTS

BD-P patients presented higher kynurenic acid levels in the cerebrospinal fluid, elevated anti- S accharomyces cerevisiae antibodies levels, and lower serum levels of dehydroepiandrosterone sulfate and progesterone than BD-NP/HC. Event-related potentials abnormalities have been identified in BD-P with respect to BD-NP. BD-P patients also presented bigger ventricles but similar hippocampal volumes compared to BD-NP/HC. Although the results are contrasting, some cognitive deficits seemed to be related to the psychotic dimension of bipolar affective disorder, such as impairment in verbal/logical memory, working memory, verbal and semantic fluency and executive functioning. Finally, polymorphisms of genes, such as NRG1, 5HTTLPR (s), COMT, DAOA and some chromosome regions (16p12 and 13q), were positively associated with BD-P.

CONCLUSION

Data about the identification of specific biomarkers for BD-P are promising, but most of them have not yet been replicated. They could lead the clinicians to an early diagnosis and proper treatment, thus ameliorating outcome of BD-P and reducing the biological changes associated with a long duration of illness. Further studies with bigger samples are needed to detect more specific biological markers of the psychotic dimension of bipolar affective disorder.

On the validity of within-nuclear-family genetic association analysis in samples of extended families

Splitting extended families into their component nuclear families to apply a genetic association method designed for nuclear families is a widespread practice in familial genetic studies. Dependence among genotypes and phenotypes of nuclear families from the same extended family arises because of genetic linkage of the tested marker with a risk variant or because of familial specificity of genetic effects due to gene-environment interaction. This raises concerns about the validity of inference conducted under the assumption of independence of the nuclear families. We indeed prove theoretically that, in a conditional logistic regression analysis applicable to disease cases and their genotyped parents, the naive model-based estimator of the variance of the coefficient estimates underestimates the true variance. However, simulations with realistic effect sizes of risk variants and variation of this effect from family to family reveal that the underestimation is negligible. The simulations also show the greater efficiency of the model-based variance estimator compared to a robust empirical estimator. Our recommendation is therefore, to use the model-based estimator of variance for inference on effects of genetic variants.

The genetics of early-onset bipolar disorder: A systematic review

BACKGROUND

Early-onset bipolar disorder has been associated with a significantly worse prognosis than late-onset BD and has been hypothesized to be a genetically homogenous subset of BD. A sizeable number of studies have investigated early-onset BD through linkage-analyses, candidate-gene association studies, genome-wide association studies (GWAS), and analyses of copy number variants (CNVs), but this literature has not yet been reviewed.

METHODS

A systematic review was conducted using the PubMed database on articles published online before January 15, 2015 and after 1990. Separate searches were made for linkage studies, candidate gene-association studies, GWAS, and studies on CNVs.

RESULTS

Seventy-three studies were included in our review. There is a lack of robust positive findings on the genetics of early-onset BD in any major molecular genetics method.

LIMITATIONS

Early-onset populations were quite small in some studies. Variance in study methods hindered efforts to interpret results or conduct meta-analysis.

CONCLUSIONS

The field is still at an early phase for research on early-onset BD. The largely null findings mirror the results of most genetics research on BD. Although most studies were underpowered, the null findings could mean that early-onset BD may not be as genetically homogenous as has been hypothesized or even that early-onset BD does not differ genetically from adult-onset BD. Nevertheless, clinically the probabilistic developmental risk trajectories associated with early-onset that may not be primarily genetically determined continued to warrant scrutiny. Future research should dramatically expand sample sizes, use atheoretical research methods like GWAS, and standardize methods.

Extension of the generalized disequilibrium test to polytomous phenotypes and two-locus models

We extend the usual logistic model between a dichotomous phenotype and an allele count in two ways: a polytomous phenotype with K > 2 levels, and modeling of allele counts at two unlinked marker loci. Inference is based on within-family information to guard against potential bias due to population genetic structure. Score tests of the model coefficients taking into account the correlation between relatives in entire pedigrees are derived as an extension of the Generalized Disequilibrium Test (GDT). Simulations confirm that the tests have the expected statistical properties, and that their power exceeds that of the GDT under a favorable scenario. The score tests are illustrated with candidate genetic markers, a major psychosis phenotype and a cognitive endophenotype in large kindreds from Eastern Quebec.