Claudin-10 Expression and the Gene Expression Pattern of Thick Ascending Limb Cells

Many genomic, anatomical and functional differences exist between the medullary (MTAL) and the cortical thick ascending limb of the loop of Henle (CTAL), including a higher expression of claudin-10 (CLDN10) in the MTAL than in the CTAL. Therefore, we assessed to what extent the Cldn10 gene expression is a determinant of differential gene expression between MTAL and CTAL. RNAs extracted from CTAL and MTAL microdissected from wild type (WT) and Cldn10 knock out mice (cKO) were analyzed by RNAseq. Differential and enrichment analyses (GSEA) were performed with interactive R Shiny software. Between WT and cKO MTAL, 637 genes were differentially expressed, whereas only 76 were differentially expressed between WT and cKO CTAL. Gene expression patterns and GSEA analyses in all replicates showed that WT MTAL did not cluster with the other replicates; no hierarchical clustering could be found between WT CTAL, cKO CTAL and cKO MTAL. Compared to WT replicates, cKO replicates were enriched in Cldn16, Cldn19, Pth1r, (parathyroid hormone receptor type 1), Casr (calcium sensing receptor) and Vdr (Vitamin D Receptor) mRNA in both the cortex and medulla. Cldn10 is associated with gene expression patterns, including genes specifically involved in divalent cations reabsorption in the TAL.

Conserved plant transcriptional responses to microgravity from two consecutive spaceflight experiments

Introduction

Understanding how plants adapt to the space environment is essential, as plants will be a valuable component of long duration space missions. Several spaceflight experiments have focused on transcriptional profiling as a means of understanding plant adaptation to microgravity. However, there is limited overlap between results from different experiments. Differences in experimental conditions and hardware make it difficult to find a consistent response across experiments and to distinguish the primary effects of microgravity from other spaceflight effects.

Methods

Plant Signaling (PS) and Plant RNA Regulation (PRR) were two separate spaceflight experiments conducted on the International Space Station utilizing the European Modular Cultivation System (EMCS). The EMCS provided a lighted environment for plant growth with centrifugal capabilities providing an onboard 1 g control.

Results and discussion

An RNA-Seq analysis of shoot samples from PS and PRR revealed a significant overlap of genes differentially expressed in microgravity between the two experiments. Relative to onboard 1 g controls, genes involved in transcriptional regulation, shoot development, and response to auxin and light were upregulated in microgravity in both experiments. Conversely, genes involved in defense response, abiotic stress, Ca++ signaling, and cell wall modification were commonly downregulated in both datasets. The downregulation of stress responses in microgravity in these two experiments is interesting as these pathways have been previously observed as upregulated in spaceflight compared to ground controls. Similarly, we have observed many stress response genes to be upregulated in the 1 g onboard control compared to ground reference controls; however these genes were specifically downregulated in microgravity. In addition, we analyzed the sRNA landscape of the 1 g and microgravity (μ g) shoot samples from PRR. We identified three miRNAs (miR319c, miR398b, and miR8683) which were upregulated in microgravity, while several of their corresponding target genes were found to be downregulated in microgravity. Interestingly, the downregulated target genes are enriched in those encoding chloroplast-localized enzymes and proteins. These results uncover microgravity unique transcriptional changes and highlight the validity and importance of an onboard 1 g control.

Patterns of metastatic progression and association with clinical outcomes in canine osteosarcoma: A necropsy study of 83 dogs

Osteosarcoma is a highly metastatic primary bone tumour that occurs spontaneously in both pet dogs and humans. Patterns of metastasis to organs beyond the most common site (lung) are poorly characterised and it is unknown whether specific associations between patterns of metastatic progression and patient features exist. This retrospective study characterised the necropsy findings of 83 dogs receiving standardised therapy and clinical monitoring in a prospective clinical trial setting to document patterns of metastasis and correlate outcomes with these patterns and other patient and tumour-specific factors. A total of 20 different sites of metastasis were documented, with lung as the most common site, followed by bone, kidney, liver, and heart. Two distinct clusters of dogs were identified based on patterns of metastasis. There was no significant association between site of enrollment, trial arm, sex, serum alkaline phosphatase (ALP) activity, or tumour location and clinical outcomes. A second cancer type was identified at necropsy in 10 dogs (10/83; 12%). These data showcase the extensive nature of osteosarcoma metastasis beyond the lung and provide a benchmark for clinical monitoring of the disease. Further, this study provides insight into transcriptional features of primary tumours that may relate to a propensity for osteosarcoma metastasis to specific organs and tissues.

Elucidating the Proximal Tubule HNF4A Gene Regulatory Network in Human Kidney Organoids

SIGNIFICANCE STATEMENT

HNF4 genes promote proximal tubule differentiation in mice, but their function in human nephrogenesis is not fully defined. This study uses human pluripotent stem cell (PSC)-derived kidney organoids as a model to investigate HNF4A and HNF4G functions. The loss of HNF4A , but not HNF4G , impaired reabsorption-related molecule expression and microvilli formation in human proximal tubules. Cleavage under targets and release using nuclease (CUT&RUN) sequencing and CRISPR-mediated transcriptional activation (CRISPRa) further confirm that HNF4A directly regulates its target genes. Human kidney organoids provide a good model for studying transcriptional regulation in human kidney development.

BACKGROUND

The proximal tubule plays a major role in electrolyte homeostasis. Previous studies have shown that HNF4A regulates reabsorption-related genes and promotes proximal tubule differentiation during murine kidney development. However, the functions and gene regulatory mechanisms of HNF4 family genes in human nephrogenesis have not yet been investigated.

METHODS

We generated HNF4A -knock out (KO), HNF4G -KO, and HNF4A/4G -double KO human pluripotent stem cell lines, differentiated each into kidney organoids, and used immunofluorescence analysis, electron microscopy, and RNA-seq to analyze them. We probed HNF4A-binding sites genome-wide by cleavage under targets and release using nuclease sequencing in both human adult kidneys and kidney organoid-derived proximal tubular cells. Clustered Regularly Interspaced Short Palindromic Repeats-mediated transcriptional activation validated HNF4A and HNF4G function in proximal tubules during kidney organoid differentiation.

RESULTS

Organoids lacking HNF4A , but not HNF4G , showed reduced expression of transport-related, endocytosis-related, and brush border-related genes, as well as disorganized brush border structure in the apical lumen of the organoid proximal tubule. Cleavage under targets and release using nuclease revealed that HNF4A primarily bound promoters and enhancers of genes that were downregulated in HNF4A -KO, suggesting direct regulation. Induced expression of HNF4A or HNF4G by CRISPR-mediated transcriptional activation drove increased expression of selected target genes during kidney organoid differentiation.

CONCLUSIONS

This study reveals regulatory mechanisms of HNF4A and HNF4G during human proximal tubule differentiation. The experimental strategy can be applied more broadly to investigate transcriptional regulation in human kidney development.

Genome-wide translational response of Candida albicans to fluconazole treatment

Azoles are commonly used for the treatment of fungal infections, and the ability of human fungal pathogens to rapidly respond to azole treatment is critical for the development of antifungal resistance. While the roles of genetic mutations, chromosomal rearrangements, and transcriptional mechanisms in azole resistance have been well-characterized, very little is known about post-transcriptional and translational mechanisms that drive this process. In addition, most previous genome-wide studies have focused on transcriptional responses to azole treatment and likely serve as inaccurate proxies for changes in protein expression due to extensive post-transcriptional and translational regulation. In this study, we use ribosome profiling to provide the first picture of the global translational response of a major human fungal pathogen, Candida albicans, to treatment with fluconazole (Flu), one of the most widely used azole drugs. We identify sets of genes showing significantly altered translational efficiency, including genes associated with a variety of biological processes such as the cell cycle, DNA repair, cell wall/cell membrane biosynthesis, transport, signaling, DNA- and RNA-binding activities, and protein synthesis. We observe both similarities and differences among the most highly represented gene categories (as defined by gene ontology) that are regulated by fluconazole at the translational vs transcriptional levels. Importantly, however, very few genes that are translationally regulated by fluconazole are also controlled transcriptionally under this condition. Our findings suggest that C. albicans possesses distinct translational mechanisms that are important for the response to antifungal treatment, which could eventually be targeted by novel antifungal therapies. IMPORTANCE Azoles are one of the most commonly used drug classes to treat human fungal pathogens. While point mutations, chromosomal rearrangements, and transcriptional mechanisms that drive azole resistance have been well-characterized, we know very little about the role of translational mechanisms. In this study, we determined the global translational profile of genes that are expressed in the major human fungal pathogen Candida albicans in response to fluconazole, one of the most widely used azole drugs. We find both similarities and differences among the most highly represented categories of genes regulated by fluconazole at the transcriptional and translational levels. Interestingly, however, many of the specific genes that are regulated by fluconazole at the translational level do not appear to be controlled by transcriptional mechanisms under this condition. Our results suggest that distinct C. albicans translational mechanisms control the response to antifungals and could eventually be targeted in the development of new therapies.

Neurodegenerative Disease Associated Pathways in the Brains of Triple Transgenic Alzheimer's Model Mice Are Reversed Following Two Weeks of Peripheral Administration of Fasudil

The pan Rho-associated coiled-coil-containing protein kinase (ROCK) inhibitor fasudil acts as a vasodilator and has been used as a medication for post-cerebral stroke for the past 29 years in Japan and China. More recently, based on the involvement of ROCK inhibition in synaptic function, neuronal survival, and processes associated with neuroinflammation, it has been suggested that the drug may be repurposed for neurodegenerative diseases. Indeed, fasudil has demonstrated preclinical efficacy in many neurodegenerative disease models. To facilitate an understanding of the wider biological processes at play due to ROCK inhibition in the context of neurodegeneration, we performed a global gene expression analysis on the brains of Alzheimer's disease model mice treated with fasudil via peripheral IP injection. We then performed a comparative analysis of the fasudil-driven transcriptional profile with profiles generated from a meta-analysis of multiple neurodegenerative diseases. Our results show that fasudil tends to drive gene expression in a reverse sense to that seen in brains with post-mortem neurodegenerative disease. The results are most striking in terms of pathway enrichment analysis, where pathways perturbed in Alzheimer's and Parkinson's diseases are overwhelmingly driven in the opposite direction by fasudil treatment. Thus, our results bolster the repurposing potential of fasudil by demonstrating an anti-neurodegenerative phenotype in a disease context and highlight the potential of in vivo transcriptional profiling of drug activity.

