Genome-wide transcription analyses in rice using tiling microarrays

RNA Pol II Assembly Affects ncRNA Expression

RNA pol II assembly occurs in the cytoplasm before translocation of the enzyme to the nucleus. Affecting this assembly influences mRNA transcription in the nucleus and mRNA decay in the cytoplasm. However, very little is known about the consequences on ncRNA synthesis. In this work, we show that impairment of RNA pol II assembly leads to a decrease in cryptic non-coding RNAs (preferentially CUTs and SUTs). This alteration is partially restored upon overcoming the assembly defect. Notably, this drop in ncRNAs is only partially dependent on the nuclear exosome, which suggests a major specific effect of enzyme assembly. Our data also point out a defect in transcription termination, which leads us to propose that CTD phosphatase Rtr1 could be involved in this process.

Linking discoveries, mechanisms, and technologies to develop a clearer perspective on plant long noncoding RNAs

Long noncoding RNAs (lncRNAs) are a large and diverse class of genes in eukaryotic genomes that contribute to a variety of regulatory processes. Functionally characterized lncRNAs play critical roles in plants, ranging from regulating flowering to controlling lateral root formation. However, findings from the past decade have revealed that thousands of lncRNAs are present in plant transcriptomes, and characterization has lagged far behind identification. In this setting, distinguishing function from noise is challenging. However, the plant community has been at the forefront of discovery in lncRNA biology, providing many functional and mechanistic insights that have increased our understanding of this gene class. In this review, we examine the key discoveries and insights made in plant lncRNA biology over the past two and a half decades. We describe how discoveries made in the pregenomics era have informed efforts to identify and functionally characterize lncRNAs in the subsequent decades. We provide an overview of the functional archetypes into which characterized plant lncRNAs fit and speculate on new avenues of research that may uncover yet more archetypes. Finally, this review discusses the challenges facing the field and some exciting new molecular and computational approaches that may help inform lncRNA comparative and functional analyses.

Twenty years of rice genomics research: From sequencing and functional genomics to quantitative genomics

Since the completion of the rice genome sequencing project in 2005, we have entered the era of rice genomics, which is still in its ascendancy. Rice genomics studies can be classified into three stages: structural genomics, functional genomics, and quantitative genomics. Structural genomics refers primarily to genome sequencing for the construction of a complete map of rice genome sequence. This is fundamental for rice genetics and molecular biology research. Functional genomics aims to decode the functions of rice genes. Quantitative genomics is large-scale sequence- and statistics-based research to define the quantitative traits and genetic features of rice populations. Rice genomics has been a transformative influence on rice biological research and contributes significantly to rice breeding, making rice a good model plant for studying crop sciences.

Targeting epigenetics and lncRNAs in liver disease: From mechanisms to therapeutics

Early onset and progression of liver diseases can be driven by aberrant transcriptional regulation. Different transcriptional regulation processes, such as RNA/DNA methylation, histone modification, and ncRNA-mediated targeting, can regulate biological processes in healthy cells, as well also under various pathological conditions, especially liver disease. Numerous studies over the past decades have demonstrated that liver disease has a strong epigenetic component. Therefore, the epigenetic basis of liver disease has challenged our knowledge of epigenetics, and epigenetics field has undergone an important transformation: from a biological phenomenon to an emerging focus of disease research. Furthermore, inhibitors of different epigenetic regulators, such as m⁶A-related factors, are being explored as potential candidates for preventing and treating liver diseases. In the present review, we summarize and discuss the current knowledge of five distinct but interconnected and interdependent epigenetic processes in the context of hepatic diseases: RNA methylation, DNA methylation, histone methylation, miRNAs, and lncRNAs. Finally, we discuss the potential therapeutic implications and future challenges and ongoing research in the field. Our review also provides a perspective for identifying therapeutic targets and new hepatic biomarkers of liver disease, bringing precision research and disease therapy to the modern era of epigenetics.

Crucial Cell Signaling Compounds Crosstalk and Integrative Multi-Omics Techniques for Salinity Stress Tolerance in Plants

In the era of rapid climate change, abiotic stresses are the primary cause for yield gap in major agricultural crops. Among them, salinity is considered a calamitous stress due to its global distribution and consequences. Salinity affects plant processes and growth by imposing osmotic stress and destroys ionic and redox signaling. It also affects phytohormone homeostasis, which leads to oxidative stress and eventually imbalances metabolic activity. In this situation, signaling compound crosstalk such as gasotransmitters [nitric oxide (NO), hydrogen sulfide (H2S), hydrogen peroxide (H2O2), calcium (Ca), reactive oxygen species (ROS)] and plant growth regulators (auxin, ethylene, abscisic acid, and salicylic acid) have a decisive role in regulating plant stress signaling and administer unfavorable circumstances including salinity stress. Moreover, recent significant progress in omics techniques (transcriptomics, genomics, proteomics, and metabolomics) have helped to reinforce the deep understanding of molecular insight in multiple stress tolerance. Currently, there is very little information on gasotransmitters and plant growth regulator crosstalk and inadequacy of information regarding the integration of multi-omics technology during salinity stress. Therefore, there is an urgent need to understand the crucial cell signaling crosstalk mechanisms and integrative multi-omics techniques to provide a more direct approach for salinity stress tolerance. To address the above-mentioned words, this review covers the common mechanisms of signaling compounds and role of different signaling crosstalk under salinity stress tolerance. Thereafter, we mention the integration of different omics technology and compile recent information with respect to salinity stress tolerance.

Comparative analysis of the transcriptomes of two rice subspecies during domestication

Two subspecies of rice, Oryza sativa ssp. indica and O. sativa ssp. japonica, with reproductive isolation and differences in morphology and phenotypic differences, were established during the process of rice domestication. To understand how domestication has changed the transcriptomes of the two rice subspecies and given rise to the phenotypic differences, we obtained approximately 700 Gb RNA-Seq data from 26 indica and 25 japonica accessions, and identified 97,005 transcribed fragments and 4579 novel transcriptionally active regions. The two rice subspecies had significantly different gene expression profiles, we identified 1,357 (3.3% in all genes) differentially expressed genes (DEGs) between indica and japonica rice. Combining existing gene function studies, it is found that some of these differential genes are related to the differentiation of the two subspecies, such as grain shape and cold tolerance, etc. Functional annotation of these DEGs indicates that they are involved in cell wall biosynthesis and reproductive processes. Furthermore, compared with the non-DEGs, the DEGs from both subspecies had more 5'flanking regions with low polymorphism to divergence ratios, indicating a stronger positive selection pressure on the regulation of the DEGs. This study improves our understanding of the rice genome by comparatively analyzing the transcriptomes of indica and japonica rice and identifies DEGs those may be responsible for the reproductive isolation and phenotypic differences between the two rice subspecies.

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/.

