Expansion of Schizophrenia Gene Network Knowledge Using Machine Learning Selected Signals From Dorsolateral Prefrontal Cortex and Amygdala RNA-seq Data.

It is widely accepted, given the complex nature of schizophrenia (SCZ) gene networks, that a few or a small number of genes are unlikely to represent the underlying functional pathways responsible for SCZ pathogenesis. Several studies from large cohorts have been performed to search for key SCZ network genes using different analytical approaches, such as differential expression tests, genome-wide association study (GWAS), copy number variations, and differential methylations, or from the analysis of mutations residing in the coding regions of the genome. However, only a small portion (<10%) of candidate genes identified in these studies were considered SCZ disease-associated genes in SCZ pathways. RNA sequencing (RNA-seq) has been a powerful method to detect functional signals. In this study, we used RNA-seq data from the dorsolateral prefrontal cortex (DLPFC) from 254 individuals and RNA-seq data from the amygdala region from 46 individuals. Analysis was performed using machine learning methods, including random forest and factor analysis, to prioritize the numbers of genes from previous SCZ studies. For genes most differentially expressed between SCZ and healthy controls, 18 were added to known SCZ-associated pathways. These include three genes (GNB2, ITPR1, and PLCB2) for the glutamatergic synapse pathway, six genes (P2RX6, EDNRB, GHR, GRID2, TSPO, and S1PR1) for neuroactive ligand-receptor interaction, eight genes (CAMK2G, MAP2K1, RAF1, PDE3A, RRAS2, VAV1, ATP1B2, and GLI3) for the cAMP signaling pathway, and four genes (GNB2, CAMK2G, ITPR1, and PLCB2) for the dopaminergic synapse pathway. Besides the previously established pathways, 103 additional gene interactions were expanded to SCZ-associated networks, which were shared among both the DLPFC and amygdala regions. The novel knowledge of molecular targets gained from this study brings opportunities for a more complete picture of the SCZ pathogenesis. A noticeable fact is that hub genes, in the expanded networks, are not necessary differentially expressed or containing hotspots from GWAS studies, indicating that individual methods, such as differential expression tests, are not enough to identify the underlying SCZ pathways and that more integrative analysis is required to unfold the pathobiology of SCZ.

Heritability and genome-wide association study of benign prostatic hyperplasia (BPH) in the eMERGE network.

Benign prostatic hyperplasia (BPH) results in a significant public health burden due to the morbidity caused by the disease and many of the available remedies. As much as 70% of men over 70 will develop BPH. Few studies have been conducted to discover the genetic determinants of BPH risk. Understanding the biological basis for this condition may provide necessary insight for development of novel pharmaceutical therapies or risk prediction. We have evaluated SNP-based heritability of BPH in two cohorts and conducted a genome-wide association study (GWAS) of BPH risk using 2,656 cases and 7,763 controls identified from the Electronic Medical Records and Genomics (eMERGE) network. SNP-based heritability estimates suggest that roughly 60% of the phenotypic variation in BPH is accounted for by genetic factors. We used logistic regression to model BPH risk as a function of principal components of ancestry, age, and imputed genotype data, with meta-analysis performed using METAL. The top result was on chromosome 22 in SYN3 at rs2710383 (p-value = 4.6 × 10-7; Odds Ratio = 0.69, 95% confidence interval = 0.55-0.83). Other suggestive signals were near genes GLGC, UNCA13, SORCS1 and between BTBD3 and SPTLC3. We also evaluated genetically-predicted gene expression in prostate tissue. The most significant result was with increasing predicted expression of ETV4 (chr17; p-value = 0.0015). Overexpression of this gene has been associated with poor prognosis in prostate cancer. In conclusion, although there were no genome-wide significant variants identified for BPH susceptibility, we present evidence supporting the heritability of this phenotype, have identified suggestive signals, and evaluated the association between BPH and genetically-predicted gene expression in prostate.

Non-coding RNA dysregulation in the amygdala region of schizophrenia patients contributes to the pathogenesis of the disease.

Schizophrenia (SCZ) is a neuropsychiatric disorder with a complex genetic etiology. The redundancy of the gene networks underlying SCZ indicates that many gene combinations have the potential to cause a system dysfunction that can manifest as SCZ or a related neurodevelopmental disorder. Recent studies show that small non-coding microRNA (miRNA) and long non-coding RNA (lncRNA) are important factors in shaping these networks and are dynamically regulated by neuronal activation. We investigated the genome-wide transcription profiles of 46 human amygdala samples obtained from 22 SCZ patients and 24 healthy controls. Using RNA sequencing (RNA-seq), we determined lncRNA expression levels in all samples and generated miRNA profiles for 27 individuals (13 cases and 14 controls). Previous studies have identified differentially expressed miRNAs in SCZ, including miR-132, miR-212, and miR-34a/miR-34c. Here we report differential expression of a novel miRNA, miR1307, in SCZ. Notably, miR1307 maps to a locus previously associated with SCZ through GWAS. Additionally, one lncRNA that was overexpressed in SCZ, AC005009.2, also maps to a region previously associated with SCZ based on GWAS and overlapped SCZ-related genes. The results were replicated in a large independent data set of 254 dorsolateral prefrontal cortex samples from the CommonMind consortium. Taken together, these results suggest that miRNA and lncRNAs are important contributors to the pathogenesis of SCZ.