Activation of Angiopoietin-Tie2 Signaling Protects the Kidney from Ischemic Injury by Modulation of Endothelial-Specific Pathways

SIGNIFICANCE STATEMENT

Ischemia-reperfusion AKI (IR-AKI) is common and causes significant morbidity. Effective treatments are lacking. However, preclinical studies suggest that inhibition of angiopoietin-Tie2 vascular signaling promotes injury, whereas activation of Tie2 is protective. We show that kidney ischemia leads to increased levels of the endothelial-specific phosphatase vascular endothelial protein tyrosine phosphatase (VE-PTP; PTPRB), which inactivates Tie2. Activation of Tie2 through VE-PTP deletion, or delivery of a novel angiopoietin mimetic (Hepta-ANG1), abrogated IR-AKI in mice. Single-cell RNAseq analysis showed Tie2 activation promotes increased Entpd1 expression, downregulation of FOXO1 target genes in the kidney vasculature, and emergence of a new subpopulation of glomerular endothelial cells. Our data provide a molecular basis and identify a candidate therapeutic to improve endothelial integrity and kidney function after IR-AKI.

BACKGROUND

Ischemia-reperfusion AKI (IR-AKI) is estimated to affect 2%-7% of all hospitalized patients. The significant morbidity and mortality associated with AKI indicates urgent need for effective treatments. Previous studies have shown activation of the vascular angiopoietin-Tie2 tyrosine kinase signaling pathway abrogates ischemia-reperfusion injury (IRI). We extended previous studies to (1) determine the molecular mechanism(s) underlying kidney injury and protection related to decreased or increased activation of Tie2, respectively, and (2) to test the hypothesis that deletion of the Tie2 inhibitory phosphatase vascular endothelial protein tyrosine phosphatase (VE-PTP) or injection of a new angiopoietin mimetic protects the kidney from IRI by common molecular mechanism(s).

METHODS

Bilateral IR-AKI was performed in VE-PTP wild-type or knockout mice and in C57BL/6J mice treated with Hepta-ANG1 or vehicle. Histologic, immunostaining, and single-cell RNA sequencing analyses were performed.

RESULTS

The phosphatase VE-PTP, which negatively regulates the angiopoietin-Tie2 pathway, was upregulated in kidney endothelial cells after IRI, and genetic deletion of VE-PTP in mice protected the kidney from IR-AKI. Injection of Hepta-ANG1 potently activated Tie2 and protected the mouse kidney from IRI. Single-cell RNAseq analysis of kidneys from Hepta-ANG1-treated and vehicle-treated mice identified endothelial-specific gene signatures and emergence of a new glomerular endothelial subpopulation associated with improved kidney function. Overlap was found between endothelial-specific genes upregulated by Hepta-ANG1 treatment and those downregulated in HUVECs with constitutive FOXO1 activation, including Entpd1 / ENTPD1 that modulates purinergic receptor signaling.

CONCLUSIONS

Our data support a key role of the endothelium in the development of IR-AKI, introduce Hepta-ANG1 as a putative new therapeutic biologic, and report a model to explain how IRI reduces Tie2 signaling and how Tie2 activation protects the kidney.

PODCAST

This article contains a podcast at https://dts.podtrac.com/redirect.mp3/www.asn-online.org/media/podcast/JASN/2023_05_23_JSN_Ang_EP23_052323.mp3.

Butt-seq: a new method for facile profiling of transcription

A wide range of sequencing methods has been developed to assess nascent RNA transcription and resolve the single-nucleotide position of RNA polymerase genome-wide. These techniques are often burdened with high input material requirements and lengthy protocols. We leveraged the template-switching properties of thermostable group II intron reverse transcriptase (TGIRT) and developed Butt-seq (bulk analysis of nascent transcript termini sequencing), which can produce libraries from purified nascent RNA in 6 h and from as few as 10,000 cells-an improvement of at least 10-fold over existing techniques. Butt-seq shows that inhibition of the superelongation complex (SEC) causes promoter-proximal pausing to move upstream in a fashion correlated with subnucleosomal fragments. To address transcriptional regulation in a tissue, Butt-seq was used to measure the circadian regulation of transcription from fly heads. All the results indicate that Butt-seq is a simple and powerful technique to analyze transcription at a high level of resolution.

Lineage Tracing and Single-Nucleus Multiomics Reveal Novel Features of Adaptive and Maladaptive Repair after Acute Kidney Injury

SIGNIFICANCE STATEMENT

Understanding the mechanisms underlying adaptive and maladaptive renal repair after AKI and their long-term consequences is critical to kidney health. The authors used lineage tracing of cycling cells and single-nucleus multiomics (profiling transcriptome and chromatin accessibility) after AKI. They demonstrated that AKI triggers a cell-cycle response in most epithelial and nonepithelial kidney cell types. They also showed that maladaptive proinflammatory proximal tubule cells (PTCs) persist until 6 months post-AKI, although they decreased in abundance over time, in part, through cell death. Single-nucleus multiomics of lineage-traced cells revealed regulatory features of adaptive and maladaptive repair. These included activation of cell state-specific transcription factors and cis-regulatory elements, and effects in PTCs even after adaptive repair, weeks after the injury event.

BACKGROUND

AKI triggers a proliferative response as part of an intrinsic cellular repair program, which can lead to adaptive renal repair, restoring kidney structure and function, or maladaptive repair with the persistence of injured proximal tubule cells (PTCs) and an altered kidney structure. However, the cellular and molecular understanding of these repair programs is limited.

METHODS

To examine chromatin and transcriptional responses in the same cell upon ischemia-reperfusion injury (IRI), we combined genetic fate mapping of cycling ( Ki67+ ) cells labeled early after IRI with single-nucleus multiomics-profiling transcriptome and chromatin accessibility in the same nucleus-and generated a dataset of 83,315 nuclei.

RESULTS

AKI triggered a broad cell cycle response preceded by cell type-specific and global transcriptional changes in the nephron, the collecting and vascular systems, and stromal and immune cell types. We observed a heterogeneous population of maladaptive PTCs throughout proximal tubule segments 6 months post-AKI, with a marked loss of maladaptive cells from 4 weeks to 6 months. Gene expression and chromatin accessibility profiling in the same nuclei highlighted differences between adaptive and maladaptive PTCs in the activity of cis-regulatory elements and transcription factors, accompanied by corresponding changes in target gene expression. Adaptive repair was associated with reduced expression of genes encoding transmembrane transport proteins essential to kidney function.

CONCLUSIONS

Analysis of genome organization and gene activity with single-cell resolution using lineage tracing and single-nucleus multiomics offers new insight into the regulation of renal injury repair. Weeks to months after mild-to-moderate IRI, maladaptive PTCs persist with an aberrant epigenetic landscape, and PTCs exhibit an altered transcriptional profile even following adaptive repair.

Transcriptional Profiling of Rat Prefrontal Cortex after Acute Inescapable Footshock Stress

Stress is a primary risk factor for psychiatric disorders such as Major Depressive Disorder (MDD) and Post Traumatic Stress Disorder (PTSD). The response to stress involves the regulation of transcriptional programs, which is supposed to play a role in coping with stress. To evaluate transcriptional processes implemented after exposure to unavoidable traumatic stress, we applied microarray expression analysis to the PFC of rats exposed to acute footshock (FS) stress that were sacrificed immediately after the 40 min session or 2 h or 24 h after. While no substantial changes were observed at the single gene level immediately after the stress session, gene set enrichment analysis showed alterations in neuronal pathways associated with glia development, glia-neuron networking, and synaptic function. Furthermore, we found alterations in the expression of gene sets regulated by specific transcription factors that could represent master regulators of the acute stress response. Of note, these pathways and transcriptional programs are activated during the early stress response (immediately after FS) and are already turned off after 2 h-while at 24 h, the transcriptional profile is largely unaffected. Overall, our analysis provided a transcriptional landscape of the early changes triggered by acute unavoidable FS stress in the PFC of rats, suggesting that the transcriptional wave is fast and mild, but probably enough to activate a cellular response to acute stress.

Digital Spatial Profiling of Glomerular Gene Expression in Pauci-Immune Focal Necrotizing Glomerulonephritis

Pauci-immune focal necrotizing glomerulonephritis (piFNGN) involves asynchronous onset and progression of injurious lesions in biopsies. Pathologists can describe this heterogeneity within a biopsy, but translating the information into prognostic or expression analyses is challenging. Understanding the underlying molecular processes could improve treatment; however, bulk or single-cell transcriptomic analyses of dissociated tissue disregard the heterogeneity of glomerular injury. We characterize protein and mRNA expression of individual glomeruli in 20 biopsies from 18 patients with antineutrophil cytoplasmic antibody-associated piFNGN using the NanoString digital spatial profiling (DSP) platform. For this purpose, circular annotations of glomeruli were analyzed using protein, immuno-oncology RNA, and Cancer Transcriptome Atlas panels (n=120, 72, and 48 glomeruli, respectively). Histologic evaluation of glomerular patterns of injury was performed in adjacent serial sections. Expression data were processed by log2 transformation, quantile normalization, and batch adjustment. DSP revealed distinct but overlapping gene expression profiles relating to the morphological evolution of injurious lesions, including dynamic expression of various immune checkpoint regulators. Enrichment analysis indicated deregulated pathways that underline known and highlight novel potential mechanisms of disease. Moreover, by capturing individual glomeruli, DSP describes heterogeneity between and within biopsies. We demonstrate the benefit of spatial profiling for characterization of heterogeneous glomerular injury, indicating novel molecular correlates of glomerular injury in piFNGN.