The Role of Noncoding RNAs in Double-Strand Break Repair

Genome stability is constantly threatened by DNA lesions generated by different environmental factors as well as endogenous processes. If not properly and timely repaired, damaged DNA can lead to mutations or chromosomal rearrangements, well-known reasons for genetic diseases or cancer in mammals, or growth abnormalities and/or sterility in plants. To prevent deleterious consequences of DNA damage, a sophisticated system termed DNA damage response (DDR) detects DNA lesions and initiates DNA repair processes. In addition to many well-studied canonical proteins involved in this process, noncoding RNA (ncRNA) molecules have recently been discovered as important regulators of the DDR pathway, extending the broad functional repertoire of ncRNAs to the maintenance of genome stability. These ncRNAs are mainly connected with double-strand breaks (DSBs), the most dangerous type of DNA lesions. The possibility to intentionally generate site-specific DSBs in the genome with endonucleases constitutes a powerful tool to study, in vivo, how DSBs are processed and how ncRNAs participate in this crucial event. In this review, we will summarize studies reporting the different roles of ncRNAs in DSB repair and discuss how genome editing approaches, especially CRISPR/Cas systems, can assist DNA repair studies. We will summarize knowledge concerning the functional significance of ncRNAs in DNA repair and their contribution to genome stability and integrity, with a focus on plants.

Genome-Wide Identification of Long Non-Coding RNAs and Their Regulatory Networks Involved in Apis mellifera ligustica Response to Nosema ceranae Infection

Long non-coding RNAs (lncRNAs) are a diverse class of transcripts that structurally resemble mRNAs but do not encode proteins, and lncRNAs have been proven to play pivotal roles in a wide range of biological processes in animals and plants. However, knowledge of expression patterns and potential roles of honeybee lncRNA response to Nosema ceranae infection is completely unknown. Here, we performed whole transcriptome strand-specific RNA sequencing of normal midguts of Apis mellifera ligustica workers (Am7CK, Am10CK) and N. ceranae-inoculated midguts (Am7T, Am10T), followed by comprehensive analyses using bioinformatic and molecular approaches. A total of 6353 A. m. ligustica lncRNAs were identified, including 4749 conserved lncRNAs and 1604 novel lncRNAs. These lncRNAs had minimal sequence similarities with other known lncRNAs in other species; however, their structural features were similar to counterparts in mammals and plants, including shorter exon and intron length, lower exon number, and lower expression level, compared with protein-coding transcripts. Further, 111 and 146 N. ceranae-responsive lncRNAs were identified from midguts at 7-days post-inoculation (dpi) and 10 dpi compared with control midguts. Twelve differentially expressed lncRNAs (DElncRNAs) were shared by Am7CK vs. Am7T and Am10CK vs. Am10T comparison groups, while the numbers of unique DElncRNAs were 99 and 134, respectively. Functional annotation and pathway analysis showed that the DElncRNAs may regulate the expression of neighboring genes by acting in cis and trans fashion. Moreover, we discovered 27 lncRNAs harboring eight known miRNA precursors and 513 lncRNAs harboring 2257 novel miRNA precursors. Additionally, hundreds of DElncRNAs and their target miRNAs were found to form complex competitive endogenous RNA (ceRNA) networks, suggesting that these DElncRNAs may act as miRNA sponges. Furthermore, DElncRNA-miRNA-mRNA networks were constructed and investigated, the results demonstrated that a portion of the DElncRNAs were likely to participate in regulating the host material and energy metabolism as well as cellular and humoral immune host responses to N. ceranae invasion. Our findings revealed here offer not only a rich genetic resource for further investigation of the functional roles of lncRNAs involved in the A. m. ligustica response to N. ceranae infection, but also a novel insight into understanding the host-pathogen interaction during honeybee microsporidiosis.

The Involvement of Long Noncoding RNAs in Response to Plant Stress

Plant growth and productivity are greatly impacted by environmental stresses. Therefore, plants have evolved mechanisms which allow them to adapt to abiotic stresses through alterations in gene expression and metabolism. In recent years, studies have investigated the role of long noncoding RNA (lncRNA) in regulating gene expression in plants and characterized their involvement in various biological functions through their regulation of DNA methylation, DNA structural modifications, histone modifications, and RNA-RNA interactions. Genome-wide transcriptome analyses have identified various types of noncoding RNAs (ncRNAs) that respond to abiotic stress. These ncRNAs are in addition to the well-known housekeeping ncRNAs, such as rRNAs, tRNAs, snoRNAs, and snRNAs. In this review, recent research pertaining to the role of lncRNAs in the response of plants to abiotic stress is summarized and discussed.

Pervasive, Genome-Wide Transcription in the Organelle Genomes of Diverse Plastid-Bearing Protists

Organelle genomes are among the most sequenced kinds of chromosome. This is largely because they are small and widely used in molecular studies, but also because next-generation sequencing technologies made sequencing easier, faster, and cheaper. However, studies of organelle RNA have not kept pace with those of DNA, despite huge amounts of freely available eukaryotic RNA-sequencing (RNA-seq) data. Little is known about organelle transcription in nonmodel species, and most of the available eukaryotic RNA-seq data have not been mined for organelle transcripts. Here, we use publicly available RNA-seq experiments to investigate organelle transcription in 30 diverse plastid-bearing protists with varying organelle genomic architectures. Mapping RNA-seq data to organelle genomes revealed pervasive, genome-wide transcription, regardless of the taxonomic grouping, gene organization, or noncoding content. For every species analyzed, transcripts covered ≥85% of the mitochondrial and/or plastid genomes (all of which were ≤105 kb), indicating that most of the organelle DNA—coding and noncoding—is transcriptionally active. These results follow earlier studies of model species showing that organellar transcription is coupled and ubiquitous across the genome, requiring significant downstream processing of polycistronic transcripts. Our findings suggest that noncoding organelle DNA can be transcriptionally active, raising questions about the underlying function of these transcripts and underscoring the utility of publicly available RNA-seq data for recovering complete genome sequences. If pervasive transcription is also found in bigger organelle genomes (>105 kb) and across a broader range of eukaryotes, this could indicate that noncoding organelle RNAs are regulating fundamental processes within eukaryotic cells.

Transport and homeostasis of potassium and phosphate: limiting factors for sustainable crop production

