Identification of recurrent regulated alternative splicing events across human solid tumors

Cancer is a complex disease that involves aberrant gene expression regulation. Discriminating the modified expression patterns driving tumor biology from the many that have no or little contribution is important for understanding cancer molecular basis. Recurrent deregulation patterns observed in multiple cancer types are enriched for such driver events. Here, we studied splicing alterations in hundreds of matched tumor and normal RNA-seq samples of eight solid cancer types. We found hundreds of cassette exons for which splicing was altered in multiple cancer types and identified a set of highly frequent altered splicing events. Specific splicing regulators, including RBFOX2, MBNL1/2 and QKI, appear to account for many splicing alteration events in multiple cancer types. Together, our results provide a first global analysis of regulated splicing alterations in cancer and identify common events with a potential causative role in solid tumor development.

DNA and RNA editing of retrotransposons accelerate mammalian genome evolution

Genome evolution is commonly viewed as a gradual process that is driven by random mutations that accumulate over time. However, DNA- and RNA-editing enzymes have been identified that can accelerate evolution by actively modifying the genomically encoded information. The apolipoprotein B mRNA editing enzymes, catalytic polypeptide-like (APOBECs) are potent restriction factors that can inhibit retroelements by cytosine-to-uridine editing of retroelement DNA after reverse transcription. In some cases, a retroelement may successfully integrate into the genome despite being hypermutated. Such events introduce unique sequences into the genome and are thus a source of genomic innovation. adenosine deaminases that act on RNA (ADARs) catalyze adenosine-to-inosine editing in double-stranded RNA, commonly formed by oppositely oriented retroelements. The RNA editing confers plasticity to the transcriptome by generating many transcript variants from a single genomic locus. If the editing produces a beneficial variant, the genome may maintain the locus that produces the RNA-edited transcript for its novel function. Here, we discuss how these two powerful editing mechanisms, which both target inserted retroelements, facilitate expedited genome evolution.

The majority of transcripts in the squid nervous system are extensively recoded by A-to-I RNA editing

RNA editing by adenosine deamination alters genetic information from the genomic blueprint. When it recodes mRNAs, it gives organisms the option to express diverse, functionally distinct, protein isoforms. All eumetazoans, from cnidarians to humans, express RNA editing enzymes. However, transcriptome-wide screens have only uncovered about 25 transcripts harboring conserved recoding RNA editing sites in mammals and several hundred recoding sites in Drosophila. These studies on few established models have led to the general assumption that recoding by RNA editing is extremely rare. Here we employ a novel bioinformatic approach with extensive validation to show that the squid Doryteuthis pealeii recodes proteins by RNA editing to an unprecedented extent. We identify 57,108 recoding sites in the nervous system, affecting the majority of the proteins studied. Recoding is tissue-dependent, and enriched in genes with neuronal and cytoskeletal functions, suggesting it plays an important role in brain physiology.

DOI:http://dx.doi.org/10.7554/eLife.05198.001

For living cells to create a protein, a genetic code found in its DNA must first be ‘transcribed’ to create a corresponding molecule of messenger RNA (mRNA). DNA and RNA are both made from smaller molecules called nucleotides that are linked together into long chains; the information in both DNA and RNA is contained in the sequence of these molecules. The mRNA nucleotides coding for proteins are ‘translated’ in groups of three, and most of these nucleotide triplets instruct for a specific amino acid to be added to the newly forming protein.

DNA sequences were thought to exactly correspond with the sequence of amino acids in the resulting protein. However, it is now known that processes called RNA editing can change the nucleotide sequence of the mRNA molecules after they have been transcribed from the DNA. One such editing process, called A-to-I editing, alters the ‘A’ nucleotide so that the translation machinery reads it as a ‘G’ nucleotide instead. In some—but not all—cases, this event will change, or ‘recode’, the amino acid encoded by this stretch of mRNA, which may change how the protein behaves. This ability to create a range of proteins from a single DNA sequence could help organisms to evolve new traits.

Evidence of amino acid recoding has only been found to a very limited extent in the few species investigated so far. There has been some evidence that suggests that recoding might occur more often, and alter more proteins, in squids and octopuses. However, this could not be confirmed as the genomes of these species have not been sequenced, and these sequences were required to investigate RNA recoding using existing techniques.

Alon et al. have now developed a new approach that allows the recoding sites to be identified in organisms whose genomes have not been sequenced. Using this technique—which compares mRNA sequences with the DNA sequence they have been transcribed from—to examine the squid nervous system revealed over 57,000 recoding sites where an ‘A’ nucleotide had been modified to ‘G’ and thereby changed the coded amino acid. Many of the identified mRNA molecules had been recoded in more than one place, and many more of these than expected changed the amino acid sequence of the protein translated from them. Alon et al. therefore suggest that RNA editing may have been crucial in the evolution of the squid's nervous system, and suggest that recoding should be considered a normal part of the process used by squids to make proteins.

