Finding fusion genes resulting from chromosome rearrangement by analyzing the expressed sequence databases

Oncogenic gene fusions in cancer: from biology to therapy

Oncogenic gene fusions occur across a broad range of cancers and are a defining feature of some cancer types. Cancers driven by gene fusion products tend to respond well to targeted therapies, where available; thus, detection of potentially targetable oncogenic fusions is necessary to select optimal treatment. Detection methods include non-sequencing methods, such as fluorescence in situ hybridization and immunohistochemistry, and sequencing methods, such as DNA- and RNA-based next-generation sequencing (NGS). While NGS is an efficient way to analyze multiple genes of interest at once, economic and technical factors may preclude its use in routine care globally, despite several guideline recommendations. The aim of this review is to present a summary of oncogenic gene fusions, with a focus on fusions that affect tyrosine kinase signaling, and to highlight the importance of testing for oncogenic fusions. We present an overview of the identification of oncogenic gene fusions and therapies approved for the treatment of cancers harboring gene fusions, and summarize data regarding treating fusion-positive cancers with no current targeted therapies and clinical studies of fusion-positive cancers. Although treatment options may be limited for patients with rare alterations, healthcare professionals should identify patients most likely to benefit from oncogenic gene fusion testing and initiate the appropriate targeted therapy to achieve optimal treatment outcomes.

FGFR3-TACC3 fusion gene promotes glioblastoma malignant progression through the activation of STAT3 signaling pathway

Objective

The Fibroblast growth factor receptors 3-transforming acidic coiled-coil-containing protein 3 (FGFR3-TACC3, F3-T3) oncogenic fusion gene, identified in malignant tumors such as gliomas and bladder cancer, has been particularly noted in recurrent gliomas where it is considered to drive malignant progression, thus presenting itself as a viable therapeutic target. However, the precise mechanism by which F3-T3 facilitates the malignant progression of glioma is not fully understood.

Methods

Correction analysis of STAT3 and FGFR3 with major glioma mutation types and pan-cancer analysis was conducted using The Cancer Genome Atlas (TCGA) database. A series of phenotypic experiments, including CCK-8, EdU, colony-formation assay, wound healing assay, and transwell assay were conducted to detect the effects of F3-T3 on proliferation, invasion, and migration of glioma cells. The association between F3-T3 and epithelial-mesenchymal transition (EMT) was investigated through enrichment analysis of the E-MTAB-6037 gene chip database and confirmed by western blot. The underling mechanism were further inferred and validated through RNA sequencing, E-MTAB-6037 gene chip data, and western blot. The relationship between p-STAT3 expression and the WHO grade of glioma was evaluated using immunohistochemistry (IHC) and tissue microarray analysis. Furthermore, the results of vivo experiments and IHC has confirmed the impact of F3-T3 on glioma malignant progression and activation of the STAT3 signaling pathway.

Results

The experimental results from this study indicate that F3-T3 accelerates the epithelial-mesenchymal transition (EMT) process in glioma cells, thereby promoting their proliferation, invasion, and migration capabilities. Mechanistically, it was determined through RNA sequencing that the signal transducer and activator of transcription 3 (STAT3) signaling pathway is crucial for the malignant progression of F3-T3. This finding was further supported through follow-up experiments conducted after STAT3 knockdown. The role of the STAT3 pathway in gliomas was also reinforced through bioinformatic analysis and immunohistochemistry (IHC) on tissue microarrays (TMA). Further in vivo experiments corroborated the role of F3-T3 in enhancing glioma growth and progression.

Conclusion

F3-T3 facilitates the proliferation, invasion, migration and EMT of glioma cells, thereby promoting their malignant progression through STAT3 signaling activation. These findings highlight its potential as a therapeutic target for glioma treatment.

Repositioning of centromere-associated repeats during karyotype evolution in Oryzias fishes

The karyotype, which is the number and shape of chromosomes, is a fundamental characteristic of all eukaryotes. Karyotypic changes play an important role in many aspects of evolutionary processes, including speciation. In organisms with monocentric chromosomes, it was previously thought that chromosome number changes were mainly caused by centric fusions and fissions, whereas chromosome shape changes, that is, changes in arm numbers, were mainly due to pericentric inversions. However, recent genomic and cytogenetic studies have revealed examples of alternative cases, such as tandem fusions and centromere repositioning, found in the karyotypic changes within and between species. Here, we employed comparative genomic approaches to investigate whether centromere repositioning occurred during karyotype evolution in medaka fishes. In the medaka family (Adrianichthyidae), the three phylogenetic groups differed substantially in their karyotypes. The Oryzias latipes species group has larger numbers of chromosome arms than the other groups, with most chromosomes being metacentric. The O. javanicus species group has similar numbers of chromosomes to the O. latipes species group, but smaller arm numbers, with most chromosomes being acrocentric. The O. celebensis species group has fewer chromosomes than the other two groups and several large metacentric chromosomes that were likely formed by chromosomal fusions. By comparing the genome assemblies of O. latipes, O. javanicus, and O. celebensis, we found that repositioning of centromere-associated repeats might be more common than simple pericentric inversion. Our results demonstrated that centromere repositioning may play a more important role in karyotype evolution than previously appreciated.

Treatment for a primary multidrug-resistant B-cell acute lymphoblastic leukemia patient carrying a SSBP2-CSF1R fusion gene: a case report

SSBP2-CSF1R is an important biomarker for clinical diagnosis and prognosis of Philadelphia chromosome-like acute lymphoblastic leukemia (Ph-like ALL). This case report presents a pediatric Ph-like ALL patient carrying the SSBP2-CSF1R fusion gene. The patient was resistant to most conventional chemotherapy regimens and to dasatinib, an inhibitor that has been reported to have a therapeutic effect on SSBP2-CSF1R fusion Ph-like ALL, as she remained minimal residual disease (MRD) positive (detection by flow cytometry) and SSBP2-CSF1R fusion gene (detection by RT-PCR) positive after five rounds of such regimens. We thus conducted a large-scale in vitro screening to assess the sensitivity of the patient's leukemic cells to anti-cancer drugs. Based on the susceptibility results, we chose to combine cytarabine, homoharringtonine, dexamethasone, fludarabine, vindesine, and epirubicin for treatment. Clinical results showed that after a course of treatment, both MRD and SSBP2-CSF1R fusion gene turned negative, and there was no recurrence during an 18-month follow-up. In conclusion, our study suggests that the SSBP2-CSF1R fusion gene may be an important biomarker of primary drug resistance in Ph-like ALL, and indicate that the combination of cytarabine, homoharringtonine, dexamethasone, fludarabine, vindesine, and epirubicin can achieve optimal therapeutic results in this category of patients.

Small Molecule Inhibitors as Therapeutic Agents Targeting Oncogenic Fusion Proteins: Current Status and Clinical

Oncogenic fusion proteins, arising from chromosomal rearrangements, have emerged as prominent drivers of tumorigenesis and crucial therapeutic targets in cancer research. In recent years, the potential of small molecular inhibitors in selectively targeting fusion proteins has exhibited significant prospects, offering a novel approach to combat malignancies harboring these aberrant molecular entities. This review provides a comprehensive overview of the current state of small molecular inhibitors as therapeutic agents for oncogenic fusion proteins. We discuss the rationale for targeting fusion proteins, elucidate the mechanism of action of inhibitors, assess the challenges associated with their utilization, and provide a summary of the clinical progress achieved thus far. The objective is to provide the medicinal community with current and pertinent information and to expedite the drug discovery programs in this area.

