A smart RBS library and its prediction model for robust and accurate fine-tuning of gene expression in Bacillus species.

Bacillus species, such as Bacillus subtilis and Bacillus licheniformis, are important industrial bacteria. However, there is a lack of standardized and predictable genetic tools for convenient and reproducible assembly of genetic modules in Bacillus species to realize their full potential. In this study, we constructed a Ribosome Binding Site (RBS) library in B. licheniformis, which provides incremental regulation of expression levels over a 104-fold range. Additionally, we developed a model to quantify the resulting translation rates. We successfully demonstrated the robust expression of various target genes using the RBS library and showed that the model accurately predicts the translation rates of arbitrary coding genes. Importantly, we also extended the use of the RBS library and prediction model to B. subtilis, B. thuringiensis, and B. amyloliquefacie. The versatility of the RBS library and its prediction model enables quantification of biological behavior, facilitating reliable forward engineering of gene expression.

Abscisic acid regulates secondary cell-wall formation and lignin deposition in Arabidopsis thaliana through phosphorylation of NST1.

Plant secondary cell-wall (SCW) deposition and lignification are affected by both seasonal factors and abiotic stress, and these responses may involve the hormone abscisic acid (ABA). However, the mechanisms involved are not clear. Here we show that mutations that limit ABA synthesis or signaling reduce the extent of SCW thickness and lignification in Arabidopsis thaliana through the core ABA-signaling pathway involving SnRK2 kinases. SnRK2.2. 3 and 6 physically interact with the SCW regulator NAC SECONDARY WALL THICKENING PROMOTING FACTOR 1 (NST1), a NAC family transcription factor that orchestrates the transcriptional activation of a suite of downstream SCW biosynthesis genes, some of which are involved in the biosynthesis of cellulose and lignin. This interaction leads to phosphorylation of NST1 at Ser316, a residue that is highly conserved among NST1 proteins from dicots, but not monocots, and is required for transcriptional activation of downstream SCW-related gene promoters. Loss of function of NST1 in the snd1 mutant background results in lack of SCWs in the interfascicular fiber region of the stem, and the Ser316Ala mutant of NST1 fails to complement this phenotype and ABA-induced lignin pathway gene expression. The discovery of NST1 as a key substrate for phosphorylation by SnRK2 suggests that the ABA-mediated core-signaling cascade provided land plants with a hormone-modulated, competitive desiccation-tolerance strategy allowing them to differentiate water-conducting and supporting tissues built of cells with thicker cell walls.

Multifeature analyses of vascular cambial cells reveal longevity mechanisms in old Ginkgo biloba trees.

Aging is a universal property of multicellular organisms. Although some tree species can live for centuries or millennia, the molecular and metabolic mechanisms underlying their longevity are unclear. To address this, we investigated age-related changes in the vascular cambium from 15- to 667-y-old Ginkgo biloba trees. The ring width decreased sharply during the first 100 to 200 y, with only a slight change after 200 y of age, accompanied by decreasing numbers of cambial cell layers. In contrast, average basal area increment (BAI) continuously increased with aging, showing that the lateral meristem can retain indeterminacy in old trees. The indole-3-acetic acid (IAA) concentration in cambial cells decreased with age, whereas the content of abscisic acid (ABA) increased significantly. In addition, cell division-, cell expansion-, and differentiation-related genes exhibited significantly lower expression in old trees, especially miR166 and HD-ZIP III interaction networks involved in cambial activity. Disease resistance-associated genes retained high expression in old trees, along with genes associated with synthesis of preformed protective secondary metabolites. Comprehensive evaluation of the expression of genes related to autophagy, senescence, and age-related miRNAs, together with analysis of leaf photosynthetic efficiencies and seed germination rates, demonstrated that the old trees are still in a healthy, mature state, and senescence is not manifested at the whole-plant level. Taken together, our results reveal that long-lived trees have evolved compensatory mechanisms to maintain a balance between growth and aging processes. This involves continued cambial divisions, high expression of resistance-associated genes, and continued synthetic capacity of preformed protective secondary metabolites.

