RNA-Seq analysis identifies genes associated with differential reproductive success under drought-stress in accessions of wild barley Hordeum spontaneum

Pangenome characterization and analysis of the NAC gene family reveals genes for Sclerotinia sclerotiorum resistance in sunflower (Helianthus annuus)

BACKGROUND

Sunflower (Helianthus annuus) is one of the most important economic crops in oilseed production worldwide. The different cultivars exhibit variability in their resistance genes. The NAC transcription factor (TF) family plays diverse roles in plant development and stress responses. With the completion of the H. annuus genome sequence, the entire complement of genes coding for NACs has been identified. However, the reference genome of a single individual cannot cover all the genetic information of the species.

RESULTS

Considering only a single reference genome to study gene families will miss many meaningful genes. A pangenome-wide survey and characterization of the NAC genes in sunflower species were conducted. In total, 139 HaNAC genes are identified, of which 114 are core and 25 are variable. Phylogenetic analysis of sunflower NAC proteins categorizes these proteins into 16 subgroups. 138 HaNACs are randomly distributed on 17 chromosomes. SNP-based haplotype analysis shows haplotype diversity of the HaNAC genes in wild accessions is richer than in landraces and modern cultivars. Ten HaNAC genes in the basal stalk rot (BSR) resistance quantitative trait loci (QTL) are found. A total of 26 HaNAC genes are differentially expressed in response to Sclerotinia head rot (SHR). A total of 137 HaNAC genes are annotated in Gene Ontology (GO) and are classified into 24 functional groups. GO functional enrichment analysis reveals that HaNAC genes are involved in various functions of the biological process.

CONCLUSIONS

We identified NAC genes in H. annuus (HaNAC) on a pangenome-wide scale and analyzed S. sclerotiorum resistance-related NACs. This study provided a theoretical basis for further genomic improvement targeting resistance-related NAC genes in sunflowers.

Heat and drought induced transcriptomic changes in barley varieties with contrasting stress response phenotypes

Drought and heat stress substantially impact plant growth and productivity. When subjected to drought or heat stress, plants exhibit reduction in growth resulting in yield losses. The occurrence of these two stresses together intensifies their negative effects. Unraveling the molecular changes in response to combined abiotic stress is essential to breed climate-resilient crops. In this study, transcriptome profiles were compared between stress-tolerant (Otis), and stress-sensitive (Golden Promise) barley genotypes subjected to drought, heat, and combined heat and drought stress for five days during heading stage. The major differences that emerged from the transcriptome analysis were the overall number of differentially expressed genes was relatively higher in Golden Promise (GP) compared to Otis. The differential expression of more than 900 transcription factors in GP and Otis may aid this transcriptional reprogramming in response to abiotic stress. Secondly, combined heat and water deficit stress results in a unique and massive transcriptomic response that cannot be predicted from individual stress responses. Enrichment analyses of gene ontology terms revealed unique and stress type-specific adjustments of gene expression. Weighted Gene Co-expression Network Analysis identified genes associated with RNA metabolism and Hsp70 chaperone components as hub genes that can be useful for engineering tolerance to multiple abiotic stresses. Comparison of the transcriptomes of unstressed Otis and GP plants identified several genes associated with biosynthesis of antioxidants and osmolytes were higher in the former that maybe providing innate tolerance capabilities to effectively combat hostile conditions. Lines with different repertoire of innate tolerance mechanisms can be effectively leveraged in breeding programs for developing climate-resilient barley varieties with superior end-use traits.

Combined analysis of transcriptome and metabolome reveals the molecular mechanism and candidate genes of Haloxylon drought tolerance

Haloxylon ammodendron and Haloxylon persicum, as typical desert plants, show strong drought tolerance and environmental adaptability. They are ideal model plants for studying the molecular mechanisms of drought tolerance. Transcriptomic and metabolomic analyses were performed to reveal the response mechanisms of H. ammodendron and H. persicum to a drought environment at the levels of transcription and physiological metabolism. The results showed that the morphological structures of H. ammodendron and H. persicum showed adaptability to drought stress. Under drought conditions, the peroxidase activity, abscisic acid content, auxin content, and gibberellin content of H. ammodendron increased, while the contents of proline and malondialdehyde decreased. The amino acid content of H. persicum was increased, while the contents of proline, malondialdehyde, auxin, and gibberellin were decreased. Under drought conditions, 12,233 and 17,953 differentially expressed genes (DEGs) were identified in H. ammodendron and H. persicum , respectively, including members of multiple transcription factor families such as FAR1, AP2/ERF, C2H2, bHLH, MYB, C2C2, and WRKY that were significantly up-regulated under drought stress. In the positive ion mode, 296 and 452 differential metabolites (DEMs) were identified in H. ammodendron and H. persicum, respectively; in the negative ion mode, 252 and 354 DEMs were identified, primarily in carbohydrate and lipid metabolism. A combined transcriptome and metabolome analysis showed that drought stress promoted the glycolysis/gluconeogenesis pathways of H. ammodendron and H. persicum and increased the expression of amino acid synthesis pathways, consistent with the physiological results. In addition, transcriptome and metabolome were jointly used to analyze the expression changes of the genes/metabolites of H. ammodendron and H. persicum that were associated with drought tolerance but were regulated differently in the two plants. This study identified drought-tolerance genes and metabolites in H. ammodendron and H. persicum and has provided new ideas for studying the drought stress response of Haloxylon.

Effects of Salt Stress on Transcriptional and Physiological Responses in Barley Leaves with Contrasting Salt Tolerance

