A Breast Cancer Polygenic Risk Score Is Feasible for Risk Stratification in the Norwegian Population.

BACKGROUND: Statistical associations of numerous single nucleotide polymorphisms with breast cancer (BC) have been identified in genome-wide association studies (GWAS). Recent evidence suggests that a Polygenic Risk Score (PRS) can be a useful risk stratification instrument for a BC screening strategy, and a PRS test has been developed for clinical use. The performance of the PRS is yet unknown in the Norwegian population. AIM: To evaluate the performance of PRS models for BC in a Norwegian dataset. METHODS: We investigated a sample of 1053 BC cases and 7094 controls from different regions of Norway. PRS values were calculated using four PRS models, and their performance was evaluated by the area under the curve (AUC) and the odds ratio (OR). The effect of the PRS on the age of onset of BC was determined by a Cox regression model, and the lifetime absolute risk of developing BC was calculated using the iCare tool. RESULTS: The best performing PRS model included 3820 SNPs, which yielded an AUC = 0.625 and an OR = 1.567 per one standard deviation increase. The PRS values of the samples correlate with an increased risk of BC, with a hazard ratio of 1.494 per one standard deviation increase (95% confidence interval of 1.406-1.588). The individuals in the highest decile of the PRS have at least twice the risk of developing BC compared to the individuals with a median PRS. The results in this study with Norwegian samples are coherent with the findings in the study conducted using Estonian and UK Biobank samples. CONCLUSION: The previously validated PRS models have a similar observed accuracy in the Norwegian data as in the UK and Estonian populations. A PRS provides a meaningful association with the age of onset of BC and lifetime risk. Therefore, as suggested in Estonia, a PRS may also be integrated into the screening strategy for BC in Norway.

Estimating number of European eel (Anguilla anguilla) individuals using environmental DNA and haplotype count in small rivers.

Knowledge about population genetic data is important for effective conservation management. Genetic research traditionally requires sampling directly from the organism, for example tissue, which can be challenging, time-consuming, and harmful to the animal. Environmental DNA (eDNA) approaches offer a way to sample genetic material noninvasively. In attempts to estimate population size of aquatic species using eDNA, researchers have found positive correlations between biomass and eDNA concentrations, but the approach is debated because of variations in the production and degrading of DNA in water. Recently, a more accurate eDNA-approach has emerged, focusing on the genomic differences between individuals. In this study, we used eDNA from water samples to estimate the number of European eel (Anguilla anguilla) individuals by examining haplotypes in the mitochondrial D-loop region, both in a closed aquatic environment with 10 eels of known haplotypes and in three rivers. The results revealed that it was possible to find every eel haplotype in the eDNA sample collected from the closed environment. We also found 13 unique haplotypes in the eDNA samples from the three rivers, which probably represent 13 eel individuals. This means that it is possible to obtain genomic information from European eel eDNA in water; however, more research is needed to develop the approach into a possible future tool for population quantification.

Combining population genomics with demographic analyses highlights habitat patchiness and larval dispersal as determinants of connectivity in coastal fish species.

Gene flow shapes spatial genetic structure and the potential for local adaptation. Among marine animals with nonmigratory adults, the presence or absence of a pelagic larval stage is thought to be a key determinant in shaping gene flow and the genetic structure of populations. In addition, the spatial distribution of suitable habitats is expected to influence the distribution of biological populations and their connectivity patterns. We used whole genome sequencing to study demographic history and reduced representation (double-digest restriction associated DNA) sequencing data to analyse spatial genetic structure in broadnosed pipefish (Syngnathus typhle). Its main habitat is eelgrass beds, which are patchily distributed along the study area in southern Norway. Demographic connectivity among populations was inferred from long-term (~30-year) population counts that uncovered a rapid decline in spatial correlations in abundance with distance as short as ~2 km. These findings were contrasted with data for two other fish species that have a pelagic larval stage (corkwing wrasse, Symphodus melops; black goby, Gobius niger). For these latter species, we found wider spatial scales of connectivity and weaker genetic isolation-by-distance patterns, except where both species experienced a strong barrier to gene flow, seemingly due to lack of suitable habitat. Our findings verify expectations that a fragmented habitat and absence of a pelagic larval stage promote genetic structure, while presence of a pelagic larvae stage increases demographic connectivity and gene flow, except perhaps over extensive habitat gaps.

