Evidence for a preferential targeting of 3'-UTRs by cis-encoded natural antisense transcripts

Non-random retention of protein-coding overlapping genes in Metazoa

Background

Although the overlap of transcriptional units occurs frequently in eukaryotic genomes, its evolutionary and biological significance remains largely unclear. Here we report a comparative analysis of overlaps between genes coding for well-annotated proteins in five metazoan genomes (human, mouse, zebrafish, fruit fly and worm).

Results

For all analyzed species the observed number of overlapping genes is always lower than expected assuming functional neutrality, suggesting that gene overlap is negatively selected. The comparison to the random distribution also shows that retained overlaps do not exhibit random features: antiparallel overlaps are significantly enriched, while overlaps lying on the same strand and those involving coding sequences are highly underrepresented. We confirm that overlap is mostly species-specific and provide evidence that it frequently originates through the acquisition of terminal, non-coding exons. Finally, we show that overlapping genes tend to be significantly co-expressed in a breast cancer cDNA library obtained by 454 deep sequencing, and that different overlap types display different patterns of reciprocal expression.

Conclusion

Our data suggest that overlap between protein-coding genes is selected against in Metazoa. However, when retained it may be used as a species-specific mechanism for the reciprocal regulation of neighboring genes. The tendency of overlaps to involve non-coding regions of the genes leads to the speculation that the advantages achieved by an overlapping arrangement may be optimized by evolving regulatory non-coding transcripts.

Processing of naturally occurring sense/antisense transcripts of the vertebrate Slc34a gene into short RNAs

Overlapping sense/antisense RNAs transcribed in opposite directions from the same genomic locus are common in vertebrates. The impact of antisense transcription on gene regulation and cell biology is largely unknown. We show that sense/antisense RNAs of an evolutionarily conserved phosphate transporter gene (Slc34a2a) are coexpressed in a short time window during embryonic development of zebrafish at 48 hours postfertilization (hpf). To address the mechanism by which coexpressed sense/antisense transcripts are processed, we injected sense/antisense RNAs in various combinations into Xenopus oocytes. In the cytoplasm RNAs were stable in whatever combination expressed. In the nucleus coinjected sense/antisense transcripts were degraded into short RNAs of approximately 23 bases within 4 h. A homologous transcript from toad or another isoform (Slc34a2b) from zebrafish failed to trigger processing. In oocytes that were primed with nuclear sense/antisense RNA coinjections, a reporter RNA was rapidly degraded. We produced evidence that the observed processing of complementary transcripts was not restricted to Xenopus oocytes, because Slc34a-related short RNAs were detected in zebrafish embryos by Northern blotting. Signals were observed at stages that showed coexpression of sense/antisense transcripts. Remarkably, strand-specific probes revealed that the orientation of short RNAs was developmentally regulated. In addition, RNA from zebrafish embryos 48 hpf was able to induce degradation of reporter constructs in Xenopus oocytes. Our findings may give important clues to understanding the physiological role of the widespread antisense transcription.

Osteopontin and systemic lupus erythematosus association: a probable gene-gender interaction

Osteopontin (SPP1) is an important bone matrix mediator found to have key roles in inflammation and immunity. SPP1 genetic polymorphisms and increased osteopontin protein levels have been reported to be associated with SLE in small patient collections. The present study evaluates association between SPP1 polymorphisms and SLE in a large cohort of 1141 unrelated SLE patients [707 European-American (EA) and 434 African-American (AA)], and 2009 unrelated controls (1309 EA and 700 AA). Population-based case-control association analyses were performed. To control for potential population stratification, admixture adjusted logistic regression, genomic control (GC), structured association (STRAT), and principal components analysis (PCA) were applied. Combined analysis of 2 ethnic groups, showed the minor allele of 2 SNPs (rs1126616T and rs9138C) significantly associated with higher risk of SLE in males (P = 0.0005, OR = 1.73, 95% CI = 1.28–2.33), but not in females. Indeed, significant gene-gender interactions in the 2 SNPs, rs1126772 and rs9138, were detected (P = 0.001 and P = 0.0006, respectively). Further, haplotype analysis identified rs1126616T-rs1126772A-rs9138C which demonstrated significant association with SLE in general (P = 0.02, OR = 1.30, 95%CI 1.08–1.57), especially in males (P = 0.0003, OR = 2.42, 95%CI 1.51–3.89). Subgroup analysis with single SNPs and haplotypes also identified a similar pattern of gender-specific association in AA and EA. GC, STRAT, and PCA results within each group showed consistent associations. Our data suggest SPP1 is associated with SLE, and this association is especially stronger in males. To our knowledge, this report serves as the first association of a specific autosomal gene with human male lupus.

