From Neglecting to Including Cultivar-Specific Per Se Temperature Responses: Extending the Concept of Thermal Time in Field Crops.

Predicting plant development, a longstanding goal in plant physiology, involves 2 interwoven components: continuous growth and the progression of growth stages (phenology). Current models for winter wheat and soybean assume species-level growth responses to temperature. We challenge this assumption, suggesting that cultivar-specific temperature responses substantially affect phenology. To investigate, we collected field-based growth and phenology data in winter wheat and soybean over multiple years. We used diverse models, from linear to neural networks, to assess growth responses to temperature at various trait and covariate levels. Cultivar-specific nonlinear models best explained phenology-related cultivar-environment interactions. With cultivar-specific models, additional relations to other stressors than temperature were found. The availability of the presented field phenotyping tools allows incorporating cultivar-specific temperature response functions in future plant physiology studies, which will deepen our understanding of key factors that influence plant development. Consequently, this work has implications for crop breeding and cultivation under adverse climatic conditions.

High-resolution bulked segregant analysis enables candidate gene identification for bacterial wilt resistance in Italian ryegrass (Lolium multiflorum Lam.).

Bacterial wilt, caused by Xanthomonas translucens pv. graminis (Xtg), is a serious disease of economically important forage grasses, including Italian ryegrass (Lolium multiflorum Lam.). A major QTL for resistance to Xtg was previously identified, but the precise location as well as the genetic factors underlying the resistance are yet to be determined. To this end, we applied a bulked segregant analysis (BSA) approach, using whole-genome deep sequencing of pools of the most resistant and most susceptible individuals of a large (n = 7484) biparental F2 population segregating for resistance to Xtg. Using chromosome-level genome assemblies as references, we were able to define a ~300 kb region highly associated with resistance on pseudo-chromosome 4. Further investigation of this region revealed multiple genes with a known role in disease resistance, including genes encoding for Pik2-like disease resistance proteins, cysteine-rich kinases, and RGA4- and RGA5-like disease resistance proteins. Investigation of allele frequencies in the pools and comparative genome analysis in the grandparents of the F2 population revealed that some of these genes contain variants with allele frequencies that correspond to the expected heterozygosity in the resistant grandparent. This study emphasizes the efficacy of combining BSA studies in very large populations with whole genome deep sequencing and high-quality genome assemblies to pinpoint regions associated with a binary trait of interest and accurately define a small set of candidate genes. Furthermore, markers identified in this region hold significant potential for marker-assisted breeding strategies to breed resistance to Xtg in Italian ryegrass cultivars more efficiently.

ChromaX: a fast and scalable breeding program simulator.

SUMMARY: ChromaX is a Python library that enables the simulation of genetic recombination, genomic estimated breeding value calculations, and selection processes. By utilizing GPU processing, it can perform these simulations up to two orders of magnitude faster than existing tools with standard hardware. This offers breeders and scientists new opportunities to simulate genetic gain and optimize breeding schemes. AVAILABILITY AND IMPLEMENTATION: The documentation is available at https://chromax.readthedocs.io. The code is available at https://github.com/kora-labs/chromax.

Genome-wide association study of lipase and esterase in wholegrain wheat flour (Triticum aestivum L.).

Lipase activity is one of the main causes of the lipid rancidity in wholegrain wheat flour, leading to its short shelf life. Genetically diverse wheat germplasm offers potential for the selection of wheat cultivars with low lipase activity for stable wholegrain end use. This study evaluated 300 European wheat cultivars harvested in 2015 and 2016 on the genetic association of lipase and esterase activities in wholegrain wheat flour. Esterase and lipase activities in wholegrain flour were measured photometrically with p-nitrophenyl butyrate and p-nitrophenyl palmitate as substrates, respectively. Both enzyme activities showed wide ranges among all cultivars within each year, with differences up to 2.5-fold. The two years showed low correlations between each other, indicating a large environmental impact on the enzyme activities. Cultivars 'Julius' and 'Bueno' were suggested to be better suited for stable wholegrain products, as they had consistently low esterase and lipase activities compared to the other cultivars. A genome-wide association study revealed associations with single nucleotide polymorphisms in genes located on the high-quality wheat genome sequence of the International Wheat Genome Sequencing Consortium. Eight and four candidate genes were tentatively proposed to be associated to esterase and lipase activity, respectively, in wholegrain flour. Our work shows esterase and lipase activities from a new perspective, that combines reverse genetics to understand the underlying causes. This study outlines the possibilities and limitations to improve lipid stability of wholegrain wheat by genomics-assisted breeding methods, thereby offering new opportunities to optimize the quality of wholegrain wheat flour and wholegrain products.

