Expression profiling of antisense transcripts on DNA arrays

Natural antisense transcripts as versatile regulators of gene expression

Long non-coding RNAs (lncRNAs) are emerging as a major class of gene products that have central roles in cell and developmental biology. Natural antisense transcripts (NATs) are an important subset of lncRNAs that are expressed from the opposite strand of protein-coding and non-coding genes and are a genome-wide phenomenon in both eukaryotes and prokaryotes. In eukaryotes, a myriad of NATs participate in regulatory pathways that affect expression of their cognate sense genes. Recent developments in the study of NATs and lncRNAs and large-scale sequencing and bioinformatics projects suggest that whether NATs regulate expression, splicing, stability or translation of the sense transcript is influenced by the pattern and degrees of overlap between the sense-antisense pair. Moreover, epigenetic gene regulatory mechanisms prevail in somatic cells whereas mechanisms dependent on the formation of double-stranded RNA intermediates are prevalent in germ cells. The modulating effects of NATs on sense transcript expression make NATs rational targets for therapeutic interventions.

Non-coding antisense transcripts: fine regulation of gene expression in cancer

Natural antisense transcripts (NATs) are coding or non-coding RNA sequences transcribed on the opposite direction from the same genomic locus. NATs are widely distributed throughout the human genome and seem to play crucial roles in physiological and pathological processes, through newly described and targeted mechanisms. NATs represent the intricate complexity of the genome organization and constitute another layer of potential targets in disease. Here, we focus on the interesting and unique role of non-coding NATs in cancer, paying particular attention to those acting as miRNA sponges.

Interdependent Transcription of a Natural Sense/Antisense Transcripts Pair (SLC34A1/PFN3)

Natural antisense transcripts (NATs) constitute a significant group of regulatory, long noncoding RNAs. They are prominently expressed in testis but are also detectable in other organs. NATs are transcribed at low levels and co-expressed with related protein coding sense transcripts. Nowadays NATs are generally considered as regulatory, long noncoding RNAs without closer focus on the inevitable interference between sense and antisense expression. This work describes a cellular system where sense and antisense transcription of a specific locus (SLC34A1/PFN3) is induced using epigenetic modifiers and CRISPR-Cas9. The renal cell lines HEK293 and HKC-8 do not express SLC34A1/PFN3 under normal culture conditions. Five-day exposure to dexamethasone significantly stimulates sense transcript (SLC34A1) levels and antisense (PFN3) minimally; the effect is only seen in HEK293 cells. Enhanced expression is paralleled by reduced sense promoter methylation and an increase in activating histone marks. Expression is further modulated by cassettes that stimulate the expression of sense or antisense transcript but disrupt protein coding potential. Constitutive expression of a 5′-truncated SLC34A1 transcript increases sense expression independent of dexamethasone induction but also stimulates antisense expression. Concordant expression is confirmed with the antisense knock-in that also enhances sense expression. The antisense effect acts on transcription in cis since transient transfection with sense or antisense constructs fails to stimulate the expression of the opposite transcript. These results suggest that bi-directional transcription of the SLC34A1/PFN3 locus has a stimulatory influence on the expression of the opposite transcript involving epigenetic changes of the promoters. In perspective of extensive, previous research into bi-directionally transcribed SLC34A loci, the findings underpin a hypothesis where NATs display different biological roles in soma and germ cells. Accordingly, we propose that in somatic cells, NATs act like lncRNAs–with the benefit of close proximity to a potential target gene. In germ cells, however, recent evidence suggests different biological roles for NATs that require RNA complementarity and double-stranded RNA formation.

Widespread formation of double-stranded RNAs in testis

The testis transcriptome is highly complex and includes RNAs that potentially hybridize to form double-stranded RNA (dsRNA). We isolated dsRNA using the monoclonal J2 antibody and deep-sequenced the enriched samples from testes of juvenile Dicer1 knockout mice, age-matched controls, and adult animals. Comparison of our data set with recently published data from mouse liver revealed that the dsRNA transcriptome in testis is markedly different from liver: In testis, dsRNA-forming transcripts derive from mRNAs including promoters and immediate downstream regions, whereas in somatic cells they originate more often from introns and intergenic transcription. The genes that generate dsRNA are significantly expressed in isolated male germ cells with particular enrichment in pachytene spermatocytes. dsRNA formation is lower on the sex (X and Y) chromosomes. The dsRNA transcriptome is significantly less complex in juvenile mice as compared to adult controls and, possibly as a consequence, the knockout of Dicer1 has only a minor effect on the total number of transcript peaks associated with dsRNA. The comparison between dsRNA-associated genes in testis and liver with a reported set of genes that produce endogenous siRNAs reveals a significant overlap in testis but not in liver. Testis dsRNAs also significantly associate with natural antisense genes-again, this feature is not observed in liver. These findings point to a testis-specific mechanism involving natural antisense transcripts and the formation of dsRNAs that feed into the RNA interference pathway, possibly to mitigate the mutagenic impacts of recombination and transposon mobilization.

Small regulatory noncoding RNAs in Drosophila melanogaster: biogenesis and biological functions

