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.

Kidney organoids: accurate models or fortunate accidents.

There are now many reports of human kidney organoids generated via the directed differentiation of human pluripotent stem cells (PSCs) based on an existing understanding of mammalian kidney organogenesis. Such kidney organoids potentially represent tractable tools for the study of normal human development and disease with improvements in scale, structure, and functional maturation potentially providing future options for renal regeneration. The utility of such organotypic models, however, will ultimately be determined by their developmental accuracy. While initially inferred from mouse models, recent transcriptional analyses of human fetal kidney have provided greater insight into nephrogenesis. In this review, we discuss how well human kidney organoids model the human fetal kidney and how the remaining differences challenge their utility.

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.

Unravelling the genetic basis and regulation networks related to fibre quality improvement using chromosome segment substitution lines in cotton.

The elucidation of genetic architecture and molecular regulatory networks underlying complex traits remains a significant challenge in life science, largely due to the substantial background effects that arise from epistasis and gene-environment interactions. The chromosome segment substitution line (CSSL) is an ideal material for genetic and molecular dissection of complex traits due to its near-isogenic properties; yet a comprehensive analysis, from the basic identification of substitution segments to advanced regulatory network, is still insufficient. Here, we developed two cotton CSSL populations on the Gossypium hirsutum background, representing wide adaptation and high lint yield, with introgression from G. barbadense, representing superior fibre quality. We sequenced 99 CSSLs that demonstrated significant differences from G. hirsutum in fibre, and characterized 836 dynamic fibre transcriptomes in three crucial developmental stages. We developed a workflow for precise resolution of chromosomal substitution segments; the genome sequencing revealed substitutions collectively representing 87.25% of the G. barbadense genome. Together, the genomic and transcriptomic survey identified 18 novel fibre-quality-related quantitative trait loci with high genetic contributions and the comprehensive landscape of fibre development regulation. Furthermore, analysis determined unique cis-expression patterns in CSSLs to be the driving force for fibre quality alteration; building upon this, the co-expression regulatory network revealed biological relationships among the noted pathways and accurately described the molecular interactions of GhHOX3, GhRDL1 and GhEXPA1 during fibre elongation, along with reliable predictions for their interactions with GhTBA8A5. Our study will enhance more strategic employment of CSSL in crop molecular biology and breeding programmes.

Systematic identification and characterization of five transcription factors mediating the oxidative stress response in Candida albicans.

Candida albicans is an opportunistic human fungal pathogen that causes superficial and systemic infections, particularly in immunocompromised individuals. In response to C. albicans infection, innate immune cells of the host produce and accumulate reactive oxygen species (ROS), which can lead to irreversible damage and apoptosis of fungal cells. Several transcription factors involved in this oxidative stress response have been identified; however, a systematic study to identify the transcription factors that mediate the oxidative stress response has not yet been conducted. Here, we screened a comprehensive transcription factor mutant library consisting of 211 transcription factor deletion mutant strains in the presence and absence of hydrogen peroxide (H2O2), a potent ROS inducer, and identified five transcription factors (Skn7, Dpb4, Cap1, Dal81, and Stp2) that are sensitive to H2O2. Genome-wide transcriptional profiling revealed that H2O2 induces a discrete set of differentially regulated genes among the five identified transcription factor mutant strains. Functional enrichment analysis identified KEGG pathways pertaining to glycolysis/gluconeogenesis, amino sugar and nucleotide sugar metabolism, and ribosome synthesis as the most enriched pathways. GO term analysis of the top common differentially expressed genes among the transcription factor mutant strains identified hexose catabolism and iron transport as the most enriched GO terms upon exposure to H2O2. This study is the first to systematically identify and characterise the transcription factors involved in the response to H2O2. Based on our transcriptional profiling results, we found that exposure to H2O2 modulates several downstream genes involved in fungal virulence. Overall, this study sheds new light on the metabolism, physiological functions, and cellular processes involved in the H2O2-induced oxidative stress response in C. albicans.

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.

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.

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.

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.

Ablation of Death-Associated Protein Kinase 1 Changes the Transcriptomic Profile and Alters Neural-Related Pathways in the Brain.

