Characterization of Genetic Diversity and Genome-Wide Association Mapping of Three Agronomic Traits in Qingke Barley (Hordeum Vulgare L.) in the Qinghai-Tibet Plateau

Landscape genomics reveals adaptive genetic differentiation driven by multiple environmental variables in naked barley on the Qinghai-Tibetan Plateau

Understanding the local adaptation of crops has long been a concern of evolutionary biologists and molecular ecologists. Identifying the adaptive genetic variability in the genome is crucial not only to provide insights into the genetic mechanism of local adaptation but also to explore the adaptation potential of crops. This study aimed to identify the climatic drivers of naked barley landraces and putative adaptive loci driving local adaptation on the Qinghai-Tibetan Plateau (QTP). To this end, a total of 157 diverse naked barley accessions were genotyped using the genotyping-by-sequencing approach, which yielded 3123 high-quality SNPs for population structure analysis and partial redundancy analysis, and 37,636 SNPs for outlier analysis. The population structure analysis indicated that naked barley landraces could be divided into four groups. We found that the genomic diversity of naked barley landraces could be partly traced back to the geographical and environmental diversity of the landscape. In total, 136 signatures associated with temperature, precipitation, and ultraviolet radiation were identified, of which 13 had pleiotropic effects. We mapped 447 genes, including a known gene HvSs1. Some genes involved in cold stress and regulation of flowering time were detected near eight signatures. Taken together, these results highlight the existence of putative adaptive loci in naked barley on QTP and thus improve our current understanding of the genetic basis of local adaptation.

Genome diversity and highland-adaptative variation in Tibet barley landrace population of China

Barley landraces accumulated variation in adapting to extreme highland environments during long-term domestication in Tibet, but little is known about their population structure and genomic selection traces. In this study, tGBS (tunable genotyping by sequencing) sequencing, molecular marker and phenotypic analyses were conducted on 1,308 highland and 58 inland barley landraces in China. The accessions were divided into six sub-populations and clearly distinguished most six-rowed, naked barley accessions (Qingke in Tibet) from inland barley. Genome-wide differentiation was observed in all five sub-populations of Qingke and inland barley accessions. High genetic differentiation in the pericentric regions of chromosomes 2H and 3H contributed to formation of five types of Qingke. Ten haplotypes of the pericentric regions of 2H, 3H, 6H and 7H were further identified as associated with ecological diversification of these sub-populations. There was genetic exchange between eastern and western Qingke but they shared the same progenitor. The identification of 20 inland barley types indicated multiple origins of Qingke in Tibet. The distribution of the five types of Qingke corresponded to specific environments. Two predominant highland-adaptative variations were identified for low temperature tolerance and grain color. Our results provide new insights into the origin, genome differentiation, population structure and highland adaptation in highland barley which will benefit both germplasm enhancement and breeding of naked barley.

Single nucleotide polymorphism (SNP) markers for genetic diversity and population structure study in Ethiopian barley (Hordeum vulgare L.) germplasm

BACKGROUND

High-density single nucleotide polymorphisms (SNPs) are the most abundant and robust form of genetic variants and hence make highly favorable markers to determine the genetic diversity and relationship, enhancing the selection of breeding materials and the discovery of novel genes associated with economically important traits. In this study, a total of 105 barley genotypes were sampled from various agro-ecologies of Ethiopia and genotyped using 10 K single nucleotide polymorphism (SNP) markers. The refined dataset was used to assess genetic diversity and population structure.

RESULTS

The average gene diversity was 0.253, polymorphism information content (PIC) of 0.216, and minor allelic frequency (MAF) of 0.118 this revealed a high genetic variation in barley genotypes. The genetic differentiation also showed the existence of variations, ranging from 0.019 to 0.117, indicating moderate genetic differentiation between barley populations. Analysis of molecular variance (AMOVA) revealed that 46.43% and 52.85% of the total genetic variation occurred within the accessions and populations, respectively. The heat map, principal components and population structure analysis further confirm the presence of four distinct clusters.

CONCLUSIONS

This study confirmed that there is substantial genetic variation among the different barley genotypes. This information is useful in genomics, genetics and barley breeding.

The competitive strategies of poisonous weeds Elsholtzia densa Benth. on the Qinghai Tibet Plateau: Allelopathy and improving soil environment

Introduction

The competitive strategies of plants play a crucial role in their growth. Allelopathy is one of the weapons that plants use to improve their competitive advantage.

Methods

In order to explore the competitive strategy of a poisonous weed Elsholtzia densa Benth. (E. densa) on the Qinghai-Tibet Plateau (QTP), the effects of decomposing substances of E. densa on growth, root border cells (RBCs) characteristics of highland crop highland barley (Hordeum vulgare L.), and soil environment were determined.

Results

