Precision medicine for suicidality: from universality to subtypes and personalization.

Suicide remains a clear, present and increasing public health problem, despite being a potentially preventable tragedy. Its incidence is particularly high in people with overt or un(der)diagnosed psychiatric disorders. Objective and precise identification of individuals at risk, ways of monitoring response to treatments and novel preventive therapeutics need to be discovered, employed and widely deployed. We sought to investigate whether blood gene expression biomarkers for suicide (that is, a 'liquid biopsy' approach) can be identified that are more universal in nature, working across psychiatric diagnoses and genders, using larger cohorts than in previous studies. Such markers may reflect and/or be a proxy for the core biology of suicide. We were successful in this endeavor, using a comprehensive stepwise approach, leading to a wealth of findings. Steps 1, 2 and 3 were discovery, prioritization and validation for tracking suicidality, resulting in a Top Dozen list of candidate biomarkers comprising the top biomarkers from each step, as well as a larger list of 148 candidate biomarkers that survived Bonferroni correction in the validation step. Step 4 was testing the Top Dozen list and Bonferroni biomarker list for predictive ability for suicidal ideation (SI) and for future hospitalizations for suicidality in independent cohorts, leading to the identification of completely novel predictive biomarkers (such as CLN5 and AK2), as well as reinforcement of ours and others previous findings in the field (such as SLC4A4 and SKA2). Additionally, we examined whether subtypes of suicidality can be identified based on mental state at the time of high SI and identified four potential subtypes: high anxiety, low mood, combined and non-affective (psychotic). Such subtypes may delineate groups of individuals that are more homogenous in terms of suicidality biology and behavior. We also studied a more personalized approach, by psychiatric diagnosis and gender, with a focus on bipolar males, the highest risk group. Such a personalized approach may be more sensitive to gender differences and to the impact of psychiatric co-morbidities and medications. We compared testing the universal biomarkers in everybody versus testing by subtypes versus personalized by gender and diagnosis, and show that the subtype and personalized approaches permit enhanced precision of predictions for different universal biomarkers. In particular, LHFP appears to be a strong predictor for suicidality in males with depression. We also directly examined whether biomarkers discovered using male bipolars only are better predictors in a male bipolar independent cohort than universal biomarkers and show evidence for a possible advantage of personalization. We identified completely novel biomarkers (such as SPTBN1 and C7orf73), and reinforced previously known biomarkers (such as PTEN and SAT1). For diagnostic ability testing purposes, we also examined as predictors phenotypic measures as apps (for suicide risk (CFI-S, Convergent Functional Information for Suicidality) and for anxiety and mood (SASS, Simplified Affective State Scale)) by themselves, as well as in combination with the top biomarkers (the combination being our a priori primary endpoint), to provide context and enhance precision of predictions. We obtained area under the curves of 90% for SI and 77% for future hospitalizations in independent cohorts. Step 5 was to look for mechanistic understanding, starting with examining evidence for the Top Dozen and Bonferroni biomarkers for involvement in other psychiatric and non-psychiatric disorders, as a mechanism for biological predisposition and vulnerability. The biomarkers we identified also provide a window towards understanding the biology of suicide, implicating biological pathways related to neurogenesis, programmed cell death and insulin signaling from the universal biomarkers, as well as mTOR signaling from the male bipolar biomarkers. In particular, HTR2A increase coupled with ARRB1 and GSK3B decreases in expression in suicidality may provide a synergistic mechanistical corrective target, as do SLC4A4 increase coupled with AHCYL1 and AHCYL2 decrease. Step 6 was to move beyond diagnostics and mechanistical risk assessment, towards providing a foundation for personalized therapeutics. Items scored positive in the CFI-S and subtypes identified by SASS in different individuals provide targets for personalized (psycho)therapy. Some individual biomarkers are targets of existing drugs used to treat mood disorders and suicidality (lithium, clozapine and omega-3 fatty acids), providing a means toward pharmacogenomics stratification of patients and monitoring of response to treatment. Such biomarkers merit evaluation in clinical trials. Bioinformatics drug repurposing analyses with the gene expression biosignatures of the Top Dozen and Bonferroni-validated universal biomarkers identified novel potential therapeutics for suicidality, such as ebselen (a lithium mimetic), piracetam (a nootropic), chlorogenic acid (a polyphenol) and metformin (an antidiabetic and possible longevity promoting drug). Finally, based on the totality of our data and of the evidence in the field to date, a convergent functional evidence score prioritizing biomarkers that have all around evidence (track suicidality, predict it, are reflective of biological predisposition and are potential drug targets) brought to the fore APOE and IL6 from among the universal biomarkers, suggesting an inflammatory/accelerated aging component that may be a targetable common denominator.

Towards understanding and predicting suicidality in women: biomarkers and clinical risk assessment.

