Kidney Single-cell Transcriptomes Predict Spatial Corticomedullary Gene Expression and Tissue Osmolality Gradients.

BACKGROUND: Single-cell transcriptomes from dissociated tissues provide insights into cell types and their gene expression and may harbor additional information on spatial position and the local microenvironment. The kidney's cells are embedded into a gradient of increasing tissue osmolality from the cortex to the medulla, which may alter their transcriptomes and provide cues for spatial reconstruction. METHODS: Single-cell or single-nuclei mRNA sequencing of dissociated mouse kidneys and of dissected cortex, outer, and inner medulla, to represent the corticomedullary axis, was performed. Computational approaches predicted the spatial ordering of cells along the corticomedullary axis and quantitated expression levels of osmo-responsive genes. In situ hybridization validated computational predictions of spatial gene-expression patterns. The strategy was used to compare single-cell transcriptomes from wild-type mice to those of mice with a collecting duct-specific knockout of the transcription factor grainyhead-like 2 (Grhl2CD-/-), which display reduced renal medullary osmolality. RESULTS: Single-cell transcriptomics from dissociated kidneys provided sufficient information to approximately reconstruct the spatial position of kidney tubule cells and to predict corticomedullary gene expression. Spatial gene expression in the kidney changes gradually and osmo-responsive genes follow the physiologic corticomedullary gradient of tissue osmolality. Single-nuclei transcriptomes from Grhl2CD-/- mice indicated a flattened expression gradient of osmo-responsive genes compared with control mice, consistent with their physiologic phenotype. CONCLUSIONS: Single-cell transcriptomics from dissociated kidneys facilitated the prediction of spatial gene expression along the corticomedullary axis and quantitation of osmotically regulated genes, allowing the prediction of a physiologic phenotype.

Psychrophilic proteases dramatically reduce single-cell RNA-seq artifacts: a molecular atlas of kidney development.

Single-cell RNA-seq is a powerful technique. Nevertheless, there are important limitations, including the technical challenges of breaking down an organ or tissue into a single-cell suspension. Invariably, this has required enzymatic incubation at 37°C, which can be expected to result in artifactual changes in gene expression patterns. Here, we describe a dissociation method that uses a protease with high activity in the cold, purified from a psychrophilic microorganism. The entire procedure is carried out at 6°C or colder, at which temperature mammalian transcriptional machinery is largely inactive, thereby effectively 'freezing in' the in vivo gene expression patterns. To test this method, we carried out RNA-seq on 20,424 single cells from postnatal day 1 mouse kidneys, comparing the results of the psychrophilic protease method with procedures using 37°C incubation. We show that the cold protease method provides a great reduction in gene expression artifacts. In addition, the results produce a single-cell resolution gene expression atlas of the newborn mouse kidney, an interesting time in development when mature nephrons are present yet nephrogenesis remains extremely active.

Cross-platform single cell analysis of kidney development shows stromal cells express Gdnf.

The developing kidney provides a useful model for study of the principles of organogenesis. In this report we use three independent platforms, Drop-Seq, Chromium 10x Genomics and Fluidigm C1, to carry out single cell RNA-Seq (scRNA-Seq) analysis of the E14.5 mouse kidney. Using the software AltAnalyze, in conjunction with the unsupervised approach ICGS, we were unable to identify and confirm the presence of 16 distinct cell populations during this stage of active nephrogenesis. Using a novel integrative supervised computational strategy, we were able to successfully harmonize and compare the cell profiles across all three technological platforms. Analysis of possible cross compartment receptor/ligand interactions identified the nephrogenic zone stroma as a source of GDNF. This was unexpected because the cap mesenchyme nephron progenitors had been thought to be the sole source of GDNF, which is a key driver of branching morphogenesis of the collecting duct system. The expression of Gdnf by stromal cells was validated in several ways, including Gdnf in situ hybridization combined with immunohistochemistry for SIX2, and marker of nephron progenitors, and MEIS1, a marker of stromal cells. Finally, the single cell gene expression profiles generated in this study confirmed and extended previous work showing the presence of multilineage priming during kidney development. Nephron progenitors showed stochastic expression of genes associated with multiple potential differentiation lineages.

A gene expression atlas of early craniofacial development.

