The ancestral chromatin landscape of land plants.

Recent studies have shown that correlations between chromatin modifications and transcription vary among eukaryotes. This is the case for marked differences between the chromatin of the moss Physcomitrium patens and the liverwort Marchantia polymorpha. Mosses and liverworts diverged from hornworts, altogether forming the lineage of bryophytes that shared a common ancestor with land plants. We aimed to describe chromatin in hornworts to establish synapomorphies across bryophytes and approach a definition of the ancestral chromatin organization of land plants. We used genomic methods to define the 3D organization of chromatin and map the chromatin landscape of the model hornwort Anthoceros agrestis. We report that nearly half of the hornwort transposons were associated with facultative heterochromatin and euchromatin and formed the center of topologically associated domains delimited by protein coding genes. Transposons were scattered across autosomes, which contrasted with the dense compartments of constitutive heterochromatin surrounding the centromeres in flowering plants. Most of the features observed in hornworts are also present in liverworts or in mosses but are distinct from flowering plants. Hence, the ancestral genome of bryophytes was likely a patchwork of units of euchromatin interspersed within facultative and constitutive heterochromatin. We propose this genome organization was ancestral to land plants.

The Polycomb repressive complex 2 deposits H3K27me3 and represses transposable elements in a broad range of eukaryotes.

The mobility of transposable elements (TEs) contributes to evolution of genomes. Their uncontrolled activity causes genomic instability; therefore, expression of TEs is silenced by host genomes. TEs are marked with DNA and H3K9 methylation, which are associated with silencing in flowering plants, animals, and fungi. However, in distantly related groups of eukaryotes, TEs are marked by H3K27me3 deposited by the Polycomb repressive complex 2 (PRC2), an epigenetic mark associated with gene silencing in flowering plants and animals. The direct silencing of TEs by PRC2 has so far only been shown in one species of ciliates. To test if PRC2 silences TEs in a broader range of eukaryotes, we generated mutants with reduced PRC2 activity and analyzed the role of PRC2 in extant species along the lineage of Archaeplastida and in the diatom P. tricornutum. In this diatom and the red alga C. merolae, a greater proportion of TEs than genes were repressed by PRC2, whereas a greater proportion of genes than TEs were repressed by PRC2 in bryophytes. In flowering plants, TEs contained potential cis-elements recognized by transcription factors and associated with neighbor genes as transcriptional units repressed by PRC2. Thus, silencing of TEs by PRC2 is observed not only in Archaeplastida but also in diatoms and ciliates, suggesting that PRC2 deposited H3K27me3 to silence TEs in the last common ancestor of eukaryotes. We hypothesize that during the evolution of Archaeplastida, TE fragments marked with H3K27me3 were selected to shape transcriptional regulation, controlling networks of genes regulated by PRC2.

Comprehensive molecular characterization of complete mitogenome assemblies of 33 Eimeria isolates infecting domestic chickens.

BACKGROUND: Coccidiosis caused by Eimeria is one of the most severe chicken diseases and poses a great economic threat to the poultry industry. Understanding the evolutionary biology of chicken Eimeria parasites underpins development of new interactions toward the improved prevention and control of this poultry disease. METHODS: We presented an evolutionary blueprint of chicken coccidia by genetically characterizing complete mitogenome assemblies of 33 isolates representing all seven known Eimeria species infecting chickens in China. Further genome- and gene-level phylogenies were also achieved to better understand the evolutionary relationships of these chicken Eimeria at the species level. RESULTS: 33 mitogenomes of chicken eimerian parasites ranged from 6148 bp to 6480 bp in size and encoded typical mitochondrial compositions of apicomplexan parasites including three protein-coding genes (PCGs), seven fragmented small subunit (SSU) and 12/13 fragmented large subunit (LSU) rRNAs. Comparative genomics provided an evolutionary scenario for the genetic diversity of PCGs-cytochrome c oxidase subunits 1 and 3 (cox1 and cox3) and cytochrome b (cytb); all were under purifying selection with cox1 and cox3 being the lowest and highest evolutionary rates, respectively. Genome-wide phylogenies classified the 33 Eimeria isolates into seven subgroups, and furthermore Eimeria tenella and Eimeria necatrix were determined to be more closely related to each other than to the other eight congenic species. Single/concatenated mitochondrial protein gene-based phylogenies supported cox1 as the genetic marker for evolutionary and phylogenetic studies for avain coccidia. CONCLUSIONS: To our knowledge, these are the first population-level mitogenomic data on the genus Eimeria, and its comprehensive molecular characterization provides valuable resources for systematic, population genetic and evolutionary biological studies of apicomplexan parasites in poultry.

