Search

BackgroundTransposable elements (TEs) are major components of plant genomes. Despite being regarded as “junk DNA” at first, TEs play important roles for the organisms they are found in. The most obvious and easily recognizable effects caused by TEs result from their mobility, which can disrupt coding sequences or promoter regions. However, with the recent advances in transcriptomics, it is becoming increasingly evident that TEs can act as an additional layer of gene expression regulation through a number of processes, which can involve production of non-coding RNAs. Here, we describe how Tnt1, a stress-responsive LTR-retrotransposon, interferes with gene expression and modulate a number of developmental aspects in tobacco.ResultsThrough an RNAi approach, we generated tobacco (HP) lines knocked-down for Tnt1 expression. Quantitative RT-PCR experiments confirm that Tnt1 is downregulated in HP lines after ethylene exposure. A RNA-seq experiment was performed and through two independent bioinformatic approaches (with different stringencies) we found 932 and 97 differentially expressed genes in HP lines. A number of phenotypes were observed in such lines, namely lesion mimicry in leaves, underdevelopment of the root system, overproduction of root hairs and early loss of seed viability. Folding prediction of part of the Tnt1 mRNA reveals putative stem-loop secondary structures containing transcriptional regulation sequences, suggesting it could be a source of small RNAs. We also propose a model to explain the Tnt1 expression in both homeostatic and stress conditions, and how it could interact with stress-responsive genes.ConclusionsOur results are consistent that interferences with Tnt1 transcript levels correlate with transcriptomic and phenotypic changes, suggesting a functional role for this element during plant development and stress response.

Alterations in microRNA (miRNA) processing have been previously linked to aging. Here we used the small molecule enoxacin to pharmacologically interfere with miRNA biogenesis and study how it affects aging in C. elegans. Enoxacin extended worm lifespan and promoted survival under normal and oxidative stress conditions. Enoxacin-induced longevity required the transcription factor SKN-1/Nrf2 and was blunted by the antioxidant N-acetyl-cysteine, suggesting a prooxidant-mediated mitohormetic response. The longevity effects of enoxacin were also dependent on the miRNA pathway, consistent with changes in miRNA expression elicited by the drug. Among these differentially expressed miRNAs, the widely conserved miR-34-5p was found to play an important role in enoxacin-mediated longevity. Enoxacin treatment down-regulated miR-34-5p and did not further extend lifespan of long-lived mir-34 mutants. Moreover, N-acetyl-cysteine abrogated mir-34(gk437)-induced longevity. Evidence also points to double-stranded RNA-specific adenosine deaminases (ADARs) as new targets of enoxacin since ADAR loss-of-function abrogates enoxacin-induced lifespan extension. Thus, enoxacin increases lifespan by reducing miR-34-5p levels, interfering with the redox balance and promoting healthspan., Highlights • Enoxacin extends lifespan and promotes healthspan in C. elegans. • Enoxacin downregulates miR-34-5p to prolong lifespan. • Longevity induced by enoxacin and mir-34 loss of function involves a mitohormetic response. • Double-stranded RNA-specific adenosine deaminases (ADARs) are potential targets of enoxacin.

The risk of developing metastatic disease in breast cancer patients is traditionally predictable based on the number of positive axillary lymph nodes, complemented with additional clinicopathological factors. However, since lymph node-negative patients have a 20–30% probability of developing metastatic disease, lymph node information alone is insufficient to accurately assess individual risk. Molecular approaches, such as multigene expression panels, analyze a set of cancer-related genes that more accurately predict the early risk of metastasis and the treatment response. Here, we present N-Myc downstream-regulated gene 4 (NDRG4) epigenetic silencing as a mechanistic biomarker of metastasis in ductal invasive breast tumors. While aberrant NDRG4 DNA hypermethylation is significantly associated with the development of metastatic disease, downregulation of NDRG4 transcription and protein expression is functionally associated with enhanced lymph node adhesion and cell mobility. Here, we show that epigenetic silencing of NDRG4 modulates integrin signaling by assembling β1-integrins into large punctate clusters at the leading edge of tumor cells to promote an “adhesive switch,” decreasing cell adhesion to fibronectin and increasing cell adhesion and migration towards vitronectin, an important component of human lymph nodes. Taken together, our functional and clinical observations suggest that NDRG4 is a potential mechanistic biomarker in breast cancer that is functionally associated with metastatic disease., Biomarker: Methylated gene promoter increases risk of metastatic disease Epigenetic silencing of a putative tumor suppressor protein confers a metastatic advantage to breast cancer cells and could serve as a mechanistic biomarker of likely disease spread. Researchers in São Paulo, Brazil, led by Érico Costa from the Molecular Oncology Center of the Hospital Sírio-Libanês, showed that NDRG4 (short for N-Myc downstream-regulated gene 4) is actively expressed in normal breast tissue. But in breast tumors and in tumor cell lines, methyl tags cover the promoter region of the NDRG4-encoding gene, resulting in reduced protein expression, increased tumor size, elevated cell mobility and other molecular features indicative of a more aggressive disease state. The researchers experimentally silenced NDRG4 expression in two non-metastatic breast cancer cell lines and observed changes in the adhesive and migratory properties of cells that promote tumor invasiveness and spread.

