Micro-dissection and integration of long and short reads to create a robust catalog of kidney compartment-specific isoforms.

Studying isoform expression at the microscopic level has always been a challenging task. A classical example is kidney, where glomerular and tubulo-interstitial 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 in the glomerular and tubulo-interstitial compartments, respectively. Of note, relying solely on either short or long reads would have resulted in many erroneous identifications. We identified distinct pathways involved in glomerular and tubulo-interstitial compartments at the isoform level, creating an important experimental and computational resource for the kidney research community.

Gene expression profiles of diabetic kidney disease and neuropathy in eNOS knockout mice: Predictors of pathology and RAS blockade effects.

Diabetic kidney disease (DKD) and diabetic peripheral neuropathy (DPN) are two common diabetic complications. However, their pathogenesis remains elusive and current therapies are only modestly effective. We evaluated genome-wide expression to identify pathways involved in DKD and DPN progression in db/db eNOS-/- mice receiving renin-angiotensin-aldosterone system (RAS)-blocking drugs to mimic the current standard of care for DKD patients. Diabetes and eNOS deletion worsened DKD, which improved with RAS treatment. Diabetes also induced DPN, which was not affected by eNOS deletion or RAS blockade. Given the multiple factors affecting DKD and the graded differences in disease severity across mouse groups, an automatic data analysis method, SOM, or self-organizing map was used to elucidate glomerular transcriptional changes associated with DKD, whereas pairwise bioinformatic analysis was used for DPN. These analyses revealed that enhanced gene expression in several pro-inflammatory networks and reduced expression of development genes correlated with worsening DKD. Although RAS treatment ameliorated the nephropathy phenotype, it did not alter the more abnormal gene expression changes in kidney. Moreover, RAS exacerbated expression of genes related to inflammation and oxidant generation in peripheral nerves. The graded increase in inflammatory gene expression and decrease in development gene expression with DKD progression underline the potentially important role of these pathways in DKD pathogenesis. Since RAS blockers worsened this gene expression pattern in both DKD and DPN, it may partly explain the inadequate therapeutic efficacy of such blockers.

Genome-Wide Functional Annotation of Human Protein-Coding Splice Variants Using Multiple Instance Learning.

The vast majority of human multiexon genes undergo alternative splicing and produce a variety of splice variant transcripts and proteins, which can perform different functions. These protein-coding splice variants (PCSVs) greatly increase the functional diversity of proteins. Most functional annotation algorithms have been developed at the gene level; the lack of isoform-level gold standards is an important intellectual limitation for currently available machine learning algorithms. The accumulation of a large amount of RNA-seq data in the public domain greatly increases our ability to examine the functional annotation of genes at isoform level. In the present study, we used a multiple instance learning (MIL)-based approach for predicting the function of PCSVs. We used transcript-level expression values and gene-level functional associations from the Gene Ontology database. A support vector machine (SVM)-based 5-fold cross-validation technique was applied. Comparatively, genes with multiple PCSVs performed better than single PCSV genes, and performance also improved when more examples were available to train the models. We demonstrated our predictions using literature evidence of ADAM15, LMNA/C, and DMXL2 genes. All predictions have been implemented in a web resource called "IsoFunc", which is freely available for the global scientific community through http://guanlab.ccmb.med.umich.edu/isofunc .

Computational Inferences of the Functions of Alternative/Noncanonical Splice Isoforms Specific to HER2+/ER-/PR- Breast Cancers, a Chromosome 17 C-HPP Study.

This study was conducted as a part of the Chromosome-Centric Human Proteome Project (C-HPP) of the Human Proteome Organization. The main objective is to identify and evaluate functionality of a set of specific noncanonical isoforms expressed in HER2-neu positive, estrogen receptor negative (ER-), and progesterone receptor negative (PR-) breast cancers (HER2+/ER-/PR- BC), an aggressive subtype of breast cancers that cause significant morbidity and mortality. We identified 11 alternative splice isoforms that were differentially expressed in HER2+/ER-/PR- BC compared to normal mammary, triple negative breast cancer and triple positive breast cancer tissues (HER2+/ER+/PR+). We used a stringent criterion that differentially expressed noncanonical isoforms (adjusted p value < 0.05) and have to be expressed in all replicates of HER2+/ER-/PR- BC samples, and the trend in differential expression (up or down) is the same in all comparisons. Of the 11 protein isoforms, six were overexpressed in HER2+/ER-/PR- BC. We explored possible functional roles of these six proteins using several complementary computational tools. Biological processes including cell cycle events and glycolysis were linked to four of these proteins. For example, glycolysis was the top ranking functional process for DMXL2 isoform 3, with a fold change of 27 compared to just two for the canonical protein. No previous reports link DMXL2 with any metabolic processes; the canonical protein is known to participate in signaling pathways. Our results clearly indicate distinct functions for the six overexpressed alternative splice isoforms, and these functions could be specific to HER2+/ER-/PR- tumor progression. Further detailed analysis is warranted as these proteins could be explored as potential biomarkers and therapeutic targets for HER2+/ER-/PR- BC patients.

MI-PVT: A Tool for Visualizing the Chromosome-Centric Human Proteome.

