Genetic variants associated with the root system architecture of oilseed rape (Brassica napus L.) under contrasting phosphate supply

Genome-wide association study of Verticillium longisporum resistance in Brassica genotypes

Verticillium stripe, caused by Verticillium longisporum, presents an emerging threat to Canadian canola (Brassica napus). Initially detected in Manitoba in 2014, the presence of this pathogen has since been confirmed across western Canada. Infections by V. longisporum can result in yield losses of up to 50%, which is a cause for concern given the susceptibility of most commercial Canadian canola cultivars. The objective of this study was to screen a collection of 211 Brassica genotypes for their reactions to V. longisporum, and to use genome-wide association study (GWAS) to identify single nucleotide polymorphism (SNP) markers for resistance. The plant material consisted of 110 rutabaga (B. napus ssp. napobrassica), 35 canola, 40 Brassica rapa, and 15 Brassica oleracea accessions or cultivars, alongside 11 hosts of the European Clubroot Differential (ECD) set. These materials were screened for resistance under greenhouse conditions and were genotyped using a 19K Brassica SNP array. Three general linear models (GLM), four mixed linear models (MLM), and three GWAS methods were employed to evaluate the markers. Eleven non-commercial Brassica accessions and 9 out of 35 commercial canola cultivars displayed a low normalized area under the disease progress curve (AUDPCnorm.). The non-commercial accessions could prove valuable as potential sources of resistance against V. longisporum. Forty-five SNP markers were identified to be significantly associated with V. longisporum resistance using single-SNP based GWAS analysis. In comparison, haplotype-based GWAS analyses identified 10 to 25 haplotype blocks to be significantly associated with V. longisporum resistance. Between 20% and 56% of QTLs identified by the more conventional single-SNP based GWAS analysis were also detected by the haplotype-based GWAS analysis. The overlapping genomic regions identified by the two GWAS methods present promising hotspots for marker-assisted selection in the future development of Verticillium stripe-resistant canola.

Genetic control of root morphology in response to nitrogen across rapeseed diversity

Rapeseed (Brassica napus L.) is an oil-containing crop of great economic value but with considerable nitrogen requirement. Breeding root systems that efficiently absorb nitrogen from the soil could be a driver to ensure genetic gains for more sustainable rapeseed production. The aim of this study is to identify genomic regions that regulate root morphology in response to nitrate availability. The natural variability offered by 300 inbred lines was screened at two experimental locations. Seedlings grew hydroponically with low or elevated nitrate levels. Fifteen traits related to biomass production and root morphology were measured. On average across the panel, a low nitrate level increased the root-to-shoot biomass ratio and the lateral root length. A large phenotypic variation was observed, along with important heritability values and genotypic effects, but low genotype-by-nitrogen interactions. Genome-wide association study and bulk segregant analysis were used to identify loci regulating phenotypic traits. The first approach nominated 319 SNPs that were combined into 80 QTLs. Three QTLs identified on the A07 and C07 chromosomes were stable across nitrate levels and/or experimental locations. The second approach involved genotyping two groups of individuals from an experimental F2 population created by crossing two accessions with contrasting lateral root lengths. These individuals were found in the tails of the phenotypic distribution. Co-localized QTLs found in both mapping approaches covered a chromosomal region on the A06 chromosome. The QTL regions contained some genes putatively involved in root organogenesis and represent selection targets for redesigning the root morphology of rapeseed.

Trade-offs between root-secreted acid phosphatase and root morphology traits, and their contribution to phosphorus acquisition in Brassica napus

Oilseed rape (Brassica napus) is one of the most important oil crops in the world and shows sensitivity to low phosphorus (P) availability. In many soils, organic P (Po) is the main component of the soil P pool. Po must be mineralised to Pi through phosphatases, and then taken up by plants. However, the relationship between root-secreted acid phosphatases (APase) and root morphology traits, two important P-acquisition strategies in response to P deficiency, is unclear among B. napus genotypes. This study aimed to understand their relationship and how they affect P acquisition, which is crucial for the sustainable utilisation of agricultural P resources. This study showed significant genotypic variations in root-secreted APase activity per unit root fresh weight (SAP) and total root-secreted APase activity per plant (total SAP) among 350 B. napus genotypes. Seed yield was positively correlated with total SAP but not significantly correlated with SAP. Six root traits of 18 B. napus genotypes with contrasting root biomass were compared under normal Pi, low Pi and Po. Genotypes with longer total root length (TRL) reduced SAP, but those with shorter TRL increased SAP under P deficiency. Additionally, TRL was important in P-acquisition under three P treatments, and total SAP was also important in P-acquisition under Po treatment. In conclusion, trade-offs existed between the two P-acquisition strategies among B. napus genotypes under P-deficient conditions. Total SAP was an important root trait under Po conditions. These results might help to breed B. napus with greater P-acquisition ability under low P availability conditions.

Genome-wide association study reveals candidate genes controlling root system architecture under low phosphorus supply at seedling stage in Brassica napus

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.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01411-2.

Integrating genome-wide association studies with selective sweep reveals genetic loci associated with tolerance to low phosphate availability in Brassica napus