Imprinting of an evolutionarily conserved antisense transcript gene APeg3

APeg3 is an antisense transcript gene of Peg3, which has been recently identified from rat brain. Careful analyses of EST databases indicated that a homologous transcript also exists in other mammalian species, including mouse, cow and human. 5'-and 3'-RACE experiments have subsequently identified a 900-bp cDNA sequence of APeg3 from mouse brain. Mouse APeg3 is localized in the 3'UTR of Peg3 with an intronless genomic structure. The expression of mouse APeg3 is derived mainly from the paternal allele, indicating the imprinting of this antisense transcript gene in brain. Strand-specific RNA analyses also revealed the expression of both human and cow APEG3 in adult brains. In sum, our study confirms that the mammalian PEG3 locus harbors an antisense transcript gene displaying paternal allele-specific expression, and the evolutionary conservation further suggests potential roles of this transcript gene for the function of this imprinted domain.

Antisense transcription in the human cytomegalovirus transcriptome

Human cytomegalovirus (HCMV) infections are prevalent in human populations and can cause serious diseases, especially in those with compromised or immature immune systems. The HCMV genome of 230 kb is among the largest of the herpesvirus genomes. Although the entire sequence of the laboratory-adapted AD169 strain of HCMV has been available for 18 years, the precise number of viral genes is still in question. We undertook an analysis of the HCMV transcriptome as an approach to enumerate and analyze the gene products of HCMV. Transcripts of HCMV-infected fibroblasts were isolated at different times after infection and used to generate cDNA libraries representing different temporal classes of viral genes. cDNA clones harboring viral sequences were selected and subjected to sequence analysis. Of the 604 clones analyzed, 45% were derived from genomic regions predicted to be noncoding. Additionally, at least 55% of the cDNA clones in this study were completely or partially antisense to known or predicted HCMV genes. The remarkable accumulation of antisense transcripts during infection suggests that currently available genomic maps based on open-reading-frame and other in silico analyses may drastically underestimate the true complexity of viral gene products. These findings also raise the possibility that aspects of both the HCMV life cycle and genome organization are influenced by antisense transcription. Correspondingly, virus-derived noncoding and antisense transcripts may shed light on HCMV pathogenesis and may represent a new class of targets for antiviral therapies.

Fine-tuning of the innate immune response by microRNAs

MicroRNAs (miRNAs) are emerging as potent regulators of many biological processes, including cellular differentiation and disease. Recently, miRNA has been directly involved in innate immunity and transduction signalling by Toll-like receptors and the ensuing cytokine response. In this review, we present an overview of what is currently known of the involvement of miRNA and RNA interference components in the fine-tuning of innate immune responses.

Sense-antisense pairs in mammals: functional and evolutionary considerations

Analysis of a catalog of S-AS pairs in the human and mouse genomes revealed several putative roles for natural antisense transcripts and showed that some are artifacts of cDNA library construction.

Background

A significant number of genes in mammalian genomes are being found to have natural antisense transcripts (NATs). These sense-antisense (S-AS) pairs are believed to be involved in several cellular phenomena.

Results

Here, we generated a catalog of S-AS pairs occurring in the human and mouse genomes by analyzing different sources of expressed sequences available in the public domain plus 122 massively parallel signature sequencing (MPSS) libraries from a variety of human and mouse tissues. Using this dataset of almost 20,000 S-AS pairs in both genomes we investigated, in a computational and experimental way, several putative roles that have been assigned to NATs, including gene expression regulation. Furthermore, these global analyses allowed us to better dissect and propose new roles for NATs. Surprisingly, we found that a significant fraction of NATs are artifacts produced by genomic priming during cDNA library construction.

Conclusion

We propose an evolutionary and functional model in which alternative polyadenylation and retroposition account for the origin of a significant number of functional S-AS pairs in mammalian genomes.

Hypoxia-inducible expression of a natural cis-antisense transcript inhibits endothelial nitric-oxide synthase

