Analysis of wheat SAGE tags reveals evidence for widespread antisense transcription

Antisense Transcription in Plants: A Systematic Review and an Update on cis-NATs of Sugarcane

Initially, natural antisense transcripts (NATs, natRNAs, or asRNAs) were considered repressors; however, their functions in gene regulation are diverse. Positive, negative, or neutral correlations to the cognate gene expression have been noted. Although the first studies were published about 50 years ago, there is still much to be investigated regarding antisense transcripts in plants. A systematic review of scientific publications available in the Web of Science databases was conducted to contextualize how the studying of antisense transcripts has been addressed. Studies were classified considering three categories: “Natural antisense” (208), artificial antisense used in “Genetic Engineering” (797), or “Natural antisense and Genetic Engineering”-related publications (96). A similar string was used for a systematic search in the NCBI Gene database. Of the 1132 antisense sequences found for plants, only 0.8% were cited in PubMed and had antisense information confirmed. This value was the lowest when compared to fungi (2.9%), bacteria (2.3%), and mice (54.1%). Finally, we present an update for the cis-NATs identified in Saccharum spp. Of the 1413 antisense transcripts found in different experiments, 25 showed concordant expressions, 22 were discordant, 1264 did not correlate with the cognate genes, and 102 presented variable results depending on the experiment.

Independent origin of MIRNA genes controlling homologous target genes by partial inverted duplication of antisense-transcribed sequences

Some microRNAs (miRNAs) are key regulators of developmental processes, mainly by controlling the accumulation of transcripts encoding transcription factors that are important for morphogenesis. MADS-box genes encode a family of transcription factors which control diverse developmental processes in flowering plants. Here we study the convergent evolution of two MIRNA (MIR) gene families, named MIR444 and MIR824, targeting members of the same clade of MIKC^C -group MADS-box genes. We show that these two MIR genes most likely originated independently in monocots (MIR444) and in Brassicales (eudicots, MIR824). We provide evidence that, in both cases, the future target gene was transcribed in antisense prior to the evolution of the MIR genes. Both MIR genes then likely originated by a partial inverted duplication of their target genes, resulting in natural antisense organization of the newly evolved MIR gene and its target gene at birth. We thus propose a model for the origin of MIR genes, MEPIDAS (MicroRNA Evolution by Partial Inverted Duplication of Antisense-transcribed Sequences). MEPIDAS is a refinement of the inverted duplication hypothesis. According to MEPIDAS, a MIR gene evolves at a genomic locus at which the future target gene is also transcribed in the antisense direction. A partial inverted duplication at this locus causes the antisense transcript to fold into a stem-loop structure that is recognized by the miRNA biogenesis machinery to produce a miRNA that regulates the gene at this locus. Our analyses exemplify how to elucidate the origin of conserved miRNAs by comparative genomics and will guide future studies. OPEN RESEARCH BADGE: This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://www.ncbi.nlm.nih.gov/genbank/.

An Overview of Methodologies in Studying lncRNAs in the High-Throughput Era: When Acronyms ATTACK!

The discovery of pervasive transcription in eukaryotic genomes provided one of many surprising (and perhaps most surprising) findings of the genomic era and led to the uncovering of a large number of previously unstudied transcriptional events. This pervasive transcription leads to the production of large numbers of noncoding RNAs (ncRNAs) and thus opened the window to study these diverse, abundant transcripts of unclear relevance and unknown function. Since that discovery, recent advances in high-throughput sequencing technologies have identified a large collection of ncRNAs, from microRNAs to long noncoding RNAs (lncRNAs). Subsequent discoveries have shown that many lncRNAs play important roles in various eukaryotic processes; these discoveries have profoundly altered our understanding of the regulation of eukaryotic gene expression. Although the identification of ncRNAs has become a standard experimental approach, the functional characterization of these diverse ncRNAs remains a major challenge. In this chapter, we highlight recent progress in the methods to identify lncRNAs and the techniques to study the molecular function of these lncRNAs and the application of these techniques to the study of plant lncRNAs.