Women are under-represented in research on suicidality to date. Although women have a lower rate of suicide completion than men, due in part to the less-violent methods used, they have a higher rate of suicide attempts. Our group has previously identified genomic (blood gene expression biomarkers) and clinical information (apps) predictors for suicidality in men. We now describe pilot studies in women. We used a powerful within-participant discovery approach to identify genes that change in expression between no suicidal ideation (no SI) and high suicidal ideation (high SI) states (n=12 participants out of a cohort of 51 women psychiatric participants followed longitudinally, with diagnoses of bipolar disorder, depression, schizoaffective disorder and schizophrenia). We then used a Convergent Functional Genomics (CFG) approach to prioritize the candidate biomarkers identified in the discovery step by using all the prior evidence in the field. Next, we validated for suicidal behavior the top-ranked biomarkers for SI, in a demographically matched cohort of women suicide completers from the coroner's office (n=6), by assessing which markers were stepwise changed from no SI to high SI to suicide completers. We then tested the 50 biomarkers that survived Bonferroni correction in the validation step, as well as top increased and decreased biomarkers from the discovery and prioritization steps, in a completely independent test cohort of women psychiatric disorder participants for prediction of SI (n=33) and in a future follow-up cohort of psychiatric disorder participants for prediction of psychiatric hospitalizations due to suicidality (n=24). Additionally, we examined how two clinical instruments in the form of apps, Convergent Functional Information for Suicidality (CFI-S) and Simplified Affective State Scale (SASS), previously tested in men, perform in women. The top CFI-S item distinguishing high SI from no SI states was the chronic stress of social isolation. We then showed how the clinical information apps combined with the 50 validated biomarkers into a broad predictor (UP-Suicide), our apriori primary end point, predicts suicidality in women. UP-Suicide had a receiver-operating characteristic (ROC) area under the curve (AUC) of 82% for predicting SI and an AUC of 78% for predicting future hospitalizations for suicidality. Some of the individual components of the UP-Suicide showed even better results. SASS had an AUC of 81% for predicting SI, CFI-S had an AUC of 84% and the combination of the two apps had an AUC of 87%. The top biomarker from our sequential discovery, prioritization and validation steps, BCL2, predicted future hospitalizations due to suicidality with an AUC of 89%, and the panel of 50 validated biomarkers (BioM-50) predicted future hospitalizations due to suicidality with an AUC of 94%. The best overall single blood biomarker for predictions was PIK3C3 with an AUC of 65% for SI and an AUC of 90% for future hospitalizations. Finally, we sought to understand the biology of the biomarkers. BCL2 and GSK3B, the top CFG scoring validated biomarkers, as well as PIK3C3, have anti-apoptotic and neurotrophic effects, are decreased in expression in suicidality and are known targets of the anti-suicidal mood stabilizer drug lithium, which increases their expression and/or activity. Circadian clock genes were overrepresented among the top markers. Notably, PER1, increased in expression in suicidality, had an AUC of 84% for predicting future hospitalizations, and CSNK1A1, decreased in expression, had an AUC of 96% for predicting future hospitalizations. Circadian clock abnormalities are related to mood disorder, and sleep abnormalities have been implicated in suicide. Docosahexaenoic acid signaling was one of the top biological pathways overrepresented in validated biomarkers, which is of interest given the potential therapeutic and prophylactic benefits of omega-3 fatty acids. Some of the top biomarkers from the current work in women showed co-directionality of change in expression with our previous work in men, whereas others had changes in opposite directions, underlying the issue of biological context and differences in suicidality between the two genders. With this study, we begin to shed much needed light in the area of female suicidality, identify useful objective predictors and help understand gender commonalities and differences. During the conduct of the study, one participant committed suicide. In retrospect, when the analyses were completed, her UP-Suicide risk prediction score was at the 100 percentile of all participants tested.

Understanding and predicting suicidality using a combined genomic and clinical risk assessment approach.

