Genome-wide association study of drought tolerance traits in sugar beet germplasms at the seedling stage

Abstract

Introduction:

Sugar beets are an important crop for global sugar production. Intense drought and the increasing lack of water resources pose a great threat to sugar beet cultivation. It is a priority to investigate favourable germplasms and functional genes to improve the breeding of drought tolerant plants.

Methods:

Thus, in this study, 328 sugar beet germplasms were used in a genome-wide association study (GWAS) to identify single nucleotide polymorphism (SNP) markers and candidate genes associated with drought tolerance.

Results:

The results showed that under drought stress (9% PEG-6000), there were 11 significantly associated loci on chromosomes 2, 3, 5, 7, and 9 from the 108946 SNPs filtered using a mixed linear model (MLM). Genome-wide association analysis combined with qRT-PCR identified 13 genes that were significantly differentially expressed in drought-tolerant extreme materials.

Discussion:

These candidate genes mainly exhibited functions such as regulating sugar metabolism, maintaining internal environmental stability and participating in photosystem repair. This study provides valuable information for exploring the molecular mechanisms of drought tolerance and improvement in sugar beet.

1 Introduction

Sugar beet is an important sugar crop. It has strong drought tolerance in the middle and late growth period, but weak drought tolerance at seedling stage. Most regions in the world are facing the problem of drought, and the sugar beet plant area of Heilongjiang Province in China is no exception. In spring, there are many southerly winds, large evaporation and small precipitation. The development of drought-tolerant varieties becomes more and more important (Bloch et al., 2006), and the selection and utilization of excellent drought tolerant germplasm resources and drought tolerant genes play an important role in the utilization efficiency and the improvement of agricultural production in water-shortage areas.

Mapping of quantitative trait loci (QTLs) is an effective tool often used to reveal the genetic basis of complex quantitative traits in crops (Huang et al., 2006). However, using traditional molecular markers such as restriction fragment length polymorphism (RFLP) and simple sequence repeats (SSR), only a few QTLs associated with drought stress can be identified (Jha et al., 2021; Dwiningsih et al., 2022). With the development of functional genomics and transcriptomics in recent years, a large number of genes have been found to be involved in plant drought tolerance, including protein kinases, transcription factors and some structural genes (Wang et al., 2021b). However, only a few genes have been functional validated. It remains a major challenge to accurately and effectively screen more drought stress-related genes and use them to breed drought-tolerant varieties.

Drought tolerance in plants is a complex quantitative trait that is controlled by multiple genes and involves multiple physiological and biochemical metabolic pathways. Although QTL mapping is a powerful method for detecting genomic regions associated with complex traits, the genetic effects of QTLs may not exist or may simply not be tested in different genetic backgrounds and environments (Shariatipour et al., 2021). Compared with traditional QTL, genome-wide association study (GWAS) can utilize genome-wide single nucleotide polymorphisms (SNPs) as molecular markers to analyze the genetic basis of complex traits, assisted by high-throughput genotyping platforms (Zheng et al., 2021). With the continuous development of high-throughput sequencing and high-resolution metabolic assays, as well as the development of multiple bioinformatics techniques and statistical methods to provide a basis for precise localization of complex trait gene variants, genome-wide association analysis is playing an increasingly important role in mining genetic loci for research (Zheng et al., 2021). Significant markers associated with aboveground biomass, fruit size, lodging score and leaf elongation in beans were identified through a whole genome-wide association analysis of 96 common bean germplasms. Specifically, 7 significant markers were found under irrigated conditions, and 5 significant markers were found under water deficit conditions (Hoyos-Villegas et al., 2017). In 2019, Mathew et al. used genome-wide association analysis to identify 75 significant marker-trait associations for population structure and marker traits of wheat biomass traits under drought and non-stress conditions. A total of 37 presumptive candidate genes were screened through gene annotation on IWGSC RefSeq 1.1 (Mathew et al., 2019). In recent years, the genome-wide association method has been applied in many crop drought tolerance studies and achieved good results. However, the current research on drought tolerance of sugar beet mainly focus on physiological and biochemical responses. Although traditional breeding methods such as analysis of combinin ability of lines have been used for selection of superior sugar beet genotypes for sugar related traits (Hassani et al., 2020), GWAS provide more reliable information about markers linked with sugar trait for use in marker-assisted selection We believe that using genome-wide association analysis in drought tolerance studies of beet germplasm resources can quickly and accurately excavate drought tolerance related genes. Therefore, in this study, 328 beet germplasm resources were simulated by PEG-6000 drought stress. It is hoped to explore significant loci and potential candidate genes related to drought tolerance of beet germplasm resources by genome-wide association analysis, and to provide reference and basis for further research on the isolation of related genes and molecular marker-assisted selection of beet drought tolerance.

2 Materials and methods

2.1 Test material

The material for this study was provided by the National Beet Medium-Term Gene Bank at Heilongjiang University. The 328 sugar beet germplasm resources used in the trial were from 17 countries (Supplementary Table S1). 200 from China, 9 from Japan, 1 from North Korea, 15 from Russia, 49 from the United States, 1 from Austria, 2 from Belgium, 10 from Poland, 1 from Denmark, 10 from Germany, 3 from France, 12 from the Netherlands, 2 from Romania, 5 from Sweden, 2 from Hungary, 4 from Italy, and 2 from the United Kingdom (Figure 1). This experiment was conducted in a sophisticated climate chamber with a constant environment of 25°C and 60% humidity during the day and 18°C and 60% humidity during the night. There were 4 rows for each genotype, and 6 plants in each row were cultured in Hoagland’s solution. After growing one pair of true leaves, two of the four rows of germplasms were selected and treated with PEG-6000 (Livak and Schmittgen) at a concentration of 9% to simulate drought stress (Tan et al., 2023). The remaining two rows of germplasms continued to be cultured in Hoagland’s solution as controls (CG). The solution was changed once every 3 days. The materials were sampled after the growth of three pairs of true leaves (Cai et al., 2015).

FIGURE 1

