DiGAS: Differential gene allele spectrum as a descriptor in genetic studies.

Diagnosing individuals with complex genetic diseases is a challenging task. Computational methodologies exploit information at the genotype level by taking into account single nucleotide polymorphisms (SNPs) leveraging the results of genome-wide association studies analysis to assign a statistical significance to each SNP. Recent methodologies extend such an approach by aggregating SNP significance at the genetic level to identify genes that are related to the condition under study. However, such methodologies still suffer from the initial SNP analysis limitations. Here, we present DiGAS, a tool for diagnosing genetic conditions by computing significance, by means of SNP information, directly at the complex level of genetic regions. Such an approach is based on a generalized notion of allele spectrum, which evaluates the complete genetic alterations of the SNP set belonging to a genetic region at the population level. The statistical significance of a region is then evaluated through a differential allele spectrum analysis between the conditions of individuals belonging to the population. Tests, performed on well-established datasets regarding Alzheimer's disease, show that DiGAS outperforms the state of the art in distinguishing between sick and healthy subjects.

Traffic-related ultrafine particles impair mitochondrial functions in human olfactory mucosa cells - Implications for Alzheimer's disease.

Constituents of air pollution, the ultrafine particles (UFP) with a diameter of ≤0.1 µm, are considerably related to traffic emissions. Several studies link air pollution to Alzheimer's disease (AD), yet the exact relationship between the two remains poorly understood. Mitochondria are known targets of environmental toxicants, and their dysfunction is associated with neurodegenerative diseases. The olfactory mucosa (OM), located at the rooftop of the nasal cavity, is directly exposed to the environment and in contact with the brain. Mounting evidence suggests that the UFPs can impact the brain directly through the olfactory tract. By using primary human OM cultures established from nasal biopsies of cognitively healthy controls and individuals diagnosed with AD, we aimed to decipher the effects of traffic-related UFPs on mitochondria. The UFP samples were collected from the exhausts of a modern heavy-duty diesel engine (HDE) without aftertreatment systems, run with renewable diesel (A0) and petroleum diesel (A20), and from an engine of a 2019 model diesel passenger car (DI-E6d) equipped with state-of-the-art aftertreatment devices and run with renewable diesel (Euro6). OM cells were exposed to three different UFPs for 24-h and 72-h, after which cellular processes were assessed on the functional and transcriptomic levels. Our results show that UFPs impair mitochondrial functions in primary human OM cells by hampering oxidative phosphorylation (OXPHOS) and redox balance, and the responses of AD cells differ from cognitively healthy controls. RNA-Seq and IPA® revealed inhibition of OXPHOS and mitochondrial dysfunction in response to UFPs A0 and A20. Functional validation confirmed that A0 and A20 impair cellular respiration, decrease ATP levels, and disturb redox balance by altering NAD and glutathione metabolism, leading to increased ROS and oxidative stress. RNA-Seq and functional assessment revealed the presence of AD-related alterations in human OM cells and that different fuels and engine technologies elicit differential effects.

Exposure to source-specific air pollution in residential areas and its association with dementia incidence: a cohort study in Northern Sweden.

The aim of this study was to investigate the relationship between source-specific ambient particulate air pollution concentrations and the incidence of dementia. The study encompassed 70,057 participants from the Västerbotten intervention program cohort in Northern Sweden with a median age of 40 years at baseline. High-resolution dispersion models were employed to estimate source-specific particulate matter (PM) concentrations, such as PM10 and PM2.5 from traffic, exhaust, and biomass (mainly wood) burning, at the residential addresses of each participant. Cox regression models, adjusted for potential confounding factors, were used for the assessment. Over 884,847 person-years of follow-up, 409 incident dementia cases, identified through national registers, were observed. The study population's average exposure to annual mean total PM10 and PM2.5 lag 1-5 years was 9.50 µg/m3 and 5.61 µg/m3, respectively. Increased risks were identified for PM10-Traffic (35% [95% CI 0-82%]) and PM2.5-Exhaust (33% [95% CI - 2 to 79%]) in the second exposure tertile for lag 1-5 years, although no such risks were observed in the third tertile. Interestingly, a negative association was observed between PM2.5-Wood burning and the risk of dementia. In summary, this register-based study did not conclusively establish a strong association between air pollution exposure and the incidence of dementia. While some evidence indicated elevated risks for PM10-Traffic and PM2.5-Exhaust, and conversely, a negative association for PM2.5-Wood burning, no clear exposure-response relationships were evident.

