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.

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.

Predictors of adverse psychological experiences surrounding genome-wide profiling for disease risk.

This study aimed to identify predictors of adverse psychological experiences among direct-to-consumer (DTC) genomic test consumers. We performed a secondary analysis on data from the Scripps Genomic Health Initiative (SGHI), which studied 2037 individuals tested with commercially available tests yielding personalized risk estimates for 23 common, genetically complex diseases. As part of the original study, the participants completed baseline and follow-up survey measures assessing demographics, personal and family health history, attitudes toward genetic testing, anxiety (State-Trait Anxiety Inventory (STAI)), test-related distress (Impact of Event Scale-Revised (IES-R)), and reactions to receipt of results. To further describe the participants who had an adverse psychological outcome, this secondary analysis defined two different variables ("distress response" and "psychologically sensitive participants") and examined their relationship to various demographic variables and other survey responses. One hundred thirty participants (6.4%) were defined as having a "distress response" to receipt of results based on changes in STAI and/or IES. Four hundred thirty-one participants (21.2%) were defined as being "psychologically sensitive" based on high STAI scores both pre- and post-receipt of results. For psychologically sensitive subjects, younger age emerged as a predictor (p < 0.0005). Family history and personal history were only significant predictors for Alzheimer's disease in the psychologically sensitive participants (p = .03) and restless leg syndrome in those with a distress response (p = .03). Psychologically sensitive participants were more likely to indicate a number of pre-test concerns than were controls, but neither group of participants were any more likely to follow up with their physician or a free genetic counseling service after the return of results.

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.

Integrated Genomic Medicine: A Paradigm for Rare Diseases and Beyond.

Individualized medicine, or the tailoring of therapeutic interventions to a patient's unique genetic, biochemical, physiological, exposure and behavioral profile, has been enhanced, if not enabled, by modern biomedical technologies such as high-throughput DNA sequencing platforms, induced pluripotent stem cell assays, biomarker discovery protocols, imaging modalities, and wireless monitoring devices. Despite successes in the isolated use of these technologies, however, it is arguable that their combined and integrated use in focused studies of individual patients is the best way to not only tailor interventions for those patients, but also shed light on treatment strategies for patients with similar conditions. This is particularly true for individuals with rare diseases since, by definition, they will require study without recourse to other individuals, or at least without recourse to many other individuals. Such integration and focus will require new biomedical scientific paradigms and infrastructure, including the creation of databases harboring study results, the formation of dedicated multidisciplinary research teams and new training programs. We consider the motivation and potential for such integration, point out areas in need of improvement, and argue for greater emphasis on improving patient health via technological innovations, not merely improving the technologies themselves. We also argue that the paradigm described can, in theory, be extended to the study of individuals with more common diseases.

Single nucleotide polymorphisms in the REG-CTNNA2 region of chromosome 2 and NEIL3 associated with impulsivity in a Native American sample.

Impulsivity is a multi-faceted construct that, while characterized by a set of correlated dimensions, is centered around a core definition that involves acting suddenly in an unplanned manner without consideration for the consequences of such behavior. Several psychiatric disorders include impulsivity as a criterion, and thus it has been suggested that it may link a number of different behavioral disorders, including substance abuse. Native Americans (NA) experience some of the highest rates of substance abuse of all the US ethnic groups. The described analyses used data from a low-coverage whole genome sequence scan to conduct a genome-wide association study (GWAS) of an impulsivity phenotype in an American Indian community sample (n = 658). Demographic and clinical information were obtained using a semi-structured interview. Impulsivity was assessed using a scale derived from the Maudsley personality inventory that combines both novelty seeking and lack of planning items. The impulsivity score was tested for association with each variant adjusted for demographic variables, and corrected for ancestry and kinship, using emmax. Simulations were conducted to calculate empirical P-values. Genome-wide significant findings were observed for a variant 50-kb upstream from catenin cadherin-associated protein, alpha 2 (CTNNA2), a neuronal-specific catenin, in the REG gene cluster. A meta-analysis of GWAS had previously identified common variants in CTNNA2 as being associated with excitement seeking. A second locus upstream of nei endonuclease VIII-like 3 (NEIL3) on chromosome 4 also achieved genome-wide significance. The association between sequence variants in these regions suggests their potential roles in the genetic regulation of this phenotype in this population.

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.