2.2 Phenotype identification and physiological trait assay

The following traits were measured in this test: embryonic axis diameter (stem diameter measured using Vernier callipers) (Yang et al., 2022b), maximum root length (a straightedge was used to measure the maximum length of the root), plant height (measured from the ground to the top of the plant using a straightedge), root fresh weight (the roots were removed, dried on paper towels, and weighed on an electronic balance), leaf fresh weight (all leaves were cut off and weighed on an electronic balance), root dry weight (roots were dried in a ventilated oven at 105°C for 15 min, then transferred to 75°C for drying to a standard weight) (Polash et al., 2018), and leaf dry weight (leaves were dried in a ventilated oven at 105°C for 15 min, then transferred to 75°C for drying to a standard weight) (Yang et al., 2022a). Related tools were Vernier callipers (0.01 mm, Brand: Sheffield, Model: S071012, Origin: Mainland China), electronic balance (0.00001 g, Brand: Mettler-Toledo GmbH, Model: MS105DU/A, Origin: Switzerland), and steel ruler (0.1 cm, Brand: Sheffield, Model: S079020, Origin: China Mainland). The relative leaf water was calculated as follows:

Relative leaf water = (Leaf fresh weight-Leaf dry weight)/(Leaf saturated fresh weight-Leaf dry weight).

Superoxide dismutase activity was evaluated by the nitrogen blue tetrazolium method (Spitz and Oberley, 1989). Soluble sugar levels were evaluated by the Anthrone-colorimetric method (Maness, 2010). Proline levels were evaluated by the acidic ninhydrin method (Carillo and Gibon, 2011). Soluble protein levels were evaluated using a BCA Protein Assay Kit (Brand: biosharp, Model: BL521A-3, Origin: Beijing, China).

The drought tolerance coefficient (DTC) was calculated for each trait:where DTC_X is the drought tolerance coefficient of each trait; DT is the mean value of the traits under drought treatment; and NT is the mean value of the traits under normal treatment.

2.3 Genotype data

