GENIX enables comparative network analysis of single-cell RNA sequencing to reveal signatures of therapeutic interventions.

Single-cell RNA sequencing (scRNA-seq) has transformed our understanding of cellular responses to perturbations such as therapeutic interventions and vaccines. Gene relevance to such perturbations is often assessed through differential expression analysis (DEA), which offers a one-dimensional view of the transcriptomic landscape. This method potentially overlooks genes with modest expression changes but profound downstream effects and is susceptible to false positives. We present GENIX (gene expression network importance examination), a computational framework that transcends DEA by constructing gene association networks and employing a network-based comparative model to identify topological signature genes. We benchmark GENIX using both synthetic and experimental datasets, including analysis of influenza vaccine-induced immune responses in peripheral blood mononuclear cells (PBMCs) from recovered COVID-19 patients. GENIX successfully emulates key characteristics of biological networks and reveals signature genes that are missed by classical DEA, thereby broadening the scope of target gene discovery in precision medicine.

Genetic analyses of inflammatory polyneuropathy and chronic inflammatory demyelinating polyradiculoneuropathy identified candidate genes.

Chronic inflammatory demyelinating polyneuropathy (CIDP) is a rare, immune-mediated disorder in which an aberrant immune response causes demyelination and axonal damage of the peripheral nerves. Genetic contribution to CIDP is unclear and no genome-wide association study (GWAS) has been reported so far. In this study, we aimed to identify CIDP-related risk loci, genes, and pathways. We first focused on CIDP, and 516 CIDP cases and 403,545 controls were included in the GWAS analysis. We also investigated genetic risk for inflammatory polyneuropathy (IP), in which we performed a GWAS study using FinnGen data and combined the results with GWAS from the UK Biobank using a fixed-effect meta-analysis. A total of 1,261 IP cases and 823,730 controls were included in the analysis. Stratified analyses by gender were performed. Mendelian randomization (MR), colocalization, and transcriptome-wide association study (TWAS) analyses were performed to identify associated genes. Gene-set analyses were conducted to identify associated pathways. We identified one genome-wide significant locus at 20q13.33 for CIDP risk among women, the top variant located at the intron region of gene CDH4. Sex-combined MR, colocalization, and TWAS analyses identified three candidate pathogenic genes for CIDP and five genes for IP. MAGMA gene-set analyses identified a total of 18 pathways related to IP or CIDP. Sex-stratified analyses identified three genes for IP among males and two genes for IP among females. Our study identified suggestive risk genes and pathways for CIDP and IP. Functional analyses should be conducted to further confirm these associations.

Novel biomarkers of preterm brain injury from blood transcriptome in sheep model of intrauterine asphyxia.

BACKGROUND: Infants born preterm have a higher incidence of neurological deficits. A key step in finding effective treatments is to identify biomarkers that reliably predict outcome. METHODS: Following umbilical cord occlusion (UCO) in pregnant sheep, whole fetal blood RNA was sequenced pre- and post-UCO, brain injury outcome was determined by battery of neuropathology scoring and the transcriptome signature correlated to the degree of brain injury. Additionally, we developed a novel analytical procedure to deduce cell blood composition over time. RESULTS: Sixty-one genes were identified with significant altered expression after UCO. In pre-UCO blood, the level of three mRNAs (Trex2, Znf280b, novel miRNA) and in post-UCO, four mRNAs (Fam184a, Angptl2, novel lincRNA and an unknown protein-coding gene) were associated to brain injury (FDR < 0.01). Several of these mRNAs are related to inflammation and angiogenesis. Pathway analysis highlighted genes playing a role in perinatal death and growth failure. Results also indicate that several leukocyte populations undergo significant changes after UCO. CONCLUSION: We have used a whole transcriptomic approach to uncover novel biomarkers in fetal blood that correlate to neuropathology in the preterm sheep brain. The current data forms a basis for future studies to investigate mechanisms of these mRNAs in the injury progression. IMPACT: Trend analysis of genes following asphyxia reveal a group of genes associated with perinatal death and growth failure. Several pre-asphyxia transcripts were associated to brain injury severity suggesting genomic susceptibility to injury. Several post-asphyxia transcripts were correlated to brain injury severity, thus, serve as potential novel biomarkers of injury outcome. Successfully adaptation of cell profiling algorithms suggests significant changes in blood cell composition following asphyxia.

Scaling up single-cell RNA-seq data analysis with CellBridge workflow.