DOI:http://dx.doi.org/10.7554/eLife.05198.002

Stem cells. m6A mRNA methylation facilitates resolution of naïve pluripotency toward differentiation

Naïve and primed pluripotent states retain distinct molecular properties, yet limited knowledge exists on how their state transitions are regulated. Here, we identify Mettl3, an N(6)-methyladenosine (m(6)A) transferase, as a regulator for terminating murine naïve pluripotency. Mettl3 knockout preimplantation epiblasts and naïve embryonic stem cells are depleted for m(6)A in mRNAs, yet are viable. However, they fail to adequately terminate their naïve state and, subsequently, undergo aberrant and restricted lineage priming at the postimplantation stage, which leads to early embryonic lethality. m(6)A predominantly and directly reduces mRNA stability, including that of key naïve pluripotency-promoting transcripts. This study highlights a critical role for an mRNA epigenetic modification in vivo and identifies regulatory modules that functionally influence naïve and primed pluripotency in an opposing manner.

Does RNA editing compensate for Alu invasion of the primate genome?

One of the distinctive features of the primate genome is the Alu element, a repetitive short interspersed element, over a million highly similar copies of which account for >10% of the genome. A direct consequence of this feature is that primates' transcriptome is highly enriched in long stable dsRNA structures, the preferred target of adenosine deaminases acting on RNAs (ADARs), which are the enzymes catalyzing A-to-I RNA editing. Indeed, A-to-I editing by ADARs is extremely abundant in primates: over a hundred million editing sites exist in their genomes. However, there are few essential editing sites conserved across mammals that have maintained their editing level despite the radical change in ADAR target landscape. Here, we review and discuss the cost of having an unusual amount of dsRNA and editing in the transcriptome, as well as the opportunities it presents, which might have contributed to the accelerated evolution of the primates.

Genome-wide analysis of Alu editability

A-to-I RNA editing is apparently the most abundant post-transcriptional modification in primates. Virtually all editing sites reside within the repetitive Alu SINEs. Alu sequences are the dominant repeats in the human genome and thus are likely to pair with neighboring reversely oriented repeats and form double-stranded RNA structures that are bound by ADAR enzymes. Editing levels vary considerably between different adenosine sites within Alu repeats. Part of the variability has been explained by local sequence and structural motifs. Here, we focus on global characteristics that affect the editability at the Alu level. We use large RNA-seq data sets to analyze the editing levels in 203 798 Alu repeats residing within human genes. The most important factor affecting Alu editability is its distance to the closest reversely oriented neighbor-average editability decays exponentially with this distance, with a typical distance of ∼800 bp. This effect alone accounts for 28% of the total variance in editability. In addition, the number of Alu repeats of the same and reverse strand in the genomic vicinity, the expressed strand of the Alu, Alu's length and subfamily and the occurrence of reversely oriented neighbor in the same intron\exon all contribute, to a lesser extent, to the Alu editability.

Mammalian conserved ADAR targets comprise only a small fragment of the human editosome

BACKGROUND

ADAR proteins are among the most extensively studied RNA binding proteins. They bind to their target and deaminate specific adenosines to inosines. ADAR activity is essential, and the editing of a subset of their targets is critical for viability. Recently, a huge number of novel ADAR targets were detected by analyzing next generation sequencing data. Most of these novel editing sites are located in lineage-specific genomic repeats, probably a result of overactivity of editing enzymes, thus masking the functional sites. In this study we aim to identify the set of mammalian conserved ADAR targets.

RESULTS

We used RNA sequencing data from human, mouse, rat, cow, opossum, and platypus to define the conserved mammalian set of ADAR targets. We found that the conserved mammalian editing sites are surprisingly small in number and have unique characteristics that distinguish them from non-conserved ones. The sites that constitute the set have a distinct genomic distribution, tend to be located in genes encoding neurotransmitter receptors or other synapse related proteins, and have higher editing and expression levels. We also found a high consistency of editing levels of this set within mice strains and between human and mouse. Tight regulation of editing in these sites across strains and species implies their functional importance.

CONCLUSIONS

Despite the discovery of numerous editing targets, only a small number of them are conserved within mammalian evolution. These sites are extremely highly conserved and exhibit unique features, such as tight regulation, and probably play a pivotal role in mammalian biology.

A-to-I RNA editing occurs at over a hundred million genomic sites, located in a majority of human genes