Oilseed rape (Brassica napus L.; B. napus) is an important oil crop worldwide. However, the genetic mechanisms of B. napus adaptations to low phosphate (P) stress are largely unknown. In this study, a genome-wide association study (GWAS) identified 68 SNPs significantly associated with seed yield (SY) under low P (LP) availability, and 7 SNPs significantly associated with phosphorus efficiency coefficient (PEC) in two trials. Among these SNPs, two, chrC07__39807169 and chrC09__14194798, were co-detected in two trials, and BnaC07.ARF9 and BnaC09.PHT1;2 were identified as candidate genes of them, respectively, by combining GWAS with quantitative reverse-transcription PCR (qRT-PCR). There were significant differences in the gene expression level of BnaC07.ARF9 and BnaC09.PHT1;2 between P-efficient and -inefficiency varieties at LP. SY_LP had a significant positive correlation with the gene expression level of both BnaC07.ARF9 and BnaC09.PHT1;2. BnaC07.ARF9 and BnaA01.PHR1 could directly bind the promoters of BnaA01.PHR1 and BnaC09.PHT1;2, respectively. Selective sweep analysis was conducted between ancient and derived B. napus, and detected 1280 putative selective signals. Within the selected region, a large number of genes related to P uptake, transport, and utilization were detected, such as purple acid phosphatase (PAP) family genes and phosphate transporter (PHT) family genes. These findings provide novel insights into the molecular targets for breeding P efficiency varieties in B. napus.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01399-9.

Uncovering genetic control of primary root length variation in Brassica napus using QTL-seq. A commentary on: 'Rapid identification of a major locus qPRL-C06 affecting primary root length in Brassica napus by QTL-seq'

This article comments on:

Wei Wang, Haijiang Liu, Yiwen Xie, Graham John King, Philip John White, Jun Zou, Fangsen Xu, and Lei Shi. Rapid identification of a major locus qPRL-C06 affecting primary root length in Brassica napus by QTL-seq, Annals of Botany, Volume 131, Issue 4, 14 March 2023, Pages 569–583, https://doi.org/10.1093/aob/mcac123

Rapid identification of a major locus qPRL-C06 affecting primary root length in Brassica napus by QTL-seq

BACKGROUND AND AIMS

Brassica napus is one of the most important oilseed crops worldwide. Seed yield of B. napus significantly correlates with the primary root length (PRL). The aims of this study were to identify quantitative trait loci (QTLs) for PRL in B. napus.

METHODS

QTL-seq and conventional QTL mapping were jointly used to detect QTLs associated with PRL in a B. napus double haploid (DH) population derived from a cross between 'Tapidor' and 'Ningyou 7'. The identified major locus was confirmed and resolved by an association panel of B. napus and an advanced backcross population. RNA-seq analysis of two long-PRL lines (Tapidor and TN20) and two short-PRL lines (Ningyou 7 and TN77) was performed to identify differentially expressed genes in the primary root underlying the target QTLs.

KEY RESULTS

A total of 20 QTLs impacting PRL in B. napus grown at a low phosphorus (P) supply were found by QTL-seq. Eight out of ten QTLs affecting PRL at a low P supply discovered by conventional QTL mapping could be detected by QTL-seq. The locus qPRL-C06 identified by QTL-seq was repeatedly detected at both an optimal P supply and a low P supply by conventional QTL mapping. This major constitutive QTL was further confirmed by regional association mapping. qPRL-C06 was delimited to a 0.77 Mb genomic region on chromosome C06 using an advanced backcross population. A total of 36 candidate genes within qPRL-C06 were identified that showed variations in coding sequences and/or exhibited significant differences in mRNA abundances in primary root between the long-PRL and short-PRL lines, including five genes involved in phytohormone biosynthesis and signaling.

CONCLUSIONS

These results both demonstrate the power of the QTL-seq in rapid QTL detection for root traits and will contribute to marker-assisted selective breeding of B. napus cultivars with increased PRL.

Genetic variation in BnGRP1 contributes to low phosphorus tolerance in Brassica napus

Lack of phosphorus (P) is a major environmental factor affecting rapeseed (Brassica napus. L) root growth and development. For breeding purposes, it is crucial to identify the molecular mechanisms of root system architecture (RSA) traits underlying low P tolerance in rapeseed. The natural variations in the glycine-rich protein gene, BnGRP1, were analyzed in the natural population of 400 rapeseed cultivars under low P stress through genome-wide association study (GWAS) and transcriptome analyses. Based on 11 SNP mutations in BnGRP1 sequence, ten types of haplotypes (Hap) were formed. Compared with the other types, the cultivar of BnGRP1Hap1 type in the panel demonstrated the longest root length and heaviest root weight. BnGRP1Hap1 overexpression in rapeseed depicted the ability to enhance its resistance in response to low P tolerance. CRISPR/Cas9-derived BnGRP1Hap4 knockout mutations in rapeseed can lead to sensitivity to low P stress. Furthermore, BnGRP1Hap1 influences the expression of phosphate transporter 1 (PHT1) genes associated with P absorption. Overall, the findings of this study highlight new mechanisms of GRP1 genes in enhancing low P tolerance in rapeseed.

Genome wide association studies for acid phosphatase activity at varying phosphorous levels in Brassica juncea L