Salt stress negatively impacts crop production worldwide. Genetic diversity among barley (Hordeum vulgare) landraces adapted to adverse conditions should provide a valuable reservoir of tolerance genes for breeding programs. To identify molecular and biochemical differences between barley genotypes, transcriptomic and antioxidant enzyme profiles along with several morpho-physiological features were compared between salt-tolerant (Boulifa) and salt-sensitive (Testour) genotypes subjected to salt stress. Decreases in biomass, photosynthetic parameters, and relative water content were low in Boulifa compared to Testour. Boulifa had better antioxidant protection against salt stress than Testour, with greater antioxidant enzymes activities including catalase, superoxide dismutase, and guaiacol peroxidase. Transcriptome assembly for both genotypes revealed greater accumulation of differentially expressed transcripts in Testour compared to Boulifa, emphasizing the elevated transcriptional response in Testour following salt exposure. Various salt-responsive genes, including the antioxidant catalase 3, the osmoprotectant betaine aldehyde dehydrogenase 2, and the transcription factors MYB20 and MYB41, were induced only in Boulifa. By contrast, several genes associated with photosystems I and II, and light receptor chlorophylls A and B, were more repressed in Testour. Co-expression network analysis identified specific gene modules correlating with differences in genotypes and morpho-physiological traits. Overall, salinity-induced differential transcript accumulation underlies the differential morpho-physiological response in both genotypes and could be important for breeding salt tolerance in barley.

Are we there yet? Driving the road to evolutionary graph-pangenomics

With increase in the number of sequenced genomes, it is now recognized that graph-based pangenomes can provide a comprehensive platform to study diversity in a population or species, from point mutations to large chromosomal rearrangements. By incorporating concepts from graph theory, a graph-pangenome can be studied directly to identify genomic regions and genes that underlie important evolutionary processes and traits. Here, I discuss how basic concepts in graph theory can be implemented to address questions in evolutionary genomics and guide future breeding efforts. Despite its compelling versatility, a graph-pangenome assembly is still challenging especially in species with large complex genomes. As technology is rapidly improving, the graph-pangenome is expected to become a central platform in genomics studies and applications. Thus, development of tools and methods that exploit the graph structure are urged to pave the route to evolutionary graph-pangenomics.

Physical geography, isolation by distance and environmental variables shape genomic variation of wild barley (Hordeum vulgare L. ssp. spontaneum) in the Southern Levant

Determining the extent of genetic variation that reflects local adaptation in crop-wild relatives is of interest for the purpose of identifying useful genetic diversity for plant breeding. We investigated the association of genomic variation with geographical and environmental factors in wild barley (Hordeum vulgare L. ssp. spontaneum) populations of the Southern Levant using genotyping by sequencing (GBS) of 244 accessions in the Barley 1K+ collection. The inference of population structure resulted in four genetic clusters that corresponded to eco-geographical habitats and a significant association between lower gene flow rates and geographical barriers, e.g. the Judaean Mountains and the Sea of Galilee. Redundancy analysis (RDA) revealed that spatial autocorrelation explained 45% and environmental variables explained 15% of total genomic variation. Only 4.5% of genomic variation was solely attributed to environmental variation if the component confounded with spatial autocorrelation was excluded. A synthetic environmental variable combining latitude, solar radiation, and accumulated precipitation explained the highest proportion of genomic variation (3.9%). When conditioned on population structure, soil water capacity was the most important environmental variable explaining 1.18% of genomic variation. Genome scans with outlier analysis and genome-environment association studies were conducted to identify adaptation signatures. RDA and outlier methods jointly detected selection signatures in the pericentromeric regions, which have reduced recombination, of the chromosomes 3H, 4H, and 5H. However, selection signatures mostly disappeared after correction for population structure. In conclusion, adaptation to the highly diverse environments of the Southern Levant over short geographical ranges had a limited effect on the genomic diversity of wild barley. This highlighted the importance of nonselective forces in genetic differentiation.

Progressive Genomic Approaches to Explore Drought- and Salt-Induced Oxidative Stress Responses in Plants under Changing Climate

Drought and salinity are the major environmental abiotic stresses that negatively impact crop development and yield. To improve yields under abiotic stress conditions, drought- and salinity-tolerant crops are key to support world crop production and mitigate the demand of the growing world population. Nevertheless, plant responses to abiotic stresses are highly complex and controlled by networks of genetic and ecological factors that are the main targets of crop breeding programs. Several genomics strategies are employed to improve crop productivity under abiotic stress conditions, but traditional techniques are not sufficient to prevent stress-related losses in productivity. Within the last decade, modern genomics studies have advanced our capabilities of improving crop genetics, especially those traits relevant to abiotic stress management. This review provided updated and comprehensive knowledge concerning all possible combinations of advanced genomics tools and the gene regulatory network of reactive oxygen species homeostasis for the appropriate planning of future breeding programs, which will assist sustainable crop production under salinity and drought conditions.

Deploying Big Data to Crack the Genotype to Phenotype Code

Mechanistically connecting genotypes to phenotypes is a longstanding and central mission of biology. Deciphering these connections will unite questions and datasets across all scales from molecules to ecosystems. Although high-throughput sequencing has provided a rich platform on which to launch this effort, tools for deciphering mechanisms further along the genome to phenome pipeline remain limited. Machine learning approaches and other emerging computational tools hold the promise of augmenting human efforts to overcome these obstacles. This vision paper is the result of a Reintegrating Biology Workshop, bringing together the perspectives of integrative and comparative biologists to survey challenges and opportunities in cracking the genotype to phenotype code and thereby generating predictive frameworks across biological scales. Key recommendations include promoting the development of minimum "best practices" for the experimental design and collection of data; fostering sustained and long-term data repositories; promoting programs that recruit, train, and retain a diversity of talent; and providing funding to effectively support these highly cross-disciplinary efforts. We follow this discussion by highlighting a few specific transformative research opportunities that will be advanced by these efforts.

EORNA, a barley gene and transcript abundance database