Functional and DNA-protein binding studies of WRKY transcription factors and their expression analysis in response to biotic and abiotic stress in wheat (Triticum aestivum L.)

WRKY, a plant-specific transcription factor family, plays vital roles in pathogen defense, abiotic stress, and phytohormone signalling. Little is known about the roles and function of WRKY transcription factors in response to rust diseases in wheat. In the present study, three TaWRKY genes encoding complete protein sequences were cloned. They belonged to class II and III WRKY based on the number of WRKY domains and the pattern of zinc finger structures. Twenty-two DNA-protein binding docking complexes predicted stable interactions of WRKY domain with W-box. Quantitative real-time-PCR using wheat near-isogenic lines with or without Lr28 gene revealed differential up- or down-regulation in response to biotic and abiotic stress treatments which could be responsible for their functional divergence in wheat. TaWRKY62 was found to be induced upon treatment with JA, MJ, and SA and reduced after ABA treatments. Maximum induction of six out of seven genes occurred at 48 h post inoculation due to pathogen inoculation. Hence, TaWRKY (49, 50, 52, 55, 57, and 62) can be considered as potential candidate genes for further functional validation as well as for crop improvement programs for stress resistance. The results of the present study will enhance knowledge towards understanding the molecular basis of mode of action of WRKY transcription factor genes in wheat and their role during leaf rust pathogenesis in particular.

Deep transcriptome sequencing provides new insights into the structural and functional organization of the wheat genome

Background

Because of its size, allohexaploid nature, and high repeat content, the bread wheat genome is a good model to study the impact of the genome structure on gene organization, function, and regulation. However, because of the lack of a reference genome sequence, such studies have long been hampered and our knowledge of the wheat gene space is still limited. The access to the reference sequence of the wheat chromosome 3B provided us with an opportunity to study the wheat transcriptome and its relationships to genome and gene structure at a level that has never been reached before.

Results

By combining this sequence with RNA-seq data, we construct a fine transcriptome map of the chromosome 3B. More than 8,800 transcription sites are identified, that are distributed throughout the entire chromosome. Expression level, expression breadth, alternative splicing as well as several structural features of genes, including transcript length, number of exons, and cumulative intron length are investigated. Our analysis reveals a non-monotonic relationship between gene expression and structure and leads to the hypothesis that gene structure is determined by its function, whereas gene expression is subject to energetic cost. Moreover, we observe a recombination-based partitioning at the gene structure and function level.

Conclusions

Our analysis provides new insights into the relationships between gene and genome structure and function. It reveals mechanisms conserved with other plant species as well as superimposed evolutionary forces that shaped the wheat gene space, likely participating in wheat adaptation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0601-9) contains supplementary material, which is available to authorized users.

EXPRSS: an Illumina based high-throughput expression-profiling method to reveal transcriptional dynamics

Background

Next Generation Sequencing technologies have facilitated differential gene expression analysis through RNA-seq and Tag-seq methods. RNA-seq has biases associated with transcript lengths, lacks uniform coverage of regions in mRNA and requires 10–20 times more reads than a typical Tag-seq. Most existing Tag-seq methods either have biases or not high throughput due to use of restriction enzymes or enzymatic manipulation of 5’ ends of mRNA or use of RNA ligations.

Results

We have developed EXpression Profiling through Randomly Sheared cDNA tag Sequencing (EXPRSS) that employs acoustic waves to randomly shear cDNA and generate sequence tags at a relatively defined position (~150-200 bp) from the 3′ end of each mRNA. Implementation of the method was verified through comparative analysis of expression data generated from EXPRSS, NlaIII-DGE and Affymetrix microarray and through qPCR quantification of selected genes. EXPRSS is a strand specific and restriction enzyme independent tag sequencing method that does not require cDNA length-based data transformations. EXPRSS is highly reproducible, is high-throughput and it also reveals alternative polyadenylation and polyadenylated antisense transcripts. It is cost-effective using barcoded multiplexing, avoids the biases of existing SAGE and derivative methods and can reveal polyadenylation position from paired-end sequencing.