Transcriptional profiling of Pseudomonas aeruginosa mature single- and dual-species biofilms in response to meropenem

Pseudomonas aeruginosa is a Gram-negative opportunistic pathogen frequently isolated from chronic infections of the cystic fibrosis lung and burn wounds, and is a major cause of antimicrobial-resistant nosocomial infections. P. aeruginosa is frequently co-isolated with the opportunistic fungal pathogen Candida albicans, with the presence of C. albicans in dual-species biofilms promoting tolerance to meropenem. Here, transcription profiling of mature P. aeruginosa single- or dual-species biofilms was carried out to understand the molecular mechanism(s) by which C. albicans enhances meropenem tolerance. C. albicans appeared to have a mild impact on the transcriptome of P. aeruginosa mature biofilms, with most differentially regulated genes being involved in interkingdom interactions (i.e. quorum sensing and phenazine biosynthesis). The addition of meropenem to mature single- or dual-species biofilms resulted in a significant bacterial transcriptional response, including the induction of the beta-lactamase, ampC, genes involved in biofilm formation. P. aeruginosa elicited a similar transcriptional response to meropenem in the presence of C. albicans, but C. albicans promoted the expression of additional efflux pumps, which could play roles in increasing the tolerance of P. aeruginosa to meropenem.

Endothelial Cell-Specific Molecule-1 Inhibits Albuminuria in Diabetic Mice

Background

Diabetic kidney disease (DKD) is the most common cause of kidney failure in the world, and novel predictive biomarkers and molecular mechanisms of disease are needed. Endothelial cell-specific molecule-1 (Esm-1) is a secreted proteoglycan that attenuates inflammation. We previously identified that a glomerular deficiency of Esm-1 associates with more pronounced albuminuria and glomerular inflammation in DKD-susceptible relative to DKD-resistant mice, but its contribution to DKD remains unexplored.

Methods

Using hydrodynamic tail-vein injection, we overexpress Esm-1 in DKD-susceptible DBA/2 mice and delete Esm-1 in DKD-resistant C57BL/6 mice to study the contribution of Esm-1 to DKD. We analyze clinical indices of DKD, leukocyte infiltration, podocytopenia, and extracellular matrix production. We also study transcriptomic changes to assess potential mechanisms of Esm-1 in glomeruli.

Results

In DKD-susceptible mice, Esm-1 inversely correlates with albuminuria and glomerular leukocyte infiltration. We show that overexpression of Esm-1 reduces albuminuria and diabetes-induced podocyte injury, independent of changes in leukocyte infiltration. Using a complementary approach, we find that constitutive deletion of Esm-1 in DKD-resistant mice modestly increases the degree of diabetes-induced albuminuria versus wild-type controls. By glomerular RNAseq, we identify that Esm-1 attenuates expression of kidney disease-promoting and interferon (IFN)-related genes, including Ackr2 and Cxcl11.

Conclusions

We demonstrate that, in DKD-susceptible mice, Esm-1 protects against diabetes-induced albuminuria and podocytopathy, possibly through select IFN signaling. Companion studies in patients with diabetes suggest a role of Esm-1 in human DKD.

Article title: Transcriptional profiling efficacy to define biological activity similarity for cosmetic ingredients' safety assessment based on next-generation read-across

The objective of this work was to use transcriptional profiling to assess the biological activity of structurally related chemicals to define their biological similarity and with that, substantiate the validity of a read-across approach usable in risk assessment. Two case studies are presented, one with 4 short alkyl chain parabens: methyl (MP), ethyl (EP), butyl (BP), and propylparaben (PP), as well as their main metabolite, p-hydroxybenzoic acid (pHBA) with the assumption that propylparaben was the target chemical; and a second one with caffeine and its main metabolites theophylline, theobromine and paraxanthine where CA was the target chemical. The comprehensive transcriptional response of MCF7, HepG2, A549 and ICell cardiomyocytes was evaluated (TempO-Seq) after exposure to vehicle-control, each paraben or pHBA, CA or its metabolites, at 3 non-cytotoxic concentrations, for 6 h. Differentially expressed genes (FDR ≥0.05, and fold change ±1.2≥) were identified for each chemical, at each concentration, and used to determine similarities. Each of the chemicals is able to elicit changes in the expression of a number of genes, as compared to controls. Importantly, the transcriptional profile elicited by each of the parabens shares a high degree of similarity across the group. The highest number of genes commonly affected was between butylparaben and PP. The transcriptional profile of the parabens is similar to the one elicited by estrogen receptor agonists, with BP being the closest structural and biological analogue for PP. In the CA case, the transcriptional profile elicited of all four methylxanthines had a high degree of similarity across the cell types, with CA and theophylline being the most active. The most robust response was obtained in the cardiomyocytes with the highest transcriptional profile similarity between CA and TP. The transcriptional profile of the methylxanthines is similar to the one elicited by inhibitors of phosphatidylinositol 3-kinase as well as other kinase inhibitors. Overall, our results support the approach of incorporating transcriptional profiling in well-designed in vitro tests as one robust stream of data to support biological similarity driven read-across procedures and strengthening the traditional structure-based approaches useful in risk assessment.

Distinct gene expression dynamics in developing and regenerating crustacean limbs

Regenerating animals have the ability to reproduce body parts that were originally made in the embryo and subsequently lost due to injury. Understanding whether regeneration mirrors development is an open question in most regenerative species. Here, we take a transcriptomics approach to examine whether leg regeneration shows similar temporal patterns of gene expression as leg development in the embryo, in the crustacean Parhyale hawaiensis. We find that leg development in the embryo shows stereotypic temporal patterns of gene expression. In contrast, the dynamics of gene expression during leg regeneration show a higher degree of variation related to the physiology of individual animals. A major driver of this variation is the molting cycle. We dissect the transcriptional signals of individual physiology and regeneration to obtain clearer temporal signals marking distinct phases of leg regeneration. Comparing the transcriptional dynamics of development and regeneration we find that, although the two processes use similar sets of genes, the temporal patterns in which these genes are deployed are different and cannot be systematically aligned.

Integrative Analysis of Genes Involved in the Global Response to Potato Wart Formation

Synchytrium endobioticum, the causal agent of potato wart disease, poses a major threat to commercial potato production. Understanding the roles of transcriptionally regulated genes following pathogen infection is necessary for understanding the system-level host response to pathogen. Although some understanding of defense mechanisms against S. endobioticum infection has been gained for incompatible interactions, the genes and signaling pathways involved in the compatible interaction remain unclear. Based on the collection of wart diseased tubers of a susceptible cultivar, we performed phenotypic and dual RNA-Seq analyses of wart lesions in seven stages of disease progression. We totally detected 5,052 differentially expressed genes (DEGs) by comparing the different stages of infection to uninfected controls. The tendency toward differential gene expression was active rather than suppressed under attack by the pathogen. The number of DEGs step-up along with the development of the disease and the first, third and seventh of the disease stages showed substantially increase of DEGs in comparison of the previous stage. The important functional groups identified via Gene ontology (GO) and KEGG enrichment were those responsible for plant-pathogen interaction, fatty acid elongation and phenylpropanoid biosynthesis. Gene coexpression networks, composed of 17 distinct gene modules that contained between 25 and 813 genes, revealed high interconnectivity of the induced response and led to the identification of a number of hub genes enriched at different stages of infection. These results provide a comprehensive perspective on the global response of potato to S. endobioticum infection and identify a potential transcriptional regulatory network underlying this susceptible response, which contribute to a better understanding of the potato-S. endobioticum pathosystem.

Post-trauma behavioral phenotype predicts the degree of vulnerability to fear relapse after extinction in male rats

Current treatments for trauma-related disorders remain ineffective for many patients.^1,2 Fear extinction deficiency is a prominent feature of these diseases,³ and many behavioral treatments rely on extinction training.^4,5 However, in many patients, therapy is followed by a relapse of symptoms, and the underpinnings of such interindividual variations in vulnerability to relapse remain unknown.^6-8 Here, we modeled interindividual differences in post-therapy fear relapse with an ethologically relevant trauma recovery paradigm. After fear conditioning, male rats underwent fear extinction while foraging in a large enriched arena, permitting the expression of a wide spectrum of behaviors. An automated multidimensional behavioral assessment revealed that post-conditioning fear response profiles clustered into two groups: some animals expressed fear by freezing more, whereas others darted more, as if fleeing from danger. Remarkably, the tendency of an animal to dart or to freeze after CS presentation during the first extinction session was, respectively, associated with stronger or weaker fear renewal. Moreover, genome-wide transcriptional profiling revealed that these groups differentially regulated specific sets of genes, some of which were previously implicated in anxiety and trauma-related disorders. Our results suggest that post-trauma behavioral phenotypes and the associated gene expression landscapes can serve as markers of fear relapse susceptibility and thus may be instrumental for future development of more effective treatments for psychiatric patients.

Molecular Characterization of Membranous Nephropathy

BACKGROUND

Molecular characterization of nephropathies may facilitate pathophysiologic insight, development of targeted therapeutics, and transcriptome-based disease classification. Although membranous nephropathy (MN) is a common cause of adult-onset nephrotic syndrome, the molecular pathways of kidney damage in MN require further definition.

METHODS

We applied a machine-learning framework to predict diagnosis on the basis of gene expression from the microdissected kidney tissue of participants in the Nephrotic Syndrome Study Network (NEPTUNE) cohort. We sought to identify differentially expressed genes between participants with MN versus those of other glomerulonephropathies across the NEPTUNE and European Renal cDNA Bank (ERCB) cohorts, to find MN-specific gene modules in a kidney-specific functional network, and to identify cell-type specificity of MN-specific genes using single-cell sequencing data from reference nephrectomy tissue.

