De Novo Assembly of the Genome of the Sea Urchin Paracentrotus lividus (Lamarck 1816).

The Mediterranean purple sea urchin Paracentrotus lividus (Lamarck 1816) is a remarkable model system for molecular, evolutionary and cell biology studies, particularly in the field of developmental biology. We sequenced the genome, performed a de novo assembly, and analysed the assembly content. The genome of P. lividus was sequenced using Illumina NextSeq 500 System (Illumina) in a 2 × 150 paired-end format. More than 30,000 open reading frames (ORFs), (more than 8000 are unique), were identified and analysed to provide molecular tools accessible for the scientific community. In particular, several genes involved in complex innate immune responses, oxidative metabolism, signal transduction, and kinome, as well as genes regulating the membrane receptors, were identified in the P. lividus genome. In this way, the employment of the Mediterranean sea urchin for investigations and comparative analyses was empowered, leading to the explanation of cis-regulatory networks and their evolution in a key developmental model occupying an important evolutionary position with respect to vertebrates and humans.

A Next-Generation Sequencing Study in a Cohort of Sicilian Patients with Parkinson's Disease.

Parkinson's disease (PD) is a multisystem and multifactorial disorder and, therefore, the application of modern genetic techniques may assist in unraveling its complex pathophysiology. We conducted a clinical-demographic evaluation of 126 patients with PD, all of whom were Caucasian and of Sicilian ancestry. DNA was extracted from the peripheral blood for each patient, followed by sequencing using a Next-Generation Sequencing system. This system was based on a custom gene panel comprising 162 genes. The sample underwent further filtering, taking into account the allele frequencies of genetic variants, their presence in the Human Gene Mutation Database, and their association in the literature with PD or other movement/neurodegenerative disorders. The largest number of variants was identified in the leucine-rich repeat kinase 2 (LRRK2) gene. However, variants in other genes, such as acid beta-glucosidase (GBA), DNA polymerase gamma catalytic subunit (POLG), and parkin RBR E3 ubiquitin protein ligase (PRKN), were also discovered. Interestingly, some of these variants had not been previously associated with PD. Enhancing our understanding of the genetic basis of PD and identifying new variants possibly linked to the disease will contribute to improved diagnostic accuracy, therapeutic developments, and prognostic insights for affected individuals.

NGS study in a sicilian case series with a genetic diagnosis for Gerstmann-Sträussler-Scheinker syndrome (PRNP, p.P102L).

BACKGROUND: Gerstmann Sträussler Scheinker (GSS) is an inherited, invariably fatal prion disease. Like other human prion diseases, GSS is caused by missense mutations in the prion protein (PrP) gene (PRNP), and by the formation and overtime accumulation of the misfolded, pathogenic scrapie PrP (PrPSc). The first mutation identified in the PRNP gene, and the one blamed as the main cause of the disease, is c.C305T:p.P102L. METHODS AND RESULTS: The Sanger sequencing method was performed on the PRNP gene for the detection of c.C305T:p.P102L mutations in a cohort of 10 subjects; moreover, a study was carried out, using Next Generation Sequencing (NGS), by sequencing a group of genes related to amyotrophic lateral sclerosis (ALS), Alzheimer's disease (AD), movement disorders and dementia which show a phenotypic profile similar to that of GSS. The results obtained from the study using NGS indicate the potential role of other genetic variants which could contribute to the various GSS phenotypes. CONCLUSIONS: In conclusion, we highlight the large clinical variability in subjects presenting with GSS and p.P102L, as well as the hypothesis that the mutation in PrP codon 102 alone is not sufficient to trigger the cardinal clinical signs of the disease; furthermore, we do not exclude the possibility that further genetic variants play a decisive role in the aspects of the various phenotypes with which GSS manifests itself.

Role and Dysregulation of miRNA in Patients with Parkinson's Disease.

