Surveying mutation density patterns around specific genomic features.

Mutation density patterns reveal unique biological properties of specific genomic regions and shed light on the mechanisms of carcinogenesis. Although previous studies reported insightful mutation density patterns associated with certain genomic regions such as transcription start sites and DNA replication origins, a tool that can systematically investigate mutational spatial patterns is still lacking. Thus, we developed MutDens, a bioinformatic tool for comprehensive analysis of mutation density patterns around genomic features, namely, genomic positions, in humans and model species. By scanning the bidirectional vicinity regions of given positions, MutDens systematically characterizes the mutation density for single-base substitution mutational classes after adjusting for total mutation burden and local nucleotide proportion. Analysis results using MutDens not only verified the previously reported transcriptional strand bias around transcription start sites and replicative strand bias around DNA replication origins, but also identified novel mutation density patterns around other genomics features, such as enhancers and retrotransposon insertion polymorphism sites. To our knowledge, MutDens is the first tool that systematically calculates, examines, and compares mutation density patterns, thus providing a valuable avenue for investigating the mutational landscapes associated with important genomic features.

Somatic mutation effects diffused over microRNA dysregulation.

MOTIVATION: As an important player in transcriptome regulation, microRNAs may effectively diffuse somatic mutation impacts to broad cellular processes and ultimately manifest disease and dictate prognosis. Previous studies that tried to correlate mutation with gene expression dysregulation neglected to adjust for the disparate multitudes of false positives associated with unequal sample sizes and uneven class balancing scenarios. RESULTS: To properly address this issue, we developed a statistical framework to rigorously assess the extent of mutation impact on microRNAs in relation to a permutation-based null distribution of a matching sample structure. Carrying out the framework in a pan-cancer study, we ascertained 9008 protein-coding genes with statistically significant mutation impacts on miRNAs. Of these, the collective miRNA expression for 83 genes showed significant prognostic power in nine cancer types. For example, in lower-grade glioma, 10 genes' mutations broadly impacted miRNAs, all of which showed prognostic value with the corresponding miRNA expression. Our framework was further validated with functional analysis and augmented with rich features including the ability to analyze miRNA isoforms; aggregative prognostic analysis; advanced annotations such as mutation type, regulator alteration, somatic motif, and disease association; and instructive visualization such as mutation OncoPrint, Ideogram, and interactive mRNA-miRNA network. AVAILABILITY AND IMPLEMENTATION: The data underlying this article are available in MutMix, at http://innovebioinfo.com/Database/TmiEx/MutMix.php.

SBSA: an online service for somatic binding sequence annotation.

Efficient annotation of alterations in binding sequences of molecular regulators can help identify novel candidates for mechanisms study and offer original therapeutic hypotheses. In this work, we developed Somatic Binding Sequence Annotator (SBSA) as a full-capacity online tool to annotate altered binding motifs/sequences, addressing diverse types of genomic variants and molecular regulators. The genomic variants can be somatic mutation, single nucleotide polymorphism, RNA editing, etc. The binding motifs/sequences involve transcription factors (TFs), RNA-binding proteins, miRNA seeds, miRNA-mRNA 3'-UTR binding target, or can be any custom motifs/sequences. Compared to similar tools, SBSA is the first to support miRNA seeds and miRNA-mRNA 3'-UTR binding target, and it unprecedentedly implements a personalized genome approach that accommodates joint adjacent variants. SBSA is empowered to support an indefinite species, including preloaded reference genomes for SARS-Cov-2 and 25 other common organisms. We demonstrated SBSA by annotating multi-omics data from over 30,890 human subjects. Of the millions of somatic binding sequences identified, many are with known severe biological repercussions, such as the somatic mutation in TERT promoter region which causes a gained binding sequence for E26 transformation-specific factor (ETS1). We further validated the function of this TERT mutation using experimental data in cancer cells. Availability:http://innovebioinfo.com/Annotation/SBSA/SBSA.php.