The destabilization of endothelial nitric-oxide synthase (eNOS) mRNA in hypoxic endothelial cells may be important in the etiology of vascular diseases, such as pulmonary hypertension. Recently, an overlapping antisense transcript to eNOS/NOS3 was implicated in the post-transcriptional regulation of eNOS. We demonstrate here that expression of sONE, also known as eNOS antisense (NOS3AS) or autophagy 9-like 2 (APG9L2), is robustly induced by hypoxia or functional deficiency of von Hippel-Lindau protein. sONE is also up-regulated in the aortas of hypoxic rats. In hypoxic endothelial cells, sONE expression negatively correlates with eNOS expression. Blocking the hypoxic induction of sONE by RNA interference attenuates the fall in both eNOS RNA and protein. We provide evidence that the induction of sONE primarily involves transcript stabilization rather than increased transcriptional activity and is von Hippel-Lindaubut not hypoxia-inducible factor 2alpha-dependent. We also demonstrate that sONE transcripts are enriched in the nucleus of normoxic cells and that hypoxia promotes an increase in the level of cytoplasmic and polyribosome-associated, sONE mRNA. The finding that eNOS expression can be regulated by an overlapping cis-antisense transcript in a stimulus-dependent fashion provides evidence that sense/antisense interactions may play a previously unappreciated role in vascular disease pathogenesis.

Genome-wide natural antisense transcription: coupling its regulation to its different regulatory mechanisms

Many genomic loci contain transcription units on both strands, therefore two oppositely oriented transcripts can overlap. Often, one strand codes for a protein, whereas the transcript from the other strand is non-encoding. Such natural antisense transcripts (NATs) can negatively regulate the conjugated sense transcript. NATs are highly prevalent in a wide range of species--for example, around 15% of human protein-encoding genes have an associated NAT. The regulatory mechanisms by which NATs act are diverse, as are the means to control their expression. Here, we review the current understanding of NAT function and its mechanistic basis, which has been gathered from both individual gene cases and genome-wide studies. In parallel, we survey findings about the regulation of NAT transcription. Finally, we hypothesize that the regulation of antisense transcription might be tailored to its mode of action. According to this model, the observed relationship between the expression patterns of NATs and their targets might indicate the regulatory mechanism that is in action.

Computational identification of microRNAs and their targets

MicroRNAs (miRNAs) are one class of newly identified riboregulators of gene expression in many eukaryotic organisms. They play important roles in multiple biological and metabolic processes, including developmental timing, signal transduction, cell maintenance and differentiation, diseases and cancers. miRNAs regulate gene expression at the posttranscriptional level by directly cleaving targeted mRNAs or repressing translation. Although the founding members of miRNAs were discovered by genetic screening approaches, experimental approaches were limited by their low efficiency, time consuming, and high cost. As an alternative, computational approaches were developed. Computational approaches for identifying miRNAs are based on the following major characteristics of miRNAs: hairpin-shaped secondary structures, high conservation for some miRNAs, and high minimal folding free energy index (MFEI). Computational approaches also play an important role in identifying miRNA targets. A majority of known miRNAs and their targets were identified by computational approaches. Several web-based or non-web-based computer software programs are publicly available for predicting miRNAs and their targets.

Expression profiling of antisense transcripts on DNA arrays

The majority of mouse genes are estimated to undergo bidirectional transcription; however, their tissue-specific distribution patterns and physiological significance are largely unknown. This is in part due to the lack of methodology to routinely assess the expression of natural antisense transcripts (NATs) on a large scale. Here we tested whether commercial DNA arrays can be used to monitor antisense transcription in mouse kidney and brain. We took advantage of the reversely annotated oligonucleotides on the U74 mouse genome array from Affymetrix that hybridize to NATs overlapping with the sense transcript in the area of the probe set. In RNA samples from mouse kidney and brain, 11.9% and 10.1%, respectively, of 5,652 potential NATs returned positive and about half of the antisense RNAs were detected in both tissues, which was similar to the fraction of sense transcripts expressed in both tissues. Notably, we found that the majority of NATs are related to the sense transcriptome since corresponding sense transcripts were detected for 92.5% (kidney) and 74.5% (brain) of the detected antisense RNAs. Antisense RNA transcription was confirmed by real-time PCR and included additional RNA samples from heart, thymus, and liver. The randomly selected transcripts showed tissue specific expression patterns and varying sense/antisense ratios. The results indicate that antisense transcriptomes are tissue specific, and although pairing of sense/antisense transcripts are known to result in rapid degradation, our data provide proof of principle that the sensitivity of commercial DNA arrays is sufficient to assess NATs in total RNA of whole organs.