A high-quality, barley gene reference transcript dataset (BaRTv1.0), was used to quantify gene and transcript abundances from 22 RNA-seq experiments, covering 843 separate samples. Using the abundance data we developed a Barley Expression Database (EORNA*) to underpin a visualisation tool that displays comparative gene and transcript abundance data on demand as transcripts per million (TPM) across all samples and all the genes. EORNA provides gene and transcript models for all of the transcripts contained in BaRTV1.0, and these can be conveniently identified through either BaRT or HORVU gene names, or by direct BLAST of query sequences. Browsing the quantification data reveals cultivar, tissue and condition specific gene expression and shows changes in the proportions of individual transcripts that have arisen via alternative splicing. TPM values can be easily extracted to allow users to determine the statistical significance of observed transcript abundance variation among samples or perform meta analyses on multiple RNA-seq experiments. * Eòrna is the Scottish Gaelic word for Barley.

EORNA, a barley gene and transcript abundance database

A high-quality, barley gene reference transcript dataset (BaRTv1.0), was used to quantify gene and transcript abundances from 22 RNA-seq experiments, covering 843 separate samples. Using the abundance data we developed a Barley Expression Database (EORNA*) to underpin a visualisation tool that displays comparative gene and transcript abundance data on demand as transcripts per million (TPM) across all samples and all the genes. EORNA provides gene and transcript models for all of the transcripts contained in BaRTV1.0, and these can be conveniently identified through either BaRT or HORVU gene names, or by direct BLAST of query sequences. Browsing the quantification data reveals cultivar, tissue and condition specific gene expression and shows changes in the proportions of individual transcripts that have arisen via alternative splicing. TPM values can be easily extracted to allow users to determine the statistical significance of observed transcript abundance variation among samples or perform meta analyses on multiple RNA-seq experiments. * Eòrna is the Scottish Gaelic word for Barley.

New insights on key genes involved in drought stress response of barley: gene networks reconstruction, hub, and promoter analysis

Background

Barley (Hordeum vulgare L.) is one of the most important cereals worldwide. Although this crop is drought-tolerant, water deficiency negatively affects its growth and production. To detect key genes involved in drought tolerance in barley, a reconstruction of the related gene network and discovery of the hub genes would help. Here, drought-responsive genes in barley were collected through analysis of the available microarray datasets (− 5 ≥ Fold change ≥ 5, adjusted p value ≤ 0.05). Protein-protein interaction (PPI) networks were reconstructed.

Results

The hub genes were identified by Cytoscape software using three Cyto-hubba algorithms (Degree, Closeness, and MNC), leading to the identification of 17 and 16 non-redundant genes at vegetative and reproductive stages, respectively. These genes consist of some transcription factors such as HvVp1, HvERF4, HvFUS3, HvCBF6, DRF1.3, HvNAC6, HvCO5, and HvWRKY42, which belong to AP2, NAC, Zinc-finger, and WRKY families. In addition, the expression pattern of four hub genes was compared between the two studied cultivars, i.e., “Yousef” (drought-tolerant) and “Morocco” (susceptible). The results of real-time PCR revealed that the expression patterns corresponded well with those determined by the microarray. Also, promoter analysis revealed that some TF families, including AP2, NAC, Trihelix, MYB, and one modular (composed of two HD-ZIP TFs), had a binding site in 85% of promoters of the drought-responsive genes and of the hub genes in barley.

Conclusions

The identified hub genes, especially those from AP2 and NAC families, might be among key TFs that regulate drought-stress response in barley and are suggested as promising candidate genes for further functional analysis.

Tapping Diversity From the Wild: From Sampling to Implementation

The diversity observed among crop wild relatives (CWRs) and their ability to flourish in unfavorable and harsh environments have drawn the attention of plant scientists and breeders for many decades. However, it is also recognized that the benefit gained from using CWRs in breeding is a potential rose between thorns of detrimental genetic variation that is linked to the trait of interest. Despite the increased interest in CWRs, little attention was given so far to the statistical, analytical, and technical considerations that should guide the sampling design, the germplasm characterization, and later its implementation in breeding. Here, we review the entire process of sampling and identifying beneficial genetic variation in CWRs and the challenge of using it in breeding. The ability to detect beneficial genetic variation in CWRs is strongly affected by the sampling design which should be adjusted to the spatial and temporal variation of the target species, the trait of interest, and the analytical approach used. Moreover, linkage disequilibrium is a key factor that constrains the resolution of searching for beneficial alleles along the genome, and later, the ability to deplete linked deleterious genetic variation as a consequence of genetic drag. We also discuss how technological advances in genomics, phenomics, biotechnology, and data science can improve the ability to identify beneficial genetic variation in CWRs and to exploit it in strive for higher-yielding and sustainable crops.

Gene co-expression network analysis to identify critical modules and candidate genes of drought-resistance in wheat

AIM

To establish a gene co-expression network for identifying principal modules and hub genes that are associated with drought resistance mechanisms, analyzing their mechanisms, and exploring candidate genes.

METHODS AND FINDINGS

42 data sets including PRJNA380841 and PRJNA369686 were used to construct the co-expression network through weighted gene co-expression network analysis (WGCNA). A total of 1,896,897,901 (284.30 Gb) clean reads and 35,021 differentially expressed genes (DEGs) were obtained from 42 samples. Functional enrichment analysis indicated that photosynthesis, DNA replication, glycolysis/gluconeogenesis, starch and sucrose metabolism, arginine and proline metabolism, and cell cycle were significantly influenced by drought stress. Furthermore, the DEGs with similar expression patterns, detected by K-means clustering, were grouped into 29 clusters. Genes involved in the modules, such as dark turquoise, yellow, and brown, were found to be appreciably linked with drought resistance. Twelve central, greatly correlated genes in stage-specific modules were subsequently confirmed and validated at the transcription levels, including TraesCS7D01G417600.1 (PP2C), TraesCS5B01G565300.1 (ERF), TraesCS4A01G068200.1 (HSP), TraesCS2D01G033200.1 (HSP90), TraesCS6B01G425300.1 (RBD), TraesCS7A01G499200.1 (P450), TraesCS4A01G118400.1 (MYB), TraesCS2B01G415500.1 (STK), TraesCS1A01G129300.1 (MYB), TraesCS2D01G326900.1 (ALDH), TraesCS3D01G227400.1 (WRKY), and TraesCS3B01G144800.1 (GT).

