CAFU: a Galaxy framework for exploring unmapped RNA-Seq data.

A widely used approach in transcriptome analysis is the alignment of short reads to a reference genome. However, owing to the deficiencies of specially designed analytical systems, short reads unmapped to the genome sequence are usually ignored, resulting in the loss of significant biological information and insights. To fill this gap, we present Comprehensive Assembly and Functional annotation of Unmapped RNA-Seq data (CAFU), a Galaxy-based framework that can facilitate the large-scale analysis of unmapped RNA sequencing (RNA-Seq) reads from single- and mixed-species samples. By taking advantage of machine learning techniques, CAFU addresses the issue of accurately identifying the species origin of transcripts assembled using unmapped reads from mixed-species samples. CAFU also represents an innovation in that it provides a comprehensive collection of functions required for transcript confidence evaluation, coding potential calculation, sequence and expression characterization and function annotation. These functions and their dependencies have been integrated into a Galaxy framework that provides access to CAFU via a user-friendly interface, dramatically simplifying complex exploration tasks involving unmapped RNA-Seq reads. CAFU has been validated with RNA-Seq data sets from wheat and Zea mays (maize) samples. CAFU is freely available via GitHub: https://github.com/cma2015/CAFU.

Evolution of the RNA N 6-Methyladenosine Methylome Mediated by Genomic Duplication.

RNA N 6-methyladenosine (m6A) modification is the most abundant form of RNA epigenetic modification in eukaryotes. Given that m6A evolution is associated with the selective constraints of nucleotide sequences in mammalian genomes, we hypothesize that m6A evolution can be linked, at least in part, to genomic duplication events in complex polyploid plant genomes. To test this hypothesis, we presented the maize (Zea mays) m6A modification landscape in a transcriptome-wide manner and identified 11,968 m6A peaks carried by 5,893 and 3,811 genes from two subgenomes (maize1 and maize2, respectively). Each of these subgenomes covered over 2,200 duplicate genes. Within these duplicate genes, those carrying m6A peaks exhibited significant differences in retention rate. This biased subgenome fractionation of m6A-methylated genes is associated with multiple sequence features and is influenced by asymmetric evolutionary rates. We also characterized the coevolutionary patterns of m6A-methylated genes and transposable elements, which can be mediated by whole genome duplication and tandem duplication. We revealed the evolutionary conservation and divergence of duplicated m6A functional factors and the potential role of m6A modification in maize responses to drought stress. This study highlights complex interplays between m6A modification and gene duplication, providing a reference for understanding the mechanisms underlying m6A evolution mediated by genome duplication events.

Evolution of intron-poor clades and expression patterns of the glycosyltransferase family 47.

MAIN CONCLUSION: A large-scale bioinformatics analysis revealed the origin and evolution of GT47 gene family, and identified two clades of intron-poor genes with putative functions in drought stress responses and seed development in maize. Glycosyltransferase family 47 (GT47) genes encode ß-galactosyltransferases and ß-glucuronyltransferases that synthesize pectin, xyloglucans and xylan, which are important components of the plant cell wall. In this study, we performed a systematic and large-scale bioinformatics analysis of GT47 gene family using 352 GT47 proteins from 15 species ranging from cyanobacteria to seed plants. The analysis results showed that GT47 family may originate in cyanobacteria and expand along the evolutionary trajectory to moss. Further analysis of 47 GT47 genes in maize revealed that they can divide into five clades with diverse exon-intron structures. Among these five clades, two were mainly composed with intron-poor genes, which may originate in the moss. Gene duplication analysis revealed that the expansion of GT47 gene family in maize was significantly driven from tandem duplication events and segmental duplication events. Significantly, almost all duplicated genes are intron-poor genes. Expression analysis indicated that several intron-poor GT47 genes may be involved in the drought stress response and seed development in maize. This work provides insight into the origin and evolutionary process, expansion mechanisms and expression patterns of GT47 genes, thus facilitating their functional investigations in the future.