Parkinson's disease (PD) is a neurodegenerative synucleinopathy that has a not yet fully understood molecular pathomechanism behind it. The role of risk genes regulated by small non-coding RNAs, or microRNAs (miRNAs), has also been highlighted in PD, where they may influence disease progression and comorbidities. In this case-control study, we analyzed miRNAs on peripheral blood mononuclear cells by means of RNA-seq in 30 participants, with the aim of identifying miRNAs differentially expressed in PD compared to age-matched healthy controls. Additionally, we investigated the pathways influenced by differentially expressed miRNAs and assessed whether a specific pathway could potentially be associated with PD susceptibility (enrichment analyses performed using the Ingenuity Pathway Analysis tools). Overall, considering that the upregulation of miRNAs might be related with the downregulation of their messenger RNA targets, and vice versa, we found several putative targets of dysregulated miRNAs (i.e., upregulated: hsa-miR-1275, hsa-miR-23a-5p, hsa-miR-432-5p, hsa-miR-4433b-3p, and hsa-miR-4443; downregulated: hsa-miR-142-5p, hsa-miR-143-3p, hsa-miR-374a-3p, hsa-miR-542-3p, and hsa-miR-99a-5p). An inverse connection between cancer and neurodegeneration, called "inverse comorbidity", has also been noted, showing that some genes or miRNAs may be expressed oppositely in neurodegenerative disorders and in some cancers. Therefore, it may be reasonable to consider these miRNAs as potential diagnostic markers and outcome measures.

CCR3 gene overexpression in patients with Down syndrome.

Chromosome 21 trisomy or Down syndrome (DS) is the most common genetic cause of intellectual disability (ID). DS is also associated with hypotonia, muscle weakness, autoimmune diseases, and congenital heart disease. C-C chemokine receptor type 3 (CCR3) plays a role in inflammatory, autoimmune, and neuronal migration mechanisms. The present study aimed to evaluate the expression of the CCR3 gene by NGS and qRT-PCR in patients with DS and normal controls (NC). The CCR3 gene was over-expressed in DS patients compared to NC. These data suggest that an over-expression of the CCR3 gene is associated with the phenotype of patients with DS.

GPR56 gene down-regulation in patients with Klinefelter Syndrome: a candidate for infertility?

BACKGROUND: The etiology of azoospermia in patients with Klinefelter Syndrome (KS) is still unknown. The protein codified by the G protein-couple receptor 56 (GPR56) belongs to the adhesion family of G protein-coupled receptors (GPRs). Its mutations are involved in the pathogenesis of intellectual disability and, according to animal studies, infertility. As the expression of GPR56 in patients with KS has not been investigated so far, this study was undertaken with the purpose of evaluating its expression in peripheral blood mononuclear cells (PBMCs) of patients with KS and normal controls. METHODS: This age-matched case-control study was performed in 10 patients with KS and 10 controls. Patients and controls underwent to blood sampling for next-generation sequencing (NGS) analysis, and differentially expressed mRNAs were identified using DESeq2 v.1.12. QRT-PCR was used to validate the results obtained by NGS analysis. TaqMan Gene Expression Assay primers were used to carry out the real-time PCR analysis for GPR56. RESULTS: GPR56 was down-regulated by -2081-fold (q-value <0.05) in PBMCs of patients with KS compared to controls. NGS data were confirmed by QRT-PCR. CONCLUSIONS: The possible contribution of the GPR56 gene down-regulation in the pathogenesis of spermatogenic failure in patients with KS is worthy to be further explored.

Functional evidence of mTORß splice variant involvement in the pathogenesis of congenital heart defects.

mTOR dysregulation has been described in pathological conditions, such as cardiovascular and overgrowth disorders. Here we report on the first case of a patient with a complex congenital heart disease and an interstitial duplication in the short arm of chromosome 1, encompassing part of the mTOR gene. Our results suggest that an intragenic mTOR microduplication might play a role in the pathogenesis of non-syndromic congenital heart defects (CHDs) due to an upregulation of mTOR/Rictor and consequently an increased phosphorylation of PI3K/AKT and MEK/ERK signaling pathways in patient-derived amniocytes. This is the first report which shows a causative role of intragenic mTOR microduplication in the etiology of an isolated complex CHD.

SOX13 gene downregulation in peripheral blood mononuclear cells of patients with Klinefelter syndrome.

Klinefelter syndrome (KS) is the most common sex chromosome disorder in men. It is characterized by germ cell loss and other variable clinical features, including autoimmunity. The sex-determining region of Y (SRY)-box 13 (Sox13) gene is expressed in mouse spermatogonia. In addition, it has been identified as islet cell autoantigen 12 (ICA12), which is involved in the pathogenesis of autoimmune diseases, including type 1 diabetes mellitus (DM) and primary biliary cirrhosis. Sox13 expression has never been investigated in patients with KS. In this age-matched, case-control study performed on ten patients with KS and ten controls, we found that SOX13 is significantly downregulated in peripheral blood mononuclear cells of patients with KS compared to controls. This finding might be consistent with the germ cell loss typical of patients with KS. However, the role of Sox13 in the pathogenesis of germ cell loss and humoral autoimmunity in patients with KS deserves to be further explored.