Chimeric chromosome landscapes of human somatic cell cultures show dependence on stress and regulation of genomic repeats by CGGBP1

Genomes of somatic cells in culture are prone to spontaneous mutations due to errors in replication and DNA repair. Some of these errors, such as chromosomal fusions, are not rectifiable and subject to selection or elimination in growing cultures. Somatic cell cultures are thus expected to generate background levels of potentially stable chromosomal chimeras. A description of the landscape of such spontaneously generated chromosomal chimeras in cultured cells will help understand the factors affecting somatic mosaicism. Here we show that short homology-associated non-homologous chromosomal chimeras occur in normal human fibroblasts and HEK293T cells at genomic repeats. The occurrence of chromosomal chimeras is enhanced by heat stress and depletion of a repeat regulatory protein CGGBP1. We also present evidence of homologous chromosomal chimeras between allelic copies in repeat-rich DNA obtained by methylcytosine immunoprecipitation. The formation of homologous chromosomal chimeras at Alu and L1 repeats increases upon depletion of CGGBP1. Our data are derived from de novo sequencing from three different cell lines under different experimental conditions and our chromosomal chimera detection pipeline is applicable to long as well as short read sequencing platforms. These findings present significant information about the generation, sensitivity and regulation of somatic mosaicism in human cell cultures.

Detection of gene mutations and gene-gene fusions in circulating cell-free DNA of glioblastoma patients: an avenue for clinically relevant diagnostic analysis

Glioblastoma (GBM) is the most common type of glioma and is uniformly fatal. Currently, tumour heterogeneity and mutation acquisition are major impedances for tailoring personalized therapy. We collected blood and tumour tissue samples from 25 GBM patients and 25 blood samples from healthy controls. Cell‐free DNA (cfDNA) was extracted from the plasma of GBM patients and from healthy controls. Tumour DNA was extracted from fresh tumour samples. Extracted DNA was sequenced using a whole‐genome sequencing procedure. We also collected 180 tumour DNA datasets from GBM patients publicly available at the TCGA/PANCANCER project. These data were analysed for mutations and gene–gene fusions that could be potential druggable targets. We found that plasma cfDNA concentrations in GBM patients were significantly elevated (22.6 ± 5 ng·mL⁻¹), as compared to healthy controls (1.4 ± 0.4 ng·mL⁻¹) of the same average age. We identified unique mutations in the cfDNA and tumour DNA of each GBM patient, including some of the most frequently mutated genes in GBM according to the COSMIC database (TP53, 18.75%; EGFR, 37.5%; NF1, 12.5%; LRP1B, 25%; IRS4, 25%). Using our gene–gene fusion database, ChiTaRS 5.0, we identified gene–gene fusions in cfDNA and tumour DNA, such as KDR–PDGFRA and NCDN–PDGFRA, which correspond to previously reported alterations of PDGFRA in GBM (44% of all samples). Interestingly, the PDGFRA protein fusions can be targeted by tyrosine kinase inhibitors such as imatinib, sunitinib, and sorafenib. Moreover, we identified BCR–ABL1 (in 8% of patients), COL1A1–PDGFB (8%), NIN–PDGFRB (8%), and FGFR1–BCR (4%) in cfDNA of patients, which can be targeted by analogues of imatinib. ROS1 fusions (CEP85L–ROS1 and GOPC–ROS1), identified in 8% of patient cfDNA, might be targeted by crizotinib, entrectinib, or larotrectinib. Thus, our study suggests that integrated analysis of cfDNA plasma concentration, gene mutations, and gene–gene fusions can serve as a diagnostic modality for distinguishing GBM patients who may benefit from targeted therapy. These results open new avenues for precision medicine in GBM, using noninvasive liquid biopsy diagnostics to assess personalized patient profiles. Moreover, repeated detection of druggable targets over the course of the disease may provide real‐time information on the evolving molecular landscape of the tumour.

Fusion genes in gynecologic tumors: the occurrence, molecular mechanism and prospect for therapy

Gene fusions are thought to be driver mutations in multiple cancers and are an important factor for poor patient prognosis. Most of them appear in specific cancers, thus satisfactory strategies can be developed for the precise treatment of these types of cancer. Currently, there are few targeted drugs to treat gynecologic tumors, and patients with gynecologic cancer often have a poor prognosis because of tumor progression or recurrence. With the application of massively parallel sequencing, a large number of fusion genes have been discovered in gynecologic tumors, and some fusions have been confirmed to be involved in the biological process of tumor progression. To this end, the present article reviews the current research status of all confirmed fusion genes in gynecologic tumors, including their rearrangement mechanism and frequency in ovarian cancer, endometrial cancer, endometrial stromal sarcoma, and other types of uterine tumors. We also describe the mechanisms by which fusion genes are generated and their oncogenic mechanism. Finally, we discuss the prospect of fusion genes as therapeutic targets in gynecologic tumors.

Identification of KANSARL as the first cancer predisposition fusion gene specific to the population of European ancestry origin

Gene fusion is one of the hallmarks of cancer. Recent advances in RNA-seq of cancer transcriptomes have facilitated the discovery of fusion transcripts. In this study, we report identification of a surprisingly large number of fusion transcripts, including six KANSARL (KANSL1-ARL17A) transcripts that resulted from the fusion between the KANSL1 and ARL17A genes using a RNA splicingcode model. Five of these six KANSARL fusion transcripts are novel. By systematic analysis of RNA-seq data of glioblastoma, prostate cancer, lung cancer, breast cancer, and lymphoma from different regions of the World, we have found that KANSARL fusion transcripts were rarely detected in the tumors of individuals from Asia or Africa. In contrast, they exist in 30 - 52% of the tumors from North Americans cancer patients. Analysis of CEPH/Utah Pedigree 1463 has revealed that KANSARL is a familially-inherited fusion gene. Further analysis of RNA-seq datasets of the 1000 Genome Project has indicated that KANSARL fusion gene is specific to 28.9% of the population of European ancestry origin. In summary, we demonstrated that KANSARL is the first cancer predisposition fusion gene associated with genetic backgrounds of European ancestry origin.

Regulator of G Protein Signaling 17 as a Negative Modulator of GPCR Signaling in Multiple Human Cancers

