Dysregulation of the PRUNE2/PCA3 genetic axis in human prostate cancer: from experimental discovery to validation in two independent patient cohorts.

Background: We have previously shown that the long non-coding (lnc)RNA prostate cancer associated 3 (PCA3; formerly prostate cancer antigen 3) functions as a trans-dominant negative oncogene by targeting the previously unrecognized prostate cancer suppressor gene PRUNE2 (a homolog of the Drosophila prune gene), thereby forming a functional unit within a unique allelic locus in human cells. Here, we investigated the PCA3/PRUNE2 regulatory axis from early (tumorigenic) to late (biochemical recurrence) genetic events during human prostate cancer progression. Methods: The reciprocal PCA3 and PRUNE2 gene expression relationship in paired prostate cancer and adjacent normal prostate was analyzed in two independent retrospective cohorts of clinically annotated cases post-radical prostatectomy: a single-institutional discovery cohort (n=107) and a multi-institutional validation cohort (n=497). We compared the tumor gene expression of PCA3 and PRUNE2 to their corresponding expression in the normal prostate. We also serially examined clinical/pathological variables including time to disease recurrence. Results: We consistently observed increased expression of PCA3 and decreased expression of PRUNE2 in prostate cancer compared with the adjacent normal prostate across all tumor grades and stages. However, there was no association between the relative gene expression levels of PCA3 or PRUNE2 and time to disease recurrence, independent of tumor grades and stages. Conclusions: We concluded that upregulation of the lncRNA PCA3 and targeted downregulation of the protein-coding PRUNE2 gene in prostate cancer could be early (rather than late) molecular events in the progression of human prostate tumorigenesis but are not associated with biochemical recurrence. Further studies of PCA3/PRUNE2 dysregulation are warranted. Funding: We received support from the Human Tissue Repository and Tissue Analysis Shared Resource from the Department of Pathology of the University of New Mexico School of Medicine and a pilot award from the University of New Mexico Comprehensive Cancer Center. RP and WA were supported by awards from the Levy-Longenbaugh Donor-Advised Fund and the Prostate Cancer Foundation. EDN reports research fellowship support from the Brazilian National Council for Scientific and Technological Development (CNPq), Brazil, and the Associação Beneficente Alzira Denise Hertzog Silva (ABADHS), Brazil. This work has been funded in part by the NCI Cancer Center Support Grants (CCSG; P30) to the University of New Mexico Comprehensive Cancer Center (CA118100) and the Rutgers Cancer Institute of New Jersey (CA072720).

Exome and Tissue-Associated Microbiota as Predictive Markers of Response to Neoadjuvant Treatment in Locally Advanced Rectal Cancer.

The clinical and pathological responses to multimodal neoadjuvant therapy in locally advanced rectal cancers (LARCs) remain unpredictable, and robust biomarkers are still lacking. Recent studies have shown that tumors present somatic molecular alterations related to better treatment response, and it is also clear that tumor-associated bacteria are modulators of chemotherapy and immunotherapy efficacy, therefore having implications for long-term survivorship and a good potential as the biomarkers of outcome. Here, we performed whole exome sequencing and 16S ribosomal RNA (rRNA) amplicon sequencing from 44 pre-treatment LARC biopsies from Argentinian and Brazilian patients, treated with neoadjuvant chemoradiotherapy or total neoadjuvant treatment, searching for predictive biomarkers of response (responders, n = 17; non-responders, n = 27). In general, the somatic landscape of LARC was not capable to predict a response; however, a significant enrichment in mutational signature SBS5 was observed in non-responders (p = 0.0021), as well as the co-occurrence of APC and FAT4 mutations (p < 0.05). Microbiota studies revealed a similar alpha and beta diversity of bacteria between response groups. Yet, the linear discriminant analysis (LDA) of effect size indicated an enrichment of Hungatella, Flavonifractor, and Methanosphaera (LDA score ≥3) in the pre-treatment biopsies of responders, while non-responders had a higher abundance of Enhydrobacter, Paraprevotella (LDA score ≥3) and Finegoldia (LDA score ≥4). Altogether, the evaluation of these biomarkers in pre-treatment biopsies could eventually predict a neoadjuvant treatment response, while in post-treatment samples, it could help in guiding non-operative treatment strategies.

Targeting a cell surface vitamin D receptor on tumor-associated macrophages in triple-negative breast cancer.

