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

HIV-1 Natural Antisense Transcription and Its Role in Viral Persistence

Natural antisense transcripts (NATs) represent a class of RNA molecules that are transcribed from the opposite strand of a protein-coding gene, and that have the ability to regulate the expression of their cognate protein-coding gene via multiple mechanisms. NATs have been described in many prokaryotic and eukaryotic systems, as well as in the viruses that infect them. The human immunodeficiency virus (HIV-1) is no exception, and produces one or more NAT from a promoter within the 3’ long terminal repeat. HIV-1 antisense transcripts have been the focus of several studies spanning over 30 years. However, a complete appreciation of the role that these transcripts play in the virus lifecycle is still lacking. In this review, we cover the current knowledge about HIV-1 NATs, discuss some of the questions that are still open and identify possible areas of future research.

TSPAN31 suppresses cell proliferation in human cervical cancer through down-regulation of its antisense pairing with CDK4

Natural antisense transcripts (NAT) are prevalent phenomena in the mammalian genome and play significant regulatory roles in gene expression. While new insights into NAT continue to be revealed, their exact function and their underlying mechanisms in human cancer remain largely unclear. We identified a NAT of CDK4, referred to TSPAN31, which inhibits CDK4 mRNA and protein expression in human cervical cancer by targeting the 3'-untranslated region (3'-UTR) of the CDK4 mRNA. Furthermore, silencing the expression of the TSPAN31 mRNA rescued the TSPAN31 3'-UTR- or the TSPAN31 full-length-induced decrease in CDK4 expression. Noteworthy, we discovered that TSPAN31, as a member of the tetraspanin family, suppressed cell proliferation by down-regulating its antisense pairing with CDK4 and decreasing retinoblastoma protein phosphorylation in human cervical cancer. Therefore, the results of the present study suggest that TSPAN31 may serve as a potential molecular target for the development of novel anti-cancer agents. SIGNIFICANCE OF THE STUDY: Natural antisense transcripts are widely found in the genome and play an important role in the growth and development of cells. TSPAN31 is natural antisense transcript, and CDK4 is an important gene in the regulation of the cell cycle. Therefore, TSPAN31 and CDK4 have great significance in the study of tumour therapeutic targets.

Landscape of Overlapping Gene Expression in the Equine Placenta

Increasing evidence suggests that overlapping genes are much more common in eukaryotic genomes than previously thought. These different-strand overlapping genes are potential sense–antisense (SAS) pairs, which might have regulatory effects on each other. In the present study, we identified the SAS loci in the equine genome using previously generated stranded, paired-end RNA sequencing data from the equine chorioallantois. We identified a total of 1261 overlapping loci. The ratio of the number of overlapping regions to chromosomal length was numerically higher on chromosome 11 followed by chromosomes 13 and 12. These results show that overlapping transcription is distributed throughout the equine genome, but that distributions differ for each chromosome. Next, we evaluated the expression patterns of SAS pairs during the course of gestation. The sense and antisense genes showed an overall positive correlation between the sense and antisense pairs. We further provide a list of SAS pairs with both positive and negative correlation in their expression patterns throughout gestation. This study characterizes the landscape of sense and antisense gene expression in the placenta for the first time and provides a resource that will enable researchers to elucidate the mechanisms of sense/antisense regulation during pregnancy.

Mechanisms of Antisense Transcription Initiation with Implications in Gene Expression, Genomic Integrity and Disease Pathogenesis

Non-coding antisense transcripts arise from the strand opposite the sense strand. Over 70% of the human genome generates non-coding antisense transcripts while less than 2% of the genome codes for proteins. Antisense transcripts and/or the act of antisense transcription regulate gene expression and genome integrity by interfering with sense transcription and modulating histone modifications or DNA methylation. Hence, they have significant pathological and physiological relevance. Indeed, antisense transcripts were found to be associated with various diseases including cancer, diabetes, cardiac and neurodegenerative disorders, and, thus, have promising potentials for prognostic and diagnostic markers and therapeutic development. However, it is not clearly understood how antisense transcription is initiated and epigenetically regulated. Such knowledge would provide new insights into the regulation of antisense transcription, and hence disease pathogenesis with therapeutic development. The recent studies on antisense transcription initiation and its epigenetic regulation, which are limited, are discussed here. Furthermore, we concisely describe how antisense transcription/transcripts regulate gene expression and genome integrity with implications in disease pathogenesis and therapeutic development.

The coding and noncoding transcriptome of Neurospora crassa

Background

Long non protein coding RNAs (lncRNAs) have been identified in many different organisms and cell types. Emerging examples emphasize the biological importance of these RNA species but their regulation and functions remain poorly understood. In the filamentous fungus Neurospora crassa, the annotation and characterization of lncRNAs is incomplete.

Results

We have performed a comprehensive transcriptome analysis of Neurospora crassa by using ChIP-seq, RNA-seq and polysome fractionation datasets. We have annotated and characterized 1478 long intergenic noncoding RNAs (lincRNAs) and 1056 natural antisense transcripts, indicating that 20% of the RNA Polymerase II transcripts of Neurospora are not coding for protein. Both classes of lncRNAs accumulate at lower levels than protein-coding mRNAs and they are considerably shorter. Our analysis showed that the vast majority of lincRNAs and antisense transcripts do not contain introns and carry less H3K4me2 modifications than similarly expressed protein coding genes. In contrast, H3K27me3 modifications inversely correlate with transcription of protein coding and lincRNA genes. We show furthermore most lincRNA sequences evolve rapidly, even between phylogenetically close species.

Conclusions

Our transcriptome analyses revealed distinct features of Neurospora lincRNAs and antisense transcripts in comparison to mRNAs and showed that the prevalence of noncoding transcripts in this organism is higher than previously anticipated. The study provides a broad repertoire and a resource for further studies of lncRNAs.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4360-8) contains supplementary material, which is available to authorized users.

Decreased expression of the long non-coding RNA SLC7A11-AS1 predicts poor prognosis and promotes tumor growth in gastric cancer

Many lncRNA and mRNA sense-antisense transcripts have been systematically identified in malignant cells. However, the molecular mechanisms of most lncRNA-mRNA pairs in gastric cancer remain largely unknown. We found the gastric cancer-associated lncRNA SLC7A11-AS1 and coding transcript mRNA SLC7A11 in human gastric cancer specimens by microarray. SLC7A11-AS1, antisense to SLC7A11, is significantly down-regulated in gastric cancer and could promote tumor growth in vitro and in vivo. The effects of SLC7A11-AS1 depend on the regulation of SLC7A11 via the ASK1-p38^MAPK/JNK signaling pathway. These findings suggest that decreased expression of SLC7A11-AS1 contributes to the progression of gastric cancer and may be a novel diagnostic biomarker and effective therapeutic target in gastric cancer patients.

Mapping of leptin and its syntenic genes to chicken chromosome 1p

Background

Misidentification of the chicken leptin gene has hampered research of leptin signaling in this species for almost two decades. Recently, the genuine leptin gene with a GC-rich (~70%) repetitive-sequence content was identified in the chicken genome but without indicating its genomic position. This suggests that such GC-rich sequences are difficult to sequence and therefore substantial regions are missing from the current chicken genome assembly.

Results

