Compositional structural brain signatures capture Alzheimer's genetic risk on brain structure along the disease continuum.

INTRODUCTION: Traditional brain imaging genetics studies have primarily focused on how genetic factors influence the volume of specific brain regions, often neglecting the overall complexity of brain architecture and its genetic underpinnings. METHODS: This study analyzed data from participants across the Alzheimer's disease (AD) continuum from the ALFA and ADNI studies. We exploited compositional data analysis to examine relative brain volumetric variations that (i) differentiate cognitively unimpaired (CU) individuals, defined as amyloid-negative (A-) based on CSF profiling, from those at different AD stages, and (ii) associated with increased genetic susceptibility to AD, assessed using polygenic risk scores. RESULTS: Distinct brain signatures differentiated CU A-individuals from amyloid-positive MCI and AD. Moreover, disease stage-specific signatures were associated with higher genetic risk of AD. DISCUSSION: The findings underscore the complex interplay between genetics and disease stages in shaping brain structure, which could inform targeted preventive strategies and interventions in preclinical AD.

Single-cell genomics and regulatory networks for 388 human brains.

Single-cell genomics is a powerful tool for studying heterogeneous tissues such as the brain. Yet, little is understood about how genetic variants influence cell-level gene expression. Addressing this, we uniformly processed single-nuclei, multi-omics datasets into a resource comprising >2.8M nuclei from the prefrontal cortex across 388 individuals. For 28 cell types, we assessed population-level variation in expression and chromatin across gene families and drug targets. We identified >550K cell-type-specific regulatory elements and >1.4M single-cell expression-quantitative-trait loci, which we used to build cell-type regulatory and cell-to-cell communication networks. These networks manifest cellular changes in aging and neuropsychiatric disorders. We further constructed an integrative model accurately imputing single-cell expression and simulating perturbations; the model prioritized ~250 disease-risk genes and drug targets with associated cell types.

An Insertion Within SIRPß1 Shows a Dual Effect Over Alzheimer's Disease Cognitive Decline Altering the Microglial Response.

Background: Microglial dysfunction plays a causative role in Alzheimer's disease (AD) pathogenesis. Here we focus on a germline insertion/deletion variant mapping SIRPß1, a surface receptor that triggers amyloid-ß(Aß) phagocytosis via TYROBP. Objective: To analyze the impact of this copy-number variant in SIRPß1 expression and how it affects AD molecular etiology. Methods: Copy-number variant proxy rs2209313 was evaluated in GERALD and GR@ACE longitudinal series. Hippocampal specimens of genotyped AD patients were also examined. SIRPß1 isoform-specific phagocytosis assays were performed in HEK393T cells. Results: The insertion alters the SIRPß1 protein isoform landscape compromising its ability to bind oligomeric Aß and its affinity for TYROBP. SIRPß1 Dup/Dup patients with mild cognitive impairment show an increased cerebrospinal fluid t-Tau/Aß ratio (p = 0.018) and a higher risk to develop AD (OR = 1.678, p = 0.018). MRIs showed that Dup/Dup patients exhibited a worse initial response to AD. At the moment of diagnosis, all patients showed equivalent Mini-Mental State Examination scores. However, AD patients with the duplication had less hippocampal degeneration (p < 0.001) and fewer white matter hyperintensities. In contrast, longitudinal studies indicate that patients bearing the duplication allele show a slower cognitive decline (p = 0.013). Transcriptional analysis also shows that the SIRPß1 duplication allele correlates with higher TREM2 expression and an increased microglial activation. Conclusions: The SIRPß1 internal duplication has opposite effects over MCI-to-Dementia conversion risk and AD progression, affecting microglial response to Aß. Given the pharmacological approaches focused on the TREM2-TYROBP axis, we believe that SIRPß1 structural variant might be considered as a potential modulator of this causative pathway.

A fast non-parametric test of association for multiple traits.

The increasing availability of multidimensional phenotypic data in large cohorts of genotyped individuals requires efficient methods to identify genetic effects on multiple traits. Permutational multivariate analysis of variance (PERMANOVA) offers a powerful non-parametric approach. However, it relies on permutations to assess significance, which hinders the analysis of large datasets. Here, we derive the limiting null distribution of the PERMANOVA test statistic, providing a framework for the fast computation of asymptotic p values. Our asymptotic test presents controlled type I error and high power, often outperforming parametric approaches. We illustrate its applicability in the context of QTL mapping and GWAS.

Multivariate Analysis and Modelling of multiple Brain endOphenotypes: Let's MAMBO!