Effect of choline amino acid ionic liquids on maize seed germination and endogenous plant hormone levels.

Prior research has established choline-based ionic liquids (ILs) as safe for various organisms. However, their impact on plants has been underexplored. To identify effective eco-friendly ILs, we synthesized seven choline amino acid ([Chl][AA]) ILs and analyzed their physiological influence on maize seed germination. In contrast to the traditionally used N-octyl pyridinium bromide IL, these seven [Chl][AA] ILs exhibited substantially lower toxicity. Moreover, within a broad treatment concentration range (10-100 mg L-1), these ILs notably enhanced maize germination indices and root and shoot growth. Specifically, treatment with 100 mg L-1 choline tryptophan resulted in a 21.2% increase in germination index compared to those of control maize. Compared to the control, the application of choline serine, choline aspartic acid, choline phenylalanine, and choline tryptophan at 100 mg L-1 led to respective increases of 23.9%, 21.5%, 22.5%, and 24.5% in maize shoot length. Analysis of endogenous hormones and free amino acid contents revealed elevated levels of growth-promoting plant hormones (gibberellic acid and zeatin) in maize shoot tips, as well as increased contents of major amino acids (glutamate, glycine, and arginine) following treatment with different [Chl][AA] ILs at 100 mg L-1. These findings indicate that [Chl][AA] holds promise for the development and application of novel low-toxicity ILs.

Effects of pyridinium-based ionic liquids with different alkyl chain lengths on the growth of maize seedlings.

The effects of four water-soluble pyridinium-based ionic liquids (ILs) differing in alkyl side chain length, namely, N-ethyl pyridinium bromide ([EPy]Br), N-butyl pyridinium bromide ([BPy]Br), N-hexyl pyridinium bromide ([HPy]Br), and N-octyl pyridinium bromide ([OPy]Br), on the growth of maize seedlings were investigated for the first time. The results revealed that the phytotoxicity of these ILs was significantly correlated with their concentration and alkyl side chain length. The 8-day 50% inhibition values indicated that the toxicity increased as the length of the alkyl chain increased: [EPy]Br < [BPy]Br < [HPy]Br < [OPy]Br. In addition, root growth was found to be more sensitive to ILs than stem growth. In response to exposure to ILs of increasing concentration, we observed different trends in the pigment contents and specific antioxidant enzyme activities in maize seedlings, whereas the contents of malondialdehyde were significantly increased. In addition, RNA sequencing analysis, performed to examine the gene expression profiles of maize leaves under [HPy]Br and [OPy]Br treatments, revealed that a larger number of genes were differentially expressed in response to [OPy]Br treatment. Furthermore, pathway enrichment analysis revealed that both [HPy]Br and [OPy]Br treatments, and particularly the latter, caused a down-regulation of genes involved in photosynthesis and carbohydrate and nitrogen metabolism. Our findings thus indicate that pyridinium-based IL toxicity might be associated with oxidative stress and changes in gene expression profiles.

Strong positive selection drives rapid diversification of R-genes in Arabidopsis relatives.

Although plant resistance (R) genes are extremely diverse and evolve rapidly, little is known about the mechanisms that generate this sequence divergence. To investigate these forces, we compared all nucleotide binding sites and leucine-rich repeat R-genes between two closely related species, Arabidopsis thaliana and Arabidopsis lyrata. Our analyses revealed two distinct evolutionary patterns driven by either positive or stabilizing selection. Most R-genes (>50%) were evolving under strong positive selection characterized by high Ka/Ks ratios (>1), frequent recombination, copy number variation, and extremely high sequence divergence between the two species. The stably selected R-genes (<30%) have exactly the opposite four characters as the positively selected genes. The remaining R-genes (about 20%) are present in only one genome and absent from the other. A higher proportion of such genes were found to be part of TNL class (23.5%) compared to the non-TNL class (5.6%), suggesting different evolutionary patterns between these two groups. A significant correlation between Ka and divergence was revealed, indicating that the rapid evolution and diversification of R-genes were initiated by selectively generated, frequently shuffled and selectively maintained non-synonymous substitutions. Our genome-wide analyses confirmed an amazing mechanism by which plants to selectively accumulate and efficiently exploit these non-synonymous substitutions for their resistance to various pathogens.