A radiation hybrid panel of chicken-hamster Wg3hCl2 cells was used to map the genome location of the chicken leptin gene. Contrary to our expectations, based on comparative genome mapping and sequence characteristics, the chicken leptin was not located on a microchromosome, which are known to contain GC-rich and repetitive regions, but at the distal tip of the largest chromosome (1p). Following conserved synteny with other vertebrates, we also mapped five additional genes to this genomic region (ARF5, SND1, LRRC4, RBM28, and FLNC), bridging the genomic gap in the current Galgal5 build for this chromosome region. All of the short scaffolds containing these genes were found to consist of GC-rich (54 to 65%) sequences comparing to the average GC-content of 40% on chromosome 1. In this syntenic group, the RNA-binding protein 28 (RBM28) was in closest proximity to leptin. We deduced the full-length of the RBM28 cDNA sequence and profiled its expression patterns detecting a negative correlation (R = − 0.7) between the expression of leptin and of RBM28 across tissues that expressed at least one of the genes above the average level. This observation suggested a local regulatory interaction between these genes. In adipose tissues, we observed a significant increase in RBM28 mRNA expression in breeds with lean phenotypes.

Conclusion

Mapping chicken leptin together with a cluster of five syntenic genes provided the final proof for its identification as the true chicken ortholog. The high GC-content observed for the chicken leptin syntenic group suggests that other similar clusters of genes in GC-rich genomic regions are missing from the current genome assembly (Galgal5), which should be resolved in future assemblies of the chicken genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0543-1) contains supplementary material, which is available to authorized users.

MiR-377 reverses cancerous phenotypes of pancreatic cells via suppressing DNMT1 and demethylating tumor suppressor genes

AIM

The aim was to investigate the effect of miR-377 on DNMT1 expression and cancer phenotype in pancreatic cancer cells.

MATERIALS & METHODS

Real-time PCR, luciferase assay, MTT and Annexin-PI staining were used.

RESULTS

Decreased miR-377 and increased DNMT1 (verified as a target for mir-377) levels in pancreatic cancer tissues and cell lines in comparison with normal tissues was confirmed to be influenced by promoter methylation. Also hypermethylation of BNIP3, SPARC, TFPI2 and PENK promoters was observed in tumor samples but not in normal tissues which negatively correlated with their expression. Restoration of miR-377 resulted in a reduction of the expression of DNMT1 and reactivation of BNIP3 and SPARC genes via promoter demethylation. Furthermore, enhanced expression of miR-377 could significantly inhibit cell proliferation and induce apoptosis.

CONCLUSION

Our findings showed that miR-377 through targeting DNMT1 could reduce DNA methylation of some tumor suppressor genes and restore their expression in pancreatic cancer cells.

Genome-wide Identification and Characterization of Natural Antisense Transcripts by Strand-specific RNA Sequencing in Ganoderma lucidum

Ganoderma lucidum is a white-rot fungus best-known for its medicinal and ligninolytic activities. To discover the underlying genes responsible for these activities, we identified and characterized the natural antisense transcripts (NATs) using strand-specific (ss) RNA-seq data obtained from the mycelia, primordia and fruiting bodies. NATs were identified using a custom pipeline and then subjected to functional enrichment and differential expression analyses. A total of 1613 cis- and 244 trans- sense and antisense transcripts were identified. Mapping to GO terms and KEGG pathways revealed that NATs were frequently associated with genes of particular functional categories in particular stages. ssRT-qPCR experiments showed that the expression profiles of 30 of 50 (60%) transcripts were highly correlated with those of the RNA-seq results (r ≥ 0.9). Expression profiles of 22 of 25 (88%) pairs of NATs and STs were highly correlated (p ≤ 0.01), with 15 having r ≥ 0.8 and 4 having r ≤ -0.8. Six lignin-modifying genes and their NATs were analyzed in detail. Diverse patterns of differential expression among different stages and positive and negative correlations were observed. These results suggested that NATs were implicated in gene expression regulation in a function-group and developmental-stage specific manner through complex mechanisms.

Natural antisense transcripts are linked to the modulation of mitochondrial function and teliospore dormancy in Ustilago maydis

The basidiomycete smut fungus Ustilago maydis causes common smut of corn. This disease is spread through the production of teliospores, which are thick-walled dormant structures characterized by low rates of respiration and metabolism. Teliospores are formed when the fungus grows within the plant, and the morphological steps involved in their formation have been described, but the molecular events leading to dormancy are not known. In U. maydis, natural antisense transcripts (NATs) can function to alter gene expression and many NATs have increased levels in the teliospore. One such NAT is as-ssm1 which is complementary to the gene for the mitochondrial seryl-tRNA synthetase (ssm1), an enzyme important to mitochondrial function. The disruption of ssm1 leads to cell lysis, indicating it is also essential for cellular viability. To assess the function of as-ssm1, it was ectopically expressed in haploid cells, where it is not normally present. This expression led to reductions in growth rate, virulence, mitochondrial membrane potential and oxygen consumption. It also resulted in the formation of as-ssm1/ssm1 double-stranded RNA and increased ssm1 transcript levels, but no change in Ssm1 protein levels was detected. Together, these findings suggest a role for as-ssm1 in facilitating teliospore dormancy through dsRNA formation and reduction of mitochondrial function.

Tissue Elasticity Bridges Cancer Stem Cells to the Tumor Microenvironment Through microRNAs: Implications for a "Watch-and-Wait" Approach to Cancer

BACKGROUND

Targeting the tumor microenvironment (TME) through which cancer stem cells (CSCs) crosstalk for cancer initiation and progression, may open new treatments different from those centered on the original hallmarks of cancer genetics thereby implying a new approach for suppression of TME driven activation of CSCs. Cancer is dynamic, heterogeneous, evolving with the TME and can be influenced by tissue-specific elasticity. One of the mediators and modulators of the crosstalk between CSCs and mechanical forces is miRNA, which can be developmentally regulated, in a tissue- and cellspecific manner.

OBJECTIVE

Here, based on our previous data, we provide a framework through which such gene expression changes in response to external mechanical forces can be understood during cancer progression. Recognizing the ways mechanical forces regulate and affect intracellular signals with applications in cancer stem cell biology. Such TME-targeted pathways shed new light on strategies for attacking cancer stem cells with fewer side effects than traditional gene-based treatments for cancer, requiring a "watchand- wait" approach. We attempt to address both normal brain microenvironment and tumor microenvironment as both works together, intertwining in pathology and physiology - a balance that needs to be maintained for the "watch-and-wait" approach to cancer.

CONCLUSION

This review connected the subjects of tissue elasticity, tumor microenvironment, epigenetic of miRNAs, and stem-cell biology that are very relevant in cancer research and therapy. It attempts to unify apparently separate entities in a complex biological web, network, and system in a realistic and practical manner, i.e., to bridge basic research with clinical application.

MiRNA expression profiling in human gliomas: upregulated miR-363 increases cell survival and proliferation