CONCLUSIONS

Analyzing the response of wheat to drought stress during different growth stages, we have detected three modules and 12 hub genes that are associated with drought resistance mechanisms, and five of those genes are newly identified for drought resistance. The references provided by these modules will promote the understanding of the drought-resistance mechanism. In addition, the candidate genes can be used as a basis of transgenic or molecular marker-assisted selection for improving the drought resistance and increasing the yields of wheat.

The Effect of Drought on Transcriptome and Hormonal Profiles in Barley Genotypes With Contrasting Drought Tolerance

Like many cereal crops, barley is also negatively affected by drought stress. However, due to its simple genome as well as enhanced stress resilient nature compared to rice and wheat, barley has been considered as a model to decipher drought tolerance in cereals. In the present study, transcriptomic and hormonal profiles along with several biochemical features were compared between drought-tolerant (Otis) and drought-sensitive (Baronesse) barley genotypes subjected to drought to identify molecular and biochemical differences between the genotypes. The drought-induced decrease in the leaf relative water content, net photosynthesis, and biomass accumulation was relatively low in Otis compared to Baronesse. The hormonal profiles did not reveal significant differences for majority of the compounds other than the GA20 and the cis-zeatin-o-glucoside (c-ZOG), whose levels were greatly increased in Otis compared to Baronesse under drought. The major differences that emerged from the transcriptome analysis are; (1), the overall number of differentially expressed genes was relatively low in drought-tolerant Otis compared to drought-sensitive Baronesse; (2), a wax biosynthesis gene (CER1), and NAC transcription factors were specifically induced in Otis but not in Baronesse; (3), the degree of upregulation of betaine aldehyde dehydrogenase and a homeobox transcription factor (genes with proven roles in imparting drought tolerance), was greater in Otis compared to Baronesse; (4) the extent of downregulation of gene expression profiles for proteins of the reaction center photosystem II (PSII) (D1 and D2) was low in Otis compared to Baronesse; and, (5), alternative splicing (AS) was also found to differ between the genotypes under drought. Taken together, the overall transcriptional responses were low in drought-tolerant Otis but the genes that could confer drought tolerance were either specifically induced or greatly upregulated in the tolerant genotype and these differences could be important for drought tolerance in barley.

Seedling response to water stress in valley oak (Quercus lobata) is shaped by different gene networks across populations

Drought is a major stress for plants, creating a strong selection pressure for traits that enable plant growth and survival in dry environments. Many drought responses are conserved species-wide responses, while others vary among populations distributed across heterogeneous environments. We tested how six populations of the widely distributed California valley oak (Quercus lobata) sampled from contrasting climates would differ in their response to soil drying relative to well-watered controls in a common environment by measuring ecophysiological traits in 93 individuals and gene expression (RNA-seq) in 42 individuals. Populations did not differ in their adjustment of turgor loss point during soil drying, suggesting a generalized species-wide response. Differential expression analysis identified 689 genes with a common response to treatment across populations and 470 genes with population-specific responses. Weighted gene co-expression network analysis (WGCNA) identified groups of genes with similar expression patterns that may be regulated together (gene modules). Several gene modules responded differently to water stress among populations, suggesting regional differences in gene network regulation. Populations from sites with a high mean annual temperature responded to the imposed water stress with significantly greater changes in gene module expression, indicating that these populations may be locally adapted to respond to drought. We propose that this variation among valley oak populations provides a mechanism for differential tolerance to the increasingly frequent and severe droughts in California.

Transcriptome profiling of Gerbera hybrida reveals that stem bending is caused by water stress and regulation of abscisic acid

Background

Gerbera hybrida is one of the most popular cut flowers in the world; however, stem bending, which always happens when gerbera flower harvested from the field, greatly limits its vase life. To date the molecular mechanisms underlying stem bending remain poorly understood.

Results

In this study, we performed high-throughput transcriptome sequencing of gerbera during stem bending using the Illumina sequencing technology. Three cDNA libraries constructed from mRNAs of gerbera stem at stem bending stage 0, 2 and 4 were sequenced. More than 300 million high-quality reads were generated and assembled into 96,492 unigenes. Among them, 34,166 unigenes were functionally annotated based on similarity search with known protein. Sequences derived from plants at different stem bending stages were mapped to the assembled transcriptome, and 9,406 differentially expressed genes (DEGs) were identified. Through Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, specific pathways were identified during the stem bending process, such as phenylpropanoid biosynthesis pathway, phenylalanine metabolism pathway, starch and sucrose metabolism pathway, and plant hormone signal transduction pathway. A total of 211 transcription factors (TFs), including TF families involved in plant senescence, such as NAC, MYB, WRKY, and AP2/ERF members, as well as TFs related to water stress tolerance, were shown to be regulated during stem bending. Gene Onotology (GO) functional enrichment analysis indicated that key genes involved in responses to osmotic and oxidative stresses were also varied in expression during this process. Furthermore, analysis of DEGs involved in the hormone signaling pathways and determination of endogenous abscisic acid (ABA) content showed that stem bending may be an ethylene-independent process, but regulated by ABA. In short, our findings suggested that the stem bending of cut gerbera may be caused by the involvement of water stress and regulation of ABA during the postharvest life.

Conclusions

