Terminal modification, sequence, length, and PIWI-protein identity determine piRNA stability.

In animals, PIWI-interacting RNAs (piRNAs) silence transposons, fight viral infections, and regulate gene expression. piRNA biogenesis concludes with 3' terminal trimming and 2'-O-methylation. Both trimming and methylation influence piRNA stability. Our biochemical data show that multiple mechanisms destabilize unmethylated mouse piRNAs, depending on whether the piRNA 5' or 3' sequence is complementary to a trigger RNA. Unlike target-directed degradation of microRNAs, complementarity-dependent destabilization of piRNAs in mice and flies is blocked by 3' terminal 2'-O-methylation and does not require base pairing to both the piRNA seed and the 3' sequence. In flies, 2'-O-methylation also protects small interfering RNAs (siRNAs) from complementarity-dependent destruction. By contrast, pre-piRNA trimming protects mouse piRNAs from a degradation pathway unaffected by trigger complementarity. In testis lysate and in vivo, internal or 3' terminal uridine- or guanine-rich tracts accelerate pre-piRNA decay. Loss of both trimming and 2'-O-methylation causes the mouse piRNA pathway to collapse, demonstrating that these modifications collaborate to stabilize piRNAs.

GTSF1 accelerates target RNA cleavage by PIWI-clade Argonaute proteins.

Argonaute proteins use nucleic acid guides to find and bind specific DNA or RNA target sequences. Argonaute proteins have diverse biological functions and many retain their ancestral endoribonuclease activity, cleaving the phosphodiester bond between target nucleotides t10 and t11. In animals, the PIWI proteins-a specialized class of Argonaute proteins-use 21-35 nucleotide PIWI-interacting RNAs (piRNAs) to direct transposon silencing, protect the germline genome, and regulate gene expression during gametogenesis1. The piRNA pathway is required for fertility in one or both sexes of nearly all animals. Both piRNA production and function require RNA cleavage catalysed by PIWI proteins. Spermatogenesis in mice and other placental mammals requires three distinct, developmentally regulated PIWI proteins: MIWI (PIWIL1), MILI (PIWIL2) and MIWI22-4 (PIWIL4). The piRNA-guided endoribonuclease activities of MIWI and MILI are essential for the production of functional sperm5,6. piRNA-directed silencing in mice and insects also requires GTSF1, a PIWI-associated protein of unknown function7-12. Here we report that GTSF1 potentiates the weak, intrinsic, piRNA-directed RNA cleavage activities of PIWI proteins, transforming them into efficient endoribonucleases. GTSF1 is thus an example of an auxiliary protein that potentiates the catalytic activity of an Argonaute protein.

PIWI-interacting RNAs: small RNAs with big functions.

In animals, PIWI-interacting RNAs (piRNAs) of 21-35 nucleotides in length silence transposable elements, regulate gene expression and fight viral infection. piRNAs guide PIWI proteins to cleave target RNA, promote heterochromatin assembly and methylate DNA. The architecture of the piRNA pathway allows it both to provide adaptive, sequence-based immunity to rapidly evolving viruses and transposons and to regulate conserved host genes. piRNAs silence transposons in the germ line of most animals, whereas somatic piRNA functions have been lost, gained and lost again across evolution. Moreover, most piRNA pathway proteins are deeply conserved, but different animals employ remarkably divergent strategies to produce piRNA precursor transcripts. Here, we discuss how a common piRNA pathway allows animals to recognize diverse targets, ranging from selfish genetic elements to genes essential for gametogenesis.

A-MYB/TCFL5 regulatory architecture ensures the production of pachytene piRNAs in placental mammals.

In male mice, the transcription factor A MYB initiates the transcription of pachytene piRNA genes during meiosis. Here, we report that A MYB activates the transcription factor Tcfl5 produced in pachytene spermatocytes. Subsequently, A MYB and TCFL5 reciprocally reinforce their own transcription to establish a positive feedback circuit that triggers pachytene piRNA production. TCFL5 regulates the expression of genes required for piRNA maturation and promotes transcription of evolutionarily young pachytene piRNA genes, whereas A-MYB activates the transcription of older pachytene piRNA genes. Intriguingly, pachytene piRNAs from TCFL5-dependent young loci initiates the production of piRNAs from A-MYB-dependent older loci ensuring the self-propagation of pachytene piRNAs. A MYB and TCFL5 act via a set of incoherent feedforward loops that drive regulation of gene expression by pachytene piRNAs during spermatogenesis. This regulatory architecture is conserved in rhesus macaque, suggesting that it was present in the last common ancestor of placental mammals.