The Nodewalk assay to quantitate chromatin fiber interactomes in very small cell populations.

The chromosome conformation capture method and its derivatives, such as circularized chromosome conformation capture, carbon copy chromosome conformation capture, high-throughput chromosome conformation capture and capture high-throughput chromosome conformation capture, have pioneered our understanding of the principles of chromosome folding in the nucleus. These technical advances, however, cannot precisely quantitate interaction frequency in very small input samples. Here we describe a protocol for the Nodewalk assay, which is based on converting chromosome conformation capture DNA samples to RNA and subsequently to cDNA using strategically placed primers. This pipeline enables the quantitative analyses of chromatin fiber interactions without compromising its sensitivity down to <300 cells, making it suitable for MiSeq analyses of higher-order chromatin structures in biopsies, circulating tumor cells and transitional cell states, for example. Importantly, the quality of the Nodewalk sample can be assessed before sequencing to avoid unnecessary costs. Moreover, it enables analyses from hundreds of different restriction enzyme fragment viewpoints within the same initial small input sample to uncover complex, genome-wide networks. Following optimization of the different steps, the entire protocol can be completed within 2 weeks.

Comparative genomics of Sarcoptes scabiei provide new insights into adaptation to permanent parasitism and within-host species divergence.

Sarcoptes scabiei is the causative agent of a highly contagious skin disease in humans and more than 100 mammals. Here, we reported the first chromosome-level reference genome of S. scabiei isolated from rabbits, with a contig N50 size of 5.92 Mb, a total assembled length of 57.30 Mb, ∼12.65% repetitive sequences and 9333 predicted protein-coding genes. The phylogenetic tree based on 1338 shared high-confidence single-copy orthologous genes estimated that the mammalian ectoparasite S. scabiei and the plant-feeding mite Tetranychus urticae separated ∼340 million years ago. Both neighbour-joining tree and principal component analysis of 20 mite populations isolated from four hosts (humans, pigs, dogs and rabbits) distributed in three countries (China, Australia and the USA) consistently supported genetic subdivisions according to host species rather than geographical location. The demographic history of S. scabiei reconstructed by multiple sequentially Markovian coalescent analysis suggested that S. scabiei isolated from rabbits, humans, dogs and pigs diverged ∼5000 years ago. Investigation of the homeobox (Hox) genes revealed that S. scabiei contains 8 of 10 canonical Hox genes that are present in the arthropod ancestor, and the absence of the Abd-A gene may correlate with the long gap between their front and back legs. Comparative genomics demonstrated that genes specific to scabies mites were mainly enriched in nutrition digestive systems, whereas genes in the families that involved detoxification (cytochrome P450, carboxyl/cholinesterases and the ATP-binding cassette transporter C group) were extremely contracted compared with that of other mites analysed in this study. Selective sweep analysis of mite populations from various pairs of two out of the four host species revealed that the strongest selective sweep signals were mainly enriched in cysteine-type peptidase activity and apoptosis. The results provided clues for the mechanisms of S. scabiei adaptation to a permanent parasitic lifestyle and knowledge that would enable further control of this highly contagious skin disease.

Canonical WNT signaling-dependent gating of MYC requires a noncanonical CTCF function at a distal binding site.