ARABIDOPSIS DEHISCENCE ZONE POLYGALACTURONASE 1 (ADPG1) releases latent defense signals in stems with reduced lignin content.

There is considerable interest in engineering plant cell wall components, particularly lignin, to improve forage quality and biomass properties for processing to fuels and bioproducts. However, modifying lignin content and/or composition in transgenic plants through down-regulation of lignin biosynthetic enzymes can induce expression of defense response genes in the absence of biotic or abiotic stress. Arabidopsis thaliana lines with altered lignin through down-regulation of hydroxycinnamoyl CoA:shikimate/quinate hydroxycinnamoyl transferase (HCT) or loss of function of cinnamoyl CoA reductase 1 (CCR1) express a suite of pathogenesis-related (PR) protein genes. The plants also exhibit extensive cell wall remodeling associated with induction of multiple cell wall-degrading enzymes, a process which renders the corresponding biomass a substrate for growth of the cellulolytic thermophile Caldicellulosiruptor bescii lacking a functional pectinase gene cluster. The cell wall remodeling also results in the release of size- and charge-heterogeneous pectic oligosaccharide elicitors of PR gene expression. Genetic analysis shows that both in planta PR gene expression and release of elicitors are the result of ectopic expression in xylem of the gene ARABIDOPSIS DEHISCENCE ZONE POLYGALACTURONASE 1 (ADPG1), which is normally expressed during anther and silique dehiscence. These data highlight the importance of pectin in cell wall integrity and the value of lignin modification as a tool to interrogate the informational content of plant cell walls.

Identification of Pueraria spp. through DNA barcoding and comparative transcriptomics.

BACKGROUND: Kudzu is a term used generically to describe members of the genus Pueraria. Kudzu roots have been used for centuries in traditional Chinese medicine in view of their high levels of beneficial isoflavones including the unique 8-C-glycoside of daidzein, puerarin. In the US, kudzu is seen as a noxious weed causing ecological and economic damage. However, not all kudzu species make puerarin or are equally invasive. Kudzu remains difficult to identify due to its diverse morphology and inconsistent nomenclature. RESULTS: We have generated sequences for the internal transcribed spacer 2 (ITS2) and maturase K (matK) regions of Pueraria montana lobata, P. montana montana, and P. phaseoloides, and identified two accessions previously used for differential analysis of puerarin biosynthesis as P. lobata and P. phaseoloides. Additionally, we have generated root transcriptomes for the puerarin-producing P. m. lobata and the non-puerarin producing P. phaseoloides. Within the transcriptomes, microsatellites were identified to aid in species identification as well as population diversity. CONCLUSIONS: The barcode sequences generated will aid in fast and efficient identification of the three kudzu species. Additionally, the microsatellites identified from the transcriptomes will aid in genetic analysis. The root transcriptomes also provide a molecular toolkit for comparative gene expression analysis towards elucidation of the biosynthesis of kudzu phytochemicals.

Leaf layer-based transcriptome profiling for discovery of epidermal-selective promoters in Medicago truncatula.