CRISPR-induced exon skipping of ß-catenin reveals tumorigenic mutants driving distinct subtypes of liver cancer.

CRISPR/Cas9-driven cancer modeling studies are based on the disruption of tumor suppressor genes by small insertions or deletions (indels) that lead to frame-shift mutations. In addition, CRISPR/Cas9 is widely used to define the significance of cancer oncogenes and genetic dependencies in loss-of-function studies. However, how CRISPR/Cas9 influences gain-of-function oncogenic mutations is elusive. Here, we demonstrate that single guide RNA targeting exon 3 of Ctnnb1 (encoding ß-catenin) results in exon skipping and generates gain-of-function isoforms in vivo. CRISPR/Cas9-mediated exon skipping of Ctnnb1 induces liver tumor formation in synergy with YAPS127A in mice. We define two distinct exon skipping-induced tumor subtypes with different histological and transcriptional features. Notably, ectopic expression of two exon-skipped ß-catenin transcript isoforms together with YAPS127A phenocopies the two distinct subtypes of liver cancer. Moreover, we identify similar CTNNB1 exon-skipping events in patients with hepatocellular carcinoma. Collectively, our findings advance our understanding of ß-catenin-related tumorigenesis and reveal that CRISPR/Cas9 can be repurposed, in vivo, to study gain-of-function mutations of oncogenes in cancer. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

The MYBL1/TCFL5 transcription network: two collaborative factors with central role in male meiosis.

Male gametogenesis, spermatogenesis, is a stepwise developmental process to generate mature sperm. The most intricate process of spermatogenesis is meiosis during which two successive cell divisions ensue with dramatic cellular and molecular changes to produce haploid cells. After entry into meiosis, several forms of regulatory events control the orderly progression of meiosis and the timely entry into post-meiotic sperm differentiation. Among other mechanisms, changes to gene expression are controlled by key transcription factors. In this review, we will discuss the gene regulatory mechanisms underlying meiotic entry, meiotic progression, and post-meiotic differentiation with a particular emphasis on the MYBL1/TCFL5 regulatory architecture and how this architecture involves in various forms of transcription network motifs to regulate gene expression.

The transcription factor TCFL5 responds to A-MYB to elaborate the male meiotic program in mice.

In brief: The testis-specific transcription factor, TCFL5, expressed in pachytene spermatocytes regulates the meiotic gene expression program in collaboration with the transcription factor A-MYB. Abstract: In male mice, the transcription factors STRA8 and MEISON initiate meiosis I. We report that STRA8/MEISON activates the transcription factors A-MYB and TCFL5, which together reprogram gene expression after spermatogonia enter into meiosis. TCFL5 promotes the transcription of genes required for meiosis, mRNA turnover, miR-34/449 production, meiotic exit, and spermiogenesis. This transcriptional architecture is conserved in rhesus macaque, suggesting TCFL5 plays a central role in meiosis and spermiogenesis in placental mammals. Tcfl5em1/em1 mutants are sterile, and spermatogenesis arrests at the mid- or late-pachytene stage of meiosis. Moreover, Tcfl5+/em1 mutants produce fewer motile sperm.

YAP1 Withdrawal in Hepatoblastoma Drives Therapeutic Differentiation of Tumor Cells to Functional Hepatocyte-Like Cells.

