Genome-wide view and characterization of natural antisense transcripts in Cannabis Sativa L.

Natural Antisense Transcripts (NATs) are a kind of complex regulatory RNAs that play crucial roles in gene expression and regulation. However, the NATs in Cannabis Sativa L., a widely economic and medicinal plant rich in cannabinoids remain unknown. In this study, we comprehensively predicted C. sativa NATs genome-wide using strand-specific RNA sequencing (ssRNA-Seq) data, and validated the expression profiles by strand-specific quantitative reverse transcription PCR (ssRT-qPCR). Consequently, a total of 307 NATs were predicted in C. sativa, including 104 cis- and 203 trans- NATs. Functional enrichment analysis demonstrated the potential involvement of the C. sativa NATs in DNA polymerase activity, RNA-DNA hybrid ribonuclease activity, and nucleic acid binding. Finally, 18 cis- and 376 trans- NAT-ST pairs were predicted to produce 621 cis- and 5,679 trans- small interfering RNA (nat-siRNAs), respectively. These nat-siRNAs were potentially involved in the biosynthesis of cannabinoids and cellulose. All these results will shed light on the regulation of NATs and nat-siRNAs in C. sativa.

Detection of Zika Virus Replication in Human Semen by Reverse-Transcription Polymerase Chain Reaction Targeting of Antisense Ribonucleic Acid.

BACKGROUND: Persistence of Zika virus (ZIKV) ribonucleic acid (RNA) in semen is common after infection. METHODS: We designed a reverse-transcription polymerase chain reaction assay that targets antisense ZIKV RNA (asRNA) to assess ZIKV replication competence in ZIKV RNA-positive semen samples. RESULTS: We detected ZIKV asRNA in semen of 9 of 19 men (47.4%) diagnosed with ZIKV infection. All asRNA-positive samples had high ZIKV loads (cycle threshold values <26) and were obtained within 21 days of symptom onset. CONCLUSIONS: The sensitivity of the asRNA assay for detection of ZIKV replication was higher than that of conventional virus isolation methods (47.4% vs 21.1%, P = .032).

The Suf Iron-Sulfur Cluster Biosynthetic System Is Essential in Staphylococcus aureus, and Decreased Suf Function Results in Global Metabolic Defects and Reduced Survival in Human Neutrophils.

Staphylococcus aureus remains a causative agent for morbidity and mortality worldwide. This is in part a result of antimicrobial resistance, highlighting the need to uncover novel antibiotic targets and to discover new therapeutic agents. In the present study, we explored the possibility that iron-sulfur (Fe-S) cluster synthesis is a viable antimicrobial target. RNA interference studies established that Suf (sulfur mobilization)-dependent Fe-S cluster synthesis is essential in S. aureus We found that sufCDSUB were cotranscribed and that suf transcription was positively influenced by sigma factor B. We characterized an S. aureus strain that contained a transposon inserted in the intergenic space between sufC and sufD (sufD*), resulting in decreased transcription of sufSUB Consistent with the transcriptional data, the sufD* strain had multiple phenotypes associated with impaired Fe-S protein maturation. They included decreased activities of Fe-S cluster-dependent enzymes, decreased growth in media lacking metabolites that require Fe-S proteins for synthesis, and decreased flux through the tricarboxylic acid (TCA) cycle. Decreased Fe-S cluster synthesis resulted in sensitivity to reactive oxygen and reactive nitrogen species, as well as increased DNA damage and impaired DNA repair. The sufD* strain also exhibited perturbed intracellular nonchelated Fe pools. Importantly, the sufD* strain did not exhibit altered exoprotein production or altered biofilm formation, but it was attenuated for survival upon challenge by human polymorphonuclear leukocytes. The results presented are consistent with the hypothesis that Fe-S cluster synthesis is a viable target for antimicrobial development.

The distribution and evolution of Arabidopsis thaliana cis natural antisense transcripts.