MAIN CONCLUSION: Transcriptomics of manually dissected leaf layers from Medicago truncatula identifies genes with preferential expression in upper and/or lower epidermis. The promoters of these genes confer epidermal-specific expression of transgenes. Improving the quality and quantity of proanthocyanidins (PAs) in forage legumes has potential to improve the nitrogen nutrition of ruminant animals and protect them from the risk of pasture bloat, as well as parasites. However, ectopic constitutive accumulation of PAs in plants by genetic engineering can significantly inhibit growth. We selected the leaf epidermis as a candidate tissue for targeted engineering of PAs or other pathways. To identify gene promoters selectively expressed in epidermal tissues, we performed comparative transcriptomic analyses in the model legume Medicago truncatula, using five tissue samples representing upper epidermis, lower epidermis, whole leaf without upper epidermis, whole leaf without lower epidermis, and whole leaf. We identified 52 transcripts preferentially expressed in upper epidermis, most of which encode genes involved in flavonoid biosynthesis, and 53 transcripts from lower epidermis, with the most enriched category being anatomical structure formation. Promoters of the preferentially expressed genes were cloned from the M. truncatula genome and shown to direct tissue-selective promoter activities in transient assays. Expression of the PA pathway transcription factor TaMYB14 under control of several of the promoters in transgenic alfalfa resulted in only modest MYB14 transcript accumulation and low levels of PA production. Activity of a subset of promoters was confirmed by transcript analysis in field-grown alfalfa plants throughout the growing season, and revealed variable but consistent expression, which was generally highest 3-4 weeks after cutting. We conclude that, although the selected promoters show acceptable tissue-specificity, they may not drive high enough transcription factor expression to activate the PA pathway.

Dissecting the transcriptional regulation of proanthocyanidin and anthocyanin biosynthesis in soybean (Glycine max).

Proanthocyanidins (PAs), also known as condensed tannins, are plant natural products that are beneficial for human and livestock health. As one of the largest grown crops in the world, soybean (Glycine max) is widely used as human food and animal feed. Many cultivated soybeans with yellow seed coats lack PAs or anthocyanins, although some soybean cultivars have coloured seed coats that contain these compounds. Here, we analyse the transcriptional control of PA and anthocyanin biosynthesis in soybean. Ectopic expression of the transcription factors (TFs) GmTT2A, GmTT2B, GmMYB5A or R in soybean hairy roots induced the accumulation of PAs (primarily in phloem tissues) or anthocyanins and led to up-regulation of 1775, 856, 1411 and 1766 genes, respectively, several of which encode enzymes involved in PA biosynthesis. The genes regulated by GmTT2A and GmTT2B partially overlapped, suggesting conserved but potentially divergent roles for these two TFs in regulating PA accumulation in soybean. The two key enzymes anthocyanidin reductase and leucoanthocyanidin reductase were differentially upregulated, by GmTT2A/GmTT2B and GmMYB5A, respectively. Transgenic soybean plants overexpressing GmTT2B or MtLAP1 (a proven up-regulator of the upstream reactions for production of precursors for PA biosynthesis in legumes) showed increased accumulation of PAs and anthocyanins, respectively, associated with transcriptional reprogramming paralleling the RNA-seq data collected in soybean hairy roots. Collectively, our results show that engineered PA biosynthesis in soybean exhibits qualitative and spatial differences from the better-studied model systems Arabidopsis thaliana and Medicago truncatula, and suggest targets for engineering PAs in soybean plants.

Growth-defense trade-offs and yield loss in plants with engineered cell walls.

As a major component of plant secondary cell walls, lignin provides structural integrity and rigidity, and contributes to primary defense by providing a physical barrier to pathogen ingress. Genetic modification of lignin biosynthesis has been adopted to reduce the recalcitrance of lignified cell walls to improve biofuel production, tree pulping properties and forage digestibility. However, lignin-modification is often, but unpredictably, associated with dwarf phenotypes. Hypotheses suggested to explain this include: collapsed vessels leading to defects in water and solute transport; accumulation of molecule(s) that are inhibitory to plant growth or deficiency of metabolites that are critical for plant growth; activation of defense pathways linked to cell wall integrity sensing. However, there is still no commonly accepted underlying mechanism for the growth defects. Here, we discuss recent data on transcriptional reprogramming in plants with modified lignin content and their corresponding suppressor mutants, and evaluate growth-defense trade-offs as a factor underlying the growth phenotypes. New approaches will be necessary to estimate how gross changes in transcriptional reprogramming may quantitatively affect growth. Better understanding of the basis for yield drag following cell wall engineering is important for the biotechnological exploitation of plants as factories for fuels and chemicals.