Stabilizing selection on Atlantic cod supergenes through a millennium of extensive exploitation.

Life on Earth has been characterized by recurring cycles of ecological stasis and disruption, relating biological eras to geological and climatic transitions through the history of our planet. Due to the increasing degree of ecological abruption caused by human influences many advocate that we now have entered the geological era of the Anthropocene, or "the age of man." Considering the ongoing mass extinction and ecosystem reshuffling observed worldwide, a better understanding of the drivers of ecological stasis will be a requisite for identifying routes of intervention and mitigation. Ecosystem stability may rely on one or a few keystone species, and the loss of such species could potentially have detrimental effects. The Atlantic cod (Gadus morhua) has historically been highly abundant and is considered a keystone species in ecosystems of the northern Atlantic Ocean. Collapses of cod stocks have been observed on both sides of the Atlantic and reported to have detrimental effects that include vast ecosystem reshuffling. By whole-genome resequencing we demonstrate that stabilizing selection maintains three extensive "supergenes" in Atlantic cod, linking these genes to species persistence and ecological stasis. Genomic inference of historic effective population sizes shows continued declines for cod in the North Sea-Skagerrak-Kattegat system through the past millennia, consistent with an early onset of the marine Anthropocene through industrialization and commercialization of fisheries throughout the medieval period.

The genetic structure of Norway.

The aim of the present study was to describe the genetic structure of the Norwegian population using genotypes from 6369 unrelated individuals with detailed information about places of residence. Using standard single marker- and haplotype-based approaches, we report evidence of two regions with distinctive patterns of genetic variation, one in the far northeast, and another in the south of Norway, as indicated by fixation indices, haplotype sharing, homozygosity, and effective population size. We detect and quantify a component of Uralic Sami ancestry that is enriched in the North. On a finer scale, we find that rates of migration have been affected by topography like mountain ridges. In the broader Scandinavian context, we detect elevated relatedness between the mid- and northern border areas towards Sweden. The main finding of this study is that despite Norway's long maritime history and as a former Danish territory, the region closest to mainland Europe in the south appears to have been an isolated region in Norway, highlighting the open sea as a barrier to gene flow into Norway.

Demographic history has shaped the strongly differentiated corkwing wrasse populations in Northern Europe.

Understanding the biological processes involved in genetic differentiation and divergence between populations within species is a pivotal aim in evolutionary biology. One particular phenomenon that requires clarification is the maintenance of genetic barriers despite the high potential for gene flow in the marine environment. Such patterns have been attributed to limited dispersal or local adaptation, and to a lesser extent to the demographic history of the species. The corkwing wrasse (Symphodus melops) is an example of a marine fish species where regions of particular strong divergence are observed. One such genetic break occurred at a surprisingly small spatial scale (FST ~0.1), over a short coastline (<60 km) in the North Sea-Skagerrak transition area in southwestern Norway. Here, we investigate the observed divergence and purported reproductive isolation using genome resequencing. Our results suggest that historical events during the post-glacial recolonization route can explain the present population structure of the corkwing wrasse in the northeast Atlantic. While the divergence across the break is strong, we detected ongoing gene flow between populations over the break suggesting recent contact or negative selection against hybrids. Moreover, we found few outlier loci and no clear genomic regions potentially being under selection. We concluded that neutral processes and random genetic drift e.g., due to founder events during colonization have shaped the population structure in this species in Northern Europe. Our findings underline the need to take into account the demographic process in studies of divergence processes.

A continuous genome assembly of the corkwing wrasse (Symphodus melops).

