A Matched Molecular and Clinical Analysis of the Epithelioid Haemangioendothelioma Cohort in the Stafford Fox Rare Cancer Program and Contextual Literature Review.

BACKGROUND: Epithelioid haemangioendothelioma (EHE) is an ultra-rare malignant vascular tumour with a prevalence of 1 per 1,000,000. It is typically molecularly characterised by a WWTR1::CAMTA1 gene fusion in approximately 90% of cases, or a YAP1::TFE3 gene fusion in approximately 10% of cases. EHE cases are typically refractory to therapies, and no anticancer agents are reimbursed for EHE in Australia. METHODS: We report a cohort of nine EHE cases with comprehensive histologic and molecular profiling from the Walter and Eliza Hall Institute of Medical Research Stafford Fox Rare Cancer Program (WEHI-SFRCP) collated via nation-wide referral to the Australian Rare Cancer (ARC) Portal. The diagnoses of EHE were confirmed by histopathological and immunohistochemical (IHC) examination. Molecular profiling was performed using the TruSight Oncology 500 assay, the TruSight RNA fusion panel, whole genome sequencing (WGS), or whole exome sequencing (WES). RESULTS: Molecular analysis of RNA, DNA or both was possible in seven of nine cases. The WWTR1::CAMTA1 fusion was identified in five cases. The YAP1::TFE3 fusion was identified in one case, demonstrating unique morphology compared to cases with the more common WWTR1::CAMTA1 fusion. All tumours expressed typical endothelial markers CD31, ERG, and CD34 and were negative for pan-cytokeratin. Cases with a WWTR1::CAMTA1 fusion displayed high expression of CAMTA1 and the single case with a YAP1::TFE3 fusion displayed high expression of TFE3. Survival was highly variable and unrelated to molecular profile. CONCLUSIONS: This cohort of EHE cases provides molecular and histopathological characterisation and matching clinical information that emphasises the molecular patterns and variable clinical outcomes and adds to our knowledge of this ultra-rare cancer. Such information from multiple studies will advance our understanding, potentially improving treatment options.

Targeting homologous recombination deficiency in uterine leiomyosarcoma.

BACKGROUND: Uterine leiomyosarcoma (uLMS) is a rare and aggressive gynaecological malignancy, with individuals with advanced uLMS having a five-year survival of < 10%. Mutations in the homologous recombination (HR) DNA repair pathway have been observed in ~ 10% of uLMS cases, with reports of some individuals benefiting from poly (ADP-ribose) polymerase (PARP) inhibitor (PARPi) therapy, which targets this DNA repair defect. In this report, we screened individuals with uLMS, accrued nationally, for mutations in the HR repair pathway and explored new approaches to therapeutic targeting. METHODS: A cohort of 58 individuals with uLMS were screened for HR Deficiency (HRD) using whole genome sequencing (WGS), whole exome sequencing (WES) or NGS panel testing. Individuals identified to have HRD uLMS were offered PARPi therapy and clinical outcome details collected. Patient-derived xenografts (PDX) were generated for therapeutic targeting. RESULTS: All 13 uLMS samples analysed by WGS had a dominant COSMIC mutational signature 3; 11 of these had high genome-wide loss of heterozygosity (LOH) (> 0.2) but only two samples had a CHORD score > 50%, one of which had a homozygous pathogenic alteration in an HR gene (deletion in BRCA2). A further three samples harboured homozygous HRD alterations (all deletions in BRCA2), detected by WES or panel sequencing, with 5/58 (9%) individuals having HRD uLMS. All five individuals gained access to PARPi therapy. Two of three individuals with mature clinical follow up achieved a complete response or durable partial response (PR) with the subsequent addition of platinum to PARPi upon minor progression during initial PR on PARPi. Corresponding PDX responses were most rapid, complete and sustained with the PARP1-specific PARPi, AZD5305, compared with either olaparib alone or olaparib plus cisplatin, even in a paired sample of a BRCA2-deleted PDX, derived following PARPi therapy in the patient, which had developed PARPi-resistance mutations in PRKDC, encoding DNA-PKcs. CONCLUSIONS: Our work demonstrates the value of identifying HRD for therapeutic targeting by PARPi and platinum in individuals with the aggressive rare malignancy, uLMS and suggests that individuals with HRD uLMS should be included in trials of PARP1-specific PARPi.

Australian Genomics: Outcomes of a 5-year national program to accelerate the integration of genomics in healthcare.