RESULTS

Glomerular gene expression alone accurately separated participants with MN from those with other nephrotic syndrome etiologies. The top predictive classifier genes from NEPTUNE participants were also differentially expressed in the ERCB participants with MN. We identified a signature of 158 genes that are significantly differentially expressed in MN across both cohorts, finding 120 of these in a validation cohort. This signature is enriched in targets of transcription factor NF-κB. Clustering these MN-specific genes in a kidney-specific functional network uncovered modules with functional enrichments, including in ion transport, cell projection morphogenesis, regulation of adhesion, and wounding response. Expression data from reference nephrectomy tissue indicated 43% of these genes are most highly expressed by podocytes.

CONCLUSIONS

These results suggest that, relative to other glomerulonephropathies, MN has a distinctive molecular signature that includes upregulation of many podocyte-expressed genes, provides a molecular snapshot of MN, and facilitates insight into MN's underlying pathophysiology.

From environmental adaptation to host survival: Attributes that mediate pathogenicity of Candida auris

Candida species are a major cause of invasive fungal infections. While Candida albicans, C. glabrata, C. parapsilosis, and C. tropicalis are the most dominant species causing life-threatening candidiasis, C. auris recently emerged as a new species causing invasive infections with high rates of clinical treatment failures. To mimic initial phases of systemic Candida infections with dissemination via the bloodstream and to elucidate the pathogenic potential of C. auris, we used an ex vivo whole blood infection model. Similar to other clinically relevant Candida spp., C. auris is efficiently killed in human blood, but showed characteristic patterns of immune cell association, survival rates, and cytokine induction. Dual-species transcriptional profiling of C. auris-infected blood revealed a unique C. auris gene expression program during infection, while the host response proofed similar and conserved compared to other Candida species. C. auris-specific responses included adaptation and survival strategies, such as counteracting oxidative burst of immune cells, but also expression of potential virulence factors, (drug) transporters, and cell surface-associated genes. Despite comparable pathogenicity to other Candida species in our model, C. auris-specific transcriptional adaptations as well as its increased stress resistance and long-term environmental survival, likely contribute to the high risk of contamination and distribution in a nosocomial setting. Moreover, infections of neutrophils with pre-starved C. auris cells suggest that environmental preconditioning can have modulatory effects on the early host interaction. In summary, we present novel insights into C. auris pathogenicity, revealing adaptations to human blood and environmental niches distinctive from other Candida species.

Sea Anemones Responding to Sex Hormones, Oxybenzone, and Benzyl Butyl Phthalate: Transcriptional Profiling and in Silico Modelling Provide Clues to Decipher Endocrine Disruption in Cnidarians

Endocrine disruption is suspected in cnidarians, but questions remain how occurs. Steroid sex hormones are detected in corals and sea anemones even though these animals do not have estrogen receptors and their repertoire of steroidogenic enzymes appears to be incomplete. Pathways associated with sex hormone biosynthesis and sterol signaling are an understudied area in cnidarian biology. The objective of this study was to identify a suite of genes that can be linked to exposure of endocrine disruptors. Exaiptasia diaphana were exposed to nominal 20ppb concentrations of estradiol (E2), testosterone (T), cholesterol, oxybenzone (BP-3), or benzyl butyl phthalate (BBP) for 4 h. Eleven genes of interest (GOIs) were chosen from a previously generated EST library. The GOIs are 17β-hydroxysteroid dehydrogenases type 14 (17β HSD14) and type 12 (17β HSD12), Niemann-Pick C type 2 (NPC2), Equistatin (EI), Complement component C3 (C3), Cathepsin L (CTSL), Patched domain-containing protein 3 (PTCH3), Smoothened (SMO), Desert Hedgehog (DHH), Zinc finger protein GLI2 (GLI2), and Vitellogenin (VTG). These GOIs were selected because of functional associations with steroid hormone biosynthesis; cholesterol binding/transport; immunity; phagocytosis; or Hedgehog signaling. Quantitative Real-Time PCR quantified expression of GOIs. In silico modelling utilized protein structures from Protein Data Bank as well as creating protein structures with SWISS-MODEL. Results show transcription of steroidogenic enzymes, and cholesterol binding/transport proteins have similar transcription profiles for E2, T, and cholesterol treatments, but different profiles when BP-3 or BBP is present. C3 expression can differentiate between exposures to BP-3 versus BBP as well as exposure to cholesterol versus sex hormones. In silico modelling revealed all ligands (E2, T, cholesterol, BBP, and BP-3) have favorable binding affinities with 17β HSD14, 17β HSD12, NPC2, SMO, and PTCH proteins. VTG expression was down-regulated in the sterol treatments but up-regulated in BP-3 and BBP treatments. In summary, these eleven GOIs collectively generate unique transcriptional profiles capable of discriminating between the five chemical exposures used in this investigation. This suite of GOIs are candidate biomarkers for detecting transcriptional changes in steroidogenesis, gametogenesis, sterol transport, and Hedgehog signaling. Detection of disruptions in these pathways offers new insight into endocrine disruption in cnidarians.

Cubilin-, megalin-, and Dab2-dependent transcription revealed by CRISPR/Cas9 knockout in kidney proximal tubule cells

The multiligand receptors megalin (Lrp2) and cubilin (Cubn) and their endocytic adaptor protein Dab2 (Dab2) play essential roles in maintaining the integrity of the apical endocytic pathway of proximal tubule (PT) cells and have complex and poorly understood roles in the development of chronic kidney disease. Here, we used RNA-sequencing and CRISPR/Cas9 knockout (KO) technology in a well-differentiated cell culture model to identify PT-specific transcriptional changes that are directly consequent to the loss of megalin, cubilin, or Dab2 expression. KO of Lrp2 had the greatest transcriptional effect, and nearly all genes whose expression was affected in Cubn KO and Dab2 KO cells were also changed in Lrp2 KO cells. Pathway analysis and more granular inspection of the altered gene profiles suggested changes in pathways with immunomodulatory functions that might trigger the pathological changes observed in KO mice and patients with Donnai-Barrow syndrome. In addition, differences in transcription patterns between Lrp2 and Dab2 KO cells suggested the possibility that altered spatial signaling by aberrantly localized receptors contributes to transcriptional changes upon the disruption of PT endocytic function. A reduction in transcripts encoding sodium-glucose cotransporter isoform 2 was confirmed in Lrp2 KO mouse kidney lysates by quantitative PCR analysis. Our results highlight the role of megalin as a master regulator and coordinator of ion transport, metabolism, and endocytosis in the PT. Compared with the studies in animal models, this approach provides a means to identify PT-specific transcriptional changes that are directly consequent to the loss of these target genes.NEW & NOTEWORTHY Megalin and cubilin receptors together with their adaptor protein Dab2 represent major components of the endocytic machinery responsible for efficient uptake of filtered proteins by the proximal tubule (PT). Dab2 and megalin expression have been implicated as both positive and negative modulators of kidney disease. We used RNA sequencing to knock out CRISPR/Cas9 cubilin, megalin, and Dab2 in highly differentiated PT cells to identify PT-specific changes that are directly consequent to knockout of each component.

Transcriptional profiling of the chicken tracheal and splenic response to virulent Mycoplasma synoviae

Mycoplasma synoviae (MS), an important avian pathogen, can cause chronic respiratory disease, eggshell apex abnormalities, infectious synovitis, and arthritis in avian species, leading serious economic losses in the global poultry industry. To date, studies have shown significant different transcript profiles using various chicken cells after MS infection. However, in vitro cell models cannot fully represent the complex in vivo regulations after adventitious infection. The objective of this study was to explore the nature of the host-pathogen interaction during MS infection. The tracheal and spleen tissues of chickens were collected at d 0, 1, 3, and 5 postinoculation, and samples were analyzed for differential gene expression using Illumina RNA sequencing. A lot of significantly differentially expressed genes (DEGs) were observed in this analysis, and 861 DEGs were observed in trachea tissues and 753 DEGs were observed in spleen samples. Many of DEGs in trachea tissues participate in a variety of cellular activities, especially cellular metabolism. Immune-related DEGs were mainly enriched at d 3, and 5 postinfection in trachea tissues. While, DEGs in spleen tissues were significantly and mainly enriched into immune-related pathways. The results of this study show the direct interactions between MS and the chicken trachea and spleen for the first time. Early dysregulation of tissue-wide gene expression as observed here set the stage for persistent infection of MS.

CALR frameshift mutations in MPN patient-derived iPSCs accelerate maturation of megakaryocytes

Calreticulin (CALR) mutations are driver mutations in myeloproliferative neoplasms (MPNs), leading to activation of the thrombopoietin receptor and causing abnormal megakaryopoiesis. Here, we generated patient-derived CALRins5- or CALRdel52-positive induced pluripotent stem cells (iPSCs) to establish an MPN disease model for molecular and mechanistic studies. We demonstrated myeloperoxidase deficiency in granulocytic cells derived from homozygous CALR mutant iPSCs, rescued by repairing the mutation using CRISPR/Cas9. iPSC-derived megakaryocytes showed characteristics of primary megakaryocytes such as formation of demarcation membrane system and cytoplasmic pro-platelet protrusions. Importantly, CALR mutations led to enhanced megakaryopoiesis and accelerated megakaryocytic development in a thrombopoietin-independent manner. Mechanistically, our study identified differentially regulated pathways in mutated versus unmutated megakaryocytes, such as hypoxia signaling, which represents a potential target for therapeutic intervention. Altogether, we demonstrate key aspects of mutated CALR-driven pathogenesis dependent on its zygosity, and found novel therapeutic targets, making our model a valuable tool for clinical drug screening in MPNs.

•

CALR-mutated iPSCs allow efficient modeling of human MPN disease

•

CRISPR-mediated repair of CALR mutations rescues normal iPSC function

•

