Bioinformatic analysis of exon repetition, exon scrambling and trans-splicing in humans

An RNA-Based Theory of Natural Universal Computation

Life is confronted with computation problems in a variety of domains including animal behavior, single-cell behavior, and embryonic development. Yet we currently do not know of a naturally existing biological system that is capable of universal computation, i.e., Turing-equivalent in scope. Generic finite-dimensional dynamical systems (which encompass most models of neural networks, intracellular signaling cascades, and gene regulatory networks) fall short of universal computation, but are assumed to be capable of explaining cognition and development. I present a class of models that bridge two concepts from distant fields: combinatory logic (or, equivalently, lambda calculus) and RNA molecular biology. A set of basic RNA editing rules can make it possible to compute any computable function with identical algorithmic complexity to that of Turing machines. The models do not assume extraordinarily complex molecular machinery or any processes that radically differ from what we already know to occur in cells. Distinct independent enzymes can mediate each of the rules and RNA molecules solve the problem of parenthesis matching through their secondary structure. In the most plausible of these models all of the editing rules can be implemented with merely cleavage and ligation operations at fixed positions relative to predefined motifs. This demonstrates that universal computation is well within the reach of molecular biology. It is therefore reasonable to assume that life has evolved - or possibly began with - a universal computer that yet remains to be discovered. The variety of seemingly unrelated computational problems across many scales can potentially be solved using the same RNA-based computation system. Experimental validation of this theory may immensely impact our understanding of memory, cognition, development, disease, evolution, and the early stages of life.

LncRNAs and PRC2: Coupled Partners in Embryonic Stem Cells

The power of embryonic stem cells (ESCs) lies in their ability to self-renew and differentiate. Behind these two unique capabilities is a fine-tuned molecular network that shapes the genetic, epigenetic, and epitranscriptomic ESC plasticity. Although RNA has been shown to be functionally important in only a small minority of long non-coding RNA genes, a growing body of evidence has highlighted the pivotal and intricate role of lncRNAs in chromatin remodeling. Due to their multifaceted nature, lncRNAs interact with DNA, RNA, and proteins, and are emerging as new modulators of extensive gene expression programs through their participation in ESC-specific regulatory circuitries. Here, we review the tight cooperation between lncRNAs and Polycomb repressive complex 2 (PRC2), which is intimately involved in determining and maintaining the ESC epigenetic landscape. The lncRNA-PRC2 partnership is fundamental in securing the fully pluripotent state of ESCs, which must be primed to differentiate properly. We also reflect on the advantages brought to this field of research by the advent of single-cell analysis.

NCLcomparator: systematically post-screening non-co-linear transcripts (circular, trans-spliced, or fusion RNAs) identified from various detectors

BACKGROUND

Non-co-linear (NCL) transcripts consist of exonic sequences that are topologically inconsistent with the reference genome in an intragenic fashion (circular or intragenic trans-spliced RNAs) or in an intergenic fashion (fusion or intergenic trans-spliced RNAs). On the basis of RNA-seq data, numerous NCL event detectors have been developed and detected thousands of NCL events in diverse species. However, there are great discrepancies in the identification results among detectors, indicating a considerable proportion of false positives in the detected NCL events. Although several helpful guidelines for evaluating the performance of NCL event detectors have been provided, a systematic guideline for measurement of NCL events identified by existing tools has not been available.

RESULTS

We develop a software, NCLcomparator, for systematically post-screening the intragenic or intergenic NCL events identified by various NCL detectors. NCLcomparator first examine whether the input NCL events are potentially false positives derived from ambiguous alignments (i.e., the NCL events have an alternative co-linear explanation or multiple matches against the reference genome). To evaluate the reliability of the identified NCL events, we define the NCL score (NCLscore) based on the variation in the number of supporting NCL junction reads identified by the tools examined. Of the input NCL events, we show that the ambiguous alignment-derived events have relatively lower NCLscore values than the other events, indicating that an NCL event with a higher NCLscore has a higher level of reliability. To help selecting highly expressed NCL events, NCLcomparator also provides a series of useful measurements such as the expression levels of the detected NCL events and their corresponding host genes and the junction usage of the co-linear splice junctions at both NCL donor and acceptor sites.

CONCLUSION

NCLcomparator provides useful guidelines, with the input of identified NCL events from various detectors and the corresponding paired-end RNA-seq data only, to help users selecting potentially high-confidence NCL events for further functional investigation. The software thus helps to facilitate future studies into NCL events, shedding light on the fundamental biology of this important but understudied class of transcripts. NCLcomparator is freely accessible at https://github.com/TreesLab/NCLcomparator .