Butyrate is a short-chain fatty acid derived from the metabolism of indigestible carbohydrates by the gut microbiota. Butyrate contributes to gut homeostasis, but it may also control inflammatory responses and host physiology in other tissues. Butyrate inhibits histone deacetylases, thereby affecting gene transcription, and also signals through the metabolite-sensing G protein receptor (GPR)109a. We produced an mAb to mouse GPR109a and found high expression on podocytes in the kidney. Wild-type and Gpr109a-/- mice were induced to develop nephropathy by a single injection of Adriamycin and treated with sodium butyrate or high butyrate-releasing high-amylose maize starch diet. Butyrate improved proteinuria by preserving podocyte at glomerular basement membrane and attenuated glomerulosclerosis and tissue inflammation. This protective phenotype was associated with increased podocyte-related proteins and a normalized pattern of acetylation and methylation at promoter sites of genes essential for podocyte function. We found that GPR109a is expressed by podocytes, and the use of Gpr109a-/- mice showed that the protective effects of butyrate depended on GPR109a expression. A prebiotic diet that releases high amounts of butyrate also proved highly effective for protection against kidney disease. Butyrate and GPR109a play a role in the pathogenesis of kidney disease and provide one of the important molecular connections between diet, the gut microbiota, and kidney disease.-Felizardo, R. J. F., de Almeida, D. C., Pereira, R. L., Watanabe, I. K. M., Doimo, N. T. S., Ribeiro, W. R., Cenedeze, M. A., Hiyane, M. I., Amano, M. T., Braga, T. T., Ferreira, C. M., Parmigiani, R. B., Andrade-Oliveira, V., Volpini, R. A., Vinolo, M. A. R., Marino, E., Robert, R., Mackay, C. R., Camara, N. O. S. Gut microbial metabolite butyrate protects against proteinuric kidney disease through epigenetic- and GPR109a-mediated mechanisms.

The phenotype of primary cells in culture varies according to the donor environmental condition. We recently showed that the time of the day imposes a molecular program linked to the inflammatory response that is heritable in culture. Here we investigated whether microRNAs (miRNAs) would show differential expression according to the time when cells were obtained, namely daytime or nighttime. Cells obtained from explants of cremaster muscle and cultivated until confluence (∼20 days) presented high CD133 expression. Global miRNA expression analysis was performed through deep sequencing in order to compare both cultured cells. A total of 504 mature miRNAs were identified, with a specific miRNA signature being associated to the light versus dark phase of a circadian cycle. miR-1249 and miR-129-2-3p were highly expressed in daytime cells, while miR-182, miR-96-5p, miR-146a-3p, miR-146a-5p, and miR-223-3p were highly expressed in nighttime cells. Nighttime cells are regulated for programs involved in cell processes and development, as well as in the inflammation, cell differentiation and maturation; while daytime cells express miRNAs that control stemness and cytoskeleton remodeling. In summary, the time of the day imposes a differential profile regarding to miRNA signature on CD133(+) cells in culture. Understanding this daily profile in the phenotype of cultured cells is highly relevant for clinical outputs, including cellular therapy approaches. J. Cell. Physiol. 231: 1953-1963, 2016. © 2016 Wiley Periodicals, Inc.

