Genome-Wide Association Study Reveals Marker Trait Associations (MTA) for Waterlogging-Triggered Adventitious Roots and Aerenchyma Formation in Barley

Waterlogging is an environmental stress, which severely affects barley growth and development. Limited availability of oxygen in the root zone negatively affects the metabolism of the whole plant. Adventitious roots (AR) and root cortical aerenchyma (RCA) formation are the most important adaptive traits that contribute to a plant’s ability to survive in waterlogged soil conditions. This study used a genome-wide association (GWAS) approach using 18,132 single nucleotide polymorphisms (SNPs) in a panel of 697 barley genotypes to reveal marker trait associations (MTA) conferring the above adaptive traits. Experiments were conducted over two consecutive years in tanks filled with soil and then validated in field experiments. GWAS analysis was conducted using general linear models (GLM), mixed linear models (MLM), and fixed and random model circulating probability unification models (FarmCPU model), with the FarmCPU showing to be the best suited model. Six and five significant (approximately −log10 (p) ≥ 5.5) MTA were identified for AR and RCA formation under waterlogged conditions, respectively. The highest −log10 (p) MTA for adventitious root and aerenchyma formation were approximately 9 and 8 on chromosome 2H and 4H, respectively. The combination of different MTA showed to be more effective in forming RCA and producing more AR under waterlogging stress. Genes from major facilitator superfamily (MFS) transporter and leucine-rich repeat (LRR) families for AR formation, and ethylene responsive factor (ERF) family genes and potassium transporter family genes for RCA formation were the potential candidate genes involved under waterlogging conditions. Several genotypes, which performed consistently well under different conditions, can be used in breeding programs to develop waterlogging-tolerant varieties.


Introduction
Waterlogging is one of the prominent abiotic stresses due to climate change and one of the major constraints for barley production worldwide. Waterlogging occurs when the soil is inadequately drained, or deposited water cannot seep into it efficiently, producing the existing water fraction in the soil surface layer 20% higher than the field water capacity [1,2]. Approximately 16% of the fertile areas around the earth undergo regular waterlogging caused by excessive rainfall, lack of soil drainage, and irregular topography, which result in severe economic loss [3].
In plants, the availability of oxygen for respiration is hampered in waterlogged organs because gas diffusion in water is about 10,000 times slower than that in the air [4,5]. To continue the energy supply, it is necessary for the waterlogged organs to switch over to anaerobic mode for energy production [6]. However, the energy trade-off via glycolysis and fermentation is inefficient when waterlogging is prolonged [4,7]. To survive from longterm waterlogging events, aerobic respiration for plants must be maintained via oxygen transport to enhance internal oxygen diffusion [8]. Root morphological and anatomical traits contribute to plant adaptation to waterlogged soils. One of such morphological traits is the formation of adventitious roots (AR), a trait that exists in many different plant species [9][10][11][12][13]. New adventitious roots facilitate water and nutrient uptake under hypoxia stress [14]. The formation of AR minimizes the distance for oxygen diffusion and improves gas diffusivity [15]. AR usually originate from the waterlogged part of hypocotyls or basal stem region, and such adaptation can replace the deteriorating primary roots [15]. Therefore, the adaptive responses of AR development to waterlogging might be more important than those of primary roots for survival in waterlogged soil [16].
Root cortical aerenchyma (RCA) formation in adventitious roots is another feature that enables plant adaptation to waterlogging conditions [17,18] by enhancing the internal diffusion of atmospheric and photosynthetic oxygen from the aerial parts to the waterlogged roots, allowing the roots to maintain aerobic respiration [19]. Aerenchyma formation increases the porosity of roots above the usual levels promoted by intercellular spaces [20]. The increase in root porosity of tolerant genotypes in response to waterlogging stress could represent the adaptation to hypoxic or anaerobic conditions [21]. RCA is a developmental process that is triggered by the hormone ethylene. Its activity contributes to the opening of gas spaces within parenchymatic tissues due to programed cell death (PCD) and cell wall modifications [16]. The formation of RCA has been studied extensively in the roots of various plant species, including barley [9], rice [22], maize [23], zea [11], wheat [21], and sugarcane [24].
Waterlogging also causes a sharp decrease in soil redox potential, resulting in very significant changes to the soil elemental profile, and significant accumulation of toxic substances occurs in soil [25]. Materials potentially toxic to plants that are accumulated in waterlogged soils include reduced manganese, iron, hydrogen sulphide, various organic acids, and ethylene [26]. AR and RCA formation are most directly related to a plant's superior ability to adapt to anoxic/toxic conditions [25]. Our previous study identified that those genotypes with higher proportion of AR and RCA formation provided substantial tolerance to prolonged waterlogging conditions and improved grain yield [18].
Genome-wide association studies (GWAS) have been applied in the trait-associated genetic studies with a large number of single nucleotide polymorphisms (SNPs) in many plant species, such as Arabidopsis thaliana [27], rice [28], and wheat [29,30]. Many QTL for waterlogging tolerance in barley have been detected [31][32][33], including those for aerenchyma formation [32,34], total root dry weight [35,36], root porosity [34], and ROS formation [37]. Most of these studies focused on only double haploid (DH) populations. GWAS is an alternative way for detecting QTL/MTA, which became an approach for unraveling the molecular genetic basis underlying the natural phenotypic variation [38]. GWAS methods have been developed to overcome many of the inherent drawbacks of the traditional QTL approaches due to a deep history of recombination events [39]. Through the recent developments in genome-wide studies, hundreds of accessions encompassing thousands of gene loci can be genotyped using high-throughput markers to enhance the efficiency of existing breeding approaches [40][41][42]. Thus, GWAS can specifically locate polymorphisms and the responsible genetic loci that are accountable for phenotypic variations to allow gene-targeted explorations [43,44].
The first GWAS for waterlogging tolerance in a wide-ranging barley collection under controlled field conditions has been conducted recently, targeting different agronomic traits and yield components [45], but did not investigate AR development and RCA formation. In this study, we used a much greater number of genotypes collected worldwide, which leads to more reproducible results [38], to (1) detect the novel MTA controlling natural phenotypic and genetic variation of AR and RCA formation in response to waterlogging; and (2) identify potential candidate genes for AR and RCA formation. Our main goal was to understand the key traits contributing to waterlogging tolerance in barley and to develop molecular markers that can be used in marker-assisted selection for accelerating a breeding program to enhance waterlogging tolerance in barley. Our secondary goal was to reproduce evidence for previous findings, using different growing conditions and a larger pool of genotypes.