Abnormal WNT signaling increases MYC expression in colon cancer cells in part via oncogenic super-enhancer-(OSE)-mediated gating of the active MYC to the nuclear pore in a poorly understood process. We show here that the principal tenet of the WNT-regulated MYC gating, facilitating nuclear export of the MYC mRNA, is regulated by a CTCF binding site (CTCFBS) within the OSE to confer growth advantage in HCT-116 cells. To achieve this, the CTCFBS directs the WNT-dependent trafficking of the OSE to the nuclear pore from intra-nucleoplasmic positions in a stepwise manner. Once the OSE reaches a peripheral position, which is triggered by a CTCFBS-mediated CCAT1 eRNA activation, its final stretch (≤0.7 µm) to the nuclear pore requires the recruitment of AHCTF1, a key nucleoporin, to the CTCFBS. Thus, a WNT/ß-catenin-AHCTF1-CTCF-eRNA circuit enables the OSE to promote pathological cell growth by coordinating the trafficking of the active MYC gene within the 3D nuclear architecture.

Genome of the Giant Panda Roundworm Illuminates Its Host Shift and Parasitic Adaptation.

Baylisascaris schroederi, a roundworm (ascaridoid) parasite specific to the bamboo-feeding giant panda (Ailuropoda melanoleuca), represents a leading cause of mortality in wild giant panda populations. Here, we present a 293-megabase chromosome-level genome assembly of B. schroederi to infer its biology, including host adaptations. Comparative genomics revealed an evolutionary trajectory accompanied by host-shift events in ascaridoid parasite lineages after host separations, suggesting their potential for transmission and rapid adaptation to new hosts. Genomic and anatomical lines of evidence, including expansion and positive selection of genes related to the cuticle and basal metabolisms, indicate that B. schroederi undergoes specific adaptations to survive in the sharp-edged bamboo-enriched gut of giant pandas by structurally increasing its cuticle thickness and efficiently utilizing host nutrients through gut parasitism. Additionally, we characterized the secretome of B. schroederi and predicted potential drug and vaccine targets for new control strategies. Overall, this genome resource provides new insights into the host adaptation of B. schroederi to the giant panda as well as the host-shift events in ascaridoid parasite lineages. Our findings on the unique biology of B. schroederi will also aid in the development of prevention and treatment measures to protect giant panda populations from roundworm parasitism.

Improved 93-11 Genome and Time-Course Transcriptome Expand Resources for Rice Genomics.

In 2002, the first crop genome was published using the rice cultivar 93-11, which is the progenitor of the first super-hybrid rice. The genome sequence has served as a reference genome for the indica cultivars, but the assembly has not been updated. In this study, we update the 93-11 genome assembly to a gap-less sequence using ultra-depth single molecule real-time (SMRT) reads, Hi-C sequencing, reference-guided, and gap-closing approach. The differences in the genome collinearity and gene content between the 93-11 and the Nipponbare reference genomes confirmed to map the indica cultivar sequencing data to the 93-11 genome, instead of the reference. Furthermore, time-course transcriptome data showed that the expression pattern was consistently correlated with the stages of seed development. Alternative splicing of starch synthesis-related genes and genomic variations of waxy make it a novel resource for targeted breeding. Collectively, the updated high quality 93-11 genome assembly can improve the understanding of the genome structures and functions of Oryza groups in molecular breeding programs.

IC4R-2.0: Rice Genome Reannotation Using Massive RNA-seq Data.

Genome reannotation aims for complete and accurate characterization of gene models and thus is of critical significance for in-depth exploration of gene function. Although the availability of massive RNA-seq data provides great opportunities for gene model refinement, few efforts have been made to adopt these precious data in rice genome reannotation. Here we reannotate the rice (Oryza sativa L. ssp. japonica) genome based on integration of large-scale RNA-seq data and release a new annotation system IC4R-2.0. In general, IC4R-2.0 significantly improves the completeness of gene structure, identifies a number of novel genes, and integrates a variety of functional annotations. Furthermore, long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) are systematically characterized in the rice genome. Performance evaluation shows that compared to previous annotation systems, IC4R-2.0 achieves higher integrity and quality, primarily attributable to massive RNA-seq data applied in genome annotation. Consequently, we incorporate the improved annotations into the Information Commons for Rice (IC4R), a database integrating multiple omics data of rice, and accordingly update IC4R by providing more user-friendly web interfaces and implementing a series of practical online tools. Together, the updated IC4R, which is equipped with the improved annotations, bears great promise for comparative and functional genomic studies in rice and other monocotyledonous species. The IC4R-2.0 annotation system and related resources are freely accessible at http://ic4r.org/.

