Rajasree Menon | University of Michigan (original) (raw)

Papers by Rajasree Menon

PLoS Computational Biology, 2013

Integrating large-scale functional genomic data has significantly accelerated our understanding o... more Integrating large-scale functional genomic data has significantly accelerated our understanding of gene functions. However, no algorithm has been developed to differentiate functions for isoforms of the same gene using high-throughput genomic data. This is because standard supervised learning requires 'ground-truth' functional annotations, which are lacking at the isoform level. To address this challenge, we developed a generic framework that interrogates public RNA-seq data at the transcript level to differentiate functions for alternatively spliced isoforms. For a specific function, our algorithm identifies the 'responsible' isoform(s) of a gene and generates classifying models at the isoform level instead of at the gene level. Through cross-validation, we demonstrated that our algorithm is effective in assigning functions to genes, especially the ones with multiple isoforms, and robust to gene expression levels and removal of homologous gene pairs. We identified genes in the mouse whose isoforms are predicted to have disparate functionalities and experimentally validated the 'responsible' isoforms using data from mammary tissue. With protein structure modeling and experimental evidence, we further validated the predicted isoform functional differences for the genes Cdkn2a and Anxa6. Our generic framework is the first to predict and differentiate functions for alternatively spliced isoforms, instead of genes, using genomic data. It is extendable to any base machine learner and other species with alternatively spliced isoforms, and shifts the current genecentered function prediction to isoform-level predictions.

Kidney International, 2020

Kidney International, 2021

Increased podocyte detachment begins immediately after kidney transplantation and is associated w... more Increased podocyte detachment begins immediately after kidney transplantation and is associated with long-term allograft failure. We hypothesized that cell-specific transcriptional changes in podocytes and glomerular endothelial cells after transplantation would offer mechanistic insights into the podocyte detachment process. To test this, we evaluated cell-specific transcriptional profiles of glomerular endothelial cells and podocytes from 14 patients of their first-year surveillance biopsies with normal histology from low immune risk recipients with no post-transplant complications and compared these to biopsies of 20 healthy living donor controls. Glomerular endothelial cells from these surveillance biopsies were enriched for genes related to fluid shear stress, angiogenesis, and interferon signaling. In podocytes, pathways were enriched for genes in response to growth factor signaling and actin cytoskeletal reorganization, but also showed evidence of podocyte stress as indicated by reduced nephrin (adhesion protein) gene expression. In parallel, transcripts coding for proteins required to maintain podocyte adherence to the underlying glomerular basement membrane were down-regulated including the major glomerular podocyte integrin α3 and the actin cytoskeleton-related gene synaptopodin. The reduction in integrin α3 protein expression in surveillance biopsies was confirmed by immunoperoxidase staining. The combined growth and stress response of patient allografts post-transplantation paralleled similar changes in a rodent model of nephrectomy-induced glomerular hypertrophic stress that progress to develop proteinuria and glomerulosclerosis with shortened kidney lifespan. Thus, even among patients with apparently healthy allografts with no detectable histologic abnormality including alloimmune injury, transcriptomic changes reflecting cell stresses are already set in motion that could drive hypertrophy-associated glomerular disease progression.

Disease Markers, 2010

Alternative splicing plays an important role in protein diversity without increasing genome size.... more Alternative splicing plays an important role in protein diversity without increasing genome size. Earlier thought to be uncommon, splicing appears to affect the majority of genes. Alternative splice variants have been detected at the mRNA level in many diseases. We have designed and demonstrated a discovery pipeline for alternative splice variant (ASV) proteins from tandem MS/MS datasets. We created a modified ECgene database with entries from exhaustive three-frame translation of Ensembl transcripts and gene models from ECgene, with periodic updates. The human database has 14 million entries; the mouse database, 10 million entries. We match MS/MS findings against these potential translation products to identify and quantify known and novel ASVs. In this review, we summarize findings and systems biology implications of biomarker candidates from a mouse model of human pancreatic ductal adenocarcinoma [28] and a mouse model of human Her2/neu-induced breast cancer [27]. The same approa...

JCI Insight, Mar 26, 2020

Kidney Precision Medicine Project (KPMP) is building a spatially-specified human tissue atlas at ... more Kidney Precision Medicine Project (KPMP) is building a spatially-specified human tissue atlas at the single-cell resolution with molecular details of the kidney in health and disease. Here, we describe the construction of an integrated reference tissue map of cells, pathways and genes using unaffected regions of nephrectomy tissues and undiseased human biopsies from 55 subjects. We use single-cell and -nucleus transcriptomics, subsegmental laser microdissection bulk transcriptomics and proteomics, near-single-cell proteomics, 3-D nondestructive and CODEX imaging, and spatial metabolomics data to hierarchically identify genes, pathways and cells. Integrated data from these different technologies coherently describe cell types/subtypes within different nephron segments and interstitium. These spatial profiles identify cell-level functional organization of the kidney tissue as indicative of their physiological functions and map different cell subtypes to genes, proteins, metabolites an...