Alternative Applications of Genotyping Array Data Using Multivariant Methods.

One of the forerunners that pioneered the revolution of high-throughput genomic technologies is the genotyping microarray technology, which can genotype millions of single-nucleotide variants simultaneously. Owing to apparent benefits, such as high speed, low cost, and high throughput, the genotyping array has gained lasting applications in genome-wide association studies (GWAS) and thus accumulated an enormous amount of data. Empowered by continuous manufactural upgrades and analytical innovation, unconventional applications of genotyping array data have emerged to address more diverse genetic problems, holding promise of boosting genetic research into human diseases through the re-mining of the rich accumulated data. Here, we review several unconventional genotyping array analysis techniques that have been built on the idea of large-scale multivariant analysis and provide empirical application examples. These unconventional outcomes of genotyping arrays include polygenic score, runs of homozygosity (ROH)/heterozygosity ratio, distant pedigree computation, and mitochondrial DNA (mtDNA) copy number inference.

Deep neural network based tissue deconvolution of circulating tumor cell RNA.

Prior research has shown that the deconvolution of cell-free RNA can uncover the tissue origin. The conventional deconvolution approaches rely on constructing a reference tissue-specific gene panel, which cannot capture the inherent variation present in actual data. To address this, we have developed a novel method that utilizes a neural network framework to leverage the entire training dataset. Our approach involved training a model that incorporated 15 distinct tissue types. Through one semi-independent and two complete independent validations, including deconvolution using a semi in silico dataset, deconvolution with a custom normal tissue mixture RNA-seq data, and deconvolution of longitudinal circulating tumor cell RNA-seq (ctcRNA) data from a cancer patient with metastatic tumors, we demonstrate the efficacy and advantages of the deep-learning approach which were exerted by effectively capturing the inherent variability present in the dataset, thus leading to enhanced accuracy. Sensitivity analyses reveal that neural network models are less susceptible to the presence of missing data, making them more suitable for real-world applications. Moreover, by leveraging the concept of organotropism, we applied our approach to trace the migration of circulating tumor cell-derived RNA (ctcRNA) in a cancer patient with metastatic tumors, thereby highlighting the potential clinical significance of early detection of cancer metastasis.

Cancer-specific expression quantitative loci are affected by expression dysregulation.

Expression quantitative trait loci (eQTLs) have been touted as the missing piece that can bridge the gap between genetic variants and phenotypes. Over the past decade, we have witnessed a sharp rise of effort in the identification and application of eQTLs. The successful application of eQTLs relies heavily on their reproducibility. The current eQTL databases such as Genotype-Tissue Expression (GTEx) were populated primarily with eQTLs deriving from germline single nucleotide polymorphisms and normal tissue gene expression. The novel scenarios that employ eQTL models for prediction purposes often involve disease phenotypes characterized by altered gene expressions. To evaluate eQTL reproducibility across diverse data sources and the effect of disease-specific gene expression alteration on eQTL identification, we conducted an eQTL study using 5178 samples from The Cancer Genome Atlas (TCGA). We found that the reproducibility of eQTLs between normal and tumor tissues was low in terms of the number of shared eQTLs. However, among the shared eQTLs, the effect directions were generally concordant. This suggests that the source of the gene expression (normal or tumor tissue) has a strong effect on the detectable eQTLs and the effect direction of the eQTLs. Additional analyses demonstrated good directional concordance of eQTLs between GTEx and TCGA. Furthermore, we found that multi-tissue eQTLs may exert opposite effects across multiple tissue types. In summary, our results suggest that eQTL prediction models need to carefully address tissue and disease dependency of eQTLs. Tissue-disease-specific eQTL databases can afford more accurate prediction models for future studies.

MutEx: a multifaceted gateway for exploring integrative pan-cancer genomic data.