Potassium (K) and phosphate (Pi) are both macronutrients essential for plant growth and crop production, but the unrenewable resources of phosphorus rock and potash have become limiting factors for food security. One critical measure to help solve this problem is to improve nutrient use efficiency (NUE) in plants by understanding and engineering genetic networks for ion uptake, translocation, and storage. Plants have evolved multiple systems to adapt to various nutrient conditions for growth and production. Within the NUE networks, transport proteins and their regulators are the primary players for maintaining nutrient homeostasis and could be utilized to engineer high NUE traits in crop plants. A large number of publications have detailed K+ and Pi transport proteins in plants over the past three decades. Meanwhile, the discovery and validation of their regulatory mechanisms are fast-track topics for research. Here, we provide an overview of K+ and Pi transport proteins and their regulatory mechanisms, which participate in the uptake, translocation, storage, and recycling of these nutrients in plants.

Present Scenario of Long Non-Coding RNAs in Plants

Small non-coding RNAs have been extensively studied in plants over the last decade. In contrast, genome-wide identification of plant long non-coding RNAs (lncRNAs) has recently gained momentum. LncRNAs are now being recognized as important players in gene regulation, and their potent regulatory roles are being studied comprehensively in eukaryotes. LncRNAs were first reported in humans in 1992. Since then, research in animals, particularly in humans, has rapidly progressed, and a vast amount of data has been generated, collected, and organized using computational approaches. Additionally, numerous studies have been conducted to understand the roles of these long RNA species in several diseases. However, the status of lncRNA investigation in plants lags behind that in animals (especially humans). Efforts are being made in this direction using computational tools and high-throughput sequencing technologies, such as the lncRNA microarray technique, RNA-sequencing (RNA-seq), RNA capture sequencing, (RNA CaptureSeq), etc. Given the current scenario, significant amounts of data have been produced regarding plant lncRNAs, and this amount is likely to increase in the subsequent years. In this review we have documented brief information about lncRNAs and their status of research in plants, along with the plant-specific resources/databases for information retrieval on lncRNAs.

Coordinating Replication with Transcription

DNA topological transitions occur when replication forks encounter other DNA transactions such as transcription. Failure in resolving such conflicts leads to generation of aberrant replication and transcription intermediates that might have adverse effects on genome stability. Cells have evolved numerous surveillance mechanisms to avoid, tolerate, and resolve such replication-transcription conflicts. Defects or non-coordination in such cellular mechanisms might have catastrophic effect on cell viability. In this chapter, we review consequences of replication encounters with transcription and its associated events, topological challenges, and how these inevitable conflicts alter the genome structure and functions.

New technologies accelerate the exploration of non-coding RNAs in horticultural plants

Non-coding RNAs (ncRNAs), that is, RNAs not translated into proteins, are crucial regulators of a variety of biological processes in plants. While protein-encoding genes have been relatively well-annotated in sequenced genomes, accounting for a small portion of the genome space in plants, the universe of plant ncRNAs is rapidly expanding. Recent advances in experimental and computational technologies have generated a great momentum for discovery and functional characterization of ncRNAs. Here we summarize the classification and known biological functions of plant ncRNAs, review the application of next-generation sequencing (NGS) technology and ribosome profiling technology to ncRNA discovery in horticultural plants and discuss the application of new technologies, especially the new genome-editing tool clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated protein 9 (Cas9) systems, to functional characterization of plant ncRNAs.

Complex Regulation of Plant Phosphate Transporters and the Gap between Molecular Mechanisms and Practical Application: What Is Missing?

It has been almost 25 years since the first report of the gene encoding a high-affinity phosphate transporter (PT), PHO84, in yeast. Since then, an increasing number of yeast PHO84 homologs as well as other genes encoding proteins with phosphate (Pi) transport activities have been identified and functionally characterized in diverse plant species. Great progress has been made also in deciphering the molecular mechanism underlying the regulation of the abundance and/or activity of these genes and their products. The regulatory genes affect plant Pi homeostasis commonly through direct or indirect regulation of the abundance of PTs at different levels. However, little has been achieved in the use of PTs for developing genetically modified crops with high phosphorus use efficiency (PUE). This might be a consequence of overemphasizing Pi uptake from the rhizosphere and lack of knowledge about the roles of PTs in Pi transport and recycling within the plant that are required to optimize PUE. Here, we mainly focused on the genes encoding proteins with Pi transport activities and the emerging understanding of their regulation at the transcriptional, post-transcriptional, translational, and post-translational levels. In addition, we propose potential strategies for effective use of PTs in improving plant growth and development.

Functions of OsDof25 in regulation of OsC4PPDK

Relative little is known about the functions of the so-called Dof zinc factors in plants. Here we report on the analysis of OsDof25 and show a function in regulation of the important C4 photosynthesis gene, OsC4PPDK in rice. Over-expression of OsDof25 enhanced the expression of OsC4PPDK in transient expression experiments by binding in a specific way to a conserved Dof binding site which was confirmed by yeast and in vitro binding studies. Expression studies using promoter GUS plants as well as qPCR experiments showed that OsDof25 expressed in different tissues including both photosynthetic and non-photosynthetic organs and that expression of OsDof25 was partially overlapping with the OsC4PPDK gene. Conclusive evidence for a role of OsDof25 in regulation of C4PPDK came from loss-of-function and gain-of-function experiments with transgenic rice, which showed that down-regulation or over-expression of OsDof25 correlated with OsC4PPDK expression and that OsDof25 has functions as transcriptional activator.

Lncing epigenetic control of transcription to cardiovascular development and disease

Transcriptional and epigenetic regulation is critical for proper heart development, cardiac homeostasis, and pathogenesis. Long noncoding RNAs have emerged as key components of the transcriptional regulatory pathways that govern cardiac development as well as stress response, signaling, and remodeling in cardiac pathologies. Within the past few years, studies have identified many long noncoding RNAs in the context of cardiovascular biology and have begun to reveal the key functions of these transcripts. In this review, we discuss the growing roles of long noncoding RNAs in different aspects of cardiovascular development as well as pathological responses during injury or disease. In addition, we discuss diverse mechanisms by which long noncoding RNAs orchestrate cardiac transcriptional programs. Finally, we explore the exciting potential of this novel class of transcripts as biomarkers and novel therapeutic targets for cardiovascular diseases.

Genome-wide characterization of long intergenic non-coding RNAs (lincRNAs) provides new insight into viral diseases in honey bees Apis cerana and Apis mellifera

BACKGROUND

Long non-coding RNAs (lncRNAs) are a class of RNAs that do not encode proteins. Recently, lncRNAs have gained special attention for their roles in various biological process and diseases.

RESULTS