The transcriptome sequences provide a valuable resource in revealing the molecular mechanism underlying stem bending of cut flower and offer novel genes that can be used to guide future studies for ornamental plant breeding.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5961-1) contains supplementary material, which is available to authorized users.

Identification of Factors Linked to Higher Water-Deficit Stress Tolerance in Amaranthus hypochondriacus Compared to Other Grain Amaranths and A. hybridus, Their Shared Ancestor

Water deficit stress (WDS)-tolerance in grain amaranths (Amaranthus hypochondriacus, A. cruentus and A. caudatus), and A. hybridus, their presumed shared ancestor, was examined. A. hypochondriacus was the most WDS-tolerant species, a trait that correlated with an enhanced osmotic adjustment (OA), a stronger expression of abscisic acid (ABA) marker genes and a more robust sugar starvation response (SSR). Superior OA was supported by higher basal hexose (Hex) levels and high Hex/sucrose (Suc) ratios in A. hypochondriacus roots, which were further increased during WDS. This coincided with increased invertase, amylase and sucrose synthase activities and a strong depletion of the starch reserves in leaves and roots. The OA was complemented by the higher accumulation of proline, raffinose, and other probable raffinose-family oligosaccharides of unknown structure in leaves and/or roots. The latter coincided with a stronger expression of Galactinol synthase 1 and Raffinose synthase in leaves. Increased SnRK1 activity and expression levels of the class II AhTPS9 and AhTPS11 trehalose phosphate synthase genes, recognized as part of the SSR network in Arabidopsis, were induced in roots of stressed A. hypochondriacus. It is concluded that these physiological modifications improved WDS in A. hypochondriacus by raising its water use efficiency.

A pan-transcriptome analysis shows that disease resistance genes have undergone more selection pressure during barley domestication

Background

It has become clear in recent years that many genes in a given species may not be found in a single genotype thus using sequences from a single genotype as reference may not be adequate for various applications.

Results

In this study we constructed a pan-transcriptome for barley by de novo assembling 288 sets of RNA-seq data from 32 cultivated barley genotypes and 31 wild barley genotypes. The pan-transcriptome consists of 756,632 transcripts with an average N50 length of 1240 bp. Of these, 289,697 (38.2%) were not found in the genome of the international reference genotype Morex. The novel transcripts are enriched with genes associated with responses to different stresses and stimuli. At the pan-transcriptome level, genotypes of wild barley have a higher proportion of disease resistance genes than cultivated ones.

Conclusions

We demonstrate that the use of the pan-transcriptome dramatically improved the efficiency in detecting variation in barley. Analysing the pan-transcriptome also found that, compared with those in other categories, disease resistance genes have gone through stronger selective pressures during domestication.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5357-7) contains supplementary material, which is available to authorized users.

Comparative transcriptomic analysis reveals common molecular factors responsive to heat and drought stress in Agrostis stolonifera

Heat and drought stress are primary abiotic stresses confining growth of cool-season grass species during summer. The objective of this study was to identify common molecular factors and metabolic pathways associated with heat and drought responses in creeping bentgrass (Agrostis stolonifera) by comparative analysis of transcriptomic profiles between plants exposed to heat and drought stress. Plants were exposed to heat stress (35/30 °C day/night temperature) or drought stress by withholding irrigation for 21 d in growth chambers. Transcriptomic profiling by RNA-seq in A. stolonifera (cv. 'Penncross') found 670 commonly up-regulated and 812 commonly down-regulated genes by heat and drought stress. Transcriptional up-regulations of differentially expressed genes (DEGs) due to heat and drought stress include genes that were highly enriched in oxylipin biosynthetic process and proline biosynthetic process. Transcriptional down-regulations of genes under heat and drought stress were highly enriched and involved in thiamine metabolic process and calcium sensing receptor. These commonly-regulated genes by heat and drought stress identified in A. stolonifera suggested that drought and heat responses shared such common molecular factors and pathways, which could be potential candidate genes for genetic modification of improving plant tolerance to the combined heat and drought stress.

Optimizing RNA-seq studies to investigate herbicide resistance

Transcriptomic profiling, specifically via RNA sequencing (RNA-seq), is becoming one of the more commonly used methods for investigating non-target site resistance (NTSR) to herbicides due to its high throughput capabilities and utility in organisms with little to no previous sequence information. A review of the weed science RNA-seq literature revealed some basic principles behind generating quality data from these types of studies. First, studies that included more replicates per biotype and took steps to control for genetic background had significantly better control of false positives and, consequently, shorter lists of potential resistance genes to sift through. Pooling of biological replicates prior to sequencing was successful in some cases, but likely contributed to an overall increase in the false discovery rate. Although the inclusion of herbicide-treated samples was common across most of the studies, it ultimately introduced difficulties in interpretation of the final results due to challenges in capturing the right sampling window after treatment and to the induction of stress responses in the injured herbicide-sensitive plants. RNA-seq is an effective tool for NTSR gene discovery, but careful consideration should be given to finding the most powerful and cost-effective balance between replicate number, sequencing depth and treatment number. © 2017 Society of Chemical Industry.

Using Biotechnology-Led Approaches to Uplift Cereal and Food Legume Yields in Dryland Environments