Somatic mutation and gene expression dysregulation are considered two major tumorigenesis factors. While independent investigations of either factor pervade, studies of associations between somatic mutations and gene expression changes have been sporadic and nonsystematic. Utilizing genomic data collected from 11 315 subjects of 33 distinct cancer types, we constructed MutEx, a pan-cancer integrative genomic database. This database records the relationships among gene expression, somatic mutation and survival data for cancer patients. MutEx can be used to swiftly explore the relationship between these genomic/clinic features within and across cancer types and, more importantly, search for corroborating evidence for hypothesis inception. Our database also incorporated Gene Ontology and several pathway databases to enhance functional annotation, and elastic net and a gene expression composite score to aid in survival analysis. To demonstrate the usability of MutEx, we provide several application examples, including top somatic mutations associated with the most extensive expression dysregulation in breast cancer, differential mutational burden downstream of DNA mismatch repair gene mutations and composite gene expression score-based survival difference in breast cancer. MutEx can be accessed at http://www.innovebioinfo.com/Databases/Mutationdb_About.php.

MetaGSCA: A tool for meta-analysis of gene set differential coexpression.

Analyses of gene set differential coexpression may shed light on molecular mechanisms underlying phenotypes and diseases. However, differential coexpression analyses of conceptually similar individual studies are often inconsistent and underpowered to provide definitive results. Researchers can greatly benefit from an open-source application facilitating the aggregation of evidence of differential coexpression across studies and the estimation of more robust common effects. We developed Meta Gene Set Coexpression Analysis (MetaGSCA), an analytical tool to systematically assess differential coexpression of an a priori defined gene set by aggregating evidence across studies to provide a definitive result. In the kernel, a nonparametric approach that accounts for the gene-gene correlation structure is used to test whether the gene set is differentially coexpressed between two comparative conditions, from which a permutation test p-statistic is computed for each individual study. A meta-analysis is then performed to combine individual study results with one of two options: a random-intercept logistic regression model or the inverse variance method. We demonstrated MetaGSCA in case studies investigating two human diseases and identified pathways highly relevant to each disease across studies. We further applied MetaGSCA in a pan-cancer analysis with hundreds of major cellular pathways in 11 cancer types. The results indicated that a majority of the pathways identified were dysregulated in the pan-cancer scenario, many of which have been previously reported in the cancer literature. Our analysis with randomly generated gene sets showed excellent specificity, indicating that the significant pathways/gene sets identified by MetaGSCA are unlikely false positives. MetaGSCA is a user-friendly tool implemented in both forms of a Web-based application and an R package "MetaGSCA". It enables comprehensive meta-analyses of gene set differential coexpression data, with an optional module of post hoc pathway crosstalk network analysis to identify and visualize pathways having similar coexpression profiles.

EditPredict: Prediction of RNA editable sites with convolutional neural network.

RNA editing exerts critical impacts on numerous biological processes. While millions of RNA editings have been identified in humans, much more are expected to be discovered. In this work, we constructed Convolutional Neural Network (CNN) models to predict human RNA editing events in both Alu regions and non-Alu regions. With a validation dataset resulting from CRISPR/Cas9 knockout of the ADAR1 enzyme, the validation accuracies reached 99.5% and 93.6% for Alu and non-Alu regions, respectively. We ported our CNN models in a web service named EditPredict. EditPredict not only works on reference genome sequences but can also take into consideration single nucleotide variants in personal genomes. In addition to the human genome, EditPredict tackles other model organisms including bumblebee, fruitfly, mouse, and squid genomes. EditPredict can be used stand-alone to predict novel RNA editing and it can be used to assist in filtering for candidate RNA editing detected from RNA-Seq data.

AnnoGen: annotating genome-wide pragmatic features.