Waterlogging Tolerance of Barley Accessions
The development of adventitious roots was quantified by allocating scores to plants (on a 1 to 10 scale) after sixty days waterlogging treatment when most of the sensitive lines showed prominent visual symptoms. The most tolerant varieties were scored 1 (a large number of AR and no dead roots), and the most sensitive/intolerant plants scored 10 (very few roots or all dead) ( Figure 1). Based on these root scores, waterlogging tolerance was divided into four levels: highly tolerant (score 1-3), tolerant (score above 3-5), sensitive (score above 5-7), and highly sensitive (score above 7-10). Based on the average values from two tank and one field experiments, the AR scores show that our barley population consists of approximately 7% highly tolerant, 31% moderately tolerant, 40% sensitive, and 20% highly sensitive accessions ( Figure 2A). For RCA formation, we also used a scoring number from 1 to 10, with the interpretations of scores shown in Figure 1; as for AR scores, lower values indicate higher tolerance to waterlogging. Among the RCA formation genotypes, about 14% genotypes were highly tolerant, 22% moderately tolerant, 27% sensitive, and 34% highly sensitive ( Figure 2B). The frequency distributions, means of all the traits, and their standard deviation are presented in Figure 2A,B. Based on the combination of adventitious root development and aerenchyma formation, genotypes are clustered into four groups: (a) tolerant lines with a large amount of AR and a high proportion of RCA formation (both scores close to 1), (b) tolerant lines with a high number of AR but poor ability for RCA formation, (c) intolerant lines with a lower number of AR and a high proportion of RCA, and (d) intolerant lines with limited AR and limited proportion of RCA (both scores approaching 10) (Table S1). Overall, AR development and RCA formation were relatively consistent, with the correlation being r = 0.58 (n = 691). The broad-sense heritability [46] estimates (as a measure of repeatability, or "intraclass correlation", of results for the same genotypes across experimental conditions) for AR development and RCA formation were 0.44 and 0.83, respectively.

SNP Distribution and Principal Components in Barley GWAS Population
SNPs are distributed all over the chromosomes ( Figure S1), with more in the tail end of the chromosome. The frequency of the SNPs in the middle area is less mainly due to cross-overs in centromere area of the chromosome. Frequency of the SNPs of most of the chromosomes ranged from 0 to 25 within one million base (Mb) window size. The first and second principal components account for only 8.1% and 5.1% of the variance in the genotype data matrix ( Figure S2).

