Association between SNPs and gene expression in multiple regions of the human brain

Genome-wide association study identified candidate genes for egg production traits in the Longyan Shan-ma duck

Egg production is an important economic trait in layer ducks and understanding the genetics basis is important for their breeding. In this study, a genome-wide association study (GWAS) for egg production traits in 303 female Longyan Shan-ma ducks was performed based on a genotyping-by-sequencing strategy. Sixty-two single nucleotide polymorphisms (SNPs) associated with egg weight traits were identified (P < 9.48 × 10-5), including 8 SNPs at 5% linkage disequilibrium (LD)-based Bonferroni-corrected genome-wide significance level (P < 4.74 × 10-6). One hundred and nineteen SNPs were associated with egg number traits (P < 9.48 × 10-5), including 13 SNPs with 5% LD-based Bonferroni-corrected genome-wide significance (P < 4.74 × 10-6). These SNPs annotated 146 target genes which contained known candidate genes for egg production traits, such as prolactin and prolactin releasing hormone receptor. This study identified that these associated genes were significantly enriched in egg production-related pathways (P < 0.05), such as the oxytocin signaling, MAPK signaling, and calcium signaling pathways. It was notable that 18 genes were differentially expressed in ovarian tissues between higher and lower egg production in Shan-ma ducks. The identified potential candidate genes and pathways provide insight into the genetic basis underlying the egg production trait of layer ducks.

Screening and analysis of single nucleotide polymorphism in the 3'-UTR microRNA target regions and its implications for lung tumorigenesis

Aim: The study aims to identify high-impact single nucleotide polymorphisms (SNPs) in miRNA target sites of genes associated with lung cancer. Materials & methods: Lung cancer genes were obtained from Uniprot KB. miRNA target site SNPs were mined from MirSNP, miRdSNP and TargetScan. SNPs were shortlisted based on binding impact, minor allele frequency and conservation. Gene expression was analyzed in genes with high-impact SNPs in healthy versus lung cancer tissue. Additionally, enrichment, pathway and network analyzes were performed. Results: 19 high-impact SNPs were identified in miRNA target sites of lung cancer-associated genes. These SNPs affect miRNA binding and gene expression. The genes are involved in key cancer related pathways. Conclusion: The identified high-impact miRNA target site SNPs and associated genes provide a starting point for case-control studies in lung cancer patients in different populations.

Small molecule probes as versatile energy acceptors: A breakthrough in photoelectrochemical sensing for sulfur dioxide recording in rat brain

Microelectrode-based photoelectrochemical (PEC) sensing is a newly developed and promising analytical technique for in vivo analysis. However, the inadequate specificity in complex environment of living bodies restricted its further in vivo application. Herein, we utilized a small molecule probe as the energy acceptor to quench the photocurrent of CdTe quantum dots through energy transfer. The efficiency of energy transfer was modulated by the concentration of target SO2, resulting in changes in photocurrent. The chemical recognition reaction between small molecule probes and SO2 enhanced the specificity of PEC sensing, thus guaranteeing its in vivo applications. Furthermore, with the use of light addressing strategy, simultaneous detection in the multiple brain regions was implemented. The energy transfer based light addressable PEC microsensor achieved monitoring fluctuations of SO2 levels in multiple brain regions of rats with epilepsy.

The causal relationship of colorectal cancer on schizophrenia: A Mendelian randomization study

Comorbidities associated with psychiatric disorders often occur in patients with cancer. A causal effect of schizophrenia on cancer was observed using Mendelian randomization (MR) analysis. However, the causal effect of colorectal cancer on schizophrenia has not been studied using MR analysis. Therefore, we performed MR analysis to investigate the causal effects of colorectal cancer on schizophrenia. We performed "two-sample summary-data Mendelian randomization" using publicly available genome-wide association studies data to investigate the causal relationship between colorectal cancer (as exposure) and schizophrenia (as outcome). The inverse variance weighted method was used to calculate causal estimates. In 2 TSMR analyses, we reported that the odds ratios for schizophrenia per log odds increase in colorectal cancer risk were 6.48 (95% confidential interval [CI] of OR 1.75-24.03; P = .005) and 9.62 × 106 (95% CI of OR 1.13-8.22 × 1013; P = .048). Pleiotropic tests and sensitivity analysis demonstrated minimal horizontal pleiotropy and robustness of the causal relationship. We provide evidence for a causal relationship between the incidence of colorectal cancer and the development of schizophrenia through TSMR analysis.

A Mouse Systems Genetics Approach Reveals Common and Uncommon Genetic Modifiers of Hepatic Lysosomal Enzyme Activities and Glycosphingolipids

Identification of genetic modulators of lysosomal enzyme activities and glycosphingolipids (GSLs) may facilitate the development of therapeutics for diseases in which they participate, including Lysosomal Storage Disorders (LSDs). To this end, we used a systems genetics approach: we measured 11 hepatic lysosomal enzymes and many of their natural substrates (GSLs), followed by modifier gene mapping by GWAS and transcriptomics associations in a panel of inbred strains. Unexpectedly, most GSLs showed no association between their levels and the enzyme activity that catabolizes them. Genomic mapping identified 30 shared predicted modifier genes between the enzymes and GSLs, which are clustered in three pathways and are associated with other diseases. Surprisingly, they are regulated by ten common transcription factors, and their majority by miRNA-340p. In conclusion, we have identified novel regulators of GSL metabolism, which may serve as therapeutic targets for LSDs and may suggest the involvement of GSL metabolism in other pathologies.

Genetic control of the human brain proteome