The wrasses (Labridae) are one of the most successful and species-rich families of the Perciformes order of teleost fish. Its members display great morphological diversity, and occupy distinct trophic levels in coastal waters and coral reefs. The cleaning behaviour displayed by some wrasses, such as corkwing wrasse (Symphodus melops), is of particular interest for the salmon aquaculture industry to combat and control sea lice infestation as an alternative to chemicals and pharmaceuticals. There are still few genome assemblies available within this fish family for comparative and functional studies, despite the rapid increase in genome resources generated during the past years. Here, we present a highly continuous genome assembly of the corkwing wrasse using PacBio SMRT sequencing (x28.8) followed by error correction with paired-end Illumina data (x132.9). The present genome assembly consists of 5040 contigs (N50 = 461,652 bp) and a total size of 614 Mbp, of which 8.5% of the genome sequence encode known repeated elements. The genome assembly covers 94.21% of highly conserved genes across ray-finned fish species. We find evidence for increased copy numbers specific for corkwing wrasse possibly highlighting diversification and adaptive processes in gene families including N-linked glycosylation (ST8SIA6) and stress response kinases (HIPK1). By comparative analyses, we discover that de novo repeats, often not properly investigated during genome annotation, encode hundreds of immune-related genes. This new genomic resource, together with the ballan wrasse (Labrus bergylta), will allow for in-depth comparative genomics as well as population genetic analyses for the understudied wrasses.

Genetic Markers of Human Evolution Are Enriched in Schizophrenia.

BACKGROUND: Why schizophrenia has accompanied humans throughout our history despite its negative effect on fitness remains an evolutionary enigma. It is proposed that schizophrenia is a by-product of the complex evolution of the human brain and a compromise for humans' language, creative thinking, and cognitive abilities. METHODS: We analyzed recent large genome-wide association studies of schizophrenia and a range of other human phenotypes (anthropometric measures, cardiovascular disease risk factors, immune-mediated diseases) using a statistical framework that draws on polygenic architecture and ancillary information on genetic variants. We used information from the evolutionary proxy measure called the Neanderthal selective sweep (NSS) score. RESULTS: Gene loci associated with schizophrenia are significantly (p = 7.30 × 10(-9)) more prevalent in genomic regions that are likely to have undergone recent positive selection in humans (i.e., with a low NSS score). Variants in brain-related genes with a low NSS score confer significantly higher susceptibility than variants in other brain-related genes. The enrichment is strongest for schizophrenia, but we cannot rule out enrichment for other phenotypes. The false discovery rate conditional on the evolutionary proxy points to 27 candidate schizophrenia susceptibility loci, 12 of which are associated with schizophrenia and other psychiatric disorders or linked to brain development. CONCLUSIONS: Our results suggest that there is a polygenic overlap between schizophrenia and NSS score, a marker of human evolution, which is in line with the hypothesis that the persistence of schizophrenia is related to the evolutionary process of becoming human.

Altered Brain Activation during Emotional Face Processing in Relation to Both Diagnosis and Polygenic Risk of Bipolar Disorder.

OBJECTIVES: Bipolar disorder (BD) is a highly heritable disorder with polygenic inheritance. Among the most consistent findings from functional magnetic imaging (fMRI) studies are limbic hyperactivation and dorsal hypoactivation. However, the relation between reported brain functional abnormalities and underlying genetic risk remains elusive. This is the first cross-sectional study applying a whole-brain explorative approach to investigate potential influence of BD case-control status and polygenic risk on brain activation. METHODS: A BD polygenic risk score (PGRS) was estimated from the Psychiatric Genomics Consortium BD case-control study, and assigned to each individual in our independent sample (N=85 BD cases and 121 healthy controls (HC)), all of whom participated in an fMRI emotional faces matching paradigm. Potential differences in BOLD response across diagnostic groups were explored at whole-brain level in addition to amygdala as a region of interest. Putative effects of BD PGRS on brain activation were also investigated. RESULTS: At whole-brain level, BD cases presented with significantly lower cuneus/precuneus activation than HC during negative face processing (Z-threshold=2.3 as cluster-level correction). The PGRS was associated positively with increased right inferior frontal gyrus (rIFG) activation during negative face processing. For amygdala activation, there were no correlations with diagnostic status or PGRS. CONCLUSIONS: These findings are in line with previous reports of reduced precuneus and altered rIFG activation in BD. While these results demonstrate the ability of PGRS to reveal underlying genetic risk of altered brain activation in BD, the lack of convergence of effects at diagnostic and PGRS level suggests that this relation is a complex one.