EIciRNA-mediated gene expression: tunability and bimodality

Biological experiments have verified that EIciRNAs (a class of circRNA) produced from pre-mRNA can regulate gene expression, but the effect of regulation remains unexplored. Here, we refine a mechanistic gene model from experimental facts, in which we assume pre-mRNA synthesizes EIciRNAs and mRNAs in a probabilistic manner, with the probability called the pathway strength, and the resulting EIciRNAs positively regulate the pre-mRNA synthesis. We show that there is a critical pathway strength such that the mRNA mean and the mRNA noise reach the highest and lowest levels, respectively. The EIciRNA can induce the unimodal and bimodal mRNA expressions, as well as the transition between them. Our investigation hints that EIciRNA is a non-negligible factor affecting cell-to-cell variability in gene expression.

Motifome comparison between modern human, Neanderthal and Denisovan

BACKGROUND

The availability of the genomes of two archaic humans, Neanderthal and Denisovan, and that of modern humans provides researchers an opportunity to investigate genetic differences between these three subspecies on a genome-wide scale. Here we describe an algorithm that predicts statistically significant motifs based on the difference between a given motif's actual and expected distributions. The algorithm was previously applied to plants but was modified for this work.

RESULTS

The result of applying the algorithm to the human, Neanderthal, and Denisovan genomes is a catalog of potential regulatory motifs in these three human subspecies. We examined the distributions of these motifs in genetic elements including human retroviruses, human accelerated regions, and human accelerated conserved noncoding sequences regions. Differences in these distributions could be the origin of differences in phenotype between the three subspecies. Twenty significant motifs common to all three genomes were found; thirty-three were found in endogenous retroviruses in Neanderthal and Denisovan. Ten of these motifs mapped to the 22 bp core of MiR-1304. The core of this genetic element regulates the ENAM and AMTN genes, which take part in odontogenesis and whose 3' UTRs contained significant motifs. The introns of 20 genes were found to contain a large number of significant motifs, which were also overrepresented in 49 human accelerated regions. These genes include NAV2, SorCS2, TRAPPC9, GRID1, PRDM16, CAMTA1, and ASIC which are all involved in neuroregulation. Further analysis of these genes using the GO database indicates that many are associated with neurodevelopment. Also, varying numbers of significant motifs were found to occur in regions of the Neanderthal and Denisovan genomes that are missing from the human genome, suggesting further functional differences between modern and archaic humans.

CONCLUSION

Although Neanderthal and Denisovan are now extinct, detailed examination of elements from their genomes can shed light on possible phenotypic and cognitive differences between these two archaic human subspecies and modern humans. Genetic similarities and differences between these three subspecies and other fossil hominids would also be of interest.

Transcriptional-Readthrough RNAs Reflect the Phenomenon of "A Gene Contains Gene(s)" or "Gene(s) within a Gene" in the Human Genome, and Thus Are Not Chimeric RNAs

Tens of thousands of chimeric RNAs, i.e., RNAs with sequences of two genes, have been identified in human cells. Most of them are formed by two neighboring genes on the same chromosome and are considered to be derived via transcriptional readthrough, but a true readthrough event still awaits more evidence and trans-splicing that joins two transcripts together remains as a possible mechanism. We regard those genomic loci that are transcriptionally read through as unannotated genes, because their transcriptional and posttranscriptional regulations are the same as those of already-annotated genes, including fusion genes formed due to genetic alterations. Therefore, readthrough RNAs and fusion-gene-derived RNAs are not chimeras. Only those two-gene RNAs formed at the RNA level, likely via trans-splicing, without corresponding genes as genomic parents, should be regarded as authentic chimeric RNAs. However, since in human cells, procedural and mechanistic details of trans-splicing have never been disclosed, we doubt the existence of trans-splicing. Therefore, there are probably no authentic chimeras in humans, after readthrough and fusion-gene derived RNAs are all put back into the group of ordinary RNAs. Therefore, it should be further determined whether in human cells all two-neighboring-gene RNAs are derived from transcriptional readthrough and whether trans-splicing truly exists.

NCLscan: accurate identification of non-co-linear transcripts (fusion, trans-splicing and circular RNA) with a good balance between sensitivity and precision