RNA molecules transmit the information encoded in the genome and generally reflect its content. Adenosine-to-inosine (A-to-I) RNA editing by ADAR proteins converts a genomically encoded adenosine into inosine. It is known that most RNA editing in human takes place in the primate-specific Alu sequences, but the extent of this phenomenon and its effect on transcriptome diversity are not yet clear. Here, we analyzed large-scale RNA-seq data and detected ∼1.6 million editing sites. As detection sensitivity increases with sequencing coverage, we performed ultradeep sequencing of selected Alu sequences and showed that the scope of editing is much larger than anticipated. We found that virtually all adenosines within Alu repeats that form double-stranded RNA undergo A-to-I editing, although most sites exhibit editing at only low levels (<1%). Moreover, using high coverage sequencing, we observed editing of transcripts resulting from residual antisense expression, doubling the number of edited sites in the human genome. Based on bioinformatic analyses and deep targeted sequencing, we estimate that there are over 100 million human Alu RNA editing sites, located in the majority of human genes. These findings set the stage for exploring how this primate-specific massive diversification of the transcriptome is utilized.

Differential inhibition of LINE1 and LINE2 retrotransposition by vertebrate AID/APOBEC proteins

Background

The role of AID/APOBEC proteins in the mammalian immune response against retroviruses and retrotransposons is well established. G to A hypermutations, the hallmark of their cytidine deaminase activity, are present in several mammalian retrotransposons. However, the role of AID/APOBEC proteins in non-mammalian retroelement restriction is not completely understood.

Results

Here we provide the first evidence of anti-retroelement activity of a reptilian APOBEC protein. The green anole lizard A1 protein displayed potent DNA mutator activity and inhibited ex vivo retrotransposition of LINE1 and LINE2 ORF1 protein encoding elements, displaying a mechanism of action similar to that of the human A1 protein. In contrast, the human A3 proteins did not require ORF1 protein to inhibit LINE retrotransposition, suggesting a differential mechanism of anti-LINE action of A1 proteins, which emerged in amniotes, and A3 proteins, exclusive to placental mammals. In accordance, genomic analyses demonstrate differential G to A DNA editing of LINE retrotransposons in the lizard genome, which is also the first evidence for G to A DNA editing in non-mammalian genomes.

Conclusion

Our data suggest that vertebrate APOBEC proteins differentially inhibit the retrotransposition of LINE elements and that the anti-retroelement activity of APOBEC proteins predates mammals.

Human housekeeping genes, revisited

Housekeeping genes are involved in basic cell maintenance and, therefore, are expected to maintain constant expression levels in all cells and conditions. Identification of these genes facilitates exposure of the underlying cellular infrastructure and increases understanding of various structural genomic features. In addition, housekeeping genes are instrumental for calibration in many biotechnological applications and genomic studies. Advances in our ability to measure RNA expression have resulted in a gradual increase in the number of identified housekeeping genes. Here, we describe housekeeping gene detection in the era of massive parallel sequencing and RNA-seq. We emphasize the importance of expression at a constant level and provide a list of 3804 human genes that are expressed uniformly across a panel of tissues. Several exceptionally uniform genes are singled out for future experimental use, such as RT-PCR control genes. Finally, we discuss both ways in which current technology can meet some of past obstacles encountered, and several as yet unmet challenges.

Activity-dependent A-to-I RNA editing in rat cortical neurons

Changes in neural activity influence synaptic plasticity/scaling, gene expression, and epigenetic modifications. We present the first evidence that short-term and persistent changes in neural activity can alter adenosine-to-inosine (A-to-I) RNA editing, a post-transcriptional site-specific modification found in several neuron-specific transcripts. In rat cortical neuron cultures, activity-dependent changes in A-to-I RNA editing in coding exons are present after 6 hr of high potassium depolarization but not after 1 hr and require calcium entry into neurons. When treatments are extended from hours to days, we observe a negative feedback phenomenon: Chronic depolarization increases editing at many sites and chronic silencing decreases editing. We present several different modulations of neural activity that change the expression of different mRNA isoforms through editing.

Identification of widespread ultra-edited human RNAs

Adenosine-to-inosine modification of RNA molecules (A-to-I RNA editing) is an important mechanism that increases transciptome diversity. It occurs when a genomically encoded adenosine (A) is converted to an inosine (I) by ADAR proteins. Sequencing reactions read inosine as guanosine (G); therefore, current methods to detect A-to-I editing sites align RNA sequences to their corresponding DNA regions and identify A-to-G mismatches. However, such methods perform poorly on RNAs that underwent extensive editing ("ultra"-editing), as the large number of mismatches obscures the genomic origin of these RNAs. Therefore, only a few anecdotal ultra-edited RNAs have been discovered so far. Here we introduce and apply a novel computational method to identify ultra-edited RNAs. We detected 760 ESTs containing 15,646 editing sites (more than 20 sites per EST, on average), of which 13,668 are novel. Ultra-edited RNAs exhibit the known sequence motif of ADARs and tend to localize in sense strand Alu elements. Compared to sites of mild editing, ultra-editing occurs primarily in Alu-rich regions, where potential base pairing with neighboring, inverted Alus creates particularly long double-stranded RNA structures. Ultra-editing sites are underrepresented in old Alu subfamilies, tend to be non-conserved, and avoid exons, suggesting that ultra-editing is usually deleterious. A possible biological function of ultra-editing could be mediated by non-canonical splicing and cleavage of the RNA near the editing sites.