In an attempt to identify long intergenic non-coding RNAs (lincRNAs) and their possible involvement in honey bee development and diseases, we analyzed RNA-seq datasets generated from Asian honey bee (Apis cerana) and western honey bee (Apis mellifera). We identified 2470 lincRNAs with an average length of 1011 bp from A. cerana and 1514 lincRNAs with an average length of 790 bp in A. mellifera. Comparative analysis revealed that 5 % of the total lincRNAs derived from both species are unique in each species. Our comparative digital gene expression analysis revealed a high degree of tissue-specific expression among the seven major tissues of honey bee, different from mRNA expression patterns. A total of 863 (57 %) and 464 (18 %) lincRNAs showed tissue-dependent expression in A. mellifera and A. cerana, respectively, most preferentially in ovary and fat body tissues. Importantly, we identified 11 lincRNAs that are specifically regulated upon viral infection in honey bees, and 10 of them appear to play roles during infection with various viruses.

CONCLUSIONS

This study provides the first comprehensive set of lincRNAs for honey bees and opens the door to discover lincRNAs associated with biological and hormone signaling pathways as well as various diseases of honey bee.

Exploring the secrets of long noncoding RNAs

High-throughput sequencing has revealed that the majority of RNAs have no capacity to encode protein. Among these non-coding transcripts, recent work has focused on the roles of long noncoding RNAs (lncRNAs) of >200 nucleotides. Although many of their attributes, such as patterns of expression, remain largely unknown, lncRNAs have key functions in transcriptional, post-transcriptional, and epigenetic gene regulation; Also, new work indicates their functions in scaffolding ribonuclear protein complexes. In plants, genome-wide identification of lncRNAs has been conducted in several species, including Zea mays, and recent research showed that lncRNAs regulate flowering time in the photoperiod pathway, and function in nodulation. In this review, we discuss the basic mechanisms by which lncRNAs regulate key cellular processes, using the large body of knowledge on animal and yeast lncRNAs to illustrate the significance of emerging work on lncRNAs in plants.

Epigenetic regulation by long noncoding RNAs in plants

Many eukaryotes, including plants, produce a large number of long noncoding RNAs (lncRNAs).Growing number of lncRNAs are being reported to have regulatory roles in various developmental processes.Emerging mechanisms underlying the function of lncRNAs indicate that lncRNAs are versatile regulatory molecules. They function as potent cis- and trans-regulators of gene expression, including the formation of modular scaffolds that recruit chromatin-modifying complexes to target chromatin. LncRNAs have also been reported in plants. Here, we describe our current understanding on potential roles of lncRNA in plants.

The impact of widespread regulatory neofunctionalization on homeolog gene evolution following whole-genome duplication in maize

Whole-genome duplications are a widespread feature of plant genome evolution, having been detected in all flowering plant lineages. Despite the prevalence of these events, the extent to which duplicated genes (homeolog gene pairs) functionally diverge (neofunctionalization) is unclear. We present a genome-wide analysis of molecular evolution and regulatory neofunctionalization in maize (Zea mays L.). We demonstrate that 13% of all homeolog gene pairs in maize are regulatory neofunctionalized in leaves, and that regulatory neofunctionalized genes experience enhanced purifying selection. We show that significantly more genes have been regulatory neofunctionalized in foliar leaves than in husk leaves and that both leaf types have experienced selection for distinct functional roles. Furthermore, we demonstrate that biased subgenome expression dominance occurs only in the presence of regulatory neofunctionalization and that in nonregulatory neofunctionalized genes subgenome dominance is progressively acquired during development. Taken together, our study reveals several novel insights into the evolution of maize, genes, and gene expression, and provides a general model for gene evolution following whole-genome duplication in plants.

Patterns of homoeologous gene expression shown by RNA sequencing in hexaploid bread wheat

BACKGROUND

Bread wheat (Triticum aestivum) has a large, complex and hexaploid genome consisting of A, B and D homoeologous chromosome sets. Therefore each wheat gene potentially exists as a trio of A, B and D homoeoloci, each of which may contribute differentially to wheat phenotypes. We describe a novel approach combining wheat cytogenetic resources (chromosome substitution 'nullisomic-tetrasomic' lines) with next generation deep sequencing of gene transcripts (RNA-Seq), to directly and accurately identify homoeologue-specific single nucleotide variants and quantify the relative contribution of individual homoeoloci to gene expression.

RESULTS

We discover, based on a sample comprising ~5-10% of the total wheat gene content, that at least 45% of wheat genes are expressed from all three distinct homoeoloci. Most of these genes show strikingly biased expression patterns in which expression is dominated by a single homoeolocus. The remaining ~55% of wheat genes are expressed from either one or two homoeoloci only, through a combination of extensive transcriptional silencing and homoeolocus loss.

CONCLUSIONS

We conclude that wheat is tending towards functional diploidy, through a variety of mechanisms causing single homoeoloci to become the predominant source of gene transcripts. This discovery has profound consequences for wheat breeding and our understanding of wheat evolution.

Regulation of metabolism by long, non-coding RNAs

Our understanding of genomic regulation was revolutionized by the discovery that the genome is pervasively transcribed, giving rise to thousands of mostly uncharacterized non-coding ribonucleic acids (ncRNAs). Long, ncRNAs (lncRNAs) have thus emerged as a novel class of functional RNAs that impinge on gene regulation by a broad spectrum of mechanisms such as the recruitment of epigenetic modifier proteins, control of mRNA decay and DNA sequestration of transcription factors. We review those lncRNAs that are implicated in differentiation and homeostasis of metabolic tissues and present novel concepts on how lncRNAs might act on energy and glucose homeostasis. Finally, the control of circadian rhythm by lncRNAs is an emerging principles of lncRNA-mediated gene regulation.

Genome-wide discovery and characterization of maize long non-coding RNAs

BACKGROUND

Long non-coding RNAs (lncRNAs) are transcripts that are 200 bp or longer, do not encode proteins, and potentially play important roles in eukaryotic gene regulation. However, the number, characteristics and expression inheritance pattern of lncRNAs in maize are still largely unknown.

RESULTS

By exploiting available public EST databases, maize whole genome sequence annotation and RNA-seq datasets from 30 different experiments, we identified 20,163 putative lncRNAs. Of these lncRNAs, more than 90% are predicted to be the precursors of small RNAs, while 1,704 are considered to be high-confidence lncRNAs. High confidence lncRNAs have an average transcript length of 463 bp and genes encoding them contain fewer exons than annotated genes. By analyzing the expression pattern of these lncRNAs in 13 distinct tissues and 105 maize recombinant inbred lines, we show that more than 50% of the high confidence lncRNAs are expressed in a tissue-specific manner, a result that is supported by epigenetic marks. Intriguingly, the inheritance of lncRNA expression patterns in 105 recombinant inbred lines reveals apparent transgressive segregation, and maize lncRNAs are less affected by cis- than by trans-genetic factors.

CONCLUSIONS

We integrate all available transcriptomic datasets to identify a comprehensive set of maize lncRNAs, provide a unique annotation resource of the maize genome and a genome-wide characterization of maize lncRNAs, and explore the genetic control of their expression using expression quantitative trait locus mapping.