BACKGROUND: Natural antisense transcripts (NATs) are regulatory RNAs that contain sequence complementary to other RNAs, these other RNAs usually being messenger RNAs. In eukaryotic genomes, cis-NATs overlap the gene they complement. RESULTS: Here, our goal is to analyze the distribution and evolutionary conservation of cis-NATs for a variety of available data sets for Arabidopsis thaliana, to gain insights into cis-NAT functional mechanisms and their significance. Cis-NATs derived from traditional sequencing are largely validated by other data sets, although different cis-NAT data sets have different prevalent cis-NAT topologies with respect to overlapping protein-coding genes. A. thaliana cis-NATs have substantial conservation (28-35% in the three substantive data sets analyzed) of expression in A. lyrata. We examined evolutionary sequence conservation at cis-NAT loci in Arabidopsis thaliana across nine sequenced Brassicaceae species (picked for optimal discernment of purifying selection), focussing on the parts of their sequences not overlapping protein-coding transcripts (dubbed 'NOLPs'). We found significant NOLP sequence conservation for 28-34% NATs across different cis-NAT sets. This NAT NOLP sequence conservation versus A. lyrata is generally significantly correlated with conservation of expression. We discover a significant enrichment of transcription factor binding sites (as evidenced by CHIP-seq data) in NOLPs compared to randomly sampled near-gene NOLP-like DNA , that is linked to significant sequence conservation. Conversely, there is no such evidence for a general significant link between NOLPs and formation of small interfering RNAs (siRNAs), with the substantial majority of unique siRNAs arising from the overlapping portions of the cis-NATs. CONCLUSIONS: In aggregate, our results suggest that many cis-NAT NOLPs function in the regulation of conserved promoter/regulatory elements that they 'over-hang'.

Genome-wide comparative analysis of digital gene expression tag profiles during maize ear development.

The present study profiled and analyzed gene expression of the maize ear at four key developmental stages. Based on genome-wide profile analysis, we detected differential mRNA of maize genes. Some of the differentially expressed genes (DEGs) were predicted to be potential candidates of maize ear development. Several well-known genes were found with reported mutant analyses, such as, compact plant2 (ct2), zea AGAMOUS homolog1 (zag1), bearded ear (bde), and silky1 (si1). MicroRNAs such as microRNA156 were predicted to target genes involved in maize ear development. Antisense transcripts were widespread throughout all the four stages, and are suspected to play important roles in maize ear development. Thus, identification and characterization of important genes and regulators at all the four developmental stages will contribute to an improved understanding of the molecular mechanisms responsible for maize ear development.

Natural antisense transcripts in Plasmodium falciparum isolates from patients with complicated malaria.

Mechanisms regulating gene expression in malaria parasites are not well understood. Little is known about how the parasite regulates its gene expression during transition from one developmental stage to another and in response to various environmental conditions. Parasites in a diseased host face environments which differ from the static, well adapted in vitro conditions. Parasites thus need to adapt quickly and effectively to these conditions by establishing transcriptional states which are best suited for better survival. With the discovery of natural antisense transcripts (NATs) in this parasite and considering the various proposed mechanisms by which NATs might regulate gene expression, it has been speculated that these might be playing a critical role in gene regulation. We report here the diversity of NATs in this parasite, using isolates taken directly from patients with differing clinical symptoms caused by malaria infection. Using a custom designed strand specific whole genome microarray, a total of 797 NATs targeted against annotated loci have been detected. Out of these, 545 NATs are unique to this study. The majority of NATs were positively correlated with the expression pattern of the sense transcript. However, 96 genes showed a change in sense/antisense ratio on comparison between uncomplicated and complicated disease conditions. The antisense transcripts map to a broad range of biochemical/metabolic pathways, especially pathways pertaining to the central carbon metabolism and stress related pathways. Our data strongly suggests that a large group of NATs detected here are unannotated transcription units antisense to annotated gene models. The results reveal a previously unknown set of NATs that prevails in this parasite, their differential regulation in disease conditions and mapping to functionally well annotated genes. The results detailed here call for studies to deduce the possible mechanism of action of NATs, which would further help in understanding the in vivo pathological adaptations of these parasites.

