Abstract
Optimal root system architecture (RSA) is essential for vigorous growth and yield in crops. Plants have evolved adaptive mechanisms in response to low phosphorus (LP) stress, and one of those is changes in RSA. Here, more than five million single-nucleotide polymorphisms (SNPs) obtained from whole-genome re-sequencing data (WGR) of an association panel of 370 oilseed rape (Brassica napus L.) were used to conduct a genome-wide association study (GWAS) of RSA traits of the panel at LP in “pouch and wick” system. Fifty-two SNPs were forcefully associated with lateral root length (LRL), total root length (TRL), lateral root density (LRD), lateral root number (LRN), mean lateral root length (MLRL), and root dry weight (RDW) at LP. There were significant correlations between phenotypic variation and the number of favorable alleles of the associated loci on chromosomes A06 (chrA06_20030601), C03 (chrC03_3535483), and C07 (chrC07_42348561), respectively. Three candidate genes (BnaA06g29270D, BnaC03g07130D, and BnaC07g43230D) were detected by combining transcriptome, candidate gene association analysis, and haplotype analysis. Cultivar carrying “CCGC” at BnaA06g29270DHap1, “CAAT” at BnaC03g07130DHap1, and “ATC” at BnaC07g43230DHap1 had greater LRL, LRN, and RDW than lines carrying other haplotypes at LP supply. The RSA of a cultivar harboring the three favorable haplotypes was further confirmed by solution culture experiments. These findings define exquisite insights into genetic architectures underlying B. napus RSA at LP and provide valuable gene resources for root breeding.
Similar content being viewed by others
Data availability
The original contributions presented in this study are included in the article/Supplementary Material, and further inquiries can be directed to the corresponding author. Raw sequencing data of genome re-sequencing are available in the Genome Sequence Archive (https://bigd.big.ac.cn/gsa/) with Bio-project IDs PRJCA002835 and PRJCA002836. All the materials in this study are available upon request.
References
Alexander DH, Novembre J, Lange K (2009) Fast model-based estimation of ancestry in unrelated individuals. Genome Res 19(9):1655–1664. https://doi.org/10.1101/gr.094052.109
Arifuzzaman M, Oladzadabbasabadi A, McClean P et al (2019) Shovelomics for phenotyping root architectural traits of rapeseed/canola (Brassica napus L.) and genome-wide association mapping. Mol Genet Genomics 294(4):985–1000. https://doi.org/10.1007/s00438-019-01563-x
Barrett JC, Fry B, Maller J et al (2005) Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 21(2):263–265. https://doi.org/10.1093/bioinformatics/bth457
Bradbury PJ, Zhang Z, Kroon DE et al (2007) TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics 23(19):2633–2635. https://doi.org/10.1093/bioinformatics/btm308
Chalhoub B, Denoeud F, Liu SY et al (2014) Early allopolyploid evolution in the post-Neolithic Brassica napus oilseed genome. Science 345(6199):950–953. https://doi.org/10.1126/science.1253435
Clark RT, Famoso AN, Zhao K et al (2013) High-throughput two-dimensional root system phenotyping platform facilitates genetic analysis of root growth and development. Plant Cell Environ 36(2):454–466. https://doi.org/10.1111/j.1365-3040.2012.02587.x
Danecek P, Auton A, Abecasis G et al (2011) The variant call format and VCFtools. Bioinformatics 27(15):2156–2158. https://doi.org/10.1093/bioinformatics/btr330
Dong HL, Tan CD, Li YZ et al (2018) Genome-wide association study reveals both overlapping and independent genetic loci to control seed weight and silique length in Brassica napus. Front Plant Sci 9:921. https://doi.org/10.3389/fpls.2018.00921
Dong SS, He WM, Ji JJ et al (2021) LDBlockShow: a fast and convenient tool for visualizing linkage disequilibrium and haplotype blocks based on variant call format files. Brief Bioinform 22(4). https://doi.org/10.1093/bib/bbaa227
Du HY, Yang C, Ding GD (2017) Genome-wide identification and characterization of SPX domain-containing members and their responses to phosphate deficiency in Brassica napus. Front Plant Sci 8. https://doi.org/10.3389/fpls.2017.00035
Duan XJ, Wang XH, Jin KM et al (2021) Genetic dissection of root angle of Brassica napus in response to low phosphorus. Front Plant Sci 12. https://doi.org/10.3389/Fpls.2021.697872
Excoffier L, Lischer HE (2010) Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour 10(3):564–567. https://doi.org/10.1111/j.1755-0998.2010.02847.x
Fedosejevs ET, Liu LNC, Abergel M et al (2017) Coimmunoprecipitation of reversibly glycosylated polypeptide with sucrose synthase from developing castor oilseeds. Febs Lett 591(23):3872–3880. https://doi.org/10.1002/1873-3468.12893
He ML, Wang SL, Zhang C et al (2021) Genetic variation of BnaA3.NIP5;1 expressing in the lateral root cap contributes to boron deficiency tolerance in Brassica napus. PLoS Genet 17(7). https://doi.org/10.1371/journal.pgen.1009661
He YJ, Hu DX, You JC et al (2019) Genome-wide association study and protein network analysis for understanding candidate genes involved in root development at the rapeseed seedling stage. Plant Physiol Biochem 137:42–52. https://doi.org/10.1016/j.plaphy.2019.01.028
Huang LY, Min Y, Schiessl S et al (2021) Integrative analysis of GWAS and transcriptome to reveal novel loci regulation flowering time in semi-winter rapeseed. Plant Sci 310. https://doi.org/10.1016/j.plantsci.2021.110980
Jan HU, Guan M, Yao M et al (2019) Genome-wide haplotype analysis improves trait predictions in Brassica napus hybrids. Plant Sci 283:157–164. https://doi.org/10.1016/j.plantsci.2019.02.007
Kittipol V, He ZS, Wang LH et al (2019) Genetic architecture of glucosinolate variation in Brassica napus. J Plant Physiol 240. https://doi.org/10.1016/J.Jplph.2019.06.001
Korte A, Farlow A (2013) The advantages and limitations of trait analysis with GWAS: a review. Plant Methods 9. https://doi.org/10.1186/1746-4811-9-29
Li F, Chen B, Xu K et al (2014) Genome-wide association study dissects the genetic architecture of seed weight and seed quality in rapeseed (Brassica napus L.). DNA Res 21(4):355–367. https://doi.org/10.1093/dnares/dsu002
Li KQ, Wang J, Kuang LQ et al (2021) Genome-wide association study and transcriptome analysis reveal key genes affecting root growth dynamics in rapeseed. Biotechnol Biofuels 14(1). https://doi.org/10.1186/s13068-021-02032-7
Li Y, Wang X, Zhang H et al (2019) Molecular identification of the phosphate transporter family 1 (PHT1) genes and their expression profiles in response to phosphorus deprivation and other abiotic stresses in Brassica napus. PLoS One 14(7). https://doi.org/10.1371/journal.pone.0220374
Liu HJ, Wang JC, Zhang BB et al (2021a) Genome-wide association study dissects the genetic control of plant height and branch number in response to low-phosphorus stress in Brassica napus. Ann Bot 128(7):919–929. https://doi.org/10.1093/aob/mcab115
Liu HJ, Li XJ, Zhang QW et al (2021b) Integrating a genome-wide association study with transcriptomic data to predict candidate genes and favourable haplotypes influencing Brassica napus seed phytate. DNA Res 28(5). https://doi.org/10.1093/dnares/dsab011
Liu HJ, Wang W, Yang M et al (2022) Genome-wide association studies of important agronomic traits in Brassica napus: what we have learned and where we are headed. Annual Plant Reviews 5:1–30
Liu YX, Wang L, Deng M et al (2015) Genome-wide association study of phosphorus-deficiency-tolerance traits in Aegilops tauschii. Theor Appl Genet 128(11):2203–2212. https://doi.org/10.1007/s00122-015-2578-x
Lu K, Peng L, Zhang C et al (2017) Genome-wide association and transcriptome analyses reveal candidate genes underlying yield-determining traits in Brassica napus. Front Plant Sci 8. https://doi.org/10.3389/Fpls.2017.00206
Luo T, Zhang YT, Zhang CN et al (2021) Genome-wide association mapping unravels the genetic control of seed vigor under low-temperature conditions in rapeseed (Brassica napus L.). Plants (Basel) 10:(3). https://doi.org/10.3390/Plants10030426
Lynch J (1995) Root architecture and plant productivity. Plant Physiol 109(1):7–13. https://doi.org/10.1104/Pp.109.1.7
Mai HJ, Lindermayr C, von Toerne C et al (2015) Iron and FER-like iron deficiency-induced transcription factor-dependent regulation of proteins and genes in Arabidopsis thaliana roots. Proteomics 15(17):3030–3047. https://doi.org/10.1002/pmic.201400351
Müller J, Toev T, Heisters M et al (2015) Iron-dependent callose deposition adjusts root meristem maintenance to phosphate availability. Dev Cell 33(2):216–230. https://doi.org/10.1016/j.devcel.2015.02.007
Muralidhara P, Weiste C, Collani S et al (2021) Perturbations in plant energy homeostasis prime lateral root initiation via SnRK1-bZIP63-ARF19 signaling. Proc Natl Acad Sci USA 118(37). https://doi.org/10.1073/pnas.2106961118
Raman H, Raman R, Qiu Y et al (2019) GWAS hints at pleiotropic roles for FLOWERING LOCUS T in flowering time and yield-related traits in canola. BMC Genomics 20(1):636. https://doi.org/10.1186/s12864-019-5964-y
Shi L, Shi TX, Broadley MR et al (2013) High-throughput root phenotyping screens identify genetic loci associated with root architectural traits in Brassica napus under contrasting phosphate availabilities. Ann Bot 112(2):381–389. https://doi.org/10.1093/aob/mcs245
Singh AP, Fridman Y, Holland N et al (2018) Interdependent nutrient availability and steroid hormone signals facilitate root growth plasticity. Dev Cell 46(1):59–72.e4. https://doi.org/10.1016/j.devcel.2018.06.002
Tan ZD, Xie ZQ, Dai LH et al (2022) Genome- and transcriptome-wide association studies reveal the genetic basis and the breeding history of seed glucosinolate content in Brassica napus. Plant Biotechnol J 20(1):211–225. https://doi.org/10.1111/pbi.13707
Tang S, Zhao H, Lu SP et al (2021) Genome- and transcriptome-wide association studies provide insights into the genetic basis of natural variation of seed oil content in Brassica napus. Mol Plant 14(3):470–487. https://doi.org/10.1016/j.molp.2020.12.003
Ticconi CA, Lucero RD, Sakhonwasee S et al (2009) ER-resident proteins PDR2 and LPR1 mediate the developmental response of root meristems to phosphate availability. Proc Natl Acad Sci U S A. 106(33):14174–14179. https://doi.org/10.1073/pnas.0901778106
Vance CP, Uhde-Stone C, Allan DL (2003) Phosphorus acquisition and use: critical adaptations by plants for securing a nonrenewable resource. New Phytol 157(3):423–447. https://doi.org/10.1046/j.1469-8137.2003.00695.x
Ward JT, Lahner B, Yakubova E, Salt DE, Raghothama KG (2008) The effect of iron on the primary root elongation of Arabidopsis during phosphate deficiency. Plant Physiol 147:1181–1191. https://doi.org/10.1104/pp.108.118562
Wang XH, Chen YL, Thomas CL et al (2017) Genetic variants associated with the root system architecture of oilseed rape (Brassica napus L.) under contrasting phosphate supply. DNA Res 24(4):407–417. https://doi.org/10.1093/dnares/dsx013
Wu D, Liang Z, Yan T et al (2019) Whole-genome resequencing of a worldwide collection of rapeseed accessions reveals the genetic basis of ecotype divergence. Mol Plant 12(1):30–43. https://doi.org/10.1016/j.molp.2018.11.007
Xiao YJ, Liu HJ, Wu LJ et al (2017) Genome-wide association studies in maize: praise and stargaze. Mol Plant 10(3):359–374. https://doi.org/10.1016/j.molp.2016.12.008
Zhang C, Dong SS, Xu JY et al (2019) PopLDdecay: a fast and effective tool for linkage disequilibrium decay analysis based on variant call format files. Bioinformatics 35(10):1786–1788. https://doi.org/10.1093/bioinformatics/bty875
Zhang Y, Thomas CL, Xiang JX et al (2016) QTL meta-analysis of root traits in Brassica napus under contrasting phosphorus supply in two growth systems. Sci Rep 6. https://doi.org/10.1038/Srep33113
Funding
This work was supported by the National Nature Science Foundation of China (Grant Nos. 31972498 and 32172662). The computations in this paper were run on the bioinformatics computing platform of the National Key Laboratory of Crop Genetic Improvement, Huazhong Agricultural University.
Author information
Authors and Affiliations
Contributions
Pan Yuan, Haijiang Liu, and Lei Shi designed the research, reviewed the writing, and drafted the manuscript. Pan Yuan, Haijiang Liu, and Xiaohua Wang participated the experiments. John P. Hammond participated in the manuscript revision.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
ESM 1
Fig. S1. Correlation of eight root related traits at low phosphorus supplies. (DOCX 257 kb)
ESM 2
Fig. S2. The LD decay of an association panel of B. napus. (DOCX 45 kb)
ESM 3
Fig. S3. Population structure of an association panel of B. napus with K from 2 to 8. (DOCX 135 kb)
ESM 4
Fig. S4. The kinship of an association panel of 370 B. napus accessions. (DOCX 175 kb)
ESM 5
Fig. S5. Distribution of linkage disequilibrium block sizes across all chromosomes. (DOCX 63 kb)
ESM 6
Table S1. List of 370 accessions of B. napus used in the study. (XLSX 20 kb)
ESM 7
Table S2. Primers used for qRT-PCR. (XLSX 9 kb)
ESM 8
Table S3. Root related traits at low phosphorus supplies in an association panel of B. napus. (XLSX 30 kb)
ESM 9
Table S4. AMOVA analysis between the K = 4 assigned populations in Brassica napus. (XLSX 9 kb)
ESM 10
Table S5. Linkage disequilibrium block in this study. (XLSX 28730 kb)
ESM 11
Table S6. Significant SNP loci for root related traits of B. napus by genome wide association study at low phosphorus supplies. (XLSX 12 kb)
ESM 12
Table S7. LD blocks harboring significant SNPs associated with RSA. (XLSX 12 kb)
ESM 13
Table S8. Candidate genes within LD decay value up and down the lead SNPs (chrA06_19934701, chrC03_3535476 and chrC07_42348526) for root related traits. (XLSX 43 kb)
ESM 14
Table S9. Comparison of SNPs detected by WGR for RSA in this study with previously identified SNPs by 60 K SNP chip for RSA at a low phosphorus supply. (XLSX 10 kb)
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Yuan, P., Liu, H., Wang, X. et al. Genome-wide association study reveals candidate genes controlling root system architecture under low phosphorus supply at seedling stage in Brassica napus. Mol Breeding 43, 63 (2023). https://doi.org/10.1007/s11032-023-01411-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s11032-023-01411-2