The decomposing allelopathic effect of E. densa on the germination and seedling growth of highland barley mainly occurred in the early stage of decomposing. The allelopathic effects were mainly on seed germination and root growth of highland barley. After treatment with its decomposing solution, the RBC's mucilage layer of highland barley thickened, and the RBC's activity decreased or even apoptosis compared with the control. However, only the above-ground part of the treatment group showed a significant difference. The effects of E. densa decomposed substances on the soil environment were evaluated from soil physicochemical properties and bacterial community. The results showed that soil bacteria varied greatly in the early stage of decomposion under different concentrations of E. densa. In addition, E. densa decomposing substances increased the soil nutrient content, extracellular enzyme activities, and bacterial community diversity. In the process of decomposition, the bacterial community structure changed constantly, but Actinobacteriota was always the dominant phylum.

Discussion

These results indicated that E. densa might adopt the following two strategies to help it gain an advantage in the competition: 1. Release allelochemicals that interfere with the defense function of surrounding plants and directly inhibit the growth and development of surrounding plants. 2. By changing the physical and chemical properties of soil and extracellular enzyme activity, residual plant decomposition can stimulate soil microbial activity, improve soil nutrition status, and create a more suitable soil environment for growth.

Evaluation of grain yield performance and its stability in various spring barley accessions under condition of different agroclimatic zones of Ukraine

Two extremely urgent problems of biological and agronomic research nowadays are ensuring an optimal balance between usage of natural resources to meet rapidly growing needs for food production and preservation of biodiversity. It is also important to extend the genetic diversity of the main crop varieties in agroecosystems. At the same time, modern varieties should be characterized by a combination of high yield and preserving yield stability under variable conditions. Solving the outlined tasks requires comprehensive research and involvement in breeding process of the genetical diversity concentrated in genebanks of the world. Barley (Hordeum vulgare L.) is one of the most important crops that satisfy the various needs of humanity. In respect to this, in 2020–2022, a multi-environment trial was conducted in three agroclimatic zones of Ukraine (Forest-Steppe, Polissia, and Northern Steppe). We studied 44 spring barley collection accessions of different ecological and geographical origin, different subspecies and groups of botanical varieties which were obtained from the National Center for Plant Genetic Resources of Ukraine. Statistical indices (Hom, Sc) and graphical models (GGE biplot, AMMI) were used to interpret the yield performance and its stability. Both individual ecological sites in different years and combinations of different sites and years of trials were characterized for productivity, discriminating power and representativeness. The environments differed quite strongly among themselves in terms of these indicators. It was established that most of the genotypes were characterized by higher adaptability to individual environmental conditions (stability in different years), compared to adaptability for all agroclimatic zones (wide adaptation). A strong cross-over genotype by environment interaction was found for most studied accessions. Nevertheless, both genotypes with very high stability in only one agroclimatic zone (Amil (UKR), Gateway (CAN)) and genotypes with a combination of high adaptability to one or two ecological niches and relatively higher wide adaptability (Stymul (UKR), Ly-1064 (UKR), Rannij (KAZ), Shedevr (UKR), and Arthur (CZE)) were identified. There were also the accessions which did not show maximum performance in the individual sites, but had relatively higher wide adaptability (Ly-1059 (UKR), Ly-1120 (UKR), Diantus (UKR), and Danielle (CZE)). In general, the naked barley genotypes were inferior to the covered ones in terms of yield potential and wide adaptability, but at the same time, some of them (CDC ExPlus (CAN), CDC Gainer (CAN), and Roseland (CAN)), accordingly to the statistical indicators, had increased stability in certain ecological sites. Among naked barley accessions relatively better wide adaptability according to the graphical analysis was found in the accession CDC McGwire (CAN), and by the statistical parameters CDC ExPlus (CAN) was better than standard. The peculiarities of yield manifestation and its variability in different spring barley genotypes in the multi-environment trial revealed in this study will contribute to the complementation and deepening of existing data in terms of the genotype by environment interaction. Our results can be used in further studies for developing spring barley variety models both with specific and wide adaptation under conditions of different agroclimatic zones of Ukraine. The disitnguished accessions of different origin and botanical affiliation are recommended for creating a new breeding material with the aim of simultaneously increasing yield potential and stability, as well as widening the genetic basis of spring barley varieties.

Genetic networks underlying salinity tolerance in wheat uncovered with genome-wide analyses and selective sweeps