Consistent levels of A-to-I RNA editing across individuals in coding sequences and non-conserved Alu repeats

BACKGROUND

Adenosine to inosine (A-to-I) RNA-editing is an essential post-transcriptional mechanism that occurs in numerous sites in the human transcriptome, mainly within Alu repeats. It has been shown to have consistent levels of editing across individuals in a few targets in the human brain and altered in several human pathologies. However, the variability across human individuals of editing levels in other tissues has not been studied so far.

RESULTS

Here, we analyzed 32 skin samples, looking at A-to-I editing level in three genes within coding sequences and in the Alu repeats of six different genes. We observed highly consistent editing levels across different individuals as well as across tissues, not only in coding targets but, surprisingly, also in the non evolutionary conserved Alu repeats.

CONCLUSIONS

Our findings suggest that A-to-I RNA-editing of Alu elements is a tightly regulated process and, as such, might have been recruited in the course of primate evolution for post-transcriptional regulatory mechanisms.

Adenosine-to-inosine RNA editing shapes transcriptome diversity in primates

Human and chimpanzee genomes are almost identical, yet humans express higher brain capabilities. Deciphering the basis for this superiority is a long sought-after challenge. Adenosine-to-inosine (A-to-I) RNA editing is a widespread modification of the transcriptome. The editing level in humans is significantly higher compared with nonprimates, due to exceptional editing within the primate-specific Alu sequences, but the global editing level of nonhuman primates has not been studied so far. Here we report the sequencing of transcribed Alu sequences in humans, chimpanzees, and rhesus monkeys. We found that, on average, the editing level in the transcripts analyzed is higher in human brain compared with nonhuman primates, even where the genomic Alu structure is unmodified. Correlated editing is observed for pairs and triplets of specific adenosines along the Alu sequences. Moreover, new editable species-specific Alu insertions, subsequent to the human-chimpanzee split, are significantly enriched in genes related to neuronal functions and neurological diseases. The enhanced editing level in the human brain and the association with neuronal functions both hint at the possible contribution of A-to-I editing to the development of higher brain function. We show here that combinatorial editing is the most significant contributor to the transcriptome repertoire and suggest that Alu editing adapted by natural selection may therefore serve as an alternate information mechanism based on the binary A/I code.

Age-related gene-specific changes of A-to-I mRNA editing in the human brain

A-to-I editing is an adenosine-to-inosine modification of mRNA particularly widespread in the human brain, where it affects thousands of genes. A growing body of evidence suggests that A-to-I RNA editing is necessary for normal development and maintenance in mammals and that its deficiencies contribute to a number of pathological states. In this study, we examined whether mRNA editing levels of two mRNA species, CYFIP2 and GABRA3, change with aging. CYFIP2 has been implicated in synaptic maintenance, while GABRA3 is a GABA receptor subunit, a part of the major inhibitory neurotransmitter system in the CNS. The levels of mRNA editing were assessed in cortex samples of 20 subjects 22-102 years old. The data show an age-dependent statistically significant decrease in editing in CYFIP2. GABRA3 editing remained much more stable with age, implying that age-related decline of RNA editing is gene-specific. This is the first report of age-dependent decline in A-to-I editing. Further examination of these and other vulnerable genes may reveal specific RNA editing mechanisms that contribute to the aging phenotype.

A survey of genomic traces reveals a common sequencing error, RNA editing, and DNA editing

While it is widely held that an organism's genomic information should remain constant, several protein families are known to modify it. Members of the AID/APOBEC protein family can deaminate DNA. Similarly, members of the ADAR family can deaminate RNA. Characterizing the scope of these events is challenging. Here we use large genomic data sets, such as the two billion sequences in the NCBI Trace Archive, to look for clusters of mismatches of the same type, which are a hallmark of editing events caused by APOBEC3 and ADAR. We align 603,249,815 traces from the NCBI trace archive to their reference genomes. In clusters of mismatches of increasing size, at least one systematic sequencing error dominates the results (G-to-A). It is still present in mismatches with 99% accuracy and only vanishes in mismatches at 99.99% accuracy or higher. The error appears to have entered into about 1% of the HapMap, possibly affecting other users that rely on this resource. Further investigation, using stringent quality thresholds, uncovers thousands of mismatch clusters with no apparent defects in their chromatograms. These traces provide the first reported candidates of endogenous DNA editing in human, further elucidating RNA editing in human and mouse and also revealing, for the first time, extensive RNA editing in Xenopus tropicalis. We show that the NCBI Trace Archive provides a valuable resource for the investigation of the phenomena of DNA and RNA editing, as well as setting the stage for a comprehensive mapping of editing events in large-scale genomic datasets.