Acid phosphatases (Apases) are an important group of enzymes that hydrolyze soil and plant phosphoesters and anhydrides to release Pi (inorganic phosphate) for plant acquisition. Their activity is strongly correlated to the phosphorus use efficiency (PUE) of plants. Indian mustard (Brassica juncea L. Czern & Coss) is a major oilseed crop that also provides protein for the animal feed industry. It exhibits low PUE. Understanding the genetics of PUE and its component traits, especially Apase activity, will help to reduce Pi fertilizer application in the crop. In the present study, we evaluated 280 genotypes of the diversity fixed foundation set of Indian mustard for Apase activity in the root (RApase) and leaf (LApase) tissues at three- low (5µM), normal (250µM) and high (1mM) Pi levels in a hydroponic system. Substantial effects of genotype and Pi level were observed for Apase activity in both tissues of the evaluated lines. Low Pi stress induced higher mean RApase and LApase activities. However, mean LApase activity was relatively more than mean RApase at all three Pi levels. JM06016, IM70 and Kranti were identified as promising genotypes with higher LApase activity and increased R/S at low Pi. Genome-wide association study revealed 10 and 4 genomic regions associated with RApase and LApase, respectively. Annotation of genomic regions in the vicinity of peak associated SNPs allowed prediction of 15 candidates, including genes encoding different family members of the acid phosphatase such as PAP10 (purple acid phosphatase 10), PAP16, PNP (polynucleotide phosphorylase) and AT5G51260 (HAD superfamily gene, subfamily IIIB acid phosphatase) genes. Our studies provide an understanding of molecular mechanism of the Apase response of B. juncea at varying Pi levels. The identified SNPs and candidate genes will support marker-assisted breeding program for improving PUE in Indian mustard. This will redeem the crop with enhanced productivity under restricted Pi reserves and degrading agro-environments.

Genome-wide association studies reveal putative QTLs for physiological traits under contrasting phosphorous conditions in wheat (Triticum aestivum L.)

A Genome-wide association (GWAS) study was conducted for phosphorous (P)-use responsive physiological traits in bread wheat at the seedling stage under contrasting P regimes. A panel of 158 diverse advanced breeding lines and released varieties, and a set of 10,800 filtered single nucleotide polymorphism (SNP) markers were used to study marker-trait associations over the eight shoot traits. Principle component analysis separated the two environments (P regimes) because of the differential response of the traits indicating the essentiality of the separate breeding programmes for each environment. Significant variations for genotypic, environmental, and genotype × environment (GEI) effects were observed for all the traits in the combined analysis of variance with moderately high broad sense heritability traits (0.50–0.73). With the different algorithms of association mapping viz., BLINK, FarmCPU, and MLM, 38 unique QTLs under non-limiting P (NLP) and 45 QTLs for limiting P (LP) conditions for various shoot traits were identified. Some of these QTLs were captured by all three algorithms. Interestingly, a Q.iari.dt.sdw.1 on chromosome 1D was found to explain the significant variations in three important physiological traits under non-limiting phosphorus (NLP) conditions. We identified the putative candidate genes for QTLs namely Q.iari.dt.chl.1, Q.iari.dt.sdw.16, Q.iari.dt.sdw.9 and Q.iari.dt.tpc.1 which are potentially involved in the mechanism regulating phosphorus use efficiency through improved P absorption due to improved root architectural traits and better mobilization such as sulfotransferase involved in postembryonic root development, WALLS ARE THIN1 (WAT1), a plant-specific protein that facilitates auxin export; lectin receptor-like kinase essentially involved in plant development, stress response during germination and lateral root development and F-box component of the SKP-Cullin-F box E3 ubiquitin ligase complex and strigolactone signal perception. Expression profiling of putative genes located in identified genomic regions against the wheat expression atlas revealed their significance based on the expression of these genes for stress response and growth development processes in wheat. Our results thus provide an important insight into understanding the genetic basis for improving PUE under phosphorus stress conditions and can shape the future breeding programme by developing and integrating molecular markers for these difficult-to-score important traits.

Integrated analysis of the lncRNA/circRNA-miRNA-mRNA expression profiles reveals novel insights into potential mechanisms in response to root-knot nematodes in peanut

Background

Peanut is the most essential oil and food crop globally due to its high oil and protein content. Root-knot nematode infects peanut roots, causing poor development and severely limiting peanut yields worldwide. The discovery of peanut genome identified a considerable number of genetic loci controlling the peanut root-knot nematode; however, the molecular mechanism of root-knot nematode remains unknown.

Results

The heterogeneous response to root-knot nematode stress in peanut roots was identified using whole-transcriptome RNA-seq. A total of 430 mRNAs, 111 miRNAs, 4453 lncRNAs, and 123 circRNAs were found to have differential expression between infected and non-infected peanuts. The expression profiles of the lncRNA/circRNA-miRNA-mRNA network were developed to understand the potential pathways that lead to root-knot nematodes in peanut roots. During root-knot nematodes stress, a total of 10 lncRNAs, 4 circRNAs, 5 miRNAs, and 13 mRNAs can create competing endogenous RNA and participate in the oxidation–reduction process as well as other biological metabolism processes in peanuts. The findings will highlight the role of peanut ceRNAs in response to root-knot nematodes.

Conclusion