Modeling a whole organ using proteomics: the avian bursa of Fabricius

While advances in proteomics have improved proteome coverage and enhanced biological modeling, modeling function in multicellular organisms requires understanding how cells interact. Here we used the chicken bursa of Fabricius, a common experimental system for B cell function, to model organ function from proteomics data. The bursa has two major functional cell types: B cells and the supporting stromal cells. We used differential detergent fractionation-multidimensional protein identification technology (DDF-MudPIT) to identify 5198 proteins from all cellular compartments. Of these, 1753 were B cell specific, 1972 were stroma specific and 1473 were shared between the two. By modeling programmed cell death (PCD), cell differentiation and proliferation, and transcriptional activation, we have improved functional annotation of chicken proteins and placed chicken-specific death receptors into the PCD process using phylogenetics. We have identified 114 transcription factors (TFs); 42 of the bursal B cell TFs have not been reported before in any B cells. We have also improved the structural annotation of a newly sequenced genome by confirming the in vivo expression of 4006 "predicted", and 6623 ab initio, ORFs. Finally, we have developed a novel method for facilitating structural annotation, "expressed peptide sequence tags" (ePSTs) and demonstrate its utility by identifying 521 potential novel proteins from the chicken "unassigned chromosome".

microRNAs as oncogenes and tumor suppressors

microRNAs (miRNAs) are a new class of non-protein-coding, endogenous, small RNAs. They are important regulatory molecules in animals and plants. miRNA regulates gene expression by translational repression, mRNA cleavage, and mRNA decay initiated by miRNA-guided rapid deadenylation. Recent studies show that some miRNAs regulate cell proliferation and apoptosis processes that are important in cancer formation. By using multiple molecular techniques, which include Northern blot analysis, real-time PCR, miRNA microarray, up- or down-expression of specific miRNAs, it was found that several miRNAs were directly involved in human cancers, including lung, breast, brain, liver, colon cancer, and leukemia. In addition, some miRNAs may function as oncogenes or tumor suppressors. More than 50% of miRNA genes are located in cancer-associated genomic regions or in fragile sites, suggesting that miRNAs may play a more important role in the pathogenesis of a limited range of human cancers than previously thought. Overexpressed miRNAs in cancers, such as mir-17-92, may function as oncogenes and promote cancer development by negatively regulating tumor suppressor genes and/or genes that control cell differentiation or apoptosis. Underexpressed miRNAs in cancers, such as let-7, function as tumor suppressor genes and may inhibit cancers by regulating oncogenes and/or genes that control cell differentiation or apoptosis. miRNA expression profiles may become useful biomarkers for cancer diagnostics. In addition, miRNA therapy could be a powerful tool for cancer prevention and therapeutics.

Genome-wide in silico identification and analysis of cis natural antisense transcripts (cis-NATs) in ten species

We developed a fast, integrative pipeline to identify cis natural antisense transcripts (cis-NATs) at genome scale. The pipeline mapped mRNAs and ESTs in UniGene to genome sequences in GoldenPath to find overlapping transcripts and combining information from coding sequence, poly(A) signal, poly(A) tail and splicing sites to deduce transcription orientation. We identified cis-NATs in 10 eukaryotic species, including 7830 candidate sense–antisense (SA) genes in 3915 SA pairs in human. The abundance of SA genes is remarkably low in worm and does not seem to be caused by the prevalence of operons. Hundreds of SA pairs are conserved across different species, even maintaining the same overlapping patterns. The convergent SA class is prevalent in fly, worm and sea squirt, but not in human or mouse as reported previously. The percentage of SA genes among imprinted genes in human and mouse is 24–47%, a range between the two previous reports. There is significant shortage of SA genes on Chromosome X in human and mouse but not in fly or worm, supporting X-inactivation in mammals as a possible cause. SA genes are over-represented in the catalytic activities and basic metabolism functions. All candidate cis-NATs can be downloaded from .

Evidence for variation in abundance of antisense transcripts between multicellular animals but no relationship between antisense transcriptionand organismic complexity