RNA interference (RNAi) is an important phenomenon that has diverse genetic regulatory functions at the pre- and posttranscriptional levels. The major trigger for the RNAi pathway is double-stranded RNA (dsRNA). dsRNA is processed to generate various types of major small noncoding RNAs (ncRNAs) that include microRNAs (miRNAs), small interfering RNAs (siRNAs) and PIWI-interacting RNAs (piRNAs) in Drosophila melanogaster (D. melanogaster). Functionally, these small ncRNAs play critical roles in virtually all biological systems and developmental pathways. Identification and processing of dsRNAs and activation of RNAi machinery are the three major academic interests that surround RNAi research. Mechanistically, some of the important biological functions of RNAi are achieved through: (i) supporting genomic stability via degradation of foreign viral genomes; (ii) suppressing the movement of transposable elements and, most importantly, (iii) post-transcriptional regulation of gene expression by miRNAs that contribute to regulation of epigenetic modifications such as heterochromatin formation and genome imprinting. Here, we review various routes of small ncRNA biogenesis, as well as different RNAi-mediated pathways in D. melanogaster with a particular focus on signaling pathways. In addition, a critical discussion of the most relevant and latest findings that concern the significant contribution of small ncRNAs to the regulation of D. melanogaster physiology and pathophysiology is presented.

A cautionary tale of sense-antisense gene pairs: independent regulation despite inverse correlation of expression

Long non-coding RNAs (lncRNAs) have been proven to play important roles in diverse cellular processes including the DNA damage response. Nearly 40% of annotated lncRNAs are transcribed in antisense direction to other genes and have often been implicated in their regulation via transcript- or transcription-dependent mechanisms. However, it remains unclear whether inverse correlation of gene expression would generally point toward a regulatory interaction between the genes. Here, we profiled lncRNA and mRNA expression in lung and liver cancer cells after exposure to DNA damage. Our analysis revealed two pairs of mRNA-lncRNA sense-antisense transcripts being inversely expressed upon DNA damage. The lncRNA NOP14-AS1 was strongly upregulated upon DNA damage, while the mRNA for NOP14 was downregulated, both in a p53-dependent manner. For another pair, the lncRNA LIPE-AS1 was downregulated, while its antisense mRNA CEACAM1 was upregulated. To test whether as expected the antisense genes would regulate each other resulting in this highly significant inverse correlation, we employed antisense oligonucleotides and RNAi to study transcript-dependent effects as well as dCas9-based transcriptional modulation by CRISPRi/CRISPRa for transcription-dependent effects. Surprisingly, despite the strong stimulus-dependent inverse correlation, our data indicate that neither transcript- nor transcription-dependent mechanisms explain the inverse regulation of NOP14-AS1:NOP14 or LIPE-AS1:CEACAM1 expression. Hence, sense-antisense pairs whose expression is strongly—positively or negatively—correlated can be nonetheless regulated independently. This highlights the requirement of individual experimental studies for each antisense pair and prohibits drawing conclusions on regulatory mechanisms from expression correlations.

Natural Antisense Transcripts at the Interface between Host Genome and Mobile Genetic Elements

Non-coding RNAs are involved in epigenetic processes, playing a role in the regulation of gene expression at the transcriptional and post-transcriptional levels. A particular group of ncRNA are natural antisense transcripts (NATs); these are transcribed in the opposite direction to protein coding transcripts and are widespread in eukaryotes. Their abundance, evidence of phylogenetic conservation and an increasing number of well-characterized examples of antisense-mediated gene regulation are indicative of essential biological roles of NATs. There is evidence to suggest that they interfere with their corresponding sense transcript to elicit concordant and discordant regulation. The main mechanisms involved include transcriptional interference as well as dsRNA formation. Sense–antisense hybrid formation can trigger RNA interference, RNA editing or protein kinase R. However, the exact molecular mechanisms elicited by NATs in the context of these regulatory roles are currently poorly understood. Several examples confirm that ectopic expression of antisense transcripts trigger epigenetic silencing of the related sense transcript. Genomic approaches suggest that the antisense transcriptome carries a broader biological significance which goes beyond the physiological regulation of the directly related sense transcripts. Because NATs show evidence of conservation we speculate that they played a role in evolution, with early eukaryotes gaining selective advantage through the regulatory effects. With the surge of genome and transcriptome sequencing projects, there is promise of a more comprehensive understanding of the biological role of NATs and the regulatory mechanisms involved.

Ectopically expressed Slc34a2a sense-antisense transcripts cause a cerebellar phenotype in zebrafish embryos depending on RNA complementarity and Dicer

Natural antisense transcripts (NATs) are complementary to protein coding genes and potentially regulate their expression. Despite widespread occurrence of NATs in the genomes of higher eukaryotes, their biological role and mechanism of action is poorly understood. Zebrafish embryos offer a unique model system to study sense-antisense transcript interplay at whole organism level. Here, we investigate putative antisense transcript-mediated mechanisms by ectopically co-expressing the complementary transcripts during early zebrafish development. In zebrafish the gene Slc34a2a (Na-phosphate transporter) is bi-directionally transcribed, the NAT predominantly during early development up to 48 hours after fertilization. Declining levels of the NAT, Slc34a2a(as), coincide with an increase of the sense transcript. At that time, sense and antisense transcripts co-localize in the endoderm at near equal amounts. Ectopic expression of the sense transcript during embryogenesis leads to specific failure to develop a cerebellum. The defect is RNA-mediated and dependent on sense-antisense complementarity. Overexpression of a Slc34a2a paralogue (Slc34a2b) or the NAT itself had no phenotypic consequences. Knockdown of Dicer rescued the brain defect suggesting that RNA interference is required to mediate the phenotype. Our results corroborate previous reports of Slc34a2a-related endo-siRNAs in two days old zebrafish embryos and emphasize the importance of coordinated expression of sense-antisense transcripts. Our findings suggest that RNAi is involved in gene regulation by certain natural antisense RNAs.

Natural antisense transcription from a comparative perspective