GateMeClass: Gate Mining and Classification of cytometry data.

MOTIVATION: Cytometry comprises powerful techniques for analyzing the cell heterogeneity of a biological sample by examining the expression of protein markers. These technologies impact especially the field of oncoimmunology, where cell identification is essential to analyze the tumor microenvironment. Several classification tools have been developed for the annotation of cytometry datasets, which include supervised tools that require a training set as a reference (i.e. reference-based) and semisupervised tools based on the manual definition of a marker table. The latter is closer to the traditional annotation of cytometry data based on manual gating. However, they require the manual definition of a marker table that cannot be extracted automatically in a reference-based fashion. Therefore, we are lacking methods that allow both classification approaches while maintaining the high biological interpretability given by the marker table. RESULTS: We present a new tool called GateMeClass (Gate Mining and Classification) which overcomes the limitation of the current methods of classification of cytometry data allowing both semisupervised and supervised annotation based on a marker table that can be defined manually or extracted from an external annotated dataset. We measured the accuracy of GateMeClass for annotating three well-established benchmark mass cytometry datasets and one flow cytometry dataset. The performance of GateMeClass is comparable to reference-based methods and marker table-based techniques, offering greater flexibility and rapid execution times. AVAILABILITY AND IMPLEMENTATION: GateMeClass is implemented in R language and is publicly available at https://github.com/simo1c/GateMeClass.

Expression of the membrane tetraspanin claudin 18 on cancer cells promotes T lymphocyte infiltration and antitumor immunity in pancreatic cancer.

Tumors weakly infiltrated by T lymphocytes poorly respond to immunotherapy. We aimed to unveil malignancy-associated programs regulating T cell entrance, arrest, and activation in the tumor environment. Differential expression of cell adhesion and tissue architecture programs, particularly the presence of the membrane tetraspanin claudin (CLDN)18 as a signature gene, demarcated immune-infiltrated from immune-depleted mouse pancreatic tumors. In human pancreatic ductal adenocarcinoma (PDAC) and non-small cell lung cancer, CLDN18 expression positively correlated with more differentiated histology and favorable prognosis. CLDN18 on the cell surface promoted accrual of cytotoxic T lymphocytes (CTLs), facilitating direct CTL contacts with tumor cells by driving the mobilization of the adhesion protein ALCAM to the lipid rafts of the tumor cell membrane through actin. This process favored the formation of robust immunological synapses (ISs) between CTLs and CLDN18-positive cancer cells, resulting in increased T cell activation. Our data reveal an immune role for CLDN18 in orchestrating T cell infiltration and shaping the tumor immune contexture.

B-cell receptor signaling activity identifies patients with mantle cell lymphoma at higher risk of progression.

Mantle cell lymphoma (MCL) is an incurable B-cell malignancy characterized by a high clinical variability. Therefore, there is a critical need to define parameters that identify high-risk patients for aggressive disease and therapy resistance. B-cell receptor (BCR) signaling is crucial for MCL initiation and progression and is a target for therapeutic intervention. We interrogated BCR signaling proteins (SYK, LCK, BTK, PLCγ2, p38, AKT, NF-κB p65, and STAT5) in 30 primary MCL samples using phospho-specific flow cytometry. Anti-IgM modulation induced heterogeneous BCR signaling responses among samples allowing the identification of two clusters with differential responses. The cluster with higher response was associated with shorter progression free survival (PFS) and overall survival (OS). Moreover, higher constitutive AKT activity was predictive of inferior response to the Bruton's tyrosine kinase inhibitor (BTKi) ibrutinib. Time-to-event analyses showed that MCL international prognostic index (MIPI) high-risk category and higher STAT5 response were predictors of shorter PFS and OS whilst MIPI high-risk category and high SYK response predicted shorter OS. In conclusion, we identified BCR signaling properties associated with poor clinical outcome and resistance to ibrutinib, thus highlighting the prognostic and predictive significance of BCR activity and advancing our understanding of signaling heterogeneity underlying clinical behavior of MCL.

Identifying the joint signature of brain atrophy and gene variant scores in Alzheimer's Disease.