Structural and Functional Annotation of Transposable Elements Revealed a Potential Regulation of Genes Involved in Rubber Biosynthesis by TE-Derived siRNA Interference in Hevea brasiliensis.

The natural rubber biosynthetic pathway is well described in Hevea, although the final stages of rubber elongation are still poorly understood. Small Rubber Particle Proteins and Rubber Elongation Factors (SRPPs and REFs) are proteins with major function in rubber particle formation and stabilization. Their corresponding genes are clustered on a scaffold1222 of the reference genomic sequence of the Hevea brasiliensis genome. Apart from gene expression by transcriptomic analyses, to date, no deep analyses have been carried out for the genomic environment of SRPPs and REFs loci. By integrative analyses on transposable element annotation, small RNAs production and gene expression, we analysed their role in the control of the transcription of rubber biosynthetic genes. The first in-depth annotation of TEs (Transposable Elements) and their capacity to produce TE-derived siRNAs (small interfering RNAs) is presented, only possible in the Hevea brasiliensis clone PB 260 for which all data are available. We observed that 11% of genes are located near TEs and their presence may interfere in their transcription at both genetic and epigenetic level. We hypothesized that the genomic environment of rubber biosynthesis genes has been shaped by TE and TE-derived siRNAs with possible transcriptional interference on their gene expression. We discussed possible functionalization of TEs as enhancers and as donors of alternative transcription start sites in promoter sequences, possibly through the modelling of genetic and epigenetic landscapes.

A Reference Genome of Bursaphelenchus mucronatus Provides New Resources for Revealing Its Displacement by Pinewood Nematode.

The Bursaphelenchus mucronatus, which was highly similar with Bursaphelenchus xylophilus in terms of morphological characteristics and biological properties-but had weaker pathogenicity to forests-was a native species often displaced by B. xylophilus when occupying the same niche. Since the draft genome of the invasive B. xylophilus has been published, the absence of a reference genome of B. mucronatus still prevents us from understanding the molecular evidences behind competitive displacement. In this study, we employed Single Molecule, Real-Time (SMRT) sequencing and a Hi-C scaffolding approach to yield a near chromosome-level assembly of B. mucronatus, including six pseudo-chromosomes. The assembly size is 73 Mb, with scaffold N50 of 11.50 Mb and contig N50 of 1.48 Mb. Comparative genomics results showed high similarity between B. xylophilus and B. mucronatus. However, the losing of orphan genes and species-specific orthologous genes in B. mucronatus may indicate weaker adaptability to the environment. The gene family contractions of GPCRs (G Protein-Coupled Receptors) and cellulases in B. mucronatus may jointly contribute to its displacement by B. xylophilus. Overall, we introduced a valuable genomic resource for molecular and evolutionary studies of B. mucronatus, especially for studying the competitive displacement by the pinewood nematode, which could help us control the pathogenicity of pine wilt diseases.

Post-transcriptional regulation of several biological processes involved in latex production in Hevea brasiliensis.