Trans-splicing and RNA editing of LSU rRNA in Diplonema mitochondria.

Mitochondrial ribosomal RNAs (rRNAs) often display reduced size and deviant secondary structure, and sometimes are fragmented, as are their corresponding genes. Here we report a mitochondrial large subunit rRNA (mt-LSU rRNA) with unprecedented features. In the protist Diplonema, the rnl gene is split into two pieces (modules 1 and 2, 534- and 352-nt long) that are encoded by distinct mitochondrial chromosomes, yet the rRNA is continuous. To reconstruct the post-transcriptional maturation pathway of this rRNA, we have catalogued transcript intermediates by deep RNA sequencing and RT-PCR. Gene modules are transcribed separately. Subsequently, transcripts are end-processed, the module-1 transcript is polyuridylated and the module-2 transcript is polyadenylated. The two modules are joined via trans-splicing that retains at the junction ∼ 26 uridines, resulting in an extent of insertion RNA editing not observed before in any system. The A-tail of trans-spliced molecules is shorter than that of mono-module 2, and completely absent from mitoribosome-associated mt-LSU rRNA. We also characterize putative antisense transcripts. Antisense-mono-modules corroborate bi-directional transcription of chromosomes. Antisense-mt-LSU rRNA, if functional, has the potential of guiding concomitantly trans-splicing and editing of this rRNA. Together, these findings open a window on the investigation of complex regulatory networks that orchestrate multiple and biochemically diverse post-transcriptional events.

Integrated detection of natural antisense transcripts using strand-specific RNA sequencing data.

Pairs of RNA molecules transcribed from partially or entirely complementary loci are called cis-natural antisense transcripts (cis-NATs), and they play key roles in the regulation of gene expression in many organisms. A promising experimental tool for profiling sense and antisense transcription is strand-specific RNA sequencing (ssRNA-seq). To identify cis-NATs using ssRNA-seq, we developed a new computational method based on a model comparison framework that incorporates the inherent variable efficiency of generating perfectly strand-specific libraries. Applying the method to new ssRNA-seq data from whole-root and cell-type-specific Arabidopsis libraries confirmed most of the known cis-NAT pairs and identified 918 additional cis-NAT pairs. Newly identified cis-NAT pairs are supported by polyadenylation data, alternative splicing patterns, and RT-PCR validation. We found 209 cis-NAT pairs that have opposite expression levels in neighboring cell types, implying cell-type-specific roles for cis-NATs. By integrating a genome-wide epigenetic profile of Arabidopsis, we identified a unique chromatin signature of cis-NATs, suggesting a connection between cis-NAT transcription and chromatin modification in plants. An analysis of small-RNA sequencing data showed that ∼4% of cis-NAT pairs produce putative cis-NAT-induced siRNAs. Taken together, our data and analyses illustrate the potential for multifaceted regulatory roles of plant cis-NATs.

Identification and characterization of a cis-encoded antisense RNA associated with the replication process of Salmonella enterica serovar Typhi.

Antisense RNAs that originate from the complementary strand of protein coding genes are involved in the regulation of gene expression in all domains of life. In bacteria, some of these antisense RNAs are transcriptional noise while others play a vital role to adapt the cell to changing environmental conditions. By deep sequencing analysis of transcriptome of Salmonella enterica serovar Typhi, a partial RNA sequence encoded in-cis to the dnaA gene was revealed. Northern blot and RACE analysis confirmed the transcription of this antisense RNA which was expressed mostly in the stationary phase of the bacterial growth and also under iron limitation and osmotic stress. Pulse expression analysis showed that overexpression of the antisense RNA resulted in a significant increase in the mRNA levels of dnaA, which will ultimately enhance their translation. Our findings have revealed that antisense RNA of dnaA is indeed transcribed not merely as a by-product of the cell's transcription machinery but plays a vital role as far as stability of dnaA mRNA is concerned.