Genome-Wide Association Study in Barley under Waterlogging Conditions
We first examined three different GWAS models: generalized linear models (GLM), mixed linear models with random effects (MLM), and FarmCPU models. The quantilequantile (QQ) plots from these different models (GLM, MLM, and FarmCPU) are shown in Figure 3A,C. The GLM model produced many −log 10 (p) values beyond what would be expected (i.e., there is an abundance of SNPs that appear relevant); the MLM method identified no significant SNPs; and the FarmCPU had a reasonably modest set of extremelog 10 (p) values. Accordingly, we focus on the SNPs identified from the FarmCPU modeling. Under control conditions, adventitious roots were not affected, and no aerenchyma was formed. However, under waterlogging conditions, adventitious roots were affected, and aerenchyma was formed. Waterlogging tolerance was divided into four levels: highly tolerant (score 1-3), tolerant (score above 3-5), sensitive (score above 5-7), and highly sensitive (score above 7-10).

MTA for Adventitious Root Formation
Using the FarmCPU model, six significant marker trait associations (MTA) were identified based on the average scores of adventitious roots (AR) formation ( Figure 3B). Table  S2 lists all significant MTA. These MTA are located on chromosomes 2H (SNP3268328), 5H (SNP3265960 and SNP5247903), 6H (SNP4790180 and SNP3266637), and 7H (SNP3271458) (Figure 4). The most significant MTA was on chromosome 2H with the observed −log 10 (p) values of about 9, with the tolerant allele of the SNP3268328 increasing the tolerance score by 0.72 (Table S2).

MTA for Aerenchyma Formation
Five significant MTA were identified for RCA formation in adventitious roots ( Figure 3D, Table S2). The MTA are on chromosomes 1H (SNP3663961), 2H (SNP3987361), 4H (SNP3433774), and 6H (SNP3913494 and SNP3272551) ( Figure 4). SNP3433774 on 4H showed the most significant association with RCA, with an observed −log 10 (p) value of about 8. The tolerant allele of SNP3433774 increased the average RCA score by 1.03 (Table S2). Two MTA on 2H were located at similar positions, thus counted as one MTA (SNP3987361) ( Figure 3D). This MTA is in a close position to the MTA on 2H for AR development. Similarly, the MTA on 6H for RCA formation were also located at near position to the MTA for AR development on 6H (Figure 4).

The Combined Effect of Different Alleles
Allelic effects of the two markers showing the most significant association (i.e., the most extreme p-values from the FarmCPU analysis) with average adventitious root development and aerenchyma formation scores are shown in Figure 5. Overall, tolerant alleles ("T" in the charts) showed lower visual scoring (healthier roots or greater proportion of RCA) than susceptible alleles. The data are for the top two SNPs (i.e., the most extreme p-values from the FarmCPU analysis) for AR development and RCA formation, respectively. The "SNP combo" data is from varieties, which have either both respective tolerant alleles, or both respective susceptible alleles. For AR development: n = 694 for SNP 3268328, n = 694 for 3271458, and n = 467 for the SNP combo. For RCA formation: n = 691 for SNP3433774, n = 691 for SNP3913494, and n = 466 for the SNP combo. The data are slightly jittered (at random) to expose the distributions, and the overlaying boxplots identify the interquartile range and the median.
The two most significant MTA for AR development showed additive effects, in the sense that the average score of the genotypes with both tolerant alleles was much lower than those with only one tolerant allele (as evidenced by the "SNP combo" data in Figure 5A). Combined allelic effects of MTA for aerenchyma formation were also significant in the same sense, and a greater number of the genotypes with both tolerant alleles showed a greater proportion of RCA formation (as evidenced by the "SNP combo" data in Figure 5B).
Simple t-tests comparing the means of varieties with the two most susceptible alleles (S) versus varieties with the two most tolerant alleles (T) were as follows. The "SNP combo" difference in means for adventitious root development scores (T minus S) was −1.38, with a 95% confidence interval of (−1.86, −0.89) and −log 10 p-value of 7.46. The "SNP combo" difference in means for aerenchyma formation scores (T minus S) was −2.30, with a 95% confidence interval of (−2.94, −1.65) and −log 10 p-value of 11.15. The rest of the test details are in Table S4.