BACKGROUND: Small RNAs modulate plant gene expression at both the transcriptional and post-transcriptional level, mostly through the induction of either targeted DNA methylation or transcript cleavage, respectively. Small RNA networks are involved in specific plant developmental processes, in signaling pathways triggered by various abiotic stresses and in interactions between the plant and viral and non-viral pathogens. They are also involved in silencing maintenance of transposable elements and endogenous viral elements. Alteration in small RNA production in response to various environmental stresses can affect all the above-mentioned processes. In rubber trees, changes observed in small RNA populations in response to trees affected by tapping panel dryness, in comparison to healthy ones, suggest a shift from a transcriptional to a post-transcriptional regulatory pathway. This is the first attempt to characterise small RNAs involved in post-transcriptional silencing and their target transcripts in Hevea. METHODS: Genes producing microRNAs (MIR genes) and loci producing trans-activated small interfering RNA (ta-siRNA) were identified in the clone PB 260 re-sequenced genome. Degradome libraries were constructed with a pool of total RNA from six different Hevea tissues in stressed and non-stressed plants. The analysis of cleaved RNA data, associated with genomics and transcriptomics data, led to the identification of transcripts that are affected by 20-22 nt small RNA-mediated post-transcriptional regulation. A detailed analysis was carried out on gene families related to latex production and in response to growth regulators. RESULTS: Compared to other tissues, latex cells had a higher proportion of transcript cleavage activity mediated by miRNAs and ta-siRNAs. Post-transcriptional regulation was also observed at each step of the natural rubber biosynthesis pathway. Among the genes involved in the miRNA biogenesis pathway, our analyses showed that all of them are expressed in latex. Using phylogenetic analyses, we show that both the Argonaute and Dicer-like gene families recently underwent expansion. Overall, our study underlines the fact that important biological pathways, including hormonal signalling and rubber biosynthesis, are subject to post-transcriptional silencing in laticifers.

AtXRN4 Affects the Turnover of Chosen miRNA*s in Arabidopsis.

Small RNA (sRNA) turnover is a key but poorly understood mechanism that determines the homeostasis of sRNAs. Animal XRN genes contribute the degradation of sRNAs, AtXRN2 and AtXRN3 also contribute the pri-miRNA processing and miRNA loop degradation in plants. However, the possible functions of the plant XRN genes in sRNA degradation are far from known. Here, we find that AtXRN4 contributes the turnover of plant sRNAs in Arabidopsis thaliana mainly by sRNA-seq, qRT-PCR and Northern blot. The mutation of AtXRN4 alters the sRNA profile and the accumulation of 21 nt sRNAs was increased. Some miRNA*s levels are significantly increased in xrn4 mutant plants. However, the accumulation of the primary miRNAs (pri-miRNAs) and miRNA precursors (pre-miRNAs) were generally unchanged in xrn4 mutant plants which indicates that AtXRN4 contributes the degradation of some miRNA*s. Moreover, AtXRN4 interacts with Arabidopsis Argonaute 2 (AtAGO2). This interaction takes place in Processing bodies (P-bodies). Taken together, our observations identified the interaction between XRN4 with AtAGO2 and suggested that plant XRN4 also contributes the turnover of sRNAs.

The genome evolution and domestication of tropical fruit mango.

BACKGROUND: Mango is one of the world's most important tropical fruits. It belongs to the family Anacardiaceae, which includes several other economically important species, notably cashew, sumac and pistachio from other genera. Many species in this family produce family-specific urushiols and related phenols, which can induce contact dermatitis. RESULTS: We generate a chromosome-scale genome assembly of mango, providing a reference genome for the Anacardiaceae family. Our results indicate the occurrence of a recent whole-genome duplication (WGD) event in mango. Duplicated genes preferentially retained include photosynthetic, photorespiration, and lipid metabolic genes that may have provided adaptive advantages to sharp historical decreases in atmospheric carbon dioxide and global temperatures. A notable example of an extended gene family is the chalcone synthase (CHS) family of genes, and particular genes in this family show universally higher expression in peels than in flesh, likely for the biosynthesis of urushiols and related phenols. Genome resequencing reveals two distinct groups of mango varieties, with commercial varieties clustered with India germplasms and demonstrating allelic admixture, and indigenous varieties from Southeast Asia in the second group. Landraces indigenous in China formed distinct clades, and some showed admixture in genomes. CONCLUSIONS: Analysis of chromosome-scale mango genome sequences reveals photosynthesis and lipid metabolism are preferentially retained after a recent WGD event, and expansion of CHS genes is likely associated with urushiol biosynthesis in mango. Genome resequencing clarifies two groups of mango varieties, discovers allelic admixture in commercial varieties, and shows distinct genetic background of landraces.

