A genome-wide association study of lateral root number for Asian cotton (Gossypium arboreum L.)
Journal of Cotton Research volume 5, Article number: 19 (2022)
The lateral root is one of the most important organs that constitute the root architecture system in plants. It can directly affect the contact area between plants and soil and plays an important role in plant structural support and nutrient absorption. Optimizing root architecture systems can greatly increase crop yields. This study was designed to identify the molecular markers and candidate genes associated with lateral root development in cotton and to evaluate correlations with yield and disease traits.
The number of lateral roots for 14-day old seedlings was recorded for 215 Gossypium arboreum accessions. A correlation analysis showed that the number of lateral roots positively correlates with the sympodial branch node and seed index traits, but negatively correlates with lint percentage. A Genome-wide association study (GWAS) identified 18 significant SNPs with 19 candidate genes associated with the lateral root number. Expression analysis identified three genes (FLA12, WRKY29, and RBOHA) associated with lateral root development.
GWAS analysis identified key SNPs and candidate genes for lateral root number, and genes of FLA12, WRKY29, and RBOHA may play a pivotal role in lateral root development in Asian cotton.
The lateral root, one of the most important organs in root system architecture, plays a pivotal role in water and nutrient acquisition that can influence crop productivity. In rice (Oryza sativa L.), the lateral root defective mutant osiaa11 and the root hairs defective mutant osrhl1 were used to demonstrate that lateral roots and not root hairs played a pivotal role in high Mn and Cd uptake (Yu et al. 2021). In maize (Zea mays), Jia et al. (2018) proved that higher lateral root branch density has greater phosphorus acquisition from low phosphorus soil.
Lateral roots are initiated from the root pericycle cells by asymmetric cell division (Casero et al. 1993). Research progresses on quantitative trait loci (QTLs) of lateral root development have been reported in crops like rice (Oryza sativa L.) (Niones et al. 2015), rap (Brassica napus L.) (Ying et al. 2016), maize (Zea mays L.) (Wang et al. 2019), and soybean (Glycine max (Linn.) Merr.) (Liang et al. 2017), but little knowledge was discovered in cotton. The lateral root development is a complex process that is mainly modulated by endogenous hormones. Auxin transport plays a pivotal role in root branching with different auxin “signal modules” regulating the various stages of lateral root development (Marhavy et al. 2014; Talboys et al. 2014). A recent study found that the AUXIN RESPONSE FACTOR 7 and its inhibitor IAA18/POTENT formed a negative regulatory loop circuit to drive the root clock in Arabidopsis, which affected the spacing of lateral roots along the primary root (Perianez-Rodriguez et al. 2021). Lv et al. (2021) overexpressed the transcription factor ERF13 in Arabidopsis, resulting in erf13 mutants with greater lateral root emergence, and showed that MPK14-mediated auxin signals controlled lateral root development through ERF13-regulated very-long-chain fatty acid biosynthesis. The auxin signal modules IAA14/SLR (SOLITARY ROOT)-ARF7-ARF19 (Fukaki et al. 2010), IAA12/BDL (BODENLOS)-ARF5 (Ive et al. 2014), PIN proteins (Lewis et al. 2011; Inahashi et al. 2018) have also been reported to be involved in the lateral root initiation and development. Other plant hormones, such as gibberellin (GA) (Farquharson 2010; Hetherington et al. 2021), abscisic acid (ABA) (Duan et al. 2013; Zhao et al. 2014), ethylene (Ivanchenko et al. 2008; Prasad et al. 2010), and jasmonate (JA) (Cai et al. 2015; Kuo et al. 2020) were also associated with the lateral root formation. In cotton (Gossypium hirsutum), the plant growth regulator mepiquat chloride promoted lateral root formation by modulating GA, ABA, and auxin homeostasis (Wu et al. 2019).
Although plant hormones play an important role in lateral root development, the mechanisms of lateral root initiation in cotton have not been elucidated. Genome-wide association study (GWAS) is a typical method to evaluate the association between millions of single nucleotide polymorphisms (SNPs) markers and a phenotype of interest in large populations for animals or plants. In cotton, many important agronomic traits like yield (Du et al. 2018), fiber quality (Fang et al. 2017; He et al. 2021) and disease-resistant related traits (Li et al. 2017) have been successfully dissected their underpinnings by GWAS. Gossypium arboreum (Asian cotton) is a diploid cotton species that can be used as a model to study root architecture development. DNA sequencing data is available for 215 G. arboreum accessions (Du et al. 2018). A GWAS analysis was conducted in the present study to detect the significant loci for seedling lateral root number. Results of this study could provide targets for cotton breeders to optimize the cotton root system structure to enhance yield potential.
Materials and methods
Determination of lateral root initiation rate
Lateral root development was initially evaluated in the Asian cotton standard line Shixiya1 (Li et al. 2014) to identify the period with the highest frequency of lateral root initiation. Seeds of Shixiya1 were surface sterilized with 15% H2O2 for ~ 30 min, then rinsed five times with steriled water. Seeds were germinated on moisten filter papers and then transplanted into hydroponic containers using a half-strength of Hoagland nutrient solution. The lateral roots were manually counted 4, 8, 12, 20 days after planting for each seedling. Five replicates were set, and each replicate included six plants. The average daily growth of lateral roots was calculated by increased lateral root number/growth days. The differences among growth periods were analyzed by t-test.
Lateral root evaluation for accessions
This evaluation used 215 G. arboreum accessions that had DNA sequencing and phenotypic data for 11 yield and disease-related traits (Du et al. 2018). Seeds for the accessions were provided by the National Medium-term Gene Bank of Cotton in China (Anyang, Henan). Thirty-six seeds with full and uniform size for each accession were first selected for sowing. Then, all accessions were evaluated using the seed germination bag method according to the manual instructions (Zhao et al. 2021). Each accession included three independent replicates (one replicate per bag), and each bag contained 10∼12 plants (Shixiya1 included as a control). The experiment was conducted in a greenhouse (Temperature: 28℃; photoperiod: light/night = 16 h/8 h; relative humidity: 60% ~ 65%). The number of lateral roots for each plant was counted 14 days after planting as described for the Shixiya1 evaluation. The mean number of lateral roots and their replicates were calculated and were used for further analysis (Additional file 1: Table S1). A basic description analysis included the range, coefficient of variation, Skewness/Kurtosis was performed by GraphPad Prism 8 (Francis 1995).
Correlation analysis of lateral root development with phenotypic traits
The lateral root data for the 215 accessions were correlated with 11 yield and disease-related traits. The 11 yield and disease-related traits were (1) BW_AY (Boll weight in Anyang (Henan)), (2) BW_SY (Boll weight in Sanya (Hainan)), (3) BW_AKS (Boll weight in Aksu (Xinjiang)), (4) SBN_AY (Sympodial branch node in Anyang (Henan)), (5) SBN_AKS (Sympodial branch node in Aksu (Xinjiang)), (6) FWDI_AY (Fusarium wilt disease index in Anyang (Henan)), (7) LP_AY (Lint percentage in Anyang (Henan)), (8) LP_SY (Lint percentage in Sanya (Hainan)), (9) LP_AKS (Lint percentage in Aksu (Xinjiang)), (10) SI_AY (Seed index in Anyang (Henan)), (11) VWDI_AY (Verticillium wilt disease index in Anyang (Henan)). A pheatmap package method was used to plot the picture of the correlation matrix in R language (R Core Team 2012).
GWAS for lateral root development
The genotype data were derived from Du et al. (2018), and DNA sequences of the same 215 accessions were used for the lateral root number in this study. The Burrows-Wheeler Aligner program (BWA, ver. 0.7.10) (Li and Durbin 2009) and GATK toolkit (ver. 3.2–2) (Mckenna et al. 2010) were used to perform the reads mapping and SNP calling with the default parameters, respectively. Software ANNOVAR tool (Kai et al. 2010) was used to annotate the identified SNPs based on the newly updated G. arboreum genome (Yang et al. 2019) annotation information. After filtering, a total of 1 425 003 high-quality SNP markers (MAF (Minor allele frequency for the genetic variant of interest) > 0.05, missing rate per site < 10%) were selected to perform the GWAS for lateral root trait with the EMMAX (Efficient mixed-model association expedited) model (Kang et al. 2010). The SNP density plot was drawn by R language with “CMplot” package (Turner 2018). The significance for LRN was first set as −log10(P) > 6.15 (1/total SNPs used, –log10(P) = 6.15) (Li et al. 2017), but no effective loci were identified in all the 13 chromosome regions with this threshold. So we manually adjusted their threshold to P < 2.81 × 10–6 (4/1 425 003 used, and − log10(P) = 5.55), and those SNPs with − log10(P) > 5.55 were defined as significant SNPs. The candidate genes were identified by the physical locations of the significant SNP and the related stronger LD (linkage disequilibrium) block regions.
LD block and haplotype analysis for lateral root development
The confidence intervals for the lateral root data were identified by IGV_2.3.91 (Helga et al. 2013) and were plotted by R language with the “qqman” package (Turner 2018). The interval regions were further performed using LD block and Neighbor-Joining analysis in Tassel software (Bradbury et al. 2007). The stronger LD block interval regions were further determined using the haplotype analysis method described by Dai et al. (2020). The different haplotypes were roughly estimated by the genotypes of different groups with the color scale method EXCEL 2010 (https://www.microsoft.com/en-us/microsoft-365/excel).
RNA-seq data of G. arboreum cultivar Shixiya1 had been sequenced by Huang et al. (2020) (Accession: PRJNA594268 ID: 594268), and we downloaded the TPM (Transcripts per million) expression data of candidate genes from https://cottonfgd.org/. The root expression value of the seedling stage was utilized in this study. These candidate genes' tissue and organ expression levels were plotted in R language with a pheatmap package.
Lateral root development for Shixiya1
The number of lateral roots observed varied significantly across the four-time periods (Fig. 1). No lateral roots were recorded 4 days after planting, whereas lateral roots were clearly visible 8 days after planting with a mean of 34.20 ± 7.60 recorded. The number of lateral roots continued to rapidly increase with a mean of 68.00 ± 9.98 recorded 12 days after planting. At 20 days after planting, the rate of lateral root initiation had slowed with a mean of 85.80 ± 10.57 observed. The average daily growth of lateral roots was kept stable on 4∼8 d and 8∼12 d, but it rapidly slowed on 12∼20 d. Using these data, a 14-day time period was selected to evaluate the 215 accessions.
Lateral root data analysis for 215 accessions
To more conveniently describe the lateral root number, we abbreviate this trait as LRN. Furthermore, LRN_AV, LRN_R1, LRN_R2, LRN_R3 represented the mean and three replicates of LRN, respectively. LRN_AV ranged from 18.86 to 74.05 across the 215 accessions with a mean of 43.21 and a coefficient of variation of 26.81% (Additional file 1: Table S2). Replicates of LRN showed a similar range, mean, and variation with LRN_AV (Fig. 2a; Additional file 1: Table S2). LRN_AV and its three replicates were biased towards the normal distribution due to their frequency distribution and Skewness/Kurtosis values, and this distribution was suitable for GWAS analysis (Fig. 2b; Additional file 1: Table S2). We are curious about the correlation of the lateral root trait with yield and disease resistance-related traits. So, a correlation was performed. The result showed that LRN_AV positively correlated with the sympodial branch node and seed index traits. Notably, a negative correlation with the lint percentage was observed (Fig. 2c). However, LRN_AV showed a weak correlation with the other yield traits like boll weight and disease-related (VWDI) traits.
GWAS for the number of lateral roots
A total of 1 425 003 high-quality SNP markers were screened out to perform the GWAS in this study, and those SNP densities were shown in Additional file 2: Fig. S1. The quantile–quantile (QQ) plots of LRN_AV, LRN_R1, LRN_R2, and LRN_R3 indicated that the EMMAX could be used to identify signals for lateral root number in Asian cotton (Additional file 3: Fig. S2). The Manhattan plots of LRN_AV and its three replicates showed stronger signals on Chr01, Chr07, Chr13, and Scaffold regions (Chr14). The strongest signal was found on Scaffold regions, and a continuous, more robust signal interval on Chr01 (Chr01:78.20∼78.90 Mb) was detected in LRN_AV and all the three replicates (Fig. 3).
Identification of significant SNPs and candidate genes
Eighteen significant SNPs were detected for the lateral root trait with 7 SNPs located on the Chr01, 5 SNPs located on Scaffold regions, 3 SNPs located on Chr07, and 2 SNPs located on Chr13 (Additional file 1: Table S3). The strongest signal was detected on SNP Chr14_23233813 with a − log10(P) value of 6.18 and 6.68 for LRN_AV and LRN_R3, respectively. The allelic variant of this site was G/A, and this signal was located on the intergenic region of Ga14G0299 and Ga14G2065. The SNPs for the continuous stronger signal interval on Chr01 included SNP Chr01_78558940, SNP Chr01_78572007, SNP Chr01_78703010, SNP Chr01_78751608, SNP Chr01_78758771, and SNP Chr01_78784026. The LD block of this interval (Chr01:78.20–78.90 Mb) showed a stronger LD block on the region of Chr01:78 284 704–78 804 477, and this LD block was referred to as LD_Chr01 (Fig. 4a, b). Haplotype analysis of the LD_Chr01 region (Fig. 4c) divided the 215 accessions into three haplotypes (HapChr01_1, HapChr01_2, HapChr01_3), including 36, 68, and 111 accessions, respectively. The LRN_AV recorded for HapChr01_2 accessions was significantly higher than HapChr01_1 and HapChr01_3, and the mean number of lateral roots for HapChr01_3 was the lowest (Fig. 4d). The LRN_R1, LRN_R2, and LRN_R3 also showed the same result as LRN_AV (Fig. 4d).
Nineteen candidate genes for LRN were detected for the five chromosomal regions (Additional file 1: Table S4). Seven candidate genes were detected on Chr01, six on the scaffold region, two on Chr05, two on Chr07, and two on Chr13. The LD_Chr01 block was associated with the four candidate genes WRKY29, WDR7, PP1, and At5g41590. The strongest signals on the Scaffold region were associated with HSP23.9 and Ga14G2065. Other associated candidate genes and their significant SNPs were listed in the Additional file 1: Table S4.
TPM expression analysis of candidate genes
Ga01G0035 (S6PDH), Ga01G1695 (WRKY29), Ga01G1696 (WDR7), Ga01G1697 (PP1), Ga01G1699 (RBOHA), Ga07G2390 (MIF2) and Ga14G1993 (FLA12) were expressed higher in the root organization with the highest TPM value (111.49) recorded for Ga14G1993 (FLA12) (Fig. 5, Additional file 1: Table S5).
Cotton is a tap root crop with a deep primary root and wide distribution of lateral roots to form a strong absorption network in the soil. Plant root can directly influence the development of above-ground organs, and different root traits may have a positive/negative correlation with yield-related traits in crops. In upland cotton, the root weight and taproot length appeared a positive correlation with lint yield, whereas the tap root showed a negative correlation with lint yield (Li et al. 2003). So far, there are no reports on the relationship between the lateral root number and yield-related traits in cotton. In this study, we found that the lateral root number negatively correlates with lint percentage and displays a positive correlation with sympodial branch node and seed index traits. Verticillium wilt (VW) disease is caused by the soil-borne fungus Verticillium dahliae, and the plant roots were usually firstly influenced when the VW disease happened. Whether the lateral root number (LRN) has a correlation with VW disease index (VWDI)? Our correlation results showed a weak correlation between the LRN and VWDI. Our findings will deepen people's understanding of the relationship between lateral root development and yield-related traits in cotton.
Research progresses on lateral root development were reported in several crops (Niones et al. 2015; Ying et al. 2016; Wang et al. 2019), but so far, no QTLs/SNPs of lateral root number was discovered in Asian cotton. While, in this study, we identified 18 significant signals in total and detected a stronger LD block on Chr01 for lateral root number in Asian cotton, which will greatly enrich researchers' understanding of molecular markers that control lateral root development in cotton.
In the present study, seven candidate genes were expressed higher in root organs, and three genes (FLA12, WRKY29, and RBOHA) were closely associated with root development. FLA12 encodes for Fasciclin-like arabinogalactan protein 12 (Johnson et al. 2011). Arabinogalactan proteins FLAs are involved in plant growth and development and play an essential role in many plant species (Seifert et al. 2014; Zang et al. 2015). In Arabidopsis, At-FLA4 is necessary for the average growth of wild-type roots under moderately elevated salinity stress. It positively regulates cell wall biosynthesis and root development via modulating ABA signals (Seifert et al. 2014). Huang et al. (2013) conducted a subcellular localization experiment for the cotton gene GhFLA1, and found the GFP fluorescence signals were displayed on the peripheral of transgenic GhFLA1(SP)-eGFPGhFLA1(GPI) Arabidopsis root cell. Further, they showed that this gene was involved in fiber initiation and elongation in upland cotton. In Barley, Stephan et al. (2016) found that the transcription factor WRKY29 was one of the most important candidate genes for root system variation, and a crucial amino acid substitution was detected within the conserved domain of WRKY29 among genotypes carrying significant and minor QTL alleles. WRKY29 was also closely associated with ABA signaling. A recent study found that OsWRKY29 could act as a new ABA signaling repressor for seed dormancy in rice with RNAi mutant lines displaying greater sensitivity after ABA treatment (Zhou et al. 2020). The WRKY29 gene was negatively correlated with the ABA content, and down-regulated in the tolerant rootstocks under osmotic stress (Hezema et al. 2021). RBOHA is a burst oxidase homolog protein, which participates in diverse biological processes in higher plants (Arthikala and Quinto 2018). Arthikala et al. (2018) found that the relative expression levels of PvRbohA were correlated with the activity of its promoter. The transgenic hairy roots of PvRbohA transcript silencing showed a significantly altered lateral root phenotype, which indicated that RbohA was involved in lateral root initiation, emergence, and development for P. Vulgaris. In rice, Groom et al. (2010) identified the complete cDNA and genomic DNA sequence of rbohA. They hypothesized that RbohA is a component of an NADPH oxidase that functions in defense-related active oxygen production on the root surface. A recent study in Brassica campestris showed that expression of RbohA and RbohD was induced by ABA and H2O2 (Zhang et al. 2019).
GWAS of seedling lateral root number for 215 G. arboreum accessions identified 18 significant signals and 19 candidate genes. Three genes, FLA12, WRKY29, and RBOHA, showed higher expression levels in root samples Shixiya1, which may be associated with lateral root development. These findings will be essential to decipher the genetics behind developmental mechanisms of lateral root formation in Asian cotton.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its additional files.
Genome-wide association studies
Lateral root number
Minor allele frequency for the genetic variant of interest
Single nucleotide polymorphism
- G. arboreum :
Gossypium arboreum L.
Transcripts per million
Boll weight in Anyang (Henan)
Boll weight in Sanya (Hainan)
Boll weight in Aksu (Xinjiang)
Sympodial branch node in Anyang (Henan)
Sympodial branch node in Aksu (Xinjiang)
Fusarium wilt disease index in Anyang (Henan)
Lint percentage in Anyang (Henan)
Lint percentage in Sanya (Hainan)
Lint percentage in Aksu (Xinjiang)
Seed index in Anyang (Henan)
Verticillium wilt disease index in Anyang (Henan)
Arthikala MK, Quinto C. RbohA coordinates lateral root emergence in common bean. Commun Integr Biol. 2018;11(2):1–10. https://doi.org/10.1080/19420889.2018.1467188.
Bradbury PJ, Zhang Z, Kroon DE, et al. TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics. 2007;23(19):2633–5. https://doi.org/10.1093/bioinformatics/btm308.
Cai X, Ping X, Yao W, et al. Activated expression of AtEDT1/HDG11 promotes lateral root formation in Arabidopsis mutant edt1 by upregulating jasmonate biosynthesis. J Integr Plant Biol. 2015;57(12):1017–30. https://doi.org/10.1111/jipb.12347.
Casero PJ, Casimiro I, Rodriguez-Gallardo L, et al. Lateral root initiation by asymmetrical transverse divisions of pericycle cells in adventitious roots of Allium cepa. Protoplasma. 1993;176:138–44. https://doi.org/10.1007/BF01378950.
Dai PH, Sun GF, Jia YH, et al. Extensive haplotypes are associated with population differentiation and environmental adaptability in Upland cotton (Gossypium hirsutum). Theor Appl Genet. 2020;133(12):3273–85. https://doi.org/10.1007/s00122-020-03668-z.
Du X, Huang G, He SP, et al. Resequencing of 243 diploid cotton accessions based on an updated A genome identifies the genetic basis of key agronomic traits. Nat Genet. 2018;50(6):796–802. https://doi.org/10.1038/s41588-018-0116-x.
Duan L, Dietrich D, Ng CH, et al. Endodermal ABA signaling promotes lateral root quiescence during salt stress in Arabidopsis seedlings. Plant Cell. 2013;25(1):324–41. https://doi.org/10.1105/tpc.112.107227.
Fang L, Wang Q, Hu Y, et al. Genomic analyses in cotton identify signatures of selection and loci associated with fiber quality and yield traits. Nat Genet. 2017;49:1089–98. https://doi.org/10.1038/ng.3887.
Farquharson KL. Gibberellin-auxin crosstalk modulates lateral root formation. Plant Cell. 2010;22(3):540. https://doi.org/10.2307/25680076.
Francis JL. GraphPad Prism(TM) Version 1.03. Blood Coagul Fibrinolysis. 1995;6(2):146. https://doi.org/10.1097/00001721-199504000-00014.
Fukaki H, Taniguchi N, Tasaka M. PICKLE is required for SOLITARY-ROOT/IAA14-mediated repression of ARF7 and ARF19 activity during Arabidopsis lateral root initiation. Plant J. 2010;48(3):380–9. https://doi.org/10.1111/j.1365-313X.2006.02882.x.
Groom QJ, Torres MA, Fordham-Skelton AP, et al. rbohA, a rice homologue of the mammalian gp91phox respiratory burst oxidase gene. Plant J. 2010;10(3):515–22. https://doi.org/10.1046/j.1365-313X.1996.10030515.x.
He SP, Sun GF, Geng XL, et al. The genomic basis of geographic differentiation and fiber improvement in cultivated cotton. Nat Genet. 2021;53:916–24. https://doi.org/10.1038/s41588-021-00844-9.
Helga T, Robinson J, Mesirov J. Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinf. 2013;14(2):178–92. https://doi.org/10.1093/bib/bbs017.
Hetherington FM, Kakkar M, Topping JF, et al. Gibberellin signaling mediates lateral root inhibition in response to K+-deprivation. Plant Physiol. 2021;185(3):1198–215. https://doi.org/10.1093/plphys/kiaa093.
Hezema YS, Shukla MR, Ayyanath MM, et al. Physiological and molecular responses of six apple rootstocks to osmotic stress. Int J Mol Sci. 2021;22(15):8263. https://doi.org/10.3390/ijms22158263.
Huang GQ, Gong SY, Xu WL, et al. A fasciclin-like arabinogalactan protein, GhFLA1, is involved in fiber initiation and elongation of cotton. Plant Physiol. 2013;161(3):1278. https://doi.org/10.1104/pp.112.203760.
Huang G, Wu Z, Percy RG, et al. Genome sequence of Gossypium herbaceum and genome updates of Gossypium arboreum and Gossypium hirsutum provide insights into cotton A-genome evolution. Nat Genet. 2020;52:516–24. https://doi.org/10.1038/s41588-020-0607-4.
Inahashi H, Shelley IJ, Yamauchi T, et al. OsPIN2, which encodes a member of the auxin efflux carrier proteins, is involved in root elongation growth and lateral root formation patterns via the regulation of auxin distribution in rice. Physiol Plant. 2018;164(2):216–25. https://doi.org/10.1111/ppl.12707.
Ivanchenko MG, Muday GK, Dubrovsky JG. Ethylene-auxin interactions regulate lateral root initiation and emergence in Arabidopsis thaliana. Plant J. 2008;55(2):335–47. https://doi.org/10.1111/j.1365-313X.2008.03528.x.
Ive DS. Multimodular auxin response controls lateral root development in Arabidopsis. Plant Signal Behav. 2014;5(5):2705–10. https://doi.org/10.4161/psb.11495.
Jia XC, Liu P, Lynch JP. Greater lateral root branching density in maize improves phosphorus acquisition from low phosphorus soil. J Exp Bot. 2018;69:4961–70. https://doi.org/10.1093/jxb/ery252.
Johnson KL, Kibble N, Bacic A, et al. A Fasciclin-Like Arabinogalactan-Protein (FLA) mutant of Arabidopsis thaliana, fla1, shows defects in shoot regeneration. PLoS One. 2011;6(9):e25154. https://doi.org/10.1371/journal.pone.0025154.
Kai W, Li M, Hakon H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38(16):e164. https://doi.org/10.1093/nar/gkq603.
Kang HM, Sul JH, Service SK, et al. Variance component model to account for sample structure in genome-wide association studies. Nat Genet. 2010;42(4):348–54. https://doi.org/10.1038/ng.548.
Kuo HY, Kang FC, Wang YY. Glucosinolate Transporter1 involves in salt-induced jasmonate signaling and alleviates the repression of lateral root growth by salt in Arabidopsis. Plant Sci. 2020;297:110487. https://doi.org/10.1016/j.plantsci.2020.110487.
Lewis DR, Negi S, Sukumar P, et al. Ethylene inhibits lateral root development, increases IAA transport and expression of PIN3 and PIN7 auxin efflux carriers. Development. 2011;138(16):3485–95. https://doi.org/10.1242/dev.065102.
Li H, Durbin R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics. 2009;25:1754–60. https://doi.org/10.1093/bioinformatics/btp324.
Li Z, Fan W, Liu Y. Correlations between lint yield and root and plant main characters of upland cotton. Acta Agric Boreali-Occidentalis Sin. 2003;1:36–9.
Li MX, Yeung JMY, Cherny SS, Sham PC. Evaluating the effective numbers of independent tests and significant p-value thresholds in commercial genotyping arrays and public imputation reference datasets. Hum Genet. 2012;131(5):747–56. https://doi.org/10.1007/s00439-011-1118-2.
Li F, Fan G, Wang K, et al. Genome sequence of the cultivated cotton Gossypium arboreum. Nat Genet. 2014;46(6):567–72. https://doi.org/10.1038/ng.2987.
Li T, Ma X, Li N, et al. Genome-wide association study discovered candidate genes of Verticillium wilt resistance in upland cotton (Gossypium hirsutum L.). Plant Biotechnol J. 2017;15(12):1520–32. https://doi.org/10.1111/pbi.12734.
Liang HZ, Dong W, Xu LJ, et al. QTL mapping for main root length and lateral root number in soybean at the seedling stage in different N, P and K environments. Sci Agric Sin. 2017;50(18):3450–60. https://doi.org/10.3864/j.issn.0578-1752.2017.18.002.
Lv BS, Wei KJ, Hu KQ, et al. MPK14-mediated auxin signaling controls lateral root development via ERF13-regulated very-long-chain fatty acid biosynthesis. Mol Plant. 2021;14:285–97. https://doi.org/10.1016/j.molp.2020.11.011.
Marhavy P, Vanstraelen M, Rybel BD, et al. Auxin reflux between the endodermis and pericycle promotes lateral root initiation. EMBO J. 2014;32(1):149–58. https://doi.org/10.1038/emboj.2012.303.
Mckenna A, Hanna M, Banks E, et al. The Genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297–303. https://doi.org/10.1101/gr.107524.110.
Niones JM, Inukai Y, Suralta RR, et al. QTL associated with lateral root plasticity in response to soil moisture fluctuation stress in rice. Plant Soil. 2015;391:63–75. https://doi.org/10.1007/s11104-015-2404-x.
Perianez-Rodriguez J, Rodriguez M, Marconi M, et al. An auxin-regulable oscillatory circuit drives the root clock in Arabidopsis. Sci Adv. 2021;7(1):eabd4722. https://doi.org/10.1126/sciadv.abd4722.
Prasad ME, Schofield A, Lyzenga W, et al. Arabidopsis RING E3 Ligase XBAT32 regulates lateral root production through its role in ethylene biosynthesis. Plant Physiol. 2010;153(4):1587–96. https://doi.org/10.4161/psb.5.11.13294.
R Core Team. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2012. http://www.R-project.org/. Accessed 15 Nov 2021.
Seifert GJ, Xue H, Tuba A. The Arabidopsis thaliana FASCICLIN LIKE ARABINOGALACTAN PROTEIN 4 gene acts synergistically with abscisic acid signalling to control root growth. Ann Bot. 2014;114(6):1125–33. https://doi.org/10.1093/aob/mcu010.
Stephan R, Annika K, Jens L, et al. Genome-wide association mapping in the global diversity set reveals new QTL controlling root system and related shoot variation in barley. Front Plant Sci. 2016;7:1061–74. https://doi.org/10.3389/fpls.2016.01061.
Talboys PJ, Healey JR, Withers P, et al. Phosphate depletion modulates auxin transport in Triticum aestivum leading to altered root branching. J Exp Bot. 2014;65(17):5023–32. https://doi.org/10.1093/jxb/eru284.
Turner SD. qqman: an R package for visualizing GWAS results using Q-Q and manhattan plots. J Open Source Softw. 2018;3(25):731. https://doi.org/10.1101/005165.
Wang H, Wei J, Li P, et al. Integrating GWAS and gene expression analysis identifies candidate genes for root morphology traits in maize at the seedling stage. Genes (Basel). 2019;10(10):773. https://doi.org/10.3390/genes10100773.
Wu Q, Du MW, Wu J, et al. Mepiquat chloride promotes cotton lateral root formation by modulating plant hormone homeostasis. BMC Plant Biol. 2019;19(1):573. https://doi.org/10.1186/s12870-019-2176-1.
Yang Z, Ge X, Yang Z, et al. Extensive intraspecific gene order and gene structural variations in upland cotton cultivars. Nat Commun. 2019;10(1):1–13. https://doi.org/10.1038/s41467-019-10820-x.
Ying Z, Thomas CL, Xiang J, et al. QTL meta-analysis of root traits in Brassica napus under contrasting phosphorus supply in two growth systems. Sci Rep. 2016;6(1):33113. https://doi.org/10.1038/srep33113.
Yu E, Naoki Y, Mao C, et al. Lateral roots but not root hairs contribute to high uptake of manganese and cadmium in rice. J Exp Bot. 2021;5:1–9. https://doi.org/10.1093/jxb/erab329.
Zang L, Zheng T, Chu Y, et al. Genome-wide analysis of the fasciclin-like arabinogalactan protein gene family reveals differential expression patterns, localization, and salt stress response in Populus. Front Plant Sci. 2015;6(65):1140. https://doi.org/10.3389/fpls.2015.01140.
Zhang TG, Lai J, Ping LI, et al. Expression analysis of RbohA and RbohD genes in Brassica campestris under different treatments. Chin J Plant Ecol. 2019;38(1):173–80. https://doi.org/10.13292/j.1000-4890.201901.024.
Zhao Y, Xing L, Wang X, et al. The ABA receptor PYL8 promotes lateral root growth by enhancing MYB77-dependent transcription of auxin-responsive genes. Sci Signal. 2014;7(328):ra53. https://doi.org/10.1126/scisignal.2005051.
Zhao Z, Hu D, Azhar MT, et al. Genome-wide association and transcriptome analysis of root color-related genes in Gossypium arboreum L. Planta. 2021;253(5):95. https://doi.org/10.1007/s00425-021-03622-3.
Zhou C, Lin Q, Lan J, et al. WRKY transcription factor OsWRKY29 represses seed dormancy in rice by weakening abscisic acid response. Front Plant Sci. 2020;11:691. https://doi.org/10.3389/fpls.2020.00691.
The authors would like to thank the National Medium-term Gene Bank of Cotton in China for their kindly providing cotton seeds.
This work was supported by Central Public-interest Scientific Institution Basal Research Fund, Chinese Academy of Agricultural Science (No.1610162021012) funded by DU Xiongming.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Table S1 Phenotypic data of lateral root number in 215 Gossypium arboreum accessions; Table S2 Basic description analysis of lateral root number in 215 Gossypium arboreum accessions; Table S3 18 significant SNPs of lateral root numbers in 215 Gossypium arboreum accessions; Table S4 19 candidate genes of lateral root number identified by the GWAS method; Table S5 TPM expression of candidate genes in different tissues
Fig. S1 SNP density plot
Fig. S2 Quantile - quantile (QQ) plots of lateral root number in Asian cotton
About this article
Cite this article
HU, D., HE, S., SUN, G. et al. A genome-wide association study of lateral root number for Asian cotton (Gossypium arboreum L.). J Cotton Res 5, 19 (2022). https://doi.org/10.1186/s42397-022-00126-7