The role of microRNAs (miRNAs) in glioma biology is increasingly recognized. To investigate the regulatory mechanisms governing the malignant signature of gliomas with different grades of malignancy, we analyzed miRNA expression profiles in human grade I-IV tumor samples and primary glioma cell cultures. Multiplex real-time PCR was used to profile miRNA expression in a set of World Health Organization (WHO) grade I (pilocytic astrocytoma), II (diffuse fibrillary astrocytoma), and IV (glioblastoma multiforme) astrocytic tumors and primary glioma cell cultures. Primary glioma cell cultures were used to evaluate the effect of transfection of specific miRNAs and miRNA inhibitors. miRNA microarray showed that a set of miRNAs was consistently upregulated in all glioma samples. miR-363 was upregulated in all tumor specimens and cell lines, and its expression correlated with tumor grading. The transfection of glioma cells with the specific inhibitor of miR-363 increased the expression level of tumor suppressor growth-associated protein 43 (GAP-43). Transfection of miR-363 induced cell survival, while inhibition of miR-363 significantly reduced glioma cell viability. Furthermore, miRNA-363 inhibition induced the downregulation of AKT, cyclin-D1, matrix metalloproteinase (MMP)-2, MMP-9, and Bcl-2 and upregulation of caspase 3. Together, these data suggest that the upregulation of miR-363 may play a role in malignant glioma signature.

Association of new deletion/duplication region at chromosome 1p21 with intellectual disability, severe speech deficit and autism spectrum disorder-like behavior: an all-in approach to solving the DPYD enigma

We describe an as yet unreported neocentric small supernumerary marker chromosome (sSMC) derived from chromosome 1p21.3p21.2. It was present in 80% of the lymphocytes in a male patient with intellectual disability, severe speech deficit, mild dysmorphic features, and hyperactivity with elements of autism spectrum disorder (ASD).

Several important neurodevelopmental genes are affected by the 3.56 Mb copy number gain of 1p21.3p21.2, which may be considered reciprocal in gene content to the recently recognized 1p21.3 microdeletion syndrome. Both 1p21.3 deletions and the presented duplication display overlapping symptoms, fitting the same disorder category. Contribution of coding and non-coding genes to the phenotype is discussed in the light of cellular and intercellular homeostasis disequilibrium. In line with this the presented 1p21.3p21.2 copy number gain correlated to 1p21.3 microdeletion syndrome verifies the hypothesis of a cumulative effect of the number of deregulated genes - homeostasis disequilibrium leading to overlapping phenotypes between microdeletion and microduplication syndromes.

Although miR-137 appears to be the major player in the 1p21.3p21.2 region, deregulation of the DPYD (dihydropyrimidine dehydrogenase) gene may potentially affect neighboring genes underlying the overlapping symptoms present in both the copy number loss and copy number gain of 1p21. Namely, the all-in approach revealed that DPYD is a complex gene whose expression is epigenetically regulated by long non-coding RNAs (lncRNAs) within the locus. Furthermore, the long interspersed nuclear element-1 (LINE-1) L1MC1 transposon inserted in DPYD intronic transcript 1 (DPYD-IT1) lncRNA with its parasites, TcMAR-Tigger5b and pair of Alu repeats appears to be the “weakest link” within the DPYD gene liable to break. Identification of the precise mechanism through which DPYD is epigenetically regulated, and underlying reasons why exactly the break (FRA1E) happens, will consequently pave the way toward preventing severe toxicity to the antineoplastic drug 5-fluorouracil (5-FU) and development of the causative therapy for the dihydropyrimidine dehydrogenase deficiency.

Novel candidate key drivers in the integrative network of genes, microRNAs, methylations, and copy number variations in squamous cell lung carcinoma

The mechanisms of lung cancer are highly complex. Not only mRNA gene expression but also microRNAs, DNA methylation, and copy number variation (CNV) play roles in tumorigenesis. It is difficult to incorporate so much information into a single model that can comprehensively reflect all these lung cancer mechanisms. In this study, we analyzed the 129 TCGA (The Cancer Genome Atlas) squamous cell lung carcinoma samples with gene expression, microRNA expression, DNA methylation, and CNV data. First, we used variance inflation factor (VIF) regression to build the whole genome integrative network. Then, we isolated the lung cancer subnetwork by identifying the known lung cancer genes and their direct regulators. This subnetwork was refined by the Bayesian method, and the directed regulations among mRNA genes, microRNAs, methylations, and CNVs were obtained. The novel candidate key drivers in this refined subnetwork, such as the methylation of ARHGDIB and HOXD3, microRNA let-7a and miR-31, and the CNV of AGAP2, were identified and analyzed. On three large public available lung cancer datasets, the key drivers ARHGDIB and HOXD3 demonstrated significant associations with the overall survival of lung cancer patients. Our results provide new insights into lung cancer mechanisms.

Long Non-coding RNA

Rapid development in genome-wide transcriptional analyses has led to the discovery of a large number of non-coding transcripts, also called long non-coding RNA (lncRNA). LncRNAs harbor biological activities including regulation of protein-coding gene expression at epigenetic, transcriptional and post-transcriptional levels. They also take a part in various physiological and pathological processes, participating in cell development, immunity, disease processes and oncogenesis.

Here I discuss and summarize, current knowledge about lncRNA origin, function and involvement in human disease.

Oligonucleotides for upregulating gene expression

Over the recent decade oligonucleotides have become an important new class of molecules, allowing therapeutic intervention through targets previously thought 'undruggable'. One of the new avenues opened up by oligonucleotide-based drugs was specific gene upregulation, which, historically, has been difficult to achieve using small-molecule drugs. This article will focus on patents covering this important development in the oligonucleotide field and highlight the different mechanisms through which the oligonucleotide-mediated gene upregulation can work, including inhibition of activity of natural antisense transcripts, interaction with promoter binding sites of noncoding regulatory RNAs, blocking of regulatory and/or miRNA binding sites in 3' UTRs, blocking splice inhibitor/enhancer sites or blocking interactions with polycomb repressive complex 2. Understanding the particular mechanism through which an oligonucleotide drug exerts its effects is highly important in drug development, as it determines the design of the drug molecule.

Examining the condition-specific antisense transcription in S. cerevisiae and S. paradoxus

Background

Recent studies have demonstrated that antisense transcription is pervasive in budding yeasts and is conserved between Saccharomyces cerevisiae and S. paradoxus. While studies have examined antisense transcripts of S. cerevisiae for inverse expression in stationary phase and stress conditions, there is a lack of comprehensive analysis of the conditional specific evolutionary characteristics of antisense transcription between yeasts. Here we attempt to decipher the evolutionary relationship of antisense transcription of S. cerevisiae and S. paradoxus cultured in mid log, early stationary phase, and heat shock conditions.

Results

Massively parallel sequencing of sequence strand-specific cDNA library was performed from RNA isolated from S. cerevisiae and S. paradoxus cells at mid log, stationary phase and heat shock conditions. We performed this analysis using a stringent set of sense ORF transcripts and non-coding antisense transcripts that were expressed in all the three conditions, as well as in both species. We found the divergence of the condition-specific anti-sense transcription levels is higher than that in condition-specific sense transcription levels, suggesting that antisense transcription played a potential role in adapting to different conditions. Furthermore, 43% of sense-antisense pairs demonstrated inverse expression in either stationary phase or heat shock conditions relative to the mid log conditions. In addition, a large part of sense-antisense pairs (67%), which demonstrated inverse expression, were highly conserved between the two species. Our results were also concordant with known functional analyses from previous studies and with the evidence from mechanistic experiments of role of individual genes.

Conclusions

By performing a genome-scale computational analysis, we have tried to evaluate the role of antisense transcription in mediating sense transcription under different environmental conditions across and in two related yeast species. Our findings suggest that antisense regulation could control expression of the corresponding sense transcript via inverse expression under a range of different circumstances.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-521) contains supplementary material, which is available to authorized users.

Long non-coding RNA-dependent transcriptional regulation in neuronal development and disease