Fine-Mapping and Comparative Genomic Analysis Reveal the Gene Composition at the S and Z Self-incompatibility Loci in Grasses.

Self-incompatibility (SI) is a genetic mechanism of hermaphroditic plants to prevent inbreeding after self-pollination. Allogamous Poaceae species exhibit a unique gametophytic SI system controlled by two multi-allelic and independent loci, S and Z. Despite intense research efforts in the last decades, the genes that determine the initial recognition mechanism are yet to be identified. Here, we report the fine-mapping of the Z-locus in perennial ryegrass (Lolium perenne L.) and provide evidence that the pollen and stigma components are determined by two genes encoding DUF247 domain proteins (ZDUF247-I and ZDUF247-II) and the gene sZ, respectively. The pollen and stigma determinants are located side-by-side and were genetically linked in 10,245 individuals of two independent mapping populations segregating for Z. Moreover, they exhibited high allelic diversity as well as tissue-specific gene expression, matching the expected characteristics of SI determinants known from other systems. Revisiting the S-locus using the latest high-quality whole-genome assemblies revealed a similar gene composition and structure as found for Z, supporting the hypothesis of a duplicated origin of the two-locus SI system of grasses. Ultimately, comparative genomic analyses across a wide range of self-compatible and self-incompatible Poaceae species revealed that the absence of a functional copy of at least one of the six putative SI determinants is accompanied by a self-compatible phenotype. Our study provides new insights into the origin and evolution of the unique gametophytic SI system in one of the largest and economically most important plant families.

Callus Induction from Diverse Explants and Genotypes Enables Robust Transformation of Perennial Ryegrass (Lolium perenne L.).

Genetic transformation of perennial ryegrass (Lolium perenne L.) is critical for fundamental and translational research in this important grass species. It often relies on Agrobacterium-mediated transformation of callus tissue. However, callus induction is restricted to a few genotypes that respond well to tissue culture. Here, we report callus induction from different perennial ryegrass genotypes and explants, such as shoot tips, seeds, and anthers, which were transformed with several plasmids for functional genomics. ß-glucuronidase (GUS) histochemical staining showed the LmdsRNAbp promoter sequence was active in stigmas, spikelets, anthers, and leaves. We also transformed calli with plasmids allowing gene silencing and gene knock-out using RNA interference and CRISPR/Cas9, respectively, for which genotypic and phenotypic investigations are ongoing. Using 19 different constructs, 262 transgenic events were regenerated. Moreover, the protocol regenerated a doubled haploid transgenic event from anther-derived calli. This work provides a proof-of-concept method for expanding the range of genotypes amenable to transformation, thus, serving research and breeding initiatives to improve this important grass crop for forage and recreation.

Identification of Candidate Genes for Self-Compatibility in Perennial Ryegrass (Lolium perenne L.).

Self-incompatibility (SI) is a genetic mechanism preventing self-pollination in ~40% of plant species. Two multiallelic loci, called S and Z, control the gametophytic SI system of the grass family (Poaceae), which contains all major forage grasses. Loci independent from S and Z have been reported to disrupt SI and lead to self-compatibility (SC). A locus causing SC in perennial ryegrass (Lolium perenne L.) was previously mapped on linkage group (LG) 5 in an F2 population segregating for SC. Using a subset of the same population (n = 68), we first performed low-resolution quantitative trait locus (QTL) mapping to exclude the presence of additional, previously undetected contributors to SC. The previously reported QTL on LG 5 explained 38.4% of the phenotypic variation, and no significant contribution from other genomic regions was found. This was verified by the presence of significantly distorted markers in the region overlapping with the QTL. Second, we fine mapped the QTL to 0.26 centimorgan (cM) using additional 2,056 plants and 23 novel sequence-based markers. Using Italian ryegrass (Lolium multiflorum Lam.) genome assembly as a reference, the markers flanking SC were estimated to span a ~3 Mb region encoding for 57 predicted genes. Among these, seven genes were proposed as relevant candidate genes based on their annotation and function described in previous studies. Our study is a step forward to identify SC genes in forage grasses and provides diagnostic markers for marker-assisted introgression of SC into elite germplasm.

Ultralong Oxford Nanopore Reads Enable the Development of a Reference-Grade Perennial Ryegrass Genome Assembly.