MOTIVATION: Genome annotation is an important step for all in-depth bioinformatics analysis. It is imperative to augment quantity and diversity of genome-wide annotation data for the latest reference genome to promote its adoption by ongoing and future impactful studies. RESULTS: We developed a python toolkit AnnoGen, which at the first time, allows the annotation of three pragmatic genomic features for the GRCh38 genome in enormous base-wise quantities. The three features are chemical binding Energy, sequence information Entropy and Homology Score. The Homology Score is an exceptional feature that captures the genome-wide homology through single-base-offset tiling windows of 100 continual nucleotide bases. AnnoGen is capable of annotating the proprietary pragmatic features for variable user-interested genomic regions and optionally comparing two parallel sets of genomic regions. AnnoGen is characterized with simple utility modes and succinct HTML report of informative statistical tables and plots. AVAILABILITY AND IMPLEMENTATION: https://github.com/shengqh/annogen.

Non-canonical RNA-DNA differences and other human genomic features are enriched within very short tandem repeats.

Very short tandem repeats bear substantial genetic, evolutional, and pathological significance in genome analyses. Here, we compiled a census of tandem mono-nucleotide/di-nucleotide/tri-nucleotide repeats (MNRs/DNRs/TNRs) in GRCh38, which we term "polytracts" in general. Of the human genome, 144.4 million nucleotides (4.7%) are occupied by polytracts, and 0.47 million single nucleotides are identified as polytract hinges, i.e., break-points of tandem polytracts. Preliminary exploration of the census suggested polytract hinge sites and boundaries of AAC polytracts may bear a higher mapping error rate than other polytract regions. Further, we revealed landscapes of polytract enrichment with respect to nearly a hundred genomic features. We found MNRs, DNRs, and TNRs displayed noticeable difference in terms of locational enrichment for miscellaneous genomic features, especially RNA editing events. Non-canonical and C-to-U RNA-editing events are enriched inside and/or adjacent to MNRs, while all categories of RNA-editing events are under-represented in DNRs. A-to-I RNA-editing events are generally under-represented in polytracts. The selective enrichment of non-canonical RNA-editing events within MNR adjacency provides a negative evidence against their authenticity. To enable similar locational enrichment analyses in relation to polytracts, we developed a software Polytrap which can handle 11 reference genomes. Additionally, we compiled polytracts of four model organisms into a Track Hub which can be integrated into USCS Genome Browser as an official track for convenient visualization of polytracts.

Advancing Pan-cancer Gene Expression Survial Analysis by Inclusion of Non-coding RNA.

Non-coding RNAs occupy a significant fraction of the human genome. Their biological significance is backed up by a plethora of emerging evidence. One of the most robust approaches to demonstrate non-coding RNA's biological relevance is through their prognostic value. Using the rich gene expression data from The Cancer Genome Altas (TCGA), we designed Advanced Expression Survival Analysis (AESA), a web tool which provides several novel survival analysis approaches not offered by previous tools. In addition to the common single-gene approach, AESA computes the gene expression composite score of a set of genes for survival analysis and utilizes permutation test or cross-validation to assess the significance of log-rank statistic and the degree of over-fitting. AESA offers survival feature selection with post-selection inference and utilizes expanded TCGA clinical data including overall, disease-specific, disease-free, and progression-free survival information. Users can analyse either protein-coding or non-coding regions of the transcriptome. We demonstrated the effectiveness of AESA using several empirical examples. Our analyses showed that non-coding RNAs perform as well as messenger RNAs in predicting survival of cancer patients. These results reinforce the potential prognostic value of non-coding RNAs. AESA is developed as a module in the freely accessible analysis suite MutEx. Abbreviation: ACC: Adrenocortical Carcinoma (n = 92); BLCA: Bladder Urothelial Carcinoma (n = 412); BRCA: Breast Invasive Carcinoma (n = 1098); CESC: Cervical Squamous Cell Carcinoma and Endocervical Adenocarcinoma (n = 307); CHOL: Cholangiocarcinoma (n = 51); COAD: Colon Adenocarcinoma (n = 461); DLBC: Lymphoid Neoplasm Diffuse Large B-cell Lymphoma (n = 58); ESCA: Oesophageal Carcinoma (n = 185); GBM: Glioblastoma Multiforme (n = 617); HNSC: Head and Neck Squamous Cell Carcinoma (n = 528); KICH: Kidney Chromophobe (n = 113); KIRC: Kidney Renal Clear Cell Carcinoma (n = 537); KIRP: Kidney Renal Papillary Cell Carcinoma (n = 291); LAML: Acute Myeloid Leukaemia (n = 200); LGG: Brain Lower Grade Glioma (n = 516); LIHC: Liver Hepatocellular Carcinoma (n = 377); LUAD: Lung Adenocarcinoma (n = 585); LUSC: Lung Squamous Cell Carcinoma (n = 504); MESO: Mesothelioma (n = 87); OV: Ovarian Serous Cystadenocarcinoma (n = 608) PAAD: Pancreatic Adenocarcinoma (n = 185); PCPG: Pheochromocytoma and Paraganglioma (n = 179); PRAD: Prostate Adenocarcinoma (n = 500); READ: Rectum Adenocarcinoma (n = 172); SARC: Sarcoma (n = 261); SKCM: Skin Cutaneous Melanoma (n = 470); STAD: Stomach Adenocarcinoma (n = 443); TGCT: Testicular Germ Cell Tumours (n = 150); THCA: Thyroid Carcinoma (n = 507) THYM: Thymoma (n = 124); UCEC: Uterine Corpus Endometrial Carcinoma (n = 560); UCS: Uterine Carcinosarcoma (n = 57); UVM: Uveal Melanoma (n = 80).