Comprehensive analysis of the mammalian transcriptome has revealed that long non-coding RNAs (lncRNAs) may make up a large fraction of cellular transcripts. Recent years have seen a surge of studies aimed at functionally characterizing the role of lncRNAs in development and disease. In this review, we discuss new findings implicating lncRNAs in controlling development of the central nervous system (CNS). The evolution of the higher vertebrate brain has been accompanied by an increase in the levels and complexities of lncRNAs expressed within the developing nervous system. Although a limited number of CNS-expressed lncRNAs are now known to modulate the activity of proteins important for neuronal differentiation, the function of the vast majority of neuronal-expressed lncRNAs is still unknown. Topics of intense current interest include the mechanism by which CNS-expressed lncRNAs might function in epigenetic and transcriptional regulation during neuronal development, and how gain and loss of function of individual lncRNAs contribute to neurological diseases.

Natural antisense transcripts and long non-coding RNA in Neurospora crassa

The prevalence of long non-coding RNAs (lncRNA) and natural antisense transcripts (NATs) has been reported in a variety of organisms. While a consensus has yet to be reached on their global importance, an increasing number of examples have been shown to be functional, regulating gene expression at the transcriptional and post-transcriptional level. Here, we use RNA sequencing data from the ABI SOLiD platform to identify lncRNA and NATs obtained from samples of the filamentous fungus Neurospora crassa grown under different light and temperature conditions. We identify 939 novel lncRNAs, of which 477 are antisense to annotated genes. Across the whole dataset, the extent of overlap between sense and antisense transcripts is large: 371 sense/antisense transcripts are complementary over 500 nts or more and 236 overlap by more than 1000 nts. Most prevalent are 3′ end overlaps between convergently transcribed sense/antisense pairs, but examples of divergently transcribed pairs and nested transcripts are also present. We confirm the expression of a subset of sense/antisense transcript pairs by qPCR. We examine the size, types of overlap and expression levels under the different environmental stimuli of light and temperature, and identify 11 lncRNAs that are up-regulated in response to light. We also find differences in transcript length and the position of introns between protein-coding transcripts that have antisense expression and transcripts with no antisense expression. These results demonstrate the ability of N. crassa lncRNAs and NATs to be regulated by different environmental stimuli and provide the scope for further investigation into the function of NATs.

Chromatin-remodeling factors mediate the balance of sense-antisense transcription at the FGF2 locus

Antisense transcription is prevalent in mammalian genomes, yet the function of many antisense transcripts remains elusive. We have previously shown that the fibroblast growth factor 2 (FGF2) gene is regulated endogenously by an overlapping antisense gene called Nudix-type motif 6 (NUDT6). However, the molecular mechanisms that determine the balance of FGF2 and NUDT6 transcripts are not yet well understood. Here we demonstrate that there is a strong negative correlation between FGF2 and NUDT6 across 7 different cell lines. Small interfering RNA-mediated knockdown of NUDT6 causes an increase in nascent FGF2 transcripts, including a short FGF2 variant that lacks sequence complementarity with NUDT6, indicating the involvement of transcriptional mechanisms. In support of this, we show that changes in histone acetylation by trichostatin A treatment, histone deacetylase inhibition, or small interfering RNA knockdown of the histone acetyltransferase CSRP2BP, oppositely affect NUDT6 and FGF2 mRNA levels. A significant increase in histone acetylation with trichostatin A treatment was only detected at the genomic region where the 2 genes overlap, suggesting that this may be an important regulatory region for determining the balance of NUDT6 and FGF2. Knockdown of the histone demethylase KDM4A similarly causes a shift in the balance of NUDT6 and FGF2 transcripts. Expression of CSRP2BP and KDM4A correlates positively with NUDT6 expression and negatively with FGF2 expression. The results presented here indicate that histone acetylation and additional chromatin modifiers are important in determining the relative levels of FGF2 and NUDT6 and support a model in which epigenetic remodeling contributes to their relative expression levels.

Complex gene expression in the dragline silk producing glands of the Western black widow (Latrodectus hesperus)

Background

Orb-web and cob-web weaving spiders spin dragline silk fibers that are among the strongest materials known. Draglines are primarily composed of MaSp1 and MaSp2, two spidroins (spider fibrous proteins) expressed in the major ampullate (MA) silk glands. Prior genetic studies of dragline silk have focused mostly on determining the sequence of these spidroins, leaving other genetic aspects of silk synthesis largely uncharacterized.

Results

Here, we used deep sequencing to profile gene expression patterns in the Western black widow, Latrodectus hesperus. We sequenced millions of 3′-anchored “tags” of cDNAs derived either from MA glands or control tissue (cephalothorax) mRNAs, then associated the tags with genes by compiling a reference database from our newly constructed normalized L. hesperus cDNA library and published L. hesperus sequences. We were able to determine transcript abundance and alternative polyadenylation of each of three loci encoding MaSp1. The ratio of MaSp1:MaSp2 transcripts varied between individuals, but on average was similar to the estimated ratio of MaSp1:MaSp2 in dragline fibers. We also identified transcription of TuSp1 in MA glands, another spidroin family member that encodes the primary component of egg-sac silk, synthesized in tubuliform glands. In addition to the spidroin paralogs, we identified 30 genes that are more abundantly represented in MA glands than cephalothoraxes and represent new candidates for involvement in spider silk synthesis.

Conclusions

Modulating expression rates of MaSp1 variants as well as MaSp2 and TuSp1 could lead to differences in mechanical properties of dragline fibers. Many of the newly identified candidate genes likely encode secreted proteins, suggesting they could be incorporated into dragline fibers or assist in protein processing and fiber assembly. Our results demonstrate previously unrecognized transcript complexity in spider silk glands.

Mechanisms of antisense transcription initiation from the 3' end of the GAL10 coding sequence in vivo

In spite of the important regulatory functions of antisense transcripts in gene expression, it remains unknown how antisense transcription is initiated. Recent studies implicated RNA polymerase II in initiation of antisense transcription. However, how RNA polymerase II is targeted to initiate antisense transcription has not been elucidated. Here, we have analyzed the association of RNA polymerase II with the antisense initiation site at the 3' end of the GAL10 coding sequence in dextrose-containing growth medium that induces antisense transcription. We find that RNA polymerase II is targeted to the antisense initiation site at GAL10 by Reb1p activator as well as general transcription factors (e.g., TFIID, TFIIB, and Mediator) for antisense transcription initiation. Intriguingly, while GAL10 antisense transcription is dependent on TFIID, its sense transcription does not require TFIID. Further, the Gal4p activator that promotes GAL10 sense transcription is dispensable for antisense transcription. Moreover, the proteasome that facilitates GAL10 sense transcription does not control its antisense transcription. Taken together, our results reveal that GAL10 sense and antisense transcriptions are regulated differently and shed much light on the mechanisms of antisense transcription initiation.

Natural antisense transcripts in plants: a review and identification in soybean infected with Phakopsora pachyrhizi SuperSAGE library