The GO classification and KEGG pathway enrichment study of core regulatory networks revealed that ceRNAs are involved in oxidation–reduction, peroxidase activity, lignin synthesis in the xylem, and flavonoid synthesis. Overall, these findings may help researchers better understand the role of non-coding RNAs in response to root-knot nematodes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08470-3.

Construction of a Quantitative Genomic Map, Identification and Expression Analysis of Candidate Genes for Agronomic and Disease-Related Traits in Brassica napus

Rapeseed is the second most important oil crop in the world. Improving seed yield and seed oil content are the two main highlights of the research. Unfortunately, rapeseed development is frequently affected by different diseases. Extensive research has been made through many years to develop elite cultivars with high oil, high yield, and/or disease resistance. Quantitative trait locus (QTL) analysis has been one of the most important strategies in the genetic deciphering of agronomic characteristics. To comprehend the distribution of these QTLs and to uncover the key regions that could simultaneously control multiple traits, 4,555 QTLs that have been identified during the last 25 years were aligned in one unique map, and a quantitative genomic map which involved 128 traits from 79 populations developed in 12 countries was constructed. The present study revealed 517 regions of overlapping QTLs which harbored 2,744 candidate genes and might affect multiple traits, simultaneously. They could be selected to customize super-rapeseed cultivars. The gene ontology and the interaction network of those candidates revealed genes that highly interacted with the other genes and might have a strong influence on them. The expression and structure of these candidate genes were compared in eight rapeseed accessions and revealed genes of similar structures which were expressed differently. The present study enriches our knowledge of rapeseed genome characteristics and diversity, and it also provided indications for rapeseed molecular breeding improvement in the future.

Genome-wide association study dissects the genetic control of plant height and branch number in response to low-phosphorus stress in Brassica napus

BACKGROUND AND AIMS

Oilseed rape (Brassica napus) is one of the most important oil crops worldwide. Phosphorus (P) deficiency severely decreases the plant height and branch number of B. napus. However, the genetic bases controlling plant height and branch number in B. napus under P deficiency remain largely unknown. This study aims to mine candidate genes for plant height and branch number by genome-wide association study (GWAS) and determine low-P-tolerance haplotypes.

METHODS

An association panel of B. napus was grown in the field with a low P supply (P, 0 kg ha-1) and a sufficient P supply (P, 40 kg ha-1) across 2 years and plant height and branch number were investigated. More than five million single-nucleotide polymorphisms (SNPs) were used to conduct GWAS of plant height and branch number at two contrasting P supplies.

KEY RESULTS

A total of 2127 SNPs were strongly associated (P < 6·25 × 10-07) with plant height and branch number at two P supplies. There was significant correlation between phenotypic variation and the number of favourable alleles of associated loci on chromosomes A10 (chrA10_821671) and C08 (chrC08_27999846), which will contribute to breeding improvement by aggregating these SNPs. BnaA10g09290D and BnaC08g26640D were identified to be associated with chrA10_821671 and chrC08_27999846, respectively. Candidate gene association analysis and haplotype analysis showed that the inbred lines carrying ATT at BnaA10g09290Hap1 and AAT at BnaC08g26640Hap1 had greater plant height than lines carrying other haplotype alleles at low P supply.

CONCLUSION

Our results demonstrate the power of GWAS in identifying genes of interest in B. napus and provided insights into the genetic basis of plant height and branch number at low P supply in B. napus. Candidate genes and favourable haplotypes may facilitate marker-based breeding efforts aimed at improving P use efficiency in B. napus.

The Bnapus50K array: a quick and versatile genotyping tool for Brassica napus genomic breeding and research

Rapeseed is a globally cultivated commercial crop, primarily grown for its oil. High-density single nucleotide polymorphism (SNP) arrays are widely used as a standard genotyping tool for rapeseed research, including for gene mapping, genome-wide association studies, germplasm resource analysis, and cluster analysis. Although considerable rapeseed genome sequencing data have been released, DNA arrays are still an attractive choice for providing additional genetic data in an era of high-throughput whole-genome sequencing. Here, we integrated re-sequencing DNA array data (32,216, 304 SNPs) from 505 inbred rapeseed lines, allowing us to develop a sensitive and efficient genotyping DNA array, Bnapus50K, with a more consistent genetic and physical distribution of probes. A total of 42,090 high-quality probes were filtered and synthesized, with an average distance between adjacent SNPs of 8 kb. To improve the practical application potential of this array in rapeseed breeding, we also added 1,618 functional probes related to important agronomic traits such as oil content, disease resistance, male sterility, and flowering time. The additional probes also included those specifically for detecting genetically modified material. These probes show a good detection efficiency and are therefore useful for gene mapping, along with crop variety improvement and identification. The novel Bnapus50K DNA array developed in this study could prove to be a quick and versatile genotyping tool for B. napus genomic breeding and research.

Integrating a genome-wide association study with transcriptomic data to predict candidate genes and favourable haplotypes influencing Brassica napus seed phytate