Enzymatic basis for C-lignin monomer biosynthesis in the seed coat of Cleome hassleriana.

C-lignin is a linear polymer of caffeyl alcohol, found in the seed coats of several exotic plant species, with promising properties for generation of carbon fibers and high value chemicals. In the ornamental plant Cleome hassleriana, guaiacyl (G) lignin is deposited in the seed coat for the first 6-12 days after pollination, after which G-lignin deposition ceases and C-lignin accumulates, providing an excellent model system to study C-lignin biosynthesis. We performed RNA sequencing of seed coats harvested at 2-day intervals throughout development. Bioinformatic analysis identified a complete set of lignin biosynthesis genes for Cleome. Transcript analysis coupled with kinetic analysis of recombinant enzymes in Escherichia coli revealed that the switch to C-lignin formation was accompanied by down-regulation of transcripts encoding functional caffeoyl CoA- and caffeic acid 3-O-methyltransferases (CCoAOMT and COMT) and a form of cinnamyl alcohol dehydrogenase (ChCAD4) with preference for coniferaldehyde as substrate, and up-regulation of a form of CAD (ChCAD5) with preference for caffealdehyde. Based on these analyses, blockage of lignin monomer methylation by down-regulation of both O-methyltransferases (OMTs) and methionine synthase (for provision of C1 units) appears to be the major factor in diversion of flux to C-lignin in the Cleome seed coat, although the change in CAD specificity also contributes based on the reduction of C-lignin levels in transgenic Cleome with down-regulation of ChCAD5. Structure modeling and mutational analysis identified amino acid residues important for the preference of ChCAD5 for caffealdehyde.

Ectopic Defense Gene Expression Is Associated with Growth Defects in Medicago truncatula Lignin Pathway Mutants.

Lignin provides essential mechanical support for plant cell walls but decreases the digestibility of forage crops and increases the recalcitrance of biofuel crops. Attempts to modify lignin content and/or composition by genetic modification often result in negative growth effects. Although several studies have attempted to address the basis for such effects in individual transgenic lines, no common mechanism linking lignin modification with perturbations in plant growth and development has yet been identified. To address whether a common mechanism exists, we have analyzed transposon insertion mutants resulting in independent loss of function of five enzymes of the monolignol pathway, as well as one double mutant, in the model legume Medicago truncatula These plants exhibit growth phenotypes from essentially wild type to severely retarded. Extensive phenotypic, transcriptomic, and metabolomics analyses, including structural characterization of differentially expressed compounds, revealed diverse phenotypic consequences of lignin pathway perturbation that were perceived early in plant development but were not predicted by lignin content or composition alone. Notable phenotypes among the mutants with severe growth impairment were increased trichome numbers, accumulation of a variety of triterpene saponins, and extensive but differential ectopic expression of defense response genes. No currently proposed model explains the observed phenotypes across all lines. We propose that reallocation of resources into defense pathways is linked to the severity of the final growth phenotype in monolignol pathway mutants of M. truncatula, although it remains unclear whether this is a cause or an effect of the growth impairment.

SINGLE FLOWER TRUSS and SELF-PRUNING signal developmental and metabolic networks to guide cotton architectures.

Patterns of indeterminate and determinate growth specify plant architecture and influence crop productivity. In cotton (Gossypium hirsutum), SINGLE FLOWER TRUSS (SFT) stimulates the transition to flowering and determinate growth, while its closely related antagonist SELF-PRUNING (SP) maintains meristems in indeterminate states to favor vegetative growth. Overexpressing GhSFT while simultaneously silencing GhSP produces highly determinate cotton with reduced foliage and synchronous fruiting. These findings suggest that GhSFT, GhSP, and genes in these signaling networks hold promise for enhancing 'annualized' growth patterns and improving cotton productivity and management. To identify the molecular programs underlying cotton growth habits, we used comparative co-expression networks, differential gene expression, and phenotypic analyses in cotton varieties expressing altered levels of GhSFT or GhSP. Using multiple cotton and tomato datasets, we identified diverse genetic modules highly correlated with SFT or SP orthologs which shared related Gene Ontologies in different crop species. Notably, altering GhSFT or GhSP levels in cotton affected the expression of genes regulating meristem fate and metabolic pathways. Further phenotypic analyses of gene products involved in photosynthesis, secondary metabolism, and cell wall biosynthesis showed that early changes in GhSFT and GhSP levels profoundly impacted later development in distal tissues. Identifying the molecular underpinnings of GhSFT and GhSP activities emphasizes their broad actions in regulating cotton architecture.