Natural antisense ranscripts (NAT) are RNA molecules complementary to other endogenous RNAs. They are capable of regulating the expression of target genes at different levels (transcription, mRNA stability, translation, etc.). Such a property makes them ideal for interventions in organisms' metabolism. The present study reviewed plant NAT aspects, including features, availability and genesis, conservation and distribution, coding capacity, NAT pair expression, and functions. Besides, an in silico identification of NATs pairs was presented, using deepSuperSAGE libraries of soybean infected or not with Phakopsora pachyrhizi. Results showed that around 1/3 of the 77,903 predicted trans-NATs (by PlantsNATsDB database) detected had unitags mapped in both sequences of each pair. The same 1/3 of the 436 foreseen cis-NATs showed unitags anchored in both sequences of the related pairs. For those unitags mapped in NAT pairs, a modulation expression was assigned as upregulated, downregulated, or constitutive, based on the statistical analysis (P < 0.05). As a result, the infected treatment promoted the expression of 2,313 trans-NATs pairs comprising unitags exclusively from that library (1,326 pairs had unitags only found in the mock library). To understand the regulation of these NAT pairs could be a key aspect in the ASR plant response.

Multiple insert size paired-end sequencing for deconvolution of complex transcriptomes

Deep sequencing of transcriptomes allows quantitative and qualitative analysis of many RNA species in a sample, with parallel comparison of expression levels, splicing variants, natural antisense transcripts, RNA editing and transcriptional start and stop sites the ideal goal. By computational modeling, we show how libraries of multiple insert sizes combined with strand-specific, paired-end (SS-PE) sequencing can increase the information gained on alternative splicing, especially in higher eukaryotes. Despite the benefits of gaining SS-PE data with paired ends of varying distance, the standard Illumina protocol allows only non-strand-specific, paired-end sequencing with a single insert size. Here, we modify the Illumina RNA ligation protocol to allow SS-PE sequencing by using a custom pre-adenylated 3' adaptor. We generate parallel libraries with differing insert sizes to aid deconvolution of alternative splicing events and to characterize the extent and distribution of natural antisense transcription in C. elegans. Despite stringent requirements for detection of alternative splicing, our data increases the number of intron retention and exon skipping events annotated in the Wormbase genome annotations by 127% and 121%, respectively. We show that parallel libraries with a range of insert sizes increase transcriptomic information gained by sequencing and that by current established benchmarks our protocol gives competitive results with respect to library quality.

Increased expression of P-glycoprotein and doxorubicin chemoresistance of metastatic breast cancer is regulated by miR-298

MicroRNAs (miRNAs) are short, noncoding RNA molecules that regulate the expression of a number of genes involved in cancer; therefore, they offer great diagnostic and therapeutic targets. We have developed doxorubicin-resistant and -sensitive metastatic human breast cancer cell lines (MDA-MB-231) to study the chemoresistant mechanisms regulated by miRNAs. We found that doxorubicin localized exclusively to the cytoplasm and was unable to reach the nuclei of resistant tumor cells because of the increased nuclear expression of MDR1/P-glycoprotein (P-gp). An miRNA array between doxorubicin-sensitive and -resistant breast cancer cells showed that reduced expression of miR-298 in doxorubicin-resistant human breast cancer cells was associated with increased expression of P-gp. In a transient transfection experiment, miR-298 directly bound to the MDR1 3' untranslated region and regulated the expression of firefly luciferase reporter in a dose-dependent manner. Overexpression of miR-298 down-regulated P-gp expression, increasing nuclear accumulation of doxorubicin and cytotoxicity in doxorubicin-resistant breast cancer cells. Furthermore, down-regulation of miR-298 increased P-gp expression and induced doxorubicin resistance in sensitive breast cancer cells. In summary, these results suggest that miR-298 directly modulates P-gp expression and is associated with the chemoresistant mechanisms of metastatic human breast cancer. Therefore, miR-298 has diagnostic and therapeutic potential for predicting doxorubicin chemoresistance in human breast cancer.

Comparative analysis of neural transcriptomes and functional implication of unannotated intronic expression

Background

The transcriptome and its regulation bridge the genome and the phenome. Recent RNA-seq studies unveiled complex transcriptomes with previously unknown transcripts and functions. To investigate the characteristics of neural transcriptomes and possible functions of previously unknown transcripts, we analyzed and compared nine recent RNA-seq datasets corresponding to tissues/organs ranging from stem cell, embryonic brain cortex to adult whole brain.

Results

We found that the neural and stem cell transcriptomes share global similarity in both gene and chromosomal expression, but are quite different from those of liver or muscle. We also found an unusually high level of unannotated expression in mouse embryonic brains. The intronic unannotated expression was found to be strongly associated with genes annotated for neurogenesis, axon guidance, negative regulation of transcription, and neural transmission. These functions are the hallmarks of the late embryonic stage cortex, and crucial for synaptogenesis and neural circuit formation.

Conclusions

Our results revealed unique global and local landscapes of neural transcriptomes. It also suggested potential functional roles for previously unknown transcripts actively expressed in the developing brain cortex. Our findings provide new insights into potentially novel genes, gene functions and regulatory mechanisms in early brain development.

Wide-ranging functions of E2F4 in transcriptional activation and repression revealed by genome-wide analysis

The E2F family of transcription factors has important roles in cell cycle progression. E2F4 is an E2F family member that has been proposed to be primarily a repressor of transcription, but the scope of its binding activity and functions in transcriptional regulation is not fully known. We used ChIP sequencing (ChIP-seq) to identify around 16 000 E2F4 binding sites which potentially regulate 7346 downstream target genes with wide-ranging functions in DNA repair, cell cycle regulation, apoptosis, and other processes. While half of all E2F4 binding sites (56%) occurred near transcription start sites (TSSs), ∼20% of sites occurred more than 20 kb away from any annotated TSS. These distal sites showed histone modifications suggesting that E2F4 may function as a long-range regulator, which we confirmed by functional experimental assays on a subset. Overexpression of E2F4 and its transcriptional cofactors of the retinoblastoma (Rb) family and its binding partner DP-1 revealed that E2F4 acts as an activator as well as a repressor. E2F4 binding sites also occurred near regulatory elements for miRNAs such as let-7a and mir-17, suggestive of regulation of miRNAs by E2F4. Taken together, our genome-wide analysis provided evidence of versatile roles of E2F4 and insights into its functions.

RNA-Seq analysis of splicing in Plasmodium falciparum uncovers new splice junctions, alternative splicing and splicing of antisense transcripts

Over 50% of genes in Plasmodium falciparum, the deadliest human malaria parasite, contain predicted introns, yet experimental characterization of splicing in this organism remains incomplete. We present here a transcriptome-wide characterization of intraerythrocytic splicing events, as captured by RNA-Seq data from four timepoints of a single highly synchronous culture. Gene model-independent analysis of these data in conjunction with publically available RNA-Seq data with HMMSplicer, an in-house developed splice site detection algorithm, revealed a total of 977 new 5' GU-AG 3' and 5 new 5' GC-AG 3' junctions absent from gene models and ESTs (11% increase to the current annotation). In addition, 310 alternative splicing events were detected in 254 (4.5%) genes, most of which truncate open reading frames. Splicing events antisense to gene models were also detected, revealing complex transcriptional arrangements within the parasite's transcriptome. Interestingly, antisense introns overlap sense introns more than would be expected by chance, perhaps indicating a functional relationship between overlapping transcripts or an inherent organizational property of the transcriptome. Independent experimental validation confirmed over 30 new antisense and alternative junctions. Thus, this largest assemblage of new and alternative splicing events to date in Plasmodium falciparum provides a more precise, dynamic view of the parasite's transcriptome.

Long Nonprotein-Coding RNAs in Plants