Structural and expressional analysis of the B-hordein genes in Tibetan hull-less barley.

The B-hordein gene family was analyzed from two Tibetan hull-less barley cultivars Z09 and Z26 (Hordeum vulgare subsp. vulgare). Fourteen B-hordein genes, designated BZ09-2 to BZ09-5 (from Z09) and BZ26-1 to BZ26-10 (from Z26), were sequenced. Seven of them, similar to a previously reported BZ09-1 from Z09, were predicted to encode putative active proteins each with a signal peptide, a repetitive domain, and a C-terminal region; seven of them were predicted to be pseudogenes. The B-hordein gene family was analyzed using all known representatives of B-hordein sequences and two most similar LMW-GSs of Triticum aestivum. Alignment of these seven putative proteins with known B-hordeins and two most similar LMW-GSs of T. aestivum revealed that they shared a common motif. A large variation was observed between numbers of repeat units of predicted B-hordeins of Z26 and Z09. Phylogenetic analysis revealed that all BZ26 clones were clustered in a subgroup, and BZ09-1 formed another subgroup by itself in the putative eight active genes. In addition, six 5'-upstream regulatory sequences of the B-hordein genes were isolated from the two accessions by a single oligonucleotide nested PCR, and several different mutations were identified in the cis-acting element GLM and two distinctive sequences (Z09P-2 and Z26P-3). Phylogenetic analysis of 5'-upstream regulatory regions of the B-hordein genes showed that members from the same accession were clustered together except for two distinct members. Quantitative real time PCR analysis indicated distinct expression levels of B-hordein genes in four developing stages of developing grains in two accessions. These findings suggest B-hordein genes have intrinsic differences between accessions, and this knowledge will be useful for incorporating the B-hordeins protein in barley breeding programs.

Transcriptomic analysis of the phytotoxic effects of 1-allyl-3-methylimidazolium chloride on the growth and plant hormone metabolic pathways of maize (Zea mays L.) seedlings.

In this study, we investigated the phytotoxicity of an imidazolium-based ionic liquid, 1-allyl-3-methylimidazolium chloride ([Amim]Cl), against maize seedlings. It was found that in response to an increase in [Amim]Cl treatment concentrations, there were significant decreases in growth parameters (fresh weights and lengths) and the photosynthetic pigment contents of maize plants, whereas in contrast, the malondialdehyde content increased. In order to determine the molecular basis of [Amim]Cl-induced plant growth inhibition, an RNA-Seq analysis to examine the gene expression profiles of selected central biological pathways was performed. And a total of 4024 genes that were differentially expressed between control and 400 mg/L [Amim]Cl-treated plants were accordingly identified. Pathway enrichment analysis for the differentially expressed genes revealed that 12 of 15 genes in the porphyrin and chlorophyll metabolic pathways were down-regulated in response to [Amim]Cl treatment. Moreover, all six genes encoding key chlorophyll synthetic enzymes were down-regulated by [Amim]Cl. With regards to plant hormone metabolic pathways, the genes encoding key enzymes involved in ethybilene and abscisic acid (ABA) biosynthesis were up-regulated in response to [Amim]Cl treatment. Genes responsible for gibberellin (GA) inactivation were also stimulated by [Amim]Cl. These observations indicate that [Amim]Cl may promote the biosynthesis of senescence-related hormones (ethylene and ABA) as well as inactivation of growth-promoting hormones (GAs). It might be concluded that the observed [Amim]Cl-induced inhibition of maize seedling growth could be associated with changes in the gene expression profiles of these metabolic pathways.