Erratum: A genome-wide association study of anorexia nervosa suggests a risk locus implicated in dysregulated leptin signaling.

A correction to this article has been published and is linked from the HTML version of this paper. The error has been fixed in the paper.

A genome-wide association study of anorexia nervosa suggests a risk locus implicated in dysregulated leptin signaling.

We conducted a genome-wide association study (GWAS) of anorexia nervosa (AN) using a stringently defined phenotype. Analysis of phenotypic variability led to the identification of a specific genetic risk factor that approached genome-wide significance (rs929626 in EBF1 (Early B-Cell Factor 1); P = 2.04 × 10-7; OR = 0.7; 95% confidence interval (CI) = 0.61-0.8) with independent replication (P = 0.04), suggesting a variant-mediated dysregulation of leptin signaling may play a role in AN. Multiple SNPs in LD with the variant support the nominal association. This demonstrates that although the clinical and etiologic heterogeneity of AN is universally recognized, further careful sub-typing of cases may provide more precise genomic signals. In this study, through a refinement of the phenotype spectrum of AN, we present a replicable GWAS signal that is nominally associated with AN, highlighting a potentially important candidate locus for further investigation.

The psychosis spectrum in a young U.S. community sample: findings from the Philadelphia Neurodevelopmental Cohort.

Little is known about the occurrence and predictors of the psychosis spectrum in large non-clinical community samples of U.S. youths. We aimed to bridge this gap through assessment of psychosis spectrum symptoms in the Philadelphia Neurodevelopmental Cohort, a collaborative investigation of clinical and neurobehavioral phenotypes in a prospectively accrued cohort of youths, funded by the National Institute of Mental Health. Youths (age 11-21; N=7,054) and collateral informants (caregiver/legal guardian) were recruited through the Children's Hospital of Philadelphia and administered structured screens of psychosis spectrum symptoms, other major psychopathology domains, and substance use. Youths were also administered a computerized neurocognitive battery assessing five neurobehavioral domains. Predictors of psychosis spectrum status in physically healthy participants (N=4,848) were examined using logistic regression. Among medically healthy youths, 3.7% reported threshold psychotic symptoms (delusions and/or hallucinations). An additional 12.3% reported significant sub-psychotic positive symptoms, with odd/unusual thoughts and auditory perceptions, followed by reality confusion, being the most discriminating and widely endorsed attenuated symptoms. A minority of youths (2.3%) endorsed subclinical negative/disorganized symptoms in the absence of positive symptoms. Caregivers reported lower symptom levels than their children. Male gender, younger age, and non-European American ethnicity were significant predictors of spectrum status. Youths with spectrum symptoms had reduced accuracy across neurocognitive domains, reduced global functioning, and increased odds of depression, anxiety, behavioral disorders, substance use and suicidal ideation. These findings have public health relevance for prevention and early intervention.

Examination of genetic variants influencing lipid traits in pediatric populations.

Previous large-scale genome-wide association studies in adult populations have implicated ∽100 loci in determining high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol, or triglyceride levels. However, whether these loci also contribute to variations of lipid traits in pediatric populations remain unknown. Here we assayed a population of Philadelphia children by high-density single nucleotide polymorphism arrays, and performed association analysis on lipid traits ascertained from lipid measurements stored in electronic medical records. We examined previously reported lipid trait associations, and found that most of them show identical direction of association in our pediatric cohorts, including genome-wide significant association on cholesteryl ester transfer protein with HDL-C levels (rs3764261, P = 2.1 × 10(-8)) and other significant associations on oxysterol-binding protein-like protein 7, low-density lipoprotein receptor-related protein 4 and low-density lipoprotein receptor-related protein 1. Additionally, we identified suggestive association on low-density lipoprotein receptor-related protein 1B with HDL-C levels (rs17736712, P = 2.1 × 10(-7)), but this signal is not supported by previous meta-analysis on adult cohorts. Finally, we examined rare copy number variants and identified deletions encompassing tetratricopeptide repeat domain 39B in two children with extreme lipid measures. Our results highlight the commonalities and differences of genetic components in determining lipid traits in pediatric versus adult populations. Furthermore, our study demonstrates the unique utility of automated information retrieval from electronic medical records in facilitating the identification of genotype-phenotype associations.