Regulators of G protein signaling (RGS) proteins modulate G protein-coupled receptor (GPCR) signaling networks by terminating signals produced by active Gα subunits. RGS17, a member of the RZ subfamily of RGS proteins, is typically only expressed in appreciable amounts in the human central nervous system, but previous works have shown that RGS17 expression is selectively upregulated in a number of malignancies, including lung, breast, prostate, and hepatocellular carcinoma. In addition, this upregulation of RGS17 is associated with a more aggressive cancer phenotype, as increased proliferation, migration, and invasion are observed. Conversely, decreased RGS17 expression diminishes the response of ovarian cancer cells to agents commonly used during chemotherapy. These somewhat contradictory roles of RGS17 in cancer highlight the need for selective, high-affinity inhibitors of RGS17 to use as chemical probes to further the understanding of RGS17 biology. Based on current evidence, these compounds could potentially have clinical utility as novel chemotherapeutics in the treatment of lung, prostate, breast, and liver cancers. Recent advances in screening technologies to identify potential inhibitors coupled with increasing knowledge of the structural requirements of RGS-Gα protein-protein interaction inhibitors make the future of drug discovery efforts targeting RGS17 promising. This review highlights recent findings related to RGS17 as both a canonical and atypical RGS protein, its role in various human disease states, and offers insights on small molecule inhibition of RGS17.

Fusion of the TBL1XR1 and HMGA1 genes in splenic hemangioma with t(3;6)(q26;p21)

RNA-sequencing of a splenic hemangioma with the karyotype 45~47,XX,t(3;6)(q26;p21) showed that this translocation generated a chimeric TBL1XR1-HMGA1 gene. This is the first time that this tumor has been subjected to genetic analysis, but the finding of an acquired clonal chromosome abnormality in cells cultured from the lesion and the presence of the TBL1XR1-HMGA1 fusion in them strongly favor the conclusion that splenic hemangiomas are of a neoplastic nature. Genomic PCR confirmed the presence of the TBL1XR1-HMGA1 fusion gene, and RT-PCR together with Sanger sequencing verified the presence of the fusion transcripts. The molecular consequences of the t(3;6) would be substantial. The cells carrying the translocation would retain only one functional copy of the wild-type TBL1XR1 gene while the other, rearranged allele could produce a putative truncated form of TBL1XR1 protein containing the LiSH and F-box-like domains. In the TBL1XR1-HMGA1 fusion transcript, furthermore, untranslated exons of HMGA1 are replaced by the first 5 exons of the TBL1XR1 gene. The result is that the entire coding region of HMGA1 comes under the control of the TBL1XR1 promoter, bringing about dysregulation of HMGA1. This is reminiscent of similar pathogenetic mechanisms involving high mobility genes in benign connective tissue tumors such as lipomas and leiomyomas.

Evolutionary significance and diversification of the phosphoglucose isomerase genes in vertebrates

Background

Phosphoglucose isomerase (PGI) genes are important multifunctional proteins whose evolution has, until now, not been well elucidated because of the limited number of completely sequenced genomes. Although the multifunctionality of this gene family has been considered as an original and innate characteristic, PGI genes may have acquired novel functions through changes in coding sequences and exon/intron structure, which are known to lead to functional divergence after gene duplication. A whole-genome comparative approach was used to estimate the rates of molecular evolution of this protein family.

Results

The results confirm the presence of two isoforms in teleost fishes and only one variant in all other vertebrates. Phylogenetic reconstructions grouped the PGI genes into five main groups: lungfishes/coelacanth/cartilaginous fishes, teleost fishes, amphibians, reptiles/birds and mammals, with the teleost group being subdivided into two subclades comprising PGI1 and PGI2. This PGI partitioning into groups is consistent with the synteny and molecular evolution results based on the estimation of the ratios of nonsynonymous to synonymous changes (Ka/Ks) and divergence rates between both PGI paralogs and orthologs. Teleost PGI2 shares more similarity with the variant found in all other vertebrates, suggesting that it has less evolved than PGI1 relative to the PGI of common vertebrate ancestor.

Conclusions

The diversification of PGI genes into PGI1 and PGI2 is consistent with a teleost-specific duplication before the radiation of this lineage, and after its split from the other infraclasses of ray-finned fishes. The low average Ka/Ks ratios within teleost and mammalian lineages suggest that both PGI1 and PGI2 are functionally constrained by purifying selection and may, therefore, have the same functions. By contrast, the high average Ka/Ks ratios and divergence rates within reptiles and birds indicate that PGI may be involved in different functions. The synteny analyses show that the genomic region harbouring PGI genes has independently undergone genomic rearrangements in mammals versus the reptile/bird lineage in particular, which may have contributed to the actual functional diversification of this gene family.

Electronic supplementary material

The online version of this article (doi:10.1186/s13104-015-1683-x) contains supplementary material, which is available to authorized users.

Recurrent ESR1-CCDC170 rearrangements in an aggressive subset of oestrogen receptor-positive breast cancers

Characterizing the genetic alterations leading to the more aggressive forms of estrogen receptor positive (ER+) breast cancers are of critical significance in breast cancer management. Here we identify recurrent rearrangements between estrogen receptor gene ESR1 and its neighbor CCDC170, which are enriched in the more aggressive and endocrine-resistant luminal-B tumors, through large-scale analyses of breast cancer transcriptome and copy number alterations. Further screening of 200 ER+ breast cancers identifies eight ESR1-CCDC170 positive tumors. These fusions encode N-terminally truncated CCDC170 proteins (ΔCCDC170). When introduced into ER+ breast cancer cells, ΔCCDC170 leads to markedly increased cell motility and anchorage-independent growth, reduced endocrine sensitivity, and enhanced xenograft tumor formation. Mechanistic studies suggest that ΔCCDC170 engages Gab1 signalosome to potentiate growth factor signaling and enhance cell motility. Together, this study identifies neoplastic ESR1-CCDC170 fusions in a more aggressive subset of ER+ breast cancer, which suggests a new concept of ER pathobiology in breast cancer.

Construction of mate pair full-length cDNAs libraries and characterization of transcriptional start sites and termination sites

To identify and characterize transcript structures ranging from transcriptional start sites (TSSs) to poly(A)-addition sites (PASs), we constructed and analyzed human TSS/PAS mate pair full-length cDNA libraries from 14 tissue types and four cell lines. The collected information enabled us to define TSS cluster (TSC) and PAS cluster (PAC) relationships for a total of 8530/9400 RefSeq genes, as well as 4251/5618 of their putative alternative promoters/terminators and 4619/4605 intervening transcripts, respectively. Analyses of the putative alternative TSCs and alternative PACs revealed that their selection appeared to be mostly independent, with rare exceptions. In those exceptional cases, pairs of transcript units rarely overlapped one another and were occasionally separated by Rad21/CTCF. We also identified a total of 172 similar cases in which TSCs and PACs spanned adjacent but distinct genes. In these cases, different transcripts may utilize different functional units of a particular gene or of adjacent genes. This approach was also useful for identifying fusion gene transcripts in cancerous cells. Furthermore, we could construct cDNA libraries in which 3′-end mate pairs were distributed randomly over the transcripts. These libraries were useful for assembling the internal structure of previously uncharacterized alternative promoter products, as well as intervening transcripts.

New methods as alternative or corrective measures for the pitfalls and artifacts of reverse transcription and polymerase chain reactions (RT-PCR) in cloning chimeric or antisense-accompanied RNA