Long non-coding RNA GAS5 expression in patients with Down syndrome.

Trisomy 21, also known as Down Syndrome (DS), is the most common chromosome abnormality and causes intellectual disability. Long non-coding RNA (lncRNA) growth arrest-specific 5 (GAS5), whose differential expression has recently been reported in patients with Klinefelter syndrome, has been addressed to play a role in the development of inflammatory and autoimmune diseases, vascular endothelial cells apoptosis and atherosclerosis, all being common features in patients with DS. Therefore, the aim of this study was to assess the lncRNA GAS5 expression profile in DS patients and in controls. lncRNA GAS5 levels were evaluated by qRT-PCR assay in 23 patients with DS and 23 age-matched controls. A significant lncRNA GAS5 down-regulation was observed in patients with DS by RT-PCR analysis, The RNA sequencing experiments confirmed the qRT-PCR data. LncRNA GAS5 down-expression may play a role in the development of some typical features of the patients with DS and, particularly, in inflammatory and autoimmune diseases.

KCTD15 is overexpressed in human childhood B-cell acute lymphoid leukemia.

Leukemic cells originate from the malignant transformation of undifferentiated myeloid/lymphoid hematopoietic progenitors normally residing in bone marrow. As the precise molecular mechanisms underlying this heterogeneous disease are yet to be disclosed, the identification and the validation of novel actors in leukemia is of extreme importance. Here, we show that KCTD15, a member of the emerging class of KCTD ((K)potassium Channel Tetramerization Domain containing) proteins, is strongly upregulated in patients affected by B-cell type acute lymphoblastic leukemia (B-ALL) and in continuous cell lines (RS4;11, REH, TOM-1, SEM) derived from this form of childhood leukemia. Interestingly, KCTD15 downregulation induces apoptosis and cell death suggesting that it has a role in cellular homeostasis and proliferation. In addition, stimulation of normal lymphocytes with the pokeweed mitogen leads to increased KCTD15 levels in a fashion comparable to those observed in proliferating leukemic cells. In this way, the role of KCTD15 is likely not confined to the B-ALL pathological state and extends to activation and proliferation of normal lymphocytes. Collectively, data here presented indicate that KCTD15 is an important and hitherto unidentified player in childhood lymphoid leukemia, and its study could open a new scenario for the identification of altered and still unknown molecular pathways in leukemia.

Identification of long non­coding RNA expression patterns useful for molecular­based classification of type I endometrial cancers.

Endometrial cancer is the most frequently diagnosed gynecologic malignant disease. Although several genetic alterations have been associated with the increased risk of endometrial cancer, to date, the diagnosis and prognosis still rely on morphological features of the tumor, such as histological type, grading and invasiveness. As molecular­based classification is desirable for optimal treatment and prognosis of these cancers, we explored the potential of lncRNAs as molecular biomarkers. To this end, we first identified by RNA sequencing (RNA­Seq) a set of lncRNAs differentially expressed in cancer vs. normal endometrial tissues, a result confirmed also by analysis of normal and cancerous endometrium RNA­Seq data from TCGA (The Cancer Genome Atlas). A significant association of a subset of these differentially expressed lncRNAs with tumor grade was then determined in 405 TCGA endometrial cancer profiles. Integrating endometrial cancer­specific expression profiles of long and small non­coding RNAs, a functional association network was then identified. These results describe for the first time a functional ῾coreá¾½ network, comprising small and long RNAs, whose deregulation is associated with endometrial neoplastic transformation, representing a set of cancer biomarkers that can be monitored and targeted for diagnosis, follow­up and therapy of these tumors.

Identification of a novel truncating mutation in PALB2 gene by a multigene sequencing panel for mutational screening of breast cancer risk-associated and related genes.