Conclusions

EXPRSS Tag-seq provides sensitive and reliable gene expression data and enables high-throughput expression profiling with relatively simple downstream analysis.

Electronic supplementary material

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

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.

Prediction and identification of natural antisense transcripts and their small RNAs in soybean (Glycine max)

BACKGROUND

Natural antisense transcripts (NATs) are a class of RNAs that contain a sequence complementary to other transcripts. NATs occur widely in eukaryotes and play critical roles in post-transcriptional regulation. Soybean NAT sequences are predicted in the PlantNATsDB, but detailed analyses of these NATs remain to be performed.

RESULTS

A total of 26,216 NATs, including 994 cis-NATs and 25,222 trans-NATs, were predicted in soybean. Each sense transcript had 1-177 antisense transcripts. We identified 21 trans-NATs using RT-PCR amplification. Additionally, we identified 179 cis-NATs and 6,629 trans-NATs that gave rise to small RNAs; these were enriched in the NAT overlapping region. The most abundant small RNAs were 21, 22, and 24 nt in length. The generation of small RNAs was biased to one stand of the NATs, and the degradation of NATs was biased. High-throughput sequencing of the degradome allowed for the global identification of NAT small interfering RNAs (nat-siRNAs) targets. 446 target genes for 165 of these nat-siRNAs were identified. The nat-siRNA target could be one transcript of a given NAT, or from other gene transcripts. We identified five NAT transcripts containing a hairpin structure that is characteristic of pre-miRNA. We identified a total of 86 microRNA (miRNA) targets that had antisense transcripts in soybean.

CONCLUSIONS

We globally identified nat-siRNAs, and the targets of nat-siRNAs in soybean. It is likely that the cis-NATs, trans-NATs, nat-siRNAs, miRNAs, and miRNA targets form complex regulatory networks.

Two modified RNA extraction methods compatible with transcript profiling and gene expression analysis for cotton roots

Efficient isolation of high-quality RNA is of prime importance for optimal transcript profiling results and further gene expression analysis. However, it is difficult for cotton roots because of lower-than-average RNA content and high content of polysaccharides, polyphenols, and other secondary metabolites. To develop simple and reliable protocols for high-quality RNA extraction from cotton roots for transcript profiling and gene expression analysis, some modifications were introduced to a reported plant RNA isolation protocol and a reagent kit method. Using method A, we successfully extracted high-quality RNA for transcript profiling from cotton roots. Gel electrophoresis analysis and polymerase chain reaction (PCR) assay indicated that RNA had good integrity without protein and genomic DNA contamination. Furthermore, the A260/280 (1.9) and A260/230 (1.6) ratios indicated that the isolated RNA was of high purity. Using method B, about 7 µg total RNA of high quality could be obtained from 0.1 g samples from cotton roots, which can be used for reverse-transcription (RT)-PCR and quantitative real-time RT-PCR. The two RNA extraction methods were used to investigate different gene expression of cotton roots (Gossypium hirsutum) infected by weak pathogenic Verticillium dahliae and the results showed they can satisfy the transcript profiling and quantitative real-time RT-PCR requirements for RNA. Supplemental materials are available for this article. Go to the publisher's online edition of Preparative Biochemistry and Biotechnology to view the supplemental file.

SOLiD-SAGE of endophyte-infected red fescue reveals numerous effects on host transcriptome and an abundance of highly expressed fungal secreted proteins