Natural antisense transcripts (NATs) can interfere with the expression of complementary sense transcripts with exquisite specificity. We have previously cloned NATs of Slc34a loci (encoding Na-phosphate transporters) from fish and mouse. Here we report the cloning of a human SLC34A1-related NAT that represents an alternatively spliced PFN3 transcript (Profilin3). The transcript is predominantly expressed in testis. Phylogenetic comparison suggests two distinct mechanisms producing Slc34a-related NATs: Alternative splicing of a transcript from a protein coding downstream gene (Pfn3, human/mouse) and transcription from the bi-directional promoter (Rbpja, zebrafish). Expression analysis suggested independent regulation of the complementary Slc34a mRNAs. Analysis of randomly selected bi-directionally transcribed human/mouse loci revealed limited phylogenetic conservation and independent regulation of NATs. They were reduced on X chromosomes and clustered in regions that escape inactivation. Locus structure and expression pattern suggest a NATs-associated regulatory mechanisms in testis unrelated to the physiological role of the sense transcript encoded protein.

Endogenous siRNAs: regulators of internal affairs

Endo-siRNAs (endogenous small-interfering RNAs) have recently emerged as versatile regulators of gene expression. They derive from double-stranded intrinsic transcripts and are processed by Dicer and associate with Argonaute proteins. In Caenorhabditis elegans, endo-siRNAs are known as 22G and 26G RNAs and are involved in genome protection and gene regulation. Drosophila melanogaster endo-siRNAs are produced with the help of specific Dicer and Argonaute isoforms and play an essential role in transposon control and the protection from viral infections. Biological functions of endo-siRNAs in vertebrates include repression of transposable elements and chromatin organization, as well as gene regulation at the transcriptional and post-transcriptional levels.

Contribution of natural antisense transcription to an endogenous siRNA signature in human cells

BACKGROUND

Eukaryotic cells express a complex layer of noncoding RNAs. An intriguing family of regulatory RNAs includes transcripts from the opposite strand of protein coding genes, so called natural antisense transcripts (NATs). Here, we test the hypothesis that antisense transcription triggers RNA interference and gives rise to endogenous short RNAs (endo-siRNAs).

RESULTS

We used cloned human embryonic kidney cells (HEK293) followed by short RNAseq to investigate the small genic RNA transcriptome. 378 genes gave rise to short RNA reads that mapped to exons of RefSeq genes. The length profile of short RNAs showed a broad peak of 20-24 nucleotides, indicative of endo-siRNAs. Collapsed reads mapped predominantly to the first and the last exon of genes (74%). RNAs reads were intersected with sequences occupied by RNAPII or bound to Argonaute (AGO1 by crosslinking, ligation, and sequencing of hybrids, CLASH). In the first exon, 94% of the reads correlated with RNAPII occupancy with an average density of 130 (relative units); this decreased to 65%/20 in middle exons and 54%/12 in the last exon. CLASH reads mapping to multi-exon genes showed little distribution bias with an average of about 5 CLASH reads overlapping with 60% of the endo-siRNA reads. However, endo-siRNAs (21-25 nt) intersecting with CLASH reads were enriched at the 5'end and decreased towards the 3'end.We then investigated the 378 genes with particular focus on features indicative for short RNA production; however, found that endo-siRNA numbers did not correlate with gene structures that favor convergent transcription. In contrast, our gene set was found notably over-represented in the NATsDB sense/antisense group as compared to non-overlapping and non-bidirectional groups. Moreover, read counts showed no correlation with the steady-state levels of the related mRNAs and the pattern of endo-siRNAs proved reproducible after an induced mutagenic insult.

CONCLUSIONS

Our results suggest that antisense transcripts contribute to low levels of endo-siRNAs in fully differentiated human cells. A characteristic endo-siRNA footprint is being produced at sites of RNAPII transcription which is also related to AGO1. This endo-siRNA signature represents an intriguing finding and its reproducibility suggests that the production of endo-siRNAs is a regulated process with potential homoeostatic impact.

Non-coding RNAs: multi-tasking molecules in the cell

In the last years it has become increasingly clear that the mammalian transcriptome is highly complex and includes a large number of small non-coding RNAs (sncRNAs) and long noncoding RNAs (lncRNAs). Here we review the biogenesis pathways of the three classes of sncRNAs, namely short interfering RNAs (siRNAs), microRNAs (miRNAs) and PIWI-interacting RNAs (piRNAs). These ncRNAs have been extensively studied and are involved in pathways leading to specific gene silencing and the protection of genomes against virus and transposons, for example. Also, lncRNAs have emerged as pivotal molecules for the transcriptional and post-transcriptional regulation of gene expression which is supported by their tissue-specific expression patterns, subcellular distribution, and developmental regulation. Therefore, we also focus our attention on their role in differentiation and development. SncRNAs and lncRNAs play critical roles in defining DNA methylation patterns, as well as chromatin remodeling thus having a substantial effect in epigenetics. The identification of some overlaps in their biogenesis pathways and functional roles raises the hypothesis that these molecules play concerted functions in vivo, creating complex regulatory networks where cooperation with regulatory proteins is necessary. We also highlighted the implications of biogenesis and gene expression deregulation of sncRNAs and lncRNAs in human diseases like cancer.

The presence, role and clinical use of spermatozoal RNAs