Background: Classification of nephrotic syndrome relies on clinical presentation and descriptive ... more Background: Classification of nephrotic syndrome relies on clinical presentation and descriptive patterns of injury on kidney biopsies. This approach does not reflect underlying disease biology, limiting the ability to predict progression or treatment response. Methods: Systems biology approaches were used to categorize patients with minimal change disease (MCD) and focal segmental glomerulosclerosis (FSGS) based on kidney biopsy tissue transcriptomics across three cohorts and assessed association with clinical outcomes. Patient- level tissue pathway activation scores were generated using differential gene expression. Then, functional enrichment and non-invasive urine biomarker candidates were identified. Biomarkers were validated in kidney organoid models and single nucleus RNA-seq (snRNAseq) from kidney biopsies. Results: Transcriptome-based categorization identified three subgroups of patients with shared molecular signatures across independent North American, European and Africa...

Studying isoform expression at the microscopic level has always been a challenging task. A classi... more Studying isoform expression at the microscopic level has always been a challenging task. A classical example is kidney, where glomerular and tubulo-insterstitial compartments carry out drastically different physiological functions and thus presumably their isoform expression also differs. We aim at developing an experimental and computational pipeline for identifying isoforms at microscopic structure-level. We microdissected glomerular and tubulo-interstitial compartments from healthy human kidney tissues from two cohorts. The two compartments were separately sequenced with the PacBio RS II platform. These transcripts were then validated using transcripts of the same samples by the traditional Illumina RNA-Seq protocol, distinct Illumina RNA-Seq short reads from European Renal cDNA Bank (ERCB) samples, and annotated GENCODE transcript list, thus identifying novel transcripts. We identified 14,739 and 14,259 annotated transcripts, and 17,268 and 13,118 potentially novel transcripts i...

American Journal of Physiology-Renal Physiology

Loss-of-function mutations in phospholipase C-ε1 (PLCE1) have been detected in patients with neph... more Loss-of-function mutations in phospholipase C-ε1 (PLCE1) have been detected in patients with nephrotic syndrome, but other family members with the same mutation were asymptomatic, suggesting additional stressor are required to cause the full phenotype. Consistent with these observations, we determined that global Plce1-deficient mice have histologically normal glomeruli and no albuminuria at baseline. Angiotensin II (ANG II) is known to induce glomerular damage in genetically susceptible individuals. Therefore, we tested whether ANG II enhances glomerular damage in Plce1-deficient mice. ANG II increased blood pressure equally in Plce1-deficient and wild-type littermates. Additionally, it led to 20-fold increased albuminuria and significantly more sclerotic glomeruli in Plce1-deficient mice compared with wild-type littermates. Furthermore, Plce1-deficient mice demonstrated diffuse mesangial expansion, podocyte loss, and focal podocyte foot process effacement. To determine whether the...

Nature Communications

Heterotopic ossification (HO) is an aberrant regenerative process with ectopic bone induction in ... more Heterotopic ossification (HO) is an aberrant regenerative process with ectopic bone induction in response to musculoskeletal trauma, in which mesenchymal stem cells (MSC) differentiate into osteochondrogenic cells instead of myocytes or tenocytes. Despite frequent cases of hospitalized musculoskeletal trauma, the inflammatory responses and cell population dynamics that regulate subsequent wound healing and tissue regeneration are still unclear. Here we examine, using a mouse model of trauma-induced HO, the local microenvironment of the initial post-injury inflammatory response. Single cell transcriptome analyses identify distinct monocyte/macrophage populations at the injury site, with their dynamic changes over time elucidated using trajectory analyses. Mechanistically, transforming growth factor beta-1 (TGFβ1)-producing monocytes/macrophages are associated with HO and aberrant chondrogenic progenitor cell differentiation, while CD47-activating peptides that reduce systemic macroph...

Blood

Blood is a complex fluid that samples all tissues in the human body. Despite complete sequence de... more Blood is a complex fluid that samples all tissues in the human body. Despite complete sequence determination of the human genome, defining genes and gene products remains a challenge. Here, we apply tandem mass spectroscopy as new source of unbiased data to interrogate genomic sequence and identify novel protein coding sequences. A six-frame translation of the Human genome was used as the query database to search for novel blood proteins in the data from the HUPO PPP. Significance is assessed using a Poisson statistical model incorporating the length of the matching sequence and the frequency of spectrum matches observed in searching the database [Nat Biotech 2006 24(3):333–8]. Matches are binned by X!Tandem hyperscore, and statistics for each score class are considered independently. The overall probability that the matches to an ORF occurred at random is calculated as the product of the probability that the matches in each score category occurred at random. The expected number of ...

Journal of the American Society of Nephrology