Despite the progress made in DNA sequencing over the last decade, reconstructing telomere-to-telomere genome assemblies of large and repeat-rich eukaryotic genomes is still difficult. More accurate basecalls or longer reads could address this issue, but no current sequencing platform can provide both simultaneously. Perennial ryegrass (Lolium perenne L.) is an example of an important species for which the lack of a reference genome assembly hindered a swift adoption of genomics-based methods into breeding programs. To fill this gap, we optimized the Oxford Nanopore Technologies' sequencing protocol, obtaining sequencing reads with an N50 of 62 kb-a very high value for a plant sample. The assembly of such reads produced a highly complete (2.3 of 2.7 Gb), correct (QV 45), and contiguous (contig N50 and N90 11.74 and 3.34 Mb, respectively) genome assembly. We show how read length was key in determining the assembly contiguity. Sequence annotation revealed the dominance of transposable elements and repeated sequences (81.6% of the assembly) and identified 38,868 protein coding genes. Almost 90% of the bases could be anchored to seven pseudomolecules, providing the first high-quality haploid reference assembly for perennial ryegrass. This protocol will enable producing longer Oxford Nanopore Technology reads for more plant samples and ushering forage grasses into modern genomics-assisted breeding programs.

Rethinking temperature effects on leaf growth, gene expression and metabolism: Diel variation matters.

Plants have evolved to grow under prominently fluctuating environmental conditions. In experiments under controlled conditions, temperature is often set to artificial, binary regimes with constant values at day and at night. This study investigated how such a diel (24 hr) temperature regime affects leaf growth, carbohydrate metabolism and gene expression, compared to a temperature regime with a field-like gradual increase and decline throughout 24 hr. Soybean (Glycine max) was grown under two contrasting diel temperature treatments. Leaf growth was measured in high temporal resolution. Periodical measurements were performed of carbohydrate concentrations, carbon isotopes as well as the transcriptome by RNA sequencing. Leaf growth activity peaked at different times under the two treatments, which cannot be explained intuitively. Under field-like temperature conditions, leaf growth followed temperature and peaked in the afternoon, whereas in the binary temperature regime, growth increased at night and decreased during daytime. Differential gene expression data suggest that a synchronization of cell division activity seems to be evoked in the binary temperature regime. Overall, the results show that the coordination of a wide range of metabolic processes is markedly affected by the diel variation of temperature, which emphasizes the importance of realistic environmental settings in controlled condition experiments.

Temperature response of wheat affects final height and the timing of stem elongation under field conditions.

In wheat, temperature affects the timing and intensity of stem elongation. Genetic variation for this process is therefore important for adaptation. This study investigates the genetic response to temperature fluctuations during stem elongation and its relationship to phenology and height. Canopy height of 315 wheat genotypes (GABI wheat panel) was scanned twice weekly in the field phenotyping platform (FIP) of ETH Zurich using a LIDAR. Temperature response was modelled using linear regressions between stem elongation and mean temperature in each measurement interval. This led to a temperature-responsive (slope) and a temperature-irresponsive (intercept) component. The temperature response was highly heritable (H2=0.81) and positively related to a later start and end of stem elongation as well as final height. Genome-wide association mapping revealed three temperature-responsive and four temperature-irresponsive quantitative trait loci (QTLs). Furthermore, putative candidate genes for temperature-responsive QTLs were frequently related to the flowering pathway in Arabidopsis thaliana, whereas temperature-irresponsive QTLs corresponded to growth and reduced height genes. In combination with Rht and Ppd alleles, these loci, together with the loci for the timing of stem elongation, accounted for 71% of the variability in height. This demonstrates how high-throughput field phenotyping combined with environmental covariates can contribute to a smarter selection of climate-resilient crops.

Genome-Wide Association Study to Identify Candidate Loci for Biomass Formation Under Water Deficit in Perennial Ryegrass.