SUMMARY: Single-cell RNA sequencing (scRNA-seq) has revolutionized the study of gene expression at the individual cell level, unraveling unprecedented insights into cellular heterogeneity. However, the analysis of scRNA-seq data remains a challenging and time-consuming task, often demanding advanced computational expertise, rendering it impractical for high-volume environments and applications. We present CellBridge, an automated workflow designed to simplify the standard procedures entailed in scRNA-seq data analysis, eliminating the need for specialized computational expertise. CellBridge utilizes state-of-the-art computational methods, integrating a range of advanced functionalities, covering the entire process from raw unaligned sequencing reads to cell type annotation. Hence, CellBridge accelerates the pace of discovery by seamlessly enabling insights into vast volumes of scRNA-seq data, without compromising workflow control and reproducibility. AVAILABILITY AND IMPLEMENTATION: The source code, detailed documentation, and materials required to reproduce the results are available on GitHub and archived in Zenodo. For the CellBridge pre-processing step (v1.0.0), access the GitHub repository at https://github.com/Sanofi-Public/PMCB-ToBridge and the Zenodo archive at https://zenodo.org/records/10246161. For the CellBridge processing step (v1.0.0), visit the GitHub repository at https://github.com/Sanofi-Public/PMCB-CellBridge and the Zenodo archive at https://zenodo.org/records/10246046.

Genome-wide association study analysis of disease severity in Acne reveals novel biological insights.

Acne vulgaris is a common skin disease that affects >85% of teenage young adults among which >8% develop severe lesions that leaves permanent scars. Genetic heritability studies of acne in twin cohorts have estimated that the heritability for acne is 80%. Previous genome-wide association studies (GWAS) have identified 50 genetic loci associated with increased risk of developing acne when compared to healthy individuals. However only a few studies have investigated genetic association with disease severity. GWAS of disease progression may provide a more effective approach to unveil potential disease modifying therapeutic targets. Here, we performed a multi-ethnic GWAS analysis to capture disease severity in acne patients by using individuals with normal acne as a control. Our cohort consists of a total of 2,956 participants, including 290 severe acne cases and 930 normal acne controls from FinnGen, and 522 cases and 1,214 controls from BioVU. We also performed mendelian randomization (MR), colocalization analyses and transcriptome-wide association study (TWAS) to identify putative causal genes. Lastly, we performed gene-set enrichment analysis using MAGMA to implicate biological pathways that drive disease severity in Acne. We identified two new loci associated with acne severity at the genome-wide significance level, six novel associated genes by MR, colocalization and TWAS analyses, including genes CDC7, SLC7A1, ADAM23, TTLL10, CDK20 and DNAJA4 , and 5 novel pathways by MAGMA analyses. Our study suggests that the etiologies of acne susceptibility and severity have limited overlap, with only 26% of known acne risk loci presenting nominal association with acne severity and none of the novel severity associated genes reported as associated with acne risk in previous GWAS.

IL-33 Expression Is Lower in Current Smokers at both Transcriptomic and Protein Levels.

Rationale: IL-33 is a proinflammatory cytokine thought to play a role in the pathogenesis of asthma and chronic obstructive pulmonary disease (COPD). A recent clinical trial using an anti-IL-33 antibody showed a reduction in exacerbation and improved lung function in ex-smokers but not current smokers with COPD. Objectives: This study aimed to understand the effects of smoking status on IL-33. Methods: We investigated the association of smoking status with the level of gene expression of IL-33 in the airways in eight independent transcriptomic studies of lung airways. Additionally, we performed Western blot analysis and immunohistochemistry for IL-33 in lung tissue to assess protein levels. Measurements and Main Results: Across the bulk RNA-sequencing datasets, IL-33 gene expression and its signaling pathway were significantly lower in current versus former or never-smokers and increased upon smoking cessation (P < 0.05). Single-cell sequencing showed that IL-33 is predominantly expressed in resting basal epithelial cells and decreases during the differentiation process triggered by smoke exposure. We also found a higher transitioning of this cellular subpopulation into a more differentiated cell type during chronic smoking, potentially driving the reduction of IL-33. Protein analysis demonstrated lower IL-33 levels in lung tissue from current versus former smokers with COPD and a lower proportion of IL-33-positive basal cells in current versus ex-smoking controls. Conclusions: We provide strong evidence that cigarette smoke leads to an overall reduction in IL-33 expression in transcriptomic and protein level, and this may be due to the decrease in resting basal cells. Together, these findings may explain the clinical observation that a recent antibody-based anti-IL-33 treatment is more effective in former than current smokers with COPD.

