Jiewei Liu

Impact of a cis-associated gene expression SNP on chromosome 20q11.22 on bipolar disorder susceptibility, hippocampal structure and cognitive performance

BACKGROUND Bipolar disorder is a highly heritable polygenic disorder. Recent enrichment analyses suggest that there may be true risk variants for bipolar disorder in the expression quantitative trait loci (eQTL) in the brain. AIMS We sought to assess the impact of eQTL variants on bipolar disorder risk by combining data from both bipolar disorder genome-wide association studies (GWAS) and brain eQTL. METHOD To detect single nucleotide polymorphisms (SNPs) that influence expression levels of genes associated with bipolar disorder, we jointly analysed data from a bipolar disorder GWAS (7481 cases and 9250 controls) and a genome-wide brain (cortical) eQTL (193 healthy controls) using a Bayesian statistical method, with independent follow-up replications. The identified risk SNP was then further tested for association with hippocampal volume (n = 5775) and cognitive performance (n = 342) among healthy individuals. RESULTS Integrative analysis revealed a significant association between a brain eQTL rs6088662 on chromosome 20q11.22 and bipolar disorder (log Bayes factor = 5.48; bipolar disorder P = 5.85 × 10(-5)). Follow-up studies across multiple independent samples confirmed the association of the risk SNP (rs6088662) with gene expression and bipolar disorder susceptibility (P = 3.54 × 10(-8)). Further exploratory analysis revealed that rs6088662 is also associated with hippocampal volume and cognitive performance in healthy individuals. CONCLUSIONS Our findings suggest that 20q11.22 is likely a risk region for bipolar disorder; they also highlight the informative value of integrating functional annotation of genetic variants for gene expression in advancing our understanding of the biological basis underlying complex disorders, such as bipolar disorder.

Adaptive evolution of interleukin-3 (IL3), a gene associated with brain volume variation in general human populations.

Greatly expanded brain volume is one of the most characteristic traits that distinguish humans from other primates. Recent studies have revealed genes responsible for the dramatically enlarged human brain size (i.e., the microcephaly genes), and it has been well documented that many microcephaly genes have undergone accelerated evolution along the human lineage. In addition to being far larger than other primates, human brain volume is also highly variable in general populations. However, the genetic basis underlying human brain volume variation remains elusive and it is not known whether genes regulating human brain volume variation also have experienced positive selection. We have previously shown that genetic variants (near the IL3 gene) on 5q33 were significantly associated with brain volume in Chinese population. Here, we provide further evidence that support the significant association of genetic variants on 5q33 with brain volume. Bioinformatic analyses suggested that rs31480 is likely to be the causal variant among the studied SNPs. Molecular evolutionary analyses suggested that IL3 might have undergone positive selection in primates and humans. Neutrality tests further revealed signatures of positive selection of IL3 in Han Chinese and Europeans. Finally, extended haplotype homozygosity (EHH) and relative EHH analyses showed that the C allele of SNP rs31480 might have experienced recent positive selection in Han Chinese. Our results suggest that IL3 is an important genetic regulator for human brain volume variation and implied that IL3 might have experienced weak or modest positive selection in the evolutionary history of humans, which may be due to its contribution to human brain volume.

Allelic variation at 5-HTTLPR is associated with brain morphology in a Chinese population

Previous studies have reported significant associations of 5-HTTLPR with brain structures mainly in Europeans, but the situations in other ethnic groups remain largely unknown. Here we examined the association of 5-HTTLPR with regional gray matter volume in Han Chinese, and observed significant association in the postcentral gyrus and precuneus cortex.

Regional selection of the brain size regulating gene CASC5 provides new insight into human brain evolution