Megakaryopoiesis in CALR-mutated iPSCs is hyperplastic and accelerated

•

Transcriptome screen of mutated megakaryocytes identifies novel therapeutic options

Transcriptional profiling of paediatric ependymomas identifies prognostically significant groups

The majority of supratentorial ependymomas in children contain oncogenic fusions, such as ZFTA-RELA or YAP1-MAMLD1. In contrast, posterior fossa (PF) ependymomas lack recurrent somatic mutations and are classified based on gene expression or methylation profiling into group A (PFA) and group B (PFB). We have applied a novel method, NanoString nCounter Technology, to identify four molecular groups among 16 supratentorial and 50 PF paediatric ependymomas, using 4-5 group-specific signature genes. Clustering analysis of 16 supratentorial ependymomas revealed 9 tumours with a RELA fusion-positive signature (RELA+), 1 tumour with a YAP1 fusion-positive signature (YAP1+), and 6 not-classified tumours. Additionally, we identified one RELA+ tumour among historically diagnosed CNS primitive neuroectodermal tumour samples. Overall, 9 of 10 tumours with the RELA+ signature possessed the ZFTA-RELA fusion as detected by next-generation sequencing (p = 0.005). Similarly, the only tumour with a YAP1+ signature exhibited the YAP1-MAMLD1 fusion. Among the remaining unclassified ependymomas, which did not exhibit the ZFTA-RELA fusion, the ZFTA-MAML2 fusion was detected in one case. Notably, among nine ependymoma patients with the RELA+ signature, eight survived at least 5 years after diagnosis. Clustering analysis of PF tumours revealed 42 samples with PFA signatures and 7 samples with PFB signatures. Clinical characteristics of patients with PFA and PFB ependymomas corroborated the previous findings. In conclusion, we confirm here that the NanoString method is a useful single tool for the diagnosis of all four main molecular groups of ependymoma. The differences in reported survival rates warrant further clinical investigation of patients with the ZFTA-RELA fusion.

Aspects of the Neurospora crassa Sulfur Starvation Response Are Revealed by Transcriptional Profiling and DNA Affinity Purification Sequencing

Accurate nutrient sensing is important for rapid fungal growth and exploitation of available resources. Sulfur is an important nutrient source found in a number of biological macromolecules, including proteins and lipids. The model filamentous fungus Neurospora crassa is capable of utilizing sulfur found in a variety of sources from amino acids to sulfate. During sulfur starvation, the transcription factor CYS-3 is responsible for upregulation of genes involved in sulfur uptake and assimilation. Using a combination of RNA sequencing and DNA affinity purification sequencing, we performed a global survey of the N. crassa sulfur starvation response and the role of CYS-3 in regulating sulfur-responsive genes. The CYS-3 transcription factor bound the promoters and regulated genes involved in sulfur metabolism. Additionally, CYS-3 directly activated the expression of a number of uncharacterized transporter genes, suggesting that regulation of sulfur import is an important aspect of regulation by CYS-3. CYS-3 also directly regulated the expression of genes involved in mitochondrial electron transfer. During sulfur starvation, genes involved in nitrogen metabolism, such as amino acid and nucleic acid metabolic pathways, along with genes encoding proteases and nucleases that are necessary for scavenging nitrogen, were activated. Sulfur starvation also caused changes in the expression of genes involved in carbohydrate metabolism, such as those encoding glycosyl hydrolases. Thus, our data suggest a connection between sulfur metabolism and other aspects of cellular metabolism. IMPORTANCE Identification of nutrients present in the environment is a challenge common to all organisms. Sulfur is an important nutrient source found in proteins, lipids, and electron carriers that are required for the survival of filamentous fungi such as Neurospora crassa. Here, we transcriptionally profiled the response of N. crassa to characterize the global response to sulfur starvation. We also used DNA affinity purification sequencing to identify the direct downstream targets of the transcription factor responsible for regulating genes involved in sulfur uptake and assimilation. Along with genes involved in sulfur metabolism, this transcription factor regulated a number of uncharacterized transporter genes and genes involved in mitochondrial electron transfer. Our data also suggest a connection between sulfur, nitrogen, and carbon metabolism, indicating that the regulation of a number of metabolic pathways is intertwined.

Differential Expression of Immune Response Genes in Asymptomatic Chronic Chagas Disease Patients Versus Healthy Subjects

Infection by the Trypanosoma cruzi parasite causes Chagas disease and triggers multiple immune mechanisms in the host to combat the pathogen. Chagas disease has a variable clinical presentation and progression, producing in the chronic phase a fragile balance between the host immune response and parasite replication that keeps patients in a clinically silent asymptomatic stage for years. Since the parasite is intracellular and replicates within cells, the cell-mediated response of the host adaptive immunity plays a critical role. This function is mainly orchestrated by T lymphocytes, which recognize parasite antigens and promote specific functions to control the infection. However, little is known about the immunological markers associated with this asymptomatic stage of the disease. In this large-scale analysis, the differential expression of 106 immune system-related genes has been analyzed using high-throughput qPCR in T. cruzi antigen-stimulated PBMC from chronic Chagas disease patients with indeterminate form (IND) and healthy donors (HD) from endemic and non-endemic areas of Chagas disease. This analysis revealed that there were no differences in the expression level of most genes under study between healthy donors from endemic and non-endemic areas determined by PCA and differential gene expression analysis. Instead, PCA revealed the existence of different expression profiles between IND patients and HD (p < 0.0001), dependent on the 32 genes included in PC1. Differential gene expression analysis also revealed 23 upregulated genes (expression fold change > 2) and 11 downregulated genes (expression fold change < 0.5) in IND patients versus HD. Enrichment analysis showed that several upregulated genes in IND patients participate in relevant immunological pathways such as antigen-dependent B cell activation, stress induction of HSP regulation, NO2-dependent IL12 pathway in NK cells, and cytokine-inflammatory response. The antigen-specific differential gene expression profile detected in these patients and the relevant immunological pathways that seem to be activated could represent potential biomarkers of the asymptomatic form of Chagas disease, helpful to diagnosis and infection control.

Developing Functional Genomics Platforms for Fungi

Fungi are responsible for diseases that result in the deaths of over a million individuals each year and devastating crop infestations that threaten global food supplies. However, outside of a select few model organisms, the majority of fungal genes are uncharacterized. The roles of these genes in the biology of the organism, pathogenesis, and mediating interactions with the environment and other microbes are unknown. Historically, fungal gene characterization has primarily relied on classical genetic screens. However, advances in sequencing technology have enabled more rapid methods of gene functional characterization. Large-scale transcriptional profiling projects are one solution to generating hypotheses about fungal gene function. Together with other 'omics techniques and newer tools that enable massively parallel mutant screens, knowledge of fungal gene function will be substantially improved. Understanding the function of fungal genes will be instrumental in increasing global food security, protecting ecosystems, and improving health outcomes.

Weighted Gene Co-Expression Network Analysis Combined with Machine Learning Validation to Identify Key Modules and Hub Genes Associated with SARS-CoV-2 Infection

The coronavirus disease-2019 (COVID-19) pandemic has caused an enormous loss of lives. Various clinical trials of vaccines and drugs are being conducted worldwide; nevertheless, as of today, no effective drug exists for COVID-19. The identification of key genes and pathways in this disease may lead to finding potential drug targets and biomarkers. Here, we applied weighted gene co-expression network analysis and LIME as an explainable artificial intelligence algorithm to comprehensively characterize transcriptional changes in bronchial epithelium cells (primary human lung epithelium (NHBE) and transformed lung alveolar (A549) cells) during severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. Our study detected a network that significantly correlated to the pathogenicity of COVID-19 infection based on identified hub genes in each cell line separately. The novel hub gene signature that was detected in our study, including PGLYRP4 and HEPHL1, may shed light on the pathogenesis of COVID-19, holding promise for future prognostic and therapeutic approaches. The enrichment analysis of hub genes showed that the most relevant biological process and KEGG pathways were the type I interferon signaling pathway, IL-17 signaling pathway, cytokine-mediated signaling pathway, and defense response to virus categories, all of which play significant roles in restricting viral infection. Moreover, according to the drug-target network, we identified 17 novel FDA-approved candidate drugs, which could potentially be used to treat COVID-19 patients through the regulation of four hub genes of the co-expression network. In conclusion, the aforementioned hub genes might play potential roles in translational medicine and might become promising therapeutic targets. Further in vitro and in vivo experimental studies are needed to evaluate the role of these hub genes in COVID-19.

Data-Driven Kidney Transplant Phenotyping as a Histology-Independent Framework for Biomarker Discovery

BACKGROUND

In transplant medicine, clinical decision making largely relies on histology of biopsy specimens. However, histology suffers from low specificity, sensitivity, and reproducibility, leading to suboptimal stratification of patients. We developed a histology-independent immune framework of kidney graft homeostasis and rejection.

METHODS

We applied tailored RNA deconvolution for leukocyte enumeration and coregulated gene network analysis to published bulk human kidney transplant RNA transcriptomes as input for unsupervised, high-dimensional phenotype clustering. We used framework-based graft survival analysis to identify a biomarker that was subsequently characterized in independent transplant biopsy specimens.

RESULTS

We found seven immune phenotypes that confirm known rejection types and uncovered novel signatures. The molecular phenotypes allow for improved graft survival analysis compared with histology, and identify a high-risk group in nonrejecting transplants. Two fibrosis-related phenotypes with distinct immune features emerged with reduced graft survival. We identified lysyl oxidase-like 2 (LOXL2)-expressing peritubular CD68+ macrophages as a framework-derived biomarker of impaired allograft function. These cells precede graft fibrosis, as demonstrated in longitudinal biopsy specimens, and may be clinically useful as a biomarker for early fibrogenesis.

CONCLUSIONS

This study provides a comprehensive, data-driven atlas of human kidney transplant phenotypes and demonstrates its utility to identify novel clinical biomarkers.