Phytate is the storage form of phosphorus in angiosperm seeds and plays vitally important roles during seed development. However, in crop plants phytate decreases bioavailability of seed-sourced mineral elements for humans, livestock and poultry, and contributes to phosphate-related water pollution. However, there is little knowledge about this trait in oilseed rape (Brassica napus). Here, a panel of 505 diverse B. napus accessions was screened in a genome-wide association study (GWAS) using 3.28 × 10⁶ single-nucleotide polymorphisms (SNPs). This identified 119 SNPs significantly associated with phytate concentration (PA_Conc) and phytate content (PA_Cont) and six candidate genes were identified. Of these, BnaA9.MRP5 represented the candidate gene for the significant SNP chrA09_5198034 (27 kb) for both PA_Cont and PA_Conc. Transcription of BnaA9.MRP5 in a low-phytate variety (LPA20) was significantly elevated compared with a high-phytate variety (HPA972). Association and haplotype analysis indicated that inbred lines carrying specific SNP haplotypes within BnaA9.MRP5 were associated with high- and low-phytate phenotypes. No significant differences in seed germination and seed yield were detected between low and high phytate cultivars examined. Candidate genes, favourable haplotypes and the low phytate varieties identified in this study will be useful for low-phytate breeding of B. napus.

Genome-wide association study and transcriptome analysis reveal key genes affecting root growth dynamics in rapeseed

BACKGROUND

In terms of global demand, rapeseed is the third-largest oilseed crop after soybeans and palm, which produces vegetable oil for human consumption and biofuel for industrial production. Roots are vital organs for plant to absorb water and attain mineral nutrients, thus they are of great importance to plant productivity. However, the genetic mechanisms regulating root development in rapeseed remain unclear. In the present study, seven root-related traits and shoot biomass traits in 280 Brassica napus accessions at five continuous vegetative stages were measured to establish the genetic basis of root growth in rapeseed.

RESULTS

The persistent and stage-specific genetic mechanisms were revealed by root dynamic analysis. Sixteen persistent and 32 stage-specific quantitative trait loci (QTL) clusters were identified through genome-wide association study (GWAS). Root samples with contrasting (slow and fast) growth rates throughout the investigated stages and those with obvious stage-specific changes in growth rates were subjected to transcriptome analysis. A total of 367 differentially expressed genes (DEGs) with persistent differential expressions throughout root development were identified, and these DEGs were significantly enriched in GO terms, such as energy metabolism and response to biotic or abiotic stress. Totally, 485 stage-specific DEGs with different expressions at specific stage were identified, and these DEGs were enriched in GO terms, such as nitrogen metabolism. Four candidate genes were identified as key persistent genetic factors and eight as stage-specific ones by integrating GWAS, weighted gene co-expression network analysis (WGCNA), and differential expression analysis. These candidate genes were speculated to regulate root system development, and they were less than 100 kb away from peak SNPs of QTL clusters. The homologs of three genes (BnaA03g52990D, BnaA06g37280D, and BnaA09g07580D) out of 12 candidate genes have been reported to regulate root development in previous studies.

CONCLUSIONS

Sixteen QTL clusters and four candidate genes controlling persistently root development, and 32 QTL clusters and eight candidate genes stage-specifically regulating root growth in rapeseed were detected in this study. Our results provide new insights into the temporal genetic mechanisms of root growth by identifying key candidate QTL/genes in rapeseed.

Genetic Dissection of Root Angle of Brassica napus in Response to Low Phosphorus

Plant root angle determines the vertical and horizontal distribution of roots in the soil layer, which further influences the acquisition of phosphorus (P) in topsoil. Large genetic variability for the lateral root angle (root angle) was observed in a linkage mapping population (BnaTNDH population) and an association panel of Brassica napus whether at a low P (LP) or at an optimal P (OP). At LP, the average root angle of both populations became smaller. Nine quantitative trait loci (QTLs) at LP and three QTLs at OP for the root angle and five QTLs for the relative root angle (RRA) were identified by the linkage mapping analysis in the BnaTNDH population. Genome-wide association studies (GWASs) revealed 11 single-nucleotide polymorphisms (SNPs) significantly associated with the root angle at LP (LPRA). The interval of a QTL for LPRA on A06 (qLPRA-A06c) overlapped with the confidence region of the leading SNP (Bn-A06-p14439400) significantly associated with LPRA. In addition, a QTL cluster on chromosome C01 associated with the root angle and the primary root length (PRL) in the "pouch and wick" high-throughput phenotyping (HTP) system, the root P concentration in the agar system, and the seed yield in the field was identified in the BnaTNDH population at LP. A total of 87 genes on A06 and 192 genes on C01 were identified within the confidence interval, and 14 genes related to auxin asymmetric redistribution and root developmental process were predicted to be candidate genes. The identification and functional analyses of these genes affecting LPRA are of benefit to the cultivar selection with optimal root system architecture (RSA) under P deficiency in Brassica napus.

Bioassociation of U(VI) and Eu(III) by Plant (Brassica napus) Suspension Cell Cultures-A Spectroscopic Investigation

In this study, we investigated the interaction of U(VI) and Eu(III) with Brassica napus suspension plant cells as a model system. Concentration-dependent (0-200 μM) bioassociation experiments showed that more than 75% of U(VI) and Eu(III) were immobilized by the cells. In addition to this phenomenon, time-dependent studies for 1 to 72 h of exposure showed a multistage bioassociation process for cells that were exposed to 200 μM U(VI), where, after initial immobilization of U(VI) within 1 h of exposure, it was released back into the culture medium starting within 24 h. A remobilization to this extent has not been previously observed. The MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide) assay was used to correlate the bioassociation behavior of Eu and U with the cell vitality. Speciation studies by spectroscopy and in silico methods highlighted various U and Eu species over the course of exposure. We were able to observe a new U species, which emerged simultaneously with the remobilization of U back into the solution, which we assume to be a U(VI) phosphate species. Thus, the interaction of U(VI) and Eu(III) with released plant metabolites could be concluded.