A marker gene-based method for identifying the cell-type of origin from single-cell RNA sequencing data.

Single-cell RNA sequencing (scRNA-seq) experiments provide opportunities to peer into complex tissues at single-cell resolution. However, insightful biological interpretation of scRNA-seq data relies upon precise identification of cell types. The ability to identify the origin of a cell quickly and accurately will greatly improve downstream analyses. We present Sargent, a transformation-free, cluster-free, single-cell annotation algorithm for rapidly identifying the cell types of origin based on cell type-specific markers. We demonstrate Sargent's high accuracy by annotating simulated datasets. Further, we compare Sargent performance against expert-annotated scRNA-seq data from human organs including PBMC, heart, kidney, and lung. We demonstrate that Sargent retains both the flexibility and biological interpretability of cluster-based manual annotation. Additionally, the automation eliminates the labor intensive and potentially biased user annotation, producing robust, reproducible, and scalable outputs.•Sargent is a transformation-free, cluster-free, single-cell annotation algorithm for rapidly identifying the cell types of origin based on cell type-specific markers.•Sargent retains both the flexibility and biological interpretability of cluster-based manual annotation.•Automation eliminates the labor intensive and potentially biased user annotation, producing robust, reproducible, and scalable outputs.

DiSiR: fast and robust method to identify ligand-receptor interactions at subunit level from single-cell RNA-sequencing data.

Most cell-cell interactions and crosstalks are mediated by ligand-receptor interactions. The advent of single-cell RNA-sequencing (scRNA-seq) techniques has enabled characterizing tissue heterogeneity at single-cell level. In the past few years, several methods have been developed to study ligand-receptor interactions at cell type level using scRNA-seq data. However, there is still no easy way to query the activity of a specific user-defined signaling pathway in a targeted way or to map the interactions of the same subunit with different ligands as part of different receptor complexes. Here, we present DiSiR, a fast and easy-to-use permutation-based software framework to investigate how individual cells are interacting with each other by analyzing signaling pathways of multi-subunit ligand-activated receptors from scRNA-seq data, not only for available curated databases of ligand-receptor interactions, but also for interactions that are not listed in these databases. We show that, when utilized to infer ligand-receptor interactions from both simulated and real datasets, DiSiR outperforms other well-known permutation-based methods, e.g. CellPhoneDB and ICELLNET. Finally, to demonstrate DiSiR's utility in exploring data and generating biologically relevant hypotheses, we apply it to COVID lung and rheumatoid arthritis (RA) synovium scRNA-seq datasets and highlight potential differences between inflammatory pathways at cell type level for control versus disease samples.

The integration of large-scale public data and network analysis uncovers molecular characteristics of psoriasis.

In recent years, a growing interest in the characterization of the molecular basis of psoriasis has been observed. However, despite the availability of a large amount of molecular data, many pathogenic mechanisms of psoriasis are still poorly understood. In this study, we performed an integrated analysis of 23 public transcriptomic datasets encompassing both lesional and uninvolved skin samples from psoriasis patients. We defined comprehensive gene co-expression network models of psoriatic lesions and uninvolved skin. Moreover, we curated and exploited a wide range of functional information from multiple public sources in order to systematically annotate the inferred networks. The integrated analysis of transcriptomics data and co-expression networks highlighted genes that are frequently dysregulated and show aberrant patterns of connectivity in the psoriatic lesion compared with the unaffected skin. Our approach allowed us to also identify plausible, previously unknown, actors in the expression of the psoriasis phenotype. Finally, we characterized communities of co-expressed genes associated with relevant molecular functions and expression signatures of specific immune cell types associated with the psoriasis lesion. Overall, integrating experimental driven results with curated functional information from public repositories represents an efficient approach to empower knowledge generation about psoriasis and may be applicable to other complex diseases.

Vitiligo Skin T Cells Are Prone to Produce Type 1 and Type 2 Cytokines to Induce Melanocyte Dysfunction and Epidermal Inflammatory Response Through Jak Signaling.