BACKGROUND: Breast cancer (BC) is the most common neoplasm in women, with 5%-10% patients showing a familial predisposition, where germline mutations in BRCA1/BRCA2 genes are found in -20% of cases. Next-generation sequencing (NGS) is among the best available options for genetic screening, providing several benefits that include enhanced sensitivity and unbiased mutation detection. PALB2 (partner and localizer of BRCA2) is a cancer predisposing gene recently described that encodes a protein partner of BRCA2 involved in DNA double-strand break repair and cell cycle control. The DNA damage response represents a key cellular event, targeted by innovative anticancer therapies, including those based on poly (ADP-ribose) polymerase (PARP) inhibitors targeting PARP1 and PARP2 enzymes, activated by DNA damage and involved in single-strand break and base excision repair. METHODS: Genomic DNA was isolated from 34 patient samples and four BC cell lines, as controls, and 27 breast cancer predisposing genes belonging to the BRCA1/BRCA2 and PARP pathways were sequenced by NGS. RESULTS: The panel described here allowed identification of several sequence variations in most investigated genes, among which we found a novel truncating mutation in PALB2. CONCLUSIONS: The NGS-based strategy designed here for molecular analysis of a customized panel of BC predisposing and related genes was found to perform effectively, providing a comprehensive exploration of all genomic sequences of the investigated genes. It is thus useful for BC molecular diagnosis, in particular for familiar cases where alterations in routinely investigated genes, such as BRCAs, result to be absent.

The nuclear receptor ERß engages AGO2 in regulation of gene transcription, RNA splicing and RISC loading.

BACKGROUND: The RNA-binding protein Argonaute 2 (AGO2) is a key effector of RNA-silencing pathways It exerts a pivotal role in microRNA maturation and activity and can modulate chromatin remodeling, transcriptional gene regulation and RNA splicing. Estrogen receptor beta (ERß) is endowed with oncosuppressive activities, antagonizing hormone-induced carcinogenesis and inhibiting growth and oncogenic functions in luminal-like breast cancers (BCs), where its expression correlates with a better prognosis of the disease. RESULTS: Applying interaction proteomics coupled to mass spectrometry to characterize nuclear factors cooperating with ERß in gene regulation, we identify AGO2 as a novel partner of ERß in human BC cells. ERß-AGO2 association was confirmed in vitro and in vivo in both the nucleus and cytoplasm and is shown to be RNA-mediated. ChIP-Seq demonstrates AGO2 association with a large number of ERß binding sites, and total and nascent RNA-Seq in ERß + vs ERß - cells, and before and after AGO2 knock-down in ERß + cells, reveals a widespread involvement of this factor in ERß-mediated regulation of gene transcription rate and RNA splicing. Moreover, isolation and sequencing by RIP-Seq of ERß-associated long and small RNAs in the cytoplasm suggests involvement of the nuclear receptor in RISC loading, indicating that it may also be able to directly control mRNA translation efficiency and stability. CONCLUSIONS: These results demonstrate that AGO2 can act as a pleiotropic functional partner of ERß, indicating that both factors are endowed with multiple roles in the control of key cellular functions.

Allelic Complexity in Long QT Syndrome: A Family-Case Study.

Congenital long QT syndrome (LQTS) is associated with high genetic and allelic heterogeneity. In some cases, more than one genetic variant is identified in the same (compound heterozygosity) or different (digenic heterozygosity) genes, and subjects with multiple pathogenic mutations may have a more severe disease. Standard-of-care clinical genetic testing for this and other arrhythmia susceptibility syndromes improves the identification of complex genotypes. Therefore, it is important to distinguish between pathogenic mutations and benign rare variants. We identified four genetic variants (KCNQ1-p.R583H, KCNH2-p.C108Y, KCNH2-p.K897T, and KCNE1-p.G38S) in an LQTS family. On the basis of in silico analysis, clinical data from our family, and the evidence from previous studies, we analyzed two mutated channels, KCNQ1-p.R583H and KCNH2-p.C108Y, using the whole-cell patch clamp technique. We found that KCNQ1-p.R583H was not associated with a severe functional impairment, whereas KCNH2-p.C108Y, a novel variant, encoded a non-functional channel that exerts dominant-negative effects on the wild-type. Notably, the common variants KCNH2-p.K897T and KCNE1-p.G38S were previously reported to produce more severe phenotypes when combined with disease-causing alleles. Our results indicate that the novel KCNH2-C108Y variant can be a pathogenic LQTS mutation, whereas KCNQ1-p.R583H, KCNH2-p.K897T, and KCNE1-p.G38S could be LQTS modifiers.