Drug repurposing strategies of relevance for Parkinson's disease

Parkinson's disease is a highly disabling, progressive neurodegenerative disease that manifests as a mix of motor and non-motor signs. Although we are equipped with some symptomatic treatments, especially for the motor signs of the disease, there are still no established disease-modifying drugs so the disease progresses unchecked. Standard drug discovery programs for disease-modifying therapies have provided key insights into the pathogenesis of Parkinson's disease but, of the many positive candidates identified in pre-clinical studies, none has yet translated into a successful clinically efficacious drug. Given the huge cost of drug discovery programs, it is not surprising that much attention has turned toward repurposing strategies. The trialing of an established therapeutic has the advantage of bypassing the need for preclinical safety testing and formulation optimization, thereby cutting both time and costs involved in getting a treatment to the clinic. Additional reduced failure rates for repurposed drugs are also a potential bonus. Many different strategies for drug repurposing are open to researchers in the Parkinson's disease field. Some of these have already proven effective in identifying suitable drugs for clinical trials, lending support to such approaches. In this review, we present a summary of the different strategies for drug repurposing, from large-scale epidemiological correlation analysis through to single-gene transcriptional approaches. We provide examples of past or ongoing studies adopting each strategy, where these exist. For strategies that have yet to be applied to Parkinson's disease, their utility is illustrated using examples taken from other disorders.

Current Methodological Challenges of Single-Cell and Single-Nucleus RNA-Sequencing in Glomerular Diseases

Single-cell RNA sequencing (scRNA-seq) and single-nucleus RNA-seq (snRNA-seq) allow transcriptomic profiling of thousands of cells from a renal biopsy specimen at a single-cell resolution. Both methods are promising tools to unravel the underlying pathophysiology of glomerular diseases. This review provides an overview of the technical challenges that should be addressed when designing single-cell transcriptomics experiments that focus on glomerulopathies. The isolation of glomerular cells from core needle biopsy specimens for single-cell transcriptomics remains difficult and depends upon five major factors. First, core needle biopsies generate little tissue material, and several samples are required to identify glomerular cells. Second, both fresh and frozen tissue samples may yield glomerular cells, although every experimental pipeline has different (dis)advantages. Third, enrichment for glomerular cells in human tissue before single-cell analysis is challenging because no effective standardized pipelines are available. Fourth, the current warm cell-dissociation protocols may damage glomerular cells and induce transcriptional artifacts, which can be minimized by using cold dissociation techniques at the cost of less efficient cell dissociation. Finally, snRNA-seq methods may be superior to scRNA-seq in isolating glomerular cells; however, the efficacy of snRNA-seq on core needle biopsy specimens remains to be proven. The field of single-cell omics is rapidly evolving, and the integration of these techniques in multiomics assays will undoubtedly create new insights in the complex pathophysiology of glomerular diseases.

Early Transcriptional Liver Signatures in Experimental Visceral Leishmaniasis

Transcriptional analysis of complex biological scenarios has been used extensively, even though sometimes the results of such analysis may prove imprecise or difficult to interpret due to an overwhelming amount of information. In this study, a large-scale real-time qPCR experiment was coupled to multivariate statistical analysis in order to describe the main immunological events underlying the early L. infantum infection in livers of BALB/c mice. High-throughput qPCR was used to evaluate the expression of 223 genes related to immunological response and metabolism 1, 3, 5, and 10 days post infection. This integrative analysis showed strikingly different gene signatures at 1 and 10 days post infection, revealing the progression of infection in the experimental model based on the upregulation of particular immunological response patterns and mediators. The gene signature 1 day post infection was not only characterized by the upregulation of mediators involved in interferon signaling and cell chemotaxis, but also the upregulation of some inhibitory markers. In contrast, at 10 days post infection, the upregulation of many inflammatory and Th1 markers characterized a more defined gene signature with the upregulation of mediators in the IL-12 signaling pathway. Our results reveal a significant connection between the expression of innate immune response and metabolic and inhibitory markers in early L. infantum infection of the liver.

A 2-Dose AERAS-402 Regimen Boosts CD8+ Polyfunctionality in HIV-Negative, BCG-Vaccinated Recipients

Despite the widespread use of BCG, tuberculosis (TB) remains a global threat. Existing vaccine candidates in clinical trials are designed to replace or boost BCG which does not provide satisfying long-term protection. AERAS-402 is a replication-deficient Ad35 vaccine encoding a fusion protein of the M. tuberculosis (Mtb) antigens 85A, 85B, and TB10.4. The present phase I trial assessed the safety and immunogenicity of AERAS-402 in participants living in India – a highly TB-endemic area. Healthy male participants aged 18–45 years with a negative QuantiFERON-TB Gold in-tube test (QFT) were recruited. Enrolled participants (n=12) were randomized 2:1 to receive two intramuscular injections of either AERAS-402 (3 x 10¹⁰ viral particles [vp]); (n=8) or placebo (n=4) on study days 0 and 28. Safety and immunogenicity parameters were evaluated for up to 182 days post the second injection. Immunogenicity was assessed by a flow cytometry-based intracellular cytokine staining (ICS) assay and transcriptional profiling. The latter was examined using dual-color-Reverse-Transcriptase-Multiplex-Ligation-dependent-Probe-Amplification (dc-RT MLPA) assay. AERAS-402 was well tolerated, and no vaccine-related serious adverse events were recorded. The vaccine-induced CD8⁺ T-cell responses were dominated by cells co-expressing IFN-γ, TNF-α, and IL-2 (“polyfunctional” cells) and were more robust than CD4⁺ T-cell responses. Five genes (CXCL10, GNLY, IFI35, IL1B and PTPRCv2) were differentially expressed between the AERAS-402-group and the placebo group, suggesting vaccine-induced responses. Further, compared to pre-vaccination, three genes (CLEC7A, PTPRCv1 and TAGAP) were consistently up-regulated following two doses of vaccination in the AERAS-402-group. No safety concerns were observed for AERAS-402 in healthy Indian adult males. The vaccine-induced predominantly polyfunctional CD8⁺ T cells in response to Ag85B, humoral immunity, and altered gene expression profiles in peripheral blood mononuclear cells (PBMCs) indicative of activation of various immunologically relevant biological pathways.

Transcriptional Profiling of Aflatoxin B1-Induced Oxidative Stress and Inflammatory Response in Macrophages

Aflatoxin B1 (AFB1) is a highly toxic mycotoxin that causes severe suppression of the immune system of humans and animals, as well as enhances reactive oxygen species (ROS) formation, causing oxidative damage. However, the mechanisms underlying the ROS formation and immunotoxicity of AFB1 are poorly understood. This study used the mouse macrophage RAW264.7 cell line and whole-transcriptome sequencing (RNA-Seq) technology to address this knowledge-gap. The results show that AFB1 induced the decrease of cell viability in a dose- and time-dependent manner. AFB1 also significantly increased intracellular productions of ROS and malondialdehyde and decreased glutathione levels. These changes correlated with increased mRNA expression of NOS2, TNF-α and CXCL2 and decreased expression of CD86. In total, 783 differentially expressed genes (DEGs) were identified via RNA-Seq technology. KEGG analysis of the oxidative phosphorylation pathway revealed that mRNA levels of ND1, ND2, ND3, ND4, ND4L, ND5, ND6, Cyt b, COX2, ATPeF0A and ATPeF08 were higher in AFB1-treated cells than control cells, whereas 14 DEGs were downregulated in the AFB1 group. Furthermore, seven immune regulatory pathways mediated by oxidative stress were identified by KEGG analysis. Altogether, these data suggest that AFB1 induces oxidative stress in macrophages via affecting the respiratory chain, which leads to the activation of several signaling pathways related to the inflammatory response.

Variable Outcomes in Neural Differentiation of Human PSCs Arise from Intrinsic Differences in Developmental Signaling Pathways

Directed differentiation of human pluripotent stem cells varies in specificity and efficiency. Stochastic, genetic, intracellular, and environmental factors affect maintenance of pluripotency and differentiation into early embryonic lineages. However, factors affecting variation in in vitro differentiation to defined cell types are not well understood. To address this, we focused on a well-established differentiation process to cerebral cortex neural progenitor cells and their neuronal progeny from human pluripotent stem cells. Analysis of 162 differentiation outcomes of 61 stem cell lines derived from 37 individuals showed that most variation occurs along gene expression axes reflecting dorsoventral and rostrocaudal spatial expression during in vivo brain development. Line-independent and line-dependent variations occur, with the latter driven largely by differences in endogenous Wnt signaling activity. Tuning Wnt signaling during a specific phase early in the differentiation process reduces variability, demonstrating that cell-line/genome-specific differentiation outcome biases can be corrected by controlling extracellular signaling.

•

Analysis of 162 iPSC neural differentiation quantifies sources of variability

•

Variation in outcomes occurs along developmental brain spatial and regional axes

•

Variation depends on iPSC-line-specific differences in Wnt/β-catenin signaling

•

Effects of signaling differences can be rescued by exogenous pathway activation

Dynamics of the transcriptional landscape during human fetal testis and ovary development

STUDY QUESTION

Which transcriptional program triggers sex differentiation in bipotential gonads and downstream cellular events governing fetal testis and ovary development in humans?

SUMMARY ANSWER

The characterization of a dynamically regulated protein-coding and non-coding transcriptional landscape in developing human gonads of both sexes highlights a large number of potential key regulators that show an early sexually dimorphic expression pattern.

WHAT IS KNOWN ALREADY

Gonadal sex differentiation is orchestrated by a sexually dimorphic gene expression program in XX and XY developing fetal gonads. A comprehensive characterization of its non-coding counterpart offers promising perspectives for deciphering the molecular events underpinning gonad development and for a complete understanding of the etiology of disorders of sex development in humans.

STUDY DESIGN, SIZE, DURATION