We have developed the web-based Michigan Proteome Visualization Tool (MI-PVT) to visualize and compare protein expression and isoform-level function across human chromosomes and tissues (http://guanlab.ccmb.med.umich.edu/mipvt). As proof of principle, we have populated the tool with Human Proteome Map (HPM) data. We were able to observe many biologically interesting features. From the vantage point of our chromosome 17 team, for example, we found more than 300 proteins from chromosome 17 expressed in each of the 30 tissues and cell types studied, with the highest number of expressed proteins being 685 in testis. Comparisons of expression levels across tissues showed low numbers of proteins expressed in esophagus, but esophagus had 12 cytoskeletal proteins coded on chromosome 17 with very high expression (>1000 spectral counts). This customized MI-PVT should be helpful for biologists to browse and study specific proteins and protein data sets across tissues and chromosomes. Users can upload any data of interest in MI-PVT for visualization. Our aim is to integrate extensive mass-spectrometric proteomic data into the tool to facilitate finding chromosome-centric protein expression and correlation across tissues.

Modeling dynamic functional relationship networks and application to ex vivo human erythroid differentiation.

MOTIVATION: Functional relationship networks, which summarize the probability of co-functionality between any two genes in the genome, could complement the reductionist focus of modern biology for understanding diverse biological processes in an organism. One major limitation of the current networks is that they are static, while one might expect functional relationships to consistently reprogram during the differentiation of a cell lineage. To address this potential limitation, we developed a novel algorithm that leverages both differentiation stage-specific expression data and large-scale heterogeneous functional genomic data to model such dynamic changes. We then applied this algorithm to the time-course RNA-Seq data we collected for ex vivo human erythroid cell differentiation. RESULTS: Through computational cross-validation and literature validation, we show that the resulting networks correctly predict the (de)-activated functional connections between genes during erythropoiesis. We identified known critical genes, such as HBD and GATA1, and functional connections during erythropoiesis using these dynamic networks, while the traditional static network was not able to provide such information. Furthermore, by comparing the static and the dynamic networks, we identified novel genes (such as OSBP2 and PDZK1IP1) that are potential drivers of erythroid cell differentiation. This novel method of modeling dynamic networks is applicable to other differentiation processes where time-course genome-scale expression data are available, and should assist in generating greater understanding of the functional dynamics at play across the genome during development. AVAILABILITY AND IMPLEMENTATION: The network described in this article is available at http://guanlab.ccmb.med.umich.edu/stageSpecificNetwork.

Systematically differentiating functions for alternatively spliced isoforms through integrating RNA-seq data.

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 gene-centered function prediction to isoform-level predictions.

A Network of Splice Isoforms for the Mouse.

The laboratory mouse is the primary mammalian species used for studying alternative splicing events. Recent studies have generated computational models to predict functions for splice isoforms in the mouse. However, the functional relationship network, describing the probability of splice isoforms participating in the same biological process or pathway, has not yet been studied in the mouse. Here we describe a rich genome-wide resource of mouse networks at the isoform level, which was generated using a unique framework that was originally developed to infer isoform functions. This network was built through integrating heterogeneous genomic and protein data, including RNA-seq, exon array, protein docking and pseudo-amino acid composition. Through simulation and cross-validation studies, we demonstrated the accuracy of the algorithm in predicting isoform-level functional relationships. We showed that this network enables the users to reveal functional differences of the isoforms of the same gene, as illustrated by literature evidence with Anxa6 (annexin a6) as an example. We expect this work will become a useful resource for the mouse genetics community to understand gene functions. The network is publicly available at: http://guanlab.ccmb.med.umich.edu/isoformnetwork.

Crowdsourced assessment of common genetic contribution to predicting anti-TNF treatment response in rheumatoid arthritis.

Rheumatoid arthritis (RA) affects millions world-wide. While anti-TNF treatment is widely used to reduce disease progression, treatment fails in ∼one-third of patients. No biomarker currently exists that identifies non-responders before treatment. A rigorous community-based assessment of the utility of SNP data for predicting anti-TNF treatment efficacy in RA patients was performed in the context of a DREAM Challenge (http://www.synapse.org/RA_Challenge). An open challenge framework enabled the comparative evaluation of predictions developed by 73 research groups using the most comprehensive available data and covering a wide range of state-of-the-art modelling methodologies. Despite a significant genetic heritability estimate of treatment non-response trait (h(2)=0.18, P value=0.02), no significant genetic contribution to prediction accuracy is observed. Results formally confirm the expectations of the rheumatology community that SNP information does not significantly improve predictive performance relative to standard clinical traits, thereby justifying a refocusing of future efforts on collection of other data.

Identification of stage-specific genes associated with lupus nephritis and response to remission induction in (NZB × NZW)F1 and NZM2410 mice.

OBJECTIVE: To elucidate the molecular mechanisms involved in renal inflammation during the progression, remission, and relapse of nephritis in murine lupus models using transcriptome analysis. METHODS: Kidneys from (NZB × NZW)F1 (NZB/NZW) and NZM2410 mice were harvested at intervals during the disease course or after remission induction. Genome-wide expression profiles were obtained from microarray analysis of perfused kidneys. Real-time polymerase chain reaction (PCR) analysis for selected genes was used to validate the microarray data. Comparisons between groups using SAM, and unbiased analysis of the entire data set using singular value decomposition and self-organizing maps were performed. RESULTS: Few changes in the renal molecular profile were detected in prenephritic kidneys, but a significant shift in gene expression, reflecting inflammatory cell infiltration and complement activation, occurred at proteinuria onset. Subsequent changes in gene expression predominantly affected mitochondrial dysfunction and metabolic stress pathways. Endothelial cell activation, tissue remodeling, and tubular damage were the major pathways associated with loss of renal function. Remission induction reversed most, but not all, of the inflammatory changes, and progression toward relapse was associated with recurrence of inflammation, mitochondrial dysfunction, and metabolic stress signatures. CONCLUSION: Immune cell infiltration and activation is associated with proteinuria onset and is reversed by immunosuppressive therapy, but disease progression is associated with renal hypoxia and metabolic stress. Optimal therapy for lupus nephritis may therefore need to target both immune and nonimmune disease mechanisms. In addition, the overlap of a substantial subset of molecular markers with those expressed in the kidneys of lupus patients suggests potential new biomarkers and therapeutic targets.