Imaging genetic studies aim to test how genetic information influences brain structure and function by combining neuroimaging-based brain features and genetic data from the same individual. Most studies focus on individual correlation and association tests between genetic variants and a single measurement of the brain. Despite the great success of univariate approaches, given the capacity of neuroimaging methods to provide a multiplicity of cerebral phenotypes, the development and application of multivariate methods become crucial. In this article, we review novel methods and strategies focused on the analysis of multiple phenotypes and genetic data. We also discuss relevant aspects of multi-trait modelling in the context of neuroimaging data.

Identification and analysis of splicing quantitative trait loci across multiple tissues in the human genome.

Alternative splicing (AS) is a fundamental step in eukaryotic mRNA biogenesis. Here, we develop an efficient and reproducible pipeline for the discovery of genetic variants that affect AS (splicing QTLs, sQTLs). We use it to analyze the GTEx dataset, generating a comprehensive catalog of sQTLs in the human genome. Downstream analysis of this catalog provides insight into the mechanisms underlying splicing regulation. We report that a core set of sQTLs is shared across multiple tissues. sQTLs often target the global splicing pattern of genes, rather than individual splicing events. Many also affect the expression of the same or other genes, uncovering regulatory loci that act through different mechanisms. sQTLs tend to be located in post-transcriptionally spliced introns, which would function as hotspots for splicing regulation. While many variants affect splicing patterns by altering the sequence of splice sites, many more modify the binding sites of RNA-binding proteins. Genetic variants affecting splicing can have a stronger phenotypic impact than those affecting gene expression.

The impact of sex on gene expression across human tissues.

Many complex human phenotypes exhibit sex-differentiated characteristics. However, the molecular mechanisms underlying these differences remain largely unknown. We generated a catalog of sex differences in gene expression and in the genetic regulation of gene expression across 44 human tissue sources surveyed by the Genotype-Tissue Expression project (GTEx, v8 release). We demonstrate that sex influences gene expression levels and cellular composition of tissue samples across the human body. A total of 37% of all genes exhibit sex-biased expression in at least one tissue. We identify cis expression quantitative trait loci (eQTLs) with sex-differentiated effects and characterize their cellular origin. By integrating sex-biased eQTLs with genome-wide association study data, we identify 58 gene-trait associations that are driven by genetic regulation of gene expression in a single sex. These findings provide an extensive characterization of sex differences in the human transcriptome and its genetic regulation.

Cell type-specific genetic regulation of gene expression across human tissues.

The Genotype-Tissue Expression (GTEx) project has identified expression and splicing quantitative trait loci in cis (QTLs) for the majority of genes across a wide range of human tissues. However, the functional characterization of these QTLs has been limited by the heterogeneous cellular composition of GTEx tissue samples. We mapped interactions between computational estimates of cell type abundance and genotype to identify cell type-interaction QTLs for seven cell types and show that cell type-interaction expression QTLs (eQTLs) provide finer resolution to tissue specificity than bulk tissue cis-eQTLs. Analyses of genetic associations with 87 complex traits show a contribution from cell type-interaction QTLs and enables the discovery of hundreds of previously unidentified colocalized loci that are masked in bulk tissue.

A limited set of transcriptional programs define major cell types.

We have produced RNA sequencing data for 53 primary cells from different locations in the human body. The clustering of these primary cells reveals that most cells in the human body share a few broad transcriptional programs, which define five major cell types: epithelial, endothelial, mesenchymal, neural, and blood cells. These act as basic components of many tissues and organs. Based on gene expression, these cell types redefine the basic histological types by which tissues have been traditionally classified. We identified genes whose expression is specific to these cell types, and from these genes, we estimated the contribution of the major cell types to the composition of human tissues. We found this cellular composition to be a characteristic signature of tissues and to reflect tissue morphological heterogeneity and histology. We identified changes in cellular composition in different tissues associated with age and sex, and found that departures from the normal cellular composition correlate with histological phenotypes associated with disease.

Effect of BDNF Val66Met on hippocampal subfields volumes and compensatory interaction with APOE-ε4 in middle-age cognitively unimpaired individuals from the ALFA study.