To further investigate the protein-coding and non-coding transcriptional landscape during gonad differentiation, we used RNA-sequencing (RNA-seq) and characterized the RNA content of human fetal testis (N = 24) and ovaries (N = 24) from 6 to 17 postconceptional week (PCW), a key period in sex determination and gonad development.

PARTICIPANTS/MATERIALS, SETTING, METHODS

First trimester fetuses (6-12 PCW) and second trimester fetuses (13-14 and 17 PCW) were obtained from legally induced normally progressing terminations of pregnancy. Total RNA was extracted from whole human fetal gonads and sequenced as paired-end 2 × 50 base reads. Resulting sequences were mapped to the human genome, allowing for the assembly and quantification of corresponding transcripts.

MAIN RESULTS AND THE ROLE OF CHANCE

This RNA-seq analysis of human fetal testes and ovaries at seven key developmental stages led to the reconstruction of 22 080 transcripts differentially expressed during testicular and/or ovarian development. In addition to 8935 transcripts displaying sex-independent differential expression during gonad development, the comparison of testes and ovaries enabled the discrimination of 13 145 transcripts that show a sexually dimorphic expression profile. The latter include 1479 transcripts differentially expressed as early as 6 PCW, including 39 transcription factors, 40 long non-coding RNAs and 20 novel genes. Despite the use of stringent filtration criteria (expression cut-off of at least 1 fragment per kilobase of exon model per million reads mapped, fold change of at least 2 and false discovery rate adjusted P values of less than <1%), the possibility of assembly artifacts and of false-positive differentially expressed transcripts cannot be fully ruled out.

LARGE-SCALE DATA

Raw data files (fastq) and a searchable table (.xlss) containing information on genomic features and expression data for all refined transcripts have been submitted to the NCBI GEO under accession number GSE116278.

LIMITATIONS, REASONS FOR CAUTION

The intrinsic nature of this bulk analysis, i.e. the sequencing of transcripts from whole gonads, does not allow direct identification of the cellular origin(s) of the transcripts characterized. Potential cellular dilution effects (e.g. as a result of distinct proliferation rates in XX and XY gonads) may account for a few of the expression profiles identified as being sexually dimorphic. Finally, transcriptome alterations that would result from exposure to pre-abortive drugs cannot be completely excluded. Although we demonstrated the high quality of the sorted cell populations used for experimental validations using quantitative RT-PCR, it cannot be totally excluded that some germline expression may correspond to cell contamination by, for example, macrophages.

WIDER IMPLICATIONS OF THE FINDINGS

For the first time, this study has led to the identification of 1000 protein-coding and non-coding candidate genes showing an early, sexually dimorphic, expression pattern that have not previously been associated with sex differentiation. Collectively, these results increase our understanding of gonad development in humans, and contribute significantly to the identification of new candidate genes involved in fetal gonad differentiation. The results also provide a unique resource that may improve our understanding of the fetal origin of testicular and ovarian dysgenesis syndromes, including cryptorchidism and testicular cancers.

STUDY FUNDING/COMPETING INTEREST(S)

This work was supported by the French National Institute of Health and Medical Research (Inserm), the University of Rennes 1, the French School of Public Health (EHESP), the Swiss National Science Foundation [SNF n° CRS115_171007 to B.J.], the French National Research Agency [ANR n° 16-CE14-0017-02 and n° 18-CE14-0038-02 to F.C.], the Medical Research Council [MR/L010011/1 to P.A.F.] and the European Community's Seventh Framework Programme (FP7/2007-2013) [under grant agreement no 212885 to P.A.F.] and from the European Union's Horizon 2020 Research and Innovation Programme [under grant agreement no 825100 to P.A.F. and S.M.G.]. There are no competing interests related to this study.

A subpopulation of embryonic microglia respond to maternal stress and influence nearby neural progenitors

The interplay between hypothalamic neurons and microglia as they integrate stressors to regulate homeostasis is of growing interest. We asked if microglia in the embryonic hypothalamus were likewise stress responsive and, if so, whether their precocious activation perturbs nearby neural stem cell (NSC) programs. We performed single-cell transcriptomics to define embryonic hypothalamic microglia heterogeneity and identified four microglial subsets, including a subpopulation adjacent to NSCs that was responsive to gestational cold stress. Stress exposure elevated CCL3 and CCL4 secretion, but only in male brains, and ex vivo CCL4 treatment of hypothalamic NSCs altered proliferation and differentiation. Concomitantly, gestational stress decreased PVN oxytocin neurons only in male embryos, which was reversed by microglia depletion. Adult offspring exposed to gestational stress displayed altered social behaviors, which was likewise microglia dependent, but only in males. Collectively, immature hypothalamic microglia play an unappreciated role in translating maternal stressors to sexually dimorphic perturbation of neurodevelopmental programs.

Transcriptome Signatures Predict Phenotypic Variations of Candida auris

Health care facilities are facing serious threats by the recently emerging human fungal pathogen Candida auris owing to its pronounced antifungal multidrug resistance and poor diagnostic tools. Distinct C. auris clades evolved seemingly simultaneously at independent geographical locations and display both genetic and phenotypic diversity. Although comparative genomics and phenotypic profiling studies are increasing, we still lack mechanistic knowledge about the C. auris species diversification and clinical heterogeneity. Since gene expression variability impacts phenotypic plasticity, we aimed to characterize transcriptomic signatures of C. auris patient isolates with distinct antifungal susceptibility profiles in this study. First, we employed an antifungal susceptibility screening of clinical C. auris isolates to identify divergent intra-clade responses to antifungal treatments. Interestingly, comparative transcriptional profiling reveals large gene expression differences between clade I isolates and one clade II strain, irrespective of their antifungal susceptibilities. However, comparisons at the clade levels demonstrate that minor changes in gene expression suffice to drive divergent drug responses. Finally, we functionally validate transcriptional signatures reflecting phenotypic divergence of clinical isolates. Thus, our results suggest that large-scale transcriptional profiling allows for predicting phenotypic diversities of patient isolates, which may help choosing suitable antifungal therapies of multidrug-resistant C. auris.

Pericyte-like cells undergo transcriptional reprogramming and distinct functional adaptations in acute lung injury

We previously reported on the role of pericyte-like cells as functional sentinel immune cells in lung injury. However, much about the biological role of pericytes in lung injury remains unknown. Lung pericyte-like cells are well-positioned to sense disruption to the epithelial barrier and coordinate local inflammatory responses due to their anatomic niche within the alveoli. In this report, we characterized transcriptional responses and functional changes in pericyte-like cells following activation by alveolar components from injured and uninjured lungs in a mouse model of acute lung injury (ALI). Purified pericyte-like cells from lung digests using PDGFRβ as a selection marker were expanded in culture as previously described (1). We induced sterile acute lung injury in mice with recombinant human Fas ligand (rhFasL) instillation followed by mechanical ventilation (1). We then collected bronchoalveolar lavage fluid (BALF) from injured and uninjured mice. Purified pericyte-like cells in culture were exposed to growth media only (control), BALF from uninjured mice, and BALF from injured mice for 6 and 24 hours. RNA collected from these treatment conditions were processed for RNAseq. Targets of interest identified by pathway analysis were validated using in vitro and in vivo assays. We observed robust global transcriptional changes in pericyte-like cells following treatment with uninjured and injured BALF at 6 hours, but this response persisted for 24 hours only after exposure to injured BALF. Functional enrichment analysis of pericytes treated with injured BALF revealed the activation of pro-inflammatory, cell migration, and angiogenesis-related pathways, whereas processes associated with tissue development and cell differentiation were down-regulated. We validated select upregulated targets in the inflammatory, angiogenic, and cell migratory pathways using functional biological assays in vitro and in vivo. We conclude that lung pericyte-like cells are highly responsive to alveolar compartment content from both uninjured and injured lungs, but injured BALF elicits a more sustained response. The inflammatory, angiogenic, and migratory changes exhibited by activated pericyte-like cells underscore the phenotypic plasticity of these specialized stromal cells in the setting of acute lung injury.

Cyclin-Dependent Kinase 1 Activity Is a Driver of Cyst Growth in Polycystic Kidney Disease

BACKGROUND

Mutations in PKD1 and PKD2, which encode the transmembrane proteins polycystin-1 and polycystin-2, respectively, cause autosomal dominant polycystic kidney disease (ADPKD). Polycystins are expressed in the primary cilium, and disrupting cilia structure significantly slows ADPKD progression following inactivation of polycystins. The cellular mechanisms of polycystin- and cilia-dependent cyst progression in ADPKD remain incompletely understood.

METHODS

Unbiased transcriptional profiling in an adult-onset Pkd2 mouse model before cysts formed revealed significant differentially expressed genes (DEGs) in Pkd2 single-knockout kidneys, which were used to identify candidate pathways dysregulated in kidneys destined to form cysts. In vivo studies validated the role of the candidate pathway in the progression of ADPKD. Wild-type and Pkd2/Ift88 double-knockout mice that are protected from cyst growth served as controls.

RESULTS

The RNASeq data identified cell proliferation as the most dysregulated pathway, with 15 of 241 DEGs related to cell cycle functions. Cdk1 appeared as a central component in this analysis. Cdk1 expression was similarly dysregulated in Pkd1 models of ADPKD, and conditional inactivation of Cdk1 with Pkd1 markedly improved the cystic phenotype and kidney function compared with inactivation of Pkd1 alone. The Pkd1/Cdk1 double knockout blocked cyst cell proliferation that otherwise accompanied Pkd1 inactivation alone.

CONCLUSIONS

Dysregulation of Cdk1 is an early driver of cyst cell proliferation in ADPKD due to Pkd1 inactivation. Selective targeting of cyst cell proliferation is an effective means of slowing ADPKD progression caused by inactivation of Pkd1.

Genome-Wide Identification of Wheat WRKY Gene Family Reveals That TaWRKY75-A Is Referred to Drought and Salt Resistances