Worldwide, one person dies every 40 seconds by suicide, a potentially preventable tragedy. A limiting step in our ability to intervene is the lack of objective, reliable predictors. We have previously provided proof of principle for the use of blood gene expression biomarkers to predict future hospitalizations due to suicidality, in male bipolar disorder participants. We now generalize the discovery, prioritization, validation, and testing of such markers across major psychiatric disorders (bipolar disorder, major depressive disorder, schizoaffective disorder, and schizophrenia) in male participants, to understand commonalities and differences. We used a powerful within-participant discovery approach to identify genes that change in expression between no suicidal ideation and high suicidal ideation states (n=37 participants out of a cohort of 217 psychiatric participants followed longitudinally). We then used a convergent functional genomics (CFG) approach with existing prior evidence in the field to prioritize the candidate biomarkers identified in the discovery step. Next, we validated the top biomarkers from the prioritization step for relevance to suicidal behavior, in a demographically matched cohort of suicide completers from the coroner's office (n=26). The biomarkers for suicidal ideation only are enriched for genes involved in neuronal connectivity and schizophrenia, the biomarkers also validated for suicidal behavior are enriched for genes involved in neuronal activity and mood. The 76 biomarkers that survived Bonferroni correction after validation for suicidal behavior map to biological pathways involved in immune and inflammatory response, mTOR signaling and growth factor regulation. mTOR signaling is necessary for the effects of the rapid-acting antidepressant agent ketamine, providing a novel biological rationale for its possible use in treating acute suicidality. Similarly, MAOB, a target of antidepressant inhibitors, was one of the increased biomarkers for suicidality. We also identified other potential therapeutic targets or biomarkers for drugs known to mitigate suicidality, such as omega-3 fatty acids, lithium and clozapine. Overall, 14% of the top candidate biomarkers also had evidence for involvement in psychological stress response, and 19% for involvement in programmed cell death/cellular suicide (apoptosis). It may be that in the face of adversity (stress), death mechanisms are turned on at a cellular (apoptosis) and organismal level. Finally, we tested the top increased and decreased biomarkers from the discovery for suicidal ideation (CADM1, CLIP4, DTNA, KIF2C), prioritization with CFG for prior evidence (SAT1, SKA2, SLC4A4), and validation for behavior in suicide completers (IL6, MBP, JUN, KLHDC3) steps in a completely independent test cohort of psychiatric participants for prediction of suicidal ideation (n=108), and in a future follow-up cohort of psychiatric participants (n=157) for prediction of psychiatric hospitalizations due to suicidality. The best individual biomarker across psychiatric diagnoses for predicting suicidal ideation was SLC4A4, with a receiver operating characteristic (ROC) area under the curve (AUC) of 72%. For bipolar disorder in particular, SLC4A4 predicted suicidal ideation with an AUC of 93%, and future hospitalizations with an AUC of 70%. SLC4A4 is involved in brain extracellular space pH regulation. Brain pH has been implicated in the pathophysiology of acute panic attacks. We also describe two new clinical information apps, one for affective state (simplified affective state scale, SASS) and one for suicide risk factors (Convergent Functional Information for Suicide, CFI-S), and how well they predict suicidal ideation across psychiatric diagnoses (AUC of 85% for SASS, AUC of 89% for CFI-S). We hypothesized a priori, based on our previous work, that the integration of the top biomarkers and the clinical information into a universal predictive measure (UP-Suicide) would show broad-spectrum predictive ability across psychiatric diagnoses. Indeed, the UP-Suicide was able to predict suicidal ideation across psychiatric diagnoses with an AUC of 92%. For bipolar disorder, it predicted suicidal ideation with an AUC of 98%, and future hospitalizations with an AUC of 94%. Of note, both types of tests we developed (blood biomarkers and clinical information apps) do not require asking the individual assessed if they have thoughts of suicide, as individuals who are truly suicidal often do not share that information with clinicians. We propose that the widespread use of such risk prediction tests as part of routine or targeted healthcare assessments will lead to early disease interception followed by preventive lifestyle modifications and proactive treatment.

Influence of individual differences in disease perception on consumer response to direct-to-consumer genomic testing.

Individuals who undergo multiplex direct-to-consumer (DTC) genomic testing receive genetic risk results for multiple conditions. To date, research has not investigated the influence of individual differences in disease perceptions among consumers on testing outcomes. A total of 2037 participants received DTC genomic testing and completed baseline and follow-up surveys assessing disease perceptions and health behaviors. Participants were asked to indicate their most feared disease of those tested. Perceived seriousness and controllability of the disease via lifestyle or medical intervention were assessed. Participants most frequently reported heart attack (19.1%) and Alzheimer's disease (18.6%) as their most feared disease. Perceived seriousness and control over the feared disease both influenced response to DTC genomic testing. Greater perceived seriousness and diminished perceived control were associated with higher, but not clinically significant levels of anxiety and distress. In some cases these associations were modified by genetic risk. No significant associations were observed for diet, exercise and screening behaviors. Individual differences in disease perceptions influence psychological outcomes following DTC genomic testing. Higher perceived seriousness may make a consumer more psychologically sensitive to test results and greater perceived control may protect against adverse psychological outcomes. Findings may inform development of educational and counseling services.

A genome-wide association study of anorexia nervosa.