The apricot (Prunus armeniaca L.) genome elucidates Rosaceae evolution and beta-carotenoid synthesis.

Apricots, scientifically known as Prunus armeniaca L, are drupes that resemble and are closely related to peaches or plums. As one of the top consumed fruits, apricots are widely grown worldwide except in Antarctica. A high-quality reference genome for apricot is still unavailable, which has become a handicap that has dramatically limited the elucidation of the associations of phenotypes with the genetic background, evolutionary diversity, and population diversity in apricot. DNA from P. armeniaca was used to generate a standard, size-selected library with an average DNA fragment size of ~20 kb. The library was run on Sequel SMRT Cells, generating a total of 16.54 Gb of PacBio subreads (N50 = 13.55 kb). The high-quality P. armeniaca reference genome presented here was assembled using long-read single-molecule sequencing at approximately 70× coverage and 171× Illumina reads (40.46 Gb), combined with a genetic map for chromosome scaffolding. The assembled genome size was 221.9 Mb, with a contig NG50 size of 1.02 Mb. Scaffolds covering 92.88% of the assembled genome were anchored on eight chromosomes. Benchmarking Universal Single-Copy Orthologs analysis showed 98.0% complete genes. We predicted 30,436 protein-coding genes, and 38.28% of the genome was predicted to be repetitive. We found 981 contracted gene families, 1324 expanded gene families and 2300 apricot-specific genes. The differentially expressed gene (DEG) analysis indicated that a change in the expression of the 9-cis-epoxycarotenoid dioxygenase (NCED) gene but not lycopene beta-cyclase (LcyB) gene results in a low ß-carotenoid content in the white cultivar "Dabaixing". This complete and highly contiguous P. armeniaca reference genome will be of help for future studies of resistance to plum pox virus (PPV) and the identification and characterization of important agronomic genes and breeding strategies in apricot.

Genome-wide analysis in Hevea brasiliensis laticifers revealed species-specific post-transcriptional regulations of several redox-related genes.

MicroRNA-mediated post-transcriptional regulation has been reported on ROS production and scavenging systems. Although microRNAs first appeared highly conserved among plant species, several aspects of biogenesis, function and evolution of microRNAs were shown to differ. High throughput transcriptome and degradome analyses enable to identify small RNAs and their mRNA targets. A non-photosynthetic tissue particularly prone to redox reactions, laticifers from Hevea brasiliensis, revealed species-specific post-transcriptional regulations. This paper sets out to identify the 407 genes of the thirty main redox-related gene families harboured by the Hevea genome. There are 161 redox-related genes expressed in latex. Thirteen of these redox-related genes were targeted by 11 microRNAs. To our knowledge, this is the first report on a mutation in the miR398 binding site of the cytosolic CuZnSOD. A working model was proposed for transcriptional and post-transcriptional regulation with respect to the predicted subcellular localization of deduced proteins.

RGAAT: A Reference-based Genome Assembly and Annotation Tool for New Genomes and Upgrade of Known Genomes.

The rapid development of high-throughput sequencing technologies has led to a dramatic decrease in the money and time required for de novo genome sequencing or genome resequencing projects, with new genome sequences constantly released every week. Among such projects, the plethora of updated genome assemblies induces the requirement of version-dependent annotation files and other compatible public dataset for downstream analysis. To handle these tasks in an efficient manner, we developed the reference-based genome assembly and annotation tool (RGAAT), a flexible toolkit for resequencing-based consensus building and annotation update. RGAAT can detect sequence variants with comparable precision, specificity, and sensitivity to GATK and with higher precision and specificity than Freebayes and SAMtools on four DNA-seq datasets tested in this study. RGAAT can also identify sequence variants based on cross-cultivar or cross-version genomic alignments. Unlike GATK and SAMtools/BCFtools, RGAAT builds the consensus sequence by taking into account the true allele frequency. Finally, RGAAT generates a coordinate conversion file between the reference and query genomes using sequence variants and supports annotation file transfer. Compared to the rapid annotation transfer tool (RATT), RGAAT displays better performance characteristics for annotation transfer between different genome assemblies, strains, and species. In addition, RGAAT can be used for genome modification, genome comparison, and coordinate conversion. RGAAT is available at https://sourceforge.net/projects/rgaat/ and https://github.com/wushyer/RGAAT_v2 at no cost.