Dynamic Responses of Barley Root Succinyl-Proteome to Short-Term Phosphate Starvation and Recovery

Barley (Hordeum vulgare L.)-a major cereal crop-has low Pi demand, which is a distinct advantage for studying the tolerance mechanisms of phosphorus deficiency. We surveyed dynamic protein succinylation events in barley roots in response to and recovery from Pi starvation by firstly evaluating the impact of Pi starvation in a Pi-tolerant (GN121) and Pi-sensitive (GN42) barley genotype exposed to long-term low Pi (40 d) followed by a high-Pi recovery for 10 d. An integrated proteomics approach involving label-free, immune-affinity enrichment, and high-resolution LC-MS/MS spectrometric analysis was then used to quantify succinylome and proteome in GN121 roots under short-term Pi starvation (6, 48 h) and Pi recovery (6, 48 h). We identified 2,840 succinylation sites (Ksuc) across 884 proteins; of which, 11 representative Ksuc motifs had the preferred amino acid residue (lysine). Furthermore, there were 81 differentially abundant succinylated proteins (DFASPs) from 119 succinylated sites, 83 DFASPs from 110 succinylated sites, 93 DFASPs from 139 succinylated sites, and 91 DFASPs from 123 succinylated sites during Pi starvation for 6 and 48 h and during Pi recovery for 6 and 48 h, respectively. Pi starvation enriched ribosome pathways, glycolysis, and RNA degradation. Pi recovery enriched the TCA cycle, glycolysis, and oxidative phosphorylation. Importantly, many of the DFASPs identified during Pi starvation were significantly overexpressed during Pi recovery. These results suggest that barley roots can regulate specific Ksuc site changes in response to Pi stress as well as specific metabolic processes. Resolving the metabolic pathways of succinylated protein regulation characteristics will improve phosphate acquisition and utilization efficiency in crops.

Genome-Wide Association Analysis for Phosphorus Use Efficiency Traits in Mungbean (Vigna radiata L. Wilczek) Using Genotyping by Sequencing Approach

Mungbean (Vigna radiata L. Wilczek) is an annual grain legume crop affected by low availability of phosphorus. Phosphorus deficiency mainly affects the growth and development of plants along with changes in root morphology and increase in root-to-shoot ratio. Deciphering the genetic basis of phosphorus use efficiency (PUE) traits can benefit our understanding of mungbean tolerance to low-phosphorus condition. To address this issue, 144 diverse mungbean genotypes were evaluated for 12 PUE traits under hydroponics with optimum- and low-phosphorus levels. The broad sense heritability of traits ranged from 0.63 to 0.92 and 0.58 to 0.92 under optimum- and low-phosphorus conditions, respectively. This study, reports for the first time such a large number of genome wide Single nucleotide polymorphisms (SNPs) (76,160) in mungbean. Further, genome wide association study was conducted using 55,634 SNPs obtained by genotyping-by-sequencing method. The results indicated that total 136 SNPs shared by both GLM and MLM models were associated with tested PUE traits under different phosphorus regimes. We have identified SNPs with highest p value (-log10(p)) for some traits like, TLA and RDW with p value (-log10(p)) of more than 6.0 at LP/OP and OP condition. We have identified nine SNPs (three for TLA and six for RDW trait) which was found to be present in chromosomes 8, 4, and 7. One SNP present in Vradi07g06230 gene contains zinc finger CCCH domain. In total, 71 protein coding genes were identified, of which 13 genes were found to be putative candidate genes controlling PUE by regulating nutrient uptake and root architectural development pathways in mungbean. Moreover, we identified three potential candidate genes VRADI11G08340, VRADI01G05520, and VRADI04G10750 with missense SNPs in coding sequence region, which results in significant variation in protein structure at tertiary level. The identified SNPs and candidate genes provide the essential information for genetic studies and marker-assisted breeding program for improving low-phosphorus tolerance in mungbean.

Genome-wide identification of meiotic recombination hot spots detected by SLAF in peanut (Arachis hypogaea L.)

Recombination hot spots (RHP), caused by meiosis, are considered to play crucial roles in improvement and domestication of crop. Cultivated peanut is one of the most important rich-source of oil and protein crops. However, no direct scale of recombination events and RHP have been estimated for peanut. To examine the scale of recombination events and RHP in peanut, a RIL population with 200 lines and a natural population with 49 cultivars were evaluated. The precise integrated map comprises 4837 SLAF markers with genetic length of 2915.46 cM and density of 1.66 markers per cM in whole genome. An average of 30.0 crossover (2.06 cMMb⁻¹) events was detected per RIL plant. The crossover events (CE) showed uneven distribution among B sub-genome (2.32) and A sub-genome (1.85). There were 4.34% and 7.86% of the genome contained large numbers of CE (> 50 cMMb⁻¹) along chromosomes in F₆ and natural population, respectively. High density of CE regions called RHP, showed negative relationship to marker haplotypes conservative region but positive to heatmap of recombination. The genes located within the RHP regions by GO categories showed the responding of environmental stimuli, which suggested that recombination plays a crucial role in peanut adaptation to changing environments