BACKGROUND Spermatozoa are highly differentiated, transcriptionally inert cells characterized by a compact nucleus with minimal cytoplasm. Nevertheless they contain a suite of unique RNAs that are delivered to oocyte upon fertilization. They are likely integrated as part of many different processes including genome recognition, consolidation-confrontation, early embryonic development and epigenetic transgenerational inherence. Spermatozoal RNAs also provide a window into the developmental history of each sperm thereby providing biomarkers of fertility and pregnancy outcome which are being intensely studied. METHODS Literature searches were performed to review the majority of spermatozoal RNA studies that described potential functions and clinical applications with emphasis on Next-Generation Sequencing. Human, mouse, bovine and stallion were compared as their distribution and composition of spermatozoal RNAs, using these techniques, have been described. RESULTS Comparisons highlighted the complexity of the population of spermatozoal RNAs that comprises rRNA, mRNA and both large and small non-coding RNAs. RNA-seq analysis has revealed that only a fraction of the larger RNAs retain their structure. While rRNAs are the most abundant and are highly fragmented, ensuring a translationally quiescent state, other RNAs including some mRNAs retain their functional potential, thereby increasing the opportunity for regulatory interactions. Abundant small non-coding RNAs retained in spermatozoa include miRNAs and piRNAs. Some, like miR-34c are essential to the early embryo development required for the first cellular division. Others like the piRNAs are likely part of the genomic dance of confrontation and consolidation. Other non-coding spermatozoal RNAs include transposable elements, annotated lnc-RNAs, intronic retained elements, exonic elements, chromatin-associated RNAs, small-nuclear ILF3/NF30 associated RNAs, quiescent RNAs, mse-tRNAs and YRNAs. Some non-coding RNAs are known to act as epigenetic modifiers, inducing histone modifications and DNA methylation, perhaps playing a role in transgenerational epigenetic inherence. Transcript profiling holds considerable potential for the discovery of fertility biomarkers for both agriculture and human medicine. Comparing the differential RNA profiles of infertile and fertile individuals as well as assessing species similarities, should resolve the regulatory pathways contributing to male factor infertility. CONCLUSIONS Dad delivers a complex population of RNAs to the oocyte at fertilization that likely influences fertilization, embryo development, the phenotype of the offspring and possibly future generations. Development is continuing on the use of spermatozoal RNA profiles as phenotypic markers of male factor status for use as clinical diagnostics of the father's contribution to the birth of a healthy child.

Mechanisms of small RNA generation from cis-NATs in response to environmental and developmental cues

A large proportion of eukaryotic genomes is transcribed from both positive and negative strands of DNA and thus may generate overlapping sense and antisense transcripts. Some of these so-called natural antisense transcripts (NATs) are possibly co-regulated. When the overlapping sense and antisense transcripts are expressed at the same time in the same cell in response to various developmental and environmental cues; they may form double-stranded RNAs, which could be recognized by the small RNA biogenesis machinery and processed into small interfering RNAs (siRNAs). cis-NAT-derived siRNAs (nat-siRNAs) are present in plants, animals, and fungi. In plants, the presence of nat-siRNAs is supported not only by Northern blot and genetic analyses, but also by the fact that there is an overall sixfold enrichment of siRNAs in the overlapping regions of cis-NATs and 19%-29% of the siRNA-generating cis-NATs in plants give rise to siRNAs only in their overlapping regions. Silencing mediated by nat-siRNAs is one of the mechanisms for regulating the expression of the cis-NATs. This review focuses on challenging issues related to the biogenesis mechanisms as well as regulation and detection of nat-siRNAs. The advantages and limitations of new technologies for detecting cis-NATs, including direct RNA sequencing and strand-specific RNA sequencing, are also discussed.

High mortality of juvenile gilthead sea bream (Sparus aurata) from photobacteriosis is associated with alternative macrophage activation and anti-inflammatory response: results of gene expression profiling of early responses in the head kidney

The halophilic bacterium Photobacterium damselae subsp. piscicida (Phdp) represents a substantial health problem for several fish species in aquaculture. Bacteria that reside free and inside phagocytes cause acute and chronic forms of photobacteriosis. Infections of juveniles rapidly kill up to 90-100% fish. Factors underlying failure of the immune protection against bacteria remain largely unknown. The reported study used a transcriptomic approach to address this issue. Juvenile sea breams (0.5 g) were challenged by immersion in salt water containing 2.89 × 10(8) CFU of a virulent Phdp and the head kidney was sampled after 24- and 48-h. Analyses were performed using the second version of a 44 k oligonucleotide DNA microarray that represents 19,734 sea bream unique transcripts and covers diverse immune pathways. Expression changes of selected immune genes were validated with qPCR. Results suggested rapid recognition of the pathogen, as testified by up-regulation of lectins and antibacterial proteins (bactericidal permeability-increasing protein lectins, lysozyme, intracellular and extracellular proteases), chemokines and chemokine receptors. Increased expression of proteins involved in iron and heme metabolism also could be a response against bacteria that are dependent on iron. However, negative regulators of immune/inflammatory response were preponderant among the up-regulated genes. A remarkable finding was the increased expression of IL-10 in concert with up-regulation of arginase I and II and proteins of the polyamine biosynthesis pathway that diverts the arginine flux from the production of reactive nitrogen species. Such expression changes are characteristic for alternatively activated macrophages that do not develop acute inflammatory responses. Immune suppression can be induced by the host to reduce tissue damages or by the pathogen to evade host response.

Identification of potential antisense transcripts in rice using conventional microarray

Natural antisense transcripts (NATs) are endogenous transcripts that contain reverse complementary sequences to other RNAs (usually called sense transcripts). NATs regulate the expression of sense transcripts in a wide range of species. The identification and analysis of NATs are the prerequisite to elucidate their functions. Microarray is a genome-wide method to detect gene expression. However, conventional microarrays do not contain the specific probes of NATs; thus, they cannot be utilized to detect NATs. In this article, we developed a novel method to identify potential NATs with the conventional microarrays. In this method of our study, we labeled the first strand cDNA from one sample with Cy5 and labeled the second strand cDNA from another sample with Cy3, and then hybridized these labeled samples with oligonucleotide microarray. Using this method, we identified 920 potential NATs in rice variety Nipponbare. Among these potential NATs, 88 of them were confirmed by either full-length cDNA or orientated ESTs (expressed sequence tags). This is the first time that a conventional oligonucleotide microarray was employed to identify NATs in rice.

