Leveraging European infrastructures to access 1 million human genomes by 2022.

Human genomics is undergoing a step change from being a predominantly research-driven activity to one driven through health care as many countries in Europe now have nascent precision medicine programmes. To maximize the value of the genomic data generated, these data will need to be shared between institutions and across countries. In recognition of this challenge, 21 European countries recently signed a declaration to transnationally share data on at least 1 million human genomes by 2022. In this Roadmap, we identify the challenges of data sharing across borders and demonstrate that European research infrastructures are well-positioned to support the rapid implementation of widespread genomic data access.

Author Correction: Leveraging European infrastructures to access 1 million human genomes by 2022.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Human herpesvirus diversity is altered in HLA class I binding peptides.

Herpesviruses are ubiquitous, genetically diverse DNA viruses, with long-term presence in humans associated with infrequent but significant pathology. Human leukocyte antigen (HLA) class I presents intracellularly derived peptide fragments from infected tissue cells to CD8+ T and natural killer cells, thereby directing antiviral immunity. Allotypes of highly polymorphic HLA class I are distinguished by their peptide binding repertoires. Because this HLA class I variation is a major determinant of herpesvirus disease, we examined if sequence diversity of virus proteins reflects evasion of HLA presentation. Using population genomic data from Epstein­Barr virus (EBV), human cytomegalovirus (HCMV), and Varicella­Zoster virus, we tested whether diversity differed between the regions of herpesvirus proteins that can be recognized, or not, by HLA class I. Herpesviruses exhibit lytic and latent infection stages, with the latter better enabling immune evasion. Whereas HLA binding peptides of lytic proteins are conserved, we found that EBV and HCMV proteins expressed during latency have increased peptide sequence diversity. Similarly, latent, but not lytic, herpesvirus proteins have greater population structure in HLA binding than nonbinding peptides. Finally, we found patterns consistent with EBV adaption to the local HLA environment, with less efficient recognition of EBV isolates by high-frequency HLA class I allotypes. Here, the frequency of CD8+ T cell epitopes inversely correlated with the frequency of HLA class I recognition. Previous analyses have shown that pathogen-mediated natural selection maintains exceptional polymorphism in HLA residues that determine peptide recognition. Here, we show that HLA class I peptide recognition impacts diversity of globally widespread pathogens.

A quality control portal for sequencing data deposited at the European genome-phenome archive.