A genetic framework underpinning salinity tolerance at reproductive stage was revealed by genome-wide SNP markers and major adaptability genes in synthetic-derived wheats, and trait-associated loci were used to predict phenotypes. Using wild relatives of crops to identify genes related to improved productivity and resilience to climate extremes is a prioritized area of crop genetic improvement. High salinity is a widespread crop production constraint, and development of salt-tolerant cultivars is a sustainable solution. We evaluated a panel of 294 wheat accessions comprising synthetic-derived wheat lines (SYN-DERs) and modern bread wheat advanced lines under control and high salinity conditions at two locations. The GWAS analysis revealed a quantitative genetic framework of more than 200 loci with minor effect underlying salinity tolerance at reproductive stage. The significant trait-associated SNPs were used to predict phenotypes using a GBLUP model, and the prediction accuracy (r²) ranged between 0.57 and 0.74. The r² values for flag leaf weight, days to flowering, biomass, and number of spikes per plant were all above 0.70, validating the phenotypic effects of the loci discovered in this study. Furthermore, the germplasm sets were compared to identify selection sweeps associated with salt tolerance loci in SYN-DERs. Six loci associated with salinity tolerance were found to be differentially selected in the SYN-DERs (12.4 Mb on chromosome (chr)1B, 7.1 Mb on chr2A, 11.2 Mb on chr2D, 200 Mb on chr3D, 600 Mb on chr6B, and 700.9 Mb on chr7B). A total of 228 reported markers and genes, including 17 well-characterized genes, were uncovered using GWAS and EigenGWAS. A linkage disequilibrium (LD) block on chr5A, including the Vrn-A1 gene at 575 Mb and its homeologs on chr5D, were strongly associated with multiple yield-related traits and flowering time under salinity stress conditions. The diversity panel was screened with more than 68 kompetitive allele-specific PCR (KASP) markers of functional genes in wheat, and the pleiotropic effects of superior alleles of Rht-1, TaGASR-A1, and TaCwi-A1 were revealed under salinity stress. To effectively utilize the extensive genetic information obtained from the GWAS analysis, a genetic interaction network was constructed to reveal correlations among the investigated traits. The genetic network data combined with GWAS, selective sweeps, and the functional gene survey provided a quantitative genetic framework for identifying differentially retained loci associated with salinity tolerance in wheat.

Genetic Diversity and Selection Signatures in Synthetic-Derived Wheats and Modern Spring Wheat

Synthetic hexaploid wheats and their derived advanced lines were subject to empirical selection in developing genetically superior cultivars. To investigate genetic diversity, patterns of nucleotide diversity, population structure, and selection signatures during wheat breeding, we tested 422 wheat accessions, including 145 synthetic-derived wheats, 128 spring wheat cultivars, and 149 advanced breeding lines from Pakistan. A total of 18,589 high-quality GBS-SNPs were identified that were distributed across the A (40%), B (49%), and D (11%) genomes. Values of population diversity parameters were estimated across chromosomes and genomes. Genome-wide average values of genetic diversity and polymorphic information content were estimated to be 0.30 and 0.25, respectively. Neighbor-joining (NJ) tree, principal component analysis (PCA), and kinship analyses revealed that synthetic-derived wheats and advanced breeding lines were genetically diverse. The 422 accessions were not separated into distinct groups by NJ analysis and confirmed using the PCA. This conclusion was validated with both relative kinship and Rogers' genetic distance analyses. EigenGWAS analysis revealed that 32 unique genome regions had undergone selection. We found that 50% of the selected regions were located in the B-genome, 29% in the D-genome, and 21% in the A-genome. Previously known functional genes or QTL were found within the selection regions associated with phenology-related traits such as vernalization, adaptability, disease resistance, and yield-related traits. The selection signatures identified in the present investigation will be useful for understanding the targets of modern wheat breeding in Pakistan.

Elucidation of gene action and combining ability for productive tillering in spring barley

The purpose of the present study is to identify breeding and genetic peculiarities for productive tillering in spring barley genotypes of different origin, purposes of usage and botanical affiliation, as well as to identify effective genetic sources to further improving of the trait. There were created two complete (6 × 6) diallel crossing schemes. Into the Scheme I elite Ukrainian (MIP Tytul and Avhur) and Western European (Datcha, Quench, Gladys, and Beatrix) malting spring barley varieties were involved. Scheme II included awnless covered barley varieties Kozyr and Vitrazh bred at the Plant Production Institute named after V. Y. Yuriev of NAAS of Ukraine, naked barley varieties Condor and CDC Rattan from Canada, as well as awned feed barley variety MIP Myroslav created at MIW and malting barley variety Sebastian from Denmark. For more reliable and informative characterization of barley varieties and their progeny for productive tillering in terms of inheritance, parameters of genetic variation and general combining ability (GCA) statistical analyses of experimental data from different (2019 and 2020) growing seasons were conducted. Accordingly to the indicator of phenotypic dominance all possible modes of inheritance were detected, except for negative dominance in the Scheme I in 2020. The degree of phenotypic dominance significantly varied depending on both varieties involved in crossing schemes and conditions of the years of trials. There was overdominance in loci in both schemes in both years. The other parameters of genetic variation showed significant differences in gene action for productive tillering between crossing Schemes. In Scheme I in both years the dominance was mainly unidirectional and due to dominant effects. In the Scheme II in both years there was multidirectional dominance. In Scheme I compliance with the additive-dominant system was revealed in 2019, but in 2020 there was a strong epistasis. In Scheme II in both years non-allelic interaction was identified. In general, the mode of gene action showed a very complex gene action for productive tillering in barley and a significant role of non-genetic factors in phenotypic manifestation of the trait. Despite this, the level of heritability in the narrow sense in both Schemes pointed to the possibility of the successful selection of individuals with genetically determined increased productive tillering in the splitting generations. In Scheme I the final selection for productive tillering will be more effective in later generations, when dominant alleles become homozygous. In Scheme II it is theoretically possible to select plants with high productive tillering on both recessive and dominant basis. In both schemes the non-allelic interaction should be taken into consideration. Spring barley varieties Beatrix, Datcha, MIP Myroslav and Kozyr can be used as effective genetic sources for involvement in crossings aimed at improving the productive tillering. The results of present study contribute to further development of studies devoted to evaluation of gene action for yield-related traits in spring barley, as well as identification of new genetic sources for plant improvement.