Death-associated protein kinase 1 (DAPK1), a Ca2+/calmodulin-dependent serine/threonine kinase, mediates various neuronal functions, including cell death. Abnormal upregulation of DAPK1 is observed in human patients with neurological diseases, such as Alzheimer's disease (AD) and epilepsy. Ablation of DAPK1 expression and suppression of DAPK1 activity attenuates neuropathology and behavior impairments. However, whether DAPK1 regulates gene expression in the brain, and whether its gene profile is implicated in neuronal disorders, remains elusive. To reveal the function and pathogenic role of DAPK1 in neurological diseases in the brain, differential transcriptional profiling was performed in the brains of DAPK1 knockout (DAPK1-KO) mice compared with those of wild-type (WT) mice by RNA sequencing. We showed significantly altered genes in the cerebral cortex, hippocampus, brain stem, and cerebellum of both male and female DAPK1-KO mice compared to those in WT mice, respectively. The genes are implicated in multiple neural-related pathways, including: AD, Parkinson's disease (PD), Huntington's disease (HD), neurodegeneration, glutamatergic synapse, and GABAergic synapse pathways. Moreover, our findings imply that the potassium voltage-gated channel subfamily A member 1 (Kcna1) may be involved in the modulation of DAPK1 in epilepsy. Our study provides insight into the pathological role of DAPK1 in the regulatory networks in the brain and new therapeutic strategies for the treatment of neurological diseases.

Neurodegeneration and neuroinflammation are linked, but independent of alpha-synuclein inclusions, in a seeding/spreading mouse model of Parkinson's disease.

A key pathological process in Parkinson's disease (PD) is the transneuronal spreading of α-synuclein. Alpha-synuclein (α-syn) is a presynaptic protein that, in PD, forms pathological inclusions. Other hallmarks of PD include neurodegeneration and microgliosis in susceptible brain regions. Whether it is primarily transneuronal spreading of α-syn particles, inclusion formation, or other mechanisms, such as inflammation, that cause neurodegeneration in PD is unclear. We used a model of spreading of α-syn induced by striatal injection of α-syn preformed fibrils into the mouse striatum to address this question. We performed quantitative analysis for α-syn inclusions, neurodegeneration, and microgliosis in different brain regions, and generated gene expression profiles of the ventral midbrain, at two different timepoints after disease induction. We observed significant neurodegeneration and microgliosis in brain regions not only with, but also without α-syn inclusions. We also observed prominent microgliosis in injured brain regions that did not correlate with neurodegeneration nor with inclusion load. Using longitudinal gene expression profiling, we observed early gene expression changes, linked to neuroinflammation, that preceded neurodegeneration, indicating an active role of microglia in this process. Altered gene pathways overlapped with those typical of PD. Our observations indicate that α-syn inclusion formation is not the major driver in the early phases of PD-like neurodegeneration, but that microglia, activated by diffusible, oligomeric α-syn, may play a key role in this process. Our findings uncover new features of α-syn induced pathologies, in particular microgliosis, and point to the necessity for a broader view of the process of α-syn spreading.

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.

Transcriptome analysis of mouse and human sinoatrial node cells reveals a conserved genetic program.

The rate of contraction of the heart relies on proper development and function of the sinoatrial node, which consists of a small heterogeneous cell population, including Tbx3+ pacemaker cells. Here, we have isolated and characterized the Tbx3+ cells from Tbx3+/Venus knock-in mice. We studied electrophysiological parameters during development and found that Venus-labeled cells are genuine Tbx3+ pacemaker cells. We analyzed the transcriptomes of late fetal FACS-purified Tbx3+ sinoatrial nodal cells and Nppb-Katushka+ atrial and ventricular chamber cardiomyocytes, and identified a sinoatrial node-enriched gene program, including key nodal transcription factors, BMP signaling and Smoc2, the disruption of which in mice did not affect heart rhythm. We also obtained the transcriptomes of the sinoatrial node region, including pacemaker and other cell types, and right atrium of human fetuses, and found a gene program including TBX3, SHOX2, ISL1 and HOX family members, and BMP and NOTCH signaling components conserved between human and mouse. We conclude that a conserved gene program characterizes the sinoatrial node region and that the Tbx3+/Venus allele provides a reliable tool for visualizing the sinoatrial node, and studying its development and function.

Transcriptional and biochemical analyses of Planomicrobium strain AX6 from Qinghai-Tibetan Plateau, China, reveal hydrogen peroxide scavenging potential.