In recent years, nonprotein-coding RNAs (or npcRNAs) have emerged as a major part of the eukaryotic transcriptome. Many new regulatory npcRNAs or riboregulators riboregulators have been discovered and characterized due to the advent of new genomic approaches. This growing number suggests that npcRNAs could play a more important role than previously believed and significantly contribute to the generation of evolutionary complexity in multicellular organisms. Regulatory npcRNAs range from small RNAs (si/miRNAs) to very large transcripts (or long npcRNAs) and play diverse functions in development and/or environmental stress responses. Small RNAs include an expanding number of 20-40 nt RNAs that function in the regulation of gene expression by affecting mRNA decay and translational inhibition or lead to DNA methylation and gene silencing. They generally involve double-stranded RNA or stem loops and imply transcriptional or posttranscriptional gene silencing (PTGS). RNA silencing besides small interfering RNA and microRNA, gene silencing in plants is also mediated by tasiRNAs (trans-acting siRNAs) and nat-siRNAs (natural antisense mediated siRNAs). In contrast to small RNAs, much less is known about the large and diverse population of long npcRNAs, and only a few have been implicated in diverse functions such as abiotic stress responses, nodulation and flower development, and sex chromosome-specific expression. Moreover, many long npcRNAs act as antisense transcripts or are substrates of the small RNA pathways, thus interfering with a variety of RNA-related metabolisms. An emerging hypothesis is that long npcRNAs, as shown for small si/miRNAs, integrate into ribonucleoprotein particles (RNPs) to modulate their function, localization, or stability to act on target mRNAs. As plants show a remarkable developmental plasticity to adapt their growth to changing environmental conditions, understanding how npcRNAs work may reveal novel mechanisms involved in growth control and differentiation and help to design new tools for biotechnological applications.

Salivary gland hypertrophy viruses: a novel group of insect pathogenic viruses

Salivary gland hypertrophy viruses (SGHVs) are a unique, unclassified group of entomopathogenic, double-stranded DNA viruses that have been reported from three genera of Diptera. These viruses replicate in nuclei of salivary gland cells in adult flies, inducing gland enlargement with little obvious external disease symptoms. Viral infection inhibits reproduction by suppressing vitellogenesis, causing testicular aberrations, and/or disrupting mating behavior. Historical and present research findings support a recent proposal of a new virus family, the Hytrosaviridae. This review describes the discovery and prevalence of different SGHVs, summarizes their biochemical characterization and taxonomy, compares morphological and histopathological properties, and details transmission routes and the influence of infection on host biology and reproduction. In addition, the potential use of SGHVs as sterilizing agents for house fly control and the deleterious impact of SGHVs on colonized tsetse flies reared for sterile insect technique are discussed.

Intronic AT skew is a defendable proxy for germline transcription but does not predict crossing-over or protein evolution rates in Drosophila melanogaster

Recent evidence suggests that germline transcription may affect both protein evolutionary rates, possibly mediated by repair processes, and recombination rates, possibly mediated by chromatin and epigenetic modification. Here, we test these propositions in Drosophila melanogaster. The challenge for such analyses is to provide defendable measures of germline gene expression. Intronic AT skew is a good candidate measure as it is thought to be a consequence, at least in part, of transcription-coupled repair. Prior evidence suggests that intronic AT skew in D. melanogaster is not affected by proximity to intron extremities and differs between transcribed DNA and flanking sequence. We now also establish that intronic AT skew is a defendable proxy for germline expression as (a) it is more similar than expected by chance between introns of the same gene (which is not accounted for by physical proximity), (b) is correlated with male germline expression, and (c) is more pronounced in broadly expressed genes. Furthermore, (d) a trend for intronic skew to differ between 3' and 5' ends of genes is particular to broadly expressed genes. Finally, (e) controlling for physical distance, introns of proximate genes are most different in skew if they have different tissue specificity. We find that intronic AT skew, employed as a proxy for germline transcription, correlates neither with recombination rates nor with the rate of protein evolution. We conclude that there is no prima facie evidence that germline expression modulates recombination rates or monotonically affects protein evolution rates in D. melanogaster.

Transcriptional regulation by small RNAs at sequences downstream from 3' gene termini

Transcriptome studies reveal many noncoding transcripts overlapping 3' gene termini. The function of these transcripts is unknown. Here we have characterized transcription at the progesterone receptor (PR) locus and identified noncoding transcripts that overlap the 3' end of the gene. Small RNAs complementary to sequences beyond the 3' terminus of PR mRNA modulated expression of PR, recruited argonaute 2 to a 3' noncoding transcript, altered occupancy of RNA polymerase II, induced chromatin changes at the PR promoter and affected responses to physiological stimuli. We found that the promoter and 3' terminal regions of the PR locus are in close proximity, providing a potential mechanism for RNA-mediated control of transcription over long genomic distances. These results extend the potential for small RNAs to regulate transcription to target sequences beyond the 3' termini of mRNA.

Natural antisense transcripts regulate gene expression in an epigenetic manner

Cytosine DNA methylation, covalent histone modifications, and RNA-mediated gene regulation are the major aspects of epigenetic regulation. Natural antisense transcripts (NATs), as a new member of regulatory RNAs, occur ubiquitously in prokaryote and eukaryote, and play significant roles in physiological or pathological processes. NATs, mostly non-coding RNAs, are involved in transcriptional interference, genomic imprinting, X inactivation, RNA editing, translational regulation, RNA export, DNA methylation, histone modifications, and so on. NATs regulate gene expression through direct interaction with the sense transcripts or indirect interaction with other targets, such as DNA methyltransferases, histone acetylases and histone deacetylases. There may be a direct link among NATs, DNA methylaton and histone modifications. Through formation of sense-antisense duplex structures, NATs exert a widespread impact on conventional gene expression at the mRNA and/or protein level and regulate sense transcripts in a concordant or discordant manner. As one of the important components in epigenetics, NATs could be a potentially rich source for scientists to exploit in the therapy of cancers and other diseases.

Role of inflammation and oxidative stress mediators in gliomas

Gliomas are the most common primary brain tumors of the central nervous system. Despite relevant progress in conventional treatments, the prognosis of such tumors remains almost invariably dismal. The genesis of gliomas is a complex, multistep process that includes cellular neoplastic transformation, resistance to apoptosis, loss of control of the cell cycle, angiogenesis, and the acquisition of invasive properties. Among a number of different biomolecular events, the existence of molecular connections between inflammation and oxidative stress pathways and the development of this cancer has been demonstrated. In particular, the tumor microenvironment, which is largely orchestrated by inflammatory molecules, is an indispensable participant in the neoplastic process, promoting proliferation, survival and migration of such tumors. Proinflammatory cytokines, such as tumor necrosis factor-alpha, interleukin-1beta, and interferon-gamma, as well as chemokines and prostaglandins, are synthesized by resident brain cells and lymphocytes invading the affected brain tissue. Key mediators of cancer progression include nuclear factor-kappaB, reactive oxygen and nitrogen species, and specific microRNAs. The collective activity of these mediators is largely responsible for a pro-tumorigenic response through changes in cell proliferation, cell death, cellular senescence, DNA mutation rates, DNA methylation and angiogenesis. We provide a general overview of the connection between specific inflammation and oxidative stress pathway molecules and gliomas. The elucidation of specific effects and interactions of these factors may provide the opportunity for the identification of new target molecules leading to improved diagnosis and treatment.

Detection of novel 3' untranslated region extensions with 3' expression microarrays

Background