Anorexia nervosa (AN) is a complex and heritable eating disorder characterized by dangerously low body weight. Neither candidate gene studies nor an initial genome-wide association study (GWAS) have yielded significant and replicated results. We performed a GWAS in 2907 cases with AN from 14 countries (15 sites) and 14 860 ancestrally matched controls as part of the Genetic Consortium for AN (GCAN) and the Wellcome Trust Case Control Consortium 3 (WTCCC3). Individual association analyses were conducted in each stratum and meta-analyzed across all 15 discovery data sets. Seventy-six (72 independent) single nucleotide polymorphisms were taken forward for in silico (two data sets) or de novo (13 data sets) replication genotyping in 2677 independent AN cases and 8629 European ancestry controls along with 458 AN cases and 421 controls from Japan. The final global meta-analysis across discovery and replication data sets comprised 5551 AN cases and 21 080 controls. AN subtype analyses (1606 AN restricting; 1445 AN binge-purge) were performed. No findings reached genome-wide significance. Two intronic variants were suggestively associated: rs9839776 (P=3.01 × 10(-7)) in SOX2OT and rs17030795 (P=5.84 × 10(-6)) in PPP3CA. Two additional signals were specific to Europeans: rs1523921 (P=5.76 × 10(-)(6)) between CUL3 and FAM124B and rs1886797 (P=8.05 × 10(-)(6)) near SPATA13. Comparing discovery with replication results, 76% of the effects were in the same direction, an observation highly unlikely to be due to chance (P=4 × 10(-6)), strongly suggesting that true findings exist but our sample, the largest yet reported, was underpowered for their detection. The accrual of large genotyped AN case-control samples should be an immediate priority for the field.

Evidence for the role of EPHX2 gene variants in anorexia nervosa.

Anorexia nervosa (AN) and related eating disorders are complex, multifactorial neuropsychiatric conditions with likely rare and common genetic and environmental determinants. To identify genetic variants associated with AN, we pursued a series of sequencing and genotyping studies focusing on the coding regions and upstream sequence of 152 candidate genes in a total of 1205 AN cases and 1948 controls. We identified individual variant associations in the Estrogen Receptor-ß (ESR2) gene, as well as a set of rare and common variants in the Epoxide Hydrolase 2 (EPHX2) gene, in an initial sequencing study of 261 early-onset severe AN cases and 73 controls (P=0.0004). The association of EPHX2 variants was further delineated in: (1) a pooling-based replication study involving an additional 500 AN patients and 500 controls (replication set P=0.00000016); (2) single-locus studies in a cohort of 386 previously genotyped broadly defined AN cases and 295 female population controls from the Bogalusa Heart Study (BHS) and a cohort of 58 individuals with self-reported eating disturbances and 851 controls (combined smallest single locus P<0.01). As EPHX2 is known to influence cholesterol metabolism, and AN is often associated with elevated cholesterol levels, we also investigated the association of EPHX2 variants and longitudinal body mass index (BMI) and cholesterol in BHS female and male subjects (N=229) and found evidence for a modifying effect of a subset of variants on the relationship between cholesterol and BMI (P<0.01). These findings suggest a novel association of gene variants within EPHX2 to susceptibility to AN and provide a foundation for future study of this important yet poorly understood condition.

Pooling/bootstrap-based GWAS (pbGWAS) identifies new loci modifying the age of onset in PSEN1 p.Glu280Ala Alzheimer's disease.

The literature on GWAS (genome-wide association studies) data suggests that very large sample sizes (for example, 50,000 cases and 50,000 controls) may be required to detect significant associations of genomic regions for complex disorders such as Alzheimer's disease (AD). Because of the challenges of obtaining such large cohorts, we describe here a novel sequential strategy that combines pooling of DNA and bootstrapping (pbGWAS) in order to significantly increase the statistical power and exponentially reduce expenses. We applied this method to a very homogeneous sample of patients belonging to a unique and clinically well-characterized multigenerational pedigree with one of the most severe forms of early onset AD, carrying the PSEN1 p.Glu280Ala mutation (often referred to as E280A mutation), which originated as a consequence of a founder effect. In this cohort, we identified novel loci genome-wide significantly associated as modifiers of the age of onset of AD (CD44, rs187116, P=1.29 × 10⁻¹²; NPHP1, rs10173717, P=1.74 × 10⁻¹²; CADPS2, rs3757536, P=1.54 × 10⁻¹°; GREM2, rs12129547, P=1.69 × 10⁻¹³, among others) as well as other loci known to be associated with AD. Regions identified by pbGWAS were confirmed by subsequent individual genotyping. The pbGWAS methodology and the genes it targeted could provide important insights in determining the genetic causes of AD and other complex conditions.

Discovery and validation of blood biomarkers for suicidality.