Genetic variation in eggplant for Nitrogen Use Efficiency under contrasting NO3 - supply

Eggplant (Solanum melongena L.) yield is highly sensitive to N fertilization, the excessive use of which is responsible for environmental and human health damage. Lowering N input together with the selection of improved Nitrogen-Use-Efficiency (NUE) genotypes, more able to uptake, utilize, and remobilize N available in soils, can be challenging to maintain high crop yields in a sustainable agriculture. The aim of this study was to explore the natural variation among eggplant accessions from different origins, in response to Low (LN) and High (HN) Nitrate (NO₃ ^- ) supply, to identify NUE-contrasting genotypes and their NUE-related traits, in hydroponic and greenhouse pot experiments. Two eggplants, AM222 and AM22, were identified as N-use efficient and inefficient, respectively, in hydroponic, and these results were confirmed in a pot experiment, when crop yield was also evaluated. Overall, our results indicated the key role of N-utilization component (NUtE) to confer high NUE. The remobilization of N from leaves to fruits may be a strategy to enhance NUtE, suggesting glutamate synthase as a key enzyme. Further, omics technologies will be used for focusing on C-N metabolism interacting networks. The availability of RILs from two other selected NUE-contrasting genotypes will allow us to detect major genes/quantitative trait loci related to NUE.

Root diversity in sesame (Sesamum indicum L.): insights into the morphological, anatomical and gene expression profiles

Sesame harbors a large diversity in root morphological and anatomical traits and a high root biomass improves the plant aboveground biomass as well as the seed yield. Sesame provides one of the most nutritious and healthy vegetable oils, sparking an increasing demand of its seeds. However, with the low yield and productivity of sesame, there is still a huge gap between the seed demand and supply. Improving the root system has a high potential to increase crop productivity, but information on the diversity of the sesame root systems is still lacking. In this study, 40 diverse sesame varieties were grown in soil and hydroponics systems and the diversity of the root system was investigated. The results showed that sesame holds a large root morphological and anatomical diversity, which can be harnessed in breeding programmes. Based on the clustering of the genotypes in hydroponics and soil culture systems, we found that similar genotypes were commonly clustered either in the small-root or in the big-root group, indicating that the hydroponics system can be employed for a large-scale root phenotyping. Our results further revealed that the root biomass positively contributes to increased seed yield in sesame, based on multi-environmental trials. By comparing the root transcriptome of two contrasting genotypes, 2897 differentially expressed genes were detected and they were enriched in phenylpropanoid biosynthesis, starch and sucrose metabolism, stilbenoid, diarylheptanoid and gingerol biosynthesis, flavonoid biosynthesis, suggesting that these pathways are crucial for sesame root growth and development. Overall, this study sheds light on the diversity of sesame root system and offers the basis for improving root traits and increasing sesame seed yield.

Root plasticity and Pi recycling within plants contribute to low-P tolerance in Tibetan wild barley

BACKGROUND

Barley is a low phosphorus (P) demand cereal crop. Tibetan wild barley, as a progenitor of cultivated barley, has revealed outstanding ability of tolerance to low-P stress. However, the underlying mechanisms of low-P adaption and the relevant genetic controlling are still unclear.

RESULTS

We identified low-P tolerant barley lines in a doubled-haploid (DH) population derived from an elite Tibetan wild barley accession and a high-yield cultivar. The tolerant lines revealed greater root plasticity in the terms of lateral root length, compared to low-P sensitive lines, in response to low-P stress. By integrating the QTLs associated with root length and root transcriptomic profiling, candidate genes encoding isoflavone reductase, nitrate reductase, nitrate transporter and transcriptional factor MYB were identified. The differentially expressed genes (DEGs) involved the growth of lateral root, Pi transport within cells as well as from roots to shoots contributed to the differences between low-P tolerant line L138 and low-P sensitive lines L73 in their ability of P acquisition and utilization.

CONCLUSIONS

The plasticity of root system is an important trait for barley to tolerate low-P stress. The low-P tolerance in the elite DH line derived from a cross of Tibetan wild barley and cultivated barley is characterized by enhanced growth of lateral root and Pi recycling within plants under low-P stress.

Temporal genetic patterns of root growth in Brassica napus L. revealed by a low-cost, high-efficiency hydroponic system