The joint modeling of genetic data and brain imaging information allows for determining the pathophysiological pathways of neurodegenerative diseases such as Alzheimer's disease (AD). This task has typically been approached using mass-univariate methods that rely on a complete set of Single Nucleotide Polymorphisms (SNPs) to assess their association with selected image-derived phenotypes (IDPs). However, such methods are prone to multiple comparisons bias and, most importantly, fail to account for potential cross-feature interactions, resulting in insufficient detection of significant associations. Ways to overcome these limitations while reducing the number of traits aim at conveying genetic information at the gene level and capturing the integrated genetic effects of a set of genetic variants, rather than looking at each SNP individually. Their associations with brain IDPs are still largely unexplored in the current literature, though they can uncover new potential genetic determinants for brain modulations in the AD continuum. In this work, we explored an explainable multivariate model to analyze the genetic basis of the grey matter modulations, relying on the AD Neuroimaging Initiative (ADNI) phase 3 dataset. Cortical thicknesses and subcortical volumes derived from T1-weighted Magnetic Resonance were considered to describe the imaging phenotypes. At the same time the genetic counterpart was represented by gene variant scores extracted by the Sequence Kernel Association Test (SKAT) filtering model. Moreover, transcriptomic analysis was carried on to assess the expression of the resulting genes in the main brain structures as a form of validation. Results highlighted meaningful genotype-phenotype interactionsas defined by three latent components showing a significant difference in the projection scores between patients and controls. Among the significant associations, the model highlighted EPHX1 and BCAS1 gene variant scores involved in neurodegenerative and myelination processes, hence relevant for AD. In particular, the first was associated with decreased subcortical volumes and the second with decreasedtemporal lobe thickness. Noteworthy, BCAS1 is particularly expressed in the dentate gyrus. Overall, the proposed approach allowed capturing genotype-phenotype interactions in a restricted study cohort that was confirmed by transcriptomic analysis, offering insights into the underlying mechanisms of neurodegeneration in AD in line with previous findings and suggesting new potential disease biomarkers.

PanDelos-frags: A methodology for discovering pangenomic content of incomplete microbial assemblies.

Pangenomics was originally defined as the problem of comparing the composition of genes into gene families within a set of bacterial isolates belonging to the same species. The problem requires the calculation of sequence homology among such genes. When combined with metagenomics, namely for human microbiome composition analysis, gene-oriented pangenome detection becomes a promising method to decipher ecosystem functions and population-level evolution. Established computational tools are able to investigate the genetic content of isolates for which a complete genomic sequence is available. However, there is a plethora of incomplete genomes that are available on public resources, which only a few tools may analyze. Incomplete means that the process for reconstructing their genomic sequence is not complete, and only fragments of their sequence are currently available. However, the information contained in these fragments may play an essential role in the analyses. Here, we present PanDelos-frags, a computational tool which exploits and extends previous results in analyzing complete genomes. It provides a new methodology for inferring missing genetic information and thus for managing incomplete genomes. PanDelos-frags outperforms state-of-the-art approaches in reconstructing gene families in synthetic benchmarks and in a real use case of metagenomics. PanDelos-frags is publicly available at https://github.com/InfOmics/PanDelos-frags.

Emissions from modern engines induce distinct effects in human olfactory mucosa cells, depending on fuel and aftertreatment.

Ultrafine particles (UFP) with a diameter of ≤0.1 µm, are contributors to ambient air pollution and derived mainly from traffic emissions, yet their health effects remain poorly characterized. The olfactory mucosa (OM) is located at the rooftop of the nasal cavity and directly exposed to both the environment and the brain. Mounting evidence suggests that pollutant particles affect the brain through the olfactory tract, however, the exact cellular mechanisms of how the OM responds to air pollutants remain poorly known. Here we show that the responses of primary human OM cells are altered upon exposure to UFPs and that different fuels and engines elicit different adverse effects. We used UFPs collected from exhausts of a heavy-duty-engine run with renewable diesel (A0) and fossil diesel (A20), and from a modern diesel vehicle run with renewable diesel (Euro6) and compared their health effects on the OM cells by assessing cellular processes on the functional and transcriptomic levels. Quantification revealed all samples as UFPs with the majority of particles being ≤0.1 µm by an aerodynamic diameter. Exposure to A0 and A20 induced substantial alterations in processes associated with inflammatory response, xenobiotic metabolism, olfactory signaling, and epithelial integrity. Euro6 caused only negligible changes, demonstrating the efficacy of aftertreatment devices. Furthermore, when compared to A20, A0 elicited less pronounced effects on OM cells, suggesting renewable diesel induces less adverse effects in OM cells. Prior studies and these results suggest that PAHs may disturb the inflammatory process and xenobiotic metabolism in the OM and that UFPs might mediate harmful effects on the brain through the olfactory route. This study provides important information on the adverse effects of UFPs in a human-based in vitro model, therefore providing new insight to form the basis for mitigation and preventive actions against the possible toxicological impairments caused by UFP exposure.