Specific patterns of PIWI-interacting small noncoding RNA expression in dysplastic liver nodules and hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) is the result of a stepwise process, often beginning with development within a cirrhotic liver of premalignant lesions, morphologically characterized by low- (LGDN) and high-grade (HGDN) dysplastic nodules. PIWI-interacting RNAs (piRNAs) are small noncoding RNAs (sncRNAs), 23-35 nucleotide-long, exerting epigenetic and post-transcriptional regulation of gene expression. Recently the PIWI-piRNA pathway, best characterized in germline cells, has been identified also in somatic tissues, including stem and cancer cells, where it influences key cellular processes.Small RNA sequencing was applied to search for liver piRNAs and to profile their expression patterns in cirrhotic nodules (CNs), LGDN, HGDN, early HCC and progressed HCC (pHCC), analyzing 55 samples (14 CN, 9 LGDN, 6 HGDN, 6 eHCC and 20 pHCC) from 17 patients, aiming at identifying possible relationships between these sncRNAs and liver carcinogenesis. We identified a 125 piRNA expression signature that characterize HCC from matched CNs, correlating also to microvascular invasion in HCC. Functional analysis of the predicted RNA targets of deregulated piRNAs indicates that these can target key signaling pathways involved in hepatocarcinogenesis and HCC progression, thereby affecting their activity. Interestingly, 24 piRNAs showed specific expression patterns in dysplastic nodules, respect to cirrhotic liver and/or pHCC.The results demonstrate that the PIWI-piRNA pathway is active in human liver, where it represents a new player in the molecular events that characterize hepatocarcinogenesis, from early stages to pHCC. Furthermore, they suggest that piRNAs might be new disease biomarkers, useful for differential diagnosis of dysplastic and neoplastic liver lesions.

Leptin as a mediator of tumor-stromal interactions promotes breast cancer stem cell activity.

Breast cancer stem cells (BCSCs) play crucial roles in tumor initiation, metastasis and therapeutic resistance. A strict dependency between BCSCs and stromal cell components of tumor microenvironment exists. Thus, novel therapeutic strategies aimed to target the crosstalk between activated microenvironment and BCSCs have the potential to improve clinical outcome. Here, we investigated how leptin, as a mediator of tumor-stromal interactions, may affect BCSC activity using patient-derived samples (n = 16) and breast cancer cell lines, and determined the potential benefit of targeting leptin signaling in these model systems. Conditioned media (CM) from cancer-associated fibroblasts and breast adipocytes significantly increased mammosphere formation in breast cancer cells and depletion of leptin from CM completely abrogated this effect. Mammosphere cultures exhibited increased leptin receptor (OBR) expression and leptin exposure enhanced mammosphere formation. Microarray analyses revealed a similar expression profile of genes involved in stem cell biology among mammospheres treated with CM and leptin. Interestingly, leptin increased mammosphere formation in metastatic breast cancers and expression of OBR as well as HSP90, a target of leptin signaling, were directly correlated with mammosphere formation in metastatic samples (r = 0.68/p = 0.05; r = 0.71/p = 0.036, respectively). Kaplan-Meier survival curves indicated that OBR and HSP90 expression were associated with reduced overall survival in breast cancer patients (HR = 1.9/p = 0.022; HR = 2.2/p = 0.00017, respectively). Furthermore, blocking leptin signaling by using a full leptin receptor antagonist significantly reduced mammosphere formation in breast cancer cell lines and patient-derived samples. Our results suggest that leptin/leptin receptor signaling may represent a potential therapeutic target that can block the stromal-tumor interactions driving BCSC-mediated disease progression.

Global gene expression profile of normal and regenerating liver in young and old mice.

The ability of the liver to regenerate and adjust its size after two/third partial hepatectomy (PH) is impaired in old rodents and humans. Here, we investigated by microarray analysis the expression pattern of hepatic genes in young and old untreated mice and the differences in gene expression profile following PH. Of the 10,237 messenger RNAs that had detectable expression, only 108 displayed a greater than 2-fold modification in gene expression levels between the two groups. These genes were involved in inflammatory and immune response, xenobiotics, and lipid and glucose metabolism. To identify the genes responsible for the different regenerative response, 10-week and 18-month-old mice subjected to PH were sacrificed at different time intervals after surgery. The results showed that 2463 transcripts had significantly different expression post PH between the two groups. However, in spite of impaired liver regeneration in old mice, cell cycle genes were similarly modified in both groups, the only exception being cyclin D1 gene which was up-regulated soon after PH in young mice, but mostly down-regulated in aged animals. Surprisingly, while in young hepatectomized mice, Yap messenger RNA (mRNA) expression was not significantly enhanced and protein expression essentially reflected the progression into cell cycle, its mRNA and protein levels were robustly increased in the liver of aged animals. Furthermore, a significant change of the age-related expression of the size regulator Yes-associated protein (YAP) was observed. Unexpectedly, while in young hepatectomized mice, Yap mRNA expression was not significantly enhanced and protein expression essentially reflected the progression into cell cycle, its mRNA and protein levels were robustly increased in the liver of aged animals. Moreover, when PH was performed on mitogen-induced enlarged livers, the earlier restoration of the original liver mass compared to animals subjected to PH only led to YAP down-regulation concomitantly with cyclin D1 up-regulation. Our data suggest that YAP activation is a size-dependent homeostatic mechanism that does not necessarily reflect cell cycle progression.