We generated an online brain pQTL resource for 7,376 proteins through the analysis of genetic and proteomic data derived from post-mortem samples of the dorsolateral prefrontal cortex of 330 older adults. The identified pQTLs tend to be non-synonymous variation, are over-represented among variants associated with brain diseases, and replicate well (77%) in an independent brain dataset. Comparison to a large study of brain eQTLs revealed that about 75% of pQTLs are also eQTLs. In contrast, about 40% of eQTLs were identified as pQTLs. These results are consistent with lower pQTL mapping power and greater evolutionary constraint on protein abundance. The latter is additionally supported by observations of pQTLs with large effects' tending to be rare, deleterious, and associated with proteins that have evidence for fewer protein-protein interactions. Mediation analyses using matched transcriptomic and proteomic data provided additional evidence that pQTL effects are often, but not always, mediated by mRNA. Specifically, we identified roughly 1.6 times more mRNA-mediated pQTLs than mRNA-independent pQTLs (550 versus 341). Our pQTL resource provides insight into the functional consequences of genetic variation in the human brain and a basis for novel investigations of genetics and disease.

Calcium-binding proteins are altered in the cerebellum in schizophrenia

Alterations in the cortico-cerebellar-thalamic-cortical circuit might underlie the diversity of symptoms in schizophrenia. However, molecular changes in cerebellar neuronal circuits, part of this network, have not yet been fully determined. Using LC-MS/MS, we screened altered candidates in pooled grey matter of cerebellum from schizophrenia subjects who committed suicide (n = 4) and healthy individuals (n = 4). Further validation by immunoblotting of three selected candidates was performed in two cohorts comprising schizophrenia (n = 20), non-schizophrenia suicide (n = 6) and healthy controls (n = 21). We found 99 significantly altered proteins, 31 of them previously reported in other brain areas by proteomic studies. Transport function was the most enriched category, while cell communication was the most prevalent function. For validation, we selected the vacuolar proton pump subunit 1 (VPP1), from transport, and two EF-hand calcium-binding proteins, calmodulin and parvalbumin, from cell communication. All candidates showed significant changes in schizophrenia (n = 7) compared to controls (n = 7). VPP1 was altered in the non-schizophrenia suicide group and increased levels of parvalbumin were linked to antipsychotics. Further validation in an independent cohort of non-suicidal chronic schizophrenia subjects (n = 13) and non-psychiatric controls (n = 14) showed that parvalbumin was increased, while calmodulin was decreased in schizophrenia. Our findings provide evidence of calcium-binding protein dysregulation in the cerebellum in schizophrenia, suggesting an impact on normal calcium-dependent synaptic functioning of cerebellar circuits. Our study also links VPP1 to suicide behaviours, suggesting a possible impairment in vesicle neurotransmitter refilling and release in these phenotypes.

Ethnicity-Dependent Effects of Schizophrenia Risk Variants of the OLIG2 Gene on OLIG2 Transcription and White Matter Integrity

Previous studies have indicated associations between several OLIG2 gene single-nucleotide polymorphisms (SNPs) and susceptibility to schizophrenia among Caucasians. Consistent with these findings, postmortem brain and diffusion tensor imaging studies have indicated that the schizophrenia-risk-associated allele (A) in the OLIG2 SNP rs1059004 predicts lower OLIG2 gene expression in the dorsolateral prefrontal cortex (DLPFC) of schizophrenia patients and reduced white matter (WM) integrity of the corona radiata in normal brains among Caucasians. In an effort to replicate the association between this variant and WM integrity among healthy Japanese, we found that the number of A alleles was positively correlated with WM integrity in some fiber tracts, including the right posterior limb of the internal capsule, and with mean blood flow in a widespread area, including the inferior frontal operculum, orbital area, and triangular gyrus. Because the A allele affected WM integrity in opposite directions in Japanese and Caucasians, we investigated a possible association between the OLIG2 gene SNPs and the expression level of OLIG2 transcripts in postmortem DLPFCs. We evaluated rs1059004 and additional SNPs in the 5' upstream and 3' downstream regions of rs1059004 to cover the broader region of the OLIG2 gene. The 2 SNPs (rs1059004 and rs9653711) had opposite effects on OLIG2 gene expression in the DLPFC in Japanese and Caucasians. These findings suggest ethnicity-dependent opposite effects of OLIG2 gene SNPs on WM integrity and OLIG2 gene expression in the brain, which may partially explain the failures in replicating associations between genetic variants and psychiatric phenotypes among ethnicities.

Predicting novel genomic regions linked to genetic disorders using GWAS and chromosome conformation data - a case study of schizophrenia

Genome-wide association studies identified numerous loci harbouring single nucleotide polymorphisms (SNPs) associated with various human diseases, although the causal role of many of them remains unknown. In this paper, we postulate that co-location and shared biological function of novel genes with genes known to associate with a specific phenotype make them potential candidates linked to the same phenotype (“guilt-by-proxy”). We propose a novel network-based approach for predicting candidate genes/genomic regions utilising the knowledge of the 3D architecture of the human genome and GWAS data. As a case study we used a well-studied polygenic disorder ‒ schizophrenia ‒ for which we compiled a comprehensive dataset of SNPs. Our approach revealed 634 novel regions covering ~398 Mb of the human genome and harbouring ~9000 genes. Using various network measures and enrichment analysis, we identified subsets of genes and investigated the plausibility of these genes/regions having an association with schizophrenia using literature search and bioinformatics resources. We identified several genes/regions with previously reported associations with schizophrenia, thus providing proof-of-concept, as well as novel candidates with no prior known associations. This approach has the potential to identify novel genes/genomic regions linked to other polygenic disorders and provide means of aggregating genes/SNPs for further investigation.

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.

PlexinA2 Forward Signaling through Rap1 GTPases Regulates Dentate Gyrus Development and Schizophrenia-like Behaviors