Australian Genomics is a national collaborative partnership of more than 100 organizations piloting a whole-of-system approach to integrating genomics into healthcare, based on federation principles. In the first five years of operation, Australian Genomics has evaluated the outcomes of genomic testing in more than 5,200 individuals across 19 rare disease and cancer flagship studies. Comprehensive analyses of the health economic, policy, ethical, legal, implementation and workforce implications of incorporating genomics in the Australian context have informed evidence-based change in policy and practice, resulting in national government funding and equity of access for a range of genomic tests. Simultaneously, Australian Genomics has built national skills, infrastructure, policy, and data resources to enable effective data sharing to drive discovery research and support improvements in clinical genomic delivery.

Cancer Predisposition Sequencing Reporter (CPSR): A flexible variant report engine for high-throughput germline screening in cancer.

The value of high-throughput germline genetic testing is increasingly recognized in clinical cancer care. Disease-associated germline variants in cancer patients are important for risk management and surveillance, surgical decisions and can also have major implications for treatment strategies since many are in DNA repair genes. With the increasing availability of high-throughput DNA sequencing in cancer clinics and research, there is thus a need to provide clinically oriented sequencing reports for germline variants and their potential therapeutic relevance on a per-patient basis. To meet this need, we have developed the Cancer Predisposition Sequencing Reporter (CPSR), an open-source computational workflow that generates a structured report of germline variants identified in known cancer predisposition genes, highlighting markers of therapeutic, prognostic and diagnostic relevance. A fully automated variant classification procedure based on more than 30 refined American College of Medical Genetics and Genomics (ACMG) criteria represents an integral part of the workflow. Importantly, the set of cancer predisposition genes profiled in the report can be flexibly chosen from more than 40 virtual gene panels established by scientific experts, enabling customization of the report for different screening purposes and clinical contexts. The report can be configured to also list actionable secondary variant findings, as recommended by ACMG. CPSR demonstrates comparable sensitivity and specificity for the detection of pathogenic variants when compared to other algorithms in the field. Technically, the tool is implemented in Python/R, and is freely available through Docker technology. Source code, documentation, example reports and installation instructions are accessible via the project GitHub page: https://github.com/sigven/cpsr.

Aß-accelerated neurodegeneration caused by Alzheimer's-associated ACE variant R1279Q is rescued by angiotensin system inhibition in mice.

Recent genome-wide association studies identified the angiotensin-converting enzyme gene (ACE) as an Alzheimer's disease (AD) risk locus. However, the pathogenic mechanism by which ACE causes AD is unknown. Using whole-genome sequencing, we identified rare ACE coding variants in AD families and investigated one, ACE1 R1279Q, in knockin (KI) mice. Similar to AD, ACE1 was increased in neurons, but not microglia or astrocytes, of KI brains, which became elevated further with age. Angiotensin II (angII) and angII receptor AT1R signaling were also increased in KI brains. Autosomal dominant neurodegeneration and neuroinflammation occurred with aging in KI hippocampus, which were absent in the cortex and cerebellum. Female KI mice exhibited greater hippocampal electroencephalograph disruption and memory impairment compared to males. ACE variant effects were more pronounced in female KI mice, suggesting a mechanism for higher AD risk in women. Hippocampal neurodegeneration was completely rescued by treatment with brain-penetrant drugs that inhibit ACE1 and AT1R. Although ACE variant-induced neurodegeneration did not depend on ß-amyloid (Aß) pathology, amyloidosis in 5XFAD mice crossed to KI mice accelerated neurodegeneration and neuroinflammation, whereas Aß deposition was unchanged. KI mice had normal blood pressure and cerebrovascular functions. Our findings strongly suggest that increased ACE1/angII signaling causes aging-dependent, Aß-accelerated selective hippocampal neuron vulnerability and female susceptibility, hallmarks of AD that have hitherto been enigmatic. We conclude that repurposed brain-penetrant ACE inhibitors and AT1R blockers may protect against AD.

RAF1 rearrangements are common in pancreatic acinar cell carcinomas.