Epigenetic silencing of SOCS5 potentiates JAK-STAT signaling and progression of T-cell acute lymphoblastic leukemia.

Activating mutations in cytokine receptors and transcriptional regulators govern aberrant signal transduction in T-cell lineage acute lymphoblastic leukemia (T-ALL). However, the roles played by suppressors of cytokine signaling remain incompletely understood. We examined the regulatory roles of suppressor of cytokine signaling 5 (SOCS5) in T-ALL cellular signaling networks and leukemia progression. We found that SOCS5 was differentially expressed in primary T-ALL and its expression levels were lowered in HOXA-deregulated leukemia harboring KMT2A gene rearrangements. Here, we report that SOCS5 expression is epigenetically regulated by DNA methyltransferase-3A-mediated DNA methylation and methyl CpG binding protein-2-mediated histone deacetylation. We show that SOCS5 negatively regulates T-ALL cell growth and cell cycle progression but has no effect on apoptotic cell death. Mechanistically, SOCS5 silencing induces activation of JAK-STAT signaling, and negatively regulates interleukin-7 and interleukin-4 receptors. Using a human T-ALL murine xenograft model, we show that genetic inactivation of SOCS5 accelerates leukemia engraftment and progression, and leukemia burden. We postulate that SOCS5 is epigenetically deregulated in T-ALL and serves as an important regulator of T-ALL cell proliferation and leukemic progression. Our results link aberrant downregulation of SOCS5 expression to the enhanced activation of the JAK-STAT and cytokine receptor-signaling cascade in T-ALL.

Distinct histone methylation and transcription profiles are established during the development of cellular quiescence in yeast.

BACKGROUND: Quiescent cells have a low level of gene activity compared to growing cells. Using a yeast model for cellular quiescence, we defined the genome-wide profiles of three species of histone methylation associated with active transcription between growing and quiescent cells, and correlated these profiles with the presence of RNA polymerase II and transcripts. RESULTS: Quiescent cells retained histone methylations normally associated with transcriptionally active chromatin and had many transcripts in common with growing cells. Quiescent cells also contained significant levels of RNA polymerase II, but only low levels of the canonical initiating and elongating forms of the polymerase. The RNA polymerase II associated with genes in quiescent cells displayed a distinct occupancy profile compared to its pattern of occupancy across genes in actively growing cells. Although transcription is generally repressed in quiescent cells, analysis of individual genes identified a period of active transcription during the development of quiescence. CONCLUSIONS: The data suggest that the transcript profile and histone methylation marks in quiescent cells were established both in growing cells and during the development of quiescence and then retained in these cells. Together, this might ensure that quiescent cells can rapidly adapt to a changing environment to resume growth.

