Association mapping with a diverse population of Puccinia graminis f. sp. tritici identified avirulence loci interacting with the barley Rpg1 stem rust resistance gene.

BACKGROUND: Wheat stem rust, caused by Puccinia graminis f. sp. tritici (Pgt), is an important disease of barley and wheat. A diverse sexual Pgt population from the Pacific Northwest (PNW) region of the US contains a high proportion of individuals with virulence on the barley stem rust resistance (R) gene, Rpg1. However, the evolutionary mechanisms of this virulence on Rpg1 are mysterious considering that Rpg1 had not been deployed in the region and the gene had remained remarkably durable in the Midwestern US and prairie provinces of Canada. METHODS AND RESULTS: To identify AvrRpg1 effectors, genome wide association studies (GWAS) were performed using 113 Pgt isolates collected from the PNW (n = 89 isolates) and Midwest (n = 24 isolates) regions of the US. Disease phenotype data were generated on two barley lines Morex and the Golden Promise transgenic (H228.2c) that carry the Rpg1 gene. Genotype data was generated by whole genome sequencing (WGS) of 96 isolates (PNW = 89 isolates and Midwest = 7 isolates) and RNA sequencing (RNAseq) data from 17 Midwestern isolates. Utilizing ~1.2 million SNPs generated from WGS and phenotype data (n = 96 isolates) on the transgenic line H228.2c, 53 marker trait associations (MTAs) were identified. Utilizing ~140 K common SNPs generated from combined analysis of WGS and RNAseq data, two significant MTAs were identified using the cv Morex phenotyping data. The 55 MTAs defined two distinct avirulence loci, on supercontig 2.30 and supercontig 2.11 of the Pgt reference genome of Pgt isolate CRL 75-36-700-3. The major avirulence locus designated AvrRpg1A was identified with the GWAS using both barley lines and was delimited to a 35 kb interval on supercontig 2.30 containing four candidate genes (PGTG_10878, PGTG_10884, PGTG_10885, and PGTG_10886). The minor avirulence locus designated AvrRpg1B identified with cv Morex contained a single candidate gene (PGTG_05433). AvrRpg1A haplotype analysis provided strong evidence that a dominant avirulence gene underlies the locus. CONCLUSIONS: The association analysis identified strong candidate AvrRpg1 genes. Further analysis to validate the AvrRpg1 genes will fill knowledge gaps in our understanding of rust effector biology and the evolution and mechanism/s of Pgt virulence on Rpg1.

Identification of Highly Repetitive Enhancers with Long-range Regulation Potential in Barley via STARR-seq.

Enhancers are DNA sequences that can strengthen transcription initiation. However, the global identification of plant enhancers is complicated due to uncertainty in the distance and orientation of enhancers, especially in species with large genomes. In this study, we performed self-transcribing active regulatory region sequencing (STARR-seq) for the first time to identify enhancers across the barley genome. A total of 7323 enhancers were successfully identified, and among 45 randomly selected enhancers, over 75% were effective as validated by a dual-luciferase reporter assay system in the lower epidermis of tobacco leaves. Interestingly, up to 53.5% of the barley enhancers were repetitive sequences, especially transposable elements (TEs), thus reinforcing the vital role of repetitive enhancers in gene expression. Both the common active mark H3K4me3 and repressive mark H3K27me3 were abundant among the barley STARR-seq enhancers. In addition, the functional range of barley STARR-seq enhancers seemed much broader than that of rice or maize and extended to ±100 kb of the gene body, and this finding was consistent with the high expression levels of genes in the genome. This study specifically depicts the unique features of barley enhancers and provides available barley enhancers for further utilization.

Promoter role of putrescine for molecular and biochemical processes under drought stress in barley.