There is now evidence that gene fusions activating the MAPK pathway are relatively common in pancreatic acinar cell carcinoma with potentially actionable BRAF or RET fusions being found in ~30%. We sought to investigate the incidence of RAF1 fusions in pancreatic malignancies with acinar cell differentiation. FISH testing for RAF1 was undertaken on 30 tumors comprising 25 'pure' acinar cell carcinomas, 2 mixed pancreatic acinar-neuroendocrine carcinomas, 1 mixed acinar cell-low grade neuroendocrine tumor and 2 pancreatoblastomas. RAF1 rearrangements were identified in 5 cases and confirmed by DNA and RNA sequencing to represent oncogenic fusions (GATM-RAF1, GOLGA4-RAF1, PDZRN3-RAF1, HERPUD1-RAF1 and TRIM33-RAF1) and to be mutually exclusive with BRAF and RET fusions, as well as KRAS mutations. Large genome-wide copy number changes were common and included 1q gain and/or 1p loss in all five RAF1 FISH-positive acinar cell carcinomas. RAF1 expression by immunohistochemistry was found in 3 of 5 (60%) of fusion-positive cases and no FISH-negative cases. Phospho-ERK1/2 expression was found in 4 of 5 RAF1-fusion-positive cases. Expression of both RAF1 and phospho-ERK1/2 was heterogeneous and often only detected at the tumor-stroma interface, thus limiting their clinical utility. We conclude that RAF1 gene rearrangements are relatively common in pancreatic acinar cell carcinomas (14.3% to 18.5% of cases) and can be effectively identified by FISH with follow up molecular testing. The combined results of several studies now indicate that BRAF, RET or RAF1 fusions occur in between one third and one-half of these tumors but are extremely rare in other pancreatic malignancies. As these fusions are potentially actionable with currently available therapies, a strong argument can be made to perform FISH or molecular testing on all pancreatic acinar cell carcinomas.

Copy number signatures and mutational processes in ovarian carcinoma.

The genomic complexity of profound copy number aberrations has prevented effective molecular stratification of ovarian cancers. Here, to decode this complexity, we derived copy number signatures from shallow whole-genome sequencing of 117 high-grade serous ovarian cancer (HGSOC) cases, which were validated on 527 independent cases. We show that HGSOC comprises a continuum of genomes shaped by multiple mutational processes that result in known patterns of genomic aberration. Copy number signature exposures at diagnosis predict both overall survival and the probability of platinum-resistant relapse. Measurement of signature exposures provides a rational framework to choose combination treatments that target multiple mutational processes.

Prioritisation of structural variant calls in cancer genomes.

Sensitivity of short read DNA-sequencing for gene fusion detection is improving, but is hampered by the significant amount of noise composed of uninteresting or false positive hits in the data. In this paper we describe a tiered prioritisation approach to extract high impact gene fusion events from existing structural variant calls. Using cell line and patient DNA sequence data we improve the annotation and interpretation of structural variant calls to best highlight likely cancer driving fusions. We also considerably improve on the automated visualisation of the high impact structural variants to highlight the effects of the variants on the resulting transcripts. The resulting framework greatly improves on readily detecting clinically actionable structural variants.

Silencing of the Drosophila ortholog of SOX5 leads to abnormal neuronal development and behavioral impairment.

SOX5 encodes a transcription factor that is expressed in multiple tissues including heart, lung and brain. Mutations in SOX5 have been previously found in patients with amyotrophic lateral sclerosis (ALS) and developmental delay, intellectual disability and dysmorphic features. To characterize the neuronal role of SOX5, we silenced the Drosophila ortholog of SOX5, Sox102F, by RNAi in various neuronal subtypes in Drosophila. Silencing of Sox102F led to misorientated and disorganized michrochaetes, neurons with shorter dendritic arborization (DA) and reduced complexity, diminished larval peristaltic contractions, loss of neuromuscular junction bouton structures, impaired olfactory perception, and severe neurodegeneration in brain. Silencing of SOX5 in human SH-SY5Y neuroblastoma cells resulted in a significant repression of WNT signaling activity and altered expression of WNT-related genes. Genetic association and meta-analyses of the results in several large family-based and case-control late-onset familial Alzheimer's disease (LOAD) samples of SOX5 variants revealed several variants that show significant association with AD disease status. In addition, analysis for rare and highly penetrate functional variants revealed four novel variants/mutations in SOX5, which taken together with functional prediction analysis, suggests a strong role of SOX5 causing AD in the carrier families. Collectively, these findings indicate that SOX5 is a novel candidate gene for LOAD with an important role in neuronal function. The genetic findings warrant further studies to identify and characterize SOX5 variants that confer risk for AD, ALS and intellectual disability.

Deep targeted sequencing of 12 breast cancer susceptibility regions in 4611 women across four different ethnicities.