// Paola Carpinetti 1, 2, 3 , Elisa Donnard 1, 2, 3 , Fabiana Bettoni 2 , Paula Asprino 2 , Fernanda Koyama 1 , Andrei Rozanski 2 , Jorge Sabbaga 4 , Angelita Habr-Gama 5, 6 , Raphael B. Parmigiani 2 , Pedro A.F. Galante 2 , Rodrigo O. Perez 1, 5, 6 , Anamaria A. Camargo 1, 2 1 Ludwig Institute for Cancer Research, Sao Paulo, SP, Brazil 2 Centro de Oncologia Molecular Hospital Sirio Libanes, Sao Paulo, SP, Brazil 3 Programa de Pos Graduacao em Bioquimica, Instituto de Quimica, Universidade de Sao Paulo, SP, Brazil 4 Centro de Oncologia Clinica, Hospital Sirio Libanes, Sao Paulo, SP, Brazil 5 Angelita & Joaquim Gama Institute, Sao Paulo, SP, Brazil 6 University of Sao Paulo, School of Medicine, Sao Paulo, SP, Brazil Correspondence to: Anamaria A. Camargo, e-mail: aacamargo@mochsl.org.br Keywords: rectal cancer, liquid biopsies, ctDNA, personalized biomarkers, neoadjuvant therapy Received: June 19, 2015 Accepted: September 24, 2015 Published: October 06, 2015 ABSTRACT Neoadjuvant chemoradiotherapy (nCRT) followed by surgery is the mainstay treatment for locally advanced rectal cancer. Variable degrees of tumor regression are observed after nCRT and alternative treatment strategies, including close surveillance without immediate surgery, have been investigated to spare patients with complete tumor regression from potentially adverse outcomes of radical surgery. However, clinical and radiological assessment of response does not allow accurate identification of patients with complete response. In addition, surveillance for recurrence is similarly important for these patients, as early detection of recurrence allows salvage resections and adjuvant interventions. We report the use of liquid biopsies and personalized biomarkers for monitoring treatment response to nCRT and detecting residual disease and recurrence in patients with rectal cancer. We sequenced the whole-genome of four rectal tumors to identify patient-specific chromosomal rearrangements that were used to monitor circulating tumor DNA (ctDNA) in liquid biopsies collected at diagnosis and during nCRT and follow-up. We compared ctDNA levels to clinical, radiological and pathological response to nCRT. Our results indicate that personalized biomarkers and liquid biopsies may not be sensitive for the detection of microscopic residual disease. However, it can be efficiently used to monitor treatment response to nCRT and detect disease recurrence, preceding increases in CEA levels and radiological diagnosis. Similar good results were observed when assessing tumor response to systemic therapy and disease progression. Our study supports the use of personalized biomarkers and liquid biopsies to tailor the management of rectal cancer patients, however, replication in a larger cohort is necessary to introduce this strategy into clinical practice.

Background Genetic studies have largely concentrated on the impact of somatic mutations found in coding regions, and have neglected mutations outside of these. However, 3' untranslated regions (3' UTR) mutations can also disrupt or create miRNA target sites, and trigger oncogene activation or tumor suppressor inactivation. Methods We used next-generation sequencing to widely screen for genetic alterations within predicted miRNA target sites of oncogenes associated with colorectal cancer, and evaluated the functional impact of a new somatic mutation. Target sequencing of 47 genes was performed for 29 primary colorectal tumor samples. For 71 independent samples, Sanger methodology was used to screen for E2F1 mutations in miRNA predicted target sites, and the functional impact of these mutations was evaluated by luciferase reporter assays. Results We identified germline and somatic alterations in E2F1. Of the 100 samples evaluated, 3 had germline alterations at the MIR205-5p target site, while one had a somatic mutation at MIR136-5p target site. E2F1 gene expression was similar between normal and tumor tissues bearing the germline alteration; however, expression was increased 4-fold in tumor tissue that harbored a somatic mutation compared to that in normal tissue. Luciferase reporter assays revealed both germline and somatic alterations increased E2F1 activity relative to wild-type E2F1. Conclusions We demonstrated that somatic mutation within E2F1:MIR136-5p target site impairs miRNA-mediated regulation and leads to increased gene activity. We conclude that somatic mutations that disrupt miRNA target sites have the potential to impact gene regulation, highlighting an important mechanism of oncogene activation.

We carried out a mutational analysis of 3,594 genes coding for cell surface proteins (Surfaceome) in 23 colorectal cancer cell lines, searching for new altered pathways, druggable mutations and mutated epitopes for targeted therapy in colorectal cancer. A total of 3,944 somatic non-synonymous substitutions and 595 InDels, occurring in 2,061 (57%) Surfaceome genes were catalogued. We identified 48 genes not previously described as mutated in colorectal tumors in the TCGA database, including genes that are mutated and expressed in >10% of the cell lines (SEMA4C, FGFRL1, PKD1, FAM38A, WDR81, TMEM136, SLC36A1, SLC26A6, IGFLR1). Analysis of these genes uncovered important roles for FGF and SEMA4 signaling in colorectal cancer with possible therapeutic implications. We also found that cell lines express on average 11 druggable mutations, including frequent mutations (>20%) in the receptor tyrosine kinases AXL and EPHA2, which have not been previously considered as potential targets for colorectal cancer. Finally, we identified 82 cell surface mutated epitopes, however expression of only 30% of these epitopes was detected in our cell lines. Notwithstanding, 92% of these epitopes were expressed in cell lines with the mutator phenotype, opening new venues for the use of “general” immune checkpoint drugs in this subset of patients.