Vitiligo is a T cell-mediated inflammatory skin disorder characterized by the loss of epidermal melanocytes. However, the contribution of melanocytes to the physiopathology of the disease in response to the T-cell microenvironment remains unclear. Here, using NanoString technology and multiplex ELISA, we show that active vitiligo perilesional skin is characterized by prominent type 1 and 2 associated immune responses. The vitiligo skin T-cell secretome downregulated melanocyte function and adhesion while increasing melanocyte mitochondrial metabolism and expression of inflammatory cytokines and chemokines by epidermal cells. The Jak1/2 inhibitor ruxolitinib strongly inhibited such effects on epidermal cells. Our data highlight that vitiligo is more complex than previously thought, with prominent combined activities of both T helper type 1- and T helper type 2-related cytokines inducing inflammatory responses of epidermal cells. Melanocytes do not appear only to be a target of T cells in vitiligo but could actively contribute to perpetuate inflammation. Jak inhibitors could prevent the impact of T cells on epidermal cells and pigmentation, highlighting their potential clinical benefit in vitiligo.

The Th1 cell regulatory circuitry is largely conserved between human and mouse.

Gene expression programs controlled by lineage-determining transcription factors are often conserved between species. However, infectious diseases have exerted profound evolutionary pressure, and therefore the genes regulated by immune-specific transcription factors might be expected to exhibit greater divergence. T-bet (Tbx21) is the immune-specific, lineage-specifying transcription factor for T helper type I (Th1) immunity, which is fundamental for the immune response to intracellular pathogens but also underlies inflammatory diseases. We compared T-bet genomic targets between mouse and human CD4+ T cells and correlated T-bet binding patterns with species-specific gene expression. Remarkably, we found that the majority of T-bet target genes are conserved between mouse and human, either via preservation of binding sites or via alternative binding sites associated with transposon-linked insertion. Species-specific T-bet binding was associated with differences in transcription factor-binding motifs and species-specific expression of associated genes. These results provide a genome-wide cross-species comparison of Th1 gene regulation that will enable more accurate translation of genetic targets and therapeutics from pre-clinical models of inflammatory and infectious diseases and cancer into human clinical trials.

Enhanced NF-κB signaling in type-2 dendritic cells at baseline predicts non-response to adalimumab in psoriasis.

Biologic therapies have transformed the management of psoriasis, but clinical outcome is variable leaving an unmet clinical need for predictive biomarkers of response. Here we perform in-depth immunomonitoring of blood immune cells of 67 patients with psoriasis, before and during therapy with the anti-TNF drug adalimumab, to identify immune mediators of clinical response and evaluate their predictive value. Enhanced NF-κBp65 phosphorylation, induced by TNF and LPS in type-2 dendritic cells (DC) before therapy, significantly correlates with lack of clinical response after 12 weeks of treatment. The heightened NF-κB activation is linked to increased DC maturation in vitro and frequency of IL-17+ T cells in the blood of non-responders before therapy. Moreover, lesional skin of non-responders contains higher numbers of dermal DC expressing the maturation marker CD83 and producing IL-23, and increased numbers of IL-17+ T cells. Finally, we identify and clinically validate LPS-induced NF-κBp65 phosphorylation before therapy as a predictive biomarker of non-response to adalimumab, with 100% sensitivity and 90.1% specificity in an independent cohort. Our study uncovers important molecular and cellular mediators underpinning adalimumab mechanisms of action in psoriasis and we propose a blood biomarker for predicting clinical outcome.

Publisher Correction: An optimized workflow for single-cell transcriptomics and repertoire profiling of purified lymphocytes from clinical samples.

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

An optimized workflow for single-cell transcriptomics and repertoire profiling of purified lymphocytes from clinical samples.

Establishing clinically relevant single-cell (SC) transcriptomic workflows from cryopreserved tissue is essential to move this emerging immune monitoring technology from the bench to the bedside. Improper sample preparation leads to detrimental cascades, resulting in loss of precious time, money and finally compromised data. There is an urgent need to establish protocols specifically designed to overcome the inevitable variations in sample quality resulting from uncontrollable factors in a clinical setting. Here, we explore sample preparation techniques relevant to a range of clinically relevant scenarios, where SC gene expression and repertoire analysis are applied to a cryopreserved sample derived from a small amount of blood, with unknown or partially known preservation history. We compare a total of ten cell-counting, viability-improvement, and lymphocyte-enrichment methods to highlight a number of unexpected findings. Trypan blue-based automated counters, typically recommended for single-cell sample quantitation, consistently overestimate viability. Advanced sample clean-up procedures significantly impact total cell yield, while only modestly increasing viability. Finally, while pre-enrichment of B cells from whole peripheral blood mononuclear cells (PBMCs) results in the most reliable BCR repertoire data, comparable T-cell enrichment strategies distort the ratio of CD4+ and CD8+ cells. Furthermore, we provide high-resolution analysis of gene expression and clonotype repertoire of different B cell subtypes. Together these observations provide both qualitative and quantitative sample preparation guidelines that increase the chances of obtaining high-quality single-cell transcriptomic and repertoire data from human PBMCs in a variety of clinical settings.