Dentate gyrus (DG) development requires specification of granule cell (GC) progenitors in the hippocampal neuroepithelium, as well as their proliferation and migration into the primordial DG. We identify the Plexin family members Plxna2 and Plxna4 as important regulators of DG development. Distribution of immature GCs is regulated by Sema5A signaling through PlxnA2 and requires a functional PlxnA2 GTPase-activating protein (GAP) domain and Rap1 small GTPases. In adult Plxna2^−/− but not Plxna2-GAP-deficient mice, the dentate GC layer is severely malformed, neurogenesis is compromised, and mossy fibers form aberrant synaptic boutons within CA3. Behavioral studies with Plxna2^−/− mice revealed deficits in associative learning, sociability, and sensorimotor gating—traits commonly observed in neuropsychiatric disorder. Remarkably, while morphological defects are minimal in Plxna2-GAP-deficient brains, defects in fear memory and sensorimotor gating persist. Since allelic variants of human PLXNA2 and RAP1 associate with schizophrenia, our studies identify a biochemical pathway important for brain development and mental health.

Dyslexia risk variant rs600753 is linked with dyslexia-specific differential allelic expression of DYX1C1

An increasing number of genetic variants involved in dyslexia development were discovered during the last years, yet little is known about the molecular functional mechanisms of these SNPs. In this study we investigated whether dyslexia candidate SNPs have a direct, disease-specific effect on local expression levels of the assumed target gene by using a differential allelic expression assay. In total, 12 SNPs previously associated with dyslexia and related phenotypes were suitable for analysis. Transcripts corresponding to four SNPs were sufficiently expressed in 28 cell lines originating from controls and a family affected by dyslexia. We observed a significant effect of rs600753 on expression levels of DYX1C1 in forward and reverse sequencing approaches. The expression level of the rs600753 risk allele was increased in the respective seven cell lines from members of the dyslexia family which might be due to a disturbed transcription factor binding sites. When considering our results in the context of neuroanatomical dyslexia-specific findings, we speculate that this mechanism may be part of the pathomechanisms underlying the dyslexia-specific brain phenotype. Our results suggest that allele-specific DYX1C1 expression levels depend on genetic variants of rs600753 and contribute to dyslexia. However, these results are preliminary and need replication.

Cis-acting single nucleotide polymorphisms alter MicroRNA-mediated regulation of human brain-expressed transcripts

Substantial variability exists in the presentation of complex neurological disorders, and the study of single nucleotide polymorphisms (SNPs) has shed light on disease mechanisms and pathophysiological variability in some cases. However, the vast majority of disease-linked SNPs have unidentified pathophysiological relevance. Here, we tested the hypothesis that SNPs within the miRNA recognition element (MRE; the region of the target transcript to which the miRNA binds) can impart changes in the expression of those genes, either by enhancing or reducing transcript and protein levels. To test this, we cross-referenced 7,153 miRNA-MRE brain interactions with the SNP database (dbSNP) to identify candidates, and functionally assessed 24 SNPs located in the 3’UTR or the coding sequence (CDS) of targets. For over half of the candidates tested, SNPs either enhanced (4 genes) or disrupted (10 genes) miRNA binding and target regulation. Additionally, SNPs causing a shift from a common to rare codon within the CDS facilitated miRNA binding downstream of the SNP, dramatically repressing target gene expression. The biological activity of the SNPs on miRNA regulation was also confirmed in induced pluripotent stem cell (iPSC) lines. These studies strongly support the notion that SNPs in the 3’UTR or the coding sequence of disease-relevant genes may be important in disease pathogenesis and should be reconsidered as candidate modifiers.

PlexinA2 Forward Signaling through Rap1 GTPases Regulates Dentate Gyrus Development and Schizophrenia-like Behaviors

Dentate gyrus (DG) development requires specification of granule cell (GC) progenitors in the hippocampal neuroepithelium, as well as their proliferation and migration into the primordial DG. We identify the Plexin family members Plxna2 and Plxna4 as important regulators of DG development. Distribution of immature GCs is regulated by Sema5A signaling through PlxnA2 and requires a functional PlxnA2 GTPase-activating protein (GAP) domain and Rap1 small GTPases. In adult Plxna2^-/- but not Plxna2-GAP-deficient mice, the dentate GC layer is severely malformed, neurogenesis is compromised, and mossy fibers form aberrant synaptic boutons within CA3. Behavioral studies with Plxna2^-/- mice revealed deficits in associative learning, sociability, and sensorimotor gating-traits commonly observed in neuropsychiatric disorder. Remarkably, while morphological defects are minimal in Plxna2-GAP-deficient brains, defects in fear memory and sensorimotor gating persist. Since allelic variants of human PLXNA2 and RAP1 associate with schizophrenia, our studies identify a biochemical pathway important for brain development and mental health.

ATP2C2 and DYX1C1 are putative modulators of dyslexia-related MMR

BACKGROUND

Dyslexia is a specific learning disorder affecting reading and spelling abilities. Its prevalence is ~5% in German-speaking individuals. Although the etiology of dyslexia largely remains to be determined, comprehensive evidence supports deficient phonological processing as a major contributing factor. An important prerequisite for phonological processing is auditory discrimination and, thus, essential for acquiring reading and spelling skills. The event-related potential Mismatch Response (MMR) is an indicator for auditory discrimination capabilities with dyslexics showing an altered late component of MMR in response to auditory input.

METHODS

In this study, we comprehensively analyzed associations of dyslexia-specific late MMRs with genetic variants previously reported to be associated with dyslexia-related phenotypes in multiple studies comprising 25 independent single-nucleotide polymorphisms (SNPs) within 10 genes.

RESULTS

