High-quality genome assembly enables prediction of allele-specific gene expression in hybrid poplar.

Poplar (Populus) is a well-established model system for tree genomics and molecular breeding, and hybrid poplar is widely used in forest plantations. However, distinguishing its diploid homologous chromosomes is difficult, complicating advanced functional studies on specific alleles. In this study, we applied a trio-binning design and PacBio high-fidelity long-read sequencing to obtain haplotype-phased telomere-to-telomere genome assemblies for the 2 parents of the well-studied F1 hybrid "84K" (Populus alba × Populus tremula var. glandulosa). Almost all chromosomes, including the telomeres and centromeres, were completely assembled for each haplotype subgenome apart from 2 small gaps on one chromosome. By incorporating information from these haplotype assemblies and extensive RNA-seq data, we analyzed gene expression patterns between the 2 subgenomes and alleles. Transcription bias at the subgenome level was not uncovered, but extensive-expression differences were detected between alleles. We developed machine-learning (ML) models to predict allele-specific expression (ASE) with high accuracy and identified underlying genome features most highly influencing ASE. One of our models with 15 predictor variables achieved 77% accuracy on the training set and 74% accuracy on the testing set. ML models identified gene body CHG methylation, sequence divergence, and transposon occupancy both upstream and downstream of alleles as important factors for ASE. Our haplotype-phased genome assemblies and ML strategy highlight an avenue for functional studies in Populus and provide additional tools for studying ASE and heterosis in hybrids.

A simulation study comparing advanced marker-assisted selection with genomic selection in tree breeding programs.

Advances in DNA sequencing technologies allow the sequencing of whole genomes of thousands of individuals and provide several million single nucleotide polymorphisms (SNPs) per individual. These data combined with precise and high-throughput phenotyping enable genome-wide association studies (GWAS) and the identification of SNPs underlying traits with complex genetic architectures. The identified causal SNPs and estimated allelic effects could then be used for advanced marker-assisted selection (MAS) in breeding programs. But could such MAS compete with the broadly used genomic selection (GS)? This question is of particular interest for the lengthy tree breeding strategies. Here, with our new software "SNPscan breeder," we simulated a simple tree breeding program and compared the impact of different selection criteria on genetic gain and inbreeding. Further, we assessed different genetic architectures and different levels of kinship among individuals of the breeding population. Interestingly, apart from progeny testing, GS using gBLUP performed best under almost all simulated scenarios. MAS based on GWAS results outperformed GS only if the allelic effects were estimated in large populations (ca. 10,000 individuals) of unrelated individuals. Notably, GWAS using 3,000 extreme phenotypes performed as good as the use of 10,000 phenotypes. GS increased inbreeding and thus reduced genetic diversity more strongly compared to progeny testing and GWAS-based selection. We discuss the practical implications for tree breeding programs. In conclusion, our analyses further support the potential of GS for forest tree breeding and improvement, although MAS may gain relevance with decreasing sequencing costs in the future.

Preparation of Poly(vinyl Alcohol) Microparticles for Freeze Protection of Sensitive Fruit Crops.

Poly(vinyl alcohol) (PVA) displays ice recrystallization inhibition (IRI) properties as many antifreeze proteins found in cold tolerant organisms. The molecular architecture and composition (molecular weight and distribution of pendant OH and acetate groups) have been studied to improve the antifreezing properties of PVA, suggesting that the molecular architecture of PVA plays an important role in IRI activity. The present work deals with the preparation of PVA microparticles using an alkaline treatment. The effect of PVA molecular weight on the morphology and antifreezeing properties of PVA microparticles was investigated. The antifreezeing property of PVA microparticles on the susceptibility of flower bud tissues to freeze damage was also evaluated. The alkaline treatment of an aqueous PVA solution produced stable polymer chain aggregates with spherical shapes. The average size of the PVA microparticles increased significantly with the increasing molecular weight of the PVA macromolecule precursor. The PVA microparticles inhibited the growth of ice crystals and blocked ice growth at concentrations as low as 0.01 % w/v. The effect of impeding ice crystal growth by preventing the joining of adjacent ice crystals is attributed to the larger size of the PVA particles adsorbed on the ice surface compared to the aggregated PVA macromolecules in saline solution. The thermal hysteresis activity of PVA macromolecules and microparticles was not detected by differential scanning calorimetry analysis. The PVA microparticles reduced the incidence of freeze injuries in flower bud tissues by 55% and their application, considering the low toxicity of PVA, has a high potential for freeze protection in fruit crops.