Sequence based identification of RNA editing sites

RNA editing diversifies the human transcriptome beyond the genomic repertoire. Recent years have proven a strategy based on genomics and computational sequence analysis as a powerful tool for identification and characterization of RNA editing. In particular, analysis of the human transcriptome has resulted in the identification of thousands of A-to-I editing sites within genomic repeats, as well as a few hundred sites located outside repeats. We review these recent advancements, emphasizing the principles underlying the various methods used. Possible directions for extending these methods are discussed.

Evidence for large diversity in the human transcriptome created by Alu RNA editing

Adenosine-to-inosine (A-to-I) RNA editing alters the original genomic content of the human transcriptome and is essential for maintenance of normal life in mammals. A-to-I editing in Alu repeats is abundant in the human genome, with many thousands of expressed Alu sequences undergoing editing. Little is known so far about the contribution of Alu editing to transcriptome complexity. Transcripts derived from a single edited Alu sequence can be edited in multiple sites, and thus could theoretically generate a large number of different transcripts. Here we explored whether the combinatorial potential nature of edited Alu sequences is actually fulfilled in the human transcriptome. We analyzed datasets of editing sites and performed an analysis of a detailed transcript set of one edited Alu sequence. We found that editing appears at many more sites than detected by earlier genomic screens. To a large extent, editing of different sites within the same transcript is only weakly correlated. Thus, rather than finding a few versions of each transcript, a large number of edited variants arise, resulting in immense transcript diversity that eclipses alternative splicing as mechanism of transcriptome diversity, although with less impact on the proteome.

Genome-wide identification of human RNA editing sites by parallel DNA capturing and sequencing

Adenosine-to-inosine (A-to-I) RNA editing leads to transcriptome diversity and is important for normal brain function. To date, only a handful of functional sites have been identified in mammals. We developed an unbiased assay to screen more than 36,000 computationally predicted nonrepetitive A-to-I sites using massively parallel target capture and DNA sequencing. A comprehensive set of several hundred human RNA editing sites was detected by comparing genomic DNA with RNAs from seven tissues of a single individual. Specificity of our profiling was supported by observations of enrichment with known features of targets of adenosine deaminases acting on RNA (ADAR) and validation by means of capillary sequencing. This efficient approach greatly expands the repertoire of RNA editing targets and can be applied to studies involving RNA editing-related human diseases.

Efficiency of complex production in changing environment

BACKGROUND

Cell function necessitates the assemblage of proteins into complexes, a process which requires further regulation on top of the fairly understood mechanisms used to control the transcription and translation of a single protein. However, not much is known about how protein levels are controlled to realize that regulation.

RESULTS

We integrated data on the composition of yeast protein complexes and the dynamics of their protein building-blocks concentrations to show how the cell regulates protein levels to optimize complex formation. We find that proteins which are subunits of the same complex tend to have similar levels which change similarly following a change in growth conditions, and that abundant proteins undergo larger decrease in their copy number when grown in minimal media. We also study the fluctuations in protein levels and find them to be significantly smaller in large complexes, and in the least abundant subunit of each complex. We use a mathematical model of complex synthesis to explain how all these observations increase the efficiency of complex synthesis, in terms of better utilization of the available molecules and better resilience to stochastic variations.

CONCLUSION

In conclusion, these results indicate an intricate regulation at all levels of protein production for the purpose of optimizing complex formation.

Intronic Alus influence alternative splicing

Examination of the human transcriptome reveals higher levels of RNA editing than in any other organism tested to date. This is indicative of extensive double-stranded RNA (dsRNA) formation within the human transcriptome. Most of the editing sites are located in the primate-specific retrotransposed element called Alu. A large fraction of Alus are found in intronic sequences, implying extensive Alu-Alu dsRNA formation in mRNA precursors. Yet, the effect of these intronic Alus on splicing of the flanking exons is largely unknown. Here, we show that more Alus flank alternatively spliced exons than constitutively spliced ones; this is especially notable for those exons that have changed their mode of splicing from constitutive to alternative during human evolution. This implies that Alu insertions may change the mode of splicing of the flanking exons. Indeed, we demonstrate experimentally that two Alu elements that were inserted into an intron in opposite orientation undergo base-pairing, as evident by RNA editing, and affect the splicing patterns of a downstream exon, shifting it from constitutive to alternative. Our results indicate the importance of intronic Alus in influencing the splicing of flanking exons, further emphasizing the role of Alus in shaping of the human transcriptome.

Altered adenosine-to-inosine RNA editing in human cancer