The 3' untranslated regions (UTRs) of transcripts are not well characterized for many genes and often extend beyond the annotated regions. Since Affymetrix 3' expression arrays were designed based on expressed sequence tags, many probesets map to intergenic regions downstream of genes. We used expression information from these probesets to predict transcript extension beyond currently known boundaries.

Results

Based on our dataset encompassing expression in 22 different murine tissues, we identified 845 genes with predicted 3'UTR extensions. These extensions have a similar conservation as known 3'UTRs, which is distinctly higher than intergenic regions. We verified 8 of the predictions by PCR and found all of the predicted regions to be expressed. The method can be extended to other 3' expression microarray platforms as we demonstrate with human data. Additional confirming evidence was obtained from public paired end read data.

Conclusions

We show that many genes have 3'UTR regions extending beyond currently known gene regions and provide a method to identify such regions based on microarray expression data. Since 3' UTR contain microRNA binding sites and other stability determining regions, identification of the full length 3' UTR is important to elucidate posttranscriptional regulation.

The expression of ELK transcription factors in adult DRG: Novel isoforms, antisense transcripts and upregulation by nerve damage

ELK transcription factors are known to be expressed in a number of regions in the nervous system. We show by RT-PCR that the previously described Elk1, Elk3/Elk3b/Elk3c and Elk4 mRNAs are expressed in adult dorsal root ganglia (DRG), together with the novel alternatively spliced isoforms Elk1b, Elk3d and Elk4c/Elk4d/Elk4e. These isoforms are also expressed in brain, heart, kidney and testis. In contrast to Elk3 protein, the novel Elk3d isoform is cytoplasmic, fails to bind ETS binding sites and yet can activate transcription by an indirect mechanism. The Elk3 and Elk4 genes are overlapped by co-expressed Pctk2 (Cdk17) and Mfsd4 genes, respectively, with the potential formation of Elk3/Pctaire2 and Elk4/Mfsd4 sense-antisense mRNA heteroduplexes. After peripheral nerve injury the Elk3 mRNA isoforms are each upregulated approximately 2.3-fold in DRG (P<0.005), whereas the natural antisense Pctaire2 isoforms show only a small increase (21%, P<0.01) and Elk1 and Elk4 mRNAs are unchanged.

Regulation of fibroblast growth factor-2 by an endogenous antisense RNA and by argonaute-2

We have previously reported that elevated fibroblast growth factor-2 (FGF-2) expression is associated with tumor recurrence and reduced survival after surgical resection of esophageal cancer and that these risks are reduced in tumors coexpressing an endogenous antisense (FGF-AS) RNA. In the present study, we examined the role of the endogenous FGF-AS transcript in the regulation of FGF-2 expression in the human lung adenocarcinoma cell line Seg-1. FGF-2 and FGF-AS were temporally and spatially colocalized in the cytoplasm of individual cells, and knockdown of either FGF-2 or FGF-AS by target-specific siRNAs resulted in dose-dependent up-regulation of the complementary transcript and its encoded protein product. Using a luciferase reporter system, we show that these effects are mediated by interaction of the endogenous antisense RNA with the 3'-untranslated region of the FGF-2 mRNA. Deletion mapping identified a 392-nucleotide sequence in the 5823-nucleotide FGF-2 untranslated tail that is targeted by FGF-AS. Small interfering RNA-mediated knockdown of either FGF-AS or FGF-2 significantly increased the stability of the complementary partner mRNA, demonstrating that these mRNAs are mutually regulatory. Knockdown of FGF-AS also resulted in reduced expression of argonaute-2 (AGO-2) and a number of other elements of the endogenous micro-RNA/RNA interference pathways. Conversely, small interfering RNA-mediated knockdown of AGO-2 significantly increased the stability of the FGF-2 mRNA transcript and the steady-state levels of both FGF-2 mRNA and protein, suggesting a role for AGO-2 in the regulation of FGF-2 expression.

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

MicroRNA involvement in the pathogenesis and management of breast cancer

MicroRNAs (miRNAs) are a highly abundant class of endogenous small non-coding RNAs (18-25 nucleotides in length) that regulate gene expression by targeting protein-coding mRNAs post-transcriptionally. miRNAs have been implicated in cancer development and progression. As miRNAs and their regulatory functions are further revealed, the more the importance of miRNA-directed gene regulation is emphasised. In the human genome, 695 mature miRNAs have been identified, although computational calculation predicts that this may increase to >1000. Deregulation of miRNA expression profiles is thought to be implicated in the pathogenesis of many human cancers including breast tumours. Breast cancer subtypes are observed to have deranged miRNA expression signatures, which makes miRNAs important targets for developing a novel molecular classification of breast cancer and opening avenues for more individualised treatment strategies for patients with breast cancer.

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.

Regulatory roles of natural antisense transcripts

Mammalian genomes encode numerous natural antisense transcripts, but the function of these transcripts is not well understood. Functional validation studies indicate that antisense transcripts are not a uniform group of regulatory RNAs but instead belong to multiple categories with some common features. Recent evidence indicates that antisense transcripts are frequently functional and use diverse transcriptional and post-transcriptional gene regulatory mechanisms to carry out a wide variety of biological roles.

Autoregulatory loop of Msx1 expression involving its antisense transcripts

The Msx1 homeogene plays an important role in epithelial-mesenchymal interactions leading organogenesis. Msx1 gene is submitted to bidirectional transcription generating a long non-coding antisense (AS) RNA potentially involved in Msx1 expression regulation. RT-Q-PCR and RNA-FISH studies indicated that transient overexpression of the Msx1 AS transcript in 705IC5 mouse odontoblasts decreased the abundance of endogenous Msx1 S mRNA at the post-transcriptional level. Conversely, Msx1 overexpression increased the AS RNA level probably by activating AS transcription. In vivo mapping by RT-PCR evidenced both Msx1 RNAs in all adult mouse tissues tested raising the issue of Msx1 function during adulthood. The expression patterns of the two RNAs were similar, confirming the tight S/AS relationship. In particular, both Msx1 mRNAs and Msx1 protein were similarly distributed in eyes, and were found in regions with a common ectodermic origin and in cells potentially involved in regeneration. In conclusion, we report that Msx1 S RNA is negatively controlled by its AS RNA at a post-transcriptional level, and that the AS RNA is retrocontrolled positively by Msx1. The tight link between Msx1 S and AS RNAs constitutes a regulatory loop resulting in a fine-tuned expression of Msx1 which appears to be significant for adult homeostasis.

Gene regulation by sense-antisense overlap of polyadenylation signals

We show here that expression of genes from convergent transcription units can be regulated by the formation of double-stranded RNA (dsRNA) in the region of overlapping polyadenylation signals. The model system employed is the mouse polyomavirus. The early and late genes of polyomavirus are transcribed from opposite strands of the circular viral genome. At early times after infection, the early genes are expressed predominantly. Late gene expression increases dramatically upon the onset of DNA replication, when a major defect in polyadenylation of the late primary transcripts generates multigenomic RNAs that are precursors to the mature late mRNAs. Embedded in these late pre-mRNAs are sequences complementary to the early RNAs that act to down-regulate early gene expression via A-to-I editing of dsRNAs. In this system, the defective polyadenylation, and consequently the production of multigenomic late RNAs, depends on the context, and perhaps also, on the A-to-I editing of the poly(A) signal that overlaps the 3'-end of early transcripts.

Analysis of transcripts from predicted open reading frames of Musca domestica salivary gland hypertrophy virus

