A tiling microarray for global analysis of chloroplast genome expression in cucumber and other plants

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.

Gene expression profiling of acute myeloid leukemia samples from adult patients with AML-M1 and -M2 through boutique microarrays, real-time PCR and droplet digital PCR

Acute myeloid leukemia (AML) is the most common and severe form of acute leukemia diagnosed in adults. Owing to its heterogeneity, AML is divided into classes associated with different treatment outcomes and specific gene expression profiles. Based on previous studies on AML, in this study, we designed and generated an AML-array containing 900 oligonucleotide probes complementary to human genes implicated in hematopoietic cell differentiation and maturation, proliferation, apoptosis and leukemic transformation. The AML-array was used to hybridize 118 samples from 33 patients with AML of the M1 and M2 subtypes of the French-American-British (FAB) classification and 15 healthy volunteers (HV). Rigorous analysis of the microarray data revealed that 83 genes were differentially expressed between the patients with AML and the HV, including genes not yet discussed in the context of AML pathogenesis. The most overexpressed genes in AML were STMN1, KITLG, CDK6, MCM5, KRAS, CEBPA, MYC, ANGPT1, SRGN, RPLP0, ENO1 and SET, whereas the most underexpressed genes were IFITM1, LTB, FCN1, BIRC3, LYZ, ADD3, S100A9, FCER1G, PTRPE, CD74 and TMSB4X. The overexpression of the CPA3 gene was specific for AML with mutated NPM1 and FLT3. Although the microarray-based method was insufficient to differentiate between any other AML subgroups, quantitative PCR approaches enabled us to identify 3 genes (ANXA3, S100A9 and WT1) whose expression can be used to discriminate between the 2 studied AML FAB subtypes. The expression levels of the ANXA3 and S100A9 genes were increased, whereas those of WT1 were decreased in the AML-M2 compared to the AML-M1 group. We also examined the association between the STMN1, CAT and ABL1 genes, and the FLT3 and NPM1 mutation status. FLT3⁺/NPM1⁻ AML was associated with the highest expression of STMN1, and ABL1 was upregulated in FLT3⁺ AML and CAT in FLT3⁻ AML, irrespectively of the NPM1 mutation status. Moreover, our results indicated that CAT and WT1 gene expression levels correlated with the response to therapy. CAT expression was highest in patients who remained longer under complete remission, whereas WT1 expression increased with treatment resistance. On the whole, this study demonstrates that the AML-array can potentially serve as a first-line screening tool, and may be helpful for the diagnosis of AML, whereas the differentiation between AML subgroups can be more successfully performed with PCR-based analysis of a few marker genes.

Cryptochrome 2 extensively regulates transcription of the chloroplast genome in tomato

Light plays a key role in the regulation of many physiological processes required for plant and chloroplast development. Plant cryptochromes (crys) play an important role in monitoring, capturing, and transmitting the light stimuli. In this study, we analyzed the effects of CRY2 overexpression on transcription of tomato chloroplast genome by a tiling array, containing about 90 000 overlapping probes (5‐nucleotide resolution). We profiled transcription in leaves of wild‐type and CRY2‐overexpressing plants grown in a diurnal cycle, to generate a comprehensive map of chloroplast transcription and to monitor potential specific modulations of the chloroplast transcriptome induced by the overexpression of CRY2. Our results demonstrate that CRY2 is a master gene of transcriptional regulation in the tomato chloroplast. In fact, it modulates the day/night mRNA abundance of about 58% of the 114 ORFs. The effect of CRY2 includes a differential extension of some transcripts at their 5′‐end, according to the period of the day. We observed that the influence of CRY2 on chloroplast transcription is not limited to coding RNA; a great number of putative noncoding micro RNA also showed differential accumulation pattern. To our knowledge, this is the first study that highlights how a photoreceptor affects the day/night transcription of the chloroplast genome.

ChloroSeq, an Optimized Chloroplast RNA-Seq Bioinformatic Pipeline, Reveals Remodeling of the Organellar Transcriptome Under Heat Stress