Adenosine-to-inosine (A-to-I) RNA editing was recently shown to be abundant in the human transcriptome, affecting thousands of genes. Employing a bioinformatic approach, we identified significant global hypoediting of Alu repetitive elements in brain, prostate, lung, kidney, and testis tumors. Experimental validation confirmed this finding, showing significantly reduced editing in Alu sequences within MED13 transcripts in brain tissues. Looking at editing of specific recoding and noncoding sites, including in cancer-related genes, a more complex picture emerged, with a gene-specific editing pattern in tumors vs. normal tissues. Additionally, we found reduced RNA levels of all three editing mediating enzymes, ADAR, ADARB1, and ADARB2, in brain tumors. The reduction of ADARB2 correlated with the grade of malignancy of glioblastoma multiforme, the most aggressive of brain tumors, displaying a 99% decrease in ADARB2 RNA levels. Consistently, overexpression of ADAR and ADARB1 in the U87 glioblastoma multiforme cell line resulted in decreased proliferation rate, suggesting that reduced A-to-I editing in brain tumors is involved in the pathogenesis of cancer. Altered epigenetic control was recently shown to play a central role in oncogenesis. We suggest that A-to-I RNA editing may serve as an additional epigenetic mechanism relevant to cancer development and progression.

RNA-editing-mediated exon evolution

BACKGROUND

Alu retroelements are specific to primates and abundant in the human genome. Through mutations that create functional splice sites within intronic Alus, these elements can become new exons in a process denoted exonization. It was recently shown that Alu elements are also heavily changed by RNA editing in the human genome.

RESULTS

Here we show that the human nuclear prelamin A recognition factor contains a primate-specific Alu-exon that exclusively depends on RNA editing for its exonization. We demonstrate that RNA editing regulates the exonization in a tissue-dependent manner, through both the creation of a functional AG 3' splice site, and alteration of functional exonic splicing enhancers within the exon. Furthermore, a premature stop codon within the Alu-exon is eliminated by an exceptionally efficient RNA editing event. The sequence surrounding this editing site is important not only for editing of that site but also for editing in other neighboring sites as well.

CONCLUSION

Our results show that the abundant RNA editing of Alu sequences can be recruited as a mechanism supporting the birth of new exons in the human genome.

RNA editing level in the mouse is determined by the genomic repeat repertoire

A-to-I RNA editing is the conversion of adenosine to inosine in double-stranded cellular and viral RNAs. Recently, abundant editing of human transcripts affecting thousands of genes has been reported. Most editing sites are confined to the primate-specific Alu repeats. Notably, the editing level in mouse was shown to be much lower. In order to find the reason for this dramatic difference, here we identify editing sites within mouse repeats and analyze the sequence properties required for RNA editing. Our results show that the overall rate of RNA editing is determined by specific properties of different repeat families such as abundance, length, and divergence. We show that the striking difference in editing levels between human and mouse is mostly due to the higher divergence of the different mouse repeats.

Letter from the editor: Adenosine-to-inosine RNA editing in Alu repeats in the human genome

Adenosine-to-inosine (A-to-I) RNA editing increases the complexity of the human transcriptome and is essential for maintenance of normal life in mammals. Most A-to-I substitutions occur within repetitive elements in the genome, mainly in Alu repeats. The phenomenon of A-to-I editing is far less abundant in mice, rats, chickens and flies than in humans, which correlates with the relative under-representation of Alu repeats in these non-primate genomes. Here, we review the recent results of bioinformatic and laboratory approaches that have estimated the extent of the editing phenomenon. We discuss the possible biological relevance of the editing pathway, its possible interaction with other cellular pathways that respond to double-stranded RNA and its possible contribution to the accelerated evolution of primates.

Connectivity and expression in protein networks: proteins in a complex are uniformly expressed

We explore the interplay between the protein-protein interactions network and the expression of the interacting proteins. It is shown that interacting proteins are expressed in significantly more similar cellular concentrations. This is largely due to interacting pairs which are part of protein complexes. We solve a generic model of complex formation and show explicitly that complexes form most efficiently when their members have roughly the same concentrations. Therefore, the observed similarity in interacting protein concentrations could be attributed to optimization for efficiency of complex formation.

Is there any sense in antisense editing?

Several recent studies have hypothesized that sense-antisense RNA-transcript pairs create dsRNA duplexes that undergo extensive A-to-I RNA editing. We studied human and mouse genomic antisense regions and found that the editing level in these areas is negligible. This observation questions the scope of sense-antisense duplexes formation in-vivo, which is the basis for several proposed regulatory mechanisms.

Algorithmic approaches for identification of RNA editing sites

Recently a number of groups have introduced computational methods for the detection of A-to-I RNA editing sites. These approaches have resulted in finding thousands of editing sites within the genomic repeats, as well as a few novel genetic recoding sites. We review these recent advancements, emphasizing the principles underlying the various methods used. Possible directions for extending these methods are discussed.

Identification of RNA editing sites in the SNP database