APDB: a database on air pollutant characterization and similarity prediction.

The World Health Organization estimates that 9 out of 10 people worldwide breathe air containing high levels of pollutants. Long-term and chronic exposure to high concentrations of air pollutants is associated with deleterious effects on vital organs, including increased inflammation in the lungs, oxidative stress in the heart and disruption of the blood-brain barrier. For this reason, in an effort to find an association between exposure to pollutants and the toxicological effects observable on human health, an online resource collecting and characterizing in detail pollutant molecules could be helpful to investigate their properties and mechanisms of action. We developed a database, APDB, collecting air-pollutant-related data from different online resources, in particular, molecules from the US Environmental Protection Agency, their associated targets and bioassays found in the PubChem chemical repository and their computed molecular descriptors and quantum mechanics properties. A web interface allows (i) to browse data by category, (ii) to navigate the database by querying molecules and targets and (iii) to visualize and download molecule and target structures as well as computed descriptors and similarities. The desired data can be freely exported in textual/tabular format and the whole database in SQL format. Database URL http://apdb.di.univr.it.

A survey on algorithms to characterize transcription factor binding sites.

Transcription factors (TFs) are key regulatory proteins that control the transcriptional rate of cells by binding short DNA sequences called transcription factor binding sites (TFBS) or motifs. Identifying and characterizing TFBS is fundamental to understanding the regulatory mechanisms governing the transcriptional state of cells. During the last decades, several experimental methods have been developed to recover DNA sequences containing TFBS. In parallel, computational methods have been proposed to discover and identify TFBS motifs based on these DNA sequences. This is one of the most widely investigated problems in bioinformatics and is referred to as the motif discovery problem. In this manuscript, we review classical and novel experimental and computational methods developed to discover and characterize TFBS motifs in DNA sequences, highlighting their advantages and drawbacks. We also discuss open challenges and future perspectives that could fill the remaining gaps in the field.

Esearch3D: propagating gene expression in chromatin networks to illuminate active enhancers.

Most cell type-specific genes are regulated by the interaction of enhancers with their promoters. The identification of enhancers is not trivial as enhancers are diverse in their characteristics and dynamic in their interaction partners. We present Esearch3D, a new method that exploits network theory approaches to identify active enhancers. Our work is based on the fact that enhancers act as a source of regulatory information to increase the rate of transcription of their target genes and that the flow of this information is mediated by the folding of chromatin in the three-dimensional (3D) nuclear space between the enhancer and the target gene promoter. Esearch3D reverse engineers this flow of information to calculate the likelihood of enhancer activity in intergenic regions by propagating the transcription levels of genes across 3D genome networks. Regions predicted to have high enhancer activity are shown to be enriched in annotations indicative of enhancer activity. These include: enhancer-associated histone marks, bidirectional CAGE-seq, STARR-seq, P300, RNA polymerase II and expression quantitative trait loci (eQTLs). Esearch3D leverages the relationship between chromatin architecture and transcription, allowing the prediction of active enhancers and an understanding of the complex underpinnings of regulatory networks. The method is available at: https://github.com/InfOmics/Esearch3D and https://doi.org/10.5281/zenodo.7737123.

Dynamic release of neuronal extracellular vesicles containing miR-21a-5p is induced by hypoxia.

Hypoxia induces changes in the secretion of extracellular vesicles (EVs) in several non-neuronal cells and pathological conditions. EVs are packed with biomolecules, such as microRNA(miR)-21-5p, which respond to hypoxia. However, the true EV association of miR-21-5p, and its functional or biomarker relevance, are inadequately characterised. Neurons are extremely sensitive cells, and it is not known whether the secretion of neuronal EVs and miR-21-5p are altered upon hypoxia. Here, we characterised the temporal EV secretion profile and cell viability of neurons under hypoxia. Hypoxia induced a rapid increase of miR-21a-5p secretion in the EVs, which preceded the elevation of hypoxia-induced tissue or cellular miR-21a-5p. Prolonged hypoxia induced cell death and the release of morphologically distinct EVs. The EVs protected miR-21a-5p from enzymatic degradation but a remarkable fraction of miR-21a-5p remained fragile and non-EV associated. The increase in miR-21a-5p secretion may have biomarker potential, as high blood levels of miR-21-5p in stroke patients were associated with significant disability at hospital discharge. Our data provides an understanding of the dynamic regulation of EV secretion from neurons under hypoxia and provides a candidate for the prediction of recovery from ischemic stroke.