Drought, which adversely affects plant growth and continuity of life and reduces product yield and quality, is one of the most common abiotic stresses at the globally. One of the polyamines that regulates plant development and reacts to abiotic stressors, including drought stress, is Putrescine (Put). This study compared the physiological and molecular effects of applying exogenous Put (10 µM) to barley (Hordeum vulgare cv. Burakbey) under drought stress (- 6.30 mPa PEG 6000). The 21-day drought stress imposed on the barley plant had a strong negative effect on plant metabolism in all experimental groups. Exogenous Put treatment under drought stress had a reformative effect on the cell cycle (transitions from G0-G1 to S and from S to G2-M), total protein content (almost 100%), endogenous polyamine content, malondialdehyde (MDA) (70%), and ascorbate peroxidase (APX) (62%) levels compared to the drought stress plants. Superoxide dismutase (SOD) (12%) and catalase (CAT) (32%) enzyme levels in the same group increased further after exogenous Put application, forming a response to drought stress. Consequently, it was discovered that the administration of exogenous Put in barley raises endogenous polyamine levels and then improves drought tolerance due to increased antioxidant capability, cell division stimulation, and total protein content.

EPAD1 Orthologs Play a Conserved Role in Pollen Exine Patterning.

The pollen wall protects pollen during dispersal and is critical for pollination recognition. In the Poaceae family, the pollen exine stereostructure exhibits a high degree of conservation with similar patterns across species. However, there remains controversy regarding the conservation of key factors involved in its formation among various Poaceae species. EPAD1, as a gene specific to the Poaceae family, and its orthologous genes play a conserved role in pollen wall formation in wheat and rice. However, they do not appear to have significant functions in maize. To further confirm the conserved function of EPAD1 in Poaceae, we performed an analysis on four EPAD1 orthologs from two distinct sub-clades within the Poaceae family. The two functional redundant barley EPAD1 genes (HvEPAD1 and HvEPAD2) from the BOP clade, along with the single copy of sorghum (SbEPAD1) and millet (SiEPAD1) from the PACMAD clade were examined. The CRISPR-Cas9-generated mutants all exhibited defects in pollen wall formation, consistent with previous findings on EPAD1 in rice and wheat. Interestingly, in barley, hvepad2 single mutant also showed apical spikelets abortion, aligning with a decreased expression level of HvEPAD1 and HvEPAD2 from the apical to the bottom of the spike. Our finding provides evidence that EPAD1 orthologs contribute to Poaceae specific pollen exine pattern formation via maintaining primexine integrity despite potential variations in copy numbers across different species.

The cap-binding complex modulates ABA-responsive transcript splicing during germination in barley (Hordeum vulgare).

To decipher the molecular bases governing seed germination, this study presents the pivotal role of the cap-binding complex (CBC), comprising CBP20 and CBP80, in modulating the inhibitory effects of abscisic acid (ABA) in barley. Using both single and double barley mutants in genes encoding the CBC, we revealed that the double mutant hvcbp20.ab/hvcbp80.b displays ABA insensitivity, in stark contrast to the hypersensitivity observed in single mutants during germination. Our comprehensive transcriptome and metabolome analysis not only identified significant alterations in gene expression and splicing patterns but also underscored the regulatory nexus among CBC, ABA, and brassinosteroid (BR) signaling pathways.

Unveiling the mysteries of HvANS: a study on anthocyanin biosynthesis in qingke (hordeum vulgare L. var. Nudum hook. f.) seeds.