Analysis of RNA-seq data often detects numerous ‘non-co-linear’ (NCL) transcripts, which comprised sequence segments that are topologically inconsistent with their corresponding DNA sequences in the reference genome. However, detection of NCL transcripts involves two major challenges: removal of false positives arising from alignment artifacts and discrimination between different types of NCL transcripts (trans-spliced, circular or fusion transcripts). Here, we developed a new NCL-transcript-detecting method (‘NCLscan’), which utilized a stepwise alignment strategy to almost completely eliminate false calls (>98% precision) without sacrificing true positives, enabling NCLscan outperform 18 other publicly-available tools (including fusion- and circular-RNA-detecting tools) in terms of sensitivity and precision, regardless of the generation strategy of simulated dataset, type of intragenic or intergenic NCL event, read depth of coverage, read length or expression level of NCL transcript. With the high accuracy, NCLscan was applied to distinguishing between trans-spliced, circular and fusion transcripts on the basis of poly(A)- and nonpoly(A)-selected RNA-seq data. We showed that circular RNAs were expressed more ubiquitously, more abundantly and less cell type-specifically than trans-spliced and fusion transcripts. Our study thus describes a robust pipeline for the discovery of NCL transcripts, and sheds light on the fundamental biology of these non-canonical RNA events in human transcriptome.

An intriguing RNA species--perspectives of circularized RNA

Circular RNAs (circRNAs), a kind of covalently closed RNA molecule, were used to be considered a type of by-products of mis-splicing events and were discovered sporadically due to the technological limits in the early years. With the great technological progress such as high-throughput next-generation sequencing, numerous circRNAs have recently been detected in many species. CircRNAs were expressed in a spatio-temporally specific manner, suggesting their regulatory functional potentials were overlooked previously. Intriguingly, some circRNAs were indeed found with critical physiological functions in certain circumstances. CircRNAs have a more stable molecular structure that can resist to exoribonuclease comparing to those linear ones, and their molecular functions include microRNA sponge, regulatory roles in transcription, mRNA traps that compete with linear splicing, templates for translation and possibly other presently unknown roles. Here, we review the discovery and characterization of circRNAs, the origination and formation mechanism, the physiological functions and the molecular roles, along with the methods for detection of circRNAs. We further look into the future and propose key questions to be answered for these magical RNA molecules.

Exon-intron circular RNAs regulate transcription in the nucleus

Noncoding RNAs (ncRNAs) have numerous roles in development and disease, and one of the prominent roles is to regulate gene expression. A vast number of circular RNAs (circRNAs) have been identified, and some have been shown to function as microRNA sponges in animal cells. Here, we report a class of circRNAs associated with RNA polymerase II in human cells. In these circRNAs, exons are circularized with introns 'retained' between exons; we term them exon-intron circRNAs or EIciRNAs. EIciRNAs predominantly localize in the nucleus, interact with U1 snRNP and promote transcription of their parental genes. Our findings reveal a new role for circRNAs in regulating gene expression in the nucleus, in which EIciRNAs enhance the expression of their parental genes in cis, and highlight a regulatory strategy for transcriptional control via specific RNA-RNA interaction between U1 snRNA and EIciRNAs.

Is an observed non-co-linear RNA product spliced in trans, in cis or just in vitro?

Global transcriptome investigations often result in the detection of an enormous number of transcripts composed of non-co-linear sequence fragments. Such ‘aberrant’ transcript products may arise from post-transcriptional events or genetic rearrangements, or may otherwise be false positives (sequencing/alignment errors or in vitro artifacts). Moreover, post-transcriptionally non-co-linear (‘PtNcl’) transcripts can arise from trans-splicing or back-splicing in cis (to generate so-called ‘circular RNA’). Here, we collected previously-predicted human non-co-linear RNA candidates, and designed a validation procedure integrating in silico filters with multiple experimental validation steps to examine their authenticity. We showed that >50% of the tested candidates were in vitro artifacts, even though some had been previously validated by RT-PCR. After excluding the possibility of genetic rearrangements, we distinguished between trans-spliced and circular RNAs, and confirmed that these two splicing forms can share the same non-co-linear junction. Importantly, the experimentally-confirmed PtNcl RNA events and their corresponding PtNcl splicing types (i.e. trans-splicing, circular RNA, or both sharing the same junction) were all expressed in rhesus macaque, and some were even expressed in mouse. Our study thus describes an essential procedure for confirming PtNcl transcripts, and provides further insight into the evolutionary role of PtNcl RNA events, opening up this important, but understudied, class of post-transcriptional events for comprehensive characterization.