Application of resequencing to rice genomics, functional genomics and evolutionary analysis

Rice is a model system used for crop genomics studies. The completion of the rice genome draft sequences in 2002 not only accelerated functional genome studies, but also initiated a new era of resequencing rice genomes. Based on the reference genome in rice, next-generation sequencing (NGS) using the high-throughput sequencing system can efficiently accomplish whole genome resequencing of various genetic populations and diverse germplasm resources. Resequencing technology has been effectively utilized in evolutionary analysis, rice genomics and functional genomics studies. This technique is beneficial for both bridging the knowledge gap between genotype and phenotype and facilitating molecular breeding via gene design in rice. Here, we also discuss the limitation, application and future prospects of rice resequencing.

Arabidopsis non-coding RNA regulation in abiotic stress responses

Plant growth and productivity are largely affected by environmental stresses. Therefore, plants have evolved unique adaptation mechanisms to abiotic stresses through fine-tuned adjustment of gene expression and metabolism. Recent advanced technologies, such as genome-wide transcriptome analysis, have revealed that a vast amount of non-coding RNAs (ncRNAs) apart from the well-known housekeeping ncRNAs such as rRNAs, tRNAs, small nuclear RNAs (snRNAs) and small nucleolar RNAs (snoRNAs) are expressed under abiotic stress conditions. These various types of ncRNAs are involved in chromatin regulation, modulation of RNA stability and translational repression during abiotic stress response. In this review, we summarize recent progress that has been made on ncRNA research in plant abiotic stress response.

Genome-wide survey of cold stress regulated alternative splicing in Arabidopsis thaliana with tiling microarray

Alternative splicing plays a major role in expanding the potential informational content of eukaryotic genomes. It is an important post-transcriptional regulatory mechanism that can increase protein diversity and affect mRNA stability. Alternative splicing is often regulated in a tissue-specific and stress-responsive manner. Cold stress, which adversely affects plant growth and development, regulates the transcription and splicing of plant splicing factors. This can affect the pre-mRNA processing of many genes. To identify cold regulated alternative splicing we applied Affymetrix Arabidopsis tiling arrays to survey the transcriptome under cold treatment conditions. A novel algorithm was used for detection of statistically relevant changes in intron expression within a transcript between control and cold growth conditions. A reverse transcription polymerase chain reaction (RT-PCR) analysis of a number of randomly selected genes confirmed the changes in splicing patterns under cold stress predicted by tiling array. Our analysis revealed new types of cold responsive genes. While their expression level remains relatively unchanged under cold stress their splicing pattern shows detectable changes in the relative abundance of isoforms. The majority of cold regulated alternative splicing introduced a premature termination codon (PTC) into the transcripts creating potential targets for degradation by the nonsense mediated mRNA decay (NMD) process. A number of these genes were analyzed in NMD-defective mutants by RT-PCR and shown to evade NMD. This may result in new and truncated proteins with altered functions or dominant negative effects. The results indicate that cold affects both quantitative and qualitative aspects of gene expression.

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.

Analysis of allele-specific gene expression using a target-oriented tiling microarray assay

Advances in molecular and computational biology in recent years have led to the development or the improvement of methods for analyzing global gene expression. In most of these efforts, it is assumed that alleles of different origins contribute equally to the transcript pool. However, accumulating evidence suggests that many genes are not equally expressed from the paternal and maternal chromosomes. In addition to imprinting, the phenomenon of imbalanced allelic expression is widespread in heterozygous individuals. To distinguish transcript pools derived from different alleles present in the same organism, a number of methods have been developed. Here, we describe an oligonucleotide tiling microarray-based assay for analyzing allele-specific gene expression. Specifically targeting single-nucleotide polymorphisms, this two-step assay offers a high-throughput and multiplex method for detecting and quantifying unequal allelic expression that is readily applicable to many experimental systems.

Integrative analysis of chromatin states in Arabidopsis identified potential regulatory mechanisms for natural antisense transcript production

Genome-wide analyses of epigenomic and transcriptomic profiles provide extensive resources for discovering epigenetic regulatory mechanisms. However, the construction of functionally relevant hypotheses from correlative patterns and the rigorous testing of these hypotheses may be challenging. We combined bioinformatics-driven hypothesis building with mutant analyses to identify potential epigenetic mechanisms using the model plant Arabidopsis thaliana. Genome-wide maps of nine histone modifications produced by ChIP-seq were used together with a strand-specific RNA-seq dataset to profile the epigenome and transcriptome of Arabidopsis. Combinatorial chromatin patterns were described by 42 major chromatin states with selected states validated using the re-ChIP assay. The functional relevance of chromatin modifications was analyzed using the ANchored CORrelative Pattern (ANCORP) method and a newly developed state-specific effects analysis (SSEA) method, which interrogates individual chromatin marks in the context of combinatorial chromatin states. Based on results from these approaches, we propose the hypothesis that cytosine methylation (5mC) and histone methylation H3K36me may synergistically repress production of natural antisense transcripts (NATs) in the context of actively expressed genes. Mutant analyses supported this proposed model at a significant proportion of the tested loci. We further identified polymerase-associated factor as a potential repressor for NAT abundance. Although the majority of tested NATs were found to localize to the nucleus, we also found evidence for cytoplasmically partitioned NATs. The significance of the subcellular localization of NATs and their biological functions remain to be defined.

Visual analysis of transcriptome data in the context of anatomical structures and biological networks

The complexity and temporal as well as spatial resolution of transcriptome datasets is constantly increasing due to extensive technological developments. Here we present methods for advanced visualization and intuitive exploration of transcriptomics data as necessary prerequisites in order to facilitate the gain of biological knowledge. Color-coding of structural images based on the expression level enables a fast visual data analysis in the background of the examined biological system. The network-based exploration of these visualizations allows for comparative analysis of genes with specific transcript patterns and supports the extraction of functional relationships even from large datasets. In order to illustrate the presented methods, the tool HIVE was applied for visualization and exploration of database-retrieved expression data for master regulators of Arabidopsis thaliana flower and seed development in the context of corresponding tissue-specific regulatory networks.

The grapevine expression atlas reveals a deep transcriptome shift driving the entire plant into a maturation program