Suicides are a leading cause of death in psychiatric patients, and in society at large. Developing more quantitative and objective ways (biomarkers) for predicting and tracking suicidal states would have immediate practical applications and positive societal implications. We undertook such an endeavor. First, building on our previous blood biomarker work in mood disorders and psychosis, we decided to identify blood gene expression biomarkers for suicidality, looking at differential expression of genes in the blood of subjects with a major mood disorder (bipolar disorder), a high-risk population prone to suicidality. We compared no suicidal ideation (SI) states and high SI states using a powerful intrasubject design, as well as an intersubject case-case design, to generate a list of differentially expressed genes. Second, we used a comprehensive Convergent Functional Genomics (CFG) approach to identify and prioritize from the list of differentially expressed gene biomarkers of relevance to suicidality. CFG integrates multiple independent lines of evidence-genetic and functional genomic data-as a Bayesian strategy for identifying and prioritizing findings, reducing the false-positives and false-negatives inherent in each individual approach. Third, we examined whether expression levels of the blood biomarkers identified by us in the live bipolar subject cohort are actually altered in the blood in an age-matched cohort of suicide completers collected from the coroner's office, and report that 13 out of the 41 top CFG scoring biomarkers (32%) show step-wise significant change from no SI to high SI states, and then to the suicide completers group. Six out of them (15%) remained significant after strict Bonferroni correction for multiple comparisons. Fourth, we show that the blood levels of SAT1 (spermidine/spermine N1-acetyltransferase 1), the top biomarker identified by us, at the time of testing for this study, differentiated future as well as past hospitalizations with suicidality, in a live cohort of bipolar disorder subjects, and exhibited a similar but weaker pattern in a live cohort of psychosis (schizophrenia/schizoaffective disorder) subjects. Three other (phosphatase and tensin homolog (PTEN), myristoylated alanine-rich protein kinase C substrate (MARCKS), and mitogen-activated protein kinase kinase kinase 3 (MAP3K3)) of the six biomarkers that survived Bonferroni correction showed similar but weaker effects. Taken together, the prospective and retrospective hospitalization data suggests SAT1, PTEN, MARCKS and MAP3K3 might be not only state biomarkers but trait biomarkers as well. Fifth, we show how a multi-dimensional approach using SAT1 blood expression levels and two simple visual-analog scales for anxiety and mood enhances predictions of future hospitalizations for suicidality in the bipolar cohort (receiver-operating characteristic curve with area under the curve of 0.813). Of note, this simple approach does not directly ask about SI, which some individuals may deny or choose not to share with clinicians. Lastly, we conducted bioinformatic analyses to identify biological pathways, mechanisms and medication targets. Overall, suicidality may be underlined, at least in part, by biological mechanisms related to stress, inflammation and apoptosis.

Enrichment of cis-regulatory gene expression SNPs and methylation quantitative trait loci among bipolar disorder susceptibility variants.

We conducted a systematic study of top susceptibility variants from a genome-wide association (GWA) study of bipolar disorder to gain insight into the functional consequences of genetic variation influencing disease risk. We report here the results of experiments to explore the effects of these susceptibility variants on DNA methylation and mRNA expression in human cerebellum samples. Among the top susceptibility variants, we identified an enrichment of cis regulatory loci on mRNA expression (eQTLs), and a significant excess of quantitative trait loci for DNA CpG methylation, hereafter referred to as methylation quantitative trait loci (mQTLs). Bipolar disorder susceptibility variants that cis regulate both cerebellar expression and methylation of the same gene are a very small proportion of bipolar disorder susceptibility variants. This finding suggests that mQTLs and eQTLs provide orthogonal ways of functionally annotating genetic variation within the context of studies of pathophysiology in brain. No lymphocyte mQTL enrichment was found, suggesting that mQTL enrichment was specific to the cerebellum, in contrast to eQTLs. Separately, we found that using mQTL information to restrict the number of single-nucleotide polymorphisms studied enhances our ability to detect a significant association. With this restriction a priori informed by the observed functional enrichment, we identified a significant association (rs12618769, P(bonferroni)<0.05) from two other GWA studies (TGen+GAIN; 2191 cases and 1434 controls) of bipolar disorder, which we replicated in an independent GWA study (WTCCC). Collectively, our findings highlight the importance of integrating functional annotation of genetic variants for gene expression and DNA methylation to advance the biological understanding of bipolar disorder.

Perceptions of genetic counseling services in direct-to-consumer personal genomic testing.

To describe consumers' perceptions of genetic counseling services in the context of direct-to-consumer personal genomic testing is the purpose of this research. Utilizing data from the Scripps Genomic Health Initiative, we assessed direct-to-consumer genomic test consumers' utilization and perceptions of genetic counseling services. At long-term follow-up, approximately 14 months post-testing, participants were asked to respond to several items gauging their interactions, if any, with a Navigenics genetic counselor, and their perceptions of those interactions. Out of 1325 individuals who completed long-term follow-up, 187 (14.1%) indicated that they had spoken with a genetic counselor. The most commonly given reason for not utilizing the counseling service was a lack of need due to the perception of already understanding one's results (55.6%). The most common reasons for utilizing the service included wanting to take advantage of a free service (43.9%) and wanting more information on risk calculations (42.2%). Among those who utilized the service, a large fraction reported that counseling improved their understanding of their results (54.5%) and genetics in general (43.9%). A relatively small proportion of participants utilized genetic counseling after direct-to-consumer personal genomic testing. Among those individuals who did utilize the service, however, a large fraction perceived it to be informative, and thus presumably beneficial.