Subtle Perturbations of the Maize Methylome Reveal Genes and Transposons Silenced by Chromomethylase or RNA-Directed DNA Methylation Pathways.

DNA methylation is a chromatin modification that can provide epigenetic regulation of gene and transposon expression. Plants utilize several pathways to establish and maintain DNA methylation in specific sequence contexts. The chromomethylase (CMT) genes maintain CHG (where H = A, C or T) methylation. The RNA-directed DNA methylation (RdDM) pathway is important for CHH methylation. Transcriptome analysis was performed in a collection of Zea mays lines carrying mutant alleles for CMT or RdDM-associated genes. While the majority of the transcriptome was not affected, we identified sets of genes and transposon families sensitive to context-specific decreases in DNA methylation in mutant lines. Many of the genes that are up-regulated in CMT mutant lines have high levels of CHG methylation, while genes that are differentially expressed in RdDM mutants are enriched for having nearby mCHH islands, implicating context-specific DNA methylation in the regulation of expression for a small number of genes. Many genes regulated by CMTs exhibit natural variation for DNA methylation and transcript abundance in a panel of diverse inbred lines. Transposon families with differential expression in the mutant genotypes show few defining features, though several families up-regulated in RdDM mutants show enriched expression in endosperm tissue, highlighting the potential importance for this pathway during reproduction. Taken together, our findings suggest that while the number of genes and transposon families whose expression is reproducibly affected by mild perturbations in context-specific methylation is small, there are distinct patterns for loci impacted by RdDM and CMT mutants.

Heritable Epigenomic Changes to the Maize Methylome Resulting from Tissue Culture.

DNA methylation can contribute to the maintenance of genome integrity and regulation of gene expression. In most situations, DNA methylation patterns are inherited quite stably. However, changes in DNA methylation can occur at some loci as a result of tissue culture resulting in somaclonal variation. To investigate heritable epigenetic changes as a consequence of tissue culture, a sequence-capture bisulfite sequencing approach was implemented to monitor context-specific DNA methylation patterns in ∼15 Mb of the maize genome for a population of plants that had been regenerated from tissue culture. Plants that have been regenerated from tissue culture exhibit gains and losses of DNA methylation at a subset of genomic regions. There was evidence for a high rate of homozygous changes to DNA methylation levels that occur consistently in multiple independent tissue culture lines, suggesting that some loci are either targeted or hotspots for epigenetic variation. The consistent changes inherited following tissue culture include both gains and losses of DNA methylation and can affect CG, CHG, or both contexts within a region. Only a subset of the tissue culture changes observed in callus plants are observed in the primary regenerants, but the majority of DNA methylation changes present in primary regenerants are passed onto offspring. This study provides insights into the susceptibility of some loci and potential mechanisms that could contribute to altered DNA methylation and epigenetic state that occur during tissue culture in plant species.

A systems approach to a spatio-temporal understanding of the drought stress response in maize.

Crops are often subjected to periods of drought stress during their life cycle. However, how stress response mechanisms contribute to the crosstalk between stress signaling pathways and developmental signaling pathways is still unknown. We built a gene co-expression network from a spatio-temporal transcriptomic map of the drought stress response in maize (Zea mays), profiled from three tissues and four developmental stages and characterized hub genes associated with duplication events, selection, and regulatory networks. Co-expression analysis grouped drought-response genes into ten modules, covering 844 highly connected genes (hub genes). Of these, 15.4% hub genes had diverged by whole-genome duplication events and 2.5% might then have been selected during natural domestication and artificial improvement processes, successively. We identified key transcription factor hubs in a transcriptional regulatory network, which may function as a crosstalk mechanism between drought stress and developmental signalling pathways in maize. Understanding the evolutionary biases that have evolved to enhance drought adaptation lays the foundation for further dissection of crosstalk between stress signalling pathways and developmental signalling pathways in maize, towards molecular design of new cultivars with desirable yield and greater stress tolerance.