BACKGROUND AND AIMS: Despite surgical and chemotherapeutic advances, the 5-year survival rate for stage IV hepatoblastoma (HB), the predominant pediatric liver tumor, remains at 27%. Yes-associated protein 1 (YAP1) and ß-catenin co-activation occurs in 80% of children's HB; however, a lack of conditional genetic models precludes tumor maintenance exploration. Thus, the need for a targeted therapy remains unmet. Given the predominance of YAP1 and ß-catenin activation in HB, we sought to evaluate YAP1 as a therapeutic target in HB. APPROACH AND RESULTS: We engineered the conditional HB murine model using hydrodynamic injection to deliver transposon plasmids encoding inducible YAP1S127A , constitutive ß-cateninDelN90 , and a luciferase reporter to murine liver. Tumor regression was evaluated using bioluminescent imaging, tumor landscape characterized using RNA and ATAC sequencing, and DNA footprinting. Here we show that YAP1S127A withdrawal mediates more than 90% tumor regression with survival for 230+ days in mice. YAP1S127A withdrawal promotes apoptosis in a subset of tumor cells, and in remaining cells induces a cell fate switch that drives therapeutic differentiation of HB tumors into Ki-67-negative hepatocyte-like HB cells ("HbHeps") with hepatocyte-like morphology and mature hepatocyte gene expression. YAP1S127A withdrawal drives the formation of hbHeps by modulating liver differentiation transcription factor occupancy. Indeed, tumor-derived hbHeps, consistent with their reprogrammed transcriptional landscape, regain partial hepatocyte function and rescue liver damage in mice. CONCLUSIONS: YAP1S127A withdrawal, without silencing oncogenic ß-catenin, significantly regresses hepatoblastoma, providing in vivo data to support YAP1 as a therapeutic target for HB. YAP1S127A withdrawal alone sufficiently drives long-term regression in HB, as it promotes cell death in a subset of tumor cells and modulates transcription factor occupancy to reverse the fate of residual tumor cells to mimic functional hepatocytes.

Depletion of TRRAP Induces p53-Independent Senescence in Liver Cancer by Down-Regulating Mitotic Genes.

Hepatocellular carcinoma (HCC) is an aggressive subtype of liver cancer with few effective treatments, and the underlying mechanisms that drive HCC pathogenesis remain poorly characterized. Identifying genes and pathways essential for HCC cell growth will aid the development of new targeted therapies for HCC. Using a kinome CRISPR screen in three human HCC cell lines, we identified transformation/transcription domain-associated protein (TRRAP) as an essential gene for HCC cell proliferation. TRRAP has been implicated in oncogenic transformation, but how it functions in cancer cell proliferation is not established. Here, we show that depletion of TRRAP or its co-factor, histone acetyltransferase KAT5, inhibits HCC cell growth through induction of p53-independent and p21-independent senescence. Integrated cancer genomics analyses using patient data and RNA sequencing identified mitotic genes as key TRRAP/KAT5 targets in HCC, and subsequent cell cycle analyses revealed that TRRAP-depleted and KAT5-depleted cells are arrested at the G2/M phase. Depletion of topoisomerase II alpha (TOP2A), a mitotic gene and TRRAP/KAT5 target, was sufficient to recapitulate the senescent phenotype of TRRAP/KAT5 knockdown. Conclusion: Our results uncover a role for TRRAP/KAT5 in promoting HCC cell proliferation by activating mitotic genes. Targeting the TRRAP/KAT5 complex is a potential therapeutic strategy for HCC.

Genetic disruption of oncogenic Kras sensitizes lung cancer cells to Fas receptor-mediated apoptosis.

Genetic lesions that activate KRAS account for ∼30% of the 1.6 million annual cases of lung cancer. Despite clinical need, KRAS is still undruggable using traditional small-molecule drugs/inhibitors. When oncogenic Kras is suppressed by RNA interference, tumors initially regress but eventually recur and proliferate despite suppression of Kras Here, we show that tumor cells can survive knockout of oncogenic Kras, indicating the existence of Kras-independent survival pathways. Thus, even if clinical KRAS inhibitors were available, resistance would remain an obstacle to treatment. Kras-independent cancer cells exhibit decreased colony formation in vitro but retain the ability to form tumors in mice. Comparing the transcriptomes of oncogenic Kras cells and Kras knockout cells, we identified 603 genes that were specifically up-regulated in Kras knockout cells, including the Fas gene, which encodes a cell surface death receptor involved in physiological regulation of apoptosis. Antibodies recognizing Fas receptor efficiently induced apoptosis of Kras knockout cells but not oncogenic Kras-expressing cells. Increased Fas expression in Kras knockout cells was attributed to decreased association of repressive epigenetic marks at the Fas promoter. Concordant with this observation, treating oncogenic Kras cells with histone deacetylase inhibitor and Fas-activating antibody efficiently induced apoptosis, thus bypassing the need to inhibit Kras. Our results suggest that activation of Fas could be exploited as an Achilles' heel in tumors initiated by oncogenic Kras.

DEBrowser: interactive differential expression analysis and visualization tool for count data.