Circular RNA is expressed across the eukaryotic tree of life

An unexpectedly large fraction of genes in metazoans (human, mouse, zebrafish, worm, fruit fly) express high levels of circularized RNAs containing canonical exons. Here we report that circular RNA isoforms are found in diverse species whose most recent common ancestor existed more than one billion years ago: fungi (Schizosaccharomyces pombe and Saccharomyces cerevisiae), a plant (Arabidopsis thaliana), and protists (Plasmodium falciparum and Dictyostelium discoideum). For all species studied to date, including those in this report, only a small fraction of the theoretically possible circular RNA isoforms from a given gene are actually observed. Unlike metazoans, Arabidopsis, D. discoideum, P. falciparum, S. cerevisiae, and S. pombe have very short introns (∼100 nucleotides or shorter), yet they still produce circular RNAs. A minority of genes in S. pombe and P. falciparum have documented examples of canonical alternative splicing, making it unlikely that all circular RNAs are by-products of alternative splicing or ‘piggyback’ on signals used in alternative RNA processing. In S. pombe, the relative abundance of circular to linear transcript isoforms changed in a gene-specific pattern during nitrogen starvation. Circular RNA may be an ancient, conserved feature of eukaryotic gene expression programs.

Integrative transcriptome sequencing identifies trans-splicing events with important roles in human embryonic stem cell pluripotency

Trans-splicing is a post-transcriptional event that joins exons from separate pre-mRNAs. Detection of trans-splicing is usually severely hampered by experimental artifacts and genetic rearrangements. Here, we develop a new computational pipeline, TSscan, which integrates different types of high-throughput long-/short-read transcriptome sequencing of different human embryonic stem cell (hESC) lines to effectively minimize false positives while detecting trans-splicing. Combining TSscan screening with multiple experimental validation steps revealed that most chimeric RNA products were platform-dependent experimental artifacts of RNA sequencing. We successfully identified and confirmed four trans-spliced RNAs, including the first reported trans-spliced large intergenic noncoding RNA (“tsRMST”). We showed that these trans-spliced RNAs were all highly expressed in human pluripotent stem cells and differentially expressed during hESC differentiation. Our results further indicated that tsRMST can contribute to pluripotency maintenance of hESCs by suppressing lineage-specific gene expression through the recruitment of NANOG and the PRC2 complex factor, SUZ12. Taken together, our findings provide important insights into the role of trans-splicing in pluripotency maintenance of hESCs and help to facilitate future studies into trans-splicing, opening up this important but understudied class of post-transcriptional events for comprehensive characterization.

Dexamethasone partially rescues ataxia telangiectasia-mutated (ATM) deficiency in ataxia telangiectasia by promoting a shortened protein variant retaining kinase activity

Ataxia telangiectasia (AT) is a rare genetic disease, still incurable, resulting from biallelic mutations in the ataxia telangiectasia-mutated (ATM) gene. Recently, short term treatment with glucocorticoid analogues improved neurological symptoms characteristic of this syndrome. Nevertheless, the molecular mechanism involved in glucocorticoid action in AT patients is not yet known. Here we describe, for the first time in mammalian cells, a short direct repeat-mediated noncanonical splicing event induced by dexamethasone, which leads to the skipping of mutations upstream of nucleotide residue 8450 of ATM coding sequence. The resulting transcript provides an alternative ORF translated in a new ATM variant with the complete kinase domain. This miniATM variant was also highlighted in lymphoblastoid cell lines from AT patients and was shown to be likely active. In conclusion, dexamethasone treatment may partly restore ATM activity in ataxia telangiectasia cells by a new molecular mechanism that overcomes most of the mutations so far described within this gene.

Origin of a novel protein-coding gene family with similar signal sequence in Schistosoma japonicum

BACKGROUND

Evolution of novel protein-coding genes is the bedrock of adaptive evolution. Recently, we identified six protein-coding genes with similar signal sequence from Schistosoma japonicum egg stage mRNA using signal sequence trap (SST). To find the mechanism underlying the origination of these genes with similar core promoter regions and signal sequence, we adopted an integrated approach utilizing whole genome, transcriptome and proteome database BLAST queries, other bioinformatics tools, and molecular analyses.

RESULTS