Cloning and characterization of four B-hordein genes from Tibetan hull-less barley (Hordeum vulgare subsp. vulgare).

Four B-hordein genes, designated BH1-BH4, were cloned using PCR amplification from two hull-less barley cultivars, ZQ7239 and ZQ148, collected from Tibet. The results of sequencing indicated that BH1-BH4 contained complete open reading frames (ORFs). Comparison of their predicted polypeptide sequences with the published sequences suggested that they all share the same basic protein structure. Phylogenetic analysis indicated that the deduced amino-acid sequences of BH1-BH4 genes were more closely related to B-hordeins from cultivated barley (Hordeum vulgare L.) than to any other prolamins from wild barley and Aegilops tauschii. Comparison of the coding regions of BH1-BH4 genes showed that BH1 had a lower sequence identity to other previously published B-hordeins than the other three B-hordeins obtained in this study. BH1 was then cloned in a bacterial expression vector based on bacteriophage T7 RNA polymerase. The resulting plasmid produced a 28.15 kDa protein in Escherichia coli. The potential value of B-hordein genes in grain quality improvement of hull-less barley has been discussed.

Identification of maize long non-coding RNAs responsive to drought stress.

Long non-coding RNAs (lncRNAs) represent a class of riboregulators that either directly act in long form or are processed to shorter miRNAs and siRNAs. Emerging evidence shows that lncRNAs participate in stress responsive regulation. In this study, to identify the putative maize lncRNAs responsive to drought stress, 8449 drought responsive transcripts were first uploaded to the Coding Potential Calculator website for classification as protein coding or non-coding RNAs, and 1724 RNAs were identified as potential non-coding RNAs. A Perl script was written to screen these 1724 ncRNAs and 664 transcripts were ultimately identified as drought-responsive lncRNAs. Of these 664 transcripts, 126 drought-responsive lncRNAs were highly similar to known maize lncRNAs; the remaining 538 transcripts were considered as novel lncRNAs. Among the 664 lncRNAs identified as drought responsive, 567 were upregulated and 97 were downregulated in drought-stressed leaves of maize. 8 lncRNAs were identified as miRNA precursor lncRNAs, 62 were classified as both shRNA and siRNA precursors, and 279 were classified as siRNA precursors. The remaining 315 lncRNAs were classified as other lncRNAs that are likely to function as longer molecules. Among these 315 lncRNAs, 10 are identified as antisense lncRNAs and 7 could pair with 17 CDS sequences with near-perfect matches. Finally, RT-qPCR results confirmed that all selected lncRNAs could respond to drought stress. These findings extend the current view on lncRNAs as ubiquitous regulators under stress conditions.

Functional analysis of the 5' regulatory region of the maize GALACTINOL SYNTHASE2 gene.

Maize (Zea mays) GALACTINOL SYNTHASE (GolS) is a key enzyme in the raffinose biosynthetic pathway. We have previously characterized the maize GolS2 (ZmGolS2) gene as heat shock induced in maize germinating seeds and cultured cells. Here we report the identification, isolation and characterization of the 1.5Kb, 5' regulatory region of the ZmGolS2 gene. The 1.5kb fragment and its deletions were tested for promoter activity by their regulation of the Renilla (Renilla reniformis) luciferase reporter gene expression in maize protoplasts cultured at either 25°C or 42°C for 24h. The expression of a constitutively expressed firefly (Photinus ssp.) luciferase gene in the same vector was used as a reference. One heat shock element (HSE) was identified by comparing the promoter activity of each fragment under normal and heat shock conditions. Deletion or triplication of this HSE motif, abolished or enhanced the heat shock response of the ZmGolS2 promoter, respectively. This HSE motif is specifically bound by proteins in the nuclear extracts of heat shock stressed, but not unstressed maize cells as confirmed by DNA-EMSA. This work helps to understand the regulatory mechanism of the ZmGolS2 gene under stress conditions.