We present a gene expression atlas of early mouse craniofacial development. Laser capture microdissection (LCM) was used to isolate cells from the principal critical microregions, whose development, differentiation and signaling interactions are responsible for the construction of the mammalian face. At E8.5, as migrating neural crest cells begin to exit the neural fold/epidermal ectoderm boundary, we examined the cranial mesenchyme, composed of mixed neural crest and paraxial mesoderm cells, as well as cells from adjacent neuroepithelium. At E9.5 cells from the cranial mesenchyme, overlying olfactory placode/epidermal ectoderm, and underlying neuroepithelium, as well as the emerging mandibular and maxillary arches were sampled. At E10.5, as the facial prominences form, cells from the medial and lateral prominences, the olfactory pit, multiple discrete regions of underlying neuroepithelium, the mandibular and maxillary arches, including both their mesenchymal and ectodermal components, as well as Rathke's pouch, were similarly sampled and profiled using both microarray and RNA-seq technologies. Further, we performed single cell studies to better define the gene expression states of the early E8.5 pioneer neural crest cells and paraxial mesoderm. Taken together, and analyzable by a variety of biological network approaches, these data provide a complementing and cross validating resource capable of fueling discovery of novel compartment specific markers and signatures whose combinatorial interactions of transcription factors and growth factors/receptors are responsible for providing the master genetic blueprint for craniofacial development.

Center volume effect on acute cellular rejection and outcomes in pediatric lung transplant recipients.

BACKGROUND: Acute cellular rejection (ACR) is common after lung transplant (LTx). We sought to determine if transplant center volume affected ACR-related outcomes in children after LTx. METHODS: The United Network for Organ Sharing (UNOS) Registry was queried for patients <18-years-of-age who underwent LTx 1987-2020. Cohorts were children who survived the first-year post transplant and were treated for ACR within that first year (ACR group) and those not treated for ACR (non-ACR). LTx center volume was defined as: high volume center (HVC) (>5LTxs/year), medium volume center (MVC) (>1≤5 LTxs/year), and low volume center (LVC) (≤1LTxs/year). RESULTS: 1320 patients were enrolled into the study; 269 (20.4%) did not experience ACR. The ACR cohort was older (median 14 [11-16] vs 13 [7-16] years, p < 0.001), female (65.3% vs 57.3%, p = 0.016), had cystic fibrosis (62.3% vs 45.5%, p < 0.001), and had a higher lung allocation score (37.3 [34.6-47.8] vs 35.8 [33-42.6], p = 0.029). The ACR cohort trended (p = 0.06) towards lower survival at 5-year (37% vs 47%) and 10-year (25% vs 34%) post-LTx. Among children at HVCs, ACR occurred in 17% of recipients (n = 98/574), compared to 18.5% (n = 73/395) at MVCs and 27% (n = 100/369) at LVCs. Children treated for ACR at HVCs had higher survival than LVCs at 5-years (52% vs 29%) and 10-years (36% vs 15%) (p < 0.001) but similar survival to MVCs at 5-years (52% vs 43%) and 10-years (36% vs 24%) (p = 0.081). No survival differences were detected in MVCs vs LVCs (p = 0.14). CONCLUSIONS: ACR treated within the first post-LTx year influence survival of children. ACR incidence was lowest at higher volume centers whereas post-ACR treatment survival outcomes were also superior.

Drinking carrot juice increases total antioxidant status and decreases lipid peroxidation in adults.

BACKGROUND: High prevalence of obesity and cardiovascular disease is attributable to sedentary lifestyle and eating diets high in fat and refined carbohydrate while eating diets low in fruit and vegetables. Epidemiological studies have confirmed a strong association between eating diets rich in fruits and vegetables and cardiovascular health. The aim of this pilot study was to determine whether drinking fresh carrot juice influences antioxidant status and cardiovascular risk markers in subjects not modifying their eating habits. METHODS: An experiment was conducted to evaluate the effects of consuming 16 fl oz of daily freshly squeezed carrot juice for three months on cardiovascular risk markers, C-reactive protein, insulin, leptin, interleukin-1α, body fat percentage, body mass index (BMI), blood pressure, antioxidant status, and malondialdehyde production. Fasting blood samples were collected pre-test and 90 days afterward to conclude the study. RESULTS: Drinking carrot juice did not affect (P > 0.1) the plasma cholesterol, triglycerides, Apo A, Apo B, LDL, HDL, body fat percentage, insulin, leptin, interleukin-1α, or C-reactive protein. Drinking carrot juice decreased (P = 0.06) systolic pressure, but did not influence diastolic pressure. Drinking carrot juice significantly (P < 0.05) increased the plasma total antioxidant capacity and decreased (P < 0.05) the plasma malondialdehyde production. CONCLUSION: Drinking carrot juice may protect the cardiovascular system by increasing total antioxidant status and by decreasing lipid peroxidation independent of any of the cardiovascular risk markers measured in the study.