BACKGROUND: Based on our previous research, a full-length cDNA sequence of HvANS gene was isolated from purple and white Qingke. The open reading frame (ORF) in the purple variety Nierumuzha was 1320 base pairs (bp), encoding 439 amino acids, while the ORF in the white variety Kunlun 10 was 1197 bp, encoding 398 amino acids. A nonsynonymous mutation was found at the position of 1195 bp (T/C) in the coding sequence (CDS) of the HvANS gene. We carried out a series of studies to further clarify the relationship between the HvANS gene and anthocyanin synthesis in Qingke. RESULTS: The conservative structural domain prediction results showed that the encoded protein belonged to the PLN03178 superfamily. Multiple comparisons showed that this protein had the highest homology with Hordeum vulgare, at 88.61%. The approximately 2000 bp promoter sequence of the HvANS gene was identical in both varieties. The real-time fluorescence PCR (qRT-PCR) results revealed that HvANS expression was either absent or very low in the roots, stems, leaves, and awns of Nierumuzha. In contrast, the HvANS expression was high in the seed coats and seeds of Nierumuzha. Likewise, in Kunlun 10, HvANS expression was either absent or very low, indicating a tissue-specific and variety-specific pattern for HvANS expression. The subcellular localization results indicated that HvANS was in the cell membrane. Metabolomic results indicated that the HvANS gene is closely related to the synthesis of three anthocyanin substances (Idaein chloride, Kinetin 9-riboside, and Cyanidin O-syringic acid). Yeast single hybridization experiments showed that the HvANS promoter interacted with HvANT1, which is the key anthocyanin regulatory protein. In a yeast two-hybrid experiment, we obtained two significantly different proteins (ZWY2020 and POMGNT2-like) and verified the results by qRT-PCR. CONCLUSIONS: These results provide a basis for further studies on the regulatory mechanism of HvANS in the synthesis of anthocyanins in Qingke purple grains.

Detection of closely linked QTLs and candidate genes controlling germination indices in response to drought and salinity stresses in barley.

The aim of current study was to identify closely linked QTLs and candidate genes related to germination indices under control, salinity and drought conditions in barley. A total of nine (a major), 28 (eight major) and 34 (five major) closely linked QTLs were mapped on the seven chromosomes in response to control, drought and salinity conditions using genome-wide composite interval mapping, respectively. The major QTLs can be used in marker-assisted selection (MAS) projects to increase tolerance to drought and salinity stresses during the germination. Overall, 422 unique candidate genes were associated with most major QTLs. Moreover, gene ontology analysis showed that candidate genes mostly involved in biological process related to signal transduction and response to stimulus in the pathway of resistance to drought and salinity stresses. Also, the protein-protein interaction network was identified 10 genes. Furthermore, 10 genes were associated with receptor-like kinase family. In addition, 16 transcription factors were detected. Three transcription factors including B3, bHLH, and FAR1 had the most encoding genes. Totally, 60 microRNAs were traced to regulate the target genes. Finally, the key genes are a suitable and reliable source for future studies to improve resistance to abiotic stress during the germination of barley.

Classification and prediction of drought and salinity stress tolerance in barley using GenPhenML.

Genetic and agronomic advances consistently lead to an annual increase in global barley yield. Since abiotic stresses (physical environmental factors that negatively affect plant growth) reduce barley yield, it is necessary to predict barley resistance. Artificial intelligence and machine learning (ML) models are new and powerful tools for predicting product resilience. Considering the research gap in the use of molecular markers in predicting abiotic stresses, this paper introduces a new approach called GenPhenML that combines molecular markers and phenotypic traits to predict the resistance of barley genotypes to drought and salinity stresses by ML models. GenPhenML uses feature selection algorithms to determine the most important molecular markers. It then identifies the best model that predicts atmospheric resistance with lower MAE, RMSE, and higher R2. The results showed that GenPhenML with a neural network model predicted the salinity stress resistance score with MAE, RMSE and R2 values of 0.1206, 0.0308 and 0.9995, respectively. Also, the NN model predicted drought stress scores with MAE, RMSE and R2 values of 0.0727, 0.0105 and 0.9999, respectively. The GenPhenML approach was also used to classify barley genotypes as resistant and stress-sensitive. The results showed that the accuracy, accuracy and F1 score of the proposed approach for salinity and drought stress classification were higher than 97%.

The CDPK-related protein kinase HvCRK2 and HvCRK4 interact with HvCML32 to negatively regulate drought tolerance in transgenic Arabidopsis thaliana.