Drought and heat in dryland agriculture challenge the enhancement of crop productivity and threaten global food security. This review is centered on harnessing genetic variation through biotechnology-led approaches to select for increased productivity and stress tolerance that will enhance crop adaptation in dryland environments. Peer-reviewed literature, mostly from the last decade and involving experiments with at least two seasons' data, form the basis of this review. It begins by highlighting the adverse impact of the increasing intensity and duration of drought and heat stress due to global warming on crop productivity and its impact on food and nutritional security in dryland environments. This is followed by (1) an overview of the physiological and molecular basis of plant adaptation to elevated CO2 (eCO2), drought, and heat stress; (2) the critical role of high-throughput phenotyping platforms to study phenomes and genomes to increase breeding efficiency; (3) opportunities to enhance stress tolerance and productivity in food crops (cereals and grain legumes) by deploying biotechnology-led approaches [pyramiding quantitative trait loci (QTL), genomic selection, marker-assisted recurrent selection, epigenetic variation, genome editing, and transgene) and inducing flowering independent of environmental clues to match the length of growing season; (4) opportunities to increase productivity in C3 crops by harnessing novel variations (genes and network) in crops' (C3, C4) germplasm pools associated with increased photosynthesis; and (5) the adoption, impact, risk assessment, and enabling policy environments to scale up the adoption of seed-technology to enhance food and nutritional security. This synthesis of technological innovations and insights in seed-based technology offers crop genetic enhancers further opportunities to increase crop productivity in dryland environments.

RNA-Seq analysis revealed genes associated with drought stress response in kabuli chickpea (Cicer arietinum L.)

Drought is the most important constraint that effects chickpea production globally. RNA-Seq has great potential to dissect the molecular mechanisms of tolerance to environmental stresses. Transcriptome profiles in roots and shoots of two contrasting Iranian kabuli chickpea genotypes (Bivanij and Hashem) were investigated under water-limited conditions at early flowering stage using RNA-Seq approach. A total of 4,572 differentially expressed genes (DEGs) were identified. Of these, 261 and 169 drought stress responsive genes were identified in the shoots and the roots, respectively, and 17 genes were common in the shoots and the roots. Gene Ontology (GO) analysis revealed several sub-categories related to the stress, including response to stress, defense response and response to stimulus in the tolerant genotype Bivanij as compared to the sensitive genotype Hashem under drought stress. In addition, several Transcription factors (TFs) were identified in major metabolic pathways such as, ABA, proline and flavonoid biosynthesis. Furthermore, a number of the DEGs were observed in "QTL-hotspot" regions which were reported earlier in chickpea. Drought tolerance dissection in the genotypes revealed that the genes and the pathways involved in shoots of Bivanij were the most important factor to make a difference between the genotypes for drought tolerance. The identified TFs in the experiment, particularly those which were up-regulated in shoots of Bivanij during drought stress, were potential candidates for enhancing tolerance to drought.

Genome-wide expression profiling of leaves and roots of watermelon in response to low nitrogen

BACKGROUND

Nitrogen (N) is a key macronutrient required for plant growth and development. In this study, watermelon plants were grown under hydroponic conditions at 0.2 mM N, 4.5 mM N, and 9 mM N for 14 days.

RESULTS

Dry weight and photosynthetic assimilation at low N (0.2 mM) was reduced by 29 and 74% compared with high N (9 mM). The photochemical activity (Fv/Fm) was also reduced from 0.78 at high N to 0.71 at low N. The N concentration in the leaf, stem, and root of watermelon under low N conditions was reduced by 68, 104, and 108%, respectively compared with 9 mM N treatment after 14 days of N treatment. In the leaf tissues of watermelon grown under low N conditions, 9598 genes were differentially expressed, out of which 4533 genes (47.22%) were up-regulated whereas, 5065 genes (52.78%) were down-regulated compared with high N. Similarly in the root tissues, 3956 genes were differentially expressed, out of which 1605 genes were up-regulated (40.57%) and 2351 genes were down-regulated (59.43%), compared with high N. Our results suggest that leaf tissues are more sensitive to N deficiency compared with root tissues. The gene ontology (GO) analysis showed that the availability of N significantly affected 19 biological processes, 8 cell component metabolic pathways, and 3 molecular functions in the leaves; and 13 biological processes, 12 molecular functions, and 5 cell component metabolic pathways in the roots of watermelon. The low affinity nitrate transporters, high affinity nitrate transporters, ammonium transporters, genes related with nitrogen assimilation, and chlorophyll and photosynthesis were expressed differentially in response to low N. Three nitrate transporters (Cla010066, Cla009721, Cla012765) substantially responded to low nitrate supply in the root and leaf tissues. Additionally, a large number of transcription factors (1365) were involved in adaptation to low N availability. The major transcription factor families identified in this study includes MYB, AP2-EREBP, bHLH, C2H2 and NAC.

CONCLUSION

Candidate genes identified in this study for nitrate uptake and transport can be targeted and utilized for further studies in watermelon breeding and improvement programs to improve N uptake and utilization efficiency.

Prominent alterations of wild barley leaf transcriptome in response to individual and combined drought acclimation and heat shock conditions

Under field conditions, drought and heat stress typically happen simultaneously and their negative impact on the agricultural production is expected to increase worldwide under the climate change scenario. In this study, we performed RNA-sequencing analysis on leaves of wild barley (Hordeum spontaneum) originated from the northern coastal region of Egypt following individual drought acclimation (DA) and heat shock (HS) treatments and their combination (CS, combined stresses) to distinguish the unique and shared differentially expressed genes (DEG). Results indicated that the number of unique genes that were differentially expressed following HS treatment exceeded the number of those expressed following DA. In addition, the number of genes that were uniquely differentially expressed in response to CS treatment exceeded the number of those of shared responses to individual DA and HS treatments. These results indicate a better adaptation of the Mediterranean wild barley to drought conditions when compared with heat stress. It also manifests that the wild barley response to CS tends to be unique rather than common. Annotation of DEG showed that metabolic processes were the most influenced biological function in response to the applied stresses.

Genomic adaptation to drought in wild barley is driven by edaphic natural selection at the Tabigha Evolution Slope