MicroRNAs enrichment in GWAS of complex human phenotypes.

BACKGROUND: The genotype information carried by Genome-wide association studies (GWAS) seems to have the potential to explain more of the 'missing heritability' of complex human phenotypes, given improved statistical approaches. Several lines of evidence support the involvement of microRNA (miRNA) and other non-coding RNA in complex human traits and diseases. We employed a novel, genetic annotation-informed enrichment method for GWAS that captures more polygenic effects than standard GWAS analysis, to investigate if miRNA-tagging Single Nucleotide Polymorphisms (SNPs) are enriched of associations with 15 complex human phenotypes. We then leveraged the enrichment using a conditional False Discovery Rate (condFDR) approach to assess any improvement in the detection of individual miRNA SNPs associated with the disorders. RESULTS: We found SNPs tagging miRNA transcription regions to be significantly enriched of associations with 10 of 15 phenotypes. The enrichment remained significant after controlling for affiliation to other genomic categories, and was confirmed by replication. Albeit only nominally significant, enrichment was found also in miRNA binding sites for 10 phenotypes out of 15. Leveraging the enrichment in the condFDR framework, we observed a 2-4-fold increase in discovery of SNPs tagging miRNA regions. CONCLUSIONS: Our results suggest that miRNAs play an important role in the polygenic architecture of complex human disorders and traits, and therefore that miRNAs are a genomic category that can and should be used to improve gene discovery.

Abundant genetic overlap between blood lipids and immune-mediated diseases indicates shared molecular genetic mechanisms.