We developed a genome-wide transcriptomic atlas of grapevine (Vitis vinifera) based on 54 samples representing green and woody tissues and organs at different developmental stages as well as specialized tissues such as pollen and senescent leaves. Together, these samples expressed ∼91% of the predicted grapevine genes. Pollen and senescent leaves had unique transcriptomes reflecting their specialized functions and physiological status. However, microarray and RNA-seq analysis grouped all the other samples into two major classes based on maturity rather than organ identity, namely, the vegetative/green and mature/woody categories. This division represents a fundamental transcriptomic reprogramming during the maturation process and was highlighted by three statistical approaches identifying the transcriptional relationships among samples (correlation analysis), putative biomarkers (O2PLS-DA approach), and sets of strongly and consistently expressed genes that define groups (topics) of similar samples (biclustering analysis). Gene coexpression analysis indicated that the mature/woody developmental program results from the reiterative coactivation of pathways that are largely inactive in vegetative/green tissues, often involving the coregulation of clusters of neighboring genes and global regulation based on codon preference. This global transcriptomic reprogramming during maturation has not been observed in herbaceous annual species and may be a defining characteristic of perennial woody plants.

Molecular Functions of Long Non-Coding RNAs in Plants

The past decade has seen dramatic changes in our understanding of the scale and complexity of eukaryotic transcriptome owing to the discovery of diverse types of short and long non-protein-coding RNAs (ncRNAs). While short ncRNA-mediated gene regulation has been extensively studied and the mechanisms well understood, the function of long ncRNAs remains largely unexplored, especially in plants. Nevertheless, functional insights generated in recent studies with mammalian systems have indicated that long ncRNAs are key regulators of a variety of biological processes. They have been shown to act as transcriptional regulators and competing endogenous RNAs (ceRNAs), to serve as molecular cargos for protein re-localization and as modular scaffolds to recruit the assembly of multiple protein complexes for chromatin modifications. Some of these functions have been found to be conserved in plants. Here, we review our current understanding of long ncRNA functions in plants and discuss the challenges in functional characterization of plant long ncRNAs.

MicroRNA-135a contributes to the development of portal vein tumor thrombus by promoting metastasis in hepatocellular carcinoma

BACKGROUND & AIMS

Portal vein tumor thrombus (PVTT) has previously been demonstrated to correlate with poor prognosis of hepatocellular carcinoma. Approximately 50-80% of HCC is accompanied by portal or hepatic vein invasion. The underlying mechanisms of PVTT development remain unclear. This study aimed to elucidate the role of miR-135a in PVTT tumorigenesis.

METHODS

In the present study, we investigated the expression of microRNAs and mRNAs in PVTT tissues using advanced microRNA and cDNA microarray techniques. MicroRNA (miR)-135a was noted to be highly over-expressed in PVTT and the cell line CSQT-2 and was selected for further study. We characterized the function of miR-135a in vitro and in vivo. We also analyzed the clinical relevance of miR-135a in relation to the prognosis and survival of HCC patients with PVTT.

RESULTS

Our analyses found that the miRNA and mRNA expression profiles of PVTT were distinct from the parenchyma tumor. Overexpression of miR-135a favors invasive and metastatic behavior in vitro. Furthermore, in a CSQT-2 orthotopic transplantation nude mouse model, blockade of miR-135a significantly reduced PVTT incidence. We also found that miR-135a was transcribed by forkhead box M1 (FOXM1), and metastasis suppressor 1 (MTSS1) was identified as the direct and functional target of miR-135a. Additionally, the cohort analysis revealed the relevance of miR-135a with respect to the prognosis and survival of HCC patients with PVTT.

CONCLUSIONS

Our data suggest an important role for miR-135a in promoting PVTT tumorigenesis and indicate the potential application of miR-135a in PVTT therapy.

"Hypothesis for the modern RNA world": a pervasive non-coding RNA-based genetic regulation is a prerequisite for the emergence of multicellular complexity

The transitions to multicellularity mark the most pivotal and distinctive events in life's history on Earth. Although several transitions to "simple" multicellularity (SM) have been recorded in both bacterial and eukaryotic clades, transitions to complex multicellularity (CM) have only happened a few times in eukaryotes. A large number of cell types (associated with large body size), increased energy consumption per gene expressed, and an increment of non-protein-coding DNA positively correlate with CM. These three factors can indeed be understood as the causes and consequences of the regulation of gene expression. Here, we discuss how a vast expansion of non-protein-coding RNA (ncRNAs) regulators rather than large numbers of novel protein regulators can easily contribute to the emergence of CM. We also propose that the evolutionary advantage of RNA-based gene regulation derives from the robustness of the RNA structure that makes it easy to combine genetic drift with functional exploration. We describe a model which aims to explain how the evolutionary dynamic of ncRNAs becomes dominated by the accessibility of advantageous mutations to innovate regulation in complex multicellular organisms. The information and models discussed here outline the hypothesis that pervasive ncRNA-based regulatory systems, only capable of being expanded and explored in higher eukaryotes, are prerequisite to complex multicellularity. Thereby, regulatory RNA molecules in Eukarya have allowed intensification of morphological complexity by stabilizing critical phenotypes and controlling developmental precision. Although the origin of RNA on early Earth is still controversial, it is becoming clear that once RNA emerged into a protocellular system, its relevance within the evolution of biological systems has been greater than we previously thought.

On the opportunistic nature of transcription and replication initiation in the metazoan genome

Cellular identity and its response to external or internal signalling variations are encoded in a cell's genome as regulatory information. The genomic regions that specify this type of information are highly variable and degenerated in their sequence determinants, as it is becoming increasingly evident through the application of genome-scale methods to study gene expression. Here, we speculate that the same scenario applies to the regulatory regions controlling where DNA replication starts in the metazoan genome. We propose that replication origins cannot be defined as unique genomic features, but rather that DNA synthesis initiates opportunistically from accessible DNA sites, making cells highly robust and adaptable to environmental or developmental changes.

Contrasting patterns of evolution following whole genome versus tandem duplication events in Populus

Comparative analysis of multiple angiosperm genomes has implicated gene duplication in the expansion and diversification of many gene families. However, empirical data and theory suggest that whole-genome and small-scale duplication events differ with respect to the types of genes preserved as duplicate pairs. We compared gene duplicates resulting from a recent whole genome duplication to a set of tandemly duplicated genes in the model forest tree Populus trichocarpa. We used a combination of microarray expression analyses of a diverse set of tissues and functional annotation to assess factors related to the preservation of duplicate genes of both types. Whole genome duplicates are 700 bp longer and are expressed in 20% more tissues than tandem duplicates. Furthermore, certain functional categories are over-represented in each class of duplicates. In particular, disease resistance genes and receptor-like kinases commonly occur in tandem but are significantly under-retained following whole genome duplication, while whole genome duplicate pairs are enriched for members of signal transduction cascades and transcription factors. The shape of the distribution of expression divergence for duplicated pairs suggests that nearly half of the whole genome duplicates have diverged in expression by a random degeneration process. The remaining pairs have more conserved gene expression than expected by chance, consistent with a role for selection under the constraints of gene balance. We hypothesize that duplicate gene preservation in Populus is driven by a combination of subfunctionalization of duplicate pairs and purifying selection favoring retention of genes encoding proteins with large numbers of interactions.