Ecological divergence at a microsite suggests adaptive evolution, and this study examined two abutting wild barley populations, each 100 m across, differentially adapted to drought tolerance on two contrasting soil types, Terra Rossa and basalt at the Tabigha Evolution Slope, Israel. We resequenced the genomes of seven and six wild barley genotypes inhabiting the Terra Rossa and basalt soils, respectively, and identified a total of 69,192,653 single-nucleotide variants (SNVs) and insertions/deletions in comparison with a reference barley genome. Comparative genomic analysis between these abutting wild barley populations involved 19,615,087 high-quality SNVs. The results revealed dramatically different selection sweep regions relevant to drought tolerance driven by edaphic natural selection within 2,577 selected genes in these regions, including key drought-responsive genes associated with ABA synthesis and degradation (such as Cytochrome P450 protein) and ABA receptor complex (such as PYL2, SNF1-related kinase). The genetic diversity of the wild barley population inhabiting Terra Rossa soil is much higher than that from the basalt soil. Additionally, we identified different sets of genes for drought adaptation in the wild barley populations from Terra Rossa soil and from wild barley populations from Evolution Canyon I at Mount Carmel. These genes are associated with abscisic acid signaling, signaling and metabolism of reactive oxygen species, detoxification and antioxidative systems, rapid osmotic adjustment, and deep root morphology. The unique mechanisms for drought adaptation of the wild barley from the Tabigha Evolution Slope may be useful for crop improvement, particularly for breeding of barley cultivars with high drought tolerance.

No Time to Waste: Transcriptome Study Reveals that Drought Tolerance in Barley May Be Attributed to Stressed-Like Expression Patterns that Exist before the Occurrence of Stress

Plant survival in adverse environmental conditions requires a substantial change in the metabolism, which is reflected by the extensive transcriptome rebuilding upon the occurrence of the stress. Therefore, transcriptomic studies offer an insight into the mechanisms of plant stress responses. Here, we present the results of global gene expression profiling of roots and leaves of two barley genotypes with contrasting ability to cope with drought stress. Our analysis suggests that drought tolerance results from a certain level of transcription of stress-influenced genes that is present even before the onset of drought. Genes that predispose the plant to better drought survival play a role in the regulatory network of gene expression, including several transcription factors, translation regulators and structural components of ribosomes. An important group of genes is involved in signaling mechanisms, with significant contribution of hormone signaling pathways and an interplay between ABA, auxin, ethylene and brassinosteroid homeostasis. Signal transduction in a drought tolerant genotype may be more efficient through the expression of genes required for environmental sensing that are active already during normal water availability and are related to actin filaments and LIM domain proteins, which may function as osmotic biosensors. Better survival of drought may also be attributed to more effective processes of energy generation and more efficient chloroplasts biogenesis. Interestingly, our data suggest that several genes involved in a photosynthesis process are required for the establishment of effective drought response not only in leaves, but also in roots of barley. Thus, we propose a hypothesis that root plastids may turn into the anti-oxidative centers protecting root macromolecules from oxidative damage during drought stress. Specific genes and their potential role in building up a drought-tolerant barley phenotype is extensively discussed with special emphasis on processes that take place in barley roots. When possible, the interconnections between particular factors are emphasized to draw a broader picture of the molecular mechanisms of drought tolerance in barley.

Dehydration induced transcriptomic responses in two Tibetan hulless barley (Hordeum vulgare var. nudum) accessions distinguished by drought tolerance

Background

The harsh environment on the Qinghai-Tibetan Plateau gives Tibetan hulless barley (Hordeum vulgare var. nudum) great ability to resist adversities such as drought, salinity, and low temperature, and makes it a good subject for the analysis of drought tolerance mechanism. To elucidate the specific gene networks and pathways that contribute to its drought tolerance, and for identifying new candidate genes for breeding purposes, we performed a transcriptomic analysis using two accessions of Tibetan hulless barley, namely Z772 (drought-tolerant) and Z013 (drought-sensitive).

Results

There were more up-regulated genes of Z772 than Z013 under both mild (5439-VS-2604) and severe (7203-VS-3359) dehydration treatments. Under mild dehydration stress, the pathways exclusively enriched in drought-tolerance genotype Z772 included Protein processing in endoplasmic reticulum, tricarboxylic acid (TCA) cycle, Wax biosynthesis, and Spliceosome. Under severe dehydration stress, the pathways that were mainly enriched in Z772 included Carbon fixation in photosynthetic organisms, Pyruvate metabolism, Porphyrin and chlorophyll metabolism. The main differentially expressed genes (DEGs) in response to dehydration stress and genes whose expression was different between tolerant and sensitive genotypes were presented in this study, respectively. The candidate genes for drought tolerance were selected based on their expression patterns.

Conclusions

The RNA-Seq data obtained in this study provided an initial overview on global gene expression patterns and networks that related to dehydration shock in Tibetan hulless barley. Furthermore, these data provided pathways and a targeted set of candidate genes that might be essential for deep analyzing the molecular mechanisms of plant tolerance to drought stress.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4152-1) contains supplementary material, which is available to authorized users.

Transcriptomics Analyses Reveal Wheat Responses to Drought Stress during Reproductive Stages under Field Conditions

Drought is a major abiotic stress that limits wheat production worldwide. To ensure food security for the rapidly increasing world population, improving wheat yield under drought stress is urgent and relevant. In this study, an RNA-seq analysis was conducted to study the effect of drought on wheat transcriptome changes during reproductive stages under field conditions. Our results indicated that drought stress during early reproductive periods had a more severe impact on wheat development, gene expression and yield than drought stress during flowering. In total, 115,656 wheat genes were detected, including 309 differentially expressed genes (DEGs) which responded to drought at various developmental stages. These DEGs were involved in many critical processes including floral development, photosynthetic activity and stomatal movement. At early developmental stages, the proteins of drought-responsive DEGs were mainly located in the nucleus, peroxisome, mitochondria, plasma membrane and chloroplast, indicating that these organelles play critical roles in drought tolerance in wheat. Furthermore, the validation of five DEGs confirmed their responsiveness to drought under different genetic backgrounds. Functional verification of DEGs of interest will occur in our subsequent research. Collectively, the results of this study not only advanced our understanding of wheat transcriptome changes under drought stress during early reproductive stages but also provided useful targets to manipulate drought tolerance in wheat at different development stages.