BACKGROUND: Current evidence supports the involvement of brain-derived neurotrophic factor (BDNF) Val66Met polymorphism, and the ε4 allele of APOE gene in hippocampal-dependent functions. Previous studies on the association of Val66Met with whole hippocampal volume included patients of a variety of disorders. However, it remains to be elucidated whether there is an impact of BDNF Val66Met polymorphism on the volumes of the hippocampal subfield volumes (HSv) in cognitively unimpaired (CU) individuals, and the interactive effect with the APOE-ε4 status. METHODS: BDNF Val66Met and APOE genotypes were determined in a sample of 430 CU late/middle-aged participants from the ALFA study (ALzheimer and FAmilies). Participants underwent a brain 3D-T1-weighted MRI scan, and volumes of the HSv were determined using Freesurfer (v6.0). The effects of the BDNF Val66Met genotype on the HSv were assessed using general linear models corrected by age, gender, education, number of APOE-ε4 alleles and total intracranial volume. We also investigated whether the association between APOE-ε4 allele and HSv were modified by BDNF Val66Met genotypes. RESULTS: BDNF Val66Met carriers showed larger bilateral volumes of the subiculum subfield. In addition, HSv reductions associated with APOE-ε4 allele were significantly moderated by BDNF Val66Met status. BDNF Met carriers who were also APOE-ε4 homozygous showed patterns of higher HSv than BDNF Val carriers. CONCLUSION: To our knowledge, the present study is the first to show that carrying the BDNF Val66Met polymorphisms partially compensates the decreased on HSv associated with APOE-ε4 in middle-age cognitively unimpaired individuals.

ggsashimi: Sashimi plot revised for browser- and annotation-independent splicing visualization.

We present ggsashimi, a command-line tool for the visualization of splicing events across multiple samples. Given a specified genomic region, ggsashimi creates sashimi plots for individual RNA-seq experiments as well as aggregated plots for groups of experiments, a feature unique to this software. Compared to the existing versions of programs generating sashimi plots, it uses popular bioinformatics file formats, it is annotation-independent, and allows the visualization of splicing events even for large genomic regions by scaling down the genomic segments between splice sites. ggsashimi is freely available at https://github.com/guigolab/ggsashimi. It is implemented in python, and internally generates R code for plotting.

Effect of the sequence data deluge on the performance of methods for detecting protein functional residues.

BACKGROUND: The exponential accumulation of new sequences in public databases is expected to improve the performance of all the approaches for predicting protein structural and functional features. Nevertheless, this was never assessed or quantified for some widely used methodologies, such as those aimed at detecting functional sites and functional subfamilies in protein multiple sequence alignments. Using raw protein sequences as only input, these approaches can detect fully conserved positions, as well as those with a family-dependent conservation pattern. Both types of residues are routinely used as predictors of functional sites and, consequently, understanding how the sequence content of the databases affects them is relevant and timely. RESULTS: In this work we evaluate how the growth and change with time in the content of sequence databases affect five sequence-based approaches for detecting functional sites and subfamilies. We do that by recreating historical versions of the multiple sequence alignments that would have been obtained in the past based on the database contents at different time points, covering a period of 20 years. Applying the methods to these historical alignments allows quantifying the temporal variation in their performance. Our results show that the number of families to which these methods can be applied sharply increases with time, while their ability to detect potentially functional residues remains almost constant. CONCLUSIONS: These results are informative for the methods' developers and final users, and may have implications in the design of new sequencing initiatives.

Screening and Validation of Novel Biomarkers in Osteoarticular Pathologies by Comprehensive Combination of Protein Array Technologies.

Osteoarthritis (OA) is one of the most prevalent articular diseases. The identification of proteins closely associated with the diagnosis, progression, prognosis, and treatment response is dramatically required for this pathology. In this work, differential serum protein profiles have been identified in OA and rheumatoid arthritis (RA) by antibody arrays containing 151 antibodies against 121 antigens in a cohort of 36 samples. Then the identified differential serum protein profiles have been validated in a larger cohort of 282 samples. The overall immunoreactivity is higher in the pathological situations in comparison with the controls. Several proteins have been identified as biomarker candidates for OA and RA. Most of these biomarker candidates are proteins related to inflammatory response, lipid metabolism, or bone and extracellular matrix formation, degradation, or remodeling.

Genetic Drivers of Epigenetic and Transcriptional Variation in Human Immune Cells.

Characterizing the multifaceted contribution of genetic and epigenetic factors to disease phenotypes is a major challenge in human genetics and medicine. We carried out high-resolution genetic, epigenetic, and transcriptomic profiling in three major human immune cell types (CD14+ monocytes, CD16+ neutrophils, and naive CD4+ T cells) from up to 197 individuals. We assess, quantitatively, the relative contribution of cis-genetic and epigenetic factors to transcription and evaluate their impact as potential sources of confounding in epigenome-wide association studies. Further, we characterize highly coordinated genetic effects on gene expression, methylation, and histone variation through quantitative trait locus (QTL) mapping and allele-specific (AS) analyses. Finally, we demonstrate colocalization of molecular trait QTLs at 345 unique immune disease loci. This expansive, high-resolution atlas of multi-omics changes yields insights into cell-type-specific correlation between diverse genomic inputs, more generalizable correlations between these inputs, and defines molecular events that may underpin complex disease risk.