Small non-coding RNA deregulation in endometrial carcinogenesis.

Small non-coding RNAs (sncRNAs) represent a heterogeneous group of <200nt-long transcripts comprising microRNAs, PIWI-interacting RNAs (piRNAs) and small-nucleolar-RNAs (snoRNAs) involved in physiological and pathological processes such as carcinogenesis and tumor progression. Aberrant sncRNA expression in cancer has been associated with specific clinical phenotypes, grading, staging, metastases development and resistance to therapy.Aim of the present work is to study the role of sncRNAs in endometrial carcinogenesis. Changes in sncRNA expression were identified by high-throughput genomic analysis of paired normal, hyperplastic and cancerous endometrial tissues obtained by endometrial biopsies (n = 10). Using smallRNA sequencing and microarrays we identified significant differences in sncRNA expression pattern between normal, hyperplastic and neoplastic endometrium. This led to the definition of a sncRNA signature (129 microRNAs, 2 of which not previously described, 10 piRNAs and 3 snoRNAs) of neoplastic transformation. Functional bioinformatics analysis identified as downstream targets multiple signaling pathways potentially involved in the hyperplastic and neoplastic tissue responses, including Wnt/ß-catenin, and ERK/MAPK and TGF-ß-Signaling.Considering the regulatory role of sncRNAs, this newly identified sncRNA signature is likely to reflect the events leading to endometrial cancer, which can be exploited to dissect the carcinogenic process including novel biomarkers for early and non-invasive diagnosis of these tumors.

Identification of cytoplasmic proteins interacting with unliganded estrogen receptor α and ß in human breast cancer cells.

Estrogen receptor subtypes (ERα and ERß) are transcription factors sharing a similar structure but exerting opposite roles in breast cancer cells. Besides the well-characterized genomic actions of nuclear ERs upon ligand binding, specific actions of ligand-free ERs in the cytoplasm also affect cellular functions. The identification of cytoplasmic interaction partners of unliganded ERα and ERß may help characterize the molecular basis of the extra-nuclear mechanism of action of these receptors, revealing novel mechanisms to explain their role in breast cancer response or resistance to endocrine therapy. To this aim, cytoplasmic extracts from human breast cancer MCF-7 cells stably expressing tandem affinity purification-tagged ERα and ERß and maintained in estrogen-free medium were subject to affinity-purification and MS analysis, leading to the identification of 84 and 142 proteins associated with unliganded ERα and ERß, respectively. Functional analyses of ER subtype-specific interactomes revealed significant differences in the molecular pathways targeted by each receptor in the cytoplasm. This work, reporting the first identification of the unliganded ERα and ERß cytoplasmic interactomes in breast cancer cells, provides novel experimental evidence on the nongenomic effects of ERs in the absence of hormonal stimulus. All MS data have been deposited in the ProteomeXchange with identifier PXD001202 (http://proteomecentral.proteomexchange.org/dataset/PXD001202).

RNA sequencing identifies specific PIWI-interacting small non-coding RNA expression patterns in breast cancer.

PIWI-interacting small non-coding RNAs (piRNAs) are genetic and epigenetic regulatory factors in germline cells, where they maintain genome stability, are involved in RNA silencing and regulate gene expression. We found that the piRNA biogenesis and effector pathway are present in human breast cancer (BC) cells and, analyzing smallRNA-Seq data generated from BC cell lines and tumor biopsies, we identified >100 BC piRNAs, including some very abundant and/or differentially expressed in mammary epithelial compared to BC cells, where this was influenced by estrogen or estrogen receptor ß, and in cancer respect to normal breast tissues. A search for mRNAs targeted by the BC piRNome revealed that eight piRNAs showing a specific expression pattern in breast tumors target key cancer cell pathways. Evidence of an active piRNA pathway in BC suggests that these small non-coding RNAs do exert transcriptional and post-transcriptional gene regulatory actions also in cancer cells.