Human evolution is marked by a continued enlargement of the brain. Previous studies on human brain evolution focused on identifying sequence divergences of brain size regulating genes between humans and nonhuman primates. However, the evolutionary pattern of the brain size regulating genes during recent human evolution is largely unknown. We conducted a comprehensive analysis of the brain size regulating gene CASC5 and found that in recent human evolution, CASC5 has accumulated many modern human specific amino acid changes, including two fixed changes and six polymorphic changes. Among human populations, 4 of the 6 amino acid polymorphic sites have high frequencies of derived alleles in East Asians, but are rare in Europeans and Africans. We proved that this between-population allelic divergence was caused by regional Darwinian positive selection in East Asians. Further analysis of brain image data of Han Chinese showed significant associations of the amino acid polymorphic sites with gray matter volume. Hence, CASC5 may contribute to the morphological and structural changes of the human brain during recent evolution. The observed between-population divergence of CASC5 variants was driven by natural selection that tends to favor a larger gray matter volume in East Asians.

GWAS-identified schizophrenia risk SNPs at TSPAN18 are highly diverged between Europeans and East Asians.

Genome‐wide association studies (GWASs) have identified multiple schizophrenia (SCZ) risk variants for samples of European and East Asian descent, but most of the identified susceptibility variants are population‐specific to either Europeans or East Asians. This strong genetic heterogeneity suggests that differential population histories may play a role in SCZ susceptibility. Here, we explored this possibility by examining the allele frequency divergence of 136 previously reported genome‐wide SCZ risk SNPs between European and East Asian populations. Our results showed that two SNPs (rs11038167 and rs11038172) at TSPAN18, reported as genome‐wide significant SCZ risk variants in Han Chinese, were entirely monomorphic in Europeans, indicating a deep between‐population divergence at this gene locus. To explore the evolutionary history of TSPAN18 in East Asians, we conducted population genetic analyses including multiple neutrality tests, the haplotype‐based iHS and EHH tests, as well as haplotype bifurcation map and network constructions. We found that the protective allele of rs11038172 (G allele) had a long extended haplotype with much slower decay compared to the A allele. The star‐like shape of the G‐allele‐carrying haplotypes indicates a recent enrichment in East Asians. Together, the evidences suggest that the protective allele of rs11038172 has experienced recent Darwinian positive selection in East Asians. These findings provide new insights that may help explain the strong genetic heterogeneity in SCZ risk and previous inconsistent association results for SCZ among both Europeans and East Asians.

Failure of replicating the association between hippocampal volume and 3 single-nucleotide polymorphisms identified from the European genome-wide association study in Asian populations

Hippocampal volume is a key brain structure for learning ability and memory process, and hippocampal atrophy is a recognized biological marker of Alzheimers disease. However, the genetic bases of hippocampal volume are still unclear although it is a heritable trait. Genome-wide association studies (GWASs) on hippocampal volume have implicated several significantly associated genetic variants in Europeans. Here, to test the contributions of these GWASs identified genetic variants to hippocampal volume in different ethnic populations, we screened the GWAS-identified candidate single-nucleotide polymorphisms in 3 independent healthy Asian brain imaging samples (a total of 990 subjects). The results showed that none of these single-nucleotide polymorphisms were associated with hippocampal volume in either individual or combined Asian samples. The replication results suggested a complexity of genetic architecture for hippocampal volume and potential genetic heterogeneity between different ethnic populations.

An evaluation of association between a novel hippocampal biology related SNP (rs7294919) and schizophrenia.

Recent genetic analyses have implicated several candidate susceptibility variants for schizophrenia. The single nucleotide polymorphism (SNP) rs7294919 is likely a schizophrenia-susceptibility variant according to its significant association with hippocampal volume, hippocampus function, and cognitive performance as well as the nominal association with schizophrenia. However, all previous analyses were conducted only in Europeans, and whether rs7294919 is associated with schizophrenia in other populations are yet to be tested. Here, we conducted a case-control analysis of rs7294919 with schizophrenia in six independent Chinese (N = 3) and Japanese (N = 3) samples, including a total of 7,352 cases and 10,824 controls. The results of our association analysis were not able to confirm the association of rs7294919 with schizophrenia (p = 0.51 in total samples, odds ratio = 1.02 for allele[C]). The absence of rs7294919’s association in Chinese and Japanese suggest a potential genetic heterogeneity in the susceptibility of schizophrenia on this locus and also demonstrate the difficulties in replicating associations of schizophrenia across different ethnic populations.