BACKGROUND: Although genome-wide association studies (GWASs) have identified thousands of disease susceptibility regions, the underlying causal mechanism in these regions is not fully known. It is likely that the GWAS signal originates from one or many as yet unidentified causal variants. METHODS: Using next-generation sequencing, we characterized 12 breast cancer susceptibility regions identified by GWASs in 2288 breast cancer cases and 2323 controls across four populations of African American, European, Japanese, and Hispanic ancestry. RESULTS: After genotype calling and quality control, we identified 137,530 single-nucleotide variants (SNVs); of those, 87.2 % had a minor allele frequency (MAF) <0.005. For SNVs with MAF >0.005, we calculated the smallest number of SNVs needed to obtain a posterior probability set (PPS) such that there is 90 % probability that the causal SNV is included. We found that the PPS for two regions, 2q35 and 11q13, contained less than 5 % of the original SNVs, dramatically decreasing the number of potentially causal SNVs. However, we did not find strong evidence supporting a causal role for any individual SNV. In addition, there were no significant gene-based rare SNV associations after correcting for multiple testing. CONCLUSIONS: This study illustrates some of the challenges faced in fine-mapping studies in the post-GWAS era, most importantly the large sample sizes needed to identify rare-variant associations or to distinguish the effects of strongly correlated common SNVs.

Iterative sorting reveals CD133+ and CD133- melanoma cells as phenotypically distinct populations.

BACKGROUND: The heterogeneity and tumourigenicity of metastatic melanoma is attributed to a cancer stem cell model, with CD133 considered to be a cancer stem cell marker in melanoma as well as other tumours, but its role has remained controversial. METHODS: We iteratively sorted CD133+ and CD133- cells from 3 metastatic melanoma cell lines, and observed tumourigenicity and phenotypic characteristics over 7 generations of serial xeno-transplantation in NOD/SCID mice. RESULTS: We demonstrate that iterative sorting is required to make highly pure populations of CD133+ and CD133- cells from metastatic melanoma, and that these two populations have distinct characteristics not related to the cancer stem cell phenotype. In vitro, gene set enrichment analysis indicated CD133+ cells were related to a proliferative phenotype, whereas CD133- cells were of an invasive phenotype. However, in vivo, serial transplantation of CD133+ and CD133- tumours over 7 generations showed that both populations were equally able to initiate and propagate tumours. Despite this, both populations remained phenotypically distinct, with CD133- cells only able to express CD133 in vivo and not in vitro. Loss of CD133 from the surface of a CD133+ cell was observed in vitro and in vivo, however CD133- cells derived from CD133+ retained the CD133+ phenotype, even in the presence of signals from the tumour microenvironment. CONCLUSION: We show for the first time the necessity of iterative sorting to isolate pure marker-positive and marker-negative populations for comparative studies, and present evidence that despite CD133+ and CD133- cells being equally tumourigenic, they display distinct phenotypic differences, suggesting CD133 may define a distinct lineage in melanoma.

Phosphoproteomic profiling of mouse primary HSPCs reveals new regulators of HSPC mobilization.

Protein phosphorylation is a central mechanism of signal transduction that both positively and negatively regulates protein function. Large-scale studies of the dynamic phosphorylation states of cell signaling systems have been applied extensively in cell lines and whole tissues to reveal critical regulatory networks, and candidate-based evaluations of phosphorylation in rare cell populations have also been informative. However, application of comprehensive profiling technologies to adult stem cell and progenitor populations has been challenging, due in large part to the scarcity of such cells in adult tissues. Here, we combine multicolor flow cytometry with highly efficient 3-dimensional high performance liquid chromatography/mass spectrometry to enable quantitative phosphoproteomic analysis from 200 000 highly purified primary mouse hematopoietic stem and progenitor cells (HSPCs). Using this platform, we identify ARHGAP25 as a novel regulator of HSPC mobilization and demonstrate that ARHGAP25 phosphorylation at serine 363 is an important modulator of its function. Our approach provides a robust platform for large-scale phosphoproteomic analyses performed with limited numbers of rare progenitor cells. Data from our study comprises a new resource for understanding the molecular signaling networks that underlie hematopoietic stem cell mobilization.

VarDict: a novel and versatile variant caller for next-generation sequencing in cancer research.