Although RNA-Seq has revolutionized transcript analysis, organellar transcriptomes are rarely assessed even when present in published datasets. Here, we describe the development and application of a rapid and convenient method, ChloroSeq, to delineate qualitative and quantitative features of chloroplast RNA metabolism from strand-specific RNA-Seq datasets, including processing, editing, splicing, and relative transcript abundance. The use of a single experiment to analyze systematically chloroplast transcript maturation and abundance is of particular interest due to frequent pleiotropic effects observed in mutants that affect chloroplast gene expression and/or photosynthesis. To illustrate its utility, ChloroSeq was applied to published RNA-Seq datasets derived from Arabidopsis thaliana grown under control and abiotic stress conditions, where the organellar transcriptome had not been examined. The most appreciable effects were found for heat stress, which induces a global reduction in splicing and editing efficiency, and leads to increased abundance of chloroplast transcripts, including genic, intergenic, and antisense transcripts. Moreover, by concomitantly analyzing nuclear transcripts that encode chloroplast gene expression regulators from the same libraries, we demonstrate the possibility of achieving a holistic understanding of the nucleus-organelle system. ChloroSeq thus represents a unique method for streamlining RNA-Seq data interpretation of the chloroplast transcriptome and its regulators.

A gene expression microarray for Nicotiana benthamiana based on de novo transcriptome sequence assembly

BACKGROUND

Nicotiana benthamiana has been widely used in laboratories around the world for studying plant-pathogen interactions and posttranscriptional gene expression silencing. Yet the exploration of its transcriptome has lagged behind due to the lack of both adequate sequence information and genome-wide analysis tools, such as DNA microarrays. Despite the increasing use of high-throughput sequencing technologies, the DNA microarrays still remain a popular gene expression tool, because they are cheaper and less demanding regarding bioinformatics skills and computational effort.

RESULTS

We designed a gene expression microarray with 103,747 60-mer probes, based on two recently published versions of N. benthamiana transcriptome (v.3 and v.5). Both versions were reconstructed from RNA-Seq data of non-strand-specific pooled-tissue libraries, so we defined the sense strand of the contigs prior to designing the probe. To accomplish this, we combined a homology search against Arabidopsis thaliana proteins and hybridization to a test 244k microarray containing pairs of probes, which represented individual contigs. We identified the sense strand in 106,684 transcriptome contigs and used this information to design an Nb-105k microarray on an Agilent eArray platform. Following hybridization of RNA samples from N. benthamiana roots and leaves we demonstrated that the new microarray had high specificity and sensitivity for detection of differentially expressed transcripts. We also showed that the data generated with the Nb-105k microarray may be used to identify incorrectly assembled contigs in the v.5 transcriptome, by detecting inconsistency in the gene expression profiles, which is indicated using multiple microarray probes that match the same v.5 primary transcripts.

CONCLUSIONS

We provided a complete design of an oligonucleotide microarray that may be applied to the research of N. benthamiana transcriptome. This, in turn, will allow the N. benthamiana research community to take full advantage of microarray capabilities for studying gene expression in this plant. Additionally, by defining the sense orientation of over 106,000 contigs, we substantially improved the functional information on the N. benthamiana transcriptome. The simple hybridization-based approach for detecting the sense orientation of computationally assembled sequences can be used for updating the transcriptomes of other non-model organisms, including cases where no significant homology to known proteins exists.

Next generation sequencing and omics in cucumber (Cucumis sativus L.) breeding directed research

In the post-genomic era the availability of genomic tools and resources is leading us to novel generation methods in plant breeding, as they facilitate the study of the genotype and its relationship with the phenotype, in particular for complex traits. In this study we have mainly concentrated on the Cucumis sativus and (but much less) Cucurbitaceae family several important vegetable crops. There are many reports on research conducted in Cucurbitaceae plant breeding programs on the ripening process, phloem transport, disease resistance, cold tolerance and fruit quality traits. This paper presents the role played by new omic technologies in the creation of knowledge on the mechanisms of the formation of the breeding features. The analysis of NGS (NGS-next generation sequencing) data allows the discovery of new genes and regulatory sequences, their positions, and makes available large collections of molecular markers. Genome-wide expression studies provide breeders with an understanding of the molecular basis of complex traits. Firstly a high density map should be created for the reference genome, then each re-sequencing data could be mapped and new markers brought out into breeding populations. The paper also presents methods that could be used in the future for the creation of variability and genomic modification of the species in question. It has been shown also the state and usefulness in breeding the chloroplastomic and mitochondriomic study.

Stress responsive miRNAs and isomiRs in cereals