ARR17 controls dioecy in Populus by repressing B-class MADS-box gene expression.

The number of dioecious species for which the genetic basis of sex determination has been resolved is rapidly increasing. Nevertheless, the molecular mechanisms downstream of the sex determinants remain largely elusive. Here, by RNA-sequencing early-flowering isogenic aspen (Populus tremula) lines differing exclusively for the sex switch gene ARR17, we show that a narrowly defined genetic network controls differential development of female and male flowers. Although ARR17 encodes a type-A response regulator supposedly involved in cytokinin (CK) hormone signalling, clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9-mediated arr17 knockout only affected the expression of a strikingly small number of genes, indicating a specific role in the regulation of floral development rather than a generic function in hormone signalling. Notably, the UNUSUAL FLORAL ORGANS (UFO) gene, encoding an F-box protein acting as a transcriptional cofactor with LEAFY (LFY) to activate B-class MADS-box gene expression, and the B-class gene PISTILLATA (PI), necessary for male floral organ development, were strongly de-repressed in the arr17 CRISPR mutants. Our data highlight a CK-independent role of the poplar response regulator ARR17 and further emphasize the minimal differences between female and male individuals. This article is part of the theme issue 'Sex determination and sex chromosome evolution in land plants'.

Sex determination and sex chromosome evolution in land plants.

Linnaeus's very first opus, written when he was 22 years old, dealt with the analogy that exists between plants and animals in how they 'propagate their species', and a revised version with a plate depicting the union of male and female Mercurialis annua plants became a foundational text on the sexuality of plants. The question how systems with separate males and females have evolved in sedentary organisms that appear ancestrally bisexual has fascinated biologists ever since. The phenomenon, termed dioecy, has important consequences for plant reproductive success and is of commercial interest since it affects seed quality and fruit production. This theme issue presents a series of articles that synthesize and challenge the current understanding of how plants achieve dioecy. The articles deal with a broad set of taxa, including Coccinia, Ginkgo, Mercurialis, Populus, Rumex and Silene, as well as overarching topics, such as the field's terminology, analogies with animal sex determination systems, evolutionary pathways to dioecy, dosage compensation, and the longevity of the two sexes. In this introduction, we focus on four topics, each addressed by several articles from different angles and with different conclusions. Our highlighting of unclear or controversial issues may help future studies to build on the current understanding and to ask new questions that will expand our knowledge of plant sexual systems. This article is part of the theme issue 'Sex determination and sex chromosome evolution in land plants'.

The Hda1 histone deacetylase limits divergent non-coding transcription and restricts transcription initiation frequency.

Nucleosome-depleted regions (NDRs) at gene promoters support initiation of RNA polymerase II transcription. Interestingly, transcription often initiates in both directions, resulting in an mRNA and a divergent non-coding (DNC) transcript of unclear purpose. Here, we characterized the genetic architecture and molecular mechanism of DNC transcription in budding yeast. Using high-throughput reverse genetic screens based on quantitative single-cell fluorescence measurements, we identified the Hda1 histone deacetylase complex (Hda1C) as a repressor of DNC transcription. Nascent transcription profiling showed a genome-wide role of Hda1C in repression of DNC transcription. Live-cell imaging of transcription revealed that mutations in the Hda3 subunit increased the frequency of DNC transcription. Hda1C contributed to decreased acetylation of histone H3 in DNC transcription regions, supporting DNC transcription repression by histone deacetylation. Our data support the interpretation that DNC transcription results as a consequence of the NDR-based architecture of eukaryotic promoters, but that it is governed by locus-specific repression to maintain genome fidelity.

Whole-genome sequencing and genome regions of special interest: Lessons from major histocompatibility complex, sex determination, and plant self-incompatibility.