Sequencing: SNPs were identified by NOVOGENE (Beijing, China) and sequenced using the Illumina HiSeq sequencing platform with double-end (Paired-End) 150 sequencing. BWA (Li and Durbin, 2009) (0.7.17): Interference information was filtered out of the raw data to obtain high-quality clean data, and effective high-quality sequencing data were aligned to the reference genome by BWA software (parameter: mem-t 4-k 32-M) (reference genome download address: ftp://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/000/511/025/GCF_000511025.2_RefBeet-1.2.2/GCF_000511025.2_RefBeet-1.2.2_genomic.fna.gz). Samtools (Soskic et al., 2019) (1.9): Samtools was used to transform the format of the sam file and build an index to generate the bai file. GATK (Sun et al., 2019) (4.2.6.1): GATK was used to detect mutations and generate group modification.vcf file. VCFtools (Ratnapriya et al., 2019) (0.1.17): A total of 108946 high-quality loci were screened for subsequent analysis using sequencing depth (dp ≥ 4), minor allele frequency (maf ≥ 0.05) and miss rate (miss ≤ 80%).

2.4 Population structure and LD analysis

VCFtools (0.1.17) was first used to convert the.vcf files ped files and then plink (1.9) software was used to convert ped format to bed format (Cadzow et al., 2014). Population structure analysis of bed files was performed using admixture (1.3.0) software, the population size K value was preset to 1–10 for classification, and the minimum value of cross entropy was used to determine the optimal number of classifications (Byrne et al., 2020). Principal component analysis of 108946 high-quality SNP loci was performed using Tassel (5.2.82) software (Breria et al., 2020), and the PCA results were visualized using R (4.1.0) software. The LD of 328 sugar beet germplasms was analysed by PopLDdecay (3.4.2). The LD intensity (r2) between two SNPs within a certain distance was calculated (Xu et al., 2018).

2.5 Genome-wide association study

The kinship analysis was performed using Tassel (5.2.82) software to obtain the kinship matrix (Zeng et al., 2017) and the results were visualized using the pheatmap package in R (4.1.0). The mixed linear model (MLM) (Kang et al., 2015) of Tassel software was used to analyse the PCA and kinship analysis results as covariates. This combination (Q + K) is often considered the most powerful, and this model can classify crops into different subgroups based on growth habit or geographical origin. Mixed models corrected with population structure and kinship can effectively reduce false positive results (Wang et al., 2007). Manhattan and Q-Q plots were visualized using the CMplot package of R (4.0.1) (Kumar et al., 2021). The mixed linear model (MLM) equation was as follows:

In the equation, Y is the vector of measured phenotypes, X is the fixed effect of SNPs, P is the fixed effect of population structure, K is the random effect of kinship, and e is the random error (Dramadri et al., 2021).

We used the powerful annotation function of the internationally recognized UniProt database (Yang and Wang, 2015) to annotate the SNPs. After screening fake SNPs according to genome annotation, we obtained high-quality and significant SNPs with a threshold of −log10p> 6.

The SNPs were screened to obtain high-quality and significantly correlated SNPs with a threshold of −log10p> 6. The size of the screened SNP interval was estimated based on LD block analysis, and genes in the region near the locus were obtained and functionally annotated. Based on the functional annotations, we could roughly determine which genes might be associated with the traits of interest and use in subsequent analysis (Bai et al., 2022).

2.6 Real-time fluorescence quantitative PCR

Four drought-tolerant (Mono HY53, Neitang578, Gansunongjiazhong, BGRC16137) and four drought-sensitive germplasms [7501A/BCE, 92017/1-8/1, Tianyansanhao (X), 92017/1-4/1] were used as extreme materials to extract total RNA according to the Trizol method. cDNA was synthesized according to the TranScript One-Step gDNA Removal and cDNA Synthesis SuperMix protocol (AT311, Beijing All Style Gold Biotechnology Co., Ltd.). qPCR was performed using the SuperReal PreMix Plus kit (Beijing Tiangen Biochemical Technology Co., Ltd., version FP210831) on a real-time fluorescent quantitative PCR instrument (Thermo Fisher Scientific Instruments, Shanghai. QuantStudio™ 1 Plus) (Liu et al., 2023). Expression level of 14 genes was analyzed with the 2^−ΔΔCT method (Livak and Schmittgen, 2001). BvGAPDH (NC_ 024800) was used as an internal control to standardize the expression levels of the different samples. All assays were performed in two independent experiments and replicated three times. Primers were summarized in Supplementary Table S5.

3 Results

3.1 Phenotypic and physiological analysis

The statistical results showed that the measured data were all close to a normal distribution, indicating that the 12 traits were all typical quantitative traits. Eight traits, including embryonic axis diameter (EAD), plant height (PH), root length (RL), leaf fresh weight (LFW), root fresh weight (RFW), root dry weight (RDW), leaf dry weight (Lander et al., 2001) and leaf relative water content (RLW), were significantly lower in the drought-treated plants than in the control plants. Three indicators, including soluble sugar levels (SS), soluble protein levels (SP) and free proline levels (Pro), were significantly higher in the drought-treated plants than in the control plants (Figure 2). Drought tolerance coefficients for each trait showed that these traits all exhibited suitable variation, indicating that they were controlled by small effect polygenes with coefficients of variation (CV) ranging from 16.29% to 79.58% (Table 1). The CV is an absolute value reflecting the degree of dispersion of each trait in the population. To better estimate the correlations between markers and traits, it is desirable to use highly dispersed trait data for statistical analysis. The traits in this test exhibited a high degree of dispersion and were suitable for genome-wide association analysis.

FIGURE 2

3.2 Population structure and genetic diversity

Population stratification refers to the existence of subgroups within a population, where the diversity among individuals within the subgroups is greater than the average diversity among individuals within the entire population (Wang et al., 2021a). Different allele frequencies at certain loci between different subpopulations can lead to false positive results when using two subpopulations for association analysis (Zheng et al., 2021). Therefore, before performing association analysis, we used a subset of 108946 SNPs without close linkage for population structure prediction. Using Tassel (5.2.82) software to construct phylogenetic trees by the neighbor-joining (NJ) method, 328 sugar beet germplasms were divided into three subgroups, but the differences between the subgroups were small (Figures 3A, D). The results from a principal component analysis (PCA) of the above three subgroups were consistent with the phylogenetic analysis, and no significant grouping was found (Figure 3E). The results indicated that the 328 germplasms without significant grouping were suitable for a GWAS. Using admixture software, the number of subgroups (K) was set to 1–10, and each K value was repeated three times to select the optimal K value according to the maximum likelihood method. The cross-validation error (CV error) value was the smallest when K = 7 (Figure 3B), so the Q-matrix with K = 7 was used as a covariate to improve the accuracy of the GWAS (Figure 3C).

FIGURE 3

A linkage disequilibrium (LD) analysis was carried out, and the LD distance decreased as the physical location of the SNPs on chromosomes gradually increased. The LD distance of 50 kb was equal to half of the maximum value (Figure 4).

FIGURE 4

The visualization of the kinship matrix showed that the variation between most of the 328 genotypes was low (light blue), and the genetic similarity between individual genotypes was high. The low genetic variation in this study would reduce the occurrence of false positive results (Figure 5).

FIGURE 5

3.3 GWAS of drought tolerance coefficients for each trait

The GWAS results showed (Figure 6; Supplementary Table S2) that a total of 183 SNPs (−log10p > 4.5) were identified in the 12 drought tolerance-related traits, including 11 SNPs with peaks >6.0, which the study mainly focused on. It was found that 1 SNP on chromosome 7 was associated with DTC_RLW; 2 SNPs associated with DTC_LFW were located on chromosomes 2 and 5; 1 SNP on chromosome 3 was associated with DTC_LDW; 2 SNPs on chromosome 9 were associated with DTC_RDW; and there were 5 SNPs associated with DTC_SS on chromosomes 2 and 7. It is worth noting that the two significant loci associated with DTC_RDW were close to each other with a distance of only 216 bp.

FIGURE 6

3.4 Candidate genes associated with sugar beet drought tolerance

In the five traits (−log10p> 6) of DTC_RLW, DTC_LFW, DTC_LDW, DTC_RDW, and DTC_SS, 50 kb upstream and downstream of the significant loci were analyzed by LD decay. A total of 24 genes in sugar beet associated with drought tolerance were identified and annotated (Supplementary Table S3).

Two genes (BVRB_7g160020 and BVRB_7g160030) were identified on chromosome 7 at position 4,745,224 may be associated with DTC_RLW (Figure 7; Supplementary Table S3). BVRB_7g160020 encodes a tyrosine-protein phosphatase, RLPH2, which is reported to be a new phospho-tyrosine-specific phosphatase belonging to the phosphoprotein phosphatase (PPP) family (Uhrig et al., 2016). BVRB_7g160030, as 1,4-alpha-glucan-branching enzyme 2-2, is mainly involved in amylose synthesis in gluconeogenesis (Han et al., 2022). It was reported that BVRB_7g160030 catalytically regulates α (1–6) glycosidic bond branching synthesis. We believe that BVRB_7g160030 is a potential target gene for the drought stress response in sugar beet.

FIGURE 7

Associated with DTC_LFW, three genes (BVRB_2g033710, BVRB_2g033720 and BVRB_2g033730) were found on chromosome 2. Two genes (BVRB_5g114800 and BVRB_5g114810) were found on chromosome 5. BVRB_2g033720, BVRB_2g033730 and BVRB_5g114800 were of unknown function (Figure 8; Supplementary Table S3). Through annotation, we found that BVRB_2g033710 functions as myosin-15; BVRB_5g114810 functions as MYB-related protein 306, the MYB structural domain is a peptide segment of approximately 51-52 amino acids containing a series of highly conserved amino acid residues and spacer sequences, and the MYB gene family is widely involved in plant metabolic regulation (Jiao et al., 2021; Johnson et al., 2022).

FIGURE 8

Two genes (BVRB_3g048910 and BVRB_3g048900) were identified on chromosome 3 position: 998,024, may be associated with DTC_LDW (Figure 9; Supplementary Table S3). By annotation, we found that BVRB_3g048910 functions as the probable methyltransferase PMT3. Methyltransferases use S-adenosylmethionine, betaine and dimethylthetin as methyl donors to generate methionine. BVRB_3g048900 functions as serine/threonine-protein kinase tricornered (TOR), which functions as a phosphate donor using ATP to phosphorylate serine and threonine residues on proteins to regulate abscisic acid, growth hormone, glucose and sucrose-mediated signalling (Ahn et al., 2011; Xiong et al., 2013; Deng et al., 2016).

FIGURE 9

Two close significant loci (position: 1,172,900, position: 1,173,116) on chromosome 9 were found to be associated with DTC_RDW (Figure 10; Supplementary Table S3). Six gene were found in the region near these two significant loci (BVRB_9g203050, BVRB_9g203040, BVRB_9g203030, BVRB_9g203020, BVRB_9g203010 and BVRB_9g20300). BVRB_9g203050 functions as photosystem II (PSII) stability/assembly factor HCF136 and is essential for PSII biogenesis. BVRB_9g203030 functions as eukaryotic initiation factor 4A-9, which promotes the hydrolysis of ATP and drives RNA deconvolution against abiotic stress; BVRB_9g203000 functions as enhancer of AG-4 protein 2, which possesses the dual ability to induce apoptosis and autophagy (Kumar et al., 2015), and is a transcription factor that acts as a flowering repressor by enhancing the expression of genes that delay flowering and suppressing nutritional growth to inhibit reproductive development (Wang et al., 2007).

FIGURE 10

Two gene regions (BVRB_2g036960 and BVRB_2g036950) were found to be associated with DTC_SS (Figure 11; Supplementary Table S3) in the region near positions 20,553,998 on chromosome 2. BVRB_2g036950 is an uncharacterized protein, and BVRB_2g036960 functions as a probable membrane-associated kinase regulator 4. Tyrosine phosphorylation controls the activation of brassinosteroid receptors through membrane-released protein kinase inhibitors.

FIGURE 11

Four SNP loci significantly associated with DTC_SS (Figure 12; Supplementary Table S3) were found on chromosome 7. No genes were found near position 15,215,427. Two genes (BVRB_7g163640 and BVRB_7g163650) were found near position 15,170,171 and they are hypothetical and uncharacterized proteins. Three genes (BVRB_7g166470, BVRB_7g166480 and BVRB_7g166490) were found near position 27,948,364, and BVRB_7g166470 functions as glycine-rich cell wall structural protein 1, a process that results in the assembly, arrangement of constituent parts, or disassembly of the cell wall. BVRB_7g166480 functions as PROTEIN LOW PSII ACCUMULATION 1, which is mainly involved in the efficient assembly of photosynthetic system II. BVRB_7g166490 functions as photosystem I reaction centre subunit IV. Two gene (BVRB_7g166670 and BVRB_7g166680) were found near positions 28,950 and BVRB_7g166670 is an uncharacterized protein, and BVRB_7g166680 functions as a ras GTPase-activating protein-binding protein.

FIGURE 12

3.5 qRT-PCR validation of candidate gene expression

By examining the overlap of the identified candidate genes with known drought stress response genes in other crop species, we identified 14 genes that may be involved in the drought tolerance process in sugar beet. In order to validate the GWAS results, qRT-PCR was performed on these genes (As shown in Figure 13). The results showed that, except for BVRB_2g036960, the other 13 genes were all significantly differently expressed between the drought tolerant and sensitive groups. The expression levels of BVRB_7g160020 and BVRB_5g114810 were 3.30-fold and 2.57-fold higher in the drought-sensitive group than them in the drought-tolerant group, respectively. BVRB_7g160030, BVRB_2g033710, BVRB_3g048910, BVRB_3g048900, BVRB_9g203050, BVRB_9g203030, BVRB_9g203000, BVRB_7g166680, BVRB_7g166470 and BVRB_7g166490 in drought-tolerant group were expressed 1.47-4.23-fold of them in the drought-sensitive group. Notably, the expression level of BVRB_7g166480 in the drought-tolerant group was 46.69-fold higher than that of the drought-sensitive group.

FIGURE 13

4 Discussion

4.1 Drought tolerance-related traits in sugar beet

Drought stress is the main abiotic stress factor that causes crop yield reduction worldwide. The sugar beet plant is more drought tolerant in the middle and late growth stages but is drought sensitive in the seedling stage, which causes great taproot loss (Fugate et al., 2018). Enhancing the drought tolerance of seedlings will greatly help to improve the sugar beet crop yield. Sadeghian and Yavari conducted mannitol-simulated drought stress on 9 sugar beet varieties, assessed the dry weight of cotyledons, fresh weight of cotyledons, RFW, RDW, and RL, and concluded that the drought tolerance traits of sugar beet seedlings can be used for selective breeding of the strain with the strongest tolerance (Sadeghian and Yavari, 2004). Islam et al. (2020) screened 7 beet varieties suitable for drought conditions by carrying out drought stress experiments on 11 varieties and measuring related physiological and biochemical traits. Through stepwise regression analysis, Moosavi et al. (2017) found that root and petiole dry weight were suitable as standard traits for drought tolerance of sugar beet, and 7221 drought tolerance genotypes were screened by these traits.

Although breeding by screening drought-tolerant varieties of sugar beet can improve drought tolerance in new varieties, this method has some limitations. The transfer of agronomic traits is easily limited by intermediate reproductive isolation, which is not conducive to genetic improvement of selected crops using genetic resources of closely or distantly related species. This study aims to improve the efficiency and precision of drought-tolerant breeding of sugar beet through research methods at the molecular level.

4.2 Genome-wide association analysis of drought tolerance-related traits in sugar beet

With the development of genotyping technology, SNP markers are characterized by high throughput, low cost and high accuracy. The application of genome-wide association analysis is universal, and it can be used to analyze both qualitative and quantitative traits (Alqudah et al., 2020). In 175 rice samples, 13 significant SNPs were screened out and 50 genes were identified, and 10 of these genes were associated with drought and/or other abiotic stress tolerance (Pantaliao et al., 2016). It was detected 132 single nucleotide polymorphisms in sesame under drought stress, and 13 potential candidate genes were found to be related with drought tolerance (Li et al., 2018). At present, GWAS has been reported in the drought tolerance studies of many crops, but it is rarely reported in sugar beet. In this study, we explored the functional genes related to drought tolerance of sugar beet through genome-wide association analysis. 328 sugar beet germplasms were analyzed, and phylogenetic tree and PCA analysis showed that there was no obvious population stratification, indicating that the heredity of sugar beet used here was less affected by geographical factors. Considering that MLM model can better control false positive results, we used this model to excavate significant sites (Zhang et al., 2010).

4.3 Potential candidate genes associated with sugar beet drought tolerance

A total of 11 SNPs loci were identified by genome-wide association analysis and 25 potential candidate genes were screened out. It was found that 14 genes had been reported to be related with drought tolerance. We conducted qRT-PCR to verify gene expression levels in the two groups of drought tolerance extreme materials and found that only BVRB_2g036960 had no difference in expression between the two groups of materials, and the expression levels of 2 genes in the drought sensitive group were significantly higher than those in the drought tolerant group. The expression levels of 11 genes in drought-tolerant group were significantly higher than those in drought-sensitive group. Therefore, these 13 differentially expressed genes might be key genes which can be used in our next research on drought tolerance in sugar beet.

BVRB_7g160020 regulated RLPH2 is a tyrosine phosphorylase of the PPP family, and tyrosine phosphorylation is essential for signaling at the cytoplasmic structural domain (Mühlenbeck et al., 2021). We also found that the gene BVRB_7g160030, which is adjacent to the gene BVRB_7g160020, is involved in straight-chain starch synthesis in sugar metabolism. α (1–6) glycosidic bond branching increases the water solubility of glycogen for storage and increases the number of nonreducing ends, making it easier for biological organisms to aggregate when glycogen is needed for energy supply, and the 1,4-alpha-glucan-branching enzyme regulates glycogen metabolism and participates in the abiotic stress response (Li et al., 2019).

We found that BVRB_7g166470 is involved in the cell wall composition. BVRB_3g048910 can regulate methionine, which can mitigate membrane damage by peroxides. It has been shown that the probable methyltransferase PMT3 is involved in the synthesis of phosphatidylcholine (PC), which is an important precursor of lipid signalling (Baldessari and Iglesias, 2012) and a ligand for regulatory proteins (Inatsugi et al., 2009). PC appears to be a major contributor to the adaptive response of plants to environmental stress, and it plays an important role in maintaining cell membrane integrity in complex plant species. The ras GTPase-activating protein-binding protein 2 associated with BVRB_7g166680 regulates ascorbic acid, which has strong reducing properties that effectively mitigate peroxide damage (Khazaei et al., 2020). BVRB_9g203030 is eukaryotic initiation factor 4A-9, which promotes ATP hydrolysis to provide energy for material transport and signaling in plants, and myosin, which is a component of the cytoskeleton to provide energy for organelle movement, material transport and apical growth (Han et al., 2021). BVRB_2g033710 is involved in the components of the cytoskeleton that provide energy for organelle movement, material transport and apex growth (Han et al., 2022). BVRB_3g048900 on chromosome 3 was found to regulate abscisic acid, and BVRB_5g114810 on chromosome 5 was related to MYB transcription factors. Abscisic acid can regulate MYB, which can promote lateral root growth, enhance root water absorption, and induce old leaf death to reduce water loss under drought conditions (Moosavi et al., 2017). The ENHANCER OF AG-4 protein regulated by BVRB_9g203000 can induce apoptosis and prolong nutritional growth (Moosavi et al., 2017; Alqudah et al., 2020).

BVRB_9g203050 on chromosome 9 regulates HCF136, which is essential for PSII biogenesis and stability (Meurer et al., 1998); the HCF136 protein is essential for the formation of early intermediates in PSII assembly, including D2 (psbD) and cytochrome b559 (Plücken et al., 2002). D2 can repair PSII and D1 expression dysfunction and maintain the stability of photosynthesis (Kanamaru and Tanaka, 2004). Protein LOW PSII ACCUMULATION (BVRB_7g166480), a required chaperone for efficient PSII assembly, is a component of the chloroplast cyst-like membrane; it binds to psbA during the biogenesis of PSII and plays an important role in repairing PSII damage to maintain photosystem stability (Peng et al., 2006). Photosystem I (PSI) reaction centre subunit IV (BVRB_7g166490) is capable of stabilizing the interaction between PsaC and the PSI core, assisting in the docking of iron oxidation reduction protein to PSI and interacting with iron oxidation reduction protein-NADP oxidoreductase (Dutta et al., 2014).

4.4 Potential impact on sugar beet breeding

Although conventional breeding can improve the performance of sugar beet in drought tolerance, artificial populations are limited by population size, insufficient recombination, and often large locus intervals, making selection both time-consuming and inefficient (He et al., 2010). Most of the natural populations used here have been reproduced over long periods of time and multiple generations, with adequate genomic recombination and accurate linkage localization, and these new genomic resources will facilitate genomic evolution in sugar beet and provide a scientific and theoretical basis for relevant research in drought tolerance.

5 Conclusion

In this study, 11 loci potentially associated with drought tolerance in sugar beet were screened using a mixed linear model (MLM). Genome-wide association analysis combined with qRT-PCR revealed that 13 genes were significantly differentially expressed between the two extreme groups of materials, with BVRB_7g166480 being 46.69-fold more expressed in the drought-tolerant group than in the drought-sensitive group. It prompted us to further investigate the role of these interesting genes in sugar beet drought tolerance improvement in the future.

References

• 1

  AhnC. S.HanJ.-A.LeeH.-S.LeeS.PaiH.-S. (2011). The PP2A regulatory subunit Tap46, a component of the TOR signaling pathway, modulates growth and metabolism in plants. Plant Cell.23 (1), 185–209. 10.1105/tpc.110.074005

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  AlqudahA. M.SallamA.BaenzigerP. S.BörnerA. (2020). Gwas: Fast-forwarding gene identification and characterization in temperate cereals: Lessons from barley–a review. J. Adv. Res.22, 119–135. 10.1016/j.jare.2019.10.013

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 3

  BaiX.WangX.WangY.WeiY.FuY.RaoJ.et al (2022). Genome-wide association study of six forage traits in ramie (boehmeria nivea L. Gaud). Plants11 (11), 1443. 10.3390/plants11111443

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 4

  BaldessariA.IglesiasL. (2012). in Lipases and phospholipases, methods and protocols. Editor SandovalG. (New York: Humana Press), 457–470.

  • Google Scholar
  • View reference in article

• 5

  BlochD.HoffmannC. M.MärländerB. (2006). Impact of water supply on photosynthesis, water use and carbon isotope discrimination of sugar beet genotypes. Eur. J. Agron.24 (3), 218–225. 10.1016/j.eja.2005.08.004

  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  BreriaC. M.HsiehC. H.YenJ.-Y.NairR.LinC.-Y.HuangS.-M.et al (2020). Population structure of the world vegetable center mungbean mini core collection and genome-wide association mapping of loci associated with variation of seed coat luster. Trop. Plant Biol.13 (1), 1–12. 10.1007/s12042-019-09236-0

  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  ByrneR. P.van RheenenW.van den BergL. H.VeldinkJ. H.McLaughlinR. L. (2020). Dutch population structure across space, time and GWAS design. Nat. Commun.11 (1), 4556. 10.1038/s41467-020-18418-4

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 8

  CadzowM.BoocockJ.NguyenH. T.WilcoxP.MerrimanT. R.BlackM. A. (2014). A bioinformatics workflow for detecting signatures of selection in genomic data. Front. Genet.5, 293. 10.3389/fgene.2014.00293

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 9

  CaiS.JiangG.YeN.ChuZ.XuX.ZhangJ.et al (2015). A key ABA catabolic gene, OsABA8ox3, is involved in drought stress resistance in rice. PLoS One10 (2), e0116646. 10.1371/journal.pone.0116646

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  CarilloP.GibonY. (2011). Protocol: Extraction and determination of proline. PrometheusWiki2011, 1–5.

  • Google Scholar
  • View reference in article

• 11

  DengK.YuL.ZhengX.ZhangK.WangW.DongP.et al (2016). Target of rapamycin is a key player for auxin signaling transduction in Arabidopsis. Front. Plant Sci.7, 291. 10.3389/fpls.2016.00291

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 12

  DramadriI. O.NkaluboS. T.KramerD. M.KellyJ. D. (2021). Genome‐wide association analysis of drought adaptive traits in common bean. Crop Sci.61 (5), 3232–3253. 10.1002/csc2.20484

  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  DuttaS.TeresinskiH. J.SmithM. D. (2014). A split-ubiquitin yeast two-hybrid screen to examine the substrate specificity of atToc159 and atToc132, two Arabidopsis chloroplast preprotein import receptors. PloS One9 (4), e95026. 10.1371/journal.pone.0095026

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  DwiningsihY.ThomasJ.KumarA.GuptaC.GillN.RuizC.et al (2022). Identification of QTLs and candidate loci associated with drought-related traits of the K/Z RIL rice population.

  • Google Scholar
  • View reference in article

• 15

  FugateK. K.LaftaA. M.EideJ. D.LiG.LulaiE. C.OlsonL. L.et al (2018). Methyl jasmonate alleviates drought stress in young sugar beet (Beta vulgaris L) plants. J. Agron. Crop Sci.204 (6), 566–576. 10.1111/jac.12286

  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  HanH.VerstraetenI.RoosjenM.MazurE.RýdzaN.HajnýJ.et al (2021). Rapid auxin-mediated phosphorylation of Myosin regulates trafficking and polarity in Arabidopsis. BioRxiv. 10.1101/2021.04.13.439603

  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  HanM.XuM.SuT.WangS.WuL.FengJ.et al (2022). Transcriptome analysis reveals critical genes and pathways in carbon metabolism and ribosome biogenesis in poplar fertilized with glutamine. Int. J. Mol. Sci.23 (17), 9998. 10.3390/ijms23179998

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 18

  HassaniM.HeidariB.StevanatoP. (2020). Combining abilities of sugar beet genotypes for root‐and sugar‐related traits under multi‐environment trials. Plant Breed.139 (1), 192–206. 10.1111/pbr.12755

  • CrossRef
  • Google Scholar
  • View reference in article

• 19

  HeZ.XiaX.BonjeanA. P.LiF.WangB.WangJ. (2010). The potential of carcinoembryonic antigen, p53, Ki-67 and glutathion Stransferase-π as clinico-histopathological markers for colorectal cancer. Cereals China688 (787), 51–57. 10.1016/S1674-8301(10)60008-5

  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  Hoyos‐VillegasV.SongQ.KellyJ. D. (2017). Genome‐wide association analysis for drought tolerance and associated traits in common bean. Plant Genome10 (1), 1–17. 10.3835/plantgenome2015.12.0122

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 21

  HuangX.CloutierS.LycarL.RadovanovicN.HumphreysD.NollJ.et al (2006). Molecular detection of QTLs for agronomic and quality traits in a doubled haploid population derived from two Canadian wheats (Triticum aestivum L). Theor. Appl. Genet.113, 753–766. 10.1007/s00122-006-0346-7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  InatsugiR.KawaiH.YamaokaY.YuY.SekiguchiA.NakamuraM.et al (2009). Isozyme-specific modes of activation of CTP: Phosphorylcholine cytidylyltransferase in Arabidopsis thaliana at low temperature. Plant Cell. Physiology50 (10), 1727–1735. 10.1093/pcp/pcp115

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  IslamM. J.KimJ. W.BegumM. K.SohelM. A. T.LimY.-S. (2020). Physiological and biochemical changes in sugar beet seedlings to confer stress adaptability under drought condition. Plants9 (11), 1511. 10.3390/plants9111511

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 24

  JhaU. C.NayyarH.PalakurthiR.JhaR.ValluriV.BajajP.et al (2021). Major QTLs and potential candidate genes for heat stress tolerance identified in chickpea (Cicer arietinum L). Front. Plant Sci.12, 655103. 10.3389/fpls.2021.655103

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  JiaoZ.XuW.NongQ.ZhangM.JianS.LuH.et al (2021). An integrative transcriptomic and metabolomic analysis of red pitaya (Hylocereus polyrhizus) seedlings in response to heat stress. Genes.12 (11), 1714. 10.3390/genes12111714

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 26

  JohnsonA. A.VenkateshB.JayammaN.ReddyB. M.PandurangaiahM.SudhakarC. (2022). R2R3 MYB transcription factor, AhMYB94 plays a crucial role in stress tolerance of a salt susceptible groundnut cultivar.

  • Google Scholar
  • View reference in article

• 27

  KanamaruK.TanakaK. (2004). Roles of chloroplast RNA polymerase sigma factors in chloroplast development and stress response in higher plants. Biosci. Biotechnol. Biochem.68 (11), 2215–2223. 10.1271/bbb.68.2215

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 28

  KangY.SakirogluM.KromN.Stanton‐GeddesJ.WangM.LeeY. C.et al (2015). Genome‐wide association of drought‐related and biomass traits with HapMap SNPs in Medicago truncatula. Plant, Cell. & Environ.38 (10), 1997–2011. 10.1111/pce.12520

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 29

  KhazaeiZ.EsmaielpourB.EstajiA. (2020). Ameliorative effects of ascorbic acid on tolerance to drought stress on pepper (Capsicum annuum L) plants. Physiology Mol. Biol. Plants26 (8), 1649–1662. 10.1007/s12298-020-00846-7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 30

  KumarD.DasB.SenR.KunduP.MannaA.SarkarA.et al (2015). Andrographolide analogue induces apoptosis and autophagy mediated cell death in U937 cells by inhibition of PI3K/Akt/mTOR pathway. PLoS One10 (10), e0139657. 10.1371/journal.pone.0139657

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 31

  KumarD.SharmaS.SharmaR.PundirS.SinghV. K.ChaturvediD.et al (2021). Genome-wide association study in hexaploid wheat identifies novel genomic regions associated with resistance to root lesion nematode (Pratylenchus thornei). Sci. Rep.11 (1), 3572. 10.1038/s41598-021-80996-0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 32

  LanderE. S.LintonL. M.BirrenB.NusbaumC.ZodyM. C.BaldwinJ.et al (2001). Initial sequencing and analysis of the human genome. Nature409 (6822), 860–921. 10.1038/35057062

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 33

  LiH.DurbinR. (2009). Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics25 (14), 1754–1760. 10.1093/bioinformatics/btp324

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 34

  LiD.DossaK.ZhangY.WeiX.WangL.ZhangY.et al (2018). GWAS uncovers differential genetic bases for drought and salt tolerances in sesame at the germination stage. Genes.9 (2), 87. 10.3390/genes9020087

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 35

  LiX.DingM.WangM.YanF.ZhangZ.LiangW.et al (2019). Gene cloning of 1,4-α-glucan branching enzyme (GBE) from nostoc flagelliforme and its differential expression in response to drought stress. Mol. Plant Breed.17 (19), 6305–6313.

  • Google Scholar
  • View reference in article

• 36

  LiuD.GaoZ.LiJ.YaoQ.TanW.XingW.et al (2023). Effects of cadmium stress on the morphology, physiology, cellular ultrastructure, and BvHIPP24 gene expression of sugar beet (Beta vulgaris L). Int. J. Phytoremediation25 (4), 455–465. 10.1080/15226514.2022.2090496

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 37

  LivakK. J.SchmittgenT. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods25 (4), 402–408. 10.1006/meth.2001.1262

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 38

  ManessN. (2010). “Extraction and analysis of soluble carbohydrates,” in Plant stress tolerance (Springer), 341–370.

  • Google Scholar
  • View reference in article

• 39

  MathewI.ShimelisH.ShayanowakoA. I. T.LaingM.ChaplotV. (2019). Genome-wide association study of drought tolerance and biomass allocation in wheat. PloS One14 (12), e0225383. 10.1371/journal.pone.0225383

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 40

  MeurerJ.PlückenH.KowallikK. V.WesthoffP. (1998). A nuclear-encoded protein of prokaryotic origin is essential for the stability of photosystem II in Arabidopsis thaliana. EMBO J.17 (18), 5286–5297. 10.1093/emboj/17.18.5286

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 41

  MoosaviS. G. R.RamazaniS. H. R.HemayatiS. S.GholizadeH. (2017). Effect of drought stress on root yield and some morpho-physiological traits in different genotypes of sugar beet (Beta Vulgaris L). J. Crop Sci. Biotechnol.20 (3), 167–174. 10.1007/s12892-017-0009-0

  • CrossRef
  • Google Scholar
  • View reference in article

• 42

  MühlenbeckH.BenderK. W.ZipfelC. (2021). Importance of tyrosine phosphorylation for transmembrane signaling in plants. Biochem. J.478 (14), 2759–2774. 10.1042/BCJ20210202

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 43

  PantaliaoG. F.NarcisoM.GuimarãesC.CastroA.ColombariJ. M.BreseghelloF.et al (2016). Genome wide association study (GWAS) for grain yield in rice cultivated under water deficit. Genetica144 (6), 651–664. 10.1007/s10709-016-9932-z

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 44

  PengL.MaJ.ChiW.GuoJ.ZhuS.LuQ.et al (2006). LOW PSII ACCUMULATION1 is involved in efficient assembly of photosystem II in Arabidopsis thaliana. Plant Cell.18 (4), 955–969. 10.1105/tpc.105.037689

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 45

  PlückenH.MüllerB.GrohmannD.WesthoffP.EichackerL. A. (2002). The HCF136 protein is essential for assembly of the photosystem II reaction center in Arabidopsis thaliana. FEBS Lett.532 (1-2), 85–90. 10.1016/s0014-5793(02)03634-7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 46

  PolashM. A. S.SakilM. A.Tahjib-Ul-ArifM.HossainM. A. (2018). Effect of salinity on osmolytes and relative water content of selected rice genotypes. Trop. Plant Res.5 (2), 227–232. 10.22271/tpr.2018.v5.i2.029

  • CrossRef
  • Google Scholar
  • View reference in article

• 47

  RatnapriyaR.SosinaO. A.StarostikM. R.KwicklisM.KapphahnR. J.FritscheL. G.et al (2019). Retinal transcriptome and eQTL analyses identify genes associated with age-related macular degeneration. Nat. Genet.51 (4), 606–610. 10.1038/s41588-019-0351-9

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 48

  SadeghianS.YavariN. (2004). Effect of water‐deficit stress on germination and early seedling growth in sugar beet. J. Agron. Crop Sci.190 (2), 138–144. 10.1111/j.1439-037x.2004.00087.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 49

  ShariatipourN.HeidariB.TahmasebiA.RichardsC. (2021). Comparative genomic analysis of quantitative trait loci associated with micronutrient contents, grain quality, and agronomic traits in wheat (Triticum aestivum L). Front. Plant Sci.12, 709817. 10.3389/fpls.2021.709817

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 50

  SoskicB.Cano-GamezE.SmythD. J.RowanW. C.NakicN.Esparza-GordilloJ.et al (2019). Chromatin activity at GWAS loci identifies T cell states driving complex immune diseases. Nat. Genet.51 (10), 1486–1493. 10.1038/s41588-019-0493-9

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 51

  SpitzD. R.OberleyL. W. (1989). An assay for superoxide dismutase activity in mammalian tissue homogenates. Anal. Biochem.179 (1), 8–18. 10.1016/0003-2697(89)90192-9

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 52

  SunD.CenH.WengH.WanL.AbdallaA.El-ManawyA. I.et al (2019). Using hyperspectral analysis as a potential high throughput phenotyping tool in GWAS for protein content of rice quality. Plant Methods15 (1), 54–16. 10.1186/s13007-019-0432-x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 53

  TanW.LiW.LiJ.LiuD.XingW. (2023). Drought resistance evaluation of sugar beet germplasms by response of phenotypic indicators. Plant Signal. Behav.18 (1), 2192570. 10.1080/15592324.2023.2192570

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 54

  UhrigR. G.LabanderaA.-M.MuhammadJ.SamuelM.MoorheadG. B. (2016). Rhizobiales-like phosphatase 2 from Arabidopsis thaliana is a novel phospho-tyrosine-specific phospho-protein phosphatase (PPP) family protein phosphatase. J. Biol. Chem.291 (11), 5926–5934. 10.1074/jbc.M115.683656

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 55

  WangQ.SajjaU.RosloskiS.HumphreyT.KimM. C.BombliesK.et al (2007). HUA2 caused natural variation in shoot morphology of A. thaliana. Curr. Biol.17 (17), 1513–1519. 10.1016/j.cub.2007.07.059

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 56

  WangM.FangZ.YooB.BejeranoG.PeltzG. (2021a). The effect of population structure on murine genome-wide association studies. Front. Genet.12, 745361. 10.3389/fgene.2021.745361

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 57

  WangX.NiuY.ZhengY. (2021b). Multiple functions of MYB transcription factors in abiotic stress responses. Int. J. Mol. Sci.22 (11), 6125. 10.3390/ijms22116125

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 58

  XiongY.McCormackM.LiL.HallQ.XiangC.SheenJ. (2013). Glucose–TOR signalling reprograms the transcriptome and activates meristems. Nature496 (7444), 181–186. 10.1038/nature12030

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 59

  XuY.YangT.ZhouY.YinS.LiP.LiuJ.et al (2018). Genome-wide association mapping of starch pasting properties in maize using single-locus and multi-locus models. Front. Plant Sci.9, 1311. 10.3389/fpls.2018.01311

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 60

  YangH.WangK. (2015). Genomic variant annotation and prioritization with ANNOVAR and wANNOVAR. Nat. Protoc.10 (10), 1556–1566. 10.1038/nprot.2015.105

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 61

  YangK.ChenG.XianJ.ChenW. (2022a). Varying relationship between vascular plant leaf area and leaf biomass along an elevational gradient on the eastern qinghai-tibet plateau. Front. Plant Sci.13, 824461. 10.3389/fpls.2022.824461

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 62

  YangY.DongS.MiaoH.LiuX.DaiZ.LiX.et al (2022b). Genome-wide association studies reveal candidate genes related to stem diameter in cucumber (cucumis sativus L). Genes.13 (6), 1095. 10.3390/genes13061095

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 63

  ZengA.ChenP.KorthK.HancockF.PereiraA.BryeK.et al (2017). Genome-wide association study (GWAS) of salt tolerance in worldwide soybean germplasm lines. Mol. Breed.37 (3), 30–14. 10.1007/s11032-017-0634-8

  • CrossRef
  • Google Scholar
  • View reference in article

• 64

  ZhangZ.ErsozE.LaiC.-Q.TodhunterR. J.TiwariH. K.GoreM. A.et al (2010). Mixed linear model approach adapted for genome-wide association studies. Nat. Genet.42 (4), 355–360. 10.1038/ng.546

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 65

  ZhengJ.ZhangZ.GongZ.LiangY.SangZ.XuY.et al (2021). Genome-wide association analysis of salt-tolerant traits in terrestrial cotton at seedling stage. Plants11 (1), 97. 10.3390/plants11010097

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article