Single-cell RNA-seq analysis reveals compartment-specific heterogeneity and plasticity of microglia.

Microglia are ubiquitous central nervous system (CNS)-resident macrophages that maintain homeostasis of neural tissues and protect them from pathogen attacks. Yet, their differentiation in different compartments remains elusive. We performed single-cell RNA-seq to compare microglial subtypes in the cortex and the spinal cord. A multi-way comparative analysis was carried out on samples from C57/BL and HIV gp120 transgenic mice at two, four, and eight months of age. The results revealed overlapping but distinct microglial populations in the cortex and the spinal cord. The differential heterogeneity of microglia in these CNS regions was further suggested by their disparity of plasticity in response to life span progression and HIV-1 pathogenic protein gp120. Our findings indicate that microglia in different CNS compartments are adapted to their local environments to fulfill region-specific biological functions.

Dissociation of Tissues for Single-Cell Analysis.

It is now a routine to carry out single-cell RNA-Seq to define the gene expression patterns of thousands of cells, thereby revolutionizing many areas of research. Projects are underway to use these techniques to create an atlas of the expressed genes in all cell types of the human body. Here we describe cold-active protease methods for single-cell dissociation of organs and tissues that better preserve the in vivo gene expression patterns.

A bigenic mouse model of FSGS reveals perturbed pathways in podocytes, mesangial cells and endothelial cells.

Focal segmental glomerulosclerosis is a major cause of end stage renal disease. Many patients prove unresponsive to available therapies. An improved understanding of the molecular basis of the disease process could provide insights leading to novel therapeutic approaches. In this study we carried out an RNA-seq analysis of the altered gene expression patterns of podocytes, mesangial cells and glomerular endothelial cells of the bigenic Cd2ap+/-, Fyn-/- mutant mouse model of FSGS. In the podocytes we observed upregulation of many genes related to the Tgfß family/pathway, including Gdnf, Tgfß1, Tgfß2, Snai2, Vegfb, Bmp4, and Tnc. The mutant podocytes also showed upregulation of Acta2, a marker of smooth muscle and associated with myofibroblasts, which are implicated in driving fibrosis. GO analysis of the podocyte upregulated genes identified elevated protein kinase activity, increased expression of growth factors, and negative regulation of cell adhesion, perhaps related to the observed podocyte loss. Both podocytes and mesangial cells showed strong upregulation of aldehyde dehydrogenase genes involved in the synthesis of retinoic acid. Similarly, the Cd2ap+/-, Fyn-/- mesangial cells, as well as podocytes in other genetic models, and the glomeruli of human FSGS patients, all show upregulation of the serine protease Prss23, with the common thread suggesting important functionality. Another gene with strong upregulation in the Cd2ap+/-, Fyn-/- mutant mesangial cells as well as multiple other mutant mouse models of FSGS was thrombospondin, which activates the secreted inactive form of Tgfß. The Cd2ap+/-, Fyn-/- mutant endothelial cells showed elevated expression of genes involved in cell proliferation, angioblast migration, angiogenesis, and neovasculature, all consistent with the formation of new blood vessels in the diseased glomerulus. The resulting global definition of the perturbed molecular pathways in the three major cell types of the mutant glomerulus provide deeper understanding of the molecular pathogenic pathways.

Single cell RNA-seq study of wild type and Hox9,10,11 mutant developing uterus.