New methods as alternative or corrective measures for the pitfalls and artifacts of reverse transcription and polymerase chain reactions (RT-PCR) in cloning chimeric or antisense-accompanied RNA.

We established new methods for cloning cDNA ends that start with reverse transcription (RT) and soon proceed with the synthesis of the second cDNA strand, avoiding manipulations of fragile RNA. Our 3'-end cloning method does not involve poly-dT primers and polymerase chain reactions (PCR), is low in efficiency but high in fidelity and can clone those RNAs without a poly-A tail. We also established a cDNA protection assay to supersede RNA protection assay. The protected cDNA can be amplified, cloned and sequenced, enhancing sensitivity and fidelity. We report that RT product using gene-specific primer (GSP) cannot be gene- or strand-specific because RNA sample contains endogenous random primers (ERP). The gene-specificity may be improved by adding a linker sequence at the 5'-end of the GSP to prime RT and using the linker as a primer in the ensuing PCR. The strand-specificity may be improved by using strand-specific DNA oligos in our protection assay. The CDK4 mRNA and TSPAN31 mRNA are transcribed from the opposite DNA strands and overlap at their 3' ends. Using this relationship as a model, we found that the overlapped sequence might serve as a primer with its antisense as the template to create a wrong-template extension in RT or PCR. We infer that two unrelated RNAs or cDNAs overlapping at the 5'- or 3'-end might create a spurious chimera in this way, and many chimeras with a homologous sequence may be such artifacts. The ERP and overlapping antisense together set complex pitfalls, which one should be aware of.

Association of RNAs with Bacillus subtilis Hfq.

The prevalence and characteristics of small regulatory RNAs (sRNAs) have not been well characterized for Bacillus subtilis, an important model system for Gram-positive bacteria. However, B. subtilis was recently found to synthesize many candidate sRNAs during stationary phase. In the current study, we performed deep sequencing on Hfq-associated RNAs and found that a small subset of sRNAs associates with Hfq, an enigmatic RNA-binding protein that stabilizes sRNAs in Gram-negatives, but whose role is largely unknown in Gram-positive bacteria. We also found that Hfq associated with antisense RNAs, antitoxin transcripts, and many mRNA leaders. Several new candidate sRNAs and mRNA leader regions were also discovered by this analysis. Additionally, mRNA fragments overlapping with start or stop codons associated with Hfq, while, in contrast, relatively few full-length mRNAs were recovered. Deletion of hfq reduced the intracellular abundance of several representative sRNAs, suggesting that B. subtilis Hfq-sRNA interactions may be functionally significant in vivo. In general, we anticipate this catalog of Hfq-associated RNAs to serve as a resource in the functional characterization of Hfq in B. subtilis.

Exosomes secreted from human colorectal cancer cell lines contain mRNAs, microRNAs and natural antisense RNAs, that can transfer into the human hepatoma HepG2 and lung cancer A549 cell lines.

Exosomes are microvesicles that are released from various cells into the extracellular space. It has been reported that the components within exosomes vary according to the type of secreted cell. In the present study, we investigated the tetraspanin family proteins CD63, CD9 and CD81 as useful collection markers of exosomes derived from the three colorectal cancer (CRC) cell lines HCT-15, SW480 and WiDr. In addition, we aimed to detect the mRNAs, microRNAs and natural antisense RNAs within the exosomes secreted from the three CRC cell lines. Furthermore, we examined whether exosomes containing their RNAs were transferred into the hepatoma cell line HepG2 and lung cancer cell line A549. CD81 was detected in exosomes secreted from the three CRC cell lines. This result indicates that CD81 can be a collection marker of exosomes derived from the three CRC cell lines. When the RNA species within exosomes derived from the three CRC cell lines were examined, the mRNAs of housekeeping genes such as ACTB and GAPDH, the microRNAs such as miR-21, miR-192 and miR-221, and the natural antisense RNAs of LRRC24, MDM2 and CDKN1A genes, were detected. We discovered their natural antisense RNAs within exosomes for the first time in the present study. Furthermore, PKH67-labeled exosomes derived from the CRC cell lines were taken up into HepG2 and A549 cells. These findings indicate that the intracellular RNAs enclosed within exosomes are secreted to the outside, and exosomes derived from the CRC cell lines are transferred into HepG2 and A549 cells. In conclusion, we reveal that exosomes derived from the CRC cell lines contain mRNAs, microRNAs and natural antisense RNAs, and can be delivered into HepG2 and A549 cells. These findings indicate that exosomal RNAs can shuttle between cells, and may be involved in the regulation of gene expression in recipient cells.