Global warming is predicted to impact many agricultural areas, which will suffer from reduced water availability. Due to precipitation changes, mild summer droughts are expected to become more frequent, even in temperate regions. For perennial ryegrass (Lolium perenne L.), an important forage grass of the Poaceae family, leaf growth is a crucial factor determining biomass accumulation and hence forage yield. Although leaf elongation has been shown to be temperature-dependent under normal conditions, the genetic regulation of leaf growth under water deficit in perennial ryegrass is poorly understood. Herein, we evaluated the response to water deprivation in a diverse panel of perennial ryegrass genotypes, employing a high-precision phenotyping platform. The study revealed phenotypic variation for growth-related traits and significant (P < 0.05) differences in leaf growth under normal conditions within the subgroups of turf and forage type cultivars. The phenotypic data was combined with genotypic variants identified using genotyping-by-sequencing to conduct a genome-wide association study (GWAS). Using GWAS, we identified DNA polymorphisms significantly associated with leaf growth reduction under water deprivation. These polymorphisms were adjacent to genes predicted to encode for phytochrome B and a MYB41 transcription factor. The result obtained in the present study will increase our understanding on the complex molecular mechanisms involved in plant growth under water deficit. Moreover, the single nucleotide polymorphism (SNP) markers identified will serve as a valuable resource in future breeding programs to select for enhanced biomass formation under mild summer drought conditions.

Phenotyping a Dynamic Trait: Leaf Growth of Perennial Ryegrass Under Water Limiting Conditions.

Water limitation is one of the major factors reducing crop productivity worldwide. In order to develop efficient breeding strategies to improve drought tolerance, accurate methods to identify when a plant reduces growth as a consequence of water deficit have yet to be established. In perennial ryegrass (Lolium perenne L.), an important forage grass of the Poaceae family, leaf elongation is a key factor determining plant growth and hence forage yield. Although leaf elongation has been shown to be temperature-dependent under non-stress conditions, the impact of water limitation on leaf elongation in perennial ryegrass is poorly understood. We describe a method for quantifying tolerance to water deficit based on leaf elongation in relation to temperature and soil moisture in perennial ryegrass. With decreasing soil moisture, three growth response phases were identified: first, a "normal" phase where growth is mainly determined by temperature, second a "slow" phase where leaf elongation decreases proportionally to soil water potential and third an "arrest" phase where leaf growth terminates. A custom R function was able to quantify the points which demarcate these phases and can be used to describe the response of plants to water deficit. Applied to different perennial ryegrass genotypes, this function revealed significant genotypic variation in the response of leaf growth to temperature and soil moisture. Dynamic phenotyping of leaf elongation can be used as a tool to accurately quantify tolerance to water deficit in perennial ryegrass and to improve this trait by breeding. Moreover, the tools presented here are applicable to study the plant response to other stresses in species with linear, graminoid leaf morphology.

Precision Phenotyping Reveals Novel Loci for Quantitative Resistance to Septoria Tritici Blotch.

Accurate, high-throughput phenotyping for quantitative traits is a limiting factor for progress in plant breeding. We developed an automated image analysis to measure quantitative resistance to septoria tritici blotch (STB), a globally important wheat disease, enabling identification of small chromosome intervals containing plausible candidate genes for STB resistance. 335 winter wheat cultivars were included in a replicated field experiment that experienced natural epidemic development by a highly diverse but fungicide-resistant pathogen population. More than 5.4 million automatically generated phenotypes were associated with 13,648 SNP markers to perform the GWAS. We identified 26 chromosome intervals explaining 1.9-10.6% of the variance associated with four independent resistance traits. Sixteen of the intervals overlapped with known STB resistance intervals, suggesting that our phenotyping approach can identify simultaneously (i.e., in a single experiment) many previously defined STB resistance intervals. Seventeen of the intervals were less than 5 Mbp in size and encoded only 173 genes, including many genes associated with disease resistance. Five intervals contained four or fewer genes, providing high priority targets for functional validation. Ten chromosome intervals were not previously associated with STB resistance, perhaps representing resistance to pathogen strains that had not been tested in earlier experiments. The SNP markers associated with these chromosome intervals can be used to recombine different forms of quantitative STB resistance that are likely to be more durable than pyramids of major resistance genes. Our experiment illustrates how high-throughput automated phenotyping can accelerate breeding for quantitative disease resistance.

Pooled DNA sequencing to identify SNPs associated with a major QTL for bacterial wilt resistance in Italian ryegrass (Lolium multiflorum Lam.).