The relationship between human inherited genomic variations and phenotypic differences has been the focus of much research effort in recent years. These studies benefit from millions of single-nucleotide polymorphism (SNP) records available in public databases, such as dbSNP. The importance of identifying false dbSNP records increases with the growing role played by SNPs in linkage analysis for disease traits. In particular, the emerging understanding of the abundance of DNA and RNA editing calls for a careful distinction between inherited SNPs and somatic DNA and RNA modifications. In order to demonstrate that some of the SNP database records are actually somatic modification, we focus on one type of these modifications, namely A-to-I RNA editing, and present evidence for hundreds of dbSNP records that are actually editing sites. We provide a list of 102 RNA editing sites previously annotated in dbSNP database as SNPs, and experimentally validate seven of these. Interestingly, we show how dbSNP can serve as a starting point to look for new editing sites. Our results, for this particular type of RNA editing, demonstrate the need for a careful analysis of SNP databases in light of the increasing recognition of the significance of somatic sequence modifications.

From genome to antivirals: SARS as a test tube

The severe acute respiratory syndrome (SARS) epidemic brought into the spotlight the need for rapid development of effective anti-viral drugs against newly emerging viruses. Researchers have leveraged the 20-year battle against AIDS into a variety of possible treatments for SARS. Most prominently, based solely on viral genome information, silencers of viral genes, viral-enzyme blockers and viral-entry inhibitors were suggested as potential therapeutic agents for SARS. In particular, inhibitors of viral entry, comprising therapeutic peptides, were based on the recently launched anti-HIV drug enfuvirtide. This could represent one of the most direct routes from genome sequencing to the discovery of antiviral drugs.

Evidence for abundant transcription of non-coding regions in the Saccharomyces cerevisiae genome

Background

Recent studies in a growing number of organisms have yielded accumulating evidence that a significant portion of the non-coding region in the genome is transcribed. We address this issue in the yeast Saccharomyces cerevisiae.

Results

Taking into account the absence of a significantly large yeast EST database, we use microarray expression data collected for genomic regions erroneously believed to be coding to study the expression pattern of non-coding regions in the Saccharomyces cerevisiae genome. We find that at least 164 out of 589 (28%) such regions are expressed under specific biological conditions. In particular, looking at the probes that are located opposing other known genes at the same genomic locus, we find that 88 out of 341 (26%) of these genes support antisense transcription. The expression patterns of these antisense genes are positively correlated. We validate these results using RT-PCR on a sample of 6 non-coding transcripts.

Conclusion

1. The yeast genome is transcribed on a scale larger than previously assumed. 2. Correlated transcription of antisense genes is abundant in the yeast genome. 3. Antisense genes in yeast are non-coding.

Is abundant A-to-I RNA editing primate-specific?

A-to-I RNA editing is common in all eukaryotes, and is associated with various neurological functions. Recently, A-to-I editing was found to occur frequently in the human transcriptome. In this article, we show that the frequency of A-to-I editing in humans is at least an order of magnitude higher than in the mouse, rat, chicken or fly genomes. The extraordinary frequency of RNA editing in human is explained by the dominance of the primate-specific Alu element in the human transcriptome, which increases the number of double-stranded RNA substrates.

Evolutionarily conserved human targets of adenosine to inosine RNA editing

A-to-I RNA editing by ADARs is a post-transcriptional mechanism for expanding the proteomic repertoire. Genetic recoding by editing was so far observed for only a few mammalian RNAs that are predominantly expressed in nervous tissues. However, as these editing targets fail to explain the broad and severe phenotypes of ADAR1 knockout mice, additional targets for editing by ADARs were always expected. Using comparative genomics and expressed sequence analysis, we identified and experimentally verified four additional candidate human substrates for ADAR-mediated editing: FLNA, BLCAP, CYFIP2 and IGFBP7. Additionally, editing of three of these substrates was verified in the mouse while two of them were validated in chicken. Interestingly, none of these substrates encodes a receptor protein but two of them are strongly expressed in the CNS and seem important for proper nervous system function. The editing pattern observed suggests that some of the affected proteins might have altered physiological properties leaving the possibility that they can be related to the phenotypes of ADAR1 knockout mice.

Evolution of multicellularity in Metazoa: comparative analysis of the subcellular localization of proteins in Saccharomyces, Drosophila and Caenorhabditis