One of the most important plant-fungal symbiotic relationships is that of cool season grasses with endophytic fungi of the genera Epichloë and Neotyphodium. These associations often confer benefits, such as resistance to herbivores and improved drought tolerance, to the hosts. One benefit that appears to be unique to fine fescue grasses is disease resistance. As a first step towards understanding the basis of the endophyte-mediated disease resistance in Festuca rubra we carried out a SOLiD-SAGE quantitative transcriptome comparison of endophyte-free and Epichloë festucae-infected F. rubra. Over 200 plant genes involved in a wide variety of physiological processes were statistically significantly differentially expressed between the two samples. Many of the endophyte expressed genes were surprisingly abundant, with the most abundant fungal tag representing over 10% of the fungal mapped tags. Many of the abundant fungal tags were for secreted proteins. The second most abundantly expressed fungal gene was for a secreted antifungal protein and is of particular interest regarding the endophyte-mediated disease resistance. Similar genes in Penicillium and Aspergillus spp. have been demonstrated to have antifungal activity. Of the 10 epichloae whole genome sequences available, only one isolate of E. festucae and Neotyphodium gansuense var inebrians have an antifungal protein gene. The uniqueness of this gene in E. festucae from F. rubra, its transcript abundance, and the secreted nature of the protein, all suggest it may be involved in the disease resistance conferred to the host, which is a unique feature of the fine fescue-endophyte symbiosis.

A genome-wide investigation of expression characteristics of natural antisense transcripts in liver and muscle samples of pigs

Natural antisense transcripts are endogenous transcripts that are complementary to the sense-strand of DNA. These transcripts have been identified in various eukaryotic species and are involved in a broad range of regulatory events and biological processes. However, their general biological functions, expression characteristics and regulatory mechanisms are still unclear. In this study, 497 liver and 586 muscle samples were harvested from a White Duroc×Erhualian F₂ resource population. The expression profiles of sense and antisense transcripts were determined by tag-based RNA sequencing. We identified 33.7% and 20.4% of transcripts having both sense and antisense expression, and 12.5% and 6.1% of transcripts only expressing antisense transcripts in liver and muscle, respectively. More than 32.2% of imprinting or predicted imprinting genes in the geneimprint database were detected with both sense and antisense expression. The correlations between sense and antisense expression in sense-antisense pairs were diverse in both liver and muscle, showing positive, negative or absent correlation. Antisense expression increases gene expression variability. More interestingly, compared to eQTL mapping of sense transcripts in which more than one eQTL was mapped for a transcript, only one eQTL was identified for each antisense transcript, and the percentage of cis-eQTL in antisense eQTL was higher than that in sense eQTL. This suggests that the expressions of antisense transcripts tend to be cis-regulated by a single genomic locus. To our knowledge, this study is the first systematical investigation of antisense transcription in pigs. The findings improve our understanding of the complexity of porcine transcriptome.

Flowering time control: another window to the connection between antisense RNA and chromatin

A high proportion of all eukaryotic genes express antisense RNA (asRNA), which accumulates to varying degrees at different loci. Whether there is a general function for asRNA is unknown, but its widespread occurrence and frequent regulation by stress suggest an important role. The best-characterized plant gene exhibiting a complex antisense transcript pattern is the Arabidopsis floral regulator FLOWERING LOCUS C (FLC). Changes occur in the accumulation, splicing, and polyadenylation of this antisense transcript, termed COOLAIR, in different environments and genotypes. These changes are associated with altered chromatin regulation and differential FLC expression, provoking mechanistic comparisons with many well-studied loci in yeast and mammals. Detailed analysis of these specific examples may shed light on the complex interplay between asRNA and chromatin modifications in different genomes.

Transcript profiling during the early development of the maize brace root via Solexa sequencing

Currently, the molecular regulation mechanisms involved in the early development of maize brace root are poorly known. To gain insight into the transcriptome dynamics that are associated with its development, genome-wide gene expression profiling was conducted by Solexa sequencing (Illumina Inc., San Diego, CA, USA). More than five million tags were generated from the stem node tissues without and with just-emerged brace roots, including 149,524 and 178,131 clean tags in the two libraries, respectively. Of these, 82,864 (55.4%) and 91,678 (51.5%) tags were matched to the reference genes. The most differentially expressed tags with a log(2) ratio > 2 or < -2 (P < 0.001) were analyzed further, representing 143 up-regulated and 152 down-regulated genes, except for unknown transcripts, which were classified into 11 functional categories. The most enriched categories were those of metabolism, signal transduction and cellular transport. Many genes or biological pathways were found to be commonly shared between brace root and lateral or adventitious root development, such as genes participating in cell wall degradation and synthesis, auxin transport and signaling, ethylene signaling, etc. Next, the expression patterns of 20 genes were assessed by quantitative real-time PCR, and the results obtained showed general agreement with the Solexa analysis. Furthermore, a comparison of the brace root transcriptome with that of maize primary root revealed substantial differences in the categories and abundances of expressed transcripts. In conclusion, we first reveal the complex changes in the transcriptome during the early development of maize brace root and provide a comprehensive set of data that are essential for understanding its molecular regulation.

