Gut microbial metabolite butyrate protects against proteinuric kidney disease through epigenetic- and GPR109a-mediated mechanisms.

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., Mariño, 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.

Light/Dark Environmental Cycle Imposes a Daily Profile in the Expression of microRNAs in Rat CD133(+) 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.

Gene copy-number polymorphism caused by retrotransposition in humans.

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.

ICRmax: an optimized approach to detect tumor-specific interchromosomal rearrangements for clinical application.

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.

Distinct patterns of somatic alterations in a lymphoblastoid and a tumor genome derived from the same individual.

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 < 0.0001). KEGG analysis showed that in the tumor genome most mutated genes were organized into signaling pathways related to tumorigenesis. No such organization or synergy was observed in the lymphoblastoid genome. Our results indicate that endogenous mutagens and replication errors can generate the overall number of mutations required to drive tumorigenesis and that it is the combination rather than the frequency of mutations that is crucial to complete tumorigenic transformation.

Identification of FAM46D as a novel cancer/testis antigen using EST data and serological analysis.

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.

NDRG4 promoter hypermethylation is a mechanistic biomarker associated with metastatic progression in breast cancer patients.

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.

Antibodies against the cancer-testis antigen CTSP-1 are frequently found in prostate cancer patients and are an independent prognostic factor for biochemical-recurrence.

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.

Histone deacetylase inhibitors selectively suppress expression of HDAC7.

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.

Enoxacin extends lifespan of C. elegans by inhibiting miR-34-5p and promoting mitohormesis.

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.

E2F1 somatic mutation within miRNA target site impairs gene regulation in colorectal cancer.

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.

Alternative polyadenylation allows differential negative feedback of human miRNA miR-579 on its host gene ZFR.

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.

The use of personalized biomarkers and liquid biopsies to monitor treatment response and disease recurrence in locally advanced rectal cancer after neoadjuvant chemoradiation.

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.

Comprehensive evaluation of the effectiveness of gene expression signatures to predict complete response to neoadjuvant chemoradiotherapy and guide surgical intervention in rectal cancer.

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.

A novel human G protein-coupled receptor is over-expressed in prostate cancer.

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.

Identification and complete sequencing of novel human transcripts through the use of mouse orthologs and testis cDNA sequences.

The correct identification of all human genes, and their derived transcripts, has not yet been achieved, and it remains one of the major aims of the worldwide genomics community. Computational programs suggest the existence of 30,000 to 40,000 human genes. However, definitive gene identification can only be achieved by experimental approaches. We used two distinct methodologies, one based on the alignment of mouse orthologous sequences to the human genome, and another based on the construction of a high-quality human testis cDNA library, in an attempt to identify new human transcripts within the human genome sequence. We generated 47 complete human transcript sequences, comprising 27 unannotated and 20 annotated sequences. Eight of these transcripts are variants of previously known genes. These transcripts were characterized according to size, number of exons, and chromosomal localization, and a search for protein domains was undertaken based on their putative open reading frames. In silico expression analysis suggests that some of these transcripts are expressed at low levels and in a restricted set of tissues.

Mutational analysis of genes coding for cell surface proteins in colorectal cancer cell lines reveal novel altered pathways, druggable mutations and mutated epitopes for targeted therapy.

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.

Systematic detection of putative tumor suppressor genes through the combined use of exome and transcriptome sequencing.

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.

Characterization of a cancer/testis (CT) antigen gene family capable of eliciting humoral response in cancer patients.

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.

Nineteen additional unpredicted transcripts from human chromosome 21.

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.