Whole-genome sequencing of non-model organisms is now widely accessible and has allowed a range of questions in the field of molecular ecology to be investigated with greater power. However, some genomic regions that are of high biological interest remain problematic for assembly and data-handling. Three such regions are the major histocompatibility complex (MHC), sex-determining regions (SDRs) and the plant self-incompatibility locus (S-locus). Using these as examples, we illustrate the challenges of both assembling and resequencing these highly polymorphic regions and how bioinformatic and technological developments are enabling new approaches to their study. Mapping short-read sequences against multiple alternative references improves genotyping comprehensiveness at the S-locus thereby contributing to more accurate assessments of allelic frequencies. Long-read sequencing, producing reads of several tens to hundreds of kilobase pairs in length, facilitates the assembly of such regions as single sequences can span the multiple duplicated gene copies of the MHC region, and sequence through repetitive stretches and translocations in SDRs and S-locus haplotypes. These advances are adding value to short-read genome resequencing approaches by allowing, for example, more accurate haplotype phasing across longer regions. Finally, we assessed further technical improvements, such as nanopore adaptive sequencing and bioinformatic tools using pangenomes, which have the potential to further expand our knowledge of a number of genomic regions that remain challenging to study with classical resequencing approaches.

Plant sex chromosomes defy evolutionary models of expanding recombination suppression and genetic degeneration.

Hundreds of land plant lineages have independently evolved separate sexes in either gametophytes (dioicy) or sporophytes (dioecy), but 43% of all dioecious angiosperms are found in just 34 entirely dioecious clades, suggesting that their mode of sex determination evolved a long time ago. Here, we review recent insights on the molecular mechanisms that underlie the evolutionary change from individuals that each produce male and female gametes to individuals specializing in the production of just one type of gamete. The canonical model of sex chromosome evolution in plants predicts that two sex-determining genes will become linked in a sex-determining region (SDR), followed by expanding recombination suppression, chromosome differentiation and, ultimately, degeneration. Experimental work, however, is showing that single genes function as master regulators in model systems, such as the liverwort Marchantia and the angiosperms Diospyros and Populus. In Populus, this type of regulatory function has been demonstrated by genome editing. In other systems, including Actinidia, Asparagus and Vitis, two coinherited factors appear to independently regulate female and male function, yet sex chromosome differentiation has remained low. We discuss the best-understood systems and evolutionary pathways to dioecy, and present a meta-analysis of the sizes and ages of SDRs. We propose that limited sexual conflict explains why most SDRs are small and sex chromosomes remain homomorphic. It appears that models of increasing recombination suppression with age do not apply because selection favours mechanisms in which sex determination depends on minimal differences, keeping it surgically precise.

Default Sex and Single Gene Sex Determination in Dioecious Plants.

A well-established hypothesis for the evolution of dioecy involves two genes linked at a sex-determining region (SDR). Recently there has been increased interest in possible single gene sex determination. Work in Populus has finally provided direct experimental evidence for single gene sex determination in plants using CRISPR-Cas9 to knock out a single gene and convert individuals from female to male. In poplar, the feminizing factor popARR17 acts as a "master regulator", analogous to the mammalian masculinizing factor SRY. The production of fully functional males from females by a simple single gene knockout is experimental evidence that an antagonistic male-determining factor does not exist in Populus. Mammals have a "default sex" (female), as do poplar trees (Populus), although the default sex in poplars is male. The occurrence of single gene sex determination with a default sex may be much commoner in plants than hitherto expected, especially when dioecy evolved via monoecy. The master regulator does not even need to be at the SDR (although it may be). In most poplars the feminizing factor popARR17 is not at the SDR, but instead a negative regulator of it. So far there is little information on how high-level regulators are connected to floral phenotype. A model is presented of how sex-determining genes could lead to different floral morphologies via MADS-box floral developmental genes.

A single gene underlies the dynamic evolution of poplar sex determination.

Although hundreds of plant lineages have independently evolved dioecy (that is, separation of the sexes), the underlying genetic basis remains largely elusive1. Here we show that diverse poplar species carry partial duplicates of the ARABIDOPSIS RESPONSE REGULATOR 17 (ARR17) orthologue in the male-specific region of the Y chromosome. These duplicates give rise to small RNAs apparently causing male-specific DNA methylation and silencing of the ARR17 gene. CRISPR-Cas9-induced mutations demonstrate that ARR17 functions as a sex switch, triggering female development when on and male development when off. Despite repeated turnover events, including a transition from the XY system to a ZW system, the sex-specific regulation of ARR17 is conserved across the poplar genus and probably beyond. Our data reveal how a single-gene-based mechanism of dioecy can enable highly dynamic sex-linked regions and contribute to maintaining recombination and integrity of sex chromosomes.