Integration of ruxolitinib into dose-intensified therapy targeted against a novel JAK2 F694L mutation in B-precursor acute lymphoblastic leukemia.

A 17-year-old girl with B-cell precursor acute lymphoblastic leukemia (BCP-ALL) with persistent minimal residual disease (MRD) who underwent standard chemotherapy was found to have a BCR-ABL1-like gene expression pattern. Genome sequencing revealed a JAK2 mutation not previously described in BCP-ALL and a potential therapeutic target. Due to concern for an on-therapy relapse, the JAK2 inhibitor ruxolitinib was incorporated into a modified chemotherapy backbone to achieve complete remission prior to stem cell transplant. Treatment was well tolerated and she had undetectable MRD prior to a matched allogeneic stem cell transplant and remained in remission at day +100.

RNA editing events in mitochondrial genes by ultra-deep sequencing methods: a comparison of cytoplasmic male sterile, fertile and restored genotypes in cotton.

Cytoplasmic male sterility (CMS) is a maternally inherited trait resulting in failure to produce functional pollen and is widely used in the production of hybrid seed. Improper RNA editing is implicated as the molecular basis for some CMS systems. However, the mechanism of CMS in cotton is unknown. This study compared RNA editing events in eight mitochondrial genes (atp1, 4, 6, 8, 9, and cox1, 2, 3) among three lines (maintainer B, CMS A, and restorer R). These events were quantified by ultra-deep sequencing of mitochondrial transcripts and sequencing of cloned versions of these genes as cDNAs. A comparison of genomic PCR and RT-PCR products detected 72 editing sites in coding sequences in the eight genes and four partial editing sites in the 3'-untranslated region of atp6. The most frequent alteration (61.4 %) resulted in changes of hydrophilic amino acids to hydrophobic amino acids and the most common alteration was proline (P) to leucine (L) (26.7 %). In atp6, RNA editing created a stop codon from a glutamine in the genomic sequence. Statistical analysis of the frequencies of RNA editing events detected differences between mtDNA genes, but no differences between cotton cytoplasms that could account for the CMS phenotype or restoration. This study represents the first work to use next-generation sequencing to identify RNA editing positions and efficiency, and possible association with CMS and restoration in plants.

The inflammatory response of human pancreatic cancer samples compared to normal controls.

Pancreatic ductal adenocarcinoma (PDAC) is a poor prognosis cancer with an aggressive growth profile that is often diagnosed at late stage and that has few curative or therapeutic options. PDAC growth has been linked to alterations in the pancreas microbiome, which could include the presence of the fungus Malassezia. We used RNA-sequencing to compare 14 matched tumor and normal (tumor adjacent) pancreatic cancer samples and found Malassezia RNA in both the PDAC and normal tissues. Although the presence of Malassezia was not correlated with tumor growth, a set of immune- and inflammatory-related genes were up-regulated in the PDAC compared to the normal samples, suggesting that they are involved in tumor progression. Gene set enrichment analysis suggests that activation of the complement cascade pathway and inflammation could be involved in pro PDAC growth.

Humanized Patient-derived Xenograft Models of Disseminated Ovarian Cancer Recapitulate Key Aspects of the Tumor Immune Environment within the Peritoneal Cavity.

