COSGAP: COntainerized Statistical Genetics Analysis Pipelines.

Summary: The collection and analysis of sensitive data in large-scale consortia for statistical genetics is hampered by multiple challenges, due to their non-shareable nature. Time-consuming issues in installing software frequently arise due to different operating systems, software dependencies, and limited internet access. For federated analysis across sites, it can be challenging to resolve different problems, including format requirements, data wrangling, setting up analysis on high-performance computing (HPC) facilities, etc. Easier, more standardized, automated protocols and pipelines can be solutions to overcome these issues. We have developed one such solution for statistical genetic data analysis using software container technologies. This solution, named COSGAP: "COntainerized Statistical Genetics Analysis Pipelines," consists of already established software tools placed into Singularity containers, alongside corresponding code and instructions on how to perform statistical genetic analyses, such as genome-wide association studies, polygenic scoring, LD score regression, Gaussian Mixture Models, and gene-set analysis. Using provided helper scripts written in Python, users can obtain auto-generated scripts to conduct the desired analysis either on HPC facilities or on a personal computer. COSGAP is actively being applied by users from different countries and projects to conduct genetic data analyses without spending much effort on software installation, converting data formats, and other technical requirements. Availability and implementation: COSGAP is freely available on GitHub (https://github.com/comorment/containers) under the GPLv3 license.

Pathway analysis of genetic markers associated with a functional MRI faces paradigm implicates polymorphisms in calcium responsive pathways.

Several lines of evidence suggest that common polygenic variation influences brain function in humans. Combining high-density genetic markers with brain imaging techniques is constricted by the practicalities of collecting sufficiently large brain imaging samples. Pathway analysis promises to leverage knowledge on function of genes to detect recurring signals of moderate effect. We adapt this approach, exploiting the deep information collected on brain function by fMRI methods, to identify molecular pathways containing genetic variants which influence brain activation during a commonly applied experiment based on a face matching task (n=246) which was developed to study neural processing of faces displaying negative emotions. Genetic markers moderately associated (p<10(-4)) with whole brain activation phenotypes constructed by applying principal components to contrast maps, were tested for pathway enrichment using permutation based methods. The most significant pathways are related to post NMDA receptor activation events, driven by genetic variants in calcium/calmodulin-dependent protein kinase II (CAMK2G, CAMK2D) and a calcium-regulated nucleotide exchange factor (RASGRF2) in which all are activated by intracellular calcium/calmodulin. The most significant effect of the combined polygenic model were localized to the left inferior frontal gyrus (p=1.03 × 10(-9)), a region primarily involved in semantic processing but also involved in processing negative emotions. These findings suggest that pathway analysis of GWAS results derived from principal component analysis of fMRI data is a promising method, to our knowledge, not previously described.

Deficit of homozygosity among 1.52 million individuals and genetic causes of recessive lethality.

Genotypes causing pregnancy loss and perinatal mortality are depleted among living individuals and are therefore difficult to find. To explore genetic causes of recessive lethality, we searched for sequence variants with deficit of homozygosity among 1.52 million individuals from six European populations. In this study, we identified 25 genes harboring protein-altering sequence variants with a strong deficit of homozygosity (10% or less of predicted homozygotes). Sequence variants in 12 of the genes cause Mendelian disease under a recessive mode of inheritance, two under a dominant mode, but variants in the remaining 11 have not been reported to cause disease. Sequence variants with a strong deficit of homozygosity are over-represented among genes essential for growth of human cell lines and genes orthologous to mouse genes known to affect viability. The function of these genes gives insight into the genetics of intrauterine lethality. We also identified 1077 genes with homozygous predicted loss-of-function genotypes not previously described, bringing the total set of genes completely knocked out in humans to 4785.

A family study of symbolic learning and synaptic plasticity in autism spectrum disorder.

The current study presents a male with autism spectrum disorder (ASD) and a 3q29 deletion, and three healthy first-degree relatives. Our magnetic resonance imaging (MRI) dataset included a healthy control subset. We describe a comprehensive multimodal approach, including equivalence class formation, neurocognitive testing, MRI, and electroencephalography (EEG)-based cortical plasticity, which can provide new insights into socio-communicative and learning impairments and neural underpinnings in ASD. On neurocognitive testing, the proband showed reduced processing speed, attending behavior, and executive function. He required more training trials in equivalence class training compared with family members and exhibited impaired priming of words compared with priming with images. The proband had smaller intracranial volume and surface area and a larger visual evoked potential (VEP) C1 amplitude than family members and intact long-term potentiation (LTP)-like visual cortex plasticity. Together, these results suggest that 3q29 deletion-related ASD is associated with impaired problem-solving strategies in complex socio-communicative and learning tasks, smaller intracranial and surface area, altered VEP amplitude, and normal LTP-like visual cortex plasticity. Further studies are needed to clarify whether this multimodal approach can be used to identify ASD subgroups with distinct neurobiological alterations and to uncover mechanisms underlying socio-communicative and learning impairments. Lay Summary: We studied learning, brain activity, and brain structure in a person with autism and a genetic aberration, and his close relatives. Compared with relatives, the person with autism required more training for learning, and visual learning was better than verbal learning. This person had some changes in the activity of the visual cortex, and the size and the surface area of the brain were reduced. Knowledge about learning and brain mechanisms is valuable for the development of training programs for individuals with autism.

A complex interplay of three R2R3 MYB transcription factors determines the profile of aliphatic glucosinolates in Arabidopsis.

While R2R3 MYB transcription factors are a large gene family of transcription factors within plants, comprehensive functional data in planta are still scarce. A model for studying R2R3 MYB control of metabolic networks is the glucosinolates (GLSs), secondary metabolites that control plant resistance against insects and pathogens and carry cancer-preventive properties. Three related members of the R2R3 MYB transcription factor family within Arabidopsis (Arabidopsis thaliana), MYB28, MYB29, and MYB76, are the commonly defined regulators of aliphatic GLS biosynthesis. We utilized new genotypes and systems analysis techniques to test the existing regulatory model in which MYB28 is the dominant regulator, MYB29 plays a minor rheostat role, and MYB76 is largely uninvolved. We unequivocally show that MYB76 is not dependent on MYB28 and MYB29 for induction of aliphatic GLSs and that MYB76 plays a role in determining the spatial distribution of aliphatic GLSs within the leaf, pointing at a potential role of MYB76 in transport regulation. Transcriptional profiling of knockout mutants revealed that GLS metabolite levels are uncoupled from the level of transcript accumulation for aliphatic GLS biosynthetic genes. This uncoupling of chemotypes from biosynthetic transcripts suggests revising our view of the regulation of GLS metabolism from a simple linear transcription factor-promoter model to a more modular system in which transcription factors cause similar chemotypes via nonoverlapping regulatory patterns. Similar regulatory networks might exist in other secondary pathways.

Attention-deficit hyperactivity disorder shares copy number variant risk with schizophrenia and autism spectrum disorder.

Attention-deficit/hyperactivity disorder (ADHD) is a highly heritable common childhood-onset neurodevelopmental disorder. Some rare copy number variations (CNVs) affect multiple neurodevelopmental disorders such as intellectual disability, autism spectrum disorders (ASD), schizophrenia and ADHD. The aim of this study is to determine to what extent ADHD shares high risk CNV alleles with schizophrenia and ASD. We compiled 19 neuropsychiatric CNVs and test 14, with sufficient power, for association with ADHD in Icelandic and Norwegian samples. Eight associate with ADHD; deletions at 2p16.3 (NRXN1), 15q11.2, 15q13.3 (BP4 & BP4.5-BP5) and 22q11.21, and duplications at 1q21.1 distal, 16p11.2 proximal, 16p13.11 and 22q11.21. Six of the CNVs have not been associated with ADHD before. As a group, the 19 CNVs associate with ADHD (OR = 2.43, P = 1.6 × 10-21), even when comorbid ASD and schizophrenia are excluded from the sample. These results highlight the pleiotropic effect of the neuropsychiatric CNVs and add evidence for ADHD, ASD and schizophrenia being related neurodevelopmental disorders rather than distinct entities.

Effects of autozygosity and schizophrenia polygenic risk on cognitive and brain developmental trajectories.

Cognitive and brain development are determined by dynamic interactions between genes and environment across the lifespan. Aside from marker-by-marker analyses of polymorphisms, biologically meaningful features of the whole genome (derived from the combined effect of individual markers) have been postulated to inform on human phenotypes including cognitive traits and their underlying biological substrate. Here, estimates of inbreeding and genetic susceptibility for schizophrenia calculated from genome-wide data-runs of homozygosity (ROH) and schizophrenia polygenic risk score (PGRS)-are analyzed in relation to cognitive abilities (n = 4183) and brain structure (n = 516) in a general-population sample of European-ancestry participants aged 8-22, from the Philadelphia Neurodevelopmental Cohort. The findings suggest that a higher ROH burden and higher schizophrenia PGRS are associated with higher intelligence. Cognition-ROH and cognition-PGRS associations obtained in this cohort may, respectively, evidence that assortative mating influences intelligence, and that individuals with high schizophrenia genetic risk who do not transition to disease status are cognitively resilient. Neuroanatomical data showed that the effects of schizophrenia PGRS on cognition could be modulated by brain structure, although larger imaging datasets are needed to accurately disentangle the underlying neural mechanisms linking IQ with both inbreeding and the genetic burden for schizophrenia.

Investigation of the genetic interaction between BDNF and DRD3 genes in suicidical behaviour in psychiatric disorders.

OBJECTIVES: Suicide is a serious public health concern, and it is partly genetic. The brain-derived neurotrophic factor (BDNF) gene has been a strong candidate in genetic studies of suicide (Dwivedi et al., Arch Gen Psychiatry 2010;60:804-815; Zai et al., Prog Neuropsychopharmacol Biol Psychiatry 2012;34:1412-1418) and BDNF regulates the expression of the dopamine D3 receptor. OBJECTIVE: We examined the role of the BDNF and DRD3 genes in suicide. METHODS: We analysed four tag single-nucleotide polymorphisms (SNPs) in BDNF and 15 SNPs in the D3 receptor gene DRD3 for possible association with suicide attempt history in our Canadian sample of Schizophrenia (SCZ) patients of European ancestry (N = 188). RESULTS: In this sample, we found a possible interaction between the BDNF Val66Met and DRD3 Ser9Gly SNPs in increasing the risk of suicide attempt(s) in our SCZ sample. Specifically, a larger proportion of SCZ patients who were carrying at least one copy of the minor allele at each of the Val66Met and Ser9Gly functional markers have attempted suicides compared to patients with other genotypes (Bonferroni P < 0.05). However, we could not replicate this finding in samples from other psychiatric populations. CONCLUSIONS: Taken together, the results from the present study suggest that an interaction between BDNF and DRD3 may not play a major role in the risk for suicide attempt, though further studies, especially in SCZ, are required.

Association analysis between suicidal behaviour and candidate genes of bipolar disorder and schizophrenia.

BACKGROUND: The present study investigated associations between the strongest joint genetic risk variants for bipolar disorder (BD) and schizophrenia (SCZ) and a history of suicide attempt in patients with BD, SCZ and related psychiatric disorders. METHODS: A history of suicide attempt was assessed in a sample of 1009 patients with BD, SCZ and related psychosis spectrum disorders, and associations with the joint genetic risk variants for BD and SCZ (rs2239547 (ITIH3/4-region), rs10994359 (ANK3) and rs4765905 (CACNA1C)) were investigated. Previously reported susceptibility loci for suicide attempt in BD were also investigated. Associations were tested by logistic regression with Bonferroni correction for multiple testing. RESULTS: The risk allele in rs2239547 (ITIH3/4-region) was significantly associated with a history of suicide attempt (p=0.01) after multiple testing correction (p threshold<0.017). The previous suicide attempt susceptibility loci were only nominally associated, but had the same direction of risk in the replication sample (sign test, p=0.02). LIMITATIONS: Relatively small sample size and retrospective clinical assessment. CONCLUSIONS: We detected a novel association between suicide attempt and the ITIH3/4-region in a combined group of patients with BD, SCZ and related psychosis spectrum disorders. This may be useful in understanding molecular mechanisms of suicidal behaviour in severe mental disorders, although replication is warranted.

A systems biology approach identifies a R2R3 MYB gene subfamily with distinct and overlapping functions in regulation of aliphatic glucosinolates.

BACKGROUND: Glucosinolates are natural metabolites in the order Brassicales that defend plants against both herbivores and pathogens and can attract specialized insects. Knowledge about the genes controlling glucosinolate regulation is limited. Here, we identify three R2R3 MYB transcription factors regulating aliphatic glucosinolate biosynthesis in Arabidopsis by combining several systems biology tools. METHODOLOGY/PRINCIPAL FINDINGS: MYB28 was identified as a candidate regulator of aliphatic glucosinolates based on its co-localization within a genomic region controlling variation both in aliphatic glucosinolate content (metabolite QTL) and in transcript level for genes involved in the biosynthesis of aliphatic glucosinolates (expression QTL), as well as its co-expression with genes in aliphatic glucosinolate biosynthesis. A phylogenetic analysis with the R2R3 motif of MYB28 showed that it and two homologues, MYB29 and MYB76, were members of an Arabidopsis-specific clade that included three characterized regulators of indole glucosinolates. Over-expression of the individual MYB genes showed that they all had the capacity to increase the production of aliphatic glucosinolates in leaves and seeds and induce gene expression of aliphatic biosynthetic genes within leaves. Analysis of leaves and seeds of single knockout mutants showed that mutants of MYB29 and MYB76 have reductions in only short-chained aliphatic glucosinolates whereas a mutant in MYB28 has reductions in both short- and long-chained aliphatic glucosinolates. Furthermore, analysis of a double knockout in MYB28 and MYB29 identified an emergent property of the system since the absence of aliphatic glucosinolates in these plants could not be predicted by the chemotype of the single knockouts. CONCLUSIONS/SIGNIFICANCE: It seems that these cruciferous-specific MYB regulatory genes have evolved both overlapping and specific regulatory capacities. This provides a unique system within which to study the evolution of MYB regulatory factors and their downstream targets.