KEY MESSAGE: SNPs and candidate genes associated with bacterial wilt resistance in Italian ryegrass were identified by sequencing the parental plants and pooled F1 progeny of a segregating population. Italian ryegrass (Lolium multiflorum Lam.) is one of the most important forage grass species in temperate regions. Its yield, quality and persistency can significantly be reduced by bacterial wilt, a serious disease caused by Xanthomonas translucens pv. graminis. Although a major QTL for bacterial wilt resistance has previously been reported, detailed knowledge on underlying genes and DNA markers to allow for efficient resistance breeding strategies is currently not available. We used pooled DNA sequencing to characterize a major QTL for bacterial wilt resistance of Italian ryegrass and to develop inexpensive sequence-based markers to efficiently target resistance alleles for marker-assisted recurrent selection. From the mapping population segregating for the QTL, DNA of 44 of the most resistant and 44 of the most susceptible F1 individuals was pooled and sequenced using the Illumina HiSeq 2000 platform. Allele frequencies of 18 × 106 single nucleotide polymorphisms (SNP) were determined in the resistant and susceptible pool. A total of 271 SNPs on 140 scaffold sequences of the reference parental genome showed significantly different allele frequencies in both pools. We converted 44 selected SNPs to KASP™ markers, genetically mapped these proximal to the major QTL and thus validated their association with bacterial wilt resistance. This study highlights the power of pooled DNA sequencing to efficiently target binary traits in biparental mapping populations. It delivers genome sequence data, SNP markers and potential candidate genes which will allow to implement marker-assisted strategies to fix bacterial wilt resistance in outcrossing breeding populations of Italian ryegrass.

Genetic Loci Governing Androgenic Capacity in Perennial Ryegrass (Lolium perenne L.).

Immature pollen can be induced to switch developmental pathways from gametogenesis to embryogenesis and subsequently regenerate into homozygous, diploid plants. Such androgenic production of doubled haploids is particularly useful for species where inbreeding is hampered by effective self-incompatibility systems. Therefore, increasing the generally low androgenic capacity of perennial ryegrass (Lolium perenne L.) germplasm would enable the efficient production of homozygous plant material, so that a more effective exploitation of heterosis through hybrid breeding schemes can be realized. Here, we present the results of a genome-wide association study in a heterozygous, multiparental population of perennial ryegrass (n = 391) segregating for androgenic capacity. Genotyping-by-sequencing was used to interrogate gene- dense genomic regions and revealed over 1,100 polymorphic sites. Between one and 10 quantitative trait loci (QTL) were identified for anther response, embryo and total plant production, green and albino plant production and regeneration. Most traits were under polygenic control, although a major QTL on linkage group 5 was associated with green plant regeneration. Distinct genetic factors seem to affect green and albino plant recovery. Two intriguing candidate genes, encoding chromatin binding domains of the developmental phase transition regulator, Polycomb Repressive Complex 2, were identified. Our results shed the first light on the molecular mechanisms behind perennial ryegrass microspore embryogenesis and enable marker-assisted introgression of androgenic capacity into recalcitrant germplasm of this forage crop of global significance.

A Novel Multivariate Approach to Phenotyping and Association Mapping of Multi-Locus Gametophytic Self-Incompatibility Reveals S, Z, and Other Loci in a Perennial Ryegrass (Poaceae) Population.

Self-incompatibility (SI) is a mechanism that many flowering plants employ to prevent fertilisation by self- and self-like pollen ensuring heterozygosity and hybrid vigour. Although a number of single locus mechanisms have been characterised in detail, no multi-locus systems have been fully elucidated. Historically, examples of the genetic analysis of multi-locus SI, to make analysis tractable, are either made on the progeny of bi-parental crosses, where the number of alleles at each locus is restricted, or on crosses prepared in such a way that only one of the SI loci segregates. Perennial ryegrass (Lolium perenne L.) possesses a well-documented two locus (S and Z) gametophytic incompatibility system. A more universal, realistic proof of principle study was conducted in a perennial ryegrass population in which allelic and non-allelic diversity was not artificially restricted. A complex pattern of pollinations from a diallel cross was revealed which could not possibly be interpreted easily per se, even with an already established genetic model. Instead, pollination scores were distilled into principal component scores described as Compatibility Components (CC1-CC3). These were then subjected to a conventional genome-wide association analysis. CC1 associated with markers on linkage groups (LGs) 1, 2, 3, and 6, CC2 exclusively with markers in a genomic region on LG 2, and CC3 with markers on LG 1. BLAST alignment with the Brachypodium physical map revealed highly significantly associated markers with peak associations with genes adjacent and four genes away from the chromosomal locations of candidate SI genes, S- and Z-DUF247, respectively. Further significant associations were found in a Brachypodium distachyon chromosome 3 region, having shared synteny with Lolium LG 1, suggesting further SI loci linked to S or extensive micro-re-arrangement of the genome between B. distachyon and L. perenne. Significant associations with gene sequences aligning with marker sequences on Lolium LGs 3 and 6 were also identified. We therefore demonstrate the power of a novel association genetics approach to identify the genes controlling multi-locus gametophytic SI systems and to identify novel loci potentially involved in already established SI systems.