Although patterns of somatic alterations have been reported for tumor genomes, little is known on how they compare with alterations present in non-tumor genomes. A comparison of the two would be crucial to better characterize the genetic alterations driving tumorigenesis. We sequenced the genomes of a lymphoblastoid (HCC1954BL) and a breast tumor (HCC1954) cell line derived from the same patient and compared the somatic alterations present in both. The lymphoblastoid genome presents a comparable number and similar spectrum of nucleotide substitutions to that found in the tumor genome. However, a significant difference in the ratio of non-synonymous to synonymous substitutions was observed between both genomes (P = 0.031). Protein–protein interaction analysis revealed that mutations in the tumor genome preferentially affect hub-genes (P = 0.0017) and are co-selected to present synergistic functions (P

Cancer/testis Antigens (CTAs) are immunogenic proteins with a restricted expression pattern in normal tissues and aberrant expression in different types of tumors being considered promising candidates for immunotherapy. We used the alignment between EST sequences and the human genome sequence to identify novel CT genes. By examining the EST tissue composition of known CT clusters we defined parameters for the selection of 1184 EST clusters corresponding to putative CT genes. The expression pattern of 70 CT gene candidates was evaluated by RT-PCR in 21 normal tissues, 17 tumor cell lines and 160 primary tumors. We were able to identify 4 CT genes expressed in different types of tumors. The presence of antibodies against the protein encoded by 1 of these 4 CT genes (FAM46D) was exclusively detected in plasma samples from cancer patients. Due to its restricted expression pattern and immunogenicity FAM46D represents a novel target for cancer immunotherapy.

Cancer-testis (CT) antigens are immunogenic proteins expressed in normal gametogenic tissues and in different types of tumors. CTSP-1 is a CT antigen frequently expressed in prostate tumors, and capable of eliciting humoral response in prostate cancer patients. Here, we analyzed the presence of anti-CTSP-1 antibodies in 147 patients with localized prostate cancer and determined its prognostic value for predicting biochemical-recurrence after radical prostatectomy. Anti-CTSP-1 antibodies were detected in 25% of the patients and a significant correlation (p = 0.017) between CTSP-1 protein expression and the presence of specific humoral response was observed. No association was found between the presence of antibodies and the pathological variables analyzed. On univariate analysis, patients without antibodies against CTSP-1 had a lower biochemical-recurrence free survival than did those with anti-CTSP-1 antibodies, although the difference between the groups was not statistically significant (57 vs. 75%, p = 0.075). However, the presence of antibodies against CTSP-1 was significantly associated with a better prognosis in patients with higher Gleason score (36 vs. 80%, p = 0.028). On multivariate analysis, antibodies against CTSP-1 were associated with a better prognosis (Hazard ratio = 0.41, 95% IC 0.18-0.90 p = 0.039), being the third most powerful prognostic factor among Gleason score and preoperative PSA levels. CTSP-1 should be considered a promising candidate for prostate cancer immunotherapy, since it is frequently expressed in prostate tumors and capable of eliciting humoral immune response in prostate cancer patients. Our results also suggest that humoral response against CTSP-1 could be used as a prognostic marker, especially among patients with a high Gleason score.

There are 18 histone deacetylases (HDAC) generally divided into four classes based on homology to yeast HDACs. HDACs have many protein substrates in addition to histones that are involved in regulation of gene expression, cell proliferation, and cell death. Inhibition of HDACs can cause accumulation of acetylated forms of these proteins, thus altering their function. HDAC inhibitors (HDACi), such as the hydroxamic acid–based vorinostat (suberoylanilide hydroxamic acid), inhibit the zinc-containing classes I, II, and IV, but not the NAD+-dependent class III, enzymes. HDACis are a group of novel anticancer agents. Vorinostat is the first HDACi approved for clinical use in the treatment of the cancer cutaneous T-cell lymphoma. Factors affecting expression of HDACs are not well understood. This study focuses on the effect of the HDACi vorinostat on the expression of class I and class II HDACs. We found that vorinostat selectively down-regulates HDAC7 with little or no effect on the expression of other class I or class II HDACs. Fourteen cell lines were examined, including normal, immortalized, genetically transformed, and human cancer-derived cell lines. Down-regulation of HDAC7 by vorinostat is more pronounced in transformed cells sensitive to inhibitor-induced cell death than in normal cells or cancer cells resistant to induced cell death. Modulation of HDAC7 levels by small interfering RNA–mediated knockdown or by HDAC7 overexpression is associated with growth arrest but without detectable changes in acetylation of histones or p21 gene expression. Selective down-regulation of HDAC7 protein may serve as a marker of response of tumors to HDACi. [Mol Cancer Ther 2007;6(9):2525–34]

This review focuses on the mechanisms of action of histone deacetylase (HDAC) inhibitors (HDACi), a group of recently discovered 'targeted' anticancer agents. There are 18 HDACs, which are generally divided into four classes, based on sequence homology to yeast counterparts. Classical HDACi such as the hydroxamic acid-based vorinostat (also known as SAHA and Zolinza) inhibits classes I, II and IV, but not the NAD+-dependent class III enzymes. In clinical trials, vorinostat has activity against hematologic and solid cancers at doses well tolerated by patients. In addition to histones, HDACs have many other protein substrates involved in regulation of gene expression, cell proliferation and cell death. Inhibition of HDACs causes accumulation of acetylated forms of these proteins, altering their function. Thus, HDACs are more properly called 'lysine deacetylases.' HDACi induces different phenotypes in various transformed cells, including growth arrest, activation of the extrinsic and/or intrinsic apoptotic pathways, autophagic cell death, reactive oxygen species (ROS)-induced cell death, mitotic cell death and senescence. In comparison, normal cells are relatively more resistant to HDACi-induced cell death. The plurality of mechanisms of HDACi-induced cell death reflects both the multiple substrates of HDACs and the heterogeneous patterns of molecular alterations present in different cancer cells.

The phenotype of primary cells in culture varies according to the donor environmental condition. We recently showed that the time of the day imposes a molecular program linked to the inflammatory response that is heritable in culture. Here we investigated whether microRNAs (miRNAs) would show differential expression according to the time when cells were obtained, namely daytime or nighttime. Cells obtained from explants of cremaster muscle and cultivated until confluence (∼20 days) presented high CD133 expression. Global miRNA expression analysis was performed through deep sequencing in order to compare both cultured cells. A total of 504 mature miRNAs were identified, with a specific miRNA signature being associated to the light versus dark phase of a circadian cycle. miR-1249 and miR-129-2-3p were highly expressed in daytime cells, while miR-182, miR-96-5p, miR-146a-3p, miR-146a-5p, and miR-223-3p were highly expressed in nighttime cells. Nighttime cells are regulated for programs involved in cell processes and development, as well as in the inflammation, cell differentiation and maturation; while daytime cells express miRNAs that control stemness and cytoskeleton remodeling. In summary, the time of the day imposes a differential profile regarding to miRNA signature on CD133(+) cells in culture. Understanding this daily profile in the phenotype of cultured cells is highly relevant for clinical outputs, including cellular therapy approaches. J. Cell. Physiol. 231: 1953-1963, 2016. © 2016 Wiley Periodicals, Inc.

Neoadjuvant chemoradiotherapy (nCRT) may lead to complete tumor regression in rectal cancer patients. Prediction of complete response to nCRT may allow a personalized management of rectal cancer and spare patients from unnecessary radical total mesorectal excision with or without sphincter preservation. To identify a gene expression signature capable of predicting complete pathological response (pCR) to nCRT, we performed a gene expression analysis in 25 pretreatment biopsies from patients who underwent 5FU-based nCRT using RNA-Seq. A supervised learning algorithm was used to identify expression signatures capable of predicting pCR, and the predictive value of these signatures was validated using independent samples. We also evaluated the utility of previously published signatures in predicting complete response in our cohort. We identified 27 differentially expressed genes between patients with pCR and patients with incomplete responses to nCRT. Predictive gene signatures using subsets of these 27 differentially expressed genes peaked at 81.8% accuracy. However, signatures with the highest sensitivity showed poor specificity, and vice-versa, when applied in an independent set of patients. Testing previously published signatures on our cohort also showed poor predictive value. Our results indicate that currently available predictive signatures are highly dependent on the sample set from which they are derived, and their accuracy is not superior to current imaging and clinical parameters used to assess response to nCRT and guide surgical intervention.

Cancer/testis (CT) antigens are immunogenic proteins expressed in normal gametogenic tissues and in different types of tumors. CT antigens are promising candidates for cancer immunotherapy, and the identification of novel CT antigens is a prerequisite for the development of cancer vaccines. We have identified a CT antigen, named CTSP-1, with partial similarity to the breast differentiation antigen NY-BR-1. CTSP-1 presents several splicing and polyadenylation variants and has a very restricted expression pattern among normal tissues. CTSP-1 is exclusively expressed in normal testis and is aberrantly expressed in 47.6% (10 of 21) of tumor cell lines and in 44.4% (75 of 169) of tumors from different histological types. The highest percentages of positive expression were observed in melanomas (59.0%) followed by prostate (58.0%) and lung (57.0%) tumors. CTSP-1 is part of a highly conserved gene family, and members of this family also have a restricted expression pattern and similar protein structure. Antibodies against members of this gene family were detected in 10% (14 of 141) of plasma samples from patients with a wide spectrum of tumors. The highest percentages of antibody response were observed in patients with prostate (20.8%), thyroid (20.0%), and breast (16.6%) tumors. Because of its very restricted expression pattern in normal tissues and immunogenicity in different types of tumors, CTSP-1 should be considered a promising candidate for cancer immunotherapy.

Background Neoadjuvant chemoradiotherapy (nCRT) followed by radical surgery is the preferred treatment strategy for locally advanced rectal cancer. However, complete tumor regression is observed in a significant proportion of patients after nCRT, making them ideal candidates for alternative treatment strategies to this considerably morbid procedure. Identification of such patients based on clinical findings (complete clinical response - cCR) is difficult mainly because it relies on subjective clinical and imaging studies. Our goal was to identify biomarkers capable of predicting complete response to nCRT. Methods We analyzed miRNA expression profile using deep sequencing in rectal tumor biopsies prior to nCRT. Differential expression was investigated by EdgeR for a training (n = 27) and a validation (n = 16) set of patients to identify miRNAs associated with treatment response (complete vs. incomplete). In vitro experiments with two cancer cell lines were also performed in order to evaluate the possible role of miRNAs on response to nCRT. Results We found 4 miRNAs differentially expressed between complete and incomplete responders to nCRT. In addition, validation was performed using an independent group of patients and miR-21-5p was confirmed as being overexpressed in complete responders. Overall sensitivity and specificity of miR-21-5p expression in predicting complete response to nCRT was 78% and 86% respectively. Interestingly, in a subset of patients with cCR followed by early local recurrence, the expression level of miR-21-5p was considerably low, similarly to incomplete responders. We also found SATB1, a miR-21-5p target gene and known multidrug resistance gene, whose expression was inversely correlated with miR-21-5p expression. Finally, we performed functional experiments and showed that miR-21-5p and SATB1 may be directly involved with poor response to nCRT in rectal cancer patients. Conclusions This study suggests miR-21-5p as a promising predictive biomarker, which should aid in the selection of patients with cCR to nCRT that potentially could be spared from radical surgery. Electronic supplementary material The online version of this article (doi:10.1186/s12920-014-0068-7) contains supplementary material, which is available to authorized users.

The era of whole-genome sequencing has revealed that gene copy-number changes caused by duplication and deletion events have important evolutionary, functional, and phenotypic consequences. Recent studies have therefore focused on revealing the extent of variation in copy-number within natural populations of humans and other species. These studies have found a large number of copy-number variants (CNVs) in humans, many of which have been shown to have clinical or evolutionary importance. For the most part, these studies have failed to detect an important class of gene copy-number polymorphism: gene duplications caused by retrotransposition, which result in a new intron-less copy of the parental gene being inserted into a random location in the genome. Here we describe a computational approach leveraging next-generation sequence data to detect gene copy-number variants caused by retrotransposition (retroCNVs), and we report the first genome-wide analysis of these variants in humans. We find that retroCNVs account for a substantial fraction of gene copy-number differences between any two individuals. Moreover, we show that these variants may often result in expressed chimeric transcripts, underscoring their potential for the evolution of novel gene functions. By locating the insertion sites of these duplicates, we are able to show that retroCNVs have had an important role in recent human adaptation, and we also uncover evidence that positive selection may currently be driving multiple retroCNVs toward fixation. Together these findings imply that retroCNVs are an especially important class of polymorphism, and that future studies of copy-number variation should search for these variants in order to illuminate their potential evolutionary and functional relevance., Author Summary Recent studies of human genetic variation have revealed that, in addition to differing at single nucleotide polymorphisms, individuals differ in copy-number at many regions of the genome. These copy-number variants (CNVs) are caused by duplication or deletion events and often affect functional sequences such as genes. Efforts to reveal the functional impact of CNVs have identified many variants increasing the risk of various disorders, and some that are adaptive. However, these studies mostly fail to detect gene duplications caused by retrotransposition, in which an mRNA transcript is reverse-transcribed and reinserted into the genome, yielding a new intron-less gene copy. Here we describe a method leveraging next-generation sequence data to accurately detect gene copy-number variants caused by retrotransposition, or retroCNVs, and apply this method to hundreds of whole-genome sequences from three different human subpopulations. We find that these variants account for a substantial number of gene copy-number differences between individuals, and that gene retrotransposition may often result in both deleterious and beneficial mutations. Indeed, we present evidence that two of these new gene duplications may be adaptive. These results imply that retroCNVs are an especially important class of CNV and should be included in future studies of human copy-number variation.

Eighteen histone deacetylases (HDACs) are present in humans, categorized into two groups: zinc-dependent enzymes (HDAC1–11) and NAD + -dependent enzymes (sirtuins 1–7). Among zinc-dependent HDACs, HDAC6 is unique. It has a cytoplasmic localization, two catalytic sites, a ubiquitin-binding site, and it selectively deacetylases α-tubulin and Hsp90. Here, we report the discovery that the redox regulatory proteins, peroxiredoxin (Prx) I and Prx II are specific targets of HDAC6. Prx are antioxidants enzymes whose main function is H 2 O 2 reduction. Prx are elevated in many cancers and neurodegenerative diseases. The acetylated form of Prx accumulates in the absence of an active HDAC6. Acetylation of Prx increases its reducing activity, its resistance to superoxidation, and its resistance to transition to high-molecular-mass complexes. Thus, HDAC6 and Prx are targets for modulating intracellular redox status in therapeutic strategies for disorders as disparate as cancers and neurodegenerative diseases.

About half of the known miRNA genes are located within protein-coding host genes, and are thus subject to co-transcription. Accumulating data indicate that this coupling may be an intrinsic mechanism to directly regulate the host gene’s expression, constituting a negative feedback loop. Inevitably, the cell requires a yet largely unknown repertoire of methods to regulate this control mechanism. We propose APA as one possible mechanism by which negative feedback of intronic miRNA on their host genes might be regulated. Using in-silico analyses, we found that host genes that contain seed matching sites for their intronic miRNAs yield longer 32UTRs with more polyadenylation sites. Additionally, the distribution of polyadenylation signals differed significantly between these host genes and host genes of miRNAs that do not contain potential miRNA binding sites. We then transferred these in-silico results to a biological example and investigated the relationship between ZFR and its intronic miRNA miR-579 in a U87 cell line model. We found that ZFR is targeted by its intronic miRNA miR-579 and that alternative polyadenylation allows differential targeting. We additionally used bioinformatics analyses and RNA-Seq to evaluate a potential cross-talk between intronic miRNAs and alternative polyadenylation. CPSF2, a gene previously associated with alternative polyadenylation signal recognition, might be linked to intronic miRNA negative feedback by altering polyadenylation signal utilization.

Alternative splicing is one of the major sources of the large transcriptional diversity found in human cells. Splicing variants have been shown to be associated with features like spreading and progression in several human tumors. Therefore, such variants may be of great importance as both diagnostic and therapeutic tools. Here, by using a set of criteria regarding the expression pattern of splicing variants and statistical analyses, we were able to screen the genome for exons overexpressed in tumors of specific tissues. However, as in other analyses attempting to identify tumor-associated variants, our list of candidates was seriously inflated with cases of genes differentially expressed in tumors. To exclude these cases and increase the probability of finding bona fide regulated splicing variants, we performed a serial analysis of gene expression (SAGE), excluding those genes that were shown to be upregulated in tumors. This allowed us to predict the overexpression of single exons in specific tumors. Our final group of candidates includes 1,386 exons belonging to 638 genes. Experimental validation of a few candidates in normal tissue, tumor cell lines, and patient samples suggests that most of these candidates are indeed tumor-associated exons. Further functional classification of our candidate genes shows that our final list is slightly inflated with cancer-related genes.

G protein-coupled receptors (GPCRs) are involved in a large variety of physiological functions. The number of known members that belong to this large family of receptors has been rapidly increasing. Now, with the availability of the human genome sequence databases, further family members are being identified. We describe the identification of a novel GPCR that shows no significant amino acid identity to any one of the known members of the GPCR superfamily. The gene expression pattern of this receptor is restricted: in normal tissues it is confined to the nervous system and testis, but we also detected gene expression in several tumor types, most notably prostate cancer, suggesting a potential role for this gene as a marker for this disease.

Neoadjuvant chemoradiation (nCRT) is the preferred initial treatment for locally advanced rectal cancer due to improved local disease control, lower toxicity and significant tumor regression. Response to nCRT varies substantially among patients and complete tumor regression can be observed in up to 42% of the patients. To avoid unnecessary postoperative morbidity associated with radical surgery, patients exhibiting significant tumor regression have been referred for conservative surgical approaches or even no immediate surgery with strict follow-up (Wait and Watch approach). However, assessment of treatment response remains a significant challenge due to the subjectivity of the clinical assessment and limitations of available radiological studies. In this context, identification of biomarkers capable of accessing tumor regression and complete response to nCRT would allow accurate selection of patients for alternative treatment strategies without immediate radical surgery after nCRT. Circulating DNA carrying tumor-specific genetic alterations can be found in the cell-free fraction of the blood (circulating tumor DNA - ctDNA) and has been successfully used to monitor tumor dynamics in both hematological and solid tumors. Recent advances in sequencing technologies have enabled the rapid and cost effective identification of genetic alterations in individual tumors. These alterations, in special tumor-specific chromosomal rearrangements, can then be used as personalized biomarkers to monitor treatment response and detect minimal residual disease, overcoming the problem imposed by the absence of a universally recurrent genetic alteration in solid tumors. In the present work, we sequenced the whole genome of 6 rectal tumors using mate-pair reads and identified tumor-specific chromosomal rearrangements (TSCR) that were subsequently validated by Sanger sequencing. Personalized assays were then designed to monitor these TSCR in serum samples collected at different time points during nCRT. TSCRs were detected in all pre-treatment serum samples and in post-treatment serum samples from patients presenting marginal tumor regression after nCRT. Post-treatment samples from patients with significant or complete tumor regression were mainly negative for the presence of TSCR. Our preliminary results suggest that personalized biomarkers can be used in the clinics for accessing tumor regression and defining the best surgical approach in rectal cancer patients treated with neoadjuvant chemoradiation. Citation Information: Mol Cancer Ther 2013;12(11 Suppl):A33. Citation Format: Paola Carpinetti, Elisa Donnard, Paula Asprino, Fabiana Bettoni, Fernanda Koyama, Angelita Habr-Gama, Raphael Parmigiani, Pedro Galante, Rodrigo Perez, Anamaria A. Camargo. The use of personalized biomarkers for accessing tumor regression and defining surgical approach in rectal cancer patients treated with neoadjuvant chemoradiation. [abstract]. In: Proceedings of the AACR-NCI-EORTC International Conference: Molecular Targets and Cancer Therapeutics; 2013 Oct 19-23; Boston, MA. Philadelphia (PA): AACR; Mol Cancer Ther 2013;12(11 Suppl):Abstract nr A33.

Somatically acquired chromosomal rearrangements occur at early stages during tumorigenesis and can be used to indirectly detect tumor cells, serving as highly sensitive and tumor-specific biomarkers. Advances in high-throughput sequencing have allowed the genome-wide identification of patient-specific chromosomal rearrangements to be used as personalized biomarkers to efficiently assess response to treatment, detect residual disease and monitor disease recurrence. However, sequencing and data processing costs still represent major obstacles for the widespread application of personalized biomarkers in oncology. We developed a computational pipeline (ICRmax) for the cost-effective identification of a minimal set of tumor-specific interchromosomal rearrangements (ICRs). We examined ICRmax performance on sequencing data from rectal tumors and simulated data achieving an average accuracy of 68% for ICR identification. ICRmax identifies ICRs from low-coverage sequenced tumors, eliminates the need to sequence a matched normal tissue and significantly reduces the costs that limit the utilization of personalized biomarkers in the clinical setting.

The identification of all human chromosome 21 (HC21) genes is a necessary step in understanding the molecular pathogenesis of trisomy 21 (Down syndrome). The first analysis of the sequence of 21q included 127 previously characterized genes and predicted an additional 98 novel anonymous genes. Recently we evaluated the quality of this annotation by characterizing a set of HC21 open reading frames (C21orfs) identified by mapping spliced expressed sequence tags (ESTs) and predicted genes (PREDs), identified only in silico. This study underscored the limitations of in silico-only gene prediction, as many PREDs were incorrectly predicted. To refine the HC21 annotation, we have developed a reliable algorithm to extract and stringently map sequences that contain bona fide 3' transcript ends to the genome. We then created a specific 21q graphical display allowing an integrated view of the data that incorporates new ESTs as well as features such as CpG islands, repeats, and gene predictions. Using these tools we identified 27 new putative genes. To validate these, we sequenced previously cloned cDNAs and carried out RT-PCR, 5'- and 3'-RACE procedures, and comparative mapping. These approaches substantiated 19 new transcripts, thus increasing the HC21 gene count by 9.5%. These transcripts were likely not previously identified because they are small and encode small proteins. We also identified four transcriptional units that are spliced but contain no obvious open reading frame. The HC21 data presented here further emphasize that current gene prediction algorithms miss a substantial number of transcripts that nevertheless can be identified using a combination of experimental approaches and multiple refined algorithms.

Background To identify potential tumor suppressor genes, genome-wide data from exome and transcriptome sequencing were combined to search for genes with loss of heterozygosity and allele-specific expression. The analysis was conducted on the breast cancer cell line HCC1954, and a lymphoblast cell line from the same individual, HCC1954BL. Results By comparing exome sequences from the two cell lines, we identified loss of heterozygosity events at 403 genes in HCC1954 and at one gene in HCC1954BL. The combination of exome and transcriptome sequence data also revealed 86 and 50 genes with allele specific expression events in HCC1954 and HCC1954BL, which comprise 5.4% and 2.6% of genes surveyed, respectively. Many of these genes identified by loss of heterozygosity and allele-specific expression are known or putative tumor suppressor genes, such as BRCA1, MSH3 and SETX, which participate in DNA repair pathways. Conclusions Our results demonstrate that the combined application of high throughput sequencing to exome and allele-specific transcriptome analysis can reveal genes with known tumor suppressor characteristics, and a shortlist of novel candidates for the study of tumor suppressor activities.