Directional RNA deep sequencing sheds new light on the transcriptional response of Anabaena sp. strain PCC 7120 to combined-nitrogen deprivation.

BACKGROUND: Cyanobacteria are potential sources of renewable chemicals and biofuels and serve as model organisms for bacterial photosynthesis, nitrogen fixation, and responses to environmental changes. Anabaena (Nostoc) sp. strain PCC 7120 (hereafter Anabaena) is a multicellular filamentous cyanobacterium that can "fix" atmospheric nitrogen into ammonia when grown in the absence of a source of combined nitrogen. Because the nitrogenase enzyme is oxygen sensitive, Anabaena forms specialized cells called heterocysts that create a microoxic environment for nitrogen fixation. We have employed directional RNA-seq to map the Anabaena transcriptome during vegetative cell growth and in response to combined-nitrogen deprivation, which induces filaments to undergo heterocyst development. Our data provide an unprecedented view of transcriptional changes in Anabaena filaments during the induction of heterocyst development and transition to diazotrophic growth. RESULTS: Using the Illumina short read platform and a directional RNA-seq protocol, we obtained deep sequencing data for RNA extracted from filaments at 0, 6, 12, and 21 hours after the removal of combined nitrogen. The RNA-seq data provided information on transcript abundance and boundaries for the entire transcriptome. From these data, we detected novel antisense transcripts within the UTRs (untranslated regions) and coding regions of key genes involved in heterocyst development, suggesting that antisense RNAs may be important regulators of the nitrogen response. In addition, many 5' UTRs were longer than anticipated, sometimes extending into upstream open reading frames (ORFs), and operons often showed complex structure and regulation. Finally, many genes that had not been previously identified as being involved in heterocyst development showed regulation, providing new candidates for future studies in this model organism. CONCLUSIONS: Directional RNA-seq data were obtained that provide comprehensive mapping of transcript boundaries and abundance for all transcribed RNAs in Anabaena filaments during the response to nitrogen deprivation. We have identified genes and noncoding RNAs that are transcriptionally regulated during heterocyst development. These data provide detailed information on the Anabaena transcriptome as filaments undergo heterocyst development and begin nitrogen fixation.

Complete characterization of the edited transcriptome of the mitochondrion of Physarum polycephalum using deep sequencing of RNA.

RNAs transcribed from the mitochondrial genome of Physarum polycephalum are heavily edited. The most prevalent editing event is the insertion of single Cs, with Us and dinucleotides also added at specific sites. The existence of insertional editing makes gene identification difficult and localization of editing sites has relied upon characterization of individual cDNAs. We have now determined the complete mitochondrial transcriptome of Physarum using Illumina deep sequencing of purified mitochondrial RNA. We report the first instances of A and G insertions and sites of partial and extragenic editing in Physarum mitochondrial RNAs, as well as an additional 772 C, U and dinucleotide insertions. The notable lack of antisense RNAs in our non-size selected, directional library argues strongly against an RNA-guided editing mechanism. Also of interest are our findings that sites of C to U changes are unedited at a significantly higher frequency than insertional editing sites and that substitutional editing of neighboring sites appears to be coupled. Finally, in addition to the characterization of RNAs from 17 predicted genes, our data identified nine new mitochondrial genes, four of which encode proteins that do not resemble other proteins in the database. Curiously, one of the latter mRNAs contains no editing sites.