Application of a low-cost and high-efficiency hydroponic system in a rapeseed population verified two types of genetic factors ("persistent" and "stage-specific") that control root development. The root system is a vital plant component for nutrient and water acquisition and is targeted to enhance plant productivity. Genetic dissection of the root system generally focuses on a single stage, but roots grow continuously during plant development. To reveal the temporal genetic patterns of root development, we measured nine root-related traits in a rapeseed recombinant inbred line population at six continuous stages during vegetative growth, using a modified hydroponic system with low-cost and high-efficiency features that could synchronize plant growth under controlled conditions. Phenotypic correlation and growth dynamic analysis suggested the existence of two types of genetic factors ("persistent" and "stage-specific") that control root development. Dynamic (unconditional and conditional) quantitative trait loci (QTL) mapping detected 28 stage-specific and 23 persistent QTLs related to root growth. Among them, 13 early stage-specific, 19 persistent and 8 later stage-specific QTLs were detected at 7 DAS (days after sowing), 16 DAS and 5 EL (expanding leaf stage), respectively, providing efficient and adaptable stages for QTL identification. The effective prediction of biomass accumulation using root morphological traits (up to 96.6% or 92.64% at a specific stage or the final stage, respectively) verified that root growth allocation with maximum root uptake area facilitated biomass accumulation. Furthermore, marker-assistant selection, which combined the "persistent" and "stage-specific" QTLs, proved their effectiveness for root improvement with an excellent uptake area. Our results highlight the potential of high-throughput and precise phenotyping to assess the dynamic genetics of root growth and provide new insights into ideotype root system-based biomass breeding.

Shovelomics for phenotyping root architectural traits of rapeseed/canola (Brassica napus L.) and genome-wide association mapping

Root system in plants plays an important role in mining moisture and nutrients from the soil and is positively correlated to yield in many crops including rapeseed/canola (Brassica napus L.). Substantial phenotypic diversity in root architectural traits among the B. napus growth types leads to a scope of root system improvement in breeding populations. In this study, 216 diverse genotypes were phenotyped for five different root architectural traits following shovelomics approach in the field condition during 2015 and 2016. A single nucleotide polymorphism (SNP) marker panel consisting of 30,262 SNPs was used to conduct genome-wide association study to detect marker/trait association. A total of 31 significant marker loci were identified at 0.01 percentile tail P value cutoff for different root traits. Six marker loci for soil-level taproot diameter (R₁Dia), six loci for belowground taproot diameter (R₂Dia), seven loci for number of primary root branches (PRB), eight loci for root angle, and eight loci for root score (RS) were detected in this study. Several markers associated with root diameters R₁Dia and R₂Dia were also associated with PRB and RS. Significant phenotypic correlation between these traits was observed in both environments. Therefore, taproot diameter appears to be a major determinant of the canola root system architecture and can be used as proxy for other root traits. Fifteen candidate genes related to root traits and root development were detected within 100 kbp upstream and downstream of different significant markers. The identified markers associated with different root architectural traits can be considered for marker-assisted selection for root traits in canola in future.

Genomic and Transcriptomic Analysis Identified Gene Clusters and Candidate Genes for Oil Content in Peanut (Arachis hypogaea L.)

Peanut (Arachis hypogaea), a major source of vegetable oil in many Asian countries, has become an integral part of human diet globally due to its high nutritional properties and option to consume in different forms. In order to meet the demand of vegetable oil, many peanut breeding programs of China have intensified their efforts in increasing oil content in newly bred varieties for reducing the import of edible oils in China. In this context, transcriptome sequencing data generated on 49 peanut cultivars were analyzed to identify candidate genes and develop molecular markers for seed oil content across multiple environments. Transcriptome analysis identified 5458 differentially expressed genes (DEGs) including 2243 positive DEGs and 3215 negative DEGs involved in oil synthesis process. Genome-wide association study identified 48 significant insertion/deletion (InDel) markers associated with seed oil content across five environments. A comparative genomics and transcriptomics analysis detected a total of 147 common gene clusters located in 17 chromosomes. Interestingly, an InDel cluster associated with seed oil content on A03 chromosome was detected in three different environments. Candidate genes identified on A03 form a haplotype, in which variable alleles were found to be different in oil content in an independent population. This locus is important for understanding the genetic control of peanut oil content and may be useful for marker-assisted selection in peanut breeding programs.

Genetic Dissection of Root Morphological Traits Related to Nitrogen Use Efficiency in Brassica napus L. under Two Contrasting Nitrogen Conditions

As the major determinant for nutrient uptake, root system architecture (RSA) has a massive impact on nitrogen use efficiency (NUE). However, little is known the molecular control of RSA as related to NUE in rapeseed. Here, a rapeseed recombinant inbred line population (BnaZNRIL) was used to investigate root morphology (RM, an important component for RSA) and NUE-related traits under high-nitrogen (HN) and low-nitrogen (LN) conditions by hydroponics. Data analysis suggested that RM-related traits, particularly root size had significantly phenotypic correlations with plant dry biomass and N uptake irrespective of N levels, but no or little correlation with N utilization efficiency (NUtE), providing the potential to identify QTLs with pleiotropy or specificity for RM- and NUE-related traits. A total of 129 QTLs (including 23 stable QTLs, which were repeatedly detected at least two environments or different N levels) were identified and 83 of them were integrated into 22 pleiotropic QTL clusters. Five RM-NUE, ten RM-specific and three NUE-specific QTL clusters with same directions of additive-effect implied two NUE-improving approaches (RM-based and N utilization-based directly) and provided valuable genomic regions for NUE improvement in rapeseed. Importantly, all of four major QTLs and most of stable QTLs (20 out of 23) detected here were related to RM traits under HN and/or LN levels, suggested that regulating RM to improve NUE would be more feasible than regulating N efficiency directly. These results provided the promising genomic regions for marker-assisted selection on RM-based NUE improvement in rapeseed.