Phase I/II open-label trial of intravenous allogeneic mesenchymal stromal cell therapy in adults with recessive dystrophic epidermolysis bullosa.

BACKGROUND: Recessive dystrophic epidermolysis bullosa (RDEB) is a hereditary blistering disorder due to a lack of type VII collagen. At present, treatment is mainly supportive. OBJECTIVES: To determine whether intravenous allogeneic bone marrow-derived mesenchymal stromal/stem cells (BM-MSCs) are safe in RDEB adults and if the cells improve wound healing and quality of life. METHODS: We conducted a prospective, phase I/II, open-label study recruiting 10 RDEB adults to receive 2 intravenous infusions of BM-MSCs (on day 0 and day 14; each dose 2-4 × 106 cells/kg). RESULTS: BM-MSCs were well tolerated with no serious adverse events to 12 months. Regarding efficacy, there was a transient reduction in disease activity scores (8/10 subjects) and a significant reduction in itch. One individual showed a transient increase in type VII collagen. LIMITATIONS: Open-label trial with no placebo. CONCLUSIONS: MSC infusion is safe in RDEB adults and can have clinical benefits for at least 2 months.

Interplay between the human gut microbiome and host metabolism.

The human gut is inhabited by a complex and metabolically active microbial ecosystem. While many studies focused on the effect of individual microbial taxa on human health, their overall metabolic potential has been under-explored. Using whole-metagenome shotgun sequencing data in 1,004 twins, we first observed that unrelated subjects share, on average, almost double the number of metabolic pathways (82%) than species (43%). Then, using 673 blood and 713 faecal metabolites, we found metabolic pathways to be associated with 34% of blood and 95% of faecal metabolites, with over 18,000 significant associations, while species showed less than 3,000 associations. Finally, we estimated that the microbiome was involved in a dialogue between 71% of faecal, and 15% of blood, metabolites. This study underlines the importance of studying the microbial metabolic potential rather than focusing purely on taxonomy to find therapeutic and diagnostic targets, and provides a unique resource describing the interplay between the microbiome and the systemic and faecal metabolic environments.

Genome-wide association study in frontal fibrosing alopecia identifies four susceptibility loci including HLA-B*07:02.

Frontal fibrosing alopecia (FFA) is a recently described inflammatory and scarring type of hair loss affecting almost exclusively women. Despite a dramatic recent increase in incidence the aetiopathogenesis of FFA remains unknown. We undertake genome-wide association studies in females from a UK cohort, comprising 844 cases and 3,760 controls, a Spanish cohort of 172 cases and 385 controls, and perform statistical meta-analysis. We observe genome-wide significant association with FFA at four genomic loci: 2p22.2, 6p21.1, 8q24.22 and 15q2.1. Within the 6p21.1 locus, fine-mapping indicates that the association is driven by the HLA-B*07:02 allele. At 2p22.1, we implicate a putative causal missense variant in CYP1B1, encoding the homonymous xenobiotic- and hormone-processing enzyme. Transcriptomic analysis of affected scalp tissue highlights overrepresentation of transcripts encoding components of innate and adaptive immune response pathways. These findings provide insight into disease pathogenesis and characterise FFA as a genetically predisposed immuno-inflammatory disorder driven by HLA-B*07:02.

The cholesterol biosynthesis pathway regulates IL-10 expression in human Th1 cells.

The mechanisms controlling CD4+ T cell switching from an effector to an anti-inflammatory (IL-10+) phenotype play an important role in the persistence of chronic inflammatory diseases. Here, we identify the cholesterol biosynthesis pathway as a key regulator of this process. Pathway analysis of cultured cytokine-producing human T cells reveals a significant association between IL-10 and cholesterol metabolism gene expression. Inhibition of the cholesterol biosynthesis pathway with atorvastatin or 25-hydroxycholesterol during switching from IFNγ+ to IL-10+ shows a specific block in immune resolution, defined as a significant decrease in IL-10 expression. Mechanistically, the master transcriptional regulator of IL10 in T cells, c-Maf, is significantly decreased by physiological levels of 25-hydroxycholesterol. Strikingly, progression to rheumatoid arthritis is associated with altered expression of cholesterol biosynthesis genes in synovial biopsies of predisposed individuals. Our data reveal a link between sterol metabolism and the regulation of the anti-inflammatory response in human CD4+ T cells.