Triple-negative breast cancer (TNBC) is an aggressive tumor with limited treatment options and poor prognosis. We applied the in vivo phage display technology to isolate peptides homing to the immunosuppressive cellular microenvironment of TNBC as a strategy for non-malignant target discovery. We identified a cyclic peptide (CSSTRESAC) that specifically binds to a vitamin D receptor, protein disulfide-isomerase A3 (PDIA3) expressed on the cell surface of tumor-associated macrophages (TAM), and targets breast cancer in syngeneic TNBC, non-TNBC xenograft, and transgenic mouse models. Systemic administration of CSSTRESAC to TNBC-bearing mice shifted the cytokine profile toward an antitumor immune response and delayed tumor growth. Moreover, CSSTRESAC enabled ligand-directed theranostic delivery to tumors and a mathematical model confirmed our experimental findings. Finally, in silico analysis showed PDIA3-expressing TAM in TNBC patients. This work uncovers a functional interplay between a cell surface vitamin D receptor in TAM and antitumor immune response that could be therapeutically exploited.

A novel MYT1L mutation in a boy with syndromic obesity: Case report and literature review.

BACKGROUND: Pathogenic variants involving the MYT1L gene lead to an autosomal dominant form of syndromic obesity, characterized by polyphagia, intellectual disability/developmental delay, and behavioral problems, and that a characteristic facial phenotype does not seem to be recognizable. METHODS: Trio whole exome sequencing was performed in a 10-year-old Brazilian male presenting polyphagia, severe early-onset obesity, intellectual disability, speech delay, macrocephaly, frontal bossing, telecanthus, strabismus, and hypogenitalism. Additionally, we performed a literature review of patients carrying non-copy number MYT1L variants. RESULTS: A de novo genetic variant not previously reported in MYT1L (NM_015025.4:c.2990C>A) was identified in the proband and classified as pathogenic. From a literature search, 22 further patients carrying non-copy number MYT1L variants were identified, evidencing that although the associated phenotype is quite variable, intellectual disability/developmental and speech delays are always present. Further, most patients have obesity or overweight due to polyphagia. Macrocephaly, strabismus, behavioral problems, and hand/feet malformations are also recurrent features. CONCLUSIONS: We described the first Brazilian case of MYT1L related syndrome and highlighted clinical characteristics based on the literature. Other syndromic forms of obesity such as Prader-Willi, Bardet-Biedl, Börjeson-Forssman-Lehmann, MORM, Cohen, Alstrom, and Kleefstra type 1 syndromes should be considered in the differential diagnosis. Further, although obesity is frequent, it is not an obligatory feature of all carriers of MYT1L mutations.

APOBEC-mediated DNA alterations: A possible new mechanism of carcinogenesis in EBV-positive gastric cancer.

Mechanisms of viral oncogenesis are diverse and include the off-target activity of enzymes expressed by the infected cells, which evolved to target viral genomes for controlling their infection. Among these enzymes, the single-strand DNA editing capability of APOBECs represent a well-conserved viral infection response that can also cause untoward mutations in the host DNA. Here we show, after evaluating somatic single-nucleotide variations and transcriptome data in 240 gastric cancer samples, a positive correlation between APOBEC3s mRNA-expression and the APOBEC-mutation signature, both increased in EBV+ tumors. The correlation was reinforced by the observation of APOBEC mutations preferentially occurring in the genomic loci of the most active transcripts. This EBV infection and APOBEC3 mutation-signature axis were confirmed in a validation cohort of 112 gastric cancer patients. Our findings suggest that APOBEC3 upregulation in EBV+ cancer may boost the mutation load, providing further clues to the mechanisms of EBV-induced gastric carcinogenesis. After further validation, this EBV-APOBEC axis may prove to be a secondary driving force in the mutational evolution of EBV+ gastric tumors, whose consequences in terms of prognosis and treatment implications should be vetted.

A Systems-level Characterization of the Differentiation of Human Embryonic Stem Cells into Mesenchymal Stem Cells.

Mesenchymal stem/stromal cells (MSCs) are self-renewing multipotent cells with regenerative, secretory and immunomodulatory capabilities that are beneficial for the treatment of various diseases. To avoid the issues that come with using tissue-derived MSCs in therapy, MSCs may be generated by the differentiation of human embryonic stems cells (hESCs) in culture. However, the changes that occur during the differentiation process have not been comprehensively characterized. Here, we combined transcriptome, proteome and phosphoproteome profiling to perform an in-depth, multi-omics study of the hESCs-to-MSCs differentiation process. Based on RNA-to-protein correlation, we determined a set of high confidence genes that are important to differentiation. Among the earliest and strongest induced proteins with extensive differential phosphorylation was AHNAK, which we hypothesized to be a defining factor in MSC biology. We observed two distinct expression waves of developmental HOX genes and an AGO2-to-AGO3 switch in gene silencing. Exploring the kinetic of noncoding ORFs during differentiation, we mapped new functions to well annotated long noncoding RNAs (CARMN, MALAT, NEAT1, LINC00152) as well as new candidates which we identified to be important to the differentiation process. Phosphoproteome analysis revealed ESC and MSC-specific phosphorylation motifs with PAK2 and RAF1 as top predicted upstream kinases in MSCs. Our data represent a rich systems-level resource on ESC-to-MSC differentiation that will be useful for the study of stem cell biology.