Calcium-dependent protein kinase (CDPK) as one of calcium sensors were play important roles in stress responses. CDPK-related protein kinase (CRK) was identified as subgroup III of CDPK has been characterized in many plants, but the members and functions of CRK genes in hulless barley (Hordeum vulgare L.) has not been clarified. Here, we identified four HvCRK genes and named HvCRK1-4 according to chromosomes localization. Moreover, the physiological function of highly induced genes of HvCRK2 and HvCRK4 were investigated in drought stress tolerance by examining their overexpression transgenic lines functions generated in Arabidopsis thaliana. Under drought stress, both overexpression HvCRK2 and HvCRK4 displayed reduced drought resistance, and accompanied by higher accumulation levels of ROS. Notably, overexpression of HvCRK2 and HvCRK4 reduced sensitivity to exogenous ABA, meanwhile the expression of ABA-responsive genes in transgenic plants were down-regulated compared to the wild type in response to drought stress. Furthermore, the physically interaction of HvCRK2 and HvCRK4 with calmodulin (CaM) and calmodulin-like (CML) proteins were determined in vivo, the further results showed that HvCML32 binds to HvCRK2/4 S_TKC structural domains and negatively regulates drought tolerance. In summary, this study identified HvCRK members and indicated that HvCRK2 and HvCRK4 genes play negative roles in drought tolerance, and provide insight into potential molecular mechanism of HvCRK2 and HvCRK4 in response to drought stress.

Investigating the regulatory role of HvANT2 in anthocyanin biosynthesis through protein-motif interaction in Qingke.