Novel retrotransposed imprinted locus identified at human 6p25.

Differentially methylated regions (DMRs) are stable epigenetic features within or in proximity to imprinted genes. We used this feature to identify candidate human imprinted loci by quantitative DNA methylation analysis. We discovered a unique DMR at the 5'-end of FAM50B at 6p25.2. We determined that sense transcripts originating from the FAM50B locus are expressed from the paternal allele in all human tissues investigated except for ovary, in which expression is biallelic. Furthermore, an antisense transcript, FAM50B-AS, was identified to be monoallelically expressed from the paternal allele in a variety of tissues. Comparative phylogenetic analysis showed that FAM50B orthologs are absent in chicken and platypus, but are present and biallelically expressed in opossum and mouse. These findings indicate that FAM50B originated in Therians after divergence from Prototherians via retrotransposition of a gene on the X chromosome. Moreover, our data are consistent with acquisition of imprinting during Eutherian evolution after divergence of Glires from the Euarchonta mammals. FAM50B expression is deregulated in testicular germ cell tumors, and loss of imprinting occurs frequently in testicular seminomas, suggesting an important role for FAM50B in spermatogenesis and tumorigenesis. These results also underscore the importance of accounting for parental origin in understanding the mechanism of 6p25-related diseases.

A simple and efficient method to search for selected primary transcripts: non-coding and antisense RNAs in the human pathogen Enterococcus faecalis.

Enterococcus faecalis is a commensal bacterium and a major opportunistic human pathogen. In this study, we combined in silico predictions with a novel 5'RACE-derivative method coined '5'tagRACE', to perform the first search for non-coding RNAs (ncRNAs) encoded on the E. faecalis chromosome. We used the 5'tagRACE to simultaneously probe and characterize primary transcripts, and demonstrate here the simplicity, the reliability and the sensitivity of the method. The 5'tagRACE is complementary to tiling arrays or RNA-sequencing methods, and is also directly applicable to deep RNA sequencing and should significantly improve functional studies of bacterial RNA landscapes. From 45 selected loci of the E. faecalis chromosome, we discovered and mapped 29 novel ncRNAs, 10 putative novel mRNAs and 16 antisense transcriptional organizations. We describe in more detail the oxygen-dependent expression of one ncRNA located in an E. faecalis pathogenicity island, the existence of an ncRNA that is antisense to the ncRNA modulator of the RNA polymerase, SsrS and provide evidences for the functional interplay between two distinct toxin-antitoxin modules.

Strand-specific transcriptome profiling with directly labeled RNA on genomic tiling microarrays.

BACKGROUND: With lower manufacturing cost, high spot density, and flexible probe design, genomic tiling microarrays are ideal for comprehensive transcriptome studies. Typically, transcriptome profiling using microarrays involves reverse transcription, which converts RNA to cDNA. The cDNA is then labeled and hybridized to the probes on the arrays, thus the RNA signals are detected indirectly. Reverse transcription is known to generate artifactual cDNA, in particular the synthesis of second-strand cDNA, leading to false discovery of antisense RNA. To address this issue, we have developed an effective method using RNA that is directly labeled, thus by-passing the cDNA generation. This paper describes this method and its application to the mapping of transcriptome profiles. RESULTS: RNA extracted from laboratory cultures of Porphyromonas gingivalis was fluorescently labeled with an alkylation reagent and hybridized directly to probes on genomic tiling microarrays specifically designed for this periodontal pathogen. The generated transcriptome profile was strand-specific and produced signals close to background level in most antisense regions of the genome. In contrast, high levels of signal were detected in the antisense regions when the hybridization was done with cDNA. Five antisense areas were tested with independent strand-specific RT-PCR and none to negligible amplification was detected, indicating that the strong antisense cDNA signals were experimental artifacts. CONCLUSIONS: An efficient method was developed for mapping transcriptome profiles specific to both coding strands of a bacterial genome. This method chemically labels and uses extracted RNA directly in microarray hybridization. The generated transcriptome profile was free of cDNA artifactual signals. In addition, this method requires fewer processing steps and is potentially more sensitive in detecting small amount of RNA compared to conventional end-labeling methods due to the incorporation of more fluorescent molecules per RNA fragment.