Our data, in combination with database analyses showed evidences of expression of these genes both at the mRNA and protein levels exclusively in all developmental stages of S. japonicum. The signal sequence motif was identified in 27 distinct S. japonicum UniGene entries with multiple mRNA transcripts, and in 34 genome contigs distributed within 18 scaffolds with evidence of genome-wide dispersion. No homolog of these genes or similar domain was found in deposited data from any other organism. We observed preponderance of flanking repetitive elements (REs), albeit partial copies, especially of the RTE-like and Perere class at either side of the duplication source locus. The role of REs as major mediators of DNA-level recombination leading to dispersive duplication is discussed with evidence from our analyses. We also identified a stepwise pathway towards functional selection in evolving genes by alternative splicing. Equally, the possible transcription models of some protein-coding representatives of the duplicons are presented with evidence of expression in vitro.

CONCLUSION

Our findings contribute to the accumulating evidence of the role of REs in the generation of evolutionary novelties in organisms' genomes.

Circular RNAs are the predominant transcript isoform from hundreds of human genes in diverse cell types

Most human pre-mRNAs are spliced into linear molecules that retain the exon order defined by the genomic sequence. By deep sequencing of RNA from a variety of normal and malignant human cells, we found RNA transcripts from many human genes in which the exons were arranged in a non-canonical order. Statistical estimates and biochemical assays provided strong evidence that a substantial fraction of the spliced transcripts from hundreds of genes are circular RNAs. Our results suggest that a non-canonical mode of RNA splicing, resulting in a circular RNA isoform, is a general feature of the gene expression program in human cells.

Genes, mutations, and human inherited disease at the dawn of the age of personalized genomics

The number of reported germline mutations in human nuclear genes, either underlying or associated with inherited disease, has now exceeded 100,000 in more than 3,700 different genes. The availability of these data has both revolutionized the study of the morbid anatomy of the human genome and facilitated "personalized genomics." With approximately 300 new "inherited disease genes" (and approximately 10,000 new mutations) being identified annually, it is pertinent to ask how many "inherited disease genes" there are in the human genome, how many mutations reside within them, and where such lesions are likely to be located? To address these questions, it is necessary not only to reconsider how we define human genes but also to explore notions of gene "essentiality" and "dispensability."Answers to these questions are now emerging from recent novel insights into genome structure and function and through complete genome sequence information derived from multiple individual human genomes. However, a change in focus toward screening functional genomic elements as opposed to genes sensu stricto will be required if we are to capitalize fully on recent technical and conceptual advances and identify new types of disease-associated mutation within noncoding regions remote from the genes whose function they disrupt.

A conserved unusual posttranscriptional processing mediated by short, direct repeated (SDR) sequences in plants

In several stress responsive gene loci of monocot cereal crops, we have previously identified an unusual posttranscriptional processing mediated by paired presence of short direct repeated (SDR) sequences at 5' and 3' splicing junctions that are distinct from conventional (U2/U12-type) splicing boundaries. By using the known SDR-containing sequences as probes, 24 plant candidate genes involved in diverse functional pathways from both monocots and dicots that potentially possess SDR-mediated posttranscriptional processing were predicted in the GenBank database. The SDRs-mediated posttranscriptional processing events including cis- and trans-actions were experimentally detected in majority of the predicted candidates. Extensive sequence analysis demonstrates several types of SDR-associated splicing peculiarities including partial exon deletion, exon fragment repetition, exon fragment scrambling and trans-splicing that result in either loss of partial exon or unusual exonic sequence rearrangements within or between RNA molecules. In addition, we show that the paired presence of SDR is necessary but not sufficient in SDR-mediated splicing in transient expression and stable transformation systems. We also show prokaryote is incapable of SDR-mediated premRNA splicing.

Entropic contributions to the splicing process

It has been recently argued that depletion attraction may play an important role in different aspects of cellular organization, ranging from the organization of transcriptional activity in transcription factories to the formation of nuclear bodies. In this paper, we suggest a new application of these ideas in the context of the splicing process, a crucial step of messenger RNA maturation in eukaryotes. We shall show that entropy effects and the resulting depletion attraction may explain the relevance of the aspecific intron length variable in the choice of splice-site recognition modality. On top of that, some qualitative features of the genome architecture of higher eukaryotes can find evolutionary realistic motivation in the light of our model.

Bioinformatics detection of alternative splicing