Genomics and epidemiology for gastric adenocarcinomas (GE4GAC): a Brazilian initiative to study gastric cancer

Gastric cancer (GC) is the fifth most common type of cancer worldwide with high incidences in Asia, Central, and South American countries. This patchy distribution means that GC studies are neglected by large research centers from developed countries. The need for further understanding of this complex disease, including the local importance of epidemiological factors and the rich ancestral admixture found in Brazil, stimulated the implementation of the GE4GAC project. GE4GAC aims to embrace epidemiological, clinical, molecular and microbiological data from Brazilian controls and patients with malignant and pre-malignant gastric disease. In this letter, we summarize the main goals of the project, including subject and sample accrual and current findings

A total transcriptome profiling method for plasma-derived extracellular vesicles: applications for liquid biopsies.

Extracellular vesicles (EVs) are key mediators of intercellular communication. Part of their biological effects can be attributed to the transfer of cargos of diverse types of RNAs, which are promising diagnostic and prognostic biomarkers. EVs found in human biofluids are a valuable source for the development of minimally invasive assays. However, the total transcriptional landscape of EVs is still largely unknown. Here we develop a new method for total transcriptome profiling of plasma-derived EVs by next generation sequencing (NGS) from limited quantities of patient-derived clinical samples, which enables the unbiased characterization of the complete RNA cargo, including both small- and long-RNAs, in a single library preparation step. This approach was applied to RNA extracted from EVs isolated by ultracentrifugation from the plasma of five healthy volunteers. Among the most abundant RNAs identified we found small RNAs such as tRNAs, miRNAs and miscellaneous RNAs, which have largely unknown functions. We also identified protein-coding and long noncoding transcripts, as well as circular RNA species that were also experimentally validated. This method enables, for the first time, the full spectrum of transcriptome data to be obtained from minute patient-derived samples, and will therefore potentially allow the identification of cell-to-cell communication mechanisms and biomarkers.

signeR: an empirical Bayesian approach to mutational signature discovery.

MOTIVATION: Mutational signatures can be used to understand cancer origins and provide a unique opportunity to group tumor types that share the same origins and result from similar processes. These signatures have been identified from high throughput sequencing data generated from cancer genomes by using non-negative matrix factorisation (NMF) techniques. Current methods based on optimization techniques are strongly sensitive to initial conditions due to high dimensionality and nonconvexity of the NMF paradigm. In this context, an important question consists in the determination of the actual number of signatures that best represent the data. The extraction of mutational signatures from high-throughput data still remains a daunting task. RESULTS: Here we present a new method for the statistical estimation of mutational signatures based on an empirical Bayesian treatment of the NMF model. While requiring minimal intervention from the user, our method addresses the determination of the number of signatures directly as a model selection problem. In addition, we introduce two new concepts of significant clinical relevance for evaluating the mutational profile. The advantages brought by our approach are shown by the analysis of real and synthetic data. The later is used to compare our approach against two alternative methods mostly used in the literature and with the same NMF parametrization as the one considered here. Our approach is robust to initial conditions and more accurate than competing alternatives. It also estimates the correct number of signatures even when other methods fail. Results on real data agree well with current knowledge. AVAILABILITY AND IMPLEMENTATION: signeR is implemented in R and C ++, and is available as a R package at http://bioconductor.org/packages/signeR CONTACT: itojal@cipe.accamargo.org.brSupplementary information: Supplementary data are available at Bioinformatics online.

Tissue-Associated Bacterial Alterations in Rectal Carcinoma Patients Revealed by 16S rRNA Community Profiling.