ENO regulates tomato fruit size through the floral meristem development network.

A dramatic evolution of fruit size has accompanied the domestication and improvement of fruit-bearing crop species. In tomato (Solanum lycopersicum), naturally occurring cis-regulatory mutations in the genes of the CLAVATA-WUSCHEL signaling pathway have led to a significant increase in fruit size generating enlarged meristems that lead to flowers with extra organs and bigger fruits. In this work, by combining mapping-by-sequencing and CRISPR/Cas9 genome editing methods, we isolated EXCESSIVE NUMBER OF FLORAL ORGANS (ENO), an AP2/ERF transcription factor which regulates floral meristem activity. Thus, the ENO gene mutation gives rise to plants that yield larger multilocular fruits due to an increased size of the floral meristem. Genetic analyses indicate that eno exhibits synergistic effects with mutations at the LOCULE NUMBER (encoding SlWUS) and FASCIATED (encoding SlCLV3) loci, two central players in the evolution of fruit size in the domestication of cultivated tomatoes. Our findings reveal that an eno mutation causes a substantial expansion of SlWUS expression domains in a flower-specific manner. In vitro binding results show that ENO is able to interact with the GGC-box cis-regulatory element within the SlWUS promoter region, suggesting that ENO directly regulates SlWUS expression domains to maintain floral stem-cell homeostasis. Furthermore, the study of natural allelic variation of the ENO locus proved that a cis-regulatory mutation in the promoter of ENO had been targeted by positive selection during the domestication process, setting up the background for significant increases in fruit locule number and fruit size in modern tomatoes.

A Novel Hydrocolloid Film Based on Pectin, Starch and Gunnera tinctoria and Ugni molinae Plant Extracts for Wound Dressing Applications.

BACKGROUND: The biodegradable and biocompatible nature of pectin-based films is of particular interest in wound dressing applications, due to its non-toxicity, pH-sensitivity and gelling activity. An approach to improve the mechanical properties, the release profile of bioactive compounds as well as the performance in wet environments of pectin-based films is mixing with other biopolymers. OBJECTIVE: To prepare hydrocolloid films based on crosslinked pectin / starch blend loaded with bioactive extracts from leaves of G. tinctoria and U. molinae with controlled release of bioactive compounds and healing property. METHODS: The hydrocolloid films were characterized by FTIR, SEM, and TGA-FTIR techniques and their tensile properties, water uptake, and polyphenolic release profile in aqueous media were evaluated. The dermal anti inflammatory activity of the hydrocolloid films was assessed by the mouse ear inflammation test. The wound healing property of the loaded hydrocolloid films was explored in a rat model and in a clinical trial (sacrum pressure ulcer). RESULTS: The films showed an adequate water-uptake capacity between 100-160%. The release of active compounds from the hydrocolloid films followed the Korsmeyer-Peppas equation. The mechanical properties of hydrocolloid films were not affected by the plant extracts within the concentration range used. The incorporation of the bioactive extracts in the polysaccharide films inhibited the topical edematous response by about 50%. The topical application of the loaded hydrocolloid film on the pressure ulcer is completely closed after 17 days without showing any adverse reaction. CONCLUSION: A novel hydrocolloid matrix was produced from crosslinked starch-pectin, which exhibited suitable chemical-physical properties to be used as a carrier of plant extracts with wound healing properties.

The Diversity and Dynamics of Sex Determination in Dioecious Plants.