First, we demonstrated validity of these SNPs for dyslexia in our sample by showing that additional inclusion of a polygenic risk score improved prediction of impaired writing compared with a model that used MMR alone. Secondly, a multifactorial regression analysis was conducted to uncover the subset of the 25 SNPs that is associated with the dyslexia-specific late component of MMR. In total, four independent SNPs within DYX1C1 and ATP2C2 were found to be associated with MMR stronger than expected from multiple testing. To explore potential pathomechanisms, we annotated these variants with functional data including tissue-specific expression analysis and eQTLs.

CONCLUSION

Our findings corroborate the late component of MMR as a potential endophenotype for dyslexia and support tripartite relationships between dyslexia-related SNPs, the late component of MMR and dyslexia.

The epigenomics of schizophrenia, in the mouse

Large-scale consortia including the Psychiatric Genomics Consortium, the Common Minds Consortium, BrainSeq and PsychENCODE, and many other studies taken together provide increasingly detailed insights into the genetic and epigenetic risk architectures of schizophrenia (SCZ) and offer vast amounts of molecular information, but with largely unexplored therapeutic potential. Here we discuss how epigenomic studies in human brain could guide animal work to test the impact of disease-associated alterations in chromatin structure and function on cognition and behavior. For example, transcription factors such as MYOCYTE-SPECIFIC ENHANCER FACTOR 2C (MEF2C), or multiple regulators of the open chromatin mark, methyl-histone H3-lysine 4, are associated with the genetic risk architectures of common psychiatric disease and alterations in chromatin structure and function in diseased brain tissue. Importantly, these molecules also affect cognition and behavior in genetically engineered mice, including virus-mediated expression changes in prefrontal cortex (PFC) and other key nodes in the circuitry underlying psychosis. Therefore, preclinical and small laboratory animal work could target genomic sequences affected by chromatin alterations in SCZ. To this end, in vivo editing of enhancer and other regulatory non-coding DNA by RNA-guided nucleases including CRISPR-Cas, and designer transcription factors, could be expected to deliver pipelines for novel therapeutic approaches aimed at improving cognitive dysfunction and other core symptoms of SCZ.

Genetic regulation of gene expression in the epileptic human hippocampus