Identification of sense and antisense transcripts regulated by drought in sugarcane

Sugarcane is an important sugar and energy crop that can be used efficiently for biofuels production. The development of sugarcane cultivars tolerant to drought could allow for the expansion of plantations to sub-prime regions. Knowledge on the mechanisms underlying drought responses and its relationship with carbon partition would greatly help to define routes to increase yield. In this work we studied sugarcane responses to drought using a custom designed oligonucleotide array with 21,901 different probes. The oligoarrays were designed to contain probes that detect transcription in both sense and antisense orientation. We validated the results obtained using quantitative real-time PCR (qPCR). A total of 987 genes were differentially expressed in at least one sample of sugarcane plants submitted to drought for 24, 72 and 120 h. Among them, 928 were sense transcripts and 59 were antisense transcripts. Genes related to Carbohydrate Metabolism, RNA Metabolism and Signal Transduction were selected for gene expression validation by qPCR that indicated a validation percentage of 90%. From the probes presented on the array, 75% of the sense probes and 11.9% of the antisense probes have signal above background and can be classified as expressed sequences. Our custom sugarcane oligonucleotide array provides sensitivity and good coverage of sugarcane transcripts for the identification of a representative proportion of natural antisense transcripts (NATs) and sense-antisense transcript pairs (SATs). The antisense transcriptome showed, in most cases, co-expression with respective sense transcripts.

Hormone-dependent expression of a steroidogenic acute regulatory protein natural antisense transcript in MA-10 mouse tumor Leydig cells

Cholesterol transport is essential for many physiological processes, including steroidogenesis. In steroidogenic cells hormone-induced cholesterol transport is controlled by a protein complex that includes steroidogenic acute regulatory protein (StAR). Star is expressed as 3.5-, 2.8-, and 1.6-kb transcripts that differ only in their 3′-untranslated regions. Because these transcripts share the same promoter, mRNA stability may be involved in their differential regulation and expression. Recently, the identification of natural antisense transcripts (NATs) has added another level of regulation to eukaryotic gene expression. Here we identified a new NAT that is complementary to the spliced Star mRNA sequence. Using 5′ and 3′ RACE, strand-specific RT-PCR, and ribonuclease protection assays, we demonstrated that Star NAT is expressed in MA-10 Leydig cells and steroidogenic murine tissues. Furthermore, we established that human chorionic gonadotropin stimulates Star NAT expression via cAMP. Our results show that sense-antisense Star RNAs may be coordinately regulated since they are co-expressed in MA-10 cells. Overexpression of Star NAT had a differential effect on the expression of the different Star sense transcripts following cAMP stimulation. Meanwhile, the levels of StAR protein and progesterone production were downregulated in the presence of Star NAT. Our data identify antisense transcription as an additional mechanism involved in the regulation of steroid biosynthesis.

What are natural antisense transcripts good for?

NATs (natural antisense transcripts) are important regulators of eukaryotic gene expression. Interference between the expression of protein-coding sense transcripts and the corresponding NAT is well documented. In the present review, we focus on an additional, higher-order role of NATs that is currently emerging. The recent discovery of endogenous siRNAs (short interfering RNAs), as well as NAT-induced transcriptional gene silencing, are key to the proposed novel function of NATs.

Development of an oligo DNA microarray for the European sea bass and its application to expression profiling of jaw deformity

Background

The European sea bass (Dicentrarchus labrax) is a marine fish of great importance for fisheries and aquaculture. Functional genomics offers the possibility to discover the molecular mechanisms underlying productive traits in farmed fish, and a step towards the application of marker assisted selection methods in this species. To this end, we report here on the development of an oligo DNA microarray for D. labrax.

Results

A database consisting of 19,048 unique transcripts was constructed, of which 12,008 (63%) could be annotated by similarity and 4,692 received a GO functional annotation. Two non-overlapping 60mer probes were designed for each unique transcript and in-situ synthesized on glass slides using Agilent SurePrint™ technology. Probe design was positively completed for 19,035 target clusters; the oligo microarray was then applied to profile gene expression in mandibles and whole-heads of fish affected by prognathism, a skeletal malformation that strongly affects sea bass production. Statistical analysis identified 242 transcripts that are significantly down-regulated in deformed individuals compared to normal fish, with a significant enrichment in genes related to nervous system development and functioning. A set of genes spanning a wide dynamic range in gene expression level were selected for quantitative RT-PCR validation. Fold change correlation between microarray and qPCR data was always significant.

Conclusions

The microarray platform developed for the European sea bass has a high level of flexibility, reliability, and reproducibility. Despite the well known limitations in achieving a proper functional annotation in non-model species, sufficient information was obtained to identify biological processes that are significantly enriched among differentially expressed genes. New insights were obtained on putative mechanisms involved on mandibular prognathism, suggesting that bone/nervous system development might play a role in this phenomenon.

Integrative analysis of the human cis-antisense gene pairs, miRNAs and their transcription regulation patterns