Dysregulation of soluble epoxide hydrolase and lipidomic profiles in anorexia nervosa.

Individuals with anorexia nervosa (AN) restrict eating and become emaciated. They tend to have an aversion to foods rich in fat. Because epoxide hydrolase 2 (EPHX2) was identified as a novel AN susceptibility gene, and because its protein product, soluble epoxide hydrolase (sEH), converts bioactive epoxides of polyunsaturated fatty acid (PUFA) to the corresponding diols, lipidomic and metabolomic targets of EPHX2 were assessed to evaluate the biological functions of EPHX2 and their role in AN. Epoxide substrates of sEH and associated oxylipins were measured in ill AN, recovered AN and gender- and race-matched controls. PUFA and oxylipin markers were tested as potential biomarkers for AN. Oxylipin ratios were calculated as proxy markers of in vivo sEH activity. Several free- and total PUFAs were associated with AN diagnosis and with AN recovery. AN displayed elevated n-3 PUFAs and may differ from controls in PUFA elongation and desaturation processes. Cytochrome P450 pathway oxylipins from arachidonic acid, linoleic acid, alpha-linolenic acid and docosahexaenoic acid PUFAs are associated with AN diagnosis. The diol:epoxide ratios suggest the sEH activity is higher in AN compared with controls. Multivariate analysis illustrates normalization of lipidomic profiles in recovered ANs. EPHX2 influences AN risk through in vivo interaction with dietary PUFAs. PUFA composition and concentrations as well as sEH activity may contribute to the pathogenesis and prognosis of AN. Our data support the involvement of EPHX2-associated lipidomic and oxylipin dysregulations in AN, and reveal their potential as biomarkers to assess responsiveness to future intervention or treatment.

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.

Genetic risk prediction and neurobiological understanding of alcoholism.

We have used a translational Convergent Functional Genomics (CFG) approach to discover genes involved in alcoholism, by gene-level integration of genome-wide association study (GWAS) data from a German alcohol dependence cohort with other genetic and gene expression data, from human and animal model studies, similar to our previous work in bipolar disorder and schizophrenia. A panel of all the nominally significant P-value SNPs in the top candidate genes discovered by CFG  (n=135 genes, 713 SNPs) was used to generate a genetic  risk prediction score (GRPS), which showed a trend towards significance (P=0.053) in separating  alcohol dependent individuals from controls in an independent German test cohort. We then validated and prioritized our top findings from this discovery work, and subsequently tested them in three independent cohorts, from two continents. A panel of all the nominally significant P-value single-nucleotide length polymorphisms (SNPs) in the top candidate genes discovered by CFG (n=135 genes, 713 SNPs) were used to generate a Genetic Risk Prediction Score (GRPS), which showed a trend towards significance (P=0.053) in separating alcohol-dependent individuals from controls in an independent German test cohort. In order to validate and prioritize the key genes that drive behavior without some of the pleiotropic environmental confounds present in humans, we used a stress-reactive animal model of alcoholism developed by our group, the D-box binding protein (DBP) knockout mouse, consistent with the surfeit of stress theory of addiction proposed by Koob and colleagues. A much smaller panel (n=11 genes, 66 SNPs) of the top CFG-discovered genes for alcoholism, cross-validated and prioritized by this stress-reactive animal model showed better predictive ability in the independent German test cohort (P=0.041). The top CFG scoring gene for alcoholism from the initial discovery step, synuclein alpha (SNCA) remained the top gene after the stress-reactive animal model cross-validation. We also tested this small panel of genes in two other independent test cohorts from the United States, one with alcohol dependence (P=0.00012) and one with alcohol abuse (a less severe form of alcoholism; P=0.0094). SNCA by itself was able to separate alcoholics from controls in the alcohol-dependent cohort (P=0.000013) and the alcohol abuse cohort (P=0.023). So did eight other genes from the panel of 11 genes taken individually, albeit to a lesser extent and/or less broadly across cohorts. SNCA, GRM3 and MBP survived strict Bonferroni correction for multiple comparisons. Taken together, these results suggest that our stress-reactive DBP animal model helped to validate and prioritize from the CFG-discovered genes some of the key behaviorally relevant genes for alcoholism. These genes fall into a series of biological pathways involved in signal transduction, transmission of nerve impulse (including myelination) and cocaine addiction. Overall, our work provides leads towards a better understanding of illness, diagnostics and therapeutics, including treatment with omega-3 fatty acids. We also examined the overlap between the top candidate genes for alcoholism from this work and the top candidate genes for bipolar disorder, schizophrenia, anxiety from previous CFG analyses conducted by us, as well as cross-tested genetic risk predictions. This revealed the significant genetic overlap with other major psychiatric disorder domains, providing a basis for comorbidity and dual diagnosis, and placing alcohol use in the broader context of modulating the mental landscape.

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.