Transcriptome analysis in Hevea brasiliensis latex revealed changes in hormone signalling pathways during ethephon stimulation and consequent Tapping Panel Dryness.

Tapping Panel Dryness (TPD) affects latex production in Hevea brasiliensis. This physiological syndrome involves the agglutination of rubber particles, which leads to partial or complete cessation of latex flow. Latex harvesting consists in tapping soft bark. Ethephon can be applied to stimulate latex flow and its regeneration in laticifers. Several studies have reported transcriptome changes in bark tissues. This study is the first report on deep RNA sequencing of latex to compare the effect of ethephon stimulation and TPD severity. Trees were carefully selected for paired-end sequencing using an Illumina HiSeq 2000. In all, 43 to 60 million reads were sequenced for each treatment in three biological replicates (slight TPD trees without ethephon stimulation, and slight and severe TPD trees with ethephon treatment). Differentially expressed genes were identified and annotated, giving 8,111 and 728 in response to ethephon in slight TPD trees and in ethephon-induced severe TPD trees, respectively. A biological network of responses to ethephon and TPD highlighted the major influence of metabolic processes and the response to stimulus, especially wounding and jasmonate depression in TPD-affected trees induced by ethephon stimulation.

REDO: RNA Editing Detection in Plant Organelles Based on Variant Calling Results.

RNA editing is a post-transcriptional or cotranscriptional process that changes the sequence of the precursor transcript by substitutions, insertions, or deletions. Almost all of the land plants undergo RNA editing in organelles (plastids and mitochondria). Although several software tools have been developed to identify RNA editing events, there has been a great challenge to distinguish true RNA editing events from genome variation, sequencing errors, and other factors. Here we introduce REDO, a comprehensive application tool for identifying RNA editing events in plant organelles based on variant call format files from RNA-sequencing data. REDO is a suite of Perl scripts that illustrate a bunch of attributes of RNA editing events in figures and tables. REDO can also detect RNA editing events in multiple samples simultaneously and identify the significant differential proportion of RNA editing loci. Comparing with similar tools, such as REDItools, REDO runs faster with higher accuracy, and more specificity at the cost of slightly lower sensitivity. Moreover, REDO annotates each RNA editing site in RNAs, whereas REDItools reports only possible RNA editing sites in genome, which need additional steps to obtain RNA editing profiles for RNAs. Overall, REDO can identify potential RNA editing sites easily and provide several functions such as detailed annotations, statistics, figures, and significantly differential proportion of RNA editing sites among different samples.

Rice Expression Database (RED): An integrated RNA-Seq-derived gene expression database for rice.

Rice is one of the most important stable food as well as a monocotyledonous model organism for the plant research community. Here, we present RED (Rice Expression Database; http://expression.ic4r.org), an integrated database of rice gene expression profiles derived entirely from RNA-Seq data. RED features a comprehensive collection of 284 high-quality RNA-Seq experiments, integrates a large number of gene expression profiles and covers a wide range of rice growth stages as well as various treatments. Based on massive expression profiles, RED provides a list of housekeeping and tissue-specific genes and dynamically constructs co-expression networks for gene(s) of interest. Besides, it provides user-friendly web interfaces for querying, browsing and visualizing expression profiles of concerned genes. Together, as a core resource in BIG Data Center, RED bears great utility for characterizing the function of rice genes and better understanding important biological processes and mechanisms underlying complex agronomic traits in rice.