We established new methods for cloning cDNA ends that start with reverse transcription (RT) and soon proceed with the synthesis of the second cDNA strand, avoiding manipulations of fragile RNA. Our 3′-end cloning method does not involve poly-dT primers and polymerase chain reactions (PCR), is low in efficiency but high in fidelity and can clone those RNAs without a poly-A tail. We also established a cDNA protection assay to supersede RNA protection assay. The protected cDNA can be amplified, cloned and sequenced, enhancing sensitivity and fidelity. We report that RT product using gene-specific primer (GSP) cannot be gene- or strand-specific because RNA sample contains endogenous random primers (ERP). The gene-specificity may be improved by adding a linker sequence at the 5′-end of the GSP to prime RT and using the linker as a primer in the ensuing PCR. The strand-specificity may be improved by using strand-specific DNA oligos in our protection assay. The CDK4 mRNA and TSPAN31 mRNA are transcribed from the opposite DNA strands and overlap at their 3′ ends. Using this relationship as a model, we found that the overlapped sequence might serve as a primer with its antisense as the template to create a wrong-template extension in RT or PCR. We infer that two unrelated RNAs or cDNAs overlapping at the 5′- or 3′-end might create a spurious chimera in this way, and many chimeras with a homologous sequence may be such artifacts. The ERP and overlapping antisense together set complex pitfalls, which one should be aware of.

Novel domain combinations in proteins encoded by chimeric transcripts

Motivation: Chimeric RNA transcripts are generated by different mechanisms including pre-mRNA trans-splicing, chromosomal translocations and/or gene fusions. It was shown recently that at least some of chimeric transcripts can be translated into functional chimeric proteins.

Results: To gain a better understanding of the design principles underlying chimeric proteins, we have analyzed 7,424 chimeric RNAs from humans. We focused on the specific domains present in these proteins, comparing their permutations with those of known human proteins. Our method uses genomic alignments of the chimeras, identification of the gene–gene junction sites and prediction of the protein domains. We found that chimeras contain complete protein domains significantly more often than in random data sets. Specifically, we show that eight different types of domains are over-represented among all chimeras as well as in those chimeras confirmed by RNA-seq experiments. Moreover, we discovered that some chimeras potentially encode proteins with novel and unique domain combinations. Given the observed prevalence of entire protein domains in chimeras, we predict that certain putative chimeras that lack activation domains may actively compete with their parental proteins, thereby exerting dominant negative effects. More generally, the production of chimeric transcripts enables a combinatorial increase in the number of protein products available, which may disturb the function of parental genes and influence their protein–protein interaction network.

Availability: our scripts are available upon request.

Contact:avalencia@cnio.es

Supplementary information:Supplementary data are available at Bioinformatics online.

The potential of carboxypeptidase M as a therapeutic target in cancer

INTRODUCTION

In the recent literature, carboxypeptidase M (CPM) emerged as a potential cancer biomarker. CPM modulates receptor signaling of kinins, anaphylatoxins, and chemokines. These CPM substrates affect proliferation, angiogenesis, and apoptosis of cancer cells. What is the evidence that CPM is a drug target for cancer therapy?

AREAS COVERED

The literature was searched using PubMed with the search terms "carboxypeptidase M" and/or "chromosome 12q13-15" eventually combined with general terms related to cancer. Information was retrieved from the GEO database and material of gene expression and proteomic studies.

EXPERT OPINION

CPM is a part of the molecular signature of many cancers. There is good evidence that it is useful for the discrimination and stratification of cancer types, possibly in combination with other markers such as EGFR and MDM2. Whether it is also a drug target remains to be determined. Lung, kidney, brain, and the reproductive system contain relatively high levels of CPM, but its functions in those tissues are largely unknown. CPM is expressed on tumor-associated macrophages. To facilitate the investigation of CPM in tumor-associated inflammation and in the other aspects of tumor biology, it is necessary to develop potent and selective CPM inhibitors.

Detection of a common chimeric transcript between human chromosomes 7 and 16

Abstract

Interchromosomal chimeric RNA molecules are often transcription products from genomic rearrangement in cancerous cells. Here we report the computational detection of an interchromosomal RNA fusion between ZC3HAV1L and CHMP1A from RNA-seq data of normal human mammary epithelial cells, and experimental confirmation of the chimeric transcript in multiple human cells and tissues. Our experimental characterization also detected three variants of the ZC3HAV1L-CHMP1A chimeric RNA, suggesting that these genes are involved in complex splicing. The fusion sequence at the novel exon-exon boundary, and the absence of corresponding DNA rearrangement suggest that this chimeric RNA is likely produced by trans-splicing in human cells.

Reviewers

This article was reviewed by Rory Johnson (nominated by Fyodor Kondrashov); Gal Avital and Itai Yanai

Using chimaeric expression sequence tag as the reference to identify three-dimensional chromosome contacts

Transcription-induced chimaeric transcripts, the potential post-transcriptional processing products, might reflect the spatial proximity of actively transcribed genes co-localized in transcription factories. A growing number of expression data deposited in databases provide us with the raw material for screening such chimaeric transcripts and using them as the probes to identify interactions between genes in cis or in trans. Based on the high-quality chimaeric transcripts gleaned from human expression sequence tag data with selection criteria, we identified the patterns of inter- and intrachromosomal gene–gene interactions. On top the contact pattern from interchromosomal interactions, we also observed an exponential behaviour of the intrachromosomal interactions within a certain length scale, which is consistent with the independent experimental results from Hi-C screening and with the Random Loop Model. A compatible result is found for mouse. Transcription-induced chimaeric transcripts, most of which might be accidental products with trivial functions, shed light on the spatial organization of chromosomes. These inter- and intrachromosomal interactions might contribute to the compaction of chromosomes, their segregation and formation of the chromosome territories, and their spatial distribution within the nucleus.

Carboxypeptidase M in apoptosis, adipogenesis and cancer

This review covers carboxypeptidase M (CPM) research that appeared in the literature since 2009. The focus is on aspects that are new or interesting from a clinical perspective. Available research tools are discussed as well as their pitfalls and limitations. Evidence is provided to suggest the potential involvement of CPM in apoptosis, adipogenesis and cancer. This evidence derives from the expression pattern of CPM and its putative substrates in cells and tissues. In recent years CPM emerged as a potential cancer biomarker, in well differentiated liposarcoma where the CPM gene is co-amplified with the oncogene MDM2; and in lung adenocarcinoma where coexpression with EGFR correlates with poor prognosis. The available data call for extended investigation of the function of CPM in tumor cells, tumor-associated macrophages, stromal cells and tumor neovascularisation. Such experiments could be instrumental to validate CPM as a therapeutic target.

ChiTaRS: a database of human, mouse and fruit fly chimeric transcripts and RNA-sequencing data