Sporadic and inflammatory forms of colorectal cancer (CRC) account for more than 80% of cases. Recent publications have shown mechanistic evidence for the involvement of gut bacteria in the development of both CRC-forms. Whereas, colon and rectal cancer have been routinely studied together as CRC, increasing evidence show these to be distinct diseases. Also, the common use of fecal samples to study microbial communities may reflect disease state but possibly not the tumor microenvironment. We performed this study to evaluate differences in bacterial communities found in tissue samples of 18 rectal-cancer subjects when compared to 18 non-cancer controls. Samples were collected during exploratory colonoscopy (non-cancer group) or during surgery for tumor excision (rectal-cancer group). High throughput 16S rRNA amplicon sequencing of the V4-V5 region was conducted on the Ion PGM platform, reads were filtered using Qiime and clustered using UPARSE. We observed significant increases in species richness and diversity in rectal cancer samples, evidenced by the total number of OTUs and the Shannon and Simpson indexes. Enterotyping analysis divided our cohort into two groups, with the majority of rectal cancer samples clustering into one enterotype, characterized by a greater abundance of Bacteroides and Dorea. At the phylum level, rectal-cancer samples had increased abundance of candidate phylum OD1 (also known as Parcubacteria) whilst non-cancer samples had increased abundance of Planctomycetes. At the genera level, rectal-cancer samples had higher abundances of Bacteroides, Phascolarctobacterium, Parabacteroides, Desulfovibrio, and Odoribacter whereas non-cancer samples had higher abundances of Pseudomonas, Escherichia, Acinetobacter, Lactobacillus, and Bacillus. Two Bacteroides fragilis OTUs were more abundant among rectal-cancer patients seen through 16S rRNA amplicon sequencing, whose presence was confirmed by immunohistochemistry and enrichment verified by digital droplet PCR. Our findings point to increased bacterial richness and diversity in rectal cancer, along with several differences in microbial community composition. Our work is the first to present evidence for a possible role of bacteria such as B. fragilis and the phylum Parcubacteria in rectal cancer, emphasizing the need to study tissue-associated bacteria and specific regions of the gastrointestinal tract in order to better understand the possible links between the microbiota and rectal cancer.

Plasmodium Infection Promotes Genomic Instability and AID-Dependent B Cell Lymphoma.

Chronic infection with Plasmodium falciparum was epidemiologically associated with endemic Burkitt's lymphoma, a mature B cell cancer characterized by chromosome translocation between the c-myc oncogene and Igh, over 50 years ago. Whether infection promotes B cell lymphoma, and if so by which mechanism, remains unknown. To investigate the relationship between parasitic disease and lymphomagenesis, we used Plasmodium chabaudi (Pc) to produce chronic malaria infection in mice. Pc induces prolonged expansion of germinal centers (GCs), unique compartments in which B cells undergo rapid clonal expansion and express activation-induced cytidine deaminase (AID), a DNA mutator. GC B cells elicited during Pc infection suffer widespread DNA damage, leading to chromosome translocations. Although infection does not change the overall rate, it modifies lymphomagenesis to favor mature B cell lymphomas that are AID dependent and show chromosome translocations. Thus, malaria infection favors mature B cell cancers by eliciting protracted AID expression in GC B cells. PAPERCLIP.

Epigenetic targeting of activation-induced cytidine deaminase.

Activation-induced cytidine deaminase (AID) initiates class switch recombination (CSR) and somatic hypermutation (SHM) by deaminating cytosine residues in immunoglobulin genes (Igh, Igκ, and Igλ). At a lower frequency, AID also causes collateral DNA damage at non-Ig loci, including genes that are rearranged or mutated in B-cell lymphoma. Precisely how AID is recruited to these off-target sites is not entirely understood. To gain further insight into how AID selects its targets, we compared AID-mediated translocations in two different cell types, B cells and mouse embryonic fibroblasts (MEFs). AID targets a distinct set of hotspots in the two cell types. In both cases, hotspots are concentrated in highly transcribed but stalled genes. However, transcription alone is insufficient to recruit AID activity. Comparison of genes similarly transcribed in B cells and MEFs but targeted in only one of the two cell types reveals a common set of epigenetic features associated with AID recruitment in both cells. AID target genes are enriched in chromatin modifications associated with active enhancers (such as H3K27Ac) and marks of active transcription (such as H3K36me3) in both fibroblasts and B cells, indicating that these features are universal mediators of AID recruitment.

Identification of chromosomal translocation hotspots via scan statistics.