Genome wide association study of plant height and tiller number in hulless barley

Hulless barley (Hordeum vulgare L. var. nudum), also called naked barley, is a unique variety of cultivated barley. The genome-wide specific length amplified fragment sequencing (SLAF-seq) method is a rapid deep sequencing technology that is used for the selection and identification of genetic loci or markers. In this study, we collected 300 hulless barley accessions and used the SLAF-seq method to identify candidate genes involved in plant height (PH) and tiller number (TN). We obtained a total of 1407 M paired-end reads, and 228,227 SLAF tags were developed. After filtering using an integrity threshold of >0.8 and a minor allele frequency of >0.05, 14,504,892 single-nucleotide polymorphisms (SNP) loci were screened out. The remaining SNPs were used for the construction of a neighbour-joining phylogenetic tree, and the three subcluster members showed no obvious differentiation among regional varieties. We used a genome wide association study approach to identify 1006 and 113 SNPs associated with TN and PH, respectively. Based on best linear unbiased predictors (BLUP), 41 and 29 SNPs associated with TN and PH, respectively. Thus, several of genes, including Hd3a and CKX5, may be useful candidates for the future genetic breeding of hulless barley. Taken together, our results provide insight into the molecular mechanisms controlling barley architecture, which is important for breeding and yield.

Labelling Selective Sweeps Used in Durum Wheat Breeding from a Diverse and Structured Panel of Landraces and Cultivars

Simple Summary

Evaluation of the genetic diversity of a crop species is a critical step for breeding. Landraces are essential to avoid genetic erosion, and Mediterranean landraces are an important group of genetic resources due to their high genetic variability, adaptation to local conditions in rainfed environments, and their resilience to pests and pathogens. This study uses a genome-wide association approach employing eigenvectors to identify selective sweeps among Mediterranean durum wheat landraces and a world panel of modern cultivars.

Abstract

A panel of 387 durum wheat genotypes including Mediterranean landraces and modern cultivars was characterized with 46,161 diversity arrays technology (DArTseq) markers. Analysis of population structure uncovered the existence of five subpopulations (SP) related to the pattern of migration of durum wheat from the domestication area to the west of the Mediterranean basin (SPs 1, 2, and 3) and further improved germplasm (SPs 4 and 5). The total genetic diversity (H_T) was 0.40 with a genetic differentiation (G_ST) of 0.08 and a mean gene flow among SPs of 6.02. The lowest gene flow was detected between SP 1 (presumably the ancient genetic pool of the panel) and SPs 4 and 5. However, gene flow from SP 2 to modern cultivars was much higher. The highest gene flow was detected between SP 3 (western Mediterranean germplasm) and SP 5 (North American and European cultivars). A genome wide association study (GWAS) approach using the top ten eigenvectors as phenotypic data revealed the presence of 89 selective sweeps, represented as quantitative trait loci (QTL) hotspots, widely distributed across the durum wheat genome. A principal component analysis (PCoA) using 147 markers with −log₁₀p > 5 identified three regions located on chromosomes 2A, 2B and 3A as the main drivers for differentiation of Mediterranean landraces. Gene flow between SPs offers clues regarding the putative use of Mediterranean old durum germplasm by the breeding programs represented in the structure analysis. EigenGWAS identified selective sweeps among landraces and modern cultivars. The analysis of the corresponding genomic regions in the ‘Zavitan’, ‘Svevo’ and ‘Chinese Spring’ genomes discovered the presence of important functional genes including Ppd, Vrn, Rht, and gene models involved in important biological processes including LRR-RLK, MADS-box, NAC, and F-box.