Chimeric RNAs that comprise two or more different transcripts have been identified in many cancers and among the Expressed Sequence Tags (ESTs) isolated from different organisms; they might represent functional proteins and produce different disease phenotypes. The ChiTaRS database of Chimeric Transcripts and RNA-Sequencing data (http://chitars.bioinfo.cnio.es/) collects more than 16 000 chimeric RNAs from humans, mice and fruit flies, 233 chimeras confirmed by RNA-seq reads and ∼2000 cancer breakpoints. The database indicates the expression and tissue specificity of these chimeras, as confirmed by RNA-seq data, and it includes mass spectrometry results for some human entries at their junctions. Moreover, the database has advanced features to analyze junction consistency and to rank chimeras based on the evidence of repeated junction sites. Finally, ‘Junction Search’ screens through the RNA-seq reads found at the chimeras’ junction sites to identify putative junctions in novel sequences entered by users. Thus, ChiTaRS is an extensive catalog of human, mouse and fruit fly chimeras that will extend our understanding of the evolution of chimeric transcripts in eukaryotes and can be advantageous in the analysis of human cancer breakpoints.

Identification and analysis of pig chimeric mRNAs using RNA sequencing data

Background

Gene fusion is ubiquitous over the course of evolution. It is expected to increase the diversity and complexity of transcriptomes and proteomes through chimeric sequence segments or altered regulation. However, chimeric mRNAs in pigs remain unclear. Here we identified some chimeric mRNAs in pigs and analyzed the expression of them across individuals and breeds using RNA-sequencing data.

Results

The present study identified 669 putative chimeric mRNAs in pigs, of which 251 chimeric candidates were detected in a set of RNA-sequencing data. The 618 candidates had clear trans-splicing sites, 537 of which obeyed the canonical GU-AG splice rule. Only two putative pig chimera variants whose fusion junction was overlapped with that of a known human chimeric mRNA were found. A set of unique chimeric events were considered middle variances in the expression across individuals and breeds, and revealed non-significant variance between sexes. Furthermore, the genomic region of the 5^′ partner gene shares a similar DNA sequence with that of the 3^′ partner gene for 458 putative chimeric mRNAs. The 81 of those shared DNA sequences significantly matched the known DNA-binding motifs in the JASPAR CORE database. Four DNA motifs shared in parental genomic regions had significant similarity with known human CTCF binding sites.

Conclusions

The present study provided detailed information on some pig chimeric mRNAs. We proposed a model that trans-acting factors, such as CTCF, induced the spatial organisation of parental genes to the same transcriptional factory so that parental genes were coordinatively transcribed to give birth to chimeric mRNAs.

Chromosomal aberrations in solid tumors

Ever since the identification of the exact number of human chromosomes in 1956, several cancer-specific chromosomal abnormalities have been identified in different tumors. Among the various genetic changes, such as alterations in oncogenes, tumor suppressor genes, and microRNA genes, recurrent chromosome translocations have been identified as an important class of mutations in hematological malignancies, soft tissue sarcomas, and more recently in prostate cancer and lung cancer. Recurrent gene fusions are used for cancer classification and as diagnostic markers, and some have been successfully targeted for drug development. Recent advances in high-throughput sequencing technology and the ambitious undertaking of "The Cancer Genome Atlas" (TCGA) project will help drive the identification of the underlying genetic aberrations in most of the solid cancers. This chapter presents an overview on the current status of the knowledge on chromosome aberrations in solid cancers, cytogenetic and noncytogenetic methods for the characterization of changes at the DNA and RNA levels, technological advancements in high-throughput characterization of the cancer genome and transcriptome, and the current understanding of the molecular mechanism involved in the formation of gene fusions in solid cancer.

Sensitive gene fusion detection using ambiguously mapping RNA-Seq read pairs

MOTIVATION

Paired-end whole transcriptome sequencing provides evidence for fusion transcripts. However, due to the repetitiveness of the transcriptome, many reads have multiple high-quality mappings. Previous methods to find gene fusions either ignored these reads or required additional longer single reads. This can obscure up to 30% of fusions and unnecessarily discards much of the data.

RESULTS

We present a method for using paired-end reads to find fusion transcripts without requiring unique mappings or additional single read sequencing. Using simulated data and data from tumors and cell lines, we show that our method can find fusions with ambiguously mapping read pairs without generating numerous spurious fusions from the many mapping locations.

AVAILABILITY

A C++ and Python implementation of the method demonstrated in this article is available at http://exon.ucsd.edu/ShortFuse.

CONTACT

mckinsel@ucsd.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Transcriptional consequences of genomic structural aberrations in breast cancer

Using a long-span, paired-end deep sequencing strategy, we have comprehensively identified cancer genome rearrangements in eight breast cancer genomes. Herein, we show that 40%-54% of these structural genomic rearrangements result in different forms of fusion transcripts and that 44% are potentially translated. We find that single segmental tandem duplication spanning several genes is a major source of the fusion gene transcripts in both cell lines and primary tumors involving adjacent genes placed in the reverse-order position by the duplication event. Certain other structural mutations, however, tend to attenuate gene expression. From these candidate gene fusions, we have found a fusion transcript (RPS6KB1-VMP1) recurrently expressed in ∼30% of breast cancers associated with potential clinical consequences. This gene fusion is caused by tandem duplication on 17q23 and appears to be an indicator of local genomic instability altering the expression of oncogenic components such as MIR21 and RPS6KB1.

Fusion genes and chromosome translocations in the common epithelial cancers

It has been known for 25 years that fusion genes play a central role in leukaemias and sarcomas but they have been neglected in the common carcinomas, largely because of technical limitations of cytogenetics. In the last few years it has emerged that gene fusions, caused by chromosome translocations, inversions, deletions, etc., are important in the common epithelial cancers, such as prostate and lung carcinoma. Most prostate cancers, for example, have an androgen-regulated fusion of one of the ETS transcription factor gene family. Early results of genome-wide searches for gene fusions in breast and other epithelial cancers suggest that most individual tumours will have several fused genes. Fusion genes are exceptionally powerful mutations. In their simplest form they can turn on expression by promoter insertion but they can also, for example, force dimerization of a protein or change its subcellular location. They are correspondingly important clinically, in classification and management and as targets for therapy. This review surveys what we know of fusion genes in the carcinomas, summarizes the technical advances that now make it possible to search systematically for such genes, and concludes by putting fusion genes into the current picture of mutation in cancers.

Detection of human interchromosomal trans-splicing in sequence databanks

Trans-splicing is a common phenomenon in nematodes and kinetoplastids, and it has also been reported in other organisms, including humans. Up to now, all in silico strategies to find evidence of trans-splicing in humans have required that the candidate sequences follow the consensus splicing site rules (spliceosome-mediated mechanism). However, this criterion is not supported by the best human experimental evidence, which, except in a single case, do not follow canonical splicing sites. Moreover, recent findings describe a novel alternative tRNA mediated trans-splicing mechanism, which prescinds the spliceosome machinery. In order to answer the question, 'Are there hybrid mRNAs in sequence databanks, whose characteristics resemble those of the best human experimental evidence?', we have developed a methodology that successfully identified 16 hybrid mRNAs which might be instances of interchromosomal trans-splicing. Each hybrid mRNA is formed by a trans-spliced region (TSR), which was successfully mapped either onto known genes or onto a human endogenous retrovirus (HERV-K) transcript which supports their transcription. The existence of these hybrid mRNAs indicates that trans-splicing may be more widespread than believed. Furthermore, non-canonical splice site patterns suggest that infrequent splicing sites may occur under special conditions, or that an alternative trans-splicing mechanism is involved. Finally, our candidates are supposedly from normal tissue, and a recent study has reported that trans-splicing may occur not only in malignant tissues, but in normal tissues as well. Our methodology can be applied to 5'-UTR, coding sequences and 3'-UTR in order to find new candidates for a posteriori experimental confirmation.

ChimerDB 2.0--a knowledgebase for fusion genes updated

Chromosome translocations and gene fusions are frequent events in the human genome and have been found to cause diverse types of tumor. ChimerDB is a knowledgebase of fusion genes identified from bioinformatics analysis of transcript sequences in the GenBank and various other public resources such as the Sanger cancer genome project (CGP), OMIM, PubMed and the Mitelman’s database. In this updated version, we significantly modified the algorithm of identifying fusion transcripts. Specifically, the new algorithm is more sensitive and has detected 2699 fusion transcripts with high confidence. Furthermore, it can identify interchromosomal translocations as well as the intrachromosomal deletions or inversions of large DNA segments. Importantly, results from the analysis of next-generation sequencing data in the short read archives are incorporated as well. We updated and integrated all contents (GenBank, Sanger CGP, OMIM, PubMed publications and the Mitelman’s database), and the user-interface has been improved to support diverse types of searches and to enhance the user convenience especially in browsing PubMed articles. We also developed a new alignment viewer that should facilitate examining reliability of fusion transcripts and inferring functional significance. We expect ChimerDB 2.0, available at http://ercsb.ewha.ac.kr/fusiongene, to be a valuable tool in identifying biomarkers and drug targets.

Chimeric transcript discovery by paired-end transcriptome sequencing

Recurrent gene fusions are a prevalent class of mutations arising from the juxtaposition of 2 distinct regions, which can generate novel functional transcripts that could serve as valuable therapeutic targets in cancer. Therefore, we aim to establish a sensitive, high-throughput methodology to comprehensively catalog functional gene fusions in cancer by evaluating a paired-end transcriptome sequencing strategy. Not only did a paired-end approach provide a greater dynamic range in comparison with single read based approaches, but it clearly distinguished the high-level "driving" gene fusions, such as BCR-ABL1 and TMPRSS2-ERG, from potential lower level "passenger" gene fusions. Also, the comprehensiveness of a paired-end approach enabled the discovery of 12 previously undescribed gene fusions in 4 commonly used cell lines that eluded previous approaches. Using the paired-end transcriptome sequencing approach, we observed read-through mRNA chimeras, tissue-type restricted chimeras, converging transcripts, diverging transcripts, and overlapping mRNA transcripts. Last, we successfully used paired-end transcriptome sequencing to detect previously undescribed ETS gene fusions in prostate tumors. Together, this study establishes a highly specific and sensitive approach for accurately and comprehensively cataloguing chimeras within a sample using paired-end transcriptome sequencing.

Mutation-associated fusion cancer genes in solid tumors

Chromosomal translocations and fusion oncogenes serve as the ultimate biomarker for clinicians as they show specificity for distinct histopathologic malignancies while simultaneously encoding an etiologic mutation and a therapeutic target. Previously considered a minor mutational event in epithelial solid tumors, new methodologies that do not rely on the detection of macroscopic cytogenetic alterations, as well as access to large series of annotated clinical material, are expanding the inventory of recurrent fusion oncogenes in both common and rare solid epithelial tumors. Unexpectedly, related assays are also revealing a high number of tandem or chimeric transcripts in normal tissues including, in one provocative case, a template for a known fusion oncogene. These observations may force us to reassess long-held views on the definition of a gene. They also raise the possibility that some rearrangements might represent constitutive forms of a physiological chimeric transcript. Defining the chimeric transcriptome in both health (transcription-induced chimerism and intergenic splicing) and disease (mutation-associated fusion oncogenes) will play an increasingly important role in the diagnosis, prognosis, and therapy of patients with cancer.

Oncogenic gene fusions in epithelial carcinomas

New discoveries regarding recurrent chromosomal aberrations in epithelial tumors have challenged the view that gene fusions play a minor role in these cancers. It is now known that recurrent fusions characterize significant subsets of prostate, breast, lung and renal-cell carcinomas, among others. This work has generated new insights into the molecular subtypes of tumors and highlighted important advances in bioinformatics, sequencing, and microarray technology as tools for gene fusion discovery. Given the ubiquity of tyrosine kinases and transcription factors in gene fusions, further interest in the potential 'druggability' of gene fusions with targeted therapeutics has also flourished. Nevertheless, the majority of chromosomal abnormalities in epithelial cancers remain uncharacterized, underscoring the limitations of our knowledge of carcinogenesis and the requirement for further research.

A sequence-level map of chromosomal breakpoints in the MCF-7 breast cancer cell line yields insights into the evolution of a cancer genome

By applying a method that combines end-sequence profiling and massively parallel sequencing, we obtained a sequence-level map of chromosomal aberrations in the genome of the MCF-7 breast cancer cell line. A total of 157 distinct somatic breakpoints of two distinct types, dispersed and clustered, were identified. A total of 89 breakpoints are evenly dispersed across the genome. A majority of dispersed breakpoints are in regions of low copy repeats (LCRs), indicating a possible role for LCRs in chromosome breakage. The remaining 68 breakpoints form four distinct clusters of closely spaced breakpoints that coincide with the four highly amplified regions in MCF-7 detected by array CGH located in the 1p13.1-p21.1, 3p14.1-p14.2, 17q22-q24.3, and 20q12-q13.33 chromosomal cytobands. The clustered breakpoints are not significantly associated with LCRs. Sequences flanking most (95%) breakpoint junctions are consistent with double-stranded DNA break repair by nonhomologous end-joining or template switching. A total of 79 known or predicted genes are involved in rearrangement events, including 10 fusions of coding exons from different genes and 77 other rearrangements. Four fusions result in novel expressed chimeric mRNA transcripts. One of the four expressed fusion products (RAD51C-ATXN7) and one gene truncation (BRIP1 or BACH1) involve genes coding for members of protein complexes responsible for homology-driven repair of double-stranded DNA breaks. Another one of the four expressed fusion products (ARFGEF2-SULF2) involves SULF2, a regulator of cell growth and angiogenesis. We show that knock-down of SULF2 in cell lines causes tumorigenic phenotypes, including increased proliferation, enhanced survival, and increased anchorage-independent growth.

Array painting reveals a high frequency of balanced translocations in breast cancer cell lines that break in cancer-relevant genes

Chromosome translocations in the common epithelial cancers are abundant, yet little is known about them. They have been thought to be almost all unbalanced and therefore dismissed as mostly mediating tumour suppressor loss. We present a comprehensive analysis by array painting of the chromosome translocations of breast cancer cell lines HCC1806, HCC1187 and ZR-75-30. In array painting, chromosomes are isolated by flow cytometry, amplified and hybridized to DNA microarrays. A total of 200 breakpoints were identified and all were mapped to 1 Mb resolution on bacterial artificial chromosome (BAC) arrays, then 40 selected breakpoints, including all balanced breakpoints, were further mapped on tiling-path BAC arrays or to around 2 kb resolution using oligonucleotide arrays. Many more of the translocations were balanced at 1 Mb resolution than expected, either reciprocal (eight in total) or balanced for at least one participating chromosome (19 paired breakpoints). Second, many of the breakpoints were at genes that are plausible targets of oncogenic translocation, including balanced breaks at CTCF, EP300/p300 and FOXP4. Two gene fusions were demonstrated, TAX1BP1-AHCY and RIF1-PKD1L1. Our results support the idea that chromosome rearrangements may play an important role in common epithelial cancers such as breast cancer.

Fusion transcripts and transcribed retrotransposed loci discovered through comprehensive transcriptome analysis using Paired-End diTags (PETs)

Identification of unconventional functional features such as fusion transcripts is a challenging task in the effort to annotate all functional DNA elements in the human genome. Paired-End diTag (PET) analysis possesses a unique capability to accurately and efficiently characterize the two ends of DNA fragments, which may have either normal or unusual compositions. This unique nature of PET analysis makes it an ideal tool for uncovering unconventional features residing in the human genome. Using the PET approach for comprehensive transcriptome analysis, we were able to identify fusion transcripts derived from genome rearrangements and actively expressed retrotransposed pseudogenes, which would be difficult to capture by other means. Here, we demonstrate this unique capability through the analysis of 865,000 individual transcripts in two types of cancer cells. In addition to the characterization of a large number of differentially expressed alternative 5' and 3' transcript variants and novel transcriptional units, we identified 70 fusion transcript candidates in this study. One was validated as the product of a fusion gene between BCAS4 and BCAS3 resulting from an amplification followed by a translocation event between the two loci, chr20q13 and chr17q23. Through an examination of PETs that mapped to multiple genomic locations, we identified 4055 retrotransposed loci in the human genome, of which at least three were found to be transcriptionally active. The PET mapping strategy presented here promises to be a useful tool in annotating the human genome, especially aberrations in human cancer genomes.

HYBRIDdb: a database of hybrid genes in the human genome

Background

Hybrid genes are candidate risk factors for human tumors by inducing mutation, translocation, inversion, or rearrangement of genes. The occurrence of hybrid genes may also have given rise to new transcripts during hominid evolution.

Description

HYBRIDdb is a database of hybrid genes in humans. This system encompasses the bioinformatics analysis of mRNA, EST, cDNA, and genomic DNA sequences in the INDC databases, and can be used to identify hybrid genes. We searched for hybrid genes among the 28,171 genes listed in the NCBI database, and analyzed their structural patterns in the human genome. The 2,344 gene pairs were detected as hybrid forms of transcriptional products. We classified the hybrid genes into two groups: chromosomal-mediated translocation fusion transcripts and transcription-mediated fusion transcripts.

Conclusion

The HYBRIDdb database will provide genome scientists with insight into potential roles for hybrid genes in human evolution and disease.

Tentative mapping of transcription-induced interchromosomal interaction using chimeric EST and mRNA data

Recent studies on chromosome conformation show that chromosomes colocalize in the nucleus, bringing together active genes in transcription factories. This spatial proximity of actively transcribing genes could provide a means for RNA interaction at the transcript level. We have screened public databases for chimeric EST and mRNA sequences with the intent of mapping transcription-induced interchromosomal interactions. We suggest that chimeric transcripts may be the result of close encounters of active genes, either as functional products or "noise" in the transcription process, and that they could be used as probes for chromosome interactions. We have found a total of 5,614 chimeric ESTs and 587 chimeric mRNAs that meet our selection criteria. Due to their higher quality, the mRNA findings are of particular interest and we hope that they may serve as food for thought for specialists in diverse areas of molecular biology.

Molecular classification of mucoepidermoid carcinomas-prognostic significance of the MECT1-MAML2 fusion oncogene

Mucoepidermoid carcinomas (MECs) of the salivary and bronchial glands are characterized by a recurrent t(11;19)(q21;p13) translocation resulting in a MECT1-MAML2 fusion in which the CREB-binding domain of the CREB coactivator MECT1 (also known as CRTC1, TORC1 or WAMTP1) is fused to the transactivation domain of the Notch coactivator MAML2. To gain further insights into the molecular pathogenesis of MECs, we cytogenetically and molecularly characterized a series of 29 MECs. A t(11;19) and/or an MECT1-MAML2 fusion was detected in more than 55% of the tumors. Several cases with cryptic rearrangements that resulted in gene fusions were detected. In fusion-negative MECs, the most common aberration was a single or multiple trisomies. Western blot and immunohistochemical studies demonstrated that the MECT1-MAML2 fusion protein was expressed in all MEC-specific cell types. In addition, cotransfection experiments showed that the fusion protein colocalized with CREB in homogeneously distributed nuclear granules. Analyses of potential downstream targets of the fusion revealed differential expression of the cAMP/CREB (FLT1 and NR4A2) and Notch (HES1 and HES5) target genes in fusion-positive and fusion-negative MECs. Moreover, clinical follow-up studies revealed that fusion-positive patients had a significantly lower risk of local recurrence, metastases, or tumor-related death compared to fusion-negative patients (P = 0.0012). When considering tumor-related deaths only, the estimated median survival for fusion-positive patients was greater than 10 years compared to 1.6 years for fusion-negative patients. These findings suggest that molecularly classifying MECs on the basis of an MECT1-MAML2 fusion is histopathologically and clinically relevant and that the fusion is a useful marker in predicting the biological behavior of MECs.

Decoding the fine-scale structure of a breast cancer genome and transcriptome

A comprehensive understanding of cancer is predicated upon knowledge of the structure of malignant genomes underlying its many variant forms and the molecular mechanisms giving rise to them. It is well established that solid tumor genomes accumulate a large number of genome rearrangements during tumorigenesis. End Sequence Profiling (ESP) maps and clones genome breakpoints associated with all types of genome rearrangements elucidating the structural organization of tumor genomes. Here we extend the ESP methodology in several directions using the breast cancer cell line MCF-7. First, targeted ESP is applied to multiple amplified loci, revealing a complex process of rearrangement and co-amplification in these regions reminiscent of breakage/fusion/bridge cycles. Second, genome breakpoints identified by ESP are confirmed using a combination of DNA sequencing and PCR. Third, in vitro functional studies assign biological function to a rearranged tumor BAC clone, demonstrating that it encodes anti-apoptotic activity. Finally, ESP is extended to the transcriptome identifying four novel fusion transcripts and providing evidence that expression of fusion genes may be common in tumors. These results demonstrate the distinct advantages of ESP including: (1) the ability to detect all types of rearrangements and copy number changes; (2) straightforward integration of ESP data with the annotated genome sequence; (3) immortalization of the genome; (4) ability to generate tumor-specific reagents for in vitro and in vivo functional studies. Given these properties, ESP could play an important role in a tumor genome project.

ChimerDB--a knowledgebase for fusion sequences

Chromosome translocation and gene fusion are frequent events in the human genome and are often the cause of many types of tumor. ChimerDB is the database of fusion sequences encompassing bioinformatics analysis of mRNA and expressed sequence tag (EST) sequences in the GenBank, manual collection of literature data and integration with other known database such as OMIM. Our bioinformatics analysis identifies the fusion transcripts that have non-overlapping alignments at multiple genomic loci. Fusion events at exon-exon borders are selected to filter out the cloning artifacts in cDNA library preparation. The result is classified into two groups--genuine chromosome translocation and fusion between neighboring genes owing to intergenic splicing. We also integrated manually collected literature and OMIM data for chromosome translocation as an aid to assess the validity of each fusion event. The database is available at http://genome.ewha.ac.kr/ChimerDB/ for human, mouse and rat genomes.

Alternative pre-mRNA processing regulates cell-type specific expression of the IL4l1 and NUP62 genes

BACKGROUND

Given the complexity of higher organisms, the number of genes encoded by their genomes is surprisingly small. Tissue specific regulation of expression and splicing are major factors enhancing the number of the encoded products. Commonly these mechanisms are intragenic and affect only one gene.

RESULTS

Here we provide evidence that the IL4I1 gene is specifically transcribed from the apparent promoter of the upstream NUP62 gene, and that the first two exons of NUP62 are also contained in the novel IL4I1_2 variant. While expression of IL4I1 driven from its previously described promoter is found mostly in B cells, the expression driven by the NUP62 promoter is restricted to cells in testis (Sertoli cells) and in the brain (e.g., Purkinje cells). Since NUP62 is itself ubiquitously expressed, the IL4I1_2 variant likely derives from cell type specific alternative pre-mRNA processing.

CONCLUSION

Comparative genomics suggest that the promoter upstream of the NUP62 gene originally belonged to the IL4I1 gene and was later acquired by NUP62 via insertion of a retroposon. Since both genes are apparently essential, the promoter had to serve two genes afterwards. Expression of the IL4I1 gene from the "NUP62" promoter and the tissue specific involvement of the pre-mRNA processing machinery to regulate expression of two unrelated proteins indicate a novel mechanism of gene regulation.

Allelic variation in gene expression identified through computational analysis of the dbEST database

Differential expression between the two alleles of an individual and between people with different genotypes has been commonly observed. Quantitative differences in gene expression between people may provide the genetic basis for the phenotypic difference between individuals and may be the primary cause of complex diseases. In this paper, we developed a computational method to identify genes that displayed allelic variation in gene expression in human EST libraries. To model allele-specific gene expression, we first identified EST libraries in which both A and B alleles were expressed and then identified allelic variation in gene expression based on the EST counts for each allele using a binomial test. Among 1107 SNPs that had a sufficient number of ESTs for the analysis, 524 (47%) displayed allelic variation in at least one cDNA library. We verified experimentally the allelic variation in gene expression for 6 of these SNPs. The frequency of allelic variation observed in EST libraries was similar to the previous studies using the SNP chip and primer extension method. We found that genes that displayed allelic variation were distributed throughout the human genome and were enriched in certain chromosome regions. The SNPs and genes identified in this study will provide a rich source for evaluating the effects of those SNPs and associated haplotypes in human health and diseases.

GMAP: a genomic mapping and alignment program for mRNA and EST sequences

MOTIVATION

We introduce GMAP, a standalone program for mapping and aligning cDNA sequences to a genome. The program maps and aligns a single sequence with minimal startup time and memory requirements, and provides fast batch processing of large sequence sets. The program generates accurate gene structures, even in the presence of substantial polymorphisms and sequence errors, without using probabilistic splice site models. Methodology underlying the program includes a minimal sampling strategy for genomic mapping, oligomer chaining for approximate alignment, sandwich DP for splice site detection, and microexon identification with statistical significance testing.

RESULTS

On a set of human messenger RNAs with random mutations at a 1 and 3% rate, GMAP identified all splice sites accurately in over 99.3% of the sequences, which was one-tenth the error rate of existing programs. On a large set of human expressed sequence tags, GMAP provided higher-quality alignments more often than blat did. On a set of Arabidopsis cDNAs, GMAP performed comparably with GeneSeqer. In these experiments, GMAP demonstrated a several-fold increase in speed over existing programs.

AVAILABILITY

Source code for gmap and associated programs is available at http://www.gene.com/share/gmap

SUPPLEMENTARY INFORMATION

http://www.gene.com/share/gmap.

Allelic Losses at 1p36 and 19q13 in Gliomas: Correlation with Histologic Classification, Definition of a 150-kb Minimal Deleted Region on 1p36, and Evaluation of CAMTA1 as a Candidate Tumor Suppressor Gene

Purpose: Allelic loss at 1p is seen in 70% to 85% of oligodendrogliomas (typically in association with 19q allelic loss) and 20-30% of astrocytomas. Because most 1p deletions in gliomas involve almost the entire chromosome arm, narrowing the region of the putative tumor suppressor gene has been difficult. To better define the histologic correlates of different patterns of 1p and 19q loss, we evaluated 1p/19q status in a large group of gliomas. This also allowed us to define a very small minimal deleted region (MDR) on 1p36.

Experimental Design: Among 205 consecutive cases of glioma studied for 1p loss of heterozygosity (LOH), 112 tumors were evaluated for both 1p and 19q LOH using at least three polymorphic markers on 1p and 19q each. The latter group included both low-grade tumors (oligodendroglioma, diffuse astrocytoma, and “oligoastrocytoma”) and high-grade tumors (anaplastic oligodendrogliomas, anaplastic astrocytomas, anaplastic oligoastrocytomas). Tumors with small segmental 1p losses (defined as LOH at some loci with retention of heterozygosity at other loci) were studied using a more extensive panel of markers to define the 1p MDR. The candidate gene was screened for mutations and its expression was studied by qualitative and quantitative reverse transcriptase-PCR and Northern blotting.

Results: Allelic losses on 1p and 19q, either separately or combined, were more common in classic oligodendrogliomas than in either astrocytomas or oligoastrocytomas (P < 0.0001). Classic oligodendrogliomas showed 1p loss in 35 of 42 (83%) cases, 19q loss in 28 of 39 (72%), and these were combined in 27 of 39 (69%) cases. There was no significant difference in 1p/19q LOH status between low-grade and anaplastic oligodendrogliomas. In contrast, no astrocytomas and only 6 of 30 (20%) oligoastrocytic tumors had combined 1p/19q loss. Although rare, 1p deletions were more often segmental in astrocytomas (5 of 6, 83%) than in oligodendrogliomas (3 of 35, 9%; P = 0.006). Eleven tumors (6 oligodendrogliomas or having oligodendroglial components, 5 purely astrocytic) with small segmental 1p losses underwent further detailed LOH mapping. All informative tumors in the oligodendroglial group and 2 of 3 informative astrocytomas showed LOH at 1p36.23, with a 150-kb MDR located between D1S2694 and D1S2666, entirely within the CAMTA1 transcription factor gene. Mutation analysis of the exons encoding conserved regions of CAMTA1 showed no somatic mutations in 10 gliomas, including 6 cases with and 4 cases without 1p LOH. CAMTA1 is normally expressed predominantly in non-neoplastic adult brain tissue. Relative to the latter, the expression level of CAMTA1 was low in oligodendroglial tumors and was further halved in cases with 1p deletion compared with those without 1p deletion (Mann-Whitney, P = 0.03).

Conclusions: Our data confirm the strong association of combined 1p/19q loss with classic oligodendroglioma histology and identify a very small segment of 1p36 located within CAMTA1 that was deleted in all oligodendroglial tumors with 1p LOH. This MDR also overlaps the neuroblastoma 1p36 MDR. CAMTA1 shows no evidence of inactivation by somatic mutations but its expression is reduced by half in cases with 1p LOH, suggesting that the functional effects of CAMTA1 haploinsufficiency warrant further investigation.