BACKGROUND: The bacterial mechanisms responsible for hydrogen peroxide (H2O2) scavenging have been well-reported, yet little is known about how bacteria isolated from cold-environments respond to H2O2 stress. Therefore, we investigated the transcriptional profiling of the Planomicrobium strain AX6 strain isolated from the cold-desert ecosystem in the Qaidam Basin, Qinghai-Tibet Plateau, China, in response to H2O2 stress aiming to uncover the molecular mechanisms associated with H2O2 scavenging potential. METHODS: We investigated the H2O2-scavenging potential of the bacterial Planomicrobium strain AX6 isolated from the cold-desert ecosystem in the Qaidam Basin, Qinghai-Tibet Plateau, China. Furthermore, we used high-throughput RNA-sequencing to unravel the molecular aspects associated with the H2O2 scavenging potential of the Planomicrobium strain AX6 isolate. RESULTS: In total, 3,427 differentially expressed genes (DEGs) were identified in Planomicrobium strain AX6 isolate in response to 4 h of H2O2 (1.5 mM) exposure. Besides, Kyoto Encyclopedia of Genes and Genomes pathway and Gene Ontology analyses revealed the down- and/or up-regulated pathways following H2O2 treatment. Our study not only identified the H2O2 scavenging capability of the strain nevertheless also a range of mechanisms to cope with the toxic effect of H2O2 through genes involved in oxidative stress response. Compared to control, several genes coding for antioxidant proteins, including glutathione peroxidase (GSH-Px), Coproporphyrinogen III oxidase, and superoxide dismutase (SOD), were relatively up-regulated in Planomicrobium strain AX6, when exposed to H2O2. CONCLUSIONS: Overall, the results suggest that the up-regulated genes responsible for antioxidant defense pathways serve as essential regulatory mechanisms for removing H2O2 in Planomicrobium strain AX6. The DEGs identified here could provide a competitive advantage for the existence of Planomicrobium strain AX6 in H2O2-polluted environments.

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.

Transcriptional Profiling of Pseudomonas aeruginosa Infections.

Pseudomonas aeruginosa is an opportunistic pathogen that causes life-devastating acute as well as chronic biofilm-associated infections with limited treatment options. Its success is largely due to its remarkable adaptability. P. aeruginosa uses different long- and short-term adaptive mechanisms to increase its fitness, both at the population level through genetic diversification and at the individual cell level by adapting gene expression. These adapted gene expression profiles can be fixed by the accumulation of patho-adaptive mutations. The latter are often found in transcriptional regulators and lead to rewiring of the regulatory network to promote survival at the infected host site. In this chapter, we review recent developments in transcriptional profiling and explain how these provide new insights into the establishment and maintenance of P. aeruginosa infections. We illustrate what can be learned from the application of advanced RNA-seq technology, such as ex vivo RNA-seq, host-pathogen crosstalk (dual RNA-seq), or recording of transcriptional heterogeneity within a bacterial population (single-cell RNA-seq). In addition, we discuss how large transcriptome datasets from a variety of clinical isolates can be used to gain an expanded understanding of bacterial adaptation during the infection process. Global genotype-phenotype correlation studies provide a unique opportunity to discover new evolutionary pathways of infection-related phenotypes and led to the discovery of different strategies of the pathogen P. aeruginosa to build a biofilm. Insights gained from large-scale, multi-layered functional -omics approaches will continue to contribute to a more comprehensive understanding of P. aeruginosa adaptation to the host habitat and promises to pave the way for novel strategies to combat recalcitrant infections.

Early sarcomere and metabolic defects in a zebrafish pitx2c cardiac arrhythmia model.

Atrial fibrillation (AF) is the most common type of cardiac arrhythmia. The major AF susceptibility locus 4q25 establishes long-range interactions with the promoter of PITX2, a transcription factor gene with critical functions during cardiac development. While many AF-linked loci have been identified in genome-wide association studies, mechanistic understanding into how genetic variants, including those at the 4q25 locus, increase vulnerability to AF is mostly lacking. Here, we show that loss of pitx2c in zebrafish leads to adult cardiac phenotypes with substantial similarities to pathologies observed in AF patients, including arrhythmia, atrial conduction defects, sarcomere disassembly, and altered cardiac metabolism. These phenotypes are also observed in a subset of pitx2c+/- fish, mimicking the situation in humans. Most notably, the onset of these phenotypes occurs at an early developmental stage. Detailed analyses of pitx2c loss- and gain-of-function embryonic hearts reveal changes in sarcomeric and metabolic gene expression and function that precede the onset of cardiac arrhythmia first observed at larval stages. We further find that antioxidant treatment of pitx2c-/- larvae significantly reduces the incidence and severity of cardiac arrhythmia, suggesting that metabolic dysfunction is an important driver of conduction defects. We propose that these early sarcomere and metabolic defects alter cardiac function and contribute to the electrical instability and structural remodeling observed in adult fish. Overall, these data provide insight into the mechanisms underlying the development and pathophysiology of some cardiac arrhythmias and importantly, increase our understanding of how developmental perturbations can predispose to functional defects in the adult heart.

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.

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.

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.