Accurate variant calling in next generation sequencing (NGS) is critical to understand cancer genomes better. Here we present VarDict, a novel and versatile variant caller for both DNA- and RNA-sequencing data. VarDict simultaneously calls SNV, MNV, InDels, complex and structural variants, expanding the detected genetic driver landscape of tumors. It performs local realignments on the fly for more accurate allele frequency estimation. VarDict performance scales linearly to sequencing depth, enabling ultra-deep sequencing used to explore tumor evolution or detect tumor DNA circulating in blood. In addition, VarDict performs amplicon aware variant calling for polymerase chain reaction (PCR)-based targeted sequencing often used in diagnostic settings, and is able to detect PCR artifacts. Finally, VarDict also detects differences in somatic and loss of heterozygosity variants between paired samples. VarDict reprocessing of The Cancer Genome Atlas (TCGA) Lung Adenocarcinoma dataset called known driver mutations in KRAS, EGFR, BRAF, PIK3CA and MET in 16% more patients than previously published variant calls. We believe VarDict will greatly facilitate application of NGS in clinical cancer research.

HIV Integration Site Analysis of Cellular Models of HIV Latency with a Probe-Enriched Next-Generation Sequencing Assay.

UNLABELLED: Antiretroviral therapy (ART) is successful in the suppression of HIV but cannot target and eradicate the latent proviral reservoir. The location of retroviral integration into the human genome is thought to play a role in the clonal expansion of infected cells and HIV persistence. We developed a high-throughput targeted sequence capture assay that uses a pool of HIV-specific probes to enrich Illumina libraries prior to deep sequencing. Using an expanded clonal population of ACH-2 cells, we demonstrate that this sequence capture assay has an extremely low false-positive rate. This assay assessed four cellular models commonly used to study HIV latency and latency-reversing agents: ACH-2 cells, J-Lat cells, the Bcl-2-transduced primary CD4(+)model, and the cultured TCM(central memory) CD4(+)model. HIV integration site characteristics and genes were compared between these cellular models and to previously reported patient data sets. Across these cellular models, there were significant differences in integration site characteristics, including orientation relative to that of the host gene, the proportion of clonally expanded sites, and the proportion located within genic regions and exons. Despite a greater diversity of minority integration sites than expected in ACH-2 cells, their integration site characteristics consistently differed from those of the other models and from the patient samples. Gene ontology analysis of highly represented genes from the patient samples found little overlap with HIV-containing genes from the cell lines. These findings show that integration site differences exist among the commonly used cellular models of HIV latency and in comparison to integration sites found in patient samples. IMPORTANCE: Despite the success of ART, currently there is no successful therapy to eradicate integrated proviruses. Cellular models of HIV latency are used to test the efficacy of latency-reversing agents, but it is unclear how well these models reflect HIV integration into the human genome in vivo We have developed a novel probe-based sequence enrichment assay to sequence and analyze integrated HIV. We compared HIV integration site characteristics between four cellular models and to previously described patient data sets. Significant differences were detected in the distribution of HIV integration sites between cellular models of HIV latency and compared to data sets from patient samples. The results from this study have implications for how well these cellular models of HIV infection truly reflect HIV integration in vivo and their applicability in drug discovery for novel latency-reversing agents.

RNA Sequencing Identifies Novel Translational Biomarkers of Kidney Fibrosis.

CKD is the gradual, asymptomatic loss of kidney function, but current tests only identify CKD when significant loss has already happened. Several potential biomarkers of CKD have been reported, but none have been approved for preclinical or clinical use. Using RNA sequencing in a mouse model of folic acid-induced nephropathy, we identified ten genes that track kidney fibrosis development, the common pathologic finding in patients with CKD. The gene expression of all ten candidates was confirmed to be significantly higher (approximately ten- to 150-fold) in three well established, mechanistically distinct mouse models of kidney fibrosis than in models of nonfibrotic AKI. Protein expression of these genes was also high in the folic acid model and in patients with biopsy-proven kidney fibrosis. mRNA expression of the ten genes increased with increasing severity of kidney fibrosis, decreased in response to therapeutic intervention, and increased only modestly (approximately two- to five-fold) with liver fibrosis in mice and humans, demonstrating specificity for kidney fibrosis. Using targeted selected reaction monitoring mass spectrometry, we detected three of the ten candidates in human urine: cadherin 11 (CDH11), macrophage mannose receptor C1 (MRC1), and phospholipid transfer protein (PLTP). Furthermore, urinary levels of each of these three proteins distinguished patients with CKD (n=53) from healthy individuals (n=53; P<0.05). In summary, we report the identification of urinary CDH11, MRC1, and PLTP as novel noninvasive biomarkers of CKD.