In Silico Identification of Candidate Genes for Fertility Restoration in Cytoplasmic Male Sterile Perennial Ryegrass (Lolium perenne L.).

Perennial ryegrass (Lolium perenne L.) is widely used for forage production in both permanent and temporary grassland systems. To increase yields in perennial ryegrass, recent breeding efforts have been focused on strategies to more efficiently exploit heterosis by hybrid breeding. Cytoplasmic male sterility (CMS) is a widely applied mechanism to control pollination for commercial hybrid seed production and although CMS systems have been identified in perennial ryegrass, they are yet to be fully characterized. Here, we present a bioinformatics pipeline for efficient identification of candidate restorer of fertility (Rf) genes for CMS. From a high-quality draft of the perennial ryegrass genome, 373 pentatricopeptide repeat (PPR) genes were identified and classified, further identifying 25 restorer of fertility-like PPR (RFL) genes through a combination of DNA sequence clustering and comparison to known Rf genes. This extensive gene family was targeted as the majority of Rf genes in higher plants are RFL genes. These RFL genes were further investigated by phylogenetic analyses, identifying three groups of perennial ryegrass RFLs. These three groups likely represent genomic regions of active RFL generation and identify the probable location of perennial ryegrass PPR-Rf genes. This pipeline allows for the identification of candidate PPR-Rf genes from genomic sequence data and can be used in any plant species. Functional markers for PPR-Rf genes will facilitate map-based cloning of Rf genes and enable the use of CMS as an efficient tool to control pollination for hybrid crop production.

Orthology Guided Transcriptome Assembly of Italian Ryegrass and Meadow Fescue for Single-Nucleotide Polymorphism Discovery.

Single-nucleotide polymorphisms (SNPs) represent natural DNA sequence variation. They can be used for various applications including the construction of high-density genetic maps, analysis of genetic variability, genome-wide association studies, and map-based cloning. Here we report on transcriptome sequencing in the two forage grasses, meadow fescue ( Huds.) and Italian ryegrass ( Lam.), and identification of various classes of SNPs. Using the Orthology Guided Assembly (OGA) strategy, we assembled and annotated a total of 18,952 and 19,036 transcripts for Italian ryegrass and meadow fescue, respectively. In addition, we used transcriptome sequence data of perennial ryegrass ( L.) from a previous study to identify 16,613 transcripts shared across all three species. Large numbers of intraspecific SNPs were identified in all three species: 248,000 in meadow fescue, 715,000 in Italian ryegrass, and 529,000 in perennial ryegrass. Moreover, we identified almost 25,000 interspecific SNPs located in 5343 genes that can distinguish meadow fescue from Italian ryegrass and 15,000 SNPs located in 3976 genes that discriminate meadow fescue from both species. All identified SNPs were positioned in silico on the seven linkage groups (LGs) of using the GenomeZipper approach. With the identification and positioning of interspecific SNPs, our study provides a valuable resource for the grass research and breeding community and will enable detailed characterization of genomic composition and gene expression analysis in prospective × hybrids.

Soil compartment is a major determinant of the impact of simulated rainfall on desert microbiota.

Although desert soils support functionally important microbial communities that affect plant growth and influence many biogeochemical processes, the impact of future changes in precipitation patterns on the microbiota and their activities is largely unknown. We performed in-situ experiments to investigate the effect of simulated rainfall on bacterial communities associated with the widespread perennial shrub, Rhazya stricta in Arabian desert soils. The bacterial community composition was distinct between three different soil compartments: surface biological crust, root-attached, and the broader rhizosphere. Simulated rainfall had no significant effect on the overall bacterial community composition, but some population-level responses were observed, especially in soil crusts where Betaproteobacteria, Sphingobacteria, and Bacilli became more abundant. Bacterial biomass in the nutrient-rich crust increased three-fold one week after watering, whereas it did not change in the rhizosphere, despite its much higher water retention. These findings indicate that between rainfall events, desert-soil microbial communities enter into stasis, with limited species turnover, and reactivate rapidly and relatively uniformly when water becomes available. However, microbiota in the crust, which was relatively enriched in nutrients and organic matter, were primarily water-limited, compared with the rhizosphere microbiota that were co-limited by nutrients and water.