Cis-antisense gene pairs (CASGPs) can transcribe mRNAs from an opposite strand of a given locus. To classify and understand diverse CASGP phenomena in the human we compiled a genome-wide catalog of CASGPs and integrated these sequences with microarray, SAGE and miRNA data. Using the concept of overlapping regions and clustering of SA transcripts by chromosome coordinates, we identified up to 9000 overlapping antisense loci. Four thousand three hundred and seventy-four of these CASGPs form 1759 complex gene architectures. We found that ∼35% (6347/18160) of RefSeq genes are overlapped with the antisense transcripts. About 30% of Affymetrix U133 microarray initial sequences map transcripts of ∼35% CASGPs and reveal mostly concordant expression in CASGPs. We found strong significant overrepresentation of human miRNA genes in loci of CASGPs. We developed a data-driven model of cross-talk between co-expressed CASGPs and DICER1-mediated miRNA pathway in normal spermatogenesis and in severe teratozoospermia. Specifically, we revealed complex SA structural–functional gene module composing the protein-coding genes, WDR6, DALRD3, NDUFAF3 and ncRNA precursors, mir-425 and mir-191, which could provide downregulation of ncRNA pathway via direct targeting DICER1 and basonuclin 2 transcripts by mir-425 and mir-191 in normal spermatogenesis, but this mechanism is switched off in severe teratozoospermia. The database is available from http://globalisland.bii.a-star.edu.sg/∼jiangtao/sas/index3.php?link =about

Experimental annotation of channel catfish virus by probabilistic proteogenomic mapping

Experimental identification of expressed proteins by proteomics constitutes the most reliable approach to identify genomic location and structure of protein-coding genes and substantially complements computational genome annotation. Channel catfish herpesvirus (CCV) is a simple comparative model for understanding herpesvirus biology and the evolution of the Herpesviridae. The canonical CCV genome has 76 predicted ORF and only 12 of these have been confirmed experimentally. We describe a modification of a statistical method, which assigns significance measures, q-values, to peptide identifications based on 2-D LC ESI MS/MS, real-decoy database searches and SEQUEST XCorr and DeltaC(n) scores. We used this approach to identify CCV proteins expressed during its replication in cell culture, to determine protein composition of mature virions and, consequently, to refine the canonical CCV genome annotation. To complement trypsin, we used partial proteinase K digestion, which yielded greater proteome coverage. At FDR <5%, for peptide identifications, we identified 25/76 previously predicted ORF using trypsin and 31/76 using proteinase K. Furthermore, we identified 17 novel protein-coding regions (7 potential ATG-initiated ORF). Most of these novel ORF encode small proteins (<100 amino acids). Directed, strand-specific reverse transcription real-time PCR confirmed RNA expression from 6/7 novel ATG-initiated ORF investigated.

Identification of novel endogenous antisense transcripts by DNA microarray analysis targeting complementary strand of annotated genes

Background

Recent transcriptomic analyses in mammals have uncovered the widespread occurrence of endogenous antisense transcripts, termed natural antisense transcripts (NATs). NATs are transcribed from the opposite strand of the gene locus and are thought to control sense gene expression, but the mechanism of such regulation is as yet unknown. Although several thousand potential sense-antisense pairs have been identified in mammals, examples of functionally characterized NATs remain limited. To identify NAT candidates suitable for further functional analyses, we performed DNA microarray-based NAT screening using mouse adult normal tissues and mammary tumors to target not only the sense orientation but also the complementary strand of the annotated genes.

Results

First, we designed microarray probes to target the complementary strand of genes for which an antisense counterpart had been identified only in human public cDNA sources, but not in the mouse. We observed a prominent expression signal from 66.1% of 635 target genes, and 58 genes of these showed tissue-specific expression. Expression analyses of selected examples (Acaa1b and Aard) confirmed their dynamic transcription in vivo. Although interspecies conservation of NAT expression was previously investigated by the presence of cDNA sources in both species, our results suggest that there are more examples of human-mouse conserved NATs that could not be identified by cDNA sources. We also designed probes to target the complementary strand of well-characterized genes, including oncogenes, and compared the expression of these genes between mammary cancerous tissues and non-pathological tissues. We found that antisense expression of 95 genes of 404 well-annotated genes was markedly altered in tumor tissue compared with that in normal tissue and that 19 of these genes also exhibited changes in sense gene expression. These results highlight the importance of NAT expression in the regulation of cellular events and in pathological conditions.

Conclusion

Our microarray platform targeting the complementary strand of annotated genes successfully identified novel NATs that could not be identified by publically available cDNA data, and as such could not be detected by the usual "sense-targeting" microarray approach. Differentially expressed NATs monitored by this platform may provide candidates for investigations of gene function. An advantage of our microarray platform is that it can be applied to any genes and target samples of interest.

Natural antisense transcript of natriuretic peptide precursor A (NPPA): structural organization and modulation of NPPA expression

BACKGROUND

Mammalian transcriptome contains a large proportion of diverse and structurally complex noncoding RNAs. One class of such RNAs, natural antisense transcripts (NATs), are derived from the opposite strand of many protein-coding genes. Although the exact structure and functional relevance of most of the NATs is unknown, their emerging role as gene expression regulators raises the hypothesis that NATs might contribute to development of complex human disorders. The goal of our study was to investigate the involvement of NATs in regulation of candidate genes for blood pressure.

RESULTS

First we analysed blood pressure candidate genes for the presence of natural antisense transcripts. In silico analysis revealed that seven genes (ADD3, NPPA, ATP1A1, NPR2, CYP17A1, ACSM3, SLC14A2) have an antisense partner transcribed from the opposite strand. We characterized NPPA and its antisense transcript (NPPA-AS) in more detail. We found that NPPA-AS is expressed in a number of human tissues as a collection of alternatively spliced isoforms and that NPPA-AS and NPPA can form RNA duplexes in vivo. We also demonstrated that a specific NPPA-AS isoform is capable of down-regulating the intron-retained NPPA mRNA variant. We studied the evolutionary conservation of NPPA-AS and were able to detect the presence of Nppa-as transcript in mouse.

CONCLUSION

Our results demonstrate functional interaction of NPPA-AS with NPPA at the RNA level and suggest that antisense transcription might be an important post-transcriptional mechanism modulating NPPA expression.

Identification of differentially expressed sense and antisense transcript pairs in breast epithelial tissues

BACKGROUND