Experimental annotation of the human pathogen Histoplasma capsulatum transcribed regions using high-resolution tiling arrays

BACKGROUND

The fungal pathogen Histoplasma capsulatum is thought to be the most common cause of fungal respiratory infections in immunocompetent humans, yet little is known about its biology. Here we provide the first genome-wide studies to experimentally validate its genome annotation. A functional interrogation of the Histoplasma genome provides critical support for continued investigation into the biology and pathogenesis of H. capsulatum and related fungi.

RESULTS

We employed a three-pronged approach to provide a functional annotation for the H. capsulatum G217B strain. First, we probed high-density tiling arrays with labeled cDNAs from cells grown under diverse conditions. These data defined 6,172 transcriptionally active regions (TARs), providing validation of 6,008 gene predictions. Interestingly, 22% of these predictions showed evidence of anti-sense transcription. Additionally, we detected transcription of 264 novel genes not present in the original gene predictions. To further enrich our analysis, we incorporated expression data from whole-genome oligonucleotide microarrays. These expression data included profiling under growth conditions that were not represented in the tiling experiment, and validated an additional 2,249 gene predictions. Finally, we compared the G217B gene predictions to other available fungal genomes, and observed that an additional 254 gene predictions had an ortholog in a different fungal species, suggesting that they represent genuine coding sequences.

CONCLUSIONS

These analyses yielded a high confidence set of validated gene predictions for H. capsulatum. The transcript sets resulting from this study are a valuable resource for further experimental characterization of this ubiquitous fungal pathogen. The data is available for interactive exploration at http://histo.ucsf.edu.

Functional consequences of bidirectional promoters

Several studies have shown that promoters of protein-coding genes are origins of pervasive non-coding RNA transcription and can initiate transcription in both directions. However, only recently have researchers begun to elucidate the functional implications of this bidirectionality and non-coding RNA production. Increasing evidence indicates that non-coding transcription at promoters influences the expression of protein-coding genes, revealing a new layer of transcriptional regulation. This regulation acts at multiple levels, from modifying local chromatin to enabling regional signal spreading and more distal regulation. Moreover, the bidirectional activity of a promoter is regulated at multiple points during transcription, giving rise to diverse types of transcripts.

Evolutionary transients in the rice transcriptome

In the canonical version of evolution by gene duplication, one copy is kept unaltered while the other is free to evolve. This process of evolutionary experimentation can persist for millions of years. Since it is so short lived in comparison to the lifetime of the core genes that make up the majority of most genomes, a substantial fraction of the genome and the transcriptome may—in principle—be attributable to what we will refer to as “evolutionary transients”, referring here to both the process and the genes that have gone or are undergoing this process. Using the rice gene set as a test case, we argue that this phenomenon goes a long way towards explaining why there are so many more rice genes than Arabidopsis genes, and why most excess rice genes show low similarity to eudicots.

Plant noncoding RNA gene discovery by "single-genome comparative genomics"

Plant genomes have undergone multiple rounds of duplications that contributed massively to the growth of gene families. The structure of resulting families has been studied in depth for protein-coding genes. However, little is known about the impact of duplications on noncoding RNA (ncRNA) genes. Here we perform a systematic analysis of duplicated regions in the rice genome in search of such ncRNA repeats. We observe that, just like their protein counterparts, most ncRNA genes have undergone multiple duplications that left visible sequence conservation footprints. The extent of ncRNA gene duplication in plants is such that these sequence footprints can be exploited for the discovery of novel ncRNA gene families on a large scale. We developed an SVM model that is able to retrieve likely ncRNA candidates among the 100,000+ repeat families in the rice genome, with a reasonably low false-positive discovery rate. Among the nearly 4000 ncRNA families predicted by this means, only 90 correspond to putative snoRNA or miRNA families. About half of the remaining families are classified as structured RNAs. New candidate ncRNAs are particularly enriched in UTR and intronic regions. Interestingly, 89% of the putative ncRNA families do not produce a detectable signal when their sequences are compared to another grass genome such as maize. Our results show that a large fraction of rice ncRNA genes are present in multiple copies and are species-specific or of recent origin. Intragenome comparison is a unique and potent source for the computational annotation of this major class of ncRNA.

Ratio-based analysis of differential mRNA processing and expression of a polyadenylation factor mutant pcfs4 using arabidopsis tiling microarray

Background

Alternative polyadenylation as a mechanism in gene expression regulation has been widely recognized in recent years. Arabidopsis polyadenylation factor PCFS4 was shown to function in leaf development and in flowering time control. The function of PCFS4 in controlling flowering time was correlated with the alternative polyadenylation of FCA, a flowering time regulator. However, genetic evidence suggested additional targets of PCFS4 that may mediate its function in both flowering time and leaf development.

Methodology/Principal Findings

To identify further targets, we investigated the whole transcriptome of a PCFS4 mutant using Affymetrix Arabidopsis genomic tiling 1.0R array and developed a data analysis pipeline, termed RADPRE (Ratio-based Analysis of Differential mRNA Processing and Expression). In RADPRE, ratios of normalized probe intensities between wild type Columbia and a pcfs4 mutant were first generated. By doing so, one of the major problems of tiling array data—variations caused by differential probe affinity—was significantly alleviated. With the probe ratios as inputs, a hierarchy of statistical tests was carried out to identify differentially processed genes (DPG) and differentially expressed genes (DEG). The false discovery rate (FDR) of this analysis was estimated by using the balanced random combinations of Col/pcfs4 and pcfs4/Col ratios as inputs. Gene Ontology (GO) analysis of the DPGs and DEGs revealed potential new roles of PCFS4 in stress responses besides flowering time regulation.

Conclusion/Significance

We identified 68 DPGs and 114 DEGs with FDR at 1% and 2%, respectively. Most of the 68 DPGs were subjected to alternative polyadenylation, splicing or transcription initiation. Quantitative PCR analysis of a set of DPGs confirmed that most of these genes were truly differentially processed in pcfs4 mutant plants. The enriched GO term “regulation of flower development” among PCFS4 targets further indicated the efficacy of the RADPRE pipeline. This simple but effective program is available upon request.

Differentiation of the maize subgenomes by genome dominance and both ancient and ongoing gene loss