Abiotic and biotic stress conditions are vital determinants in the production of cereals, the major caloric source in human nutrition. Small RNAs, miRNAs and isomiRs are central to post-transcriptional regulation of gene expression in a variety of cellular processes including development and stress responses. Several miRNAs have been identified using new technologies and have roles in stress responses in plants, including cereals. The overall knowledge about the cereal miRNA repertoire, as well as an understanding of complex miRNA mediated mechanisms of target regulation in response to stress conditions, is far from complete. Ongoing efforts that add to our understanding of complex miRNA machinery have implications in plant response to stress conditions. Additionally, sequence variants of miRNAs (isomiRNAs or isomiRs), regulation of their expression through dissection of upstream regulatory elements, the role of Processing-bodies (P-bodies) in miRNA exerted gene regulation and yet unveiled organellar plant miRNAs are newly emerging topics, which will contribute to the elucidation of the miRNA machinery and its role in cereal tolerance against abiotic and biotic stresses.

Arabidopsis chloroplast mini-ribonuclease III participates in rRNA maturation and intron recycling

RNase III proteins recognize double-stranded RNA structures and catalyze endoribonucleolytic cleavages that often regulate gene expression. Here, we characterize the functions of RNC3 and RNC4, two Arabidopsis thaliana chloroplast Mini-RNase III-like enzymes sharing 75% amino acid sequence identity. Whereas rnc3 and rnc4 null mutants have no visible phenotype, rnc3/rnc4 (rnc3/4) double mutants are slightly smaller and chlorotic compared with the wild type. In Bacillus subtilis, the RNase Mini-III is integral to 23S rRNA maturation. In Arabidopsis, we observed imprecise maturation of 23S rRNA in the rnc3/4 double mutant, suggesting that exoribonucleases generated staggered ends in the absence of specific Mini-III-catalyzed cleavages. A similar phenotype was found at the 3' end of the 16S rRNA, and the primary 4.5S rRNA transcript contained 3' extensions, suggesting that Mini-III catalyzes several processing events of the polycistronic rRNA precursor. The rnc3/4 mutant showed overaccumulation of a noncoding RNA complementary to the 4.5S-5S rRNA intergenic region, and its presence correlated with that of the extended 4.5S rRNA precursor. Finally, we found rnc3/4-specific intron degradation intermediates that are probable substrates for Mini-III and show that B. subtilis Mini-III is also involved in intron regulation. Overall, this study extends our knowledge of the key role of Mini-III in intron and noncoding RNA regulation and provides important insight into plastid rRNA maturation.

Analysis of boutique arrays: a universal method for the selection of the optimal data normalization procedure

DNA microarrays, which are among the most popular genomic tools, are widely applied in biology and medicine. Boutique arrays, which are small, spotted, dedicated microarrays, constitute an inexpensive alternative to whole-genome screening methods. The data extracted from each microarray-based experiment must be transformed and processed prior to further analysis to eliminate any technical bias. The normalization of the data is the most crucial step of microarray data pre-processing and this process must be carefully considered as it has a profound effect on the results of the analysis. Several normalization algorithms have been developed and implemented in data analysis software packages. However, most of these methods were designed for whole-genome analysis. In this study, we tested 13 normalization strategies (ten for double-channel data and three for single-channel data) available on R Bioconductor and compared their effectiveness in the normalization of four boutique array datasets. The results revealed that boutique arrays can be successfully normalized using standard methods, but not every method is suitable for each dataset. We also suggest a universal seven-step workflow that can be applied for the selection of the optimal normalization procedure for any boutique array dataset. The described workflow enables the evaluation of the investigated normalization methods based on the bias and variance values for the control probes, a differential expression analysis and a receiver operating characteristic curve analysis. The analysis of each component results in a separate ranking of the normalization methods. A combination of the ranks obtained from all the normalization procedures facilitates the selection of the most appropriate normalization method for the studied dataset and determines which methods can be used interchangeably.

Plastid non-coding RNAs: emerging candidates for gene regulation

Recent advances in transcriptomics and bioinformatics, specifically strand-specific RNA sequencing, have allowed high-throughput, comprehensive detection of low-abundance transcripts typical of the non-coding RNAs studied in bacteria and eukaryotes. Before this, few plastid non-coding RNAs (pncRNAs) had been identified, and even fewer had been investigated for any functional role in gene regulation. Relaxed plastid transcription initiation and termination result in full transcription of both chloroplast DNA strands. Following this, post-transcriptional processing produces a pool of metastable RNA species, including distinct pncRNAs. Here we review pncRNA biogenesis and possible functionality, and speculate that this RNA class may have an underappreciated role in plastid gene regulation.