The Musca domestica salivary gland hypertrophy virus (MdSGHV) is a large dsDNA virus that infects and sterilizes adult houseflies. The transcriptome of this newly described virus was analysed by rapid amplification of cDNA 3'-ends (3'-RACE) and RT-PCR. Direct sequencing of 3'-RACE products revealed 78 poly(A) transcripts containing 95 of the 108 putative ORFs. An additional six ORFs not amplified by 3'-RACE were detected by RT-PCR. Only seven of the 108 putative ORFs were not amplified by either 3'-RACE or RT-PCR. A series of 5'-RACE reactions were conducted on selected ORFs that were identified by 3'-RACE to be transcribed in tandem (tandem transcripts). In the majority of cases, the downstream ORFs were detected as single transcripts as well as components of the tandem transcripts, whereas the upstream ORFs were found only in tandem transcripts. The only exception was the upstream ORF MdSGHV084, which was differentially transcribed as a single transcript at 1 and 2 days post-infection (days p.i.) and as a tandem transcript (MdSGHV084/085) at 2 days p.i. Transcriptome analysis of MdSGHV detected splicing in the 3' untranslated region (3'-UTR) and extensive heterogeneity in the polyadenylation signals and cleavage sites. In addition, 23 overlapping antisense transcripts were found. In conclusion, sequencing the 3'-RACE products without cloning served as an alternative approach to detect both 3'-UTRs and transcript variants of this large DNA virus.

miR-21 and 221 upregulation and miR-181b downregulation in human grade II-IV astrocytic tumors

MicroRNAs (miRNAs) are small noncoding regulatory RNAs that reduce stability and/or translation of fully or partially sequence-complementary target mRNAs. Recent evidence indicates that miRNAs can function both as tumor suppressors and as oncogenes. It has been demonstrated that in glioblastoma multiforme miR-21 and 221 are upregulated whereas miR-128 and 181 are downregulated. Expression of miR-21, 221, 128a, 128b, 128c, 181a, 181b, 181c was studied using real-time quantitative reverse transcriptase polymerase chain reaction and northern blotting for human astrocytic tumors with different grade of malignancy. miR-21 and 221 were overexpressed in glioma samples, whereas miRNA 181b was downregulated compared with normal brain tissue. miRNA-21 was hyperexpressed in all tumor samples whereas higher levels of miRNA-221 were found in high-grade gliomas. This study is the first analysis of miRNAs in astrocytic tumor at different stages of malignancy. The different expression pattern observed in tumors at different stages of malignancy is probably dependent on the cell-specific repertoire of target genes of tumors sharing different molecular pathways activity and suggests miRNAs may have also a place in diagnosis and staging of brain tumors.

Variation in neighbouring genes of the dopaminergic and serotonergic systems affects feather pecking behaviour of laying hens

Feather pecking is a behavioural disorder of laying hens and has serious animal welfare and economic implications. One of the several aetiological hypotheses proposes that the disorder results from redirected exploratory behaviour. Variation in the gene encoding the dopamine D4 receptor (DRD4) has been shown to be associated with exploratory behaviour in several species, including in a passerine bird species. We therefore considered DRD4 as a candidate gene for feather pecking. We have annotated DRD4 in the chicken genome and have re-sequenced it in 140 animals belonging to: experimental layer lines divergently selected for high and low propensity to feather pecking; the unselected founder population; and two commercial lines with low and high propensity to feather pecking. We have identified two sub-haplotypes of DRD4 that are highly significantly associated with feather pecking behaviour in the experimental (P = 7.30 x 10(-7)) as well as in the commercial lines (P = 2.78 x 10(-6)). Linkage disequilibrium (LD) extends into a neighbouring gene encoding deformed epidermal autoregulatory factor 1 (DEAF1). The product of DEAF1 regulates the transcription of the gene encoding the serotonin (5-hydroxytryptamine) 1A receptor. Thus, DEAF1 represents another candidate gene for feather pecking. Re-sequencing of five animals homozygous for the 'low-pecking' sub-haplotype and of six animals homozygous for the 'high-pecking' sub-haplotype delineated an LD block of 14 833 bases spanning the two genes. None of the variants in the LD block is obviously functional. However, the haplotype information will be useful to select against the propensity to feather pecking in chicken and to elucidate the functional implications of the variants.

Sense and antisense OsDof12 transcripts in rice

Background

Antisense transcription is a widespread phenomenon in plants and mammals. Our previous data on rice gene expression analysis by microarray indicated that the sense and antisense transcripts at the OsDof12 locus were co-expressed in leaves. In current study, we analyzed the expression patterns in detail and looked for the possible mechanism related to their expression patterns.

Results

OsDof12, being a single copy gene located on rice chromosome 3, encodes a predicted Dof protein of 440 amino acids with one intron of 945 bp. The antisense transcript, OsDofl2os, overlaps with both the exonic and intronic regions of OsDof12 and encodes a functionally unknown protein of 104 amino acids with no intron. The sense-antisense OsDof12 transcripts were co-expressed within the same tissues, and their expressions were not tissue-specific in general. At different developmental stages in rice, the OsDof12 and OsDof12os transcripts exhibited reciprocal expression patterns. Interestingly, the expression of both genes was significantly induced under drought treatment, and inhibited by dark treatment. In the Pro_OsDof12-GUS and Pro_OsDof12os-GUS transgenic rice plants, the expression profiles of GUS were consistent with those of the OsDof12 and OsDof12os transcripts, respectively. In addition, the analysis of cis-regulatory elements indicated that either of the two promoters contained 74 classes of cis-regulatory elements predicted, of which the two promoter regions shared 53 classes.

Conclusion

Based on the expression profiles of OsDof12 and OsDof12os, the expression patterns of GUS in the Pro_OsDof12-GUS and Pro_OsDof12os-GUS transgenic rice plants and the predicted common cis-regulatory elements shared by the two promoters, we suggest that the co-expression patterns of OsDof12 and OsDof12os might be attributed to the basically common nature of the two promoters.

Evolution of closely linked gene pairs in vertebrate genomes

The orientation of closely linked genes in mammalian genomes is not random: there are more head-to-head (h2h) gene pairs than expected. To understand the origin of this enrichment in h2h gene pairs, we have analyzed the phylogenetic distribution of gene pairs separated by less than 600 bp of intergenic DNA (gene duos). We show here that a lack of head-to-tail (h2t) gene duos is an even more distinctive characteristic of mammalian genomes, with the platypus genome as the only exception. In nonmammalian vertebrate and in nonvertebrate genomes, the frequency of h2h, h2t, and tail-to-tail (t2t) gene duos is close to random. In tetrapod genomes, the h2t and t2t gene duos are more likely to be part of a larger gene cluster of closely spaced genes than h2h gene duos; in fish and urochordate genomes, the reverse is seen. In human and mouse tissues, the expression profiles of gene duos were skewed toward positive coexpression, irrespective of orientation. The organization of orthologs of both members of about 40% of the human gene duos could be traced in other species, enabling a prediction of the organization at the branch points of gnathostomes, tetrapods, amniotes, and euarchontoglires. The accumulation of h2h gene duos started in tetrapods, whereas that of h2t and t2t gene duos only started in amniotes. The apparent lack of evolutionary conservation of h2t and t2t gene duos relative to that of h2h gene duos is thus a result of their relatively late origin in the lineage leading to mammals; we show that once they are formed h2t and t2t gene duos are as stable as h2h gene duos.