Epidemiological studies suggest a relationship between blood lipids and immune-mediated diseases, but the nature of these associations is not well understood. We used genome-wide association studies (GWAS) to investigate shared single nucleotide polymorphisms (SNPs) between blood lipids and immune-mediated diseases. We analyzed data from GWAS (n~200,000 individuals), applying new False Discovery Rate (FDR) methods, to investigate genetic overlap between blood lipid levels [triglycerides (TG), low density lipoproteins (LDL), high density lipoproteins (HDL)] and a selection of archetypal immune-mediated diseases (Crohn's disease, ulcerative colitis, rheumatoid arthritis, type 1 diabetes, celiac disease, psoriasis and sarcoidosis). We found significant polygenic pleiotropy between the blood lipids and all the investigated immune-mediated diseases. We discovered several shared risk loci between the immune-mediated diseases and TG (n = 88), LDL (n = 87) and HDL (n = 52). Three-way analyses differentiated the pattern of pleiotropy among the immune-mediated diseases. The new pleiotropic loci increased the number of functional gene network nodes representing blood lipid loci by 40%. Pathway analyses implicated several novel shared mechanisms for immune pathogenesis and lipid biology, including glycosphingolipid synthesis (e.g. FUT2) and intestinal host-microbe interactions (e.g. ATG16L1). We demonstrate a shared genetic basis for blood lipids and immune-mediated diseases independent of environmental factors. Our findings provide novel mechanistic insights into dyslipidemia and immune-mediated diseases and may have implications for therapeutic trials involving lipid-lowering and anti-inflammatory agents.

Polygenic risk for schizophrenia associated with working memory-related prefrontal brain activation in patients with schizophrenia and healthy controls.

Schizophrenia is a highly heritable and polygenic disease, and identified common genetic variants have shown weak individual effects. Many studies have reported altered working memory (WM)-related brain activation in schizophrenia, preferentially in the frontal lobe. Such differences in brain activations could reflect inherited alterations possibly involved in the disease etiology, or rather secondary disease-related mechanisms. The use of polygenic risk scores (PGRS) based on a large number of risk polymorphisms with small effects is a valuable approach to examine the effect of cumulative genetic risk on brain functioning. This study examined the impact of cumulative genetic risk for schizophrenia on WM-related brain activations, assessed with functional magnetic resonance imaging. For each participant (63 schizophrenia patients and 118 healthy controls), we calculated a PGRS for schizophrenia based on 18 862 single-nucleotide polymorphism in a large multicenter genome-wide association study including 9146 schizophrenia patients and 12 111 controls, performed by the Psychiatric Genomics Consortium. As expected, the PGRS was significantly higher in patients compared with healthy controls. Further, the PGRS was related to differences in frontal lobe brain activation between high and low WM demand. Specifically, even in absence of main effects of diagnosis, increased PGRS was associated with decreased activation difference in the right middle-superior prefrontal cortex (BA 10/11) and the right inferior frontal gyrus (BA 45). This effect was seen in both cases and controls, and was not influenced by sex, age, or task performance. The findings support the notion of dysregulation of frontal lobe functioning as an inherited vulnerability factor in schizophrenia.

Improved detection of common variants associated with schizophrenia and bipolar disorder using pleiotropy-informed conditional false discovery rate.

Several lines of evidence suggest that genome-wide association studies (GWAS) have the potential to explain more of the "missing heritability" of common complex phenotypes. However, reliable methods to identify a larger proportion of single nucleotide polymorphisms (SNPs) that impact disease risk are currently lacking. Here, we use a genetic pleiotropy-informed conditional false discovery rate (FDR) method on GWAS summary statistics data to identify new loci associated with schizophrenia (SCZ) and bipolar disorders (BD), two highly heritable disorders with significant missing heritability. Epidemiological and clinical evidence suggest similar disease characteristics and overlapping genes between SCZ and BD. Here, we computed conditional Q-Q curves of data from the Psychiatric Genome Consortium (SCZ; n = 9,379 cases and n = 7,736 controls; BD: n = 6,990 cases and n = 4,820 controls) to show enrichment of SNPs associated with SCZ as a function of association with BD and vice versa with a corresponding reduction in FDR. Applying the conditional FDR method, we identified 58 loci associated with SCZ and 35 loci associated with BD below the conditional FDR level of 0.05. Of these, 14 loci were associated with both SCZ and BD (conjunction FDR). Together, these findings show the feasibility of genetic pleiotropy-informed methods to improve gene discovery in SCZ and BD and indicate overlapping genetic mechanisms between these two disorders.

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.

Up-regulation of NOTCH4 gene expression in bipolar disorder.

OBJECTIVE: Immunopathogenic mechanisms have been implicated in schizophrenia and bipolar disorder, and genome-wide association studies (GWAS) point to the major histocompatibility complex, a region that contains many immune-related genes. One of the strongest candidate risk genes for schizophrenia and bipolar disorder is the NOTCH4 gene within the major histocompatibility complex. The authors investigated the NOTCH4 gene expression in individuals with bipolar disorder and schizophrenia relative to healthy comparison subjects and identified putative expression quantitative trait loci in and around the NOTCH4 gene. METHOD: The authors measured and compared NOTCH4 mRNA in whole blood in 690 individuals (479 patients and 211 healthy comparison subjects) and adjusted for a range of confounders. The authors also genotyped 20 single-nucleotide polymorphisms (SNPs) and investigated possible associations between expression quantitative trait loci and NOTCH4 expression. RESULTS: The authors found a strong association between NOTCH4 expression and bipolar disorder after adjusting for a range of confounders and multiple testing. In addition, seven SNPs within the NOTCH4 gene region were associated with enhanced NOTCH4 mRNA levels. Three of these expression quantitative trait loci were independent (not in linkage disequilibrium). CONCLUSIONS: The results indicate that the association between NOTCH4 DNA markers and bipolar disorder is related to altered function at the mRNA level, supporting the notion that NOTCH4 pathways are involved in the pathophysiology of bipolar disorder.

Identification of common variants associated with human hippocampal and intracranial volumes.

Identifying genetic variants influencing human brain structures may reveal new biological mechanisms underlying cognition and neuropsychiatric illness. The volume of the hippocampus is a biomarker of incipient Alzheimer's disease and is reduced in schizophrenia, major depression and mesial temporal lobe epilepsy. Whereas many brain imaging phenotypes are highly heritable, identifying and replicating genetic influences has been difficult, as small effects and the high costs of magnetic resonance imaging (MRI) have led to underpowered studies. Here we report genome-wide association meta-analyses and replication for mean bilateral hippocampal, total brain and intracranial volumes from a large multinational consortium. The intergenic variant rs7294919 was associated with hippocampal volume (12q24.22; N = 21,151; P = 6.70 × 10(-16)) and the expression levels of the positional candidate gene TESC in brain tissue. Additionally, rs10784502, located within HMGA2, was associated with intracranial volume (12q14.3; N = 15,782; P = 1.12 × 10(-12)). We also identified a suggestive association with total brain volume at rs10494373 within DDR2 (1q23.3; N = 6,500; P = 5.81 × 10(-7)).

Genome-wide association study identifies genetic loci associated with body mass index and high density lipoprotein-cholesterol levels during psychopharmacological treatment - a cross-sectional naturalistic study.

Metabolic and cardiovascular side effects are serious clinical problems related to psychopharmacological treatment, but the underlying mechanisms are mostly unknown. We performed a genome-wide association study of metabolic and cardiovascular risk factors during pharmacological therapy. Twelve indicators of metabolic side effects as well as cardiovascular risk factors were analyzed in a naturalistic sample of 594 patients of Norwegian ancestry. We analyzed interactions between gene variants and three categories of psychopharmacological agents based on their reported potential for side effects. For body mass index (BMI), two significantly associated loci were identified on 8q21.3. There were seven markers in one 30-kb region, and the strongest signal was rs7838490. In another locus 140kb away, six markers were significant, and rs6989402 obtained the strongest signal. Both of these loci are located upstream of the gene matrix metalloproteinase 16 (MMP16). For high density lipoprotein cholesterol (HDL-C), marker rs11615274 on 12q21 was significant. The results highlight three genomic regions potentially harboring susceptibility genes for drug-induced metabolic side effects, identifying MMP16 as a candidate gene. This deserves to be replicated in additional populations to provide more evidence for molecular genetic mechanisms of side effects during psychopharmacological treatment.

Gene-based analysis of regionally enriched cortical genes in GWAS data sets of cognitive traits and psychiatric disorders.

BACKGROUND: Despite its estimated high heritability, the genetic architecture leading to differences in cognitive performance remains poorly understood. Different cortical regions play important roles in normal cognitive functioning and impairment. Recently, we reported on sets of regionally enriched genes in three different cortical areas (frontomedial, temporal and occipital cortices) of the adult rat brain. It has been suggested that genes preferentially, or specifically, expressed in one region or organ reflect functional specialisation. Employing a gene-based approach to the analysis, we used the regionally enriched cortical genes to mine a genome-wide association study (GWAS) of the Norwegian Cognitive NeuroGenetics (NCNG) sample of healthy adults for association to nine psychometric tests measures. In addition, we explored GWAS data sets for the serious psychiatric disorders schizophrenia (SCZ) (n = 3 samples) and bipolar affective disorder (BP) (n = 3 samples), to which cognitive impairment is linked. PRINCIPAL FINDINGS: At the single gene level, the temporal cortex enriched gene RAR-related orphan receptor B (RORB) showed the strongest overall association, namely to a test of verbal intelligence (Vocabulary, P = 7.7E-04). We also applied gene set enrichment analysis (GSEA) to test the candidate genes, as gene sets, for enrichment of association signal in the NCNG GWAS and in GWASs of BP and of SCZ. We found that genes differentially expressed in the temporal cortex showed a significant enrichment of association signal in a test measure of non-verbal intelligence (Reasoning) in the NCNG sample. CONCLUSION: Our gene-based approach suggests that RORB could be involved in verbal intelligence differences, while the genes enriched in the temporal cortex might be important to intellectual functions as measured by a test of reasoning in the healthy population. These findings warrant further replication in independent samples on cognitive traits.