Ancient tetraploidies are found throughout the eukaryotes. After duplication, one copy of each duplicate gene pair tends to be lost (fractionate). For all studied tetraploidies, the loss of duplicated genes, known as homeologs, homoeologs, ohnologs, or syntenic paralogs, is uneven between duplicate regions. In maize, a species that experienced a tetraploidy 5-12 million years ago, we show that in addition to uneven ancient gene loss, the two complete genomes contained within maize are differentiated by ongoing fractionation among diverse inbreds as well as by a pattern of overexpression of genes from the genome that has experienced less gene loss. These expression differences are consistent over a range of experiments quantifying RNA abundance in different tissues. We propose that the universal bias in gene loss between the genomes of this ancient tetraploid, and perhaps all tetraploids, is the result of selection against loss of the gene responsible for the majority of total expression for a duplicate gene pair. Although the tetraploidy of maize is ancient, biased gene loss and expression continue today and explain, at least in part, the remarkable genetic diversity found among modern maize cultivars.

Development of a versatile, target-oriented tiling microarray assay for measuring allele-specific gene expression

In the study of gene expression, it is often desirable to distinguish transcript pools derived from different alleles present in the same organism. We report here an oligonucleotide tiling microarray designed to specifically target 518 single nucleotide polymorphisms (SNPs) between the two sequenced rice (Oryza sativa) subspecies indica and japonica. The tiling array included all 25-mer probes interrogating each SNP by placing the polymorphic site at all 25 possible positions within the probe. Through hybridization to a titration series in which the japonica- and indica-derived cDNA templates were mixed with altering proportions, a regression model was used to screen for diagnostic probe sets for each SNP. Our result indicates that 284 (55%) SNPs have at least one diagnostic probe pair suitable for distinguishing and quantifying the relative abundance of allele-specific transcripts. As a proof-of-concept, we analyzed allele-specific expression in reciprocal indica×japonica F(1) hybrids and detected imbalanced expression at approximately one third of the SNPs. These results were validated by RNA-sequencing and allele-specific real-time PCR experiments. Together, our work demonstrates the utility and advantages of the tiling array method in interrogating large numbers of SNPs for quantifying allele-specific gene expression.

Function annotation of the rice transcriptome at single-nucleotide resolution by RNA-seq

The functional complexity of the rice transcriptome is not yet fully elucidated, despite many studies having reported the use of DNA microarrays. Next-generation DNA sequencing technologies provide a powerful approach for mapping and quantifying the transcriptome, termed RNA sequencing (RNA-seq). In this study, we applied RNA-seq to globally sample transcripts of the cultivated rice Oryza sativa indica and japonica subspecies for resolving the whole-genome transcription profiles. We identified 15,708 novel transcriptional active regions (nTARs), of which 51.7% have no homolog to public protein data and >63% are putative single-exon transcripts, which are highly different from protein-coding genes (<20%). We found that approximately 48% of rice genes show alternative splicing patterns, a percentage considerably higher than previous estimations. On the basis of the available rice gene models, 83.1% (46,472 genes) of the current rice gene models were validated by RNA-seq, and 6228 genes were identified to be extended at the 5' and/or 3' ends by at least 50 bp. Comparative transcriptome analysis demonstrated that 3464 genes exhibited differential expression patterns. The ratio of SNPs with nonsynonymous/synonymous mutations was nearly 1:1.06. In total, we interrogated and compared transcriptomes of the two rice subspecies to reveal the overall transcriptional landscape at maximal resolution.

Experimental annotation of the human pathogen Candida albicans coding and noncoding transcribed regions using high-resolution tiling arrays

BACKGROUND

Compared to other model organisms and despite the clinical relevance of the pathogenic yeast Candida albicans, no comprehensive analysis has been done to provide experimental support of its in silico-based genome annotation.

RESULTS

We have undertaken a genome-wide experimental annotation to accurately uncover the transcriptional landscape of the pathogenic yeast C. albicans using strand-specific high-density tiling arrays. RNAs were purified from cells growing under conditions relevant to C. albicans pathogenicity, including biofilm, lab-grown yeast and serum-induced hyphae, as well as cells isolated from the mouse caecum. This work provides a genome-wide experimental validation for a large number of predicted ORFs for which transcription had not been detected by other approaches. Additionally, we identified more than 2,000 novel transcriptional segments, including new ORFs and exons, non-coding RNAs (ncRNAs) as well as convincing cases of antisense gene transcription. We also characterized the 5' and 3' UTRs of expressed ORFs, and established that genes with long 5' UTRs are significantly enriched in regulatory functions controlling filamentous growth. Furthermore, we found that genomic regions adjacent to telomeres harbor a cluster of expressed ncRNAs. To validate and confirm new ncRNA candidates, we adapted an iterative strategy combining both genome-wide occupancy of the different subunits of RNA polymerases I, II and III and expression data. This comprehensive approach allowed the identification of different families of ncRNAs.

CONCLUSIONS

In summary, we provide a comprehensive expression atlas that covers relevant C. albicans pathogenic developmental stages in addition to the discovery of new ORF and non-coding genetic elements.

Analysis of antisense expression by whole genome tiling microarrays and siRNAs suggests mis-annotation of Arabidopsis orphan protein-coding genes

Background

MicroRNAs (miRNAs) and trans-acting small-interfering RNAs (tasi-RNAs) are small (20–22 nt long) RNAs (smRNAs) generated from hairpin secondary structures or antisense transcripts, respectively, that regulate gene expression by Watson-Crick pairing to a target mRNA and altering expression by mechanisms related to RNA interference. The high sequence homology of plant miRNAs to their targets has been the mainstay of miRNA prediction algorithms, which are limited in their predictive power for other kingdoms because miRNA complementarity is less conserved yet transitive processes (production of antisense smRNAs) are active in eukaryotes. We hypothesize that antisense transcription and associated smRNAs are biomarkers which can be computationally modeled for gene discovery.

Principal Findings

We explored rice (Oryza sativa) sense and antisense gene expression in publicly available whole genome tiling array transcriptome data and sequenced smRNA libraries (as well as C. elegans) and found evidence of transitivity of MIRNA genes similar to that found in Arabidopsis. Statistical analysis of antisense transcript abundances, presence of antisense ESTs, and association with smRNAs suggests several hundred Arabidopsis ‘orphan’ hypothetical genes are non-coding RNAs. Consistent with this hypothesis, we found novel Arabidopsis homologues of some MIRNA genes on the antisense strand of previously annotated protein-coding genes. A Support Vector Machine (SVM) was applied using thermodynamic energy of binding plus novel expression features of sense/antisense transcription topology and siRNA abundances to build a prediction model of miRNA targets. The SVM when trained on targets could predict the “ancient” (deeply conserved) class of validated Arabidopsis MIRNA genes with an accuracy of 84%, and 76% for “new” rapidly-evolving MIRNA genes.

Conclusions

Antisense and smRNA expression features and computational methods may identify novel MIRNA genes and other non-coding RNAs in plants and potentially other kingdoms, which can provide insight into antisense transcription, miRNA evolution, and post-transcriptional gene regulation.