Systems-level analysis of risk genes reveals the modular nature of schizophrenia

Schizophrenia (SCZ) is a complex mental disorder with high heritability. Genetic studies (especially recent genome-wide association studies) have identified many risk genes for schizophrenia. However, the physical interactions among the proteins encoded by schizophrenia risk genes remain elusive and it is not known whether the identified risk genes converge on common molecular networks or pathways. Here we systematically investigated the network characteristics of schizophrenia risk genes using the high-confidence protein-protein interactions (PPI) from the human interactome. We found that schizophrenia risk genes encode a densely interconnected PPI network (P = 4.15 × 10-31). Compared with the background genes, the schizophrenia risk genes in the interactome have significantly higher degree (P = 5.39 × 10-11), closeness centrality (P = 7.56 × 10-11), betweeness centrality (P = 1.29 × 10-11), clustering coefficient (P = 2.22 × 10-2), and shorter average shortest path length (P = 7.56 × 10-11). Based on the densely interconnected PPI network, we identified 48 hub genes and 4 modules formed by highly interconnected schizophrenia genes. We showed that the proteins encoded by schizophrenia hub genes have significantly more direct physical interactions. Gene ontology (GO) analysis revealed that cell adhesion, cell cycle, immune system response, and GABR-receptor complex categories were enriched in the modules formed by highly interconnected schizophrenia risk genes. Our study reveals that schizophrenia risk genes encode a densely interconnected molecular network and demonstrates the modular nature of schizophrenia.

Comprehensive integrative analyses identify GLT8D1 and CSNK2B as schizophrenia risk genes

Recent genome-wide association studies (GWAS) have identified multiple risk loci that show strong associations with schizophrenia. However, pinpointing the potential causal genes at the reported loci remains a major challenge. Here we identify candidate causal genes for schizophrenia using an integrative genomic approach. Sherlock integrative analysis shows that ALMS1, GLT8D1, and CSNK2B are schizophrenia risk genes, which are validated using independent brain expression quantitative trait loci (eQTL) data and integrative analysis method (SMR). Consistently, gene expression analysis in schizophrenia cases and controls further supports the potential role of these three genes in the pathogenesis of schizophrenia. Finally, we show that GLT8D1 and CSNK2B knockdown promote the proliferation and inhibit the differentiation abilities of neural stem cells, and alter morphology and synaptic transmission of neurons. These convergent lines of evidence suggest that the ALMS1, CSNK2B, and GLT8D1 genes may be involved in pathophysiology of schizophrenia.More than 100 risk loci for schizophrenia have been identified by genome-wide association studies. Here, the authors apply an integrative genomic approach to prioritize risk genes and validate GLT8D1 and CSNK2B as candidate causal genes by in vitro studies in neural stem cells.

Integrated analysis supports ATXN1 as a schizophrenia risk gene

Protein-protein interaction (PPI) is informative in identifying hidden disease risk genes that tend to interact with known risk genes usually working together in the same disease module. With the use of an integrated approach combining PPI information with pathway and expression analysis as well as genome-wide association study (GWAS), we intended to find new risk genes for schizophrenia (SCZ). We showed that ATXN1 was the only direct PPI partner of the know SCZ risk gene ZNF804A, and it also had direct PPIs with other 18 known SCZ risk genes. ATXN1 serves as one of the hub genes in the PPI network containing many known SCZ risk genes, and this network is significantly enriched for the MAPK signaling pathway. Further gene expression analysis indicated that ATXN1 is highly expressed in prefrontal cortex, and SCZ patients had significantly decreased expression compared with healthy controls. Finally, the published GWAS data supports an association of ATXN1 with SCZ as well as other psychiatric disorders though not reaching genome-wide significance. These convergent evidences support ATXN1 as a promising risk gene for SCZ, and the integrated approach serves as a useful tool for dissecting the genetic basis of psychiatric disorders.