Gene regulatory networks for lignin biosynthesis in switchgrass (Panicum virgatum).

Cell wall recalcitrance is the major challenge to improving saccharification efficiency in converting lignocellulose into biofuels. However, information regarding the transcriptional regulation of secondary cell wall biogenesis remains poor in switchgrass (Panicum virgatum), which has been selected as a biofuel crop in the United States. In this study, we present a combination of computational and experimental approaches to develop gene regulatory networks for lignin formation in switchgrass. To screen transcription factors (TFs) involved in lignin biosynthesis, we developed a modified method to perform co-expression network analysis using 14 lignin biosynthesis genes as bait (target) genes. The switchgrass lignin co-expression network was further extended by adding 14 TFs identified in this study, and seven TFs identified in previous studies, as bait genes. Six TFs (PvMYB58/63, PvMYB42/85, PvMYB4, PvWRKY12, PvSND2 and PvSWN2) were targeted to generate overexpressing and/or down-regulated transgenic switchgrass lines. The alteration of lignin content, cell wall composition and/or plant growth in the transgenic plants supported the role of the TFs in controlling secondary wall formation. RNA-seq analysis of four of the transgenic switchgrass lines revealed downstream target genes of the secondary wall-related TFs and crosstalk with other biological pathways. In vitro transactivation assays further confirmed the regulation of specific lignin pathway genes by four of the TFs. Our meta-analysis provides a hierarchical network of TFs and their potential target genes for future manipulation of secondary cell wall formation for lignin modification in switchgrass.

An Improved Brome mosaic virus Silencing Vector: Greater Insert Stability and More Extensive VIGS.

Virus-induced gene silencing (VIGS) is used extensively for gene function studies in plants. VIGS is inexpensive and rapid compared with silencing conducted through stable transformation, but many virus-silencing vectors, especially in grasses, induce only transient silencing phenotypes. A major reason for transient phenotypes is the instability of the foreign gene fragment (insert) in the vector during VIGS. Here, we report the development of a Brome mosaic virus (BMV)-based vector that better maintains inserts through modification of the original BMV vector RNA sequence. Modification of the BMV RNA3 sequence yielded a vector, BMVCP5, that better maintained phytoene desaturase and heat shock protein70-1 (HSP70-1) inserts in Nicotiana benthamiana and maize (Zea mays). Longer maintenance of inserts was correlated with greater target gene silencing and more extensive visible silencing phenotypes displaying greater tissue penetration and involving more leaves. The modified vector accumulated similarly to the original vector in N. benthamiana after agroinfiltration, thus maintaining a high titer of virus in this intermediate host used to produce virus inoculum for grass hosts. For HSP70, silencing one family member led to a large increase in the expression of another family member, an increase likely related to the target gene knockdown and not a general effect of virus infection. The cause of the increased insert stability in the modified vector is discussed in relationship to its recombination and accumulation potential. The modified vector will improve functional genomic studies in grasses, and the conceptual methods used to improve the vector may be applied to other VIGS vectors.

Co-expression networks for plant biology: why and how.