Since its launch in 2008, the European Genome-Phenome Archive (EGA) has been leading the archiving and distribution of human identifiable genomic data. In this regard, one of the community concerns is the potential usability of the stored data, as of now, data submitters are not mandated to perform any quality control (QC) before uploading their data and associated metadata information. Here, we present a new File QC Portal developed at EGA, along with QC reports performed and created for 1 694 442 files [Fastq, sequence alignment map (SAM)/binary alignment map (BAM)/CRAM and variant call format (VCF)] submitted at EGA. QC reports allow anonymous EGA users to view summary-level information regarding the files within a specific dataset, such as quality of reads, alignment quality, number and type of variants and other features. Researchers benefit from being able to assess the quality of data prior to the data access decision and thereby, increasing the reusability of data (https://ega-archive.org/blog/data-upcycling-powered-by-ega/).

CAAStools: a toolbox to identify and test Convergent Amino Acid Substitutions.

MOTIVATION: Coincidence of Convergent Amino Acid Substitutions (CAAS) with phenotypic convergences allow pinpointing genes and even individual mutations that are likely to be associated with trait variation within their phylogenetic context. Such findings can provide useful insights into the genetic architecture of complex phenotypes. RESULTS: Here we introduce CAAStools, a set of bioinformatics tools to identify and validate CAAS in orthologous protein alignments for predefined groups of species representing the phenotypic values targeted by the user. AVAILABILITY AND IMPLEMENTATION: CAAStools source code is available at http://github.com/linudz/caastools, along with documentation and examples.

The European Genome-phenome Archive in 2021.

The European Genome-phenome Archive (EGA - https://ega-archive.org/) is a resource for long term secure archiving of all types of potentially identifiable genetic, phenotypic, and clinical data resulting from biomedical research projects. Its mission is to foster hosted data reuse, enable reproducibility, and accelerate biomedical and translational research in line with the FAIR principles. Launched in 2008, the EGA has grown quickly, currently archiving over 4,500 studies from nearly one thousand institutions. The EGA operates a distributed data access model in which requests are made to the data controller, not to the EGA, therefore, the submitter keeps control on who has access to the data and under which conditions. Given the size and value of data hosted, the EGA is constantly improving its value chain, that is, how the EGA can contribute to enhancing the value of human health data by facilitating its submission, discovery, access, and distribution, as well as leading the design and implementation of standards and methods necessary to deliver the value chain. The EGA has become a key GA4GH Driver Project, leading multiple development efforts and implementing new standards and tools, and has been appointed as an ELIXIR Core Data Resource.

Genetic characterization of the ALFA study: Uncovering genetic profiles in the Alzheimer's continuum.

INTRODUCTION: In 2013, the ALzheimer's and FAmilies (ALFA) project was established to investigate pathophysiological changes in preclinical Alzheimer's disease (AD), and to foster research on early detection and preventive interventions. METHODS: We conducted a comprehensive genetic characterization of ALFA participants with respect to neurodegenerative/cerebrovascular diseases, AD biomarkers, brain endophenotypes, risk factors and aging biomarkers. We placed particular emphasis on amyloid/tau status and assessed gender differences. Multiple polygenic risk scores were computed to capture different aspects of genetic predisposition. We additionally compared AD risk in ALFA to that across the full disease spectrum from the Alzheimer's Disease Neuroimaging Initiative (ADNI). RESULTS: Results show that the ALFA project has been successful at establishing a cohort of cognitively unimpaired individuals at high genetic predisposition of AD. DISCUSSION: It is, therefore, well-suited to study early pathophysiological changes in the preclinical AD continuum. Highlights Prevalence of ε4 carriers in ALzheimer and FAmilies (ALFA) is higher than in the general European population The ALFA study is highly enriched in Alzheimer's disease (AD) genetic risk factors beyond APOE AD genetic profiles in ALFA are similar to clinical groups along the continuum ALFA has succeeded in establishing a cohort of cognitively unimpaired individuals at high genetic AD risk ALFA is well suited to study pathogenic events/early pathophysiological changes in AD.

Copy number variants and fixed duplications among 198 rhesus macaques (Macaca mulatta).

The rhesus macaque is an abundant species of Old World monkeys and a valuable model organism for biomedical research due to its close phylogenetic relationship to humans. Copy number variation is one of the main sources of genomic diversity within and between species and a widely recognized cause of inter-individual differences in disease risk. However, copy number differences among rhesus macaques and between the human and macaque genomes, as well as the relevance of this diversity to research involving this nonhuman primate, remain understudied. Here we present a high-resolution map of sequence copy number for the rhesus macaque genome constructed from a dataset of 198 individuals. Our results show that about one-eighth of the rhesus macaque reference genome is composed of recently duplicated regions, either copy number variable regions or fixed duplications. Comparison with human genomic copy number maps based on previously published data shows that, despite overall similarities in the genome-wide distribution of these regions, there are specific differences at the chromosome level. Some of these create differences in the copy number profile between human disease genes and their rhesus macaque orthologs. Our results highlight the importance of addressing the number of copies of target genes in the design of experiments and cautions against human-centered assumptions in research conducted with model organisms. Overall, we present a genome-wide copy number map from a large sample of rhesus macaque individuals representing an important novel contribution concerning the evolution of copy number in primate genomes.

Beacon v2 and Beacon networks: A "lingua franca" for federated data discovery in biomedical genomics, and beyond.

Beacon is a basic data discovery protocol issued by the Global Alliance for Genomics and Health (GA4GH). The main goal addressed by version 1 of the Beacon protocol was to test the feasibility of broadly sharing human genomic data, through providing simple "yes" or "no" responses to queries about the presence of a given variant in datasets hosted by Beacon providers. The popularity of this concept has fostered the design of a version 2, that better serves real-world requirements and addresses the needs of clinical genomics research and healthcare, as assessed by several contributing projects and organizations. Particularly, rare disease genetics and cancer research will benefit from new case level and genomic variant level requests and the enabling of richer phenotype and clinical queries as well as support for fuzzy searches. Beacon is designed as a "lingua franca" to bridge data collections hosted in software solutions with different and rich interfaces. Beacon version 2 works alongside popular standards like Phenopackets, OMOP, or FHIR, allowing implementing consortia to return matches in beacon responses and provide a handover to their preferred data exchange format. The protocol is being explored by other research domains and is being tested in several international projects.

Comparative Analysis of Mammal Genomes Unveils Key Genomic Variability for Human Life Span.

The enormous mammal's lifespan variation is the result of each species' adaptations to their own biological trade-offs and ecological conditions. Comparative genomics have demonstrated that genomic factors underlying both, species lifespans and longevity of individuals, are in part shared across the tree of life. Here, we compared protein-coding regions across the mammalian phylogeny to detect individual amino acid (AA) changes shared by the most long-lived mammals and genes whose rates of protein evolution correlate with longevity. We discovered a total of 2,737 AA in 2,004 genes that distinguish long- and short-lived mammals, significantly more than expected by chance (P = 0.003). These genes belong to pathways involved in regulating lifespan, such as inflammatory response and hemostasis. Among them, a total 1,157 AA showed a significant association with maximum lifespan in a phylogenetic test. Interestingly, most of the detected AA positions do not vary in extant human populations (81.2%) or have allele frequencies below 1% (99.78%). Consequently, almost none of these putatively important variants could have been detected by genome-wide association studies. Additionally, we identified four more genes whose rate of protein evolution correlated with longevity in mammals. Crucially, SNPs located in the detected genes explain a larger fraction of human lifespan heritability than expected, successfully demonstrating for the first time that comparative genomics can be used to enhance interpretation of human genome-wide association studies. Finally, we show that the human longevity-associated proteins are significantly more stable than the orthologous proteins from short-lived mammals, strongly suggesting that general protein stability is linked to increased lifespan.

The genetic impact of an Ebola outbreak on a wild gorilla population.

BACKGROUND: Numerous Ebola virus outbreaks have occurred in Equatorial Africa over the past decades. Besides human fatalities, gorillas and chimpanzees have also succumbed to the fatal virus. The 2004 outbreak at the Odzala-Kokoua National Park (Republic of Congo) alone caused a severe decline in the resident western lowland gorilla (Gorilla gorilla gorilla) population, with a 95% mortality rate. Here, we explore the immediate genetic impact of the Ebola outbreak in the western lowland gorilla population. RESULTS: Associations with survivorship were evaluated by utilizing DNA obtained from fecal samples from 16 gorilla individuals declared missing after the outbreak (non-survivors) and 15 individuals observed before and after the epidemic (survivors). We used a target enrichment approach to capture the sequences of 123 genes previously associated with immunology and Ebola virus resistance and additionally analyzed the gut microbiome which could influence the survival after an infection. Our results indicate no changes in the population genetic diversity before and after the Ebola outbreak, and no significant differences in microbial community composition between survivors and non-survivors. However, and despite the low power for an association analysis, we do detect six nominally significant missense mutations in four genes that might be candidate variants associated with an increased chance of survival. CONCLUSION: This study offers the first insight to the genetics of a wild great ape population before and after an Ebola outbreak using target capture experiments from fecal samples, and presents a list of candidate loci that may have facilitated their survival.

Replicability and Prediction: Lessons and Challenges from GWAS.

Since the publication of the Wellcome Trust Case Control Consortium (WTCCC) landmark study a decade ago, genome-wide association studies (GWAS) have led to the discovery of thousands of risk variants involved in disease etiology. This success story has two angles that are often overlooked. First, GWAS findings are highly replicable. This is an unprecedented phenomenon in complex trait genetics, and indeed in many areas of science, which in past decades have been plagued by false positives. At a time of increasing concerns about the lack of reproducibility, we examine the biological and methodological reasons that account for the replicability of GWAS and identify the challenges ahead. In contrast to the exemplary success of disease gene discovery, at present GWAS findings are not useful for predicting phenotypes. We close with an overview of the prospects for individualized prediction of disease risk and its foreseeable impact in clinical practice.

The shared genetic architecture of schizophrenia, bipolar disorder and lifespan.

Psychiatric disorders such as Schizophrenia (SCZ) and Bipolar Disorder (BD) represent an evolutionary paradox, as they exhibit strong negative effects on fitness, such as decreased fecundity and early mortality, yet they persist at a worldwide prevalence of approximately 1%. Molecular mechanisms affecting lifespan, which may be widely common among complex diseases with fitness effects, can be studied by the integrated analysis of data from genome-wide association studies (GWAS) of human longevity together with any disease of interest. Here, we report the first of such studies, focusing on the genetic overlap-pleiotropy-between two psychiatric disorders with shortened lifespan, SCZ and BD, and human parental lifespan (PLS) as a surrogate of life expectancy. Our results are twofold: first, we demonstrate extensive polygenic overlap between SCZ and PLS and to a lesser extent between BD and PLS. Second, we identified novel loci shared between PLS and SCZ (n = 39), and BD (n = 8). Whereas most of the identified SCZ (66%) and BD (62%) pleiotropic risk alleles were associated with reduced lifespan, we also detected some antagonistic protective alleles associated to shorter lifespans. In fact, top-associated SNPs with SCZ seems to explain longevity variance explained (LVE) better than many other life-threatening diseases, including Type 2 diabetes and most cancers, probably due to a high overlap with smoking-related pathways. Overall, our study provides evidence of a genetic burden driven through premature mortality among people with SCZ, which can have profound implications for understanding, and potentially treating, the mortality gap associated with this psychiatric disorder.

Genome-phenome explorer (GePhEx): a tool for the visualization and interpretation of phenotypic relationships supported by genetic evidence.

MOTIVATION: Association studies based on SNP arrays and Next Generation Sequencing technologies have enabled the discovery of thousands of genetic loci related to human diseases. Nevertheless, their biological interpretation is still elusive, and their medical applications limited. Recently, various tools have been developed to help bridging the gap between genomes and phenomes. To our knowledge, however none of these tools allows users to retrieve the phenotype-wide list of genetic variants that may be linked to a given disease or to visually explore the joint genetic architecture of different pathologies. RESULTS: We present the Genome-Phenome Explorer (GePhEx), a web-tool easing the visual exploration of phenotypic relationships supported by genetic evidences. GePhEx is primarily based on the thorough analysis of linkage disequilibrium between disease-associated variants and also considers relationships based on genes, pathways or drug-targets, leveraging on publicly available variant-disease associations to detect potential relationships between diseases. We demonstrate that GePhEx does retrieve well-known relationships as well as novel ones, and that, thus, it might help shedding light on the patho-physiological mechanisms underlying complex diseases. To this end, we investigate the potential relationship between schizophrenia and lung cancer, first detected using GePhEx and provide further evidence supporting a functional link between them. AVAILABILITY AND IMPLEMENTATION: GePhEx is available at: https://gephex.ega-archive.org/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Variation in predicted COVID-19 risk among lemurs and lorises.

The novel coronavirus SARS-CoV-2, which in humans leads to the disease COVID-19, has caused global disruption and more than 2 million fatalities since it first emerged in late 2019. As we write, infection rates are at their highest point globally and are rising extremely rapidly in some areas due to more infectious variants. The primary target of SARS-CoV-2 is the cellular receptor angiotensin-converting enzyme-2 (ACE2). Recent sequence analyses of the ACE2 gene predict that many nonhuman primates are also likely to be highly susceptible to infection. However, the anticipated risk is not equal across the Order. Furthermore, some taxonomic groups show high ACE2 amino acid conservation, while others exhibit high variability at this locus. As an example of the latter, analyses of strepsirrhine primate ACE2 sequences to date indicate large variation among lemurs and lorises compared to other primate clades despite low sampling effort. Here, we report ACE2 gene and protein sequences for 71 individual strepsirrhines, spanning 51 species and 19 genera. Our study reinforces previous results while finding additional variability in other strepsirrhine species, and suggests several clades of lemurs have high potential susceptibility to SARS-CoV-2 infection. Troublingly, some species, including the rare and endangered aye-aye (Daubentonia madagascariensis), as well as those in the genera Avahi and Propithecus, may be at high risk. Given that lemurs are endemic to Madagascar and among the primates at highest risk of extinction globally, further understanding of the potential threat of COVID-19 to their health should be a conservation priority. All feasible actions should be taken to limit their exposure to SARS-CoV-2.

Targeted resequencing reveals rare variants enrichment in multiple sclerosis susceptibility genes.

Although genome-wide association studies have identified a number of common variants associated with multiple sclerosis (MS) susceptibility, little is known about the relevance of rare variants. Here, we aimed to explore the role of rare variants in 14 MS risk genes (FCRL1, RGS1, TIMMDC1, HHEX, CXCR5, LTBR, TSFM, GALC, TRAF3, STAT3, TNFSF14, IFI30, CD40, and CYP24A1) by targeted resequencing in an Iberian population of 524 MS cases and 546 healthy controls. Four rare variants-enriched regions within CYP24A1, FCRL1, RGS1, and TRAF3 were identified as significantly associated with MS. Functional studies revealed significantly decreased regulator of G protein signaling 1 (RGS1) gene expression levels in peripheral blood mononuclear cells from MS patients with RGS1 rare variants compared to noncarriers, whereas no significant differences in gene expression were observed for CYP24A1, FCRL1, and TRAF3 between rare variants carriers and noncarriers. Immunophenotyping showed significant decrease in RGS1 expression in peripheral blood B lymphocytes from MS patients with RGS1 rare variants relative to noncarriers. Lastly, peripheral blood mononuclear cell from MS patients carrying RGS1 rare variants showed significantly lower induction of RGS1 gene expression by interferon-ß compared to MS patients lacking RGS1 variants. The presence of rare variants in RGS1 reinforce the ideas of high genetic heterogeneity and a role of rare variants in MS pathogenesis.

A pharmacogenetic study implicates NINJ2 in the response to Interferon-ß in multiple sclerosis.

BACKGROUND: Multiple sclerosis (MS) is a disease in which biomarker identification is fundamental to predict response to treatments and to deliver the optimal drug to patients. We previously found an association between rs7298096, a polymorphism upstream to the NINJ2 gene, and the 4-year response to interferon-ß (IFNß) treatment in MS patients. OBJECTIVES: To analyse the association between rs7298096 and time to first relapse (TTFR) during IFNß therapy in MS patients and to better investigate its functional role. METHODS: Survival analysis was applied in three MS cohorts from different countries (n = 1004). We also studied the role of the polymorphism on gene expression using GTEx portal and a luciferase assay. We interrogated GEO datasets to explore the relationship between NINJ2 expression, IFNß and TTFR. RESULTS: Rs7298096AA patients show a shorter TTFR than rs7298096G-carriers (Pmeta-analysis = 3 × 10-4, hazard ratio = 1.41). Moreover, rs7298096AA is associated with a higher NINJ2 expression in blood (p = 7.0 × 10-6), which was confirmed in vitro (p = 0.009). Finally, NINJ2 expression is downregulated by IFNß treatment and related to TTFR. CONCLUSIONS: Rs7298096 could influence MS disease activity during IFNß treatment by modulating NINJ2 expression in blood. The gene encodes for an adhesion molecule involved in inflammation and endothelial cells activation, supporting its role in MS.

Derived immune and ancestral pigmentation alleles in a 7,000-year-old Mesolithic European.

Ancient genomic sequences have started to reveal the origin and the demographic impact of farmers from the Neolithic period spreading into Europe. The adoption of farming, stock breeding and sedentary societies during the Neolithic may have resulted in adaptive changes in genes associated with immunity and diet. However, the limited data available from earlier hunter-gatherers preclude an understanding of the selective processes associated with this crucial transition to agriculture in recent human evolution. Here we sequence an approximately 7,000-year-old Mesolithic skeleton discovered at the La Braña-Arintero site in León, Spain, to retrieve a complete pre-agricultural European human genome. Analysis of this genome in the context of other ancient samples suggests the existence of a common ancient genomic signature across western and central Eurasia from the Upper Paleolithic to the Mesolithic. The La Braña individual carries ancestral alleles in several skin pigmentation genes, suggesting that the light skin of modern Europeans was not yet ubiquitous in Mesolithic times. Moreover, we provide evidence that a significant number of derived, putatively adaptive variants associated with pathogen resistance in modern Europeans were already present in this hunter-gatherer.

Biological Processes Modulating Longevity across Primates: A Phylogenetic Genome-Phenome Analysis.

Aging is a complex process affecting different species and individuals in different ways. Comparing genetic variation across species with their aging phenotypes will help understanding the molecular basis of aging and longevity. Although most studies on aging have so far focused on short-lived model organisms, recent comparisons of genomic, transcriptomic, and metabolomic data across lineages with different lifespans are unveiling molecular signatures associated with longevity. Here, we examine the relationship between genomic variation and maximum lifespan across primate species. We used two different approaches. First, we searched for parallel amino-acid mutations that co-occur with increases in longevity across the primate linage. Twenty-five such amino-acid variants were identified, several of which have been previously reported by studies with different experimental setups and in different model organisms. The genes harboring these mutations are mainly enriched in functional categories such as wound healing, blood coagulation, and cardiovascular disorders. We demonstrate that these pathways are highly enriched for pleiotropic effects, as predicted by the antagonistic pleiotropy theory of aging. A second approach was focused on changes in rates of protein evolution across the primate phylogeny. Using the phylogenetic generalized least squares, we show that some genes exhibit strong correlations between their evolutionary rates and longevity-associated traits. These include genes in the Sphingosine 1-phosphate pathway, PI3K signaling, and the Thrombin/protease-activated receptor pathway, among other cardiovascular processes. Together, these results shed light into human senescence patterns and underscore the power of comparative genomics to identify pathways related to aging and longevity.