Podocyte injury is central to many forms of kidney disease, but transcriptional signatures reflec... more Podocyte injury is central to many forms of kidney disease, but transcriptional signatures reflecting podocyte injury and compensation mechanisms are challenging to analyze in vivo. Human kidney organoids derived from pluripotent stem cells (PSCs), a new model for disease and regeneration, present an opportunity to explore the transcriptional plasticity of podocytes. Here, transcriptional profiling of over 12,000 single cells from human PSC-derived kidney organoid cultures was used to identify robust and reproducible cell-lineage gene expression signatures shared with developing human kidneys based on trajectory analysis. Surprisingly, the gene expression signature characteristic of developing glomerular epithelial cells was also observed in glomerular tissue from a kidney disease cohort. This signature correlated with proteinuria and inverse eGFR, and was confirmed in an independent podocytopathy cohort. Three genes in particular, LYPD1, PRSS23 and CDH6, were further identified as ...

Development (Cambridge, England), Aug 30, 2018

The mammalian kidney develops through reciprocal interactions between the ureteric bud and the me... more The mammalian kidney develops through reciprocal interactions between the ureteric bud and the metanephric mesenchyme to give rise to the entire collecting system and the nephrons. Most of our knowledge of the developmental regulators driving this process arises from the study of gene expression and functional genetics in mice and other animal models. In order to shed light on human kidney development, we have used single-cell transcriptomics to characterize gene expression in different cell populations, and to study individual cell dynamics and lineage trajectories during development. Single-cell transcriptome analyses of 6414 cells from five individual specimens identified 11 initial clusters of specific renal cell types as defined by their gene expression profile. Further subclustering identifies progenitors, and mature and intermediate stages of differentiation for several renal lineages. Other lineages identified include mesangium, stroma, endothelial and immune cells. Novel ma...

American journal of human genetics, Jan 2, 2018

Expression quantitative trait loci (eQTL) studies illuminate the genetics of gene expression and,... more Expression quantitative trait loci (eQTL) studies illuminate the genetics of gene expression and, in disease research, can be particularly illuminating when using the tissues directly impacted by the condition. In nephrology, there is a paucity of eQTL studies of human kidney. Here, we used whole-genome sequencing (WGS) and microdissected glomerular (GLOM) and tubulointerstitial (TI) transcriptomes from 187 individuals with nephrotic syndrome (NS) to describe the eQTL landscape in these functionally distinct kidney structures. Using MatrixEQTL, we performed cis-eQTL analysis on GLOM (n = 136) and TI (n = 166). We used the Bayesian "Deterministic Approximation of Posteriors" (DAP) to fine-map these signals, eQTLBMA to discover GLOM- or TI-specific eQTLs, and single-cell RNA-seq data of control kidney tissue to identify the cell type specificity of significant eQTLs. We integrated eQTL data with an IgA Nephropathy (IgAN) GWAS to perform a transcriptome-wide association study...

Cell stem cell, Jan 15, 2018

Organoids derived from human pluripotent stem cells are a potentially powerful tool for high-thro... more Organoids derived from human pluripotent stem cells are a potentially powerful tool for high-throughput screening (HTS), but the complexity of organoid cultures poses a significant challenge for miniaturization and automation. Here, we present a fully automated, HTS-compatible platform for enhanced differentiation and phenotyping of human kidney organoids. The entire 21-day protocol, from plating to differentiation to analysis, can be performed automatically by liquid-handling robots, or alternatively by manual pipetting. High-content imaging analysis reveals both dose-dependent and threshold effects during organoid differentiation. Immunofluorescence and single-cell RNA sequencing identify previously undetected parietal, interstitial, and partially differentiated compartments within organoids and define conditions that greatly expand the vascular endothelium. Chemical modulation of toxicity and disease phenotypes can be quantified for safety and efficacy prediction. Screening in ge...

Methods in molecular biology (Clifton, N.J.), 2017

Tens of thousands of splice isoforms of proteins have been catalogued as predicted sequences from... more Tens of thousands of splice isoforms of proteins have been catalogued as predicted sequences from transcripts in humans and other species. Relatively few have been characterized biochemically or structurally. With the extensive development of protein bioinformatics, the characterization and modeling of isoform features, isoform functions, and isoform-level networks have advanced notably. Here we present applications of the I-TASSER family of algorithms for folding and functional predictions and the IsoFunc, MIsoMine, and Hisonet data resources for isoform-level analyses of network and pathway-based functional predictions and protein-protein interactions. Hopefully, predictions and insights from protein bioinformatics will stimulate many experimental validation studies.

Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research, Jan 8, 2017

Previously, we showed that cortical mineralization is coordinately adjusted to mechanically offse... more Previously, we showed that cortical mineralization is coordinately adjusted to mechanically offset external bone size differences between A/J (narrow) and C57BL/6J (wide) mouse femora to achieve whole bone strength equivalence at adulthood. The identity of the genes and their interactions that are responsible for establishing this homeostatic state (i.e., canalization) remain unknown. We hypothesize that these inbred strains, whose inter-individual differences in bone structure and material properties mimic that observed among humans, achieve functional homeostasis by differentially adjusting key molecular pathways regulating external bone size and mineralization throughout growth. The cortices of A/J and C57BL/6J male mouse femora were phenotyped and gene expression levels were assessed across growth (i.e. 2, 4, 6, 8, 12, 16-weeks of age). A difference in total cross-sectional area (p < 0.01) and cortical tissue mineral density were apparent between mouse strains by 2-weeks of a...