Mycoplasma Infection Alters Cancer Stem Cell Properties in Vitro.

Cancer cell lines can be useful to model cancer stem cells. Infection with Mycoplasma species is an insidious problem in mammalian cell culture. While investigating stem-like properties in early passage melanoma cell lines, we noted poorly reproducible results from an aliquot of a cell line that was later found to be infected with Mycoplasma hyorhinis. Deliberate infection of other early passage melanoma cell lines aliquots induced variable and unpredictable effects on expression of putative cancer stem cell markers, clonogenicity, proliferation and global gene expression. Cell lines established in stem cell media (SCM) were equally susceptible. Mycoplasma status is rarely reported in publications using cultured cells to study the cancer stem cell hypothesis. Our work highlights the importance of surveillance for Mycoplasma infection while using any cultured cells to interrogate tumor heterogeneity.

An RNA-binding Protein, Lin28, Recognizes and Remodels G-quartets in the MicroRNAs (miRNAs) and mRNAs It Regulates.

Lin28 is an evolutionarily conserved RNA-binding protein that inhibits processing of pre-let-7 microRNAs (miRNAs) and regulates translation of mRNAs that control developmental timing, pluripotency, metabolism, and tumorigenesis. The RNA features that mediate Lin28 binding to the terminal loops of let-7 pre-miRNAs and to Lin28-responsive elements (LREs) in mRNAs are not well defined. Here we show that Lin28 target datasets are enriched for RNA sequences predicted to contain stable planar structures of 4 guanines known as G-quartets (G4s). The imino NMR spectra of pre-let-7 loops and LREs contain resonances characteristic of G4 hydrogen bonds. These sequences bind to a G4-binding fluorescent dye, N-methyl-mesoporphyrin IX (NMM). Mutations and truncations in the RNA sequence that prevent G4 formation also prevent Lin28 binding. The addition of Lin28 to a pre-let-7 loop or an LRE reduces G4 resonance intensity and NMM binding, suggesting that Lin28 may function to remodel G4s. Further, we show that NMM inhibits Lin28 binding. Incubation of a human embryonal carcinoma cell line with NMM reduces its stem cell traits. In particular it increases mature let-7 levels, decreases OCT4, HMGA1, CCNB1, CDK4, and Lin28A protein, decreases sphere formation, and inhibits colony formation. Our results suggest a previously unknown structural feature of Lin28 targets and a new strategy for manipulating Lin28 function.

Clonal dynamics of native haematopoiesis.

It is currently thought that life-long blood cell production is driven by the action of a small number of multipotent haematopoietic stem cells. Evidence supporting this view has been largely acquired through the use of functional assays involving transplantation. However, whether these mechanisms also govern native non-transplant haematopoiesis is entirely unclear. Here we have established a novel experimental model in mice where cells can be uniquely and genetically labelled in situ to address this question. Using this approach, we have performed longitudinal analyses of clonal dynamics in adult mice that reveal unprecedented features of native haematopoiesis. In contrast to what occurs following transplantation, steady-state blood production is maintained by the successive recruitment of thousands of clones, each with a minimal contribution to mature progeny. Our results demonstrate that a large number of long-lived progenitors, rather than classically defined haematopoietic stem cells, are the main drivers of steady-state haematopoiesis during most of adulthood. Our results also have implications for understanding the cellular origin of haematopoietic disease.

Sequencing of captive target transcripts identifies the network of regulated genes and functions of primate-specific miR-522.

Identifying microRNA (miRNA)-regulated genes is key to understanding miRNA function. However, many miRNA recognition elements (MREs) do not follow canonical "seed" base-pairing rules, making identification of bona fide targets challenging. Here, we apply an unbiased sequencing-based systems approach to characterize miR-522, a member of the oncogenic primate-specific chromosome 19 miRNA cluster, highly expressed in poorly differentiated cancers. To identify miRNA targets, we sequenced full-length transcripts captured by a biotinylated miRNA mimic. Within these targets, mostly noncanonical MREs were identified by sequencing RNase-resistant fragments. miR-522 overexpression reduced mRNA, protein levels, and luciferase activity of >70% of a random list of candidate target genes and MREs. Bioinformatic analysis suggested that miR-522 regulates cell proliferation, detachment, migration, and epithelial-mesenchymal transition. miR-522 induces G1 cell-cycle arrest and causes cells to detach without anoikis, become invasive, and express mesenchymal genes. Thus, our method provides a simple but effective technique for identifying miRNA-regulated genes and biological function.