BACKGROUND: Sequencing data has become a standard measure of diverse cellular activities. For example, gene expression is accurately measured by RNA sequencing (RNA-Seq) libraries, protein-DNA interactions are captured by chromatin immunoprecipitation sequencing (ChIP-Seq), protein-RNA interactions by crosslinking immunoprecipitation sequencing (CLIP-Seq) or RNA immunoprecipitation (RIP-Seq) sequencing, DNA accessibility by assay for transposase-accessible chromatin (ATAC-Seq), DNase or MNase sequencing libraries. The processing of these sequencing techniques involves library-specific approaches. However, in all cases, once the sequencing libraries are processed, the result is a count table specifying the estimated number of reads originating from each genomic locus. Differential analysis to determine which loci have different cellular activity under different conditions starts with the count table and iterates through a cycle of data assessment, preparation and analysis. Such complex analysis often relies on multiple programs and is therefore a challenge for those without programming skills. RESULTS: We developed DEBrowser as an R bioconductor project to interactively visualize every step of the differential analysis, without programming. The application provides a rich and interactive web based graphical user interface built on R's shiny infrastructure. DEBrowser allows users to visualize data with various types of graphs that can be explored further by selecting and re-plotting any desired subset of data. Using the visualization approaches provided, users can determine and correct technical variations such as batch effects and sequencing depth that affect differential analysis. We show DEBrowser's ease of use by reproducing the analysis of two previously published data sets. CONCLUSIONS: DEBrowser is a flexible, intuitive, web-based analysis platform that enables an iterative and interactive analysis of count data without any requirement of programming knowledge.

Hypermigration of macrophages through the concerted action of GRA effectors on NF-κB/p38 signaling and host chromatin accessibility potentiates Toxoplasma dissemination.

Mononuclear phagocytes facilitate the dissemination of the obligate intracellular parasite Toxoplasma gondii. Here, we report how a set of secreted parasite effector proteins from dense granule organelles (GRA) orchestrates dendritic cell-like chemotactic and pro-inflammatory activation of parasitized macrophages. These effects enabled efficient dissemination of the type II T. gondii lineage, a highly prevalent genotype in humans. We identify novel functions for effectors GRA15 and GRA24 in promoting CCR7-mediated macrophage chemotaxis by acting on NF-κB and p38 MAPK signaling pathways, respectively, with contributions of GRA16/18 and counter-regulation by effector TEEGR. Further, GRA28 boosted chromatin accessibility and GRA15/24/NF-κB-dependent transcription at the Ccr7 gene locus in primary macrophages. In vivo, adoptively transferred macrophages infected with wild-type T. gondii outcompeted macrophages infected with a GRA15/24 double mutant in migrating to secondary organs in mice. The data show that T. gondii, rather than being passively shuttled, actively promotes its dissemination by inducing a finely regulated pro-migratory state in parasitized human and murine phagocytes via co-operating polymorphic GRA effectors.

Hypermigration of macrophages through the concerted action of GRA effectors on NF-κB/p38 signaling and host chromatin accessibility potentiates Toxoplasma dissemination.

Mononuclear phagocytes facilitate the dissemination of the obligate intracellular parasite Toxoplasma gondii. Here, we report how a set of secreted parasite effector proteins from dense granule organelles (GRA) orchestrates dendritic cell-like chemotactic and pro-inflammatory activation of parasitized macrophages. These effects enabled efficient dissemination of the type II T. gondii lineage, a highly prevalent genotype in humans. We identify novel functions for effectors GRA15 and GRA24 in promoting CCR7-mediated macrophage chemotaxis by acting on NF-κB and p38 mitogen-activated protein kinase signaling pathways, respectively, with contributions by GRA16/18 and counter-regulation by effector TEEGR. Furthermore, GRA28 boosted chromatin accessibility and GRA15/24/NF-κB-dependent transcription at the Ccr7 gene locus in primary macrophages. In vivo, adoptively transferred macrophages infected with wild-type T. gondii outcompeted macrophages infected with a GRA15/24 double mutant in migrating to secondary organs in mice. The data show that T. gondii, rather than being passively shuttled, actively promotes its dissemination by inducing a finely regulated pro-migratory state in parasitized human and murine phagocytes via co-operating polymorphic GRA effectors. IMPORTANCE: Intracellular pathogens can hijack the cellular functions of infected host cells to their advantage, for example, for intracellular survival and dissemination. However, how microbes orchestrate the hijacking of complex cellular processes, such as host cell migration, remains poorly understood. As such, the common parasite Toxoplasma gondii actively invades the immune cells of humans and other vertebrates and modifies their migratory properties. Here, we show that the concerted action of a number of secreted effector proteins from the parasite, principally GRA15 and GRA24, acts on host cell signaling pathways to activate chemotaxis. Furthermore, the protein effector GRA28 selectively acted on chromatin accessibility in the host cell nucleus to selectively boost host gene expression. The joint activities of GRA effectors culminated in pro-migratory signaling within the infected phagocyte. We provide a molecular framework delineating how T. gondii can orchestrate a complex biological phenotype, such as the migratory activation of phagocytes to boost dissemination.