The diversity of inflorescences among flowering plants is captivating. Such charm is not only due to the variety of sizes, shapes, colors, and flowers displayed, but also to the range of reproductive systems. For instance, hermaphrodites occur abundantly throughout the plant kingdom with both stamens and carpels within the same flower. Nevertheless, 10% of flowering plants have separate unisexual flowers, either in different locations of the same individual (monoecy) or on different individuals (dioecy). Despite their rarity, dioecious plants provide an excellent opportunity to investigate the mechanisms involved in sex expression and the evolution of sex-determining regions (SDRs) and sex chromosomes. The SDRs and the evolution of dioecy have been studied in many species ranging from Ginkgo to important fruit crops. Some of these studies, for example in asparagus or kiwifruit, identified two sex-determining genes within the non-recombining SDR and may thus be consistent with the classical model for the evolution of dioecy from hermaphroditism via gynodioecy, that predicts two successive mutations, the first one affecting male and the second one female function, becoming linked in a region of suppressed recombination. On the other hand, aided by genome sequencing and gene editing, single factor sex determination has emerged in other species, such as persimmon or poplar. Despite the diversity of sex-determining mechanisms, a tentative comparative analysis of the known sex-determining genes and candidates in different species suggests that similar genes and pathways may be employed repeatedly for the evolution of dioecy. The cytokinin signaling pathway appears important for sex determination in several species regardless of the underlying genetic system. Additionally, tapetum-related genes often seem to act as male-promoting factors when sex is determined via two genes. We present a unified model that synthesizes the genetic networks of sex determination in monoecious and dioecious plants and will support the generation of hypothesis regarding candidate sex determinants in future studies.

Correction to: Long-term study of a subdioecious Populus ×canescens family reveals sex lability of females and reproduction behaviour of cosexual plants.

Table 4 in the original publication reports incomplete genotype names in the column "Cross" and wrong codes in the column "Generation".

Long-term study of a subdioecious Populus ×canescens family reveals sex lability of females and reproduction behaviour of cosexual plants.

KEY MESSAGE: Cosexual Populus ×canescens plants are inconstant females with life course plasticity of sex phenotype and can reproduce by selfing. Populus species are dioecious, but deviations from dioecy are reported in some cases. The objectives of this study were to investigate the phenotypic expression and the inheritance of subdioecy in a Populus ×canescens pedigree. The F1 progeny was monitored for sex during 14 years. Thirty per cent of individuals expressed deviations from dioecy and long-term plasticity of sex. Some plants started flowering as male, then became cosexual, and finally turned female. Two cosexual individuals were self-pollinated and generated a selfed progeny markedly impaired by inbreeding depression, but able to reproduce by outcrossing. Sex segregation of the F1 progeny statistically fitted the expected ratio 1:2:1 (female:male:cosexual). By analysis of DNA markers, the cosexual individuals were genetically clustered with the females. The segregation ratio and the genetic profile indicated that cosexual plants were female with altered sex phenotype. Linkage analysis identified a putative sex-determining region with suppressed recombination on chromosome 19 of the male Populus tremula parent. The male sex trait was linked to the pericentromeric region of the P. tremula chromosome 19, whereas the cosexual trait was linked to chromosome 19 of the female Populus alba parent. A genetic model is proposed to explain inheritance and phenotypic expression of sex.

RNA-seq of eight different poplar clones reveals conserved up-regulation of gene expression in response to insect herbivory.

BACKGROUND: Herbivorous insects can have a profound impact on plant growth performance. In some years, canopy damage in poplar plantations exceeds 50% of the total leaf surface, thereby possibly compromising carbon fixation and biomass yield. To assess the transcriptional response of elite poplar clones to insect feeding and to test whether this response varies between different genotypes, we performed an RNA-sequencing experiment. We deeply sequenced the transcriptomes of eight elite clones belonging to three poplar species (Populus trichocarpa, P. nigra and P. maximowiczii), under Phratora vitellinae feeding and control conditions. This allowed us to precisely quantify transcript levels of about 24,000 expressed genes. RESULTS: Our data reveal a striking overall up-regulation of gene expression under insect attack in all eight poplar clones studied. The up-regulated genes were markedly enriched for the biological process 'regulation of transcription' indicating a highly concerted restructuring of the transcriptome. A search for potential cis-regulatory elements (CREs) that may be involved in this process identified the G-box (CACGTG) as the most significant motif in the promoters of the induced genes. In line with the role of the G-box in jasmonate (JA)-mediated activation of gene expression by MYC2, several genes involved in JA biosynthesis and signaling were up-regulated in our dataset. A co-expression network analysis additionally highlighted WRKY transcription factors. Within the most prominent expression module, WRKYs were strongly overrepresented and occupied several network hubs. Finally, the insect-induced genes comprised several protein families known to be involved in plant defenses, e.g. cytochrome P450s, chitinases and protease inhibitors. CONCLUSIONS: Our data represent a comprehensive characterization of the transcriptional response of selected elite poplar clones to insect herbivory. Our results suggest that the concerted up-regulation of gene expression is controlled by JA signaling and WRKY transcription factors, and activates several defense mechanisms. Our data highlight potential targets of selection and may thus contribute to breeding insect-resistant poplar clones.