The importance of the immune microenvironment in ovarian cancer progression, metastasis, and response to therapies has become increasingly clear, especially with the new emphasis on immunotherapies. To leverage the power of patient-derived xenograft (PDX) models within a humanized immune microenvironment, three ovarian cancer PDXs were grown in humanized NBSGW (huNBSGW) mice engrafted with human CD34+ cord blood-derived hematopoietic stem cells. Analysis of cytokine levels in the ascites fluid and identification of infiltrating immune cells in the tumors demonstrated that these humanized PDX (huPDX) established an immune tumor microenvironment similar to what has been reported for patients with ovarian cancer. The lack of human myeloid cell differentiation has been a major setback for humanized mouse models, but our analysis shows that PDX engraftment increases the human myeloid population in the peripheral blood. Analysis of cytokines within the ascites fluid of huPDX revealed high levels of human M-CSF, a key myeloid differentiation factor as well as other elevated cytokines that have previously been identified in ovarian cancer patient ascites fluid including those involved in immune cell differentiation and recruitment. Human tumor-associated macrophages and tumor-infiltrating lymphocytes were detected within the tumors of humanized mice, demonstrating immune cell recruitment to tumors. Comparison of the three huPDX revealed certain differences in cytokine signatures and in the extent of immune cell recruitment. Our studies show that huNBSGW PDX models reconstitute important aspects of the ovarian cancer immune tumor microenvironment, which may recommend these models for preclinical therapeutic trials. Significance: huPDX models are ideal preclinical models for testing novel therapies. They reflect the genetic heterogeneity of the patient population, enhance human myeloid differentiation, and recruit immune cells to the tumor microenvironment.

Dominant Gene Expression Profiles Define Adenoid Cystic Carcinoma (ACC) from Different Tissues: Validation of a Gene Signature Classifier for Poor Survival in Salivary Gland ACC.

Adenoid cystic carcinoma (ACC) is an aggressive malignancy that most often arises in salivary or lacrimal glands but can also occur in other tissues. We used optimized RNA-sequencing to analyze the transcriptomes of 113 ACC tumor samples from salivary gland, lacrimal gland, breast or skin. ACC tumors from different organs displayed remarkedly similar transcription profiles, and most harbored translocations in the MYB or MYBL1 genes, which encode oncogenic transcription factors that may induce dramatic genetic and epigenetic changes leading to a dominant 'ACC phenotype'. Further analysis of the 56 salivary gland ACC tumors led to the identification of three distinct groups of patients, based on gene expression profiles, including one group with worse survival. We tested whether this new cohort could be used to validate a biomarker developed previously with a different set of 68 ACC tumor samples. Indeed, a 49-gene classifier developed with the earlier cohort correctly identified 98% of the poor survival patients from the new set, and a 14-gene classifier was almost as accurate. These validated biomarkers form a platform to identify and stratify high-risk ACC patients into clinical trials of targeted therapies for sustained clinical response.

Adult human CD133/1(+) kidney cells isolated from papilla integrate into developing kidney tubules.

Approximately 60,000 patients in the United States are waiting for a kidney transplant due to genetic, immunologic and environmentally caused kidney failure. Adult human renal stem cells could offer opportunities for autologous transplant and repair of damaged organs. Current data suggest that there are multiple progenitor types in the kidney with distinct localizations. In the present study, we characterize cells derived from human kidney papilla and show their capacity for tubulogenesis. In situ, nestin(+) and CD133/1(+) cells were found extensively intercalated between tubular epithelia in the loops of Henle of renal papilla, but not of the cortex. Populations of primary cells from the renal cortex and renal papilla were isolated by enzymatic digestion from human kidneys unsuited for transplant and immuno-enriched for CD133/1(+) cells. Isolated CD133/1(+) papillary cells were positive for nestin, as well as several human embryonic stem cell markers (SSEA4, Nanog, SOX2, and OCT4/POU5F1) and could be triggered to adopt tubular epithelial and neuronal-like phenotypes. Isolated papillary cells exhibited morphologic plasticity upon modulation of culture conditions and inhibition of asymmetric cell division. Labeled papillary cells readily associated with cortical tubular epithelia in co-culture and 3-dimensional collagen gel cultures. Heterologous organ culture demonstrated that CD133/1(+) progenitors from the papilla and cortex became integrated into developing kidney tubules. Tubular epithelia did not participate in tubulogenesis. Human renal papilla harbor cells with the hallmarks of adult kidney stem/progenitor cells that can be amplified and phenotypically modulated in culture while retaining the capacity to form new kidney tubules. This article is part of a Special Issue entitled: Polycystic Kidney Disease.