Background: Currently, there are no reports on the HvbHLH gene family in the recent barley genome (Morex_V3). Furthermore, the structural genes related to anthocyanin synthesis that interact with HvANT2 have yet to be fully identified. Methods: In this study, a bioinformatics approach was used to systematically analyze the HvbHLH gene family. The expression of this gene family was analyzed through RNA sequencing (RNA-seq), and the gene with the most significant expression level, HvANT2, was analyzed using quantitative reverse transcription polymerase chain reaction (qRT-PCR) in different tissues of two differently colored varieties. Finally, structural genes related to anthocyanin synthesis and their interactions with HvANT2 were verified using a yeast one-hybrid (Y1H) assay. Results: The study identified 161 bHLH genes, designated as HvbHLH1 to HvbHLH161, from the most recent barley genome available. Evolutionary tree analysis categorized barley bHLH TFs into 21 subfamilies, demonstrating a pronounced similarity to rice and maize. Through RNA-Seq analysis of purple and white grain Qingke, we discovered a significant transcription factor (TF), HvANT2 (HvbHLH78), associated with anthocyanin biosynthesis. Subsequently, HvANT2 protein-motifs interaction assays revealed 41 interacting motifs, three of which were validated through Y1H experiments. These validated motifs were found in the promoter regions of key structural genes (CHI, F3'H, and GT) integral to the anthocyanin synthesis pathway. These findings provide substantial evidence for the pivotal role of HvANT2 TF in anthocyanin biosynthesis.

Identification and functional analysis of long non-coding RNA (lncRNA) and metabolites response to mowing in hulless barley (Hordeum vulgare L. var. nudum hook. f.).

BACKGROUND: Hulless barley (Hordeum vulgare L. var. nudum Hook. f.) is a significant cereal crop and a substantial source of forage for livestock. Long non-coding RNAs (lncRNAs) and metabolites play crucial roles in the nutrient accumulation and regeneration of hulless barley plants following mowing. The study aimed to identify differentially expressed lncRNAs and metabolites in hulless barley plants by analyzing transcriptomic and metabolomic datasets at 2 h, 24 h, and 72 h following mowing. RESULTS: The study revealed that 190, 90, and 438 lncRNA genes were differentially expressed at the 2 h, 24 h, and 72 h time points compared to the non-mowing control. We identified 14 lncRNA genes-11 downregulated and 3 upregulated-showing consistently significant differential expression across all time points after mowing. These differentially expressed lncRNAs target genes involved in critical processes such as cytokinin signaling, cell wall degradation, storage protein accumulation, and biomass increase. In addition, we identified ten differentially expressed metabolites targeting diverse metabolic pathways, including plant hormones, alkaloids, and flavonoids, before and after mowing at various time points. Endogenous hormone analysis revealed that cytokinin most likely played a crucial role in the regeneration of hulless barley after mowing. CONCLUSIONS: This study created a comprehensive dataset of lncRNAs, metabolites, and hormones in hulless barley after mowing, revealing valuable insights into the functional characteristics of lncRNAs, metabolites, and hormones in regulating plant regeneration. The results indicated that cytokinin plays a significant role in facilitating the regeneration process of hulless barley after mowing. This comprehensive dataset is an invaluable resource for better understanding the complex mechanisms that underlie plant regeneration, with significant implications for crop improvement.

High resolution mapping of a novel non-transgressive hybrid susceptibility locus in barley exploited by P. teres f. maculata.

Hybrid genotypes can provide significant yield gains over conventional inbred varieties due to heterosis or hybrid vigor. However, hybrids can also display unintended negative attributes or phenotypes such as extreme pathogen susceptibility. The necrotrophic pathogen Pyrenophora teres f. maculata (Ptm) causes spot form net blotch, which has caused significant yield losses to barley worldwide. Here, we report on a non-transgressive hybrid susceptibility locus in barley identified between the three parental lines CI5791, Tifang and Golden Promise that are resistant to Ptm isolate 13IM.3. However, F2 progeny from CI5791 × Tifang and CI5791 × Golden Promise crosses exhibited extreme susceptibility. The susceptible phenotype segregated in a ratio of 1 resistant:1 susceptible representing a genetic segregation ratio of 1 parental (res):2 heterozygous (sus):1 parental (res) suggesting a single hybrid susceptibility locus. Genetic mapping using a total of 715 CI5791 × Tifang F2 individuals (1430 recombinant gametes) and 149 targeted SNPs delimited the hybrid susceptibility locus designated Susceptibility to Pyrenophora teres 2 (Spt2) to an ~ 198 kb region on chromosome 5H of the Morex V3 reference assembly. This single locus was independently mapped with 83 CI5791 × Golden Promise F2 individuals (166 recombinant gametes) and 180 genome wide SNPs that colocalized to the same Spt2 locus. The CI5791 genome was sequenced using PacBio Continuous Long Read technology and comparative analysis between CI5791 and the publicly available Golden Promise genome assembly determined that the delimited region contained a single high confidence Spt2 candidate gene predicted to encode a pentatricopeptide repeat-containing protein.

Genome wide association study of seedling and adult plant leaf rust resistance in two subsets of barley genetic resources.

Leaf rust (LR) caused by Puccinia hordei is a serious disease of barley worldwide, causing significant yield losses and reduced grain quality. Discovery and incorporation of new sources of resistance from gene bank accessions into barley breeding programs is essential for the development of leaf rust resistant varieties. To identify Quantitative Trait Loci (QTL) conferring LR resistance in the two barley subsets, the Generation Challenge Program (GCP) reference set of 142 accessions and the leaf rust subset constructed using the Focused Identification of Germplasm Strategy (FIGS) of 76 barley accessions, were genotyped to conduct a genome-wide association study (GWAS). The results revealed a total of 59 QTL in the 218 accessions phenotyped against barley leaf rust at the seedling stage using two P. hordei isolates (ISO-SAT and ISO-MRC), and at the adult plant stage in four environments in Morocco. Out of these 59 QTL, 10 QTL were associated with the seedling resistance (SR) and 49 QTL were associated with the adult plant resistance (APR). Four QTL showed stable effects in at least two environments for APR, whereas two common QTL associated with SR and APR were detected on chromosomes 2H and 7H. Furthermore, 39 QTL identified in this study were potentially novel. Interestingly, the sequences of 27 SNP markers encoded the candidate genes (CGs) with predicted protein functions in plant disease resistance. These results will provide new perspectives on the diversity of leaf rust resistance loci for fine mapping, isolation of resistance genes, and for marker-assisted selection for the LR resistance in barley breeding programs worldwide.

Anther Culture Protocols for Barley and Wheat.

Doubled haploid (DH) techniques remain valuable tools for wheat (Triticum aestivum L.) and barley (Hordeum vulgare L.) genetic improvement, and DH populations are used extensively in breeding and research endeavors. Several techniques are available for DH production in wheat and barley. Here, we describe two simple, robust anther culture methods used to produce more than 15,000 DH wheat and barley lines annually in Australia.

Map-Based Cloning of the Causal Gene for a Seed Dormancy Quantitative Trait Locus in Barley.

Seed dormancy is an important agronomic trait in cereal crops. Throughout the domestication of cereals, seed dormancy has been reduced to obtain uniform germination. However, grain crops must retain moderate levels of seed dormancy to prevent problems such as preharvest sprouting in wheat (Triticum aestivum) and barley (Hordeum vulgare). To produce modern cultivars with the appropriate seed dormancy levels, it is important to identify the genes responsible for seed dormancy. With recent advances in sequencing technology, several causal genes for seed dormancy quantitative trait loci (QTLs) have been identified in barley and wheat. Here, we present a method to identify causal genes for seed dormancy QTLs in barley, a method that is also applicable to other cereals.

Targeted Modification of Grain Dormancy Genes in Barley.

Transgenesis technologies, such as overexpression or RNA interference-mediated suppression, have often been used to alter the activity of target genes. More recently developed targeted genome modification methods using customizable endonucleases allow for the regulation or knockout mutation of target genes without the necessity of integrating recombinant DNA. Such approaches make it possible to create novel alleles of target genes, thereby significantly contributing to crop improvement. Among these technologies, the Cas9 endonuclease-based method is widely applied to several crops, including barley (Hordeum vulgare). In this chapter, we describe an Agrobacterium-based approach to the targeted modification of grain dormancy genes in barley using RNA-guided Cas9 nuclease.

LncRNA cis- and trans-regulation provides new insight into drought stress responses in wild barley.

Drought is one of the most common abiotic stresses that affect barley productivity. Long noncoding RNA (lncRNA) has been reported to be widely involved in abiotic stress, however, its function in the drought stress response in wild barley remains unclear. In this study, RNA sequencing was performed to identify differentially expressed lncRNAs (DElncRNA) among two wild barley and two cultivated barley genotypes. Then, the cis-regulatory networks were according to the chromosome position and the expression level correction. The GO annotation indicates that these cis-target genes are mainly involved in "ion transport transporter activity" and "metal ion transport transporter activity". Through weighted gene co-expression network analysis (WGCNA), 10 drought-related modules were identified to contract trans-regulatory networks. The KEGG annotation demonstrated that these trans-target genes were enriched for photosynthetic physiology, brassinosteroid biosynthesis, and flavonoid metabolism. In addition, we constructed the lncRNA-mediated ceRNA regulatory network by predicting the microRNA response elements (MREs). Furthermore, the expressions of lncRNAs were verified by RT-qPCR. Functional verification of a candidate lncRNA, MSTRG.32128, demonstrated its positive role in drought response and root growth and development regulation. Hormone content analysis provided insights into the regulatory mechanisms of MSTRG.32128 in root development, revealing its involvement in auxin and ethylene signal transduction pathways. These findings advance our understanding of lncRNA-mediated regulatory mechanisms in barley under drought stress. Our results will provide new insights into the functions of lncRNAs in barley responding to drought stress.

Natural selection drives emergent genetic homogeneity in a century-scale experiment with barley.

Direct observation is central to our understanding of adaptation, but evolution is rarely documented in a large, multicellular organism for more than a few generations. In this study, we observed evolution across a century-scale competition experiment, barley composite cross II (CCII). CCII was founded in 1929 in Davis, California, with thousands of genotypes, but we found that natural selection has massively reduced genetic diversity, leading to a single lineage constituting most of the population by generation 50. Selection favored alleles originating from climates similar to that of Davis and targeted loci contributing to reproductive development, including the barley diversification loci Vrs1, HvCEN, Ppd-H1, and Vrn-H2. Our findings point to selection as the predominant force shaping genomic variation in one of the world's oldest biological experiments.

Growth Properties and Metabolomic Analysis Provide Insight into Drought Tolerance in Barley (Hordeum vulgare L.).

Drought stress is a major meteorological threat to crop growth and yield. Barley (Hordeum vulgare L.) is a vital cereal crop with strong drought tolerance worldwide. However, the underlying growth properties and metabolomic regulatory module of drought tolerance remains less known. Here, we investigated the plant height, spike length, effective tiller, biomass, average spikelets, 1000-grain weight, number of seeds per plant, grain weight per plant, ash content, protein content, starch content, cellulose content, and metabolomic regulation mechanisms of drought stress in barley. Our results revealed that the growth properties were different between ZDM5430 and IL-12 under drought stress at different growth stages. We found that a total of 12,235 metabolites were identified in two barley genotype root samples with drought treatment. More than 50% of these metabolites showed significant differences between the ZDM5430 and IL-12 roots. The Kyoto Encyclopedia of Genes and Genomes pathway analysis identified 368 differential metabolites mainly involved in starch and sucrose metabolism, the pentose phosphate pathway, pyrimidine metabolism, phenylalanine, tyrosine, and tryptophan biosynthesis in ZDM5430 under drought stress, whereas the different metabolites of IL-12 under drought stress related to starch and sucrose metabolism, the pentose phosphate pathway, 2-oxocarboxylic acid metabolism, cutin, suberine and wax biosynthesis, carbon metabolism, fatty acid biosynthesis, and C5-branched dibasic acid metabolism. These metabolites have application in the tricarboxylic cycle, the urea cycle, the met salvage pathway, amino acid metabolism, unsaturated fatty acid biosynthesis, phenolic metabolism, and glycolysis. On the other hand, the expression patterns of 13 genes related to the abovementioned bioprocesses in different barley genotypes roots were proposed. These findings afford an overview for the understanding of barley roots' metabolic changes in the drought defense mechanism by revealing the differently accumulated compounds.

Disentangling sources of gene tree discordance for Hordeum species via target-enriched sequencing assays.

Hordeum is an economically and evolutionarily important genus within the Triticeae tribe of the family Poaceae, and contains 33 widely distributed and diverse species which cytologically represent four subgenomes (H, Xa, Xu and I). These wild species (except Hordeum spontaneum, which is the primary gene pool of barley) are secondary or tertiary gene-pool germplasms for barley and wheat improvement, and uncovering their complicated evolutionary relationships would benefit for future breeding programs. Here, we developed a complexity-reduced pipeline via capturing genome-wide distributed fragments via two novel target-enriched assays (HorCap v1.0 and BarPlex v1.0) in conjugation with high-throughput sequencing of the enrichments. Both assays were tested for genotyping 40 species from three genera (Hordeum, Triticum, and Aegilops) containing 82 samples 67 accessions. Either of both assays worked efficiently in genotyping, while integration of both assays can significantly improve the robustness and resolution of the Hordeum phylogenetic trees. Interestingly, the incomplete lineage sorting (ILS) was inferred for the first time as the major factor causing phylogenetic discordance among the four subgenomes, whereas in New World species (carrying I genome) post-speciation introgression events were revealed. Through revising the evolutionary relationships of the Hordeum species based on an ancestral state reconstruction for the diploids and parental donor inference for the polyploids, our results raised new queries about the Hordeum phylogeny. Moreover, both newly-developed assays are applicable in genotyping and phylogenetic analysis of Hordeum and other Triticeae wild species.