Relation between concentration fluctuations and dynamical heterogeneities in binary glass-forming liquids: A molecular dynamics simulation study.

We perform molecular dynamics simulations of binary glass-forming liquids with high dynamical contrast. In detail, we mix water-like molecules with various polarities, which exhibit strongly different dynamics but very similar structures as neat liquids. It is found that these mixtures show concentration fluctuations, which strongly grow upon cooling towards an expected mixing-demixing phase transition. Because of the developing microphase segregation, the composition and temperature ranges accessible to our approach are limited, but the studied mixtures are ideal models to ascertain the relation between concentration fluctuations and molecular dynamics, in particular, the glassy slowdown of binary liquids, including aqueous solutions. We find that the dynamics of the components decouple upon cooling. While the dynamics of the slow component resembles that of neat glass formers, that of the fast component shows different features, including sub-linear diffusion on length and time scales beyond that of local particle cages and quasi-logarithmic decays of correlation functions. We show that this "anomalous dynamics" of the fast component receives a simple explanation based on a strong dependence of the particle mobility on the local concentration in combination with microphase segregation and dynamical contrast. Explicitly, spatially heterogeneous dynamics, which are more extended and tenacious than in neat glass formers, result from growing concentration fluctuations upon cooling. Due to a high dynamical asymmetry of the components at low temperatures, the slow species sustains the microphase segregation in the mixture, while the fast component moves along the concentration and mobility gradients, leading to remarkable dynamical features.

High Level of Conservation of Mitochondrial RNA Editing Sites Among Four Populus Species.

RNA editing occurs in the endosymbiont organelles of higher plants as C-to-U conversions of defined nucleotides. The availability of large quantities of RNA sequencing data makes it possible to identify RNA editing sites and to quantify their editing extent. We have investigated RNA editing in 34 protein-coding mitochondrial transcripts of four Populus species, a genus noteworthy for its remarkably small number of RNA editing sites compared to other angiosperms. 27 of these transcripts were subject to RNA editing in at least one species. In total, 355 RNA editing sites were identified with high confidence, their editing extents ranging from 10 to 100%. The most heavily edited transcripts were ccmB with the highest density of RNA editing sites (53.7 sites / kb) and ccmFn with the highest number of sites (39 sites). Most of the editing events are at position 1 or 2 of the codons, usually altering the encoded amino acid, and are highly conserved among the species, also with regard to their editing extent. However, one SNP was found in the newly sequenced and annotated mitochondrial genome of P. alba resulting in the loss of an RNA editing site compared to P. tremula and P. davidiana This SNP causes a C-to-T transition and an amino acid exchange from Ser to Phe, highlighting the widely discussed role of RNA editing in compensating mutations.

Mutations in EID1 and LNK2 caused light-conditional clock deceleration during tomato domestication.

Circadian period and phase of cultivated tomato (Solanum lycopersicum) were changed during domestication, likely adapting the species to its new agricultural environments. Whereas the delayed circadian phase is mainly caused by allelic variation of EID1, the genetic basis of the long circadian period has remained elusive. Here we show that a partial deletion of the clock gene LNK2 is responsible for the period lengthening in cultivated tomatoes. We use resequencing data to phylogenetically classify hundreds of tomato accessions and investigate the evolution of the eid1 and lnk2 mutations along successive domestication steps. We reveal signatures of selection across the genomic region of LNK2 and different patterns of fixation of the mutant alleles. Strikingly, LNK2 and EID1 are both involved in light input to the circadian clock, indicating that domestication specifically targeted this input pathway. In line with this, we show that the clock deceleration in the cultivated tomato is light-dependent and requires the phytochrome B1 photoreceptor. Such conditional variation in circadian rhythms may be key for latitudinal adaptation in a variety of species, including crop plants and livestock.

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.