Convergent functional genomics of schizophrenia: from comprehensive understanding to genetic risk prediction.

We have used a translational convergent functional genomics (CFG) approach to identify and prioritize genes involved in schizophrenia, by gene-level integration of genome-wide association study data with other genetic and gene expression studies in humans and animal models. Using this polyevidence scoring and pathway analyses, we identify top genes (DISC1, TCF4, MBP, MOBP, NCAM1, NRCAM, NDUFV2, RAB18, as well as ADCYAP1, BDNF, CNR1, COMT, DRD2, DTNBP1, GAD1, GRIA1, GRIN2B, HTR2A, NRG1, RELN, SNAP-25, TNIK), brain development, myelination, cell adhesion, glutamate receptor signaling, G-protein-coupled receptor signaling and cAMP-mediated signaling as key to pathophysiology and as targets for therapeutic intervention. Overall, the data are consistent with a model of disrupted connectivity in schizophrenia, resulting from the effects of neurodevelopmental environmental stress on a background of genetic vulnerability. In addition, we show how the top candidate genes identified by CFG can be used to generate a genetic risk prediction score (GRPS) to aid schizophrenia diagnostics, with predictive ability in independent cohorts. The GRPS also differentiates classic age of onset schizophrenia from early onset and late-onset disease. We also show, in three independent cohorts, two European American and one African American, increasing overlap, reproducibility and consistency of findings from single-nucleotide polymorphisms to genes, then genes prioritized by CFG, and ultimately at the level of biological pathways and mechanisms. Finally, we compared our top candidate genes for schizophrenia from this analysis with top candidate genes for bipolar disorder and anxiety disorders from previous CFG analyses conducted by us, as well as findings from the fields of autism and Alzheimer. Overall, our work maps the genomic and biological landscape for schizophrenia, providing leads towards a better understanding of illness, diagnostics and therapeutics. It also reveals the significant genetic overlap with other major psychiatric disorder domains, suggesting the need for improved nosology.

The insulin gene VNTR is associated with fasting insulin levels and development of juvenile obesity.

In millions of people, obesity leads to type 2 diabetes (T2D; also known as non-insulin-dependent diabetes mellitus). During the early stages of juvenile obesity, the increase of insulin secretion in proportion to accumulated fat balances insulin resistance and protects patients from hyperglycaemia. After several decades, however,beta-cell function deteriorates and T2D develops in approximately 20% of obese patients. In modern societies, obesity has thus become the leading risk factor for T2D (ref. 5). The factors that predispose obese patients to alteration of insulin secretion upon gaining weight remain unknown. To determine which genetic factors predispose obese patients to beta-cell dysfunction, and possibly T2D, we studied single-nucleotide polymorphisms (SNPs) in the region of the insulin gene (INS) among 615 obese children. We found that, in the early phase of obesity, alleles of the INS variable number of tandem repeat (VNTR) locus are associated with different effects of body fatness on insulin secretion. Young obese patients homozygous for class I VNTR alleles secrete more insulin than those with other genotypes.

Sex-determining genes on mouse autosomes identified by linkage analysis of C57BL/6J-YPOS sex reversal.

A powerful approach for identifying mammalian primary (gonadal) sex determination genes is the molecular genetic analyses of sex reversal conditions (that is, XX individuals with testicular tissue and XY individuals with ovarian tissue). Here we determined the number and chromosomal location of autosomal and X-linked genes that cause sex reversal in C57BL/6J (B6) mice carrying a Y chromosome of Mus domesticus poschiavinus origin (YPOS). B6 XYPOS mice develop either as females with exclusively ovarian tissue or as true hermaphrodites with ovarian and testicular tissue. In contrast, the YPOS chromosome is fully masculinizing on most other inbred strain backgrounds. B6-YPOS sex reversal appears to result from the incompatibility of the Sry (sex determining region, Y chromosome) allele carried on the YPOS chromosome with B6-derived autosomal or X-linked loci. We found strong evidence for the location of one gene, designated tda1 (testis-determining, autosomal 1), at the distal end of Chromosome (Chr) 4 and a second gene, tda2, in the central region of Chr 2. A third gene, tda3, on Chr 5 is implicated, but the evidence here is not as strong. We suggest that B6 alleles at these loci predispose XYPOS fetuses to ovarian tissue development, but no single locus or combination of loci is necessary and sufficient to cause sex reversal. The TDA proteins may regulate Sry expression or form complexes with the SRY protein to regulate other genes, or the tda genes may be activated or repressed by the SRY protein.