Epilepsy is a serious and common neurological disorder. Expression quantitative loci (eQTL) analysis is a vital aid for the identification and interpretation of disease-risk loci. Many eQTLs operate in a tissue- and condition-specific manner. We have performed the first genome-wide cis-eQTL analysis of human hippocampal tissue to include not only normal (n = 22) but also epileptic (n = 22) samples. We demonstrate that disease-associated variants from an epilepsy GWAS meta-analysis and a febrile seizures (FS) GWAS are significantly more enriched with epilepsy-eQTLs than with normal hippocampal eQTLs from two larger independent published studies. In contrast, GWAS meta-analyses of two other brain diseases associated with hippocampal pathology (Alzheimer's disease and schizophrenia) are more enriched with normal hippocampal eQTLs than with epilepsy-eQTLs. These observations suggest that an eQTL analysis that includes disease-affected brain tissue is advantageous for detecting additional risk SNPs for the afflicting and closely related disorders, but not for distinct diseases affecting the same brain regions. We also show that epilepsy eQTLs are enriched within epilepsy-causing genes: an epilepsy cis-gene is significantly more likely to be a causal gene for a Mendelian epilepsy syndrome than to be a causal gene for another Mendelian disorder. Epilepsy cis-genes, compared to normal hippocampal cis-genes, are more enriched within epilepsy-causing genes. Hence, we utilize the epilepsy eQTL data for the functional interpretation of epilepsy disease-risk variants and, thereby, highlight novel potential causal genes for sporadic epilepsy. In conclusion, an epilepsy-eQTL analysis is superior to normal hippocampal tissue eQTL analyses for identifying the variants and genes underlying epilepsy.

Decreased SMG7 expression associates with lupus-risk variants and elevated antinuclear antibody production

OBJECTIVES

Following up the systemic lupus erythematosus (SLE) genome-wide association studies (GWAS) identification of NMNAT2 at rs2022013, we fine-mapped its 150 kb flanking regions containing NMNAT2 and SMG7 in a 15 292 case-control multi-ancestry population and tested functions of identified variants.

METHODS

We performed genotyping using custom array, imputation by IMPUTE 2.1.2 and allele specific functions using quantitative real-time PCR and luciferase reporter transfections. SLE peripheral blood mononuclear cells (PBMCs) were cultured with small interfering RNAs to measure antinuclear antibody (ANA) and cyto/chemokine levels in supernatants using ELISA.

RESULTS

We confirmed association at NMNAT2 in European American (EA) and Amerindian/Hispanic ancestries, and identified independent signal at SMG7 tagged by rs2702178 in EA only (p=2.4×10-8, OR=1.23 (95% CI 1.14 to 1.32)). In complete linkage disequilibrium with rs2702178, rs2275675 in the promoter region robustly associated with SMG7 mRNA levels in multiple expression quantitative trait locus (eQTL) datasets. Its risk allele was dose-dependently associated with decreased SMG7 mRNA levels in PBMCs of 86 patients with SLE and 119 controls (p=1.1×10-3 and 6.8×10-8, respectively) and conferred reduced transcription activity in transfected HEK-293 (human embryonic kidney cell line) and Raji cells (p=0.0035 and 0.0037, respectively). As a critical component in the nonsense-mediated mRNA decay pathway, SMG7 could regulate autoantigens including ribonucleoprotein (RNP) and Smith (Sm). We showed SMG7 mRNA levels in PBMCs correlated inversely with ANA titres of patients with SLE (r=-0.31, p=0.01), and SMG7 knockdown increased levels of ANA IgG and chemokine (C-C motif) ligand 19 in SLE PBMCs (p=2.0×10-5 and 2.0×10-4, respectively).

CONCLUSION

We confirmed NMNAT2 and identified independent SMG7 association with SLE. The inverse relationship between levels of the risk allele-associated SMG7 mRNAs and ANA suggested the novel contribution of mRNA surveillance pathway to SLE pathogenesis.

Four Susceptibility Loci for Gallstone Disease Identified in a Meta-analysis of Genome-Wide Association Studies

BACKGROUND & AIMS

A genome-wide association study (GWAS) of 280 cases identified the hepatic cholesterol transporter ABCG8 as a locus associated with risk for gallstone disease, but findings have not been reported from any other GWAS of this phenotype. We performed a large-scale, meta-analysis of GWASs of individuals of European ancestry with available prior genotype data, to identify additional genetic risk factors for gallstone disease.

METHODS

We obtained per-allele odds ratio (OR) and standard error estimates using age- and sex-adjusted logistic regression models within each of the 10 discovery studies (8720 cases and 55,152 controls). We performed an inverse variance weighted, fixed-effects meta-analysis of study-specific estimates to identify single-nucleotide polymorphisms that were associated independently with gallstone disease. Associations were replicated in 6489 cases and 62,797 controls.

RESULTS

We observed independent associations for 2 single-nucleotide polymorphisms at the ABCG8 locus: rs11887534 (OR, 1.69; 95% confidence interval [CI], 1.54-1.86; P = 2.44 × 10(-60)) and rs4245791 (OR, 1.27; P = 1.90 × 10(-34)). We also identified and/or replicated associations for rs9843304 in TM4SF4 (OR, 1.12; 95% CI, 1.08-1.16; P = 6.09 × 10(-11)), rs2547231 in SULT2A1 (encodes a sulfoconjugation enzyme that acts on hydroxysteroids and cholesterol-derived sterol bile acids) (OR, 1.17; 95% CI, 1.12-1.21; P = 2.24 × 10(-10)), rs1260326 in glucokinase regulatory protein (OR, 1.12; 95% CI, 1.07-1.17; P = 2.55 × 10(-10)), and rs6471717 near CYP7A1 (encodes an enzyme that catalyzes conversion of cholesterol to primary bile acids) (OR, 1.11; 95% CI, 1.08-1.15; P = 8.84 × 10(-9)). Among individuals of African American and Hispanic American ancestry, rs11887534 and rs4245791 were associated positively with gallstone disease risk, whereas the association for the rs1260326 variant was inverse.

CONCLUSIONS

In this large-scale GWAS of gallstone disease, we identified 4 loci in genes that have putative functions in cholesterol metabolism and transport, and sulfonylation of bile acids or hydroxysteroids.

Association study of MiRSNPs with schizophrenia, tardive dyskinesia and cognition

MicroRNAs (miRNAs) bind to 3'UTRs of genes and negatively regulate their expression. With ~50% of miRNAs expressing in the brain, they play an important role in neuronal development, plasticity, cognition and neurological disorders. Conserved miRNA targets are present in >60% genes in humans and are under evolutionary pressure to maintain pairing with miRNA. However, such binding may be affected by genetic variant(s) in the target sites (MiRSNPs), thereby altering gene expression. Differential expression of a large number of genes in postmortem brains of schizophrenia (SZ) patients compared to controls has been documented. Thus studying the role of MiRSNPs which are underinvestigated in SZ becomes attractive. We systematically selected 35 MiRSNPs with predicted functional relevance in 3'UTRs of genes shown previously to be associated with SZ, genotyped and tested their association with disease, using independent discovery and replication samples (total n=1017 cases; n=1073 controls). We also explored genetic associations with two sets of quantitative traits, namely tardive dyskinesia (TD) and cognitive functions disrupted in SZ in subsets of the study cohort. In the primary analysis, a significant association of MiRSNP rs7430 at PPP3CC was observed with SZ in the discovery and the replication samples [discovery: P=0.01; OR (95% CI) 1.24 (1.04-1.48); replication: P=0.03; OR (95% CI) 1.20 (1.02-1.43)]. In the exploratory analyses, five SNPs were nominally associated with TD (P values 0.04-0.004). Separately, 12 SNPs were associated with one or more of the eight cognitive domains (P values 0.05-0.003). These associations, particularly the SNP at PPP3CC merit further investigations.

Association, characterisation and meta-analysis of SNPs linked to general reading ability in a German dyslexia case-control cohort

Dyslexia is a severe disorder in the acquisition of reading and writing. Several studies investigated the role of genetics for reading, writing and spelling ability in the general population. However, many of the identified SNPs were not analysed in case-control cohorts. Here, we investigated SNPs previously linked to reading or spelling ability in the general population in a German case-control cohort. Furthermore, we characterised these SNPs for functional relevance with in silico methods and meta-analysed them with previous studies. A total of 16 SNPs within five genes were included. The total number of risk alleles was higher in cases than in controls. Three SNPs were nominally associated with dyslexia: rs7765678 within DCDC2, and rs2038137 and rs6935076 within KIAA0319. The relevance of rs2038137 and rs6935076 was further supported by the meta-analysis. Functional profiling included analysis of tissue-specific expression, annotations for regulatory elements and effects on gene expression levels (eQTLs). Thereby, we found molecular mechanistical implications for 13 of all 16 included SNPs. SNPs associated in our cohort showed stronger gene-specific eQTL effects than non-associated SNPs. In summary, our results validate SNPs previously linked to reading and spelling in the general population in dyslexics and provide insights into their putative molecular pathomechanisms.

Differential activation of immune/inflammatory response-related co-expression modules in the hippocampus across the major psychiatric disorders

The Stanley Neuropathology Consortium Integrative Database (SNCID, http://sncid.stanleyresearch.org) is a data-mining tool that includes 379 neuropathology data sets from hippocampus, as well as RNA-Seq data measured in 15 well-matched cases in each of four groups: schizophrenia, bipolar disorder (BPD), major depression (MD) and unaffected controls. We analyzed the neuropathology data from the hippocampus to identify those abnormalities that are shared between psychiatric disorders and those that are specific to each disorder. Of the 379 data sets, 20 of them showed a significant abnormality in at least one disorder as compared with unaffected controls. GABAergic markers and synaptic proteins were mainly abnormal in schizophrenia and the two mood disorders, respectively. Two immune/inflammation-related co-expression modules built from RNA-seq data from both schizophrenia and controls combined were associated with disease status, as well as negatively correlated with the GABAergic markers. The correlation between immune-related modules and schizophrenia was replicated using microarray data from an independent tissue collection. Immune/inflammation-related co-expression modules were also built from RNA-seq data from BPD cases or from MD cases but were not preserved when using data from control cases. Moreover, there was no overlap in the genes that comprise the immune/inflammation response-related modules across the different disorders. Thus, there appears to be differential activation of the immune/inflammatory response, as determined by co-expression of genes, which is associated with the major psychiatric disorders and which is also associated with the abnormal neuropathology in the disorders.

Glia-related genes and their contribution to schizophrenia

Schizophrenia, a debilitating disease with 1% prevalence in the general population, is characterized by major neuropsychiatric symptoms, including delusions, hallucinations, and deficits in emotional and social behavior. Previous studies have directed their investigations on the mechanism of schizophrenia towards neuronal dysfunction and have defined schizophrenia as a 'neuron-centric' disorder. However, along with the development of genetics and systematic biology approaches in recent years, the crucial role of glial cells in the brain has also been shown to contribute to the etiopathology of schizophrenia. Here, we summarize comprehensive data that support the involvement of glial cells (including oligodendrocytes, astrocytes, and microglial cells) in schizophrenia and list several acknowledged glia-related genes or molecules associated with schizophrenia. Instead of purely an abnormality of neurons in schizophrenia, an additional 'glial perspective' provides us a novel and promising insight into the causal mechanisms and treatment for this disease.

Understanding complex transcriptome dynamics in schizophrenia and other neurological diseases using RNA sequencing

How the human brain develops and adapts with its trillions of functionally integrated synapses remains one of the greatest mysteries of life. With tremendous advances in neuroscience, genetics, and molecular biology, we are beginning to appreciate the scope of this complexity and define some of the parameters of the systems that make it possible. These same tools are also leading to advances in our understanding of the pathophysiology of neurocognitive and neuropsychiatric disorders. Like the substrate for these problems, the etiology is usually complex-involving an array of genetic and environmental influences. To resolve these influences and derive better interventions, we need to reveal every aspect of this complexity and model their interactions and define the systems and their regulatory structure. This is particularly important at the tissue-specific molecular interface between the underlying genetic and environmental influence defined by the transcriptome. Recent advances in transcriptome analysis facilitated by RNA sequencing (RNA-Seq) can provide unprecedented insight into the functional genomics of neurological disorders. In this review, we outline the advantages of this approach and highlight some early application of this technology in the investigation of the neuropathology of schizophrenia. Recent progress of RNA-Seq studies in schizophrenia has shown that there is extraordinary transcriptome dynamics with significant levels of alternative splicing. These studies only scratch the surface of this complexity and therefore future studies with greater depth and samples size will be vital to fully explore transcriptional diversity and its underlying influences in schizophrenia and provide the basis for new biomarkers and improved treatments.

Schizophrenia as a disorder of molecular pathways

Over the last decade, transcriptome studies of postmortem tissue from subjects with schizophrenia revealed that synaptic, mitochondrial, immune system, gamma-aminobutyric acidergic, and oligodendrocytic changes are all integral parts of the disease process. The combined genetic and transcriptomic studies argue that the molecular underpinnings of the disease are even more varied than the symptomatic diversity of schizophrenia. Ultimately, to decipher the pathophysiology of human disorders in general, we will need to understand the function of hundreds of genes and regulatory elements in our genome and the consequences of their overexpression and reduced expression in a developmental context. Furthermore, integration of knowledge from various data sources remains a monumental challenge that has to be systematically addressed in the upcoming decades. In the end, our success in interpreting the molecular changes in schizophrenia will depend on our ability to understand the biology using innovative ideas and cannot depend on the hope of developing novel, more powerful technologies.

Genetic variability in the regulation of gene expression in ten regions of the human brain

Germ-line genetic control of gene expression occurs via expression quantitative trait loci (eQTLs). We present a large, exon-specific eQTL data set covering ten human brain regions. We found that cis-eQTL signals (within 1 Mb of their target gene) were numerous, and many acted heterogeneously among regions and exons. Co-regulation analysis of shared eQTL signals produced well-defined modules of region-specific co-regulated genes, in contrast to standard coexpression analysis of the same samples. We report cis-eQTL signals for 23.1% of catalogued genome-wide association study hits for adult-onset neurological disorders. The data set is publicly available via public data repositories and via http://www.braineac.org/. Our study increases our understanding of the regulation of gene expression in the human brain and will be of value to others pursuing functional follow-up of disease-associated variants.

Transcription factor 4 (TCF4) and schizophrenia: integrating the animal and the human perspective

Schizophrenia is a genetically complex disease considered to have a neurodevelopmental pathogenesis and defined by a broad spectrum of positive and negative symptoms as well as cognitive deficits. Recently, large genome-wide association studies have identified common alleles slightly increasing the risk for schizophrenia. Among the few schizophrenia-risk genes that have been consistently replicated is the basic Helix-Loop-Helix (bHLH) transcription factor 4 (TCF4). Haploinsufficiency of the TCF4 (formatting follows IUPAC nomenclature: TCF4 protein/protein function, Tcf4 rodent gene cDNA mRNA, TCF4 human gene cDNA mRNA) gene causes the Pitt-Hopkins syndrome-a neurodevelopmental disease characterized by severe mental retardation. Accordingly, Tcf4 null-mutant mice display developmental brain defects. TCF4-associated risk alleles are located in putative coding and non-coding regions of the gene. Hence, subtle changes at the level of gene expression might be relevant for the etiopathology of schizophrenia. Behavioural phenotypes obtained with a mouse model of slightly increased gene dosage and electrophysiological investigations with human risk-allele carriers revealed an overlapping spectrum of schizophrenia-relevant endophenotypes. Most prominently, early information processing and higher cognitive functions appear to be associated with TCF4 risk genotypes. Moreover, a recent human study unravelled gene × environment interactions between TCF4 risk alleles and smoking behaviour that were specifically associated with disrupted early information processing. Taken together, TCF4 is considered as an integrator ('hub') of several bHLH networks controlling critical steps of various developmental, and, possibly, plasticity-related transcriptional programs in the CNS and changes of TCF4 expression also appear to affect brain networks important for information processing. Consequently, these findings support the neurodevelopmental hypothesis of schizophrenia and provide a basis for identifying the underlying molecular mechanisms.

Effects of long-term averaging of quantitative blood pressure traits on the detection of genetic associations

Blood pressure (BP) is a heritable, quantitative trait with intraindividual variability and susceptibility to measurement error. Genetic studies of BP generally use single-visit measurements and thus cannot remove variability occurring over months or years. We leveraged the idea that averaging BP measured across time would improve phenotypic accuracy and thereby increase statistical power to detect genetic associations. We studied systolic BP (SBP), diastolic BP (DBP), mean arterial pressure (MAP), and pulse pressure (PP) averaged over multiple years in 46,629 individuals of European ancestry. We identified 39 trait-variant associations across 19 independent loci (p < 5 × 10(-8)); five associations (in four loci) uniquely identified by our LTA analyses included those of SBP and MAP at 2p23 (rs1275988, near KCNK3), DBP at 2q11.2 (rs7599598, in FER1L5), and PP at 6p21 (rs10948071, near CRIP3) and 7p13 (rs2949837, near IGFBP3). Replication analyses conducted in cohorts with single-visit BP data showed positive replication of associations and a nominal association (p < 0.05). We estimated a 20% gain in statistical power with long-term average (LTA) as compared to single-visit BP association studies. Using LTA analysis, we identified genetic loci influencing BP. LTA might be one way of increasing the power of genetic associations for continuous traits in extant samples for other phenotypes that are measured serially over time.

Overlap of expression quantitative trait loci (eQTL) in human brain and blood

BACKGROUND

Expression quantitative trait loci (eQTL) are genomic regions regulating RNA transcript expression levels. Genome-wide Association Studies (GWAS) have identified many variants, often in non-coding regions, with unknown functions and eQTL provide a possible mechanism by which these variants may influence observable phenotypes. Limited access and availability of tissues such as brain has led to the use of blood as a substitute for eQTL analyses.

METHODS

Here, we evaluate the overlap of eQTL reported in published studies conducted in blood and brain tissues to assess the utility of blood as an alternative to brain tissue in the study of neurological and psychiatric conditions. Expression QTL results from eight published brain studies were compared to blood eQTL identified in from a meta-analysis involving 5,311 individuals. We accounted for differences in SNP platforms and study design by using SNP proxies in high linkage disequilibrium with reported eQTL. The degree of overlap between studies was calculated by ascertaining if an eQTL identified in one study was also identified in the other study.

RESULTS

The percentage of eQTL overlapping for brain and blood expression after adjusting for differences in sample size ranged from 13 - 23% (mean 19.2%). Amongst pairs of brain studies eQTL overlap ranged from 0 - 35%, with higher degrees of overlap found for studies using expression data collected from the same brain region.

CONCLUSION

Our results suggest that whenever possible tissue specific to the pathophysiology of the disease being studied should be used for transcription analysis.

RNA-Seq optimization with eQTL gold standards

BACKGROUND

RNA-Sequencing (RNA-Seq) experiments have been optimized for library preparation, mapping, and gene expression estimation. These methods, however, have revealed weaknesses in the next stages of analysis of differential expression, with results sensitive to systematic sample stratification or, in more extreme cases, to outliers. Further, a method to assess normalization and adjustment measures imposed on the data is lacking.

RESULTS

To address these issues, we utilize previously published eQTLs as a novel gold standard at the center of a framework that integrates DNA genotypes and RNA-Seq data to optimize analysis and aid in the understanding of genetic variation and gene expression. After detecting sample contamination and sequencing outliers in RNA-Seq data, a set of previously published brain eQTLs was used to determine if sample outlier removal was appropriate. Improved replication of known eQTLs supported removal of these samples in downstream analyses. eQTL replication was further employed to assess normalization methods, covariate inclusion, and gene annotation. This method was validated in an independent RNA-Seq blood data set from the GTEx project and a tissue-appropriate set of eQTLs. eQTL replication in both data sets highlights the necessity of accounting for unknown covariates in RNA-Seq data analysis.

CONCLUSION

As each RNA-Seq experiment is unique with its own experiment-specific limitations, we offer an easily-implementable method that uses the replication of known eQTLs to guide each step in one's data analysis pipeline. In the two data sets presented herein, we highlight not only the necessity of careful outlier detection but also the need to account for unknown covariates in RNA-Seq experiments.

MicroRNAs as the cause of schizophrenia in 22q11.2 deletion carriers, and possible implications for idiopathic disease: a mini-review

The 22q11.2 deletion is the strongest known genetic risk factor for schizophrenia. Research has implicated microRNA-mediated dysregulation in 22q11.2 deletion syndrome (22q11.2DS) schizophrenia-risk. Primary candidate genes are DGCR8 (DiGeorge syndrome critical region gene 8), which encodes a component of the microprocessor complex essential for microRNA biogenesis, and MIR185, which encodes microRNA 185. Mouse models of 22q11.2DS have demonstrated alterations in brain microRNA biogenesis, and that DGCR8 haploinsufficiency may contribute to these alterations, e.g., via down-regulation of a specific microRNA subset. miR-185 was the top-scoring down-regulated microRNA in both the prefrontal cortex and the hippocampus, brain areas which are the key foci of schizophrenia research. This reduction in miR-185 expression contributed to dendritic and spine development deficits in hippocampal neurons. In addition, miR-185 has two validated targets (RhoA, Cdc42), both of which have been associated with altered expression levels in schizophrenia. These combined data support the involvement of miR-185 and its down-stream pathways in schizophrenia. This review summarizes evidence implicating microRNA-mediated dysregulation in schizophrenia in both 22q11.2DS-related and idiopathic cases.

Genetics of human gene expression

A steadily growing number of studies have identified and characterized expression quantitative trait loci (eQTLs) in human cell-lines, primary cells, and tissues. This class of variation has been shown to play a role in complex traits, including disease. Here, we discuss how eQTLs have the potential to accelerate discovery of disease genes and functional mechanisms underlying complex traits. We discuss how context-specificity of eQTLs is being characterized at an unprecedented scale and breadth, and how this both informs on the intricacy of human genome function, and has important ramifications for elucidating function of genetic variants of interest, particularly for those contributing to disease.

Differential network analysis reveals genetic effects on catalepsy modules

We performed short-term bi-directional selective breeding for haloperidol-induced catalepsy, starting from three mouse populations of increasingly complex genetic structure: an F₂ intercross, a heterogeneous stock (HS) formed by crossing four inbred strains (HS4) and a heterogeneous stock (HS-CC) formed from the inbred strain founders of the Collaborative Cross (CC). All three selections were successful, with large differences in haloperidol response emerging within three generations. Using a custom differential network analysis procedure, we found that gene coexpression patterns changed significantly; importantly, a number of these changes were concordant across genetic backgrounds. In contrast, absolute gene-expression changes were modest and not concordant across genetic backgrounds, in spite of the large and similar phenotypic differences. By inferring strain contributions from the parental lines, we are able to identify significant differences in allelic content between the selected lines concurrent with large changes in transcript connectivity. Importantly, this observation implies that genetic polymorphisms can affect transcript and module connectivity without large changes in absolute expression levels. We conclude that, in this case, selective breeding acts at the subnetwork level, with the same modules but not the same transcripts affected across the three selections.

Human-specific histone methylation signatures at transcription start sites in prefrontal neurons

Cognitive abilities and disorders unique to humans are thought to result from adaptively driven changes in brain transcriptomes, but little is known about the role of cis-regulatory changes affecting transcription start sites (TSS). Here, we mapped in human, chimpanzee, and macaque prefrontal cortex the genome-wide distribution of histone H3 trimethylated at lysine 4 (H3K4me3), an epigenetic mark sharply regulated at TSS, and identified 471 sequences with human-specific enrichment or depletion. Among these were 33 loci selectively methylated in neuronal but not non-neuronal chromatin from children and adults, including TSS at DPP10 (2q14.1), CNTN4 and CHL1 (3p26.3), and other neuropsychiatric susceptibility genes. Regulatory sequences at DPP10 and additional loci carried a strong footprint of hominid adaptation, including elevated nucleotide substitution rates and regulatory motifs absent in other primates (including archaic hominins), with evidence for selective pressures during more recent evolution and adaptive fixations in modern populations. Chromosome conformation capture at two neurodevelopmental disease loci, 2q14.1 and 16p11.2, revealed higher order chromatin structures resulting in physical contact of multiple human-specific H3K4me3 peaks spaced 0.5-1 Mb apart, in conjunction with a novel cis-bound antisense RNA linked to Polycomb repressor proteins and downregulated DPP10 expression. Therefore, coordinated epigenetic regulation via newly derived TSS chromatin could play an important role in the emergence of human-specific gene expression networks in brain that contribute to cognitive functions and neurological disease susceptibility in modern day humans.

Association of TLR5 sequence variants and mRNA level with cytokine transcription in pigs

The Toll-like receptor 5 (TLR5) plays a crucial role in host defense against flagellated bacteria by recognizing flagellin. Accumulating evidence suggests that single nucleotide polymorphisms (SNPs) in TLR5 have an effect on flagellin recognition and are associated with susceptibility/resistance to disease. In this study, we analyzed association of SNPs, including c.834T>G, c.1065T>C, c.1205C>T, c.1246A>T, c.1269G>A, and c.1398C>T, as well as mRNA level of TLR5 with the abundance of transcripts of cytokines in pigs. SNPs c.1246A>T and c.1269G>A were significantly associated with the transcript abundance of interleukin (IL)-2, and SNPs c.834T>G and c.1398C>T with IL-10 (P < 0.05); the haplotypes showed a tendency to affect the transcript abundance of IL-10 (P = 0.0660) and significantly associated with the transcription of TLR5 (P < 0.01); the abundance of transcripts of TLR5 and IL-10 were strongly correlated (P < 0.01). The results indicated that the SNPs, associated with the transcript abundance of cytokines, were related to immune responsiveness mediated by cytokine, which, in turn, would have a role in pig breeding for disease resistance. Furthermore, the positive correlation between the abundance of TLR5 and IL10 suggest a link between TLR5 activation and IL-10 expression in porcine.