Haploids in flowering plants: origins and exploitation

The first haploid angiosperm, a dwarf form of cotton with half the normal chromosome complement, was discovered in 1920, and in the ninety years since then such plants have been identified in many other species. They can occur either spontaneously or can be induced by modified pollination methods in vivo, or by in vitro culture of immature male or female gametophytes. Haploids represent an immediate, one-stage route to homozygous diploids and thence to F(1) hybrid production. The commercial exploitation of heterosis in such F(1) hybrids leads to the development of hybrid seed companies and subsequently to the GM revolution in agriculture. This review describes the range of techniques available for the isolation or induction of haploids and discusses their value in a range of areas, from fundamental research on mutant isolation and transformation, through to applied aspects of quantitative genetics and plant breeding. It will also focus on how molecular methods have been used recently to explore some of the underlying aspects of this fascinating developmental phenomenon.

Gene expression analysis in cadmium-stressed roots of a low cadmium-accumulating solanaceous plant, Solanum torvum

Solanum torvum Sw. cv. Torubamubiga (TB) is a low cadmium (Cd)-accumulating plant. To elucidate the molecular mechanisms of the Cd acclimation process in TB roots, transcriptional regulation was analysed in response to mild Cd treatment: 0.1 muM CdCl(2) in hydroponic solution. A unigene set consisting of 6296 unigene sequences was constructed from 18 816 TB cDNAs. The distribution of functional categories was similar to tomato, while 330 unigenes were suggested to be TB specific. For expression profiling, the SuperSAGE method was adapted for use with Illumina sequencing technology. Expression tag libraries were constructed from Cd-treated (for 3 h, 1 d, and 3 d) and untreated roots, and 34 269 species of independent tags were collected. Moreover, 6237 tags were ascribed to the TB or eggplant (aubergine) unigene sequences. Time-course changes were examined, and 2049 up- and 2022 down-regulated tags were identified. Although no tags annotated to metal transporter genes were significantly regulated, a tag annotated to AtFRD3, a xylem-loading citrate transporter, was down-regulated. In addition to induction of heavy metal chaperone proteins, antioxidative and sulphur-assimilating enzymes were induced, confirming that oxidative stress developed even using a mild Cd concentration. Rapid repression of dehydration-related transcription factors and aquaporin isoforms suggests that dehydration stress is a potential constituent of Cd-induced biochemical impediments. These transcriptional changes were also confirmed by real-time reverse transcription-PCR. Further additions of TB unigene sequences and functional analysis of the regulated tags will reveal the molecular basis of the Cd acclimation process, including the low Cd-accumulating characteristics of TB.

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

Background

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

Results

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

Conclusion

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

Transcriptomic analysis of starch biosynthesis in the developing grain of hexaploid wheat

The expression of genes involved in starch synthesis in wheat was analyzed together with the accumulation profiles of soluble sugars, starch, protein, and starch granule distribution in developing caryopses obtained from the same biological materials used for profiling of gene expression using DNA microarrays. Multiple expression patterns were detected for the different starch biosynthetic gene isoforms, suggesting their relative importance through caryopsis development. Members of the ADP-glucose pyrophosphorylase, starch synthase, starch branching enzyme, and sucrose synthase gene families showed different expression profiles; expression of some members of these gene families coincided with a period of high accumulation of starch while others did not. A biphasic pattern was observed in the rates of starch and protein accumulation which paralleled changes in global gene expression. Metabolic and regulatory genes that show a pattern of expression similar to starch accumulation and granule size distribution were identified, suggesting their coinvolvement in these biological processes.