Long first exons and epigenetic marks distinguish conserved pachytene piRNA clusters from other mammalian genes.

In the male germ cells of placental mammals, 26-30-nt-long PIWI-interacting RNAs (piRNAs) emerge when spermatocytes enter the pachytene phase of meiosis. In mice, pachytene piRNAs derive from ~100 discrete autosomal loci that produce canonical RNA polymerase II transcripts. These piRNA clusters bear 5' caps and 3' poly(A) tails, and often contain introns that are removed before nuclear export and processing into piRNAs. What marks pachytene piRNA clusters to produce piRNAs, and what confines their expression to the germline? We report that an unusually long first exon (≥ 10 kb) or a long, unspliced transcript correlates with germline-specific transcription and piRNA production. Our integrative analysis of transcriptome, piRNA, and epigenome datasets across multiple species reveals that a long first exon is an evolutionarily conserved feature of pachytene piRNA clusters. Furthermore, a highly methylated promoter, often containing a low or intermediate level of CG dinucleotides, correlates with germline expression and somatic silencing of pachytene piRNA clusters. Pachytene piRNA precursor transcripts bind THOC1 and THOC2, THO complex subunits known to promote transcriptional elongation and mRNA nuclear export. Together, these features may explain why the major sources of pachytene piRNA clusters specifically generate these unique small RNAs in the male germline of placental mammals.

Dietary suppression of MHC class II expression in intestinal epithelial cells enhances intestinal tumorigenesis.

Little is known about how interactions of diet, intestinal stem cells (ISCs), and immune cells affect early-stage intestinal tumorigenesis. We show that a high-fat diet (HFD) reduces the expression of the major histocompatibility complex class II (MHC class II) genes in intestinal epithelial cells, including ISCs. This decline in epithelial MHC class II expression in a HFD correlates with reduced intestinal microbiome diversity. Microbial community transfer experiments suggest that epithelial MHC class II expression is regulated by intestinal flora. Mechanistically, pattern recognition receptor (PRR) and interferon-gamma (IFNγ) signaling regulates epithelial MHC class II expression. MHC class II-negative (MHC-II-) ISCs exhibit greater tumor-initiating capacity than their MHC class II-positive (MHC-II+) counterparts upon loss of the tumor suppressor Apc coupled with a HFD, suggesting a role for epithelial MHC class II-mediated immune surveillance in suppressing tumorigenesis. ISC-specific genetic ablation of MHC class II increases tumor burden cell autonomously. Thus, HFD perturbs a microbiome-stem cell-immune cell interaction that contributes to tumor initiation in the intestine.

Evolutionarily conserved pachytene piRNA loci are highly divergent among modern humans.

In the fetal mouse testis, PIWI-interacting RNAs (piRNAs) guide PIWI proteins to silence transposons but, after birth, most post-pubertal pachytene piRNAs map to the genome uniquely and are thought to regulate genes required for male fertility. In the human male, the developmental classes, precise genomic origins and transcriptional regulation of postnatal piRNAs remain undefined. Here, we demarcate the genes and transcripts that produce postnatal piRNAs in human juvenile and adult testes. As in the mouse, human A-MYB drives transcription of both pachytene piRNA precursor transcripts and messenger RNAs encoding piRNA biogenesis factors. Although human piRNA genes are syntenic to those in other placental mammals, their sequences are poorly conserved. In fact, pachytene piRNA loci are rapidly diverging even among modern humans. Our findings suggest that, during mammalian evolution, pachytene piRNA genes are under few selective constraints. We speculate that pachytene piRNA diversity may provide a hitherto unrecognized driver of reproductive isolation.

The evolutionarily conserved piRNA-producing locus pi6 is required for male mouse fertility.