Transcriptomic variation among six Arabidopsis thaliana accessions identified several novel genes controlling aluminium tolerance

Differences in the expression levels of aluminium (Al) tolerance genes are a known determinant of Al tolerance among plant varieties. We combined transcriptomic analysis of six Arabidopsis thaliana accessions with contrasting Al tolerance and a reverse genetic approach to identify Al-tolerance genes responsible for differences in Al tolerance between accession groups. Gene expression variation increased in the signal transduction process under Al stress and in growth-related processes in the absence of stress. Co-expression analysis and promoter single nucleotide polymorphism searching suggested that both trans-acting polymorphisms of Al signal transduction pathway and cis-acting polymorphisms in the promoter sequences caused the variations in gene expression associated with Al tolerance. Compared with the wild type, Al sensitivity increased in T-DNA knockout (KO) lines for five genes, including TARGET OF AVRB OPERATION1 (TAO1) and an unannotated gene (At5g22530). These were identified from 53 Al-inducible genes showing significantly higher expression in tolerant accessions than in sensitive accessions. These results indicate that the difference in transcriptional signalling is partly associated with the natural variation in Al tolerance in Arabidopsis. Our study also demonstrates the feasibility of comparative transcriptome analysis by using natural genetic variation for the identification of genes responsible for Al stress tolerance.

Large Differences in Gene Expression Responses to Drought and Heat Stress between Elite Barley Cultivar Scarlett and a Spanish Landrace

Drought causes important losses in crop production every season. Improvement for drought tolerance could take advantage of the diversity held in germplasm collections, much of which has not been incorporated yet into modern breeding. Spanish landraces constitute a promising resource for barley breeding, as they were widely grown until last century and still show good yielding ability under stress. Here, we study the transcriptome expression landscape in two genotypes, an outstanding Spanish landrace-derived inbred line (SBCC073) and a modern cultivar (Scarlett). Gene expression of adult plants after prolonged stresses, either drought or drought combined with heat, was monitored. Transcriptome of mature leaves presented little changes under severe drought, whereas abundant gene expression changes were observed under combined mild drought and heat. Developing inflorescences of SBCC073 exhibited mostly unaltered gene expression, whereas numerous changes were found in the same tissues for Scarlett. Genotypic differences in physiological traits and gene expression patterns confirmed the different behavior of landrace SBCC073 and cultivar Scarlett under abiotic stress, suggesting that they responded to stress following different strategies. A comparison with related studies in barley, addressing gene expression responses to drought, revealed common biological processes, but moderate agreement regarding individual differentially expressed transcripts. Special emphasis was put in the search of co-expressed genes and underlying common regulatory motifs. Overall, 11 transcription factors were identified, and one of them matched cis-regulatory motifs discovered upstream of co-expressed genes involved in those responses.

Evolution and Adaptation of Wild Emmer Wheat Populations to Biotic and Abiotic Stresses

The genetic bottlenecks associated with plant domestication and subsequent selection in man-made agroecosystems have limited the genetic diversity of modern crops and increased their vulnerability to environmental stresses. Wild emmer wheat, the tetraploid progenitor of domesticated wheat, distributed along a wide range of ecogeographical conditions in the Fertile Crescent, has valuable "left behind" adaptive diversity to multiple diseases and environmental stresses. The biotic and abiotic stress responses are conferred by series of genes and quantitative trait loci (QTLs) that control complex resistance pathways. The study of genetic diversity, genomic organization, expression profiles, protein structure and function of biotic and abiotic stress-resistance genes, and QTLs could shed light on the evolutionary history and adaptation mechanisms of wild emmer populations for their natural habitats. The continuous evolution and adaptation of wild emmer to the changing environment provide novel solutions that can contribute to safeguarding food for the rapidly growing human population.

Transcriptome profiling analysis for two Tibetan wild barley genotypes in responses to low nitrogen

BACKGROUND

Nitrogen (N) is the most common limiting factor for crop productivity worldwide. An effective approach to solve N deficiency is to develop low N (LN) tolerant crop cultivars. Tibetan annual wild barley is well-known for its wide genetic diversity and high tolerance to poor soil fertility. Up to date, no study has been done to illustrate the mechanism of LN tolerance underlying the wild barley at transcriptional level.

RESULTS

In this study, we employed Illumina RNA-Sequencing to determine the genotypic difference in transcriptome profile using two Tibetan wild barley genotypes differing in LN tolerance (XZ149, tolerant and XZ56, sensitive). A total of 1469 differentially expressed genes (DEGs) were identified in the two genotypes at 6 h and 48 h after LN treatment. Genetic difference existed in DEGs between XZ149 and XZ56, including transporters, transcription factors (TFs), kinases, antioxidant stress and hormone signaling related genes. Meanwhile, 695 LN tolerance-associated DEGs were mainly mapped to amino acid metabolism, starch and sucrose metabolism and secondary metabolism, and involved in transporter activity, antioxidant activities, and other gene ontology (GO). XZ149 had a higher capability of N absorption and use efficiency under LN stress than XZ56. The higher expression of nitrate transporters and energy-saving assimilation pattern could be attributed to its more N uptake and higher LN tolerance. In addition, auxin (IAA) and ethylene (ETH) response pathways may be also related to the genotypic difference in LN tolerance.

CONCLUSION

The responses of XZ149 and XZ56 to LN stress differed dramatically at transcriptional level. The identified candidate genes related to LN tolerance may provide new insights into comprehensive understanding of the genotypic difference in N utilization and LN tolerance.