Exploring human biology with N-of-1 clinical trials.

Studies on humans that exploit contemporary data-intensive, high-throughput 'omic' assay technologies, such as genomics, transcriptomics, proteomics and metabolomics, have unequivocally revealed that humans differ greatly at the molecular level. These differences, which are compounded by each individual's distinct behavioral and environmental exposures, impact individual responses to health interventions such as diet and drugs. Questions about the best way to tailor health interventions to individuals based on their nuanced genomic, physiologic, behavioral, etc. profiles have motivated the current emphasis on 'precision' medicine. This review's purpose is to describe how the design and execution of N-of-1 (or personalized) multivariate clinical trials can advance the field. Such trials focus on individual responses to health interventions from a whole-person perspective, leverage emerging health monitoring technologies, and can be used to address the most relevant questions in the precision medicine era. This includes how to validate biomarkers that may indicate appropriate activity of an intervention as well as how to identify likely beneficial interventions for an individual. We also argue that multivariate N-of-1 and aggregated N-of-1 trials are ideal vehicles for advancing biomedical and translational science in the precision medicine era since the insights gained from them can not only shed light on how to treat or prevent diseases generally, but also provide insight into how to provide real-time care to the very individuals who are seeking attention for their health concerns in the first place.

Genetic regulation of catecholamine synthesis, storage and secretion in the spontaneously hypertensive rat.

Understanding catecholamine metabolism is crucial for elucidating the pathogenesis of hereditary hypertension. Here we integrated transcriptional and biochemical profiling with physiologic quantitative trait locus (eQTL and pQTL) mapping in adrenal glands of the HXB/BXH recombinant inbred (RI) strains, derived from the spontaneously hypertensive rat (SHR) and normotensive Brown Norway (BN.Lx). We found simultaneous down-regulation of five heritable transcripts in the catecholaminergic pathway in young (6 weeks) SHRs. We identified cis-acting eQTLs for Dbh, Pnmt (catecholamine biosynthesis) and Vamp1 (catecholamine secretion); enzymatic activities of Dbh and Pnmt paralleled transcripts, with pQTLs for activities mirroring eQTLs. We also detected trans-regulated expression of Vmat1 and Chga (both involved in catecholamine storage), with co-localization of these trans-eQTLs to the Pnmt locus. Pnmt re-sequencing revealed promoter polymorphisms that result in decreased response of the transfected SHR promoter to glucocorticoid, compared with BN.Lx. Of physiological pertinence, Dbh activity negatively correlated with systolic blood pressure in RI strains, whereas Pnmt activity was negatively correlated with heart rate. The finding of such cis- and trans-QTLs at an age before the onset of frank hypertension suggests that these heritable changes in biosynthetic enzyme expression represent primary genetic mechanisms for regulation of catecholamine action and blood pressure control in this widely studied model of hypertension.

Background gene expression networks significantly enhance drug response prediction by transcriptional profiling.

A central goal of gene expression studies coupled with drug response screens is to identify predictive profiles that can be exploited to stratify patients. Numerous methods have been proposed towards this end, most of them focusing on novel statistical methods and model selection techniques that attempt to uncover groups of genes, whose expression profiles are directly and robustly correlated with drug response. However, biological systems process information through the crosstalk of multiple signaling networks, whose ultimate phenotypic consequences may only be determined by the combined input of relevant interacting systems. By restricting predictive signatures to direct gene-drug correlations, biologically meaningful interactions that may serve as superior predictors are ignored. Here we demonstrate that predictive signatures, which incorporate the interaction between background gene expression patterns and individual predictive probes, can provide superior models than those that directly relate gene expression levels to pharmacological response, and thus should be more widely utilized in pharmacogenetic studies.

A genome-wide association study on common SNPs and rare CNVs in anorexia nervosa.

Anorexia nervosa (AN) is a mental illness with high mortality that most commonly afflicts adolescent female individuals. Clinical symptoms include chronic food refusal, weight loss and body image distortions. We carried out a genome-wide association study on 1033 AN cases and 3733 pediatric control subjects, all of whom were of European ancestry and were genotyped on the Illumina HumanHap610 platform (Illumina, San Diego, CA, USA). We confirmed that common single-nucleotide polymorphisms (SNPs) within OPRD1 (rs533123, P=0.0015) confer risk for AN, and obtained suggestive evidence that common SNPs near HTR1D (rs7532266, P=0.04) confer risk for restricting-type AN specifically. However, no SNPs reached genome-wide significance in our data, whereas top association signals were detected near ZNF804B, CSRP2BP, NTNG1, AKAP6 and CDH9. In parallel, we performed genome-wide analysis on copy number variations (CNVs) using the signal intensity data from the SNP arrays. We did not find evidence that AN cases have more CNVs than control subjects, nor do they have over-representation of rare or large CNVs. However, we identified several regions with rare CNVs that were only observed in AN cases, including a recurrent 13q12 deletion (1.5 Mb) disrupting SCAS in two cases, and CNVs disrupting the CNTN6/CNTN4 region in several AN cases. In conclusion, our study suggests that both common SNPs and rare CNVs may confer genetic risk to AN. These results point to intriguing genes that await further validation in independent cohorts for confirmatory roles in AN.