It is well known that WRKY transcription factors play essential roles in plants' response to diverse stress responses, especially to drought and salt stresses. However, a full comprehensive analysis of this family in wheat is still missing. Here we used in silico analysis and identified 124 WRKY genes, including 294 homeologous copies from a high-quality reference genome of wheat (Triticum aestivum). We also found that the TaWRKY gene family did not undergo gene duplication rather than gene loss during the evolutionary process. The TaWRKY family members displayed different expression profiles under several abiotic stresses, indicating their unique functions in the mediation of particular responses. Furthermore, TaWRKY75-A was highly induced after polyethylene glycol and salt treatments. The ectopic expression of TaWRKY75-A in Arabidopsis enhanced drought and salt tolerance. A comparative transcriptome analysis demonstrated that TaWRKY75-A integrated jasmonic acid biosynthetic pathway and other potential metabolic pathways to increase drought and salt resistances in transgenic Arabidopsis. Our study provides valuable insights into the WRKY family in wheat and will generate a useful genetic resource for improving wheat breeding.

Transcriptional Profiling of Immune and Inflammatory Responses in the Context of SARS-CoV-2 Fungal Superinfection in a Human Airway Epithelial Model

An increasing amount of evidence indicates a relatively high prevalence of superinfections associated with coronavirus disease 2019 (COVID-19), including invasive aspergillosis, but the underlying mechanisms remain to be characterized. In the present study, to better understand the biological impact of superinfection, we determine and compare the host transcriptional response to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) versus Aspergillus superinfection, using a model of reconstituted human airway epithelium. Our analyses reveal that both simple infection and superinfection induce strong deregulation of core components of innate immune and inflammatory responses, with a stronger response to superinfection in the bronchial epithelial model compared to its nasal counterpart. Our results also highlight unique transcriptional footprints of SARS-CoV-2 Aspergillus superinfection, such as an imbalanced type I/type III IFN, and an induction of several monocyte and neutrophil associated chemokines, that could be useful for the understanding of Aspergillus-associated COVID-19 and also the management of severe forms of aspergillosis in this specific context.

A Distinct Transcriptional Program in Human CAR T Cells Bearing the 4-1BB Signaling Domain Revealed by scRNA-Seq

T cells engineered to express chimeric antigen receptors (CARs) targeting CD19 have produced impressive outcomes for the treatment of B cell malignancies, but different products vary in kinetics, persistence, and toxicity profiles based on the co-stimulatory domains included in the CAR. In this study, we performed transcriptional profiling of bulk CAR T cell populations and single cells to characterize the transcriptional states of human T cells transduced with CD3ζ, 4-1BB-CD3ζ (BBζ), or CD28-CD3ζ (28ζ) co-stimulatory domains at rest and after activation by triggering their CAR or their endogenous T cell receptor (TCR). We identified a transcriptional signature common across CARs with the CD3ζ signaling domain, as well as a distinct program associated with the 4-1BB co-stimulatory domain at rest and after activation. CAR T cells bearing BBζ had increased expression of human leukocyte antigen (HLA) class II genes, ENPP2, and interleukin (IL)-21 axis genes, and decreased PD1 compared to 28ζ CAR T cells. Similar to previous studies, we also found BBζ CAR CD8 T cells to be enriched in a central memory cell phenotype and fatty acid metabolism genes. Our data uncovered transcriptional signatures related to costimulatory domains and demonstrated that signaling domains included in CARs uniquely shape the transcriptional programs of T cells.

Genetics, Transcriptomics and Meta-Taxonomics in Visceral Leishmaniasis

Visceral leishmaniasis (VL) caused by parasites of the Leishmania donovani complex can be fatal in susceptible individuals. Understanding the interactions between host and pathogen is one way to obtain leads to develop better drugs and for vaccine development. In recent years multiple omics-based approaches have assisted researchers to gain a more global picture of this interaction in leishmaniasis. Here we review results from studies using three omics-based approaches to study VL caused by L. donovani in India: (i) chip-based analysis of single nucleotide variants in the first genome-wide association study of host genetic risk factors for VL, followed by analysis of epitope binding to HLA DRB1 risk versus protective alleles; (ii) transcriptional profiling demonstrating pathways important in Amphotericin B treated compared to active VL cases, including demonstration that anti-interleukin-10 unleashes a storm of chemokines and cytokines in whole blood responses to soluble leishmania antigen in active cases; and (iii) a meta-taxonomic approach based on sequencing amplicons derived from regions of 16S ribosomal RNA (16S rRNA) and 18S rRNA genes that allowed us to determine composition of both prokaryotic and eukaryotic gut microflora in VL cases compared to endemic controls. Overall, our omics-based approaches demonstrate that global analyses of genetic risk factors, host responses to infection, and the interaction between host, parasite and the microbiome can point to the most critical factors that determine the outcome of infection.

Coordinated Circulating T Follicular Helper and Activated B Cell Responses Underlie the Onset of Antibody-Mediated Rejection in Kidney Transplantation

BACKGROUND

Although antibody-mediated rejection (ABMR) has been long recognized as a leading cause of allograft failure after kidney transplantation, the cellular and molecular processes underlying the induction of deleterious donor-specific antibody (DSA) responses remain poorly understood.

METHODS

Using high-dimensional flow cytometry, in vitro assays, and RNA sequencing, we concomitantly investigated the role of T follicular helper (T_FH) cells and B cells during ABMR in 105 kidney transplant recipients.

RESULTS

There were 54 patients without DSAs; of those with DSAs, ABMR emerged in 20 patients, but not in 31 patients. We identified proliferating populations of circulating T_FH cells and activated B cells emerging in blood of patients undergoing ABMR. Although these circulating T_FH cells comprised heterogeneous phenotypes, they were dominated by activated (ICOS⁺PD-1⁺) and early memory precursor (CCR7⁺CD127⁺) subsets, and were enriched for the transcription factors IRF4 and c-Maf. These circulating T_FH cells produced large amounts of IL-21 upon stimulation with donor antigen and induced B cells to differentiate into antibody-secreting cells that produced DSAs. Combined analysis of the matched circulating T_FH cell and activated B cell RNA-sequencing profiles identified highly coordinated transcriptional programs in circulating T_FH cells and B cells among patients with ABMR, which markedly differed from those of patients who did not develop DSAs or ABMR. The timing of expansion of the distinctive circulating T_FH cells and activated B cells paralleled emergence of DSAs in blood, and their magnitude was predictive of IgG3 DSA generation, more severe allograft injury, and higher rate of allograft loss.

CONCLUSIONS

Patients undergoing ABMR may benefit from monitoring and therapeutic targeting of T_FH cell-B cell interactions.

Harnessing Expressed Single Nucleotide Variation and Single Cell RNA Sequencing To Define Immune Cell Chimerism in the Rejecting Kidney Transplant

BACKGROUND

In solid organ transplantation, donor-derived immune cells are assumed to decline with time after surgery. Whether donor leukocytes persist within kidney transplants or play any role in rejection is unknown, however, in part because of limited techniques for distinguishing recipient from donor cells.

METHODS

Whole-exome sequencing of donor and recipient DNA and single-cell RNA sequencing (scRNA-seq) of five human kidney transplant biopsy cores distinguished immune cell contributions from both participants. DNA-sequence comparisons used single nucleotide variants (SNVs) identified in the exome sequences across all samples.

RESULTS

Analysis of expressed SNVs in the scRNA-seq data set distinguished recipient versus donor origin for all 81,139 cells examined. The leukocyte donor/recipient ratio varied with rejection status for macrophages and with time post-transplant for lymphocytes. Recipient macrophages displayed inflammatory activation whereas donor macrophages demonstrated antigen presentation and complement signaling. Recipient-origin T cells expressed cytotoxic and proinflammatory genes consistent with an effector cell phenotype, whereas donor-origin T cells appeared quiescent, expressing oxidative phosphorylation genes. Finally, both donor and recipient T cell clones within the rejecting kidney suggested lymphoid aggregation. The results indicate that donor-origin macrophages and T cells have distinct transcriptional profiles compared with their recipient counterparts, and that donor macrophages can persist for years post-transplantation.

CONCLUSIONS

Analysis of single nucleotide variants and their expression in single cells provides a powerful novel approach to accurately define leukocyte chimerism in a complex organ such as a transplanted kidney, coupled with the ability to examine transcriptional profiles at single-cell resolution.

Transcriptional Profiling of Normal, Stenotic, and Regurgitant Human Aortic Valves

The genetic mechanisms underlying aortic stenosis (AS) and aortic insufficiency (AI) disease progression remain unclear. We hypothesized that normal aortic valves and those with AS or AI all exhibit unique transcriptional profiles. Normal control (NC) aortic valves were collected from non-matched donor hearts that were otherwise acceptable for transplantation (n = 5). Valves with AS or AI (n = 5, each) were collected from patients undergoing surgical aortic valve replacement. High-throughput sequencing of total RNA revealed 6438 differentially expressed genes (DEGs) for AS vs. NC, 4994 DEGs for AI vs. NC, and 2771 DEGs for AS vs. AI. Among 21 DEGs of interest, APCDD1L, CDH6, COL10A1, HBB, IBSP, KRT14, PLEKHS1, PRSS35, and TDO2 were upregulated in both AS and AI compared to NC, whereas ALDH1L1, EPHB1, GPX3, HIF3A, and KCNT1 were downregulated in both AS and AI (p < 0.05). COL11A1, H19, HIF1A, KCNJ6, PRND, and SPP1 were upregulated only in AS, and NPY was downregulated only in AS (p < 0.05). The functional network for AS clustered around ion regulation, immune regulation, and lipid homeostasis, and that for AI clustered around ERK1/2 regulation. Overall, we report transcriptional profiling data for normal human aortic valves from non-matched donor hearts that were acceptable for transplantation and demonstrated that valves with AS and AI possess unique genetic signatures. These data create a roadmap for the development of novel therapeutics to treat AS and AI.