Human genetic diversity alters off-target outcomes of therapeutic gene editing.

CRISPR gene editing holds great promise to modify DNA sequences in somatic cells to treat disease. However, standard computational and biochemical methods to predict off-target potential focus on reference genomes. We developed an efficient tool called CRISPRme that considers single-nucleotide polymorphism (SNP) and indel genetic variants to nominate and prioritize off-target sites. We tested the software with a BCL11A enhancer targeting guide RNA (gRNA) showing promise in clinical trials for sickle cell disease and ß-thalassemia and found that the top candidate off-target is produced by an allele common in African-ancestry populations (MAF 4.5%) that introduces a protospacer adjacent motif (PAM) sequence. We validated that SpCas9 generates strictly allele-specific indels and pericentric inversions in CD34+ hematopoietic stem and progenitor cells (HSPCs), although high-fidelity Cas9 mitigates this off-target. This report illustrates how genetic variants should be considered as modifiers of gene editing outcomes. We expect that variant-aware off-target assessment will become integral to therapeutic genome editing evaluation and provide a powerful approach for comprehensive off-target nomination.

Combined Large Cell Neuroendocrine Carcinomas of the Lung: Integrative Molecular Analysis Identifies Subtypes with Potential Therapeutic Implications.

BACKGROUND: Combined large cell neuroendocrine carcinoma (CoLCNEC) is given by the association of LCNEC with adeno or squamous or any non-neuroendocrine carcinoma. Molecular bases of CoLCNEC pathogenesis are scant and no standardized therapies are defined. METHODS: 44 CoLCNECs: 26 with adenocarcinoma (CoADC), 7 with squamous cell carcinoma (CoSQC), 3 with small cell carcinoma (CoSCLC), 4 with atypical carcinoid (CoAC) and 4 napsin-A positive LCNEC (NapA+), were assessed for alterations in 409 genes and transcriptomic profiling of 20,815 genes. RESULTS: Genes altered included TP53 (n = 30), RB1 (n = 14) and KRAS (n = 13). Targetable alterations included six KRAS G12C mutations and ALK-EML4 fusion gene. Comparison of CoLCNEC transcriptomes with 86 lung cancers of pure histology (8 AC, 19 ADC, 19 LCNEC, 11 SCLC and 29 SQC) identified CoLCNEC as a separate entity of neuroendocrine tumours with three different molecular profiles, two of which showed a non-neuroendocrine lineage. Hypomethylation, activation of MAPK signalling and association to immunotherapy signature specifically characterized each of three CoLCNEC molecular clusters. Prognostic stratification was also provided. CONCLUSIONS: CoLCNECs are an independent histologic category. Our findings support the extension of routine evaluation of KRAS mutations, fusion genes and immune-related markers to offer new perspectives in the therapeutic management of CoLCNEC.

Stardust: improving spatial transcriptomics data analysis through space-aware modularity optimization-based clustering.

BACKGROUND: Spatial transcriptomics (ST) combines stained tissue images with spatially resolved high-throughput RNA sequencing. The spatial transcriptomic analysis includes challenging tasks like clustering, where a partition among data points (spots) is defined by means of a similarity measure. Improving clustering results is a key factor as clustering affects subsequent downstream analysis. State-of-the-art approaches group data by taking into account transcriptional similarity and some by exploiting spatial information as well. However, it is not yet clear how much the spatial information combined with transcriptomics improves the clustering result. RESULTS: We propose a new clustering method, Stardust, that easily exploits the combination of space and transcriptomic information in the clustering procedure through a manual or fully automatic tuning of algorithm parameters. Moreover, a parameter-free version of the method is also provided where the spatial contribution depends dynamically on the expression distances distribution in the space. We evaluated the proposed methods results by analyzing ST data sets available on the 10x Genomics website and comparing clustering performances with state-of-the-art approaches by measuring the spots' stability in the clusters and their biological coherence. Stability is defined by the tendency of each point to remain clustered with the same neighbors when perturbations are applied. CONCLUSIONS: Stardust is an easy-to-use methodology allowing to define how much spatial information should influence clustering on different tissues and achieving more stable results than state-of-the-art approaches.