Co-expression network analysis is one of the most powerful approaches for interpretation of large transcriptomic datasets. It enables characterization of modules of co-expressed genes that may share biological functional linkages. Such networks provide an initial way to explore functional associations from gene expression profiling and can be applied to various aspects of plant biology. This review presents the applications of co-expression network analysis in plant biology and addresses optimized strategies from the recent literature for performing co-expression analysis on plant biological systems. Additionally, we describe the combined interpretation of co-expression analysis with other genomic data to enhance the generation of biologically relevant information.

Effect of Acyl Activating Enzyme (AAE) 3 on the growth and development of Medicago truncatula.

The Acyl-Activating Enzyme (AAE) 3 gene encodes an oxalyl-CoA synthetase that catalyzes the conversion of oxalate to oxalyl-CoA in a CoA and ATP-dependent manner. Although the biochemical activity of AAE3 has been established, its biological role in plant growth and development remains unclear. To advance our understanding of the role of AAE3 in plant growth and development, we report here the characterization of two Medicago truncatula AAE3 (Mtaae3) mutants. Characterization of a Mtaae3 RNAi mutant revealed an accumulation of calcium oxalate crystals and increased seed permeability. These phenotypes were also exhibited in the Arabidopsis aae3 (Ataae3) mutants. Unlike the Ataae3 mutants, the Mtaae3 RNAi mutant did not show a reduction in vegetative growth, decreased seed germination, or increased seed calcium concentration. In an effort to clarify these phenotypic differences, a Mtaae3 Tnt1 mutant was identified and characterized. This Mtaae3 Tnt1 mutant displayed reduced vegetative growth, decreased seed germination, and increased seed calcium concentration as well as an accumulation of calcium oxalate crystals and increased seed permeability as found in Ataae3. Overall, the results presented here show the importance of AAE3 in the growth and development of plants. In addition, this study highlights the ability to separate specific growth and development phenotypes based on the level of AAE3 gene expression.

Elicitors and defense gene induction in plants with altered lignin compositions.

A reduction in the lignin content in transgenic plants induces the ectopic expression of defense genes, but the importance of altered lignin composition in such phenomena remains unclear. Two Arabidopsis lines with similar lignin contents, but strikingly different lignin compositions, exhibited different quantitative and qualitative transcriptional responses. Plants with lignin composed primarily of guaiacyl units overexpressed genes responsive to oomycete and bacterial pathogen attack, whereas plants with lignin composed primarily of syringyl units expressed a far greater number of defense genes, including some associated with cis-jasmone-mediated responses to aphids; these plants exhibited altered responsiveness to bacterial and aphid inoculation. Several of the defense genes were differentially induced by water-soluble extracts from cell walls of plants of the two lines. Glycome profiling, fractionation and enzymatic digestion studies indicated that the different lignin compositions led to differential extractability of a range of heterogeneous oligosaccharide epitopes, with elicitor activity originating from different cell wall polymers. Alteration of lignin composition affects interactions with plant cell wall matrix polysaccharides to alter the sequestration of multiple latent defense signal molecules with an impact on biotic stress responses.

Comparative cell-specific transcriptomics reveals differentiation of C4 photosynthesis pathways in switchgrass and other C4 lineages.

Almost all C4 plants require the co-ordination of the adjacent and fully differentiated cell types, mesophyll (M) and bundle sheath (BS). The C4 photosynthetic pathway operates through two distinct subtypes based on how malate is decarboxylated in BS cells; through NAD-malic enzyme (NAD-ME) or NADP-malic enzyme (NADP-ME). The diverse or unique cell-specific molecular features of M and BS cells from separate C4 subtypes of independent lineages remain to be determined. We here provide an M/BS cell type-specific transcriptome data set from the monocot NAD-ME subtype switchgrass (Panicum virgatum). A comparative transcriptomics approach was then applied to compare the M/BS mRNA profiles of switchgrass, monocot NADP-ME subtype C4 plants maize and Setaria viridis, and dicot NAD-ME subtype Cleome gynandra. We evaluated the convergence in the transcript abundance of core components in C4 photosynthesis and transcription factors to establish Kranz anatomy, as well as gene distribution of biological functions, in these four independent C4 lineages. We also estimated the divergence between NAD-ME and NADP-ME subtypes of C4 photosynthesis in the two cell types within C4 species, including differences in genes encoding decarboxylating enzymes, aminotransferases, and metabolite transporters, and differences in the cell-specific functional enrichment of RNA regulation and protein biogenesis/homeostasis. We suggest that C4 plants of independent lineages in both monocots and dicots underwent convergent evolution to establish C4 photosynthesis, while distinct C4 subtypes also underwent divergent processes for the optimization of M and BS cell co-ordination. The comprehensive data sets in our study provide a basis for further research on evolution of C4 species.