Given that humans have about the same number of genes as mice and not so many more than worm, what makes us more complex? Antisense transcripts are implicated in many aspects of gene regulation. Is there a functional connection between antisense transcription and organismic complexity, that is, is antisense regulation especially prevalent in humans? We used the same robust protocol to identify antisense transcripts in humans and five other metazoan genomes (mouse, rat, chicken, fruit fly, and nematode), and found that the estimated proportions of genes involved in antisense transcription are highly sensitive to the number of transcripts included in the analysis. By controlling for transcript abundance, we find that the probability that any given transcript is putatively involved in sense-antisense regulation is no higher in humans than in other vertebrates but appears unusually high in flies and especially low in nematodes. Similarly, there is no evidence that the proportion of sense-antisense transcripts is especially higher in humans than other vertebrates in a given subset of transcript sequences such as mRNAs, coding sequences, conserved, or nonconserved transcripts. Although antisense transcription might be enriched in mammalian brains compared with nonbrain tissues, it is no more enriched in human brain than in mouse brain. Overall, therefore, while we see striking variation between multicellular animals in the abundance of antisense transcripts, there is no evidence for a link between antisense transcription and organismic complexity. More particularly, we see no evidence that humans are in any way unusual among the vertebrates in this regard. Instead, our results suggest that antisense transcription might be prevalent in almost all metazoan genomes, nematodes being an unexplained exception.

A systematic analysis of disease-associated variants in the 3' regulatory regions of human protein-coding genes I: general principles and overview

The 3' regulatory regions (3' RRs) of human genes play an important role in regulating mRNA 3' end formation, stability/degradation, nuclear export, subcellular localization and translation and are consequently rich in regulatory elements. Although 3' RRs contain only approximately 0.2% of known disease-associated mutations, this is likely to represent a rather conservative estimate of their actual prevalence. In an attempt to catalogue 3' RR-mediated disease and also to gain a greater understanding of the functional role of regulatory elements within 3' RRs, we have performed a systematic analysis of disease-associated 3' RR variants; 121 3' RR variants in 94 human genes were collated. These included 17 mutations in the upstream core polyadenylation signal sequence (UCPAS), 81 in the upstream sequence (USS) between the translational termination codon and the UCPAS, 6 in the left arm of the 'spacer' sequence (LAS) between the UCPAS and the pre-mRNA cleavage site (CS), 3 in the right arm of the 'spacer' sequence (RAS) or downstream core polyadenylation signal sequence (DCPAS) and 7 in the downstream sequence (DSS) of the 3'-flanking region, with 7 further mutations being treated as isolated examples. All the UCPAS mutations and the rather unusual cases of the DMPK, SCA8, FCMD and GLA mutations exert a significant effect on the mRNA phenotype and are usually associated with monogenic disease. By contrast, most of the remaining variants are polymorphisms that exert a comparatively minor influence on mRNA expression, but which may nevertheless predispose to or otherwise modify complex clinical phenotypes. Considerable efforts have been made to validate/elucidate the mechanisms through which the 3' untranslated region (3' UTR) variants affect gene expression. It is hoped that the integrative approach employed here in the study of naturally occurring variants of actual or potential pathological significance will serve to complement ongoing efforts to identify all functional regulatory elements in the human genome.

Complex Loci in human and mouse genomes

Mammalian genomes harbor a larger than expected number of complex loci, in which multiple genes are coupled by shared transcribed regions in antisense orientation and/or by bidirectional core promoters. To determine the incidence, functional significance, and evolutionary context of mammalian complex loci, we identified and characterized 5,248 cis-antisense pairs, 1,638 bidirectional promoters, and 1,153 chains of multiple cis-antisense and/or bidirectionally promoted pairs from 36,606 mouse transcriptional units (TUs), along with 6,141 cis-antisense pairs, 2,113 bidirectional promoters, and 1,480 chains from 42,887 human TUs. In both human and mouse, 25% of TUs resided in cis-antisense pairs, only 17% of which were conserved between the two organisms, indicating frequent species specificity of antisense gene arrangements. A sampling approach indicated that over 40% of all TUs might actually be in cis-antisense pairs, and that only a minority of these arrangements are likely to be conserved between human and mouse. Bidirectional promoters were characterized by variable transcriptional start sites and an identifiable midpoint at which overall sequence composition changed strand and the direction of transcriptional initiation switched. In microarray data covering a wide range of mouse tissues, genes in cis-antisense and bidirectionally promoted arrangement showed a higher probability of being coordinately expressed than random pairs of genes. In a case study on homeotic loci, we observed extensive transcription of nonconserved sequences on the noncoding strand, implying that the presence rather than the sequence of these transcripts is of functional importance. Complex loci are ubiquitous, host numerous nonconserved gene structures and lineage-specific exonification events, and may have a cis-regulatory impact on the member genes.