A comparison of the subcellular assignments of proteins between the unicellular Saccharomyces cerevisiae and the multicellular Drosophila melanogaster and Caenorhabditis elegans was performed using a computational tool for the prediction of subcellular localization. Nine subcellular compartments were studied: (1) extracellular domain, (2) cell membrane, (3) cytoplasm, (4) endoplasmic reticulum, (5) Golgi apparatus, (6) lysosome, (7) peroxisome, (8) mitochondria, and (9) nucleus. The transition to multicellularity was found to be characterized by an increase in the total number of proteins encoded by the genome. Interestingly, this increase is distributed unevenly among the subcellular compartments. That is, a disproportionate increase in the number of proteins in the extracellular domain, the cell membrane, and the cytoplasm is observed in multicellular organisms, while no such increase is seen in other subcellular compartments. A possible explanation involves signal transduction. In terms of protein numbers, signal transduction pathways may be roughly described as a pyramid with an expansive base in the extracellular domain (the numerous extracellular signal proteins), progressively narrowing at the cell membrane and cytoplasmic levels, and ending in a narrow tip consisting of only a handful of transcription modulators in the nucleus. Our observations suggest that extracellular signaling interactions among metazoan cells account for the uneven increase in the numbers of proteins among subcellular compartments during the transition to multicellularity.

Systematic identification of abundant A-to-I editing sites in the human transcriptome

RNA editing by members of the ADAR (adenosine deaminases acting on RNA) family leads to site-specific conversion of adenosine to inosine (A-to-I) in precursor messenger RNAs. Editing by ADARs is believed to occur in all metazoa, and is essential for mammalian development. Currently, only a limited number of human ADAR substrates are known, whereas indirect evidence suggests a substantial fraction of all pre-mRNAs being affected. Here we describe a computational search for ADAR editing sites in the human transcriptome, using millions of available expressed sequences. We mapped 12,723 A-to-I editing sites in 1,637 different genes, with an estimated accuracy of 95%, raising the number of known editing sites by two orders of magnitude. We experimentally validated our method by verifying the occurrence of editing in 26 novel substrates. A-to-I editing in humans primarily occurs in noncoding regions of the RNA, typically in Alu repeats. Analysis of the large set of editing sites indicates the role of editing in controlling dsRNA stability.

Preferential attachment in the protein network evolution

The Saccharomyces cerevisiae protein-protein interaction map, as well as many natural and man-made networks, shares the scale-free topology. The preferential attachment model was suggested as a generic network evolution model that yields this universal topology. However, it is not clear that the model assumptions hold for the protein interaction network. Using a cross-genome comparison, we show that (a) the older a protein, the better connected it is, and (b) the number of interactions a protein gains during its evolution is proportional to its connectivity. Therefore, preferential attachment governs the protein network evolution. Evolutionary mechanisms leading to such preference and some implications are discussed.

Cloaked similarity between HIV-1 and SARS-CoV suggests an anti-SARS strategy

Background

Severe acute respiratory syndrome (SARS) is a febrile respiratory illness. The disease has been etiologically linked to a novel coronavirus that has been named the SARS-associated coronavirus (SARS-CoV), whose genome was recently sequenced. Since it is a member of the Coronaviridae, its spike protein (S2) is believed to play a central role in viral entry by facilitating fusion between the viral and host cell membranes. The protein responsible for viral-induced membrane fusion of HIV-1 (gp41) differs in length, and has no sequence homology with S2.

Results

Sequence analysis reveals that the two viral proteins share the sequence motifs that construct their active conformation. These include (1) an N-terminal leucine/isoleucine zipper-like sequence, and (2) a C-terminal heptad repeat located upstream of (3) an aromatic residue-rich region juxtaposed to the (4) transmembrane segment.

Conclusions

This study points to a similar mode of action for the two viral proteins, suggesting that anti-viral strategy that targets the viral-induced membrane fusion step can be adopted from HIV-1 to SARS-CoV. Recently the FDA approved Enfuvirtide, a synthetic peptide corresponding to the C-terminal heptad repeat of HIV-1 gp41, as an anti-AIDS agent. Enfuvirtide and C34, another anti HIV-1 peptide, exert their inhibitory activity by binding to a leucine/isoleucine zipper-like sequence in gp41, thus inhibiting a conformational change of gp41 required for its activation. We suggest that peptides corresponding to the C-terminal heptad repeat of the S2 protein may serve as inhibitors for SARS-CoV entry.

Widespread occurrence of antisense transcription in the human genome

An increasing number of eukaryotic genes are being found to have naturally occurring antisense transcripts. Here we study the extent of antisense transcription in the human genome by analyzing the public databases of expressed sequences using a set of computational tools designed to identify sense-antisense transcriptional units on opposite DNA strands of the same genomic locus. The resulting data set of 2,667 sense-antisense pairs was evaluated by microarrays containing strand-specific oligonucleotide probes derived from the region of overlap. Verification of specific cases by northern blot analysis with strand-specific riboprobes proved transcription from both DNA strands. We conclude that > or =60% of this data set, or approximately 1,600 predicted sense-antisense transcriptional units, are transcribed from both DNA strands. This indicates that the occurrence of antisense transcription, usually regarded as infrequent, is a very common phenomenon in the human genome. Therefore, antisense modulation of gene expression in human cells may be a common regulatory mechanism.