A deep transcriptomic analysis of pod development in the vanilla orchid (Vanilla planifolia).

BACKGROUND: Pods of the vanilla orchid (Vanilla planifolia) accumulate large amounts of the flavor compound vanillin (3-methoxy, 4-hydroxy-benzaldehyde) as a glucoside during the later stages of their development. At earlier stages, the developing seeds within the pod synthesize a novel lignin polymer, catechyl (C) lignin, in their coats. Genomic resources for determining the biosynthetic routes to these compounds and other flavor components in V. planifolia are currently limited. RESULTS: Using next-generation sequencing technologies, we have generated very large gene sequence datasets from vanilla pods at different times of development, and representing different tissue types, including the seeds, hairs, placental and mesocarp tissues. This developmental series was chosen as being the most informative for interrogation of pathways of vanillin and C-lignin biosynthesis in the pod and seed, respectively. The combined 454/Illumina RNA-seq platforms provide both deep sequence coverage and high quality de novo transcriptome assembly for this non-model crop species. CONCLUSIONS: The annotated sequence data provide a foundation for understanding multiple aspects of the biochemistry and development of the vanilla bean, as exemplified by the identification of candidate genes involved in lignin biosynthesis. Our transcriptome data indicate that C-lignin formation in the seed coat involves coordinate expression of monolignol biosynthetic genes with the exception of those encoding the caffeoyl coenzyme A 3-O-methyltransferase for conversion of caffeoyl to feruloyl moieties. This database provides a general resource for further studies on this important flavor species.

Plant Transcriptome Analysis with HISAT-StringTie-Ballgown and TopHat-Cufflinks Pipelines.

The development of next-generation sequencing technology has led to a burst of data in a single assay. Management of a large dataset requires high demands on bioinformatic skills and computing resources. Here we present two popular pipelines for RNA-seq data analysis, using open-source software tools HISAT-StringTie-Ballgown and TopHat-Cufflinks. To meet the need of plant scientist, we describe in detail how to perform such comprehensive analysis beginning with raw RNA-seq reads and available reference genome. It allows biologists to align short reads to a reference genome, measure the transcript abundance, and analyze gene differential expression under two or more conditions. We also discuss other RNA-seq tools that are comparable or alternative to this protocol.

DNA Methylation Dynamics in Response to Drought Stress in Crops.

Drought is one of the most hazardous environmental factors due to its severe damage on plant growth, development and productivity. Plants have evolved complex regulatory networks and resistance strategies for adaptation to drought stress. As a conserved epigenetic regulation, DNA methylation dynamically alters gene expression and chromosome interactions in plants' response to abiotic stresses. The development of omics technologies on genomics, epigenomics and transcriptomics has led to a rapid increase in research on epigenetic variation in non-model crop species. In this review, we summarize the most recent findings on the roles of DNA methylation under drought stress in crops, including methylating and demethylating enzymes, the global methylation dynamics, the dual regulation of DNA methylation on gene expression, the RNA-dependent DNA methylation (RdDM) pathway, alternative splicing (AS) events and long non-coding RNAs (lnc RNAs). We also discuss drought-induced stress memory. These epigenomic findings provide valuable potential for developing strategies to improve crop drought tolerance.