In recent years, genome-wide detection of alternative splicing based on Expressed Sequence Tag (EST) sequence alignments with mRNA and genomic sequences has dramatically expanded our understanding of the role of alternative splicing in functional regulation. This chapter reviews the data, methodology, and technical challenges of these genome-wide analyses of alternative splicing, and briefly surveys some of the uses to which such alternative splicing databases have been put. For example, with proper alternative splicing database schema design, it is possible to query genome-wide for alternative splicing patterns that are specific to particular tissues, disease states (e.g., cancer), gender, or developmental stages. EST alignments can be used to estimate exon inclusion or exclusion level of alternatively spliced exons and evolutionary changes for various species can be inferred from exon inclusion level. Such databases can also help automate design of probes for RT-PCR and microarrays, enabling high throughput experimental measurement of alternative splicing.

Patterned sequence in the transcriptome of vascular plants

BACKGROUND

Microsatellites (repeated subsequences based on motifs of one to six nucleotides) are widely used as codominant genetic markers because of their frequent polymorphism and relative selective neutrality. Minisatellites are repeats of motifs having seven or more nucleotides. The large number of EST sequences now available in public databases offers an opportunity to compare microsatellite and minisatellite properties and evaluate their evolution over a broad range of plant taxa.

RESULTS

Repeated motifs from one to 250 nucleotides long were identified in 6793306 expressed sequence tags (ESTs) from 88 genera of vascular plants, using a custom data-processing pipeline that allowed limited variation among repeats. The pipeline processed trimmed but otherwise unfiltered sequence and output nonredundant loci of at least 15 nucleotides, with degree of polymorphism and PCR primers wherever possible. Motifs that were an integral multiple of three in length were more abundant and richer in G/C than other motifs. From 80 to 85% of minisatellite motifs represented repeats within proteins, up to the 228-nucleotide repeat of ubiquitin, but not all of these repeats preserved reading frame. The remaining 15 to 20% of minisatellite motifs were associated with transcribed repetitive elements, e.g., retrotransposons. Relative microsatellite motif frequencies did not correlate tightly to phylogenetic relationship. Evolution of increased microsatellite and EST GC content was evident within the grasses. Microsatellites were less frequent in the transcriptome of genera with large genomes, but there was no evidence for greater dilution of the transcriptome with transposable element transcripts in these genera.

CONCLUSION

The relatively low correlation of microsatellite spectrum to phylogeny suggests that repeat loci evolve more rapidly than the surrounding sequence, although tissue specificity of the different EST libraries is a complicating factor. In-frame motifs are more abundant and higher in GC than frame-shifting motifs, but most EST minisatellite loci appear to represent repeats in translated sequence, regardless of whether reading frame is preserved. Motifs of four to six nucleotides are as polymorphic in EST collections as the commonly used motifs of two and three nucleotides, and they can be exploited as genetic markers with little additional effort.

The bursicon gene in mosquitoes: an unusual example of mRNA trans-splicing

The bursicon gene in Anopheles gambiae is encoded by two loci. Burs124 on chromosome arm 2L contains exons 1, 2, and 4, while burs3 on arm 2R contains exon 3. Exon 3 is efficiently spliced into position in the mature transcript. This unusual gene arrangement is ancient within mosquitoes, being shared by Aedes aegypti and Culex pipiens.

Mammalian transcription in support of hybrid mRNA and protein synthesis in testis and lung

Post-transcriptional mechanisms including differential splicing expand the protein repertoire beyond that provided by the one gene-one protein model. Trans-splicing has been observed in mammalian systems but is low level (sometimes referred to as noise), and a contribution to hybrid protein expression is unclear. In the study of rat sperm tail proteins a cDNA, called 1038, was isolated representing a hybrid mRNA derived in part from the ornithine decarboxylase antizyme 3 (Oaz3) gene located on rat chromosome 2 fused to sequences encoded by a novel gene on chromosome 4. Cytoplasmic Oaz3 mRNA is completely testis specific. However, in several tissues Oaz3 is transcribed and contributes to hybrid 1038 mRNA synthesis, without concurrent Oaz3 mRNA synthesis. 1038 mRNA directs synthesis of a hybrid 14-kDa protein, part chromosome 2- and part chromosome 4-derived as shown in vitro and in transfected cells. Antisera that recognize a chromosome 4-encoded C-terminal peptide confirm the hybrid character of endogenous 14-kDa protein and its presence in sperm tail structures and 1038-positive tissue. Our data suggest that the testis-specific OAZ3 gene may be an example of a mammalian gene that in several tissues is transcribed to contribute to a hybrid mRNA and protein. This finding expands the repertoire of known mechanisms available to cells to generate proteome diversity.