Involvement of argonaute proteins in gene silencing and activation by RNAs complementary to a non-coding transcript at the progesterone receptor promoter.

Double-stranded RNAs that are complementary to non-coding transcripts at gene promoters can activate or inhibit gene expression in mammalian cells. Understanding the mechanism for modulating gene expression by promoter-targeted antigene RNAs (agRNAs) will require identification of the proteins involved in recognition. Previous reports have implicated argonaute (AGO) proteins, but identifications have differed with involvement of AGO1, AGO2, or both AGO1 and AGO2 being reported by different studies. The roles of AGO3 and AGO4 have not been investigated. Here, we examine the role of AGO 1-4 in gene silencing and activation of the progesterone receptor (PR) gene. Expression of AGO2 is necessary for efficient gene silencing or activation and AGO2 is recruited to the non-coding transcript that overlaps the promoter during both gene silencing and activation. Expression of AGO1, AGO3 and AGO4 are not necessary for gene silencing or activation nor are AGO1, AGO3, or AGO4 recruited to the target non-coding transcript during gene activation. These data indicate that AGO2 is the primary AGO variant involved in modulating expression of PR by agRNAs.

Strand-specific RNA sequencing reveals extensive regulated long antisense transcripts that are conserved across yeast species.

BACKGROUND: Recent studies in budding yeast have shown that antisense transcription occurs at many loci. However, the functional role of antisense transcripts has been demonstrated only in a few cases and it has been suggested that most antisense transcripts may result from promiscuous bi-directional transcription in a dense genome. RESULTS: Here, we use strand-specific RNA sequencing to study anti-sense transcription in Saccharomyces cerevisiae. We detect 1,103 putative antisense transcripts expressed in mid-log phase growth, ranging from 39 short transcripts covering only the 3' UTR of sense genes to 145 long transcripts covering the entire sense open reading frame. Many of these antisense transcripts overlap sense genes that are repressed in mid-log phase and are important in stationary phase, stress response, or meiosis. We validate the differential regulation of 67 antisense transcripts and their sense targets in relevant conditions, including nutrient limitation and environmental stresses. Moreover, we show that several antisense transcripts and, in some cases, their differential expression have been conserved across five species of yeast spanning 150 million years of evolution. Divergence in the regulation of antisense transcripts to two respiratory genes coincides with the evolution of respiro-fermentation. CONCLUSIONS: Our work provides support for a global and conserved role for antisense transcription in yeast gene regulation.

Identification of regulatory RNAs in Bacillus subtilis.

Post-transcriptional regulatory mechanisms are widespread in bacteria. Interestingly, current published data hint that some of these mechanisms may be non-random with respect to their phylogenetic distribution. Although small, trans-acting regulatory RNAs commonly occur in bacterial genomes, they have been better characterized in Gram-negative bacteria, leaving the impression that they may be less important for Firmicutes. It has been presumed that Gram-positive bacteria, in particular the Firmicutes, are likely to utilize cis-acting regulatory RNAs located within the 5' mRNA leader region more often than trans-acting regulatory RNAs. In this analysis we catalog, by a deep sequencing-based approach, both classes of regulatory RNA candidates for Bacillus subtilis, the model microorganism for Firmicutes. We successfully recover most of the known small RNA regulators while also identifying a greater number of new candidate RNAs. We anticipate these data to be a broadly useful resource for analysis of post-transcriptional regulatory strategies in B. subtilis and other Firmicutes.