More than 20% of human transcripts have naturally occurring antisense products (or natural antisense transcripts--NATs), some of which may play a key role in a range of human diseases. To date, several databases of in silico defined human sense-antisense (SAS) pairs have appeared, however no study has focused on differential expression of SAS pairs in breast tissue. We therefore investigated the expression levels of sense and antisense transcripts in normal and malignant human breast epithelia using the Affymetrix HG-U133 Plus 2.0 and Almac Diagnostics Breast Cancer DSA microarray technologies as well as massively parallel signature sequencing (MPSS) data.

RESULTS

The expression of more than 2500 antisense transcripts were detected in normal breast duct luminal cells and in primary breast tumors substantially enriched for their epithelial cell content by DSA microarray. Expression of 431 NATs were confirmed by either of the other two technologies. A corresponding sense transcript could be identified on DSA for 257 antisense transcripts. Of these SAS pairs, 163 have not been previously reported. A positive correlation of differential expression between normal and malignant breast samples was observed for most SAS pairs. Orientation specific RT-QPCR of selected SAS pairs validated their expression in several breast cancer cell lines and solid breast tumours.

CONCLUSION

Disease-focused and antisense enriched microarray platforms (such as Breast Cancer DSA) confirm the assumption that antisense transcription in the human breast is more prevalent than previously anticipated. Expression of a proportion of these NATs has already been confirmed by other technologies while the true existence of the remaining ones has to be validated. Nevertheless, future studies will reveal whether the relative abundances of antisense and sense transcripts have regulatory influences on the translation of these mRNAs.

Large-scale analysis of antisense transcription in wheat using the Affymetrix GeneChip Wheat Genome Array

Background

Natural antisense transcripts (NATs) are transcripts of the opposite DNA strand to the sense-strand either at the same locus (cis-encoded) or a different locus (trans-encoded). They can affect gene expression at multiple stages including transcription, RNA processing and transport, and translation. NATs give rise to sense-antisense transcript pairs and the number of these identified has escalated greatly with the availability of DNA sequencing resources and public databases. Traditionally, NATs were identified by the alignment of full-length cDNAs or expressed sequence tags to genome sequences, but an alternative method for large-scale detection of sense-antisense transcript pairs involves the use of microarrays. In this study we developed a novel protocol to assay sense- and antisense-strand transcription on the 55 K Affymetrix GeneChip Wheat Genome Array, which is a 3' in vitro transcription (3'IVT) expression array. We selected five different tissue types for assay to enable maximum discovery, and used the 'Chinese Spring' wheat genotype because most of the wheat GeneChip probe sequences were based on its genomic sequence. This study is the first report of using a 3'IVT expression array to discover the expression of natural sense-antisense transcript pairs, and may be considered as proof-of-concept.

Results

By using alternative target preparation schemes, both the sense- and antisense-strand derived transcripts were labeled and hybridized to the Wheat GeneChip. Quality assurance verified that successful hybridization did occur in the antisense-strand assay. A stringent threshold for positive hybridization was applied, which resulted in the identification of 110 sense-antisense transcript pairs, as well as 80 potentially antisense-specific transcripts. Strand-specific RT-PCR validated the microarray observations, and showed that antisense transcription is likely to be tissue specific. For the annotated sense-antisense transcript pairs, analysis of the gene ontology terms showed a significant over-representation of transcripts involved in energy production. These included several representations of ATP synthase, photosystem proteins and RUBISCO, which indicated that photosynthesis is likely to be regulated by antisense transcripts.

Conclusion

This study demonstrated the novel use of an adapted labeling protocol and a 3'IVT GeneChip array for large-scale identification of antisense transcription in wheat. The results show that antisense transcription is relatively abundant in wheat, and may affect the expression of valuable agronomic phenotypes. Future work should select potentially interesting transcript pairs for further functional characterization to determine biological activity.

Strand selective generation of endo-siRNAs from the Na/phosphate transporter gene Slc34a1 in murine tissues

Natural antisense transcripts (NATs) are important regulators of gene expression. Recently, a link between antisense transcription and the formation of endo-siRNAs has emerged. We investigated the bi-directionally transcribed Na/phosphate cotransporter gene (Slc34a1) under the aspect of endo-siRNA processing. Mouse Slc34a1 produces an antisense transcript that represents an alternative splice product of the Pfn3 gene located downstream of Slc34a1. The antisense transcript is prominently found in testis and in kidney. Co-expression of in vitro synthesized sense/antisense transcripts in Xenopus oocytes indicated processing of the overlapping transcripts into endo-siRNAs in the nucleus. Truncation experiments revealed that an overlap of at least 29 base-pairs is required to induce processing. We detected endo-siRNAs in mouse tissues that co express Slc34a1 sense/antisense transcripts by northern blotting. The orientation of endo-siRNAs was tissue specific in mouse kidney and testis. In kidney where the Na/phosphate cotransporter fulfils its physiological function endo-siRNAs complementary to the NAT were detected, in testis both orientations were found. Considering the wide spread expression of NATs and the gene silencing potential of endo-siRNAs we hypothesized a genome-wide link between antisense transcription and monoallelic expression. Significant correlation between random imprinting and antisense transcription could indeed be established. Our findings suggest a novel, more general role for NATs in gene regulation.

Identifying protein-coding genes in genomic sequences

A review of the main computational pipelines used to generate the human reference protein-coding gene sets.

The vast majority of the biology of a newly sequenced genome is inferred from the set of encoded proteins. Predicting this set is therefore invariably the first step after the completion of the genome DNA sequence. Here we review the main computational pipelines used to generate the human reference protein-coding gene sets.

Deep sequencing-based expression analysis shows major advances in robustness, resolution and inter-lab portability over five microarray platforms