Genome-wide association study of shared components of reading disability and language impairment.

Written and verbal languages are neurobehavioral traits vital to the development of communication skills. Unfortunately, disorders involving these traits-specifically reading disability (RD) and language impairment (LI)-are common and prevent affected individuals from developing adequate communication skills, leaving them at risk for adverse academic, socioeconomic and psychiatric outcomes. Both RD and LI are complex traits that frequently co-occur, leading us to hypothesize that these disorders share genetic etiologies. To test this, we performed a genome-wide association study on individuals affected with both RD and LI in the Avon Longitudinal Study of Parents and Children. The strongest associations were seen with markers in ZNF385D (OR = 1.81, P = 5.45 × 10(-7) ) and COL4A2 (OR = 1.71, P = 7.59 × 10(-7) ). Markers within NDST4 showed the strongest associations with LI individually (OR = 1.827, P = 1.40 × 10(-7) ). We replicated association of ZNF385D using receptive vocabulary measures in the Pediatric Imaging Neurocognitive Genetics study (P = 0.00245). We then used diffusion tensor imaging fiber tract volume data on 16 fiber tracts to examine the implications of replicated markers. ZNF385D was a predictor of overall fiber tract volumes in both hemispheres, as well as global brain volume. Here, we present evidence for ZNF385D as a candidate gene for RD and LI. The implication of transcription factor ZNF385D in RD and LI underscores the importance of transcriptional regulation in the development of higher order neurocognitive traits. Further study is necessary to discern target genes of ZNF385D and how it functions within neural development of fluent language.

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.

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.

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.

Common genetic variation near the connexin-43 gene is associated with resting heart rate in African Americans: a genome-wide association study of 13,372 participants.

BACKGROUND: Genome-wide association studies have identified several genetic loci associated with variation in resting heart rate in European and Asian populations. No study has evaluated genetic variants associated with heart rate in African Americans. OBJECTIVE: To identify novel genetic variants associated with resting heart rate in African Americans. METHODS: Ten cohort studies participating in the Candidate-gene Association Resource and Continental Origins and Genetic Epidemiology Network consortia performed genome-wide genotyping of single nucleotide polymorphisms (SNPs) and imputed 2,954,965 SNPs using HapMap YRI and CEU panels in 13,372 participants of African ancestry. Each study measured the RR interval (ms) from 10-second resting 12-lead electrocardiograms and estimated RR-SNP associations using covariate-adjusted linear regression. Random-effects meta-analysis was used to combine cohort-specific measures of association and identify genome-wide significant loci (P≤2.5×10(-8)). RESULTS: Fourteen SNPs on chromosome 6q22 exceeded the genome-wide significance threshold. The most significant association was for rs9320841 (+13 ms per minor allele; P = 4.98×10(-15)). This SNP was approximately 350 kb downstream of GJA1, a locus previously identified as harboring SNPs associated with heart rate in Europeans. Adjustment for rs9320841 also attenuated the association between the remaining 13 SNPs in this region and heart rate. In addition, SNPs in MYH6, which have been identified in European genome-wide association study, were associated with similar changes in the resting heart rate as this population of African Americans. CONCLUSIONS: An intergenic region downstream of GJA1 (the gene encoding connexin 43, the major protein of the human myocardial gap junction) and an intragenic region within MYH6 are associated with variation in resting heart rate in African Americans as well as in populations of European and Asian origin.

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.

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.