Identification of blood biomarkers for psychosis using convergent functional genomics.

There are to date no objective clinical laboratory blood tests for psychotic disease states. We provide proof of principle for a convergent functional genomics (CFG) approach to help identify and prioritize blood biomarkers for two key psychotic symptoms, one sensory (hallucinations) and one cognitive (delusions). We used gene expression profiling in whole blood samples from patients with schizophrenia and related disorders, with phenotypic information collected at the time of blood draw, then cross-matched the data with other human and animal model lines of evidence. Topping our list of candidate blood biomarkers for hallucinations, we have four genes decreased in expression in high hallucinations states (Fn1, Rhobtb3, Aldh1l1, Mpp3), and three genes increased in high hallucinations states (Arhgef9, Phlda1, S100a6). All of these genes have prior evidence of differential expression in schizophrenia patients. At the top of our list of candidate blood biomarkers for delusions, we have 15 genes decreased in expression in high delusions states (such as Drd2, Apoe, Scamp1, Fn1, Idh1, Aldh1l1), and 16 genes increased in high delusions states (such as Nrg1, Egr1, Pvalb, Dctn1, Nmt1, Tob2). Twenty-five of these genes have prior evidence of differential expression in schizophrenia patients. Predictive scores, based on panels of top candidate biomarkers, show good sensitivity and negative predictive value for detecting high psychosis states in the original cohort as well as in three additional cohorts. These results have implications for the development of objective laboratory tests to measure illness severity and response to treatment in devastating disorders such as schizophrenia.

Two separate, large cohorts reveal potential modifiers of age-associated variation in visual reaction time performance.

To identify potential factors influencing age-related cognitive decline and disease, we created MindCrowd. MindCrowd is a cross-sectional web-based assessment of simple visual (sv) reaction time (RT) and paired-associate learning (PAL). svRT and PAL results were combined with 22 survey questions. Analysis of svRT revealed education and stroke as potential modifiers of changes in processing speed and memory from younger to older ages (ntotal = 75,666, nwomen = 47,700, nmen = 27,966; ages 18-85 years old, mean (M)Age = 46.54, standard deviation (SD)Age = 18.40). To complement this work, we evaluated complex visual recognition reaction time (cvrRT) in the UK Biobank (ntotal = 158,249 nwomen = 89,333 nmen = 68,916; ages 40-70 years old, MAge = 55.81, SDAge = 7.72). Similarities between the UK Biobank and MindCrowd were assessed using a subset of MindCrowd (UKBb MindCrowd) selected to mirror the UK Biobank demographics (ntotal = 39,795, nwomen = 29,640, nmen = 10,155; ages 40-70 years old, MAge = 56.59, SDAge = 8.16). An identical linear model (LM) was used to assess both cohorts. Analyses revealed similarities between MindCrowd and the UK Biobank across most results. Divergent findings from the UK Biobank included (1) a first-degree family history of Alzheimer's disease (FHAD) was associated with longer cvrRT. (2) Men with the least education were associated with longer cvrRTs comparable to women across all educational attainment levels. Divergent findings from UKBb MindCrowd included more education being associated with shorter svRTs and a history of smoking with longer svRTs from younger to older ages.

Genetic regulation of commitment to interleukin 4 production by a CD4(+) T cell-intrinsic mechanism.

The dysregulated expression of interleukin 4 (IL-4) can have deleterious effects on the outcome of infectious and allergic diseases. Despite this, the mechanisms by which naive T cells commit to IL-4 expression during differentiation into mature effector cells remain incompletely defined. As compared to cells from most strains of mice, activated CD4(+) T cells from BALB mice show a bias towards IL-4 production and T helper 2 commitment in vitro and in vivo. Here, we show that this bias arises not from an increase in the amount of IL-4 produced per cell, but rather from an increase in the proportion of CD4(+) T cells that commit to IL-4 expression. This strain-specific difference in commitment was independent of signals mediated via the IL-4 receptor and hence occurred upstream of potential autoregulatory effects of IL-4. Segregation analysis of the phenotype in an experimental backcross cohort implicated a polymorphic locus on chromosome 16. Consistent with a role in differentiation, expression of the phenotype was CD4(+) T cell intrinsic and was evident as early as 16 h after the activation of naive T cells. Probabilistic gene activation is proposed as a T cell-intrinsic mechanism capable of modulating the proportion of naive T cells that commit to IL-4 production.