Candidate Genes for Adventitious Root and Aerenchyma Formation
All six MTA for AR development and five MTA for RCA formation were considered as favorable MTA for candidate genes identification. A search for possible candidate genes for AR and RCA formation within the different regions of MTA was conducted. The putative domain family genes that were located within 2 Mb window from the MTA were considered as a possible candidate family gene. Waterlogging-responsive genes are located within the MTA region for AR, and RCA formations are listed in Table S3.
There are 28, 36, 12, 19, and 42 candidate genes that are located within 2 Mb screen of MTA.WL.AR.2H (SNP3268328), MTA.WL.AR.5H (SNP3265960 and SNP5247903), MTA.WL.AR.6H (SNP4790180 and SNP3266637), and MTA.WL.AR.7H (SNP3271458), respectively. Among all candidate genes, some of the gene families are highly induced by waterlogging or the SNPs of the MTA located inside the genes are considered decent candidate family genes. The two gene families that are most closely associated with the markers are MFS transporter superfamilies (HORVU5Hr1G093350 and HORVU5Hr1G093390) and leucine-rich repeat (LRR) family genes (HORVU5Hr1G093410 and HORVU7Hr1G107310).
The highly hypoxia-stress-responsive family genes or those closely associated with most significant MTA were considered for possible candidate family genes for RCA for-

Discussion
Barley waterlogging tolerance is influenced by many mechanisms, including RCA formation [17] and the development of AR [20]. In the past decade, many MTA/QTL have been mapped on all seven chromosomes of barley for waterlogging tolerance [31][32][33]37,47]. The value of MTA/QTL in breeding programs depends on their effects (and combined or additive effects) and whether they respond well in different backgrounds and environments. In this study, we identified some novel MTA. We also identified some other MTAs that are located at similar/near positions to those reported previously, which are associated with other traits. This element of "reproducibility" confirms the importance of the identified MTA/QTL [32]. In addition, the significant MTA for AR and RCA formation showed additive effects, thus showing the potential to be combined into more tolerant genotypes.