MOTIVATION: The detection of genomic regions unusually rich in a given pattern is an important undertaking in the analysis of next-generation sequencing data. Recent studies of chromosomal translocations in activated B lymphocytes have identified regions that are frequently translocated to c-myc oncogene. A quantitative method for the identification of translocation hotspots was crucial to this study. Here we improve this analysis by using a simple probabilistic model and the framework provided by scan statistics to define the number and location of translocation breakpoint hotspots. A key feature of our method is that it provides a global chromosome-wide nominal control level to clustering, as opposed to previous methods based on local criteria. While being motivated by a specific application, the detection of unusual clusters is a widespread problem in bioinformatics. We expect our method to be useful in the analysis of data from other experimental approaches such as of ChIP-seq and 4C-seq. RESULTS: The analysis of translocations from B lymphocytes with the method described here reveals the presence of longer hotspots when compared with those defined previously. Further, we show that the hotspot size changes substantially in the absence of DNA repair protein 53BP1. When 53BP1 deficiency is combined with overexpression of activation-induced cytidine deaminase, the hotspot length increases even further. These changes are not detected by previous methods that use local significance criteria for clustering. Our method is also able to identify several exclusive translocation hotspots located in genes of known tumor supressors. AVAILABILITY AND IMPLEMENTATION: The detection of translocation hotspots is done with hot_scan, a program implemented in R and Perl. Source code and documentation are freely available for download at https://github.com/itojal/hot_scan.

Altered expression of degranulation-related genes in CD8+ T cells in human T lymphotropic virus type I infection.

Human T lymphotropic virus type I (HTLV-1) is the etiological agent of HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM/TSP). CD8+ T cells may contribute to the protection or development of HAM/TSP. In this study we used SAGE methodology to screen for differentially expressed genes in CD8+ T cells isolated from HTLV-1 asymptomatic carriers (HAC) and from HAM/TSP patients to identify genes involved in HAM/TSP development. SAGE analysis was conducted by pooling samples according to clinical status. The comparison of gene expression profiles between HAC and HAM/TSP libraries identified 285 differentially expressed tags. We focus on cytotoxicity and cytokine-related genes due to their potential biological role in HTLV-1 infection. Our results showed that patients with HAM/TSP have high expression levels of degranulation-related genes, namely GZMH and PRF1, and of the cytoskeletal adaptor PXN. We found that GZMB and ZAP70 were overexpressed in HTLV-infected patients compared to the noninfected group. We also detected that CCL5 was higher in the HAM/TSP group compared to the HAC and CT groups. Our findings showed that CD8+ T cells of HAM/TSP patients have an inflammatory and active profile. PXN and ZAP70 overexpression in HTLV-1-infected patients was described for the first time here and reinforces this concept. However, although active and abundant, CD8+ T cells are not able to completely eliminate infected cells and prevent the development of HAM/TSP and, moreover, these cells might contribute to the pathogenesis of the disease by migrating to the central nervous system (CNS). These results should be further tested with biological functional assays to increase our understanding on the role of these molecules in the development of HTLV-1-related diseases.

53BP1 alters the landscape of DNA rearrangements and suppresses AID-induced B cell lymphoma.

Deficiencies in factors that regulate the DNA damage response enhance the incidence of malignancy by destabilizing the genome. However, the precise influence of the DNA damage response on regulation of cancer-associated rearrangements is not well defined. Here we examine the genome-wide impact of tumor protein P53-binding protein 1 (53BP1) deficiency in lymphoma and translocation. While both activation-induced cytidine deaminase (AID) and 53BP1 have been associated with cancer in humans, neither AID overexpression nor loss of 53BP1 is sufficient to produce malignancy. However, the combination of 53BP1 deficiency and AID deregulation results in B cell lymphoma. Deep sequencing of the genome of 53BP1(-/-) cancer cells and translocation capture sequencing (TC-Seq) of primary 53BP1(-/-) B cells revealed that their chromosomal rearrangements differ from those found in wild-type cells in that they show increased DNA end resection. Moreover, loss of 53BP1 alters the translocatome by increasing rearrangements to intergenic regions.

ProbFAST: Probabilistic functional analysis system tool.

BACKGROUND: The post-genomic era has brought new challenges regarding the understanding of the organization and function of the human genome. Many of these challenges are centered on the meaning of differential gene regulation under distinct biological conditions and can be performed by analyzing the Multiple Differential Expression (MDE) of genes associated with normal and abnormal biological processes. Currently MDE analyses are limited to usual methods of differential expression initially designed for paired analysis. RESULTS: We proposed a web platform named ProbFAST for MDE analysis which uses Bayesian inference to identify key genes that are intuitively prioritized by means of probabilities. A simulated study revealed that our method gives a better performance when compared to other approaches and when applied to public expression data, we demonstrated its flexibility to obtain relevant genes biologically associated with normal and abnormal biological processes. CONCLUSIONS: ProbFAST is a free accessible web-based application that enables MDE analysis on a global scale. It offers an efficient methodological approach for MDE analysis of a set of genes that are turned on and off related to functional information during the evolution of a tumor or tissue differentiation. ProbFAST server can be accessed at http://gdm.fmrp.usp.br/probfast.