The uterus is a remarkable organ that must guard against infections while maintaining the ability to support growth of a fetus without rejection. The Hoxa10 and Hoxa11 genes have previously been shown to play essential roles in uterus development and function. In this report we show that the Hoxa9,10,11, Hoxc9,10,11, Hoxd9,10,11 genes play a redundant role in the formation of uterine glands. In addition, we use single cell RNA-seq to create a high resolution gene expression atlas of the developing wild type mouse uterus. Cell types and subtypes are defined, for example dividing endothelial cells into arterial, venous, capillary, and lymphatic, while epithelial cells separate into luminal and glandular subtypes. Further, a surprising heterogeneity of stromal and myocyte cell types are identified. Transcription factor codes and ligand/receptor interactions are characterized. We also used single cell RNA-seq to globally define the altered gene expression patterns in all developing uterus cell types for two Hox mutants, with 8 or 9 mutant Hox genes. The mutants show a striking disruption of Wnt signaling as well as the Cxcl12/Cxcr4 ligand/receptor axis.

Single-Cell Transcriptomics in Medulloblastoma Reveals Tumor-Initiating Progenitors and Oncogenic Cascades during Tumorigenesis and Relapse.

Progenitor heterogeneity and identities underlying tumor initiation and relapse in medulloblastomas remain elusive. Utilizing single-cell transcriptomic analysis, we demonstrated a developmental hierarchy of progenitor pools in Sonic Hedgehog (SHH) medulloblastomas, and identified OLIG2-expressing glial progenitors as transit-amplifying cells at the tumorigenic onset. Although OLIG2+ progenitors become quiescent stem-like cells in full-blown tumors, they are highly enriched in therapy-resistant and recurrent medulloblastomas. Depletion of mitotic Olig2+ progenitors or Olig2 ablation impeded tumor initiation. Genomic profiling revealed that OLIG2 modulates chromatin landscapes and activates oncogenic networks including HIPPO-YAP/TAZ and AURORA-A/MYCN pathways. Co-targeting these oncogenic pathways induced tumor growth arrest. Together, our results indicate that glial lineage-associated OLIG2+ progenitors are tumor-initiating cells during medulloblastoma tumorigenesis and relapse, suggesting OLIG2-driven oncogenic networks as potential therapeutic targets.

Molecular Anatomy of Palate Development.

The NIH FACEBASE consortium was established in part to create a central resource for craniofacial researchers. One purpose is to provide a molecular anatomy of craniofacial development. To this end we have used a combination of laser capture microdissection and RNA-Seq to define the gene expression programs driving development of the murine palate. We focused on the E14.5 palate, soon after medial fusion of the two palatal shelves. The palate was divided into multiple compartments, including both medial and lateral, as well as oral and nasal, for both the anterior and posterior domains. A total of 25 RNA-Seq datasets were generated. The results provide a comprehensive view of the region specific expression of all transcription factors, growth factors and receptors. Paracrine interactions can be inferred from flanking compartment growth factor/receptor expression patterns. The results are validated primarily through very high concordance with extensive previously published gene expression data for the developing palate. In addition selected immunostain validations were carried out. In conclusion, this report provides an RNA-Seq based atlas of gene expression patterns driving palate development at microanatomic resolution. This FACEBASE resource is designed to promote discovery by the craniofacial research community.

Drinking orange juice increases total antioxidant status and decreases lipid peroxidation in adults.

Cardiovascular disease (CVD) is the leading cause of death in the world and is the primary cause of mortality among Americans. One of the many reasons for the pathogenesis of CVD is attributed to eating diets high in saturated fat and refined carbohydrates and low in fruits and vegetables. Epidemiological evidence has supported a strong association between eating diets rich in fruits and vegetables and cardiovascular health. An experiment was conducted utilizing 24 adults with hypercholesterolemia and hypertriglyceridemia to evaluate the impact of drinking 20 fl oz of freshly squeezed orange juice daily for 90 days on blood pressure, lipid panels, plasma antioxidant capacity, metabolic hormones, lipid peroxidation, and inflammatory markers. Except for addition of drinking orange juice, subjects did not modify their eating habits. The findings suggested that drinking orange juice does not affect (P>.1) blood pressure, lipid panels, metabolic hormones, body fat percentage, or inflammatory markers. However, total plasma antioxidant capacity was significantly increased (P<.05) and lipid peroxidation was significantly decreased (P<.05) after orange juice consumption. Drinking orange juice may protect the cardiovascular system by increasing total plasma antioxidant status and by lowering lipid peroxidation independent of other cardiovascular risk markers evaluated in this study.