Formation of ARs and RCA under Waterlogging
Roots are the first organs to be affected by waterlogging, with the death of seminal roots being observed soon after hypoxia onset [48]. AR development is a typical responsive trait of waterlogged plants, which can partially replace the damaged seminal root and create more aerenchyma to improve the capability for internal oxygen transportation [49]. The formation of AR and RCA under waterlogging conditions has been reported for many crops, such as maize, soybean, tomato, and a deep-water rice [9,13,[50][51][52]. Our previous experiments identified that AR and RCA formation are some of the most important traits associated with waterlogging tolerance in barley [9,18]. The introgression of a single major RCA-QTL into commercial barley varieties can increase the waterlogging tolerance without significant negative effects on yield and quality [18].
Major QTL for waterlogging-triggered adventitious roots have been reported in maize [53,54] and cucumber [55]. Aerenchyma formation in adventitious roots has been suggested as a reliable method for screening waterlogging tolerance in maize and barley [16,32,54]. In this experiment, a total of 11 MTA was identified for AR and RCA formation under waterlogging stress using the 697 barley accessions collected worldwide. The most significant MTA (SNP3268328) for AR development was identified on 2H, and the most significant MTA (SNP3433774) for RCA formation was on chromosome 4H (Figure 4). MTA.AR (SNP3268328) on 2H was within 10 Mb from barley waterlogging tolerance QTL identified from the Yerong/Franklin and Franklin/TX9425 populations [31,33]. MTA.RCA (SNP3433774) on 4H was within 2-5 Mb distance from the major QTL for RCA formation and waterlogging tolerance reported in several DH populations [9,32,33,48].

Candidate Genes
AR and RCA formation induced by waterlogging is a complex process that is regulated by multiple hormone signaling pathways. Ethylene, as the primary signal in response to waterlogging stress, acts in both developmental modes of adventitious rooting, promotes auxin transport, and increases the sensitivity to auxin in the rooting zones, thereby facilitating adventitious root development [60]. The production of adventitious roots is directly related to ethylene production. Ethylene has been implicated in signaling cell death in the formation of aerenchyma in the cortex of adventitious roots of maize subjected to hypoxia [61].
Potassium also plays an important role as a determinant of cell fate, with cytosolic K + acting as a trigger of programed cell death under a range of biotic and abiotic stress conditions, and the plant's ability to take up K + is significantly reduced under O 2 -deficient conditions [25]. Under waterlogging conditions, the metabolism of gamma aminobutyric acid (GABA) can bind directly to the guard cell outward rectifying K + channels, thereby improving hypoxia tolerance [25,65]. Potassium transporter family genes (HORVU4Hr1G079150, 11 mem) and (HORVU7Hr1G107400, 32 mem) are positioned within 1 Mb from MTA.WL.RCA.4H and MTA.WL.AR.7H, respectively. High-affinity nitrate transporters play an important role in nitrogen acquisition under anaerobic conditions and are the members of the major facilitator superfamily (MFS) transporter genes [66]. MFS family genes (HORVU2Hr1G085160, 10 mem) and (HORVU2Hr1G085260, 5 mem) are located within 1 Mb areas of MTA.WL.RCA.2H (SNP3987361), and (HORVU5Hr1G093350, 20 mem) and (HORVU5Hr1G093390, 4 mem) are also within 2 Mb of MTA.WL.AR.5H (SNP5247903). Leucine-rich repeat receptor (LRR)-like kinases play a fundamental role in sensing external signals and regulating gene expression responses at the cellular level and induced by abiotic stresses [67]. The responses of this gene increased the oxidative stress tolerance in the germination and early root growth in Arabidopsis [68]. One LRR gene (HORVU5Hr1G093410), six LRR genes (HORVU6Hr1G033670) and eight LRR genes (HORVU7Hr1G107310) are found within 1 Mb MTA.WL.AR.5H (SNP5247903), MTA.WL.RCA.6H (SNP3913494), and MTA.WL.AR.7H (SNP3271458), respectively.
Under waterlogging conditions, aerobic respiration is restricted, and anaerobic respiration (glycolysis and fermentation) for energy generation is an adaptive strategy of plants [7]. We also found a substantial number of genes related to glycolysis and fermentation. FHY3/FAR1 family genes operate downstream of a metabolic pathway that converts the pyruvate derived from glycolysis into alanine. The FHY3/FAR1 family genes (HORVU5Hr1G092140, one mem) are located within 2 Mb of MTA.WL.AR.5H (SNP5247903). Due to a deficiency of oxygen under waterlogging conditions, plants adjust mitochondrial respiration to glycolysis by fermentation of pyruvate to ethanol via pyruvate decarboxylase (PDC) [69]. Mitochondrial pyruvate carrier family genes (HORVU2Hr1G085120, two mem) are positioned within 1 Mb of MTA.WL.RCA.2H (SNP3987361).
To conclude, we propose that ethylene responsive factor (ERF) family genes and potassium transporter family genes for RCA formation, and MFS transporter superfamily genes and leucine-rich repeat (LRR) family genes for AR development, could both be decent candidate genes that are involved in barley plant's response to waterlogging conditions. The other potential candidate genes that have been reported to be associated with waterlogging tolerance might also play an important role for AR and RCA formation in barley.

Plant Materials
A total of 697 barley accessions, collected worldwide, were used in this study. Our pool of barley genotypes includes approximately 30% that originate from China, 5% from commercial Australian varieties, 7% from wild barley of mixed global origin, and the balance primarily from Europe. The split of genotypes between two-row and six-row barley is approximately 70% and 30%, respectively.
DNA was extracted from leaf tissue collected at the 2-leaf seedling stage from a single plant per accession and genotyped with DArTSeq (http://www.diversityarrays.com/dartapplication-dartseq, accessed on 1 February 2022). Over 33,000 DArT and 31,000 SNP markers were scored. After removing those with the same scores or with greater distortion and/or greater proportion of missing data, a total of 18,132 markers with less than 5% missing data were used for marker-trait association study.

Growing Conditions, Treatment, and Phenotyping
Experiment one (tank): The 697 genotypes were grown in tanks under natural conditions at Tasmanian Institute of Agriculture in the 2020 and 2021 growing seasons ( Figures S3 and S4). Four seeds of each genotype were sown in a tank with 0.10 m row space and 0.06 m between genotypes with two replications.
Experiment two (field): A field experiment was performed in 2021 under field waterlogging conditions at the research station at Tasmanian Institute of Agriculture, Launceston, Tasmania ( Figures S3 and S4). Fifteen seeds of each from 697 genotypes were sown in 0.06 m rows with 0.25 m row space with two replications.
Waterlogging treatment was imposed at 2-3 leaf stage in all the above experiments and continued for two months. Root cortical aerenchyma formation in adventitious roots was checked 10 days after waterlogging. Aerenchyma formation in adventitious root was scored 1 to 10 ( Figure 1). Adventitious root scoring was also scored 1 to 10 after terminating two-month waterlogging treatment (Figure 1). Redox potential (Eh) was measured with redox probe, according to the method in our previous report [18]. Mean soil Eh value was significantly lower in the waterlogged treatments (−406 mV) than that in the well-drained control (+212 mV), indicating highly saturated and hypoxic conditions under waterlogging. Mean monthly rainfall, maximum and minimum temperature were comparable between the tank experiments and field experiment ( Figure S3).

Genome-Wide Association Study
A genome-wide association study (GWAS) was conducted using rMVP package of R software [76]. In the association study, the notional p-value threshold to suggest a significant marker trait association (MTA) was set to −log 10 (p) ≥ 5.5 (which is based on a singular test Type I error rate of 0.05, a simple Bonferroni adjustment [77], and 18,132 SNPs).

Model Used for GWAS Analysis
In this study, GLM, MLM, and fixed and random model circulating probability unification (FarmCPU, [78]) models were compared using adventitious root development and root cortical aerenchyma formation scoring from the 697 accessions and 18,132 SNPs. The best model was selected based on quantile-quantile (QQ) plots. The outcome variable in each model was the average of the respective AR or RCA values across the different trials.
The FarmCPU model was found to be the most suitable one, which had a small number of extreme p-values (see the QQ plots). The FarmCPU method iteratively performs marker tests with pseudo-quantitative trait nucleotides (QTNs) as covariates in a fixed-effect model and keeps optimizing pseudo-QTNs in a random-effect model, and this process continues until no new pseudo-QTNs are added [78]. The joint capabilities of removing confounding effects between testing markers and kinship, preventing model overfitting, and controlling false positives has made the FarmCPU method one of the most efficient methods in modern GWAS [79]. Within our implementation of the rMVP model, we used the following settings: default imputation, no principal components, "max Loop" was set to 10 (this exceeded what was required for our data), and the internal "FaST-LMM" method was used. Significant markers were visualized with a Manhattan plot, and important p-value distributions from all the GLM, MLM, and FarmCPU results (expected versus observed p-values on a −log 10 scale) are shown with quantile-quantile (QQ) plots.

Evaluation of Allelic Effect of Waterlogging Tolerance
The allele at the locus of the detected SNP responsible for increasing waterlogging tolerance (lower visual scoring (healthier roots or greater proportion of RCA)) is referred to as a "positive allele (T)" while the allele associated with increasing waterlogging susceptibility is classified as a "negative allele (S)" ( Figure 5). The relative "effects" of SNPs were determined by the difference between the scores of negative alleles and those of positive alleles (Table S2), based on the rMVP output for the FarmCPU model.

Candidate Gene Associated with Waterlogging Tolerance
SNPs at the level of −log 10 (p) ≥ 5.5 were considered highly significant and were used to identify candidate genes. We used the physical positions and sequence of significant SNPs to identify the annotation of the high-confidence (HC) candidate gene. Potential candidate genes around ±2 Mb areas of the significant SNPs were searched against BLAST search-Hordeum_vulgare-Ensembl Genomes 64 (gramene.org) with default parameters. The barley reference genome assembly (IBSC v2; [80]) was used to identify the position of the possible candidate genes. The sequences of potential candidate genes were also blasted on the database of Grain Genes Class Browser: Probe (usda.gov, accessed on 17 March 2022) searching for their functional annotation. Candidate genes with a possible connection to waterlogging tolerance were taken into consideration.

Conclusions and Prospective Research
Waterlogging seriously affects barley growth. GWAS analysis identified six and five MTA in barley for adventitious root development and aerenchyma formation, respectively, under waterlogging stress. The lead MTA that we identified can be combined with other significant MTA in breeding programs for increasing the proportion of RCA and improving AR development to elevate the barley waterlogging tolerance. Genes from the MFS transporter superfamily and leucine-rich repeat (LRR) family could be suitable targets for improving adventitious root development. The AP2/ERF domain superfamily genes and potassium transporter family genes might be good candidate genes for aerenchyma formation in adventitious root under waterlogging conditions. Several genotypes, which consistently produce well-developed adventitious roots and aerenchyma under different conditions can be used in breeding programs to develop waterlogging-tolerant varieties.
Further research should focus on the following aspects. An association mapping population can be developed with barley diverse germplasm lines. Further fine mapping and positional cloning of MTA/QTL approaches can be applied to isolate the genes of interest, since the trait-loci association can be mapped with more precision and with high resolution. Pyramiding of two traits via marker-assisted backcross strategy to introgress the favorable alleles into commercial cultivated varieties will help develop advanced waterlogging-tolerant cultivars.

Conflicts of Interest:
The authors declare no conflict of interest.