The hippocampal expression profiles of wild-type mice and mice transgenic for δC-doublecortin-like kinase were compared with Solexa/Illumina deep sequencing technology and five different microarray platforms. With Illumina's digital gene expression assay, we obtained ∼2.4 million sequence tags per sample, their abundance spanning four orders of magnitude. Results were highly reproducible, even across laboratories. With a dedicated Bayesian model, we found differential expression of 3179 transcripts with an estimated false-discovery rate of 8.5%. This is a much higher figure than found for microarrays. The overlap in differentially expressed transcripts found with deep sequencing and microarrays was most significant for Affymetrix. The changes in expression observed by deep sequencing were larger than observed by microarrays or quantitative PCR. Relevant processes such as calmodulin-dependent protein kinase activity and vesicle transport along microtubules were found affected by deep sequencing but not by microarrays. While undetectable by microarrays, antisense transcription was found for 51% of all genes and alternative polyadenylation for 47%. We conclude that deep sequencing provides a major advance in robustness, comparability and richness of expression profiling data and is expected to boost collaborative, comparative and integrative genomics studies.

Arabidopsis transcriptome analysis under drought, cold, high-salinity and ABA treatment conditions using a tiling array

Plants respond and adapt to drought, cold and high-salinity stresses in order to survive. In this study, we applied Arabidopsis Affymetrix tiling arrays to study the whole genome transcriptome under drought, cold, high-salinity and ABA treatment conditions. The bioinformatic analysis using the tiling array data showed that 7,719 non-AGI transcriptional units (TUs) exist in the unannotated "intergenic" regions of Arabidopsis genome. These include 1,275 and 181 TUs that are induced and downregulated, respectively, by the stress or ABA treatments. Most of the non-AGI TUs are hypothetical non-protein-coding RNAs. About 80% of the non-AGI TUs belong to pairs of the fully overlapping sense-antisense transcripts (fSATs). Significant linear correlation between the expression ratios (treated/untreated) of the sense TUs and the ratios of the antisense TUs was observed in the SATs of AGI code/non-AGI TU. We studied the biogenesis mechanisms of the stress- or ABA-inducible antisense RNAs and found that the expression of sense TUs is necessary for the stress- or ABA-inducible expression of the antisense TUs in the fSATs (AGI code/non-AGI TU).

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.

Genome-wide analysis of antisense transcription with Affymetrix exon array

Background

A large number of natural antisense transcripts have been identified in human and mouse genomes. Study of their potential functions clearly requires cost-efficient method for expression analysis.

Results

Here we show that Affymetrix Exon arrays, which were designed to detect conventional transcripts in the sense orientation, can be used to monitor antisense expression across all exonic loci in mammalian genomes. Through modification of the cDNA synthesis protocol, we labeled single-strand cDNA in the reverse orientation as in the standard protocol, thus enabling the detection of antisense transcripts using the same array. Applying this technique to human Jurkat cells, we identified antisense transcription at 2,088 exonic loci of 1,516 UniGene clusters. Many of these antisense transcripts were not observed previously and some were validated by orientation-specific RT-PCR.

Conclusion

Our results suggest that with a modified protocol Affymetrix human, mouse and rat Exon arrays can be used as a routine method for genome-wide analysis of antisense transcription in these genomes.

Strand bias structure in mouse DNA gives a glimpse of how chromatin structure affects gene expression

Background

On a single strand of genomic DNA the number of As is usually about equal to the number of Ts (and similarly for Gs and Cs), but deviations have been noted for transcribed regions and origins of replication.

Results

The mouse genome is shown to have a segmented structure defined by strand bias. Transcription is known to cause a strand bias and numerous analyses are presented to show that the strand bias in question is not caused by transcription. However, these strand bias segments influence the position of genes and their unspliced length. The position of genes within the strand bias structure affects the probability that a gene is switched on and its expression level. Transcription has a highly directional flow within this structure and the peak volume of transcription is around 20 kb from the A-rich/T-rich segment boundary on the T-rich side, directed away from the boundary. The A-rich/T-rich boundaries are SATB1 binding regions, whereas the T-rich/A-rich boundary regions are not.

Conclusion

The direct cause of the strand bias structure may be DNA replication. The strand bias segments represent a further biological feature, the chromatin structure, which in turn influences the ease of transcription.

Uncovering information on expression of natural antisense transcripts in Affymetrix MOE430 datasets

Background

The function and significance of the widespread expression of natural antisense transcripts (NATs) is largely unknown. The ability to quantitatively assess changes in NAT expression for many different transcripts in multiple samples would facilitate our understanding of this relatively new class of RNA molecules.

Results

Here, we demonstrate that standard expression analysis Affymetrix MOE430 and HG-U133 GeneChips contain hundreds of probe sets that detect NATs. Probe sets carrying a "Negative Strand Matching Probes" annotation in NetAffx were validated using Ensembl by manual and automated approaches. More than 50 % of the 1,113 probe sets with "Negative Strand Matching Probes" on the MOE430 2.0 GeneChip were confirmed as detecting NATs. Expression of selected antisense transcripts as indicated by Affymetrix data was confirmed using strand-specific RT-PCR. Thus, Affymetrix datasets can be mined to reveal information about the regulated expression of a considerable number of NATs. In a correlation analysis of 179 sense-antisense (SAS) probe set pairs using publicly available data from 1637 MOE430 2.0 GeneChips a significant number of SAS transcript pairs were found to be positively correlated.

Conclusion

Standard expression analysis Affymetrix GeneChips can be used to measure many different NATs. The large amount of samples deposited in microarray databases represents a valuable resource for a quantitative analysis of NAT expression and regulation in different cells, tissues and biological conditions.