Pachytene PIWI-interacting RNAs (piRNAs), which comprise >80% of small RNAs in the adult mouse testis, have been proposed to bind and regulate target RNAs like microRNAs, cleave targets like short interfering RNAs or lack biological function altogether. Although piRNA pathway protein mutants are male sterile, no biological function has been identified for any mammalian piRNA-producing locus. Here, we report that males lacking piRNAs from a conserved mouse pachytene piRNA locus on chromosome 6 (pi6) produce sperm with defects in capacitation and egg fertilization. Moreover, heterozygous embryos sired by pi6-/- fathers show reduced viability in utero. Molecular analyses suggest that pi6 piRNAs repress gene expression by cleaving messenger RNAs encoding proteins required for sperm function. pi6 also participates in a network of piRNA-piRNA precursor interactions that initiate piRNA production from a second piRNA locus on chromosome 10, as well as pi6 itself. Our data establish a direct role for pachytene piRNAs in spermiogenesis and embryo viability.

CRISPR-SONIC: targeted somatic oncogene knock-in enables rapid in vivo cancer modeling.

CRISPR/Cas9 has revolutionized cancer mouse models. Although loss-of-function genetics by CRISPR/Cas9 is well-established, generating gain-of-function alleles in somatic cancer models is still challenging because of the low efficiency of gene knock-in. Here we developed CRISPR-based Somatic Oncogene kNock-In for Cancer Modeling (CRISPR-SONIC), a method for rapid in vivo cancer modeling using homology-independent repair to integrate oncogenes at a targeted genomic locus. Using a dual guide RNA strategy, we integrated a plasmid donor in the 3'-UTR of mouse ß-actin, allowing co-expression of reporter genes or oncogenes from the ß-actin promoter. We showed that knock-in of oncogenic Ras and loss of p53 efficiently induced intrahepatic cholangiocarcinoma in mice. Further, our strategy can generate bioluminescent liver cancer to facilitate tumor imaging. This method simplifies in vivo gain-of-function genetics by facilitating targeted integration of oncogenes.

CRISPR/Cas9-mediated genome editing induces exon skipping by alternative splicing or exon deletion.

CRISPR is widely used to disrupt gene function by inducing small insertions and deletions. Here, we show that some single-guide RNAs (sgRNAs) can induce exon skipping or large genomic deletions that delete exons. For example, CRISPR-mediated editing of ß-catenin exon 3, which encodes an autoinhibitory domain, induces partial skipping of the in-frame exon and nuclear accumulation of ß-catenin. A single sgRNA can induce small insertions or deletions that partially alter splicing or unexpected larger deletions that remove exons. Exon skipping adds to the unexpected outcomes that must be accounted for, and perhaps taken advantage of, in CRISPR experiments.

MicroRNA expression patterns associated with hyperfunctioning and non-hyperfunctioning phenotypes in adrenocortical adenomas.

BACKGROUND: The adrenocortical adenoma (ACA) entity includes aldosterone-producing adenoma (APA), cortisol-producing adenoma (CPA), and non-hyperfunctioning adenoma (NHFA) phenotypes. While gene mutations and mRNA expression profiles have been partly characterized, less is known about the alterations involving microRNA (miRNA) expression. AIM: To characterize miRNA expression profile in relation to the subtypes of ACAs. SUBJECTS AND METHODS: miRNA expression profiles were determined in 26 ACAs (nine APAs, ten CPAs, and seven NHFAs) and four adrenal references using microarray-based screening. Significance analysis of microarrays (SAM) was carried out to identify differentially expressed miRNAs between ACA and adrenal cortices or between tumor subtypes. Selected differentially expressed miRNAs were validated in an extended series of 43 ACAs and ten adrenal references by quantitative RT-PCR. RESULTS: An hierarchical clustering revealed separate clusters for APAs and CPAs, while the NHFAs were found spread out within the APA/CPA clusters. When NHFA was excluded, the clustering analysis showed a better separation between APA and CPA. SAM analysis identified 40 over-expressed and three under-expressed miRNAs in the adenomas as compared with adrenal references. Fourteen miRNAs were common among the three ACA subtypes. Furthermore, we found specific miRNAs associated with different tumor phenotypes. CONCLUSION: The results suggest that miRNA expression profiles can distinguish different subtypes of ACA, which may contribute to a deeper understanding of ACA development and potential therapeutics.