A systems genomics and genetics approach to identify the genetic regulatory network for lignin content in Brassica napus seeds

Seed quality traits of oilseed rape, Brassica napus (B. napus), exhibit quantitative inheritance determined by its genetic makeup and the environment via the mediation of a complex genetic architecture of hundreds to thousands of genes. Thus, instead of single gene analysis, network-based systems genomics and genetics approaches that combine genotype, phenotype, and molecular phenotypes offer a promising alternative to uncover this complex genetic architecture. In the current study, systems genetics approaches were used to explore the genetic regulation of lignin traits in B. napus seeds. Four QTL (qLignin_A09_1, qLignin_A09_2, qLignin_A09_3, and qLignin_C08) distributed on two chromosomes were identified for lignin content. The qLignin_A09_2 and qLignin_C08 loci were homologous QTL from the A and C subgenomes, respectively. Genome-wide gene regulatory network analysis identified eighty-three subnetworks (or modules); and three modules with 910 genes in total, were associated with lignin content, which was confirmed by network QTL analysis. eQTL (expression quantitative trait loci) analysis revealed four cis-eQTL genes including lignin and flavonoid pathway genes, cinnamoyl-CoA-reductase (CCR1), and TRANSPARENT TESTA genes TT4, TT6, TT8, as causal genes. The findings validated the power of systems genetics to identify causal regulatory networks and genes underlying complex traits. Moreover, this information may enable the research community to explore new breeding strategies, such as network selection or gene engineering, to rewire networks to develop climate resilience crops with better seed quality.

A major yellow-seed QTL on chromosome A09 significantly increases the oil content and reduces the fiber content of seed in Brassica napus

A major yellow-seed QTL on chromosome A09 significantly increases the oil content and reduces the fiber content of seed in Brassica napus. The yellow-seed trait (YST) has always been a main breeding objective for rapeseed because yellow-seeded B. napus generally contains higher oil contents, fewer pigments and polyphenols and lower fiber content than black-seeded B. napus, although the mechanism controlling this correlation remains unclear. In this study, QTL mapping was implemented for YST based on a KN double haploid population derived from the hybridization of yellow-seeded B. napus N53-2 with a high oil content and black-seeded Ken-C8 with a relatively low oil content. Ten QTLs were identified, including four stable QTLs that could be detected in multiple environments. A major QTL, cqSC-A09, on chromosome A09 was identified by both QTL mapping and BSR-Seq technology, and explained more than 41% of the phenotypic variance. The major QTL cqSC-A09 for YST not only controls the seed color but also affects the oil and fiber contents in seeds. More importantly, the advantageous allele could increase the oil content and reduce the pigment and fiber content at the same time. This is the first QTL reported to control seed color, oil content and fiber content simultaneously with a large effect and has great application value for breeding high oil varieties with high seed quality. Important candidate genes, including BnaA09. JAZ1, BnaA09. GH3.3 and BnaA09. LOX3, were identified for cqSC-A09 by combining sequence variation annotation, expression differences and an interaction network, which lays a foundation for further cloning and breeding applications in the future.

QTL Genetic Mapping Study for Traits Affecting Meal Quality in Winter Oilseed Rape (Brassica Napus L.)

Rapeseed (Brassica napus L.) meal is an important source of protein, but the presence of anti-nutritional compounds, such as fibre and glucosinolates, still limits its use as a livestock feed. Understanding the genetic basis of seed fibre biosynthesis would help to manipulate its content in seeds of oilseed rape. Here, we applied high-resolution skim genotyping by sequencing (SkimGBS) and characterised 187,835 single-nucleotide polymorphism (SNP) markers across a mapping population subsequently used for a genetic mapping study (R/qtl). This approach allowed the identification of 11 stable QTL related to seed quality traits and led to the identification of potential functional genes underlying these traits. Among these, key genes with a known role in carbohydrate metabolic process, cell wall, lignin, and flavonoid biosynthesis, including cellulase GH5, TT10/LAC15, TT4, and SUC2, were found. This study furthers the understanding of the molecular mechanisms underlying seed fibre content and provides new markers for molecular breeding in B. napus.

Stable and novel QTL identification and new insights into the genetic networks affecting seed fiber traits in Brassica napus

QTL mapping for fiber-related traits and elucidation of a stable and novel QTL affecting seed lignin content, cellulose content and seed oil content. Dissection of the genetic networks for fiber biosynthesis is important for improving the seed oil content and meal value of Brassica napus. In this study, the genetic basis of seed fiber biosynthesis in B. napus was investigated via quantitative trait locus (QTL) analysis using a doubled haploid population derived from 'KenC-8' crossed with 'N53-2.' Seed lignin content (LC), cellulose content (CC) and hemicellulose content (HC) were significantly negatively correlated with seed oil content (OC). Co-localization QTLs among LC, CC, HC and OC on A09 were found with contributions ranging from 9.87 to 48.50%. Seven co-localization QTLs involved in the fiber component and OC were further verified by bulked segregant analysis (BSA). The unique QTL uqA9-12 might be a real and new QTL that was commonly identified by QTL mapping and BSA and simultaneously affected LC, CC and OC with opposite additive effects. A potential regulatory network controlling seed fiber biosynthesis was constructed to dissect the complex mechanism of seed fiber and oil accumulation, and numerous candidate genes were identified in the fiber-related QTL regions. These results provided an enrichment of QTLs and potential candidates for fiber biosynthesis, as well as useful new information for understanding the complex genetic mechanism underlying rapeseed seed fiber accumulation.

Identification of stable QTLs for seed oil content by combined linkage and association mapping in Brassica napus

Seed oil content is an important agricultural trait in rapeseed breeding. Although numerous quantitative trait locus (QTL) have been identified, most of them cannot be applied in practical breeding mainly due to environmental instability or large confidence intervals. The purpose of this study was to identify and validate high quality and more stable QTLs by combining linkage mapping and genome-wide association study (GWAS). For linkage mapping, we constructed two F₂ populations from crosses of high-oil content (∼50%) lines 6F313 and 61616 with a low-oil content (∼40%) line 51070. Two high density linkage maps spanned 1987cM (1659 bins) and 1856cM (1746 bins), respectively. For GWAS, we developed more than 34,000 high-quality SNP markers based on 227 accessions. Finally, 40 QTLs and 29 associations were established by linkage and association mapping in different environments. After merging the results, 32 consensus QTLs were obtained and 7 of them were identified by both mapping methods. Seven overlapping QTLs covered an average confidence interval of 183kb and explained the phenotypic variation of 10.23 to 24.45%. We further developed allele-specific PCR primers to identify each of the seven QTLs. These stable QTLs should be useful in gene cloning and practical breeding application.

Disruption of Germination and Seedling Development in Brassica napus by Mutations Causing Severe Seed Hormonal Imbalance

The Brassica napus (oilseed rape) accession 1012-98 shows a disturbed germination phenotype that was thought to be associated with its lack of testa pigmentation and thin seed coat. Here, we demonstrate that the disturbed germination and seedling development are actually due to independent mutations that disrupt the balance of hormone metabolites and their regulators in the seeds. High-throughput UPLC-MS/MS hormone profiling of seeds and seedlings before and after germination revealed that 1012-98 has a severely disturbed hormone balance with extremely atypical, excessive quantities of auxin and ABA metabolites. The resulting hypersensitivity to abscisic acid (ABA) and a corresponding increase in dormancy often results in death of the embryo after imbibition or high frequencies of disturbed, often lethal developmental phenotypes, resembling Arabidopsis mutants for the auxin regulatory factor gene ARF10 or the auxin-overproducing transgenic line iaaM-OX. Molecular cloning of Brassica ARF10 orthologs revealed four loci in normal B. napus, two derived from the Brassica A genome and two from the C genome. On the other hand, the phenotypic mutant 1012-98 exhibited amplification of C-genome BnaC.ARF10 copy number along with a chimeric allele originating from recombination between homeologous A and C genome loci which lead to minor increase of Bna.ARF10 transcription on the critical timepoint for seed germination, the indirect regulator of ABI3, the germinative inhibitor. Bna.GH3.5 expression was upregulated to conjugate free auxin to IAA-asp between 2 and 6 DAS. Functional amino acid changes were also found in important DNA binding domains of one BnaC.ARF10 locus, suggesting that regulatory changes in Bna.ARF10 are collectively responsible for the observed phenotpyes in 1012-98. To our knowledge, this study is the first to report disruption of germination and seedling development in Brassica napus caused by the crosstalk of auxin-ABA and the corresponding regulators Bna.ARF10 and Bna.GH3.5.

Genome-Wide Analysis of Seed Acid Detergent Lignin (ADL) and Hull Content in Rapeseed (Brassica napus L.)

A stable yellow-seeded variety is the breeding goal for obtaining the ideal rapeseed (Brassica napus L.) plant, and the amount of acid detergent lignin (ADL) in the seeds and the hull content (HC) are often used as yellow-seeded rapeseed screening indices. In this study, a genome-wide association analysis of 520 accessions was performed using the Q + K model with a total of 31,839 single-nucleotide polymorphism (SNP) sites. As a result, three significant associations on the B. napus chromosomes A05, A09, and C05 were detected for seed ADL content. The peak SNPs were within 9.27, 14.22, and 20.86 kb of the key genes BnaA.PAL4, BnaA.CAD2/BnaA.CAD3, and BnaC.CCR1, respectively. Further analyses were performed on the major locus of A05, which was also detected in the seed HC examination. A comparison of our genome-wide association study (GWAS) results and previous linkage mappings revealed a common chromosomal region on A09, which indicates that GWAS can be used as a powerful complementary strategy for dissecting complex traits in B. napus. Genomic selection (GS) utilizing the significant SNP markers based on the GWAS results exhibited increased predictive ability, indicating that the predictive ability of a given model can be substantially improved by using GWAS and GS.

Molecular breeding in Brassica for salt tolerance: importance of microsatellite (SSR) markers for molecular breeding in Brassica

Salinity is one of the important abiotic factors for any crop management in irrigated as well as rainfed areas, which leads to poor harvests. This yield reduction in salt affected soils can be overcome by improving salt tolerance in crops or by soil reclamation. Salty soils can be reclaimed by leaching the salt or by cultivation of salt tolerance crops. Salt tolerance is a quantitative trait controlled by several genes. Poor knowledge about mechanism of its inheritance makes slow progress in its introgression into target crops. Brassica is known to be a good reclamation crop. Inter and intra specific variation within Brassica species shows potential of molecular breeding to raise salinity tolerant genotypes. Among the various molecular markers, SSR markers are getting high attention, since they are randomly sparsed, highly variable and show co-dominant inheritance. Furthermore, as sequencing techniques are improving and softwares to find SSR markers are being developed, SSR markers technology is also evolving rapidly. Comparative SSR marker studies targeting Arabidopsis thaliana and Brassica species which lie in the same family will further aid in studying the salt tolerance related QTLs and subsequent identification of the "candidate genes" and finding out the origin of important QTLs. Although, there are a few reports on molecular breeding for improving salt tolerance using molecular markers in Brassica species, usage of SSR markers has a big potential to improve salt tolerance in Brassica crops. In order to obtain best harvests, role of SSR marker driven breeding approaches play important role and it has been discussed in this review especially for the introgression of salt tolerance traits in crops.

Adult plant development in triticale (× triticosecale wittmack) is controlled by dynamic genetic patterns of regulation

Many biologically and agronomically important traits are dynamic and show temporal variation. In this study, we used triticale (× Triticosecale Wittmack) as a model crop to assess the genetic dynamics underlying phenotypic plasticity of adult plant development. To this end, a large mapping population with 647 doubled haploid lines derived from four partially connected families from crosses among six parents was scored for developmental stage at three different time points. Using genome-wide association mapping, we identified main effect and epistatic quantitative trait loci (QTL) at all three time points. Interestingly, some of these QTL were identified at all time points, whereas others appear to only contribute to the genetic architecture at certain developmental stages. Our results illustrate the temporal contribution of QTL to the genetic control of adult plant development and more generally, the temporal genetic patterns of regulation that underlie dynamic traits.

Genome-wide association study dissects the genetic architecture of seed weight and seed quality in rapeseed (Brassica napus L.)

Association mapping can quickly and efficiently dissect complex agronomic traits. Rapeseed is one of the most economically important polyploid oil crops, although its genome sequence is not yet published. In this study, a recently developed 60K Brassica Infinium(®) SNP array was used to analyse an association panel with 472 accessions. The single-nucleotide polymorphisms (SNPs) of the array were in silico mapped using 'pseudomolecules' representative of the genome of rapeseed to establish their hypothetical order and to perform association mapping of seed weight and seed quality. As a result, two significant associations on A8 and C3 of Brassica napus were detected for erucic acid content, and the peak SNPs were found to be only 233 and 128 kb away from the key genes BnaA.FAE1 and BnaC.FAE1. BnaA.FAE1 was also identified to be significantly associated with the oil content. Orthologues of Arabidopsis thaliana HAG1 were identified close to four clusters of SNPs associated with glucosinolate content on A9, C2, C7 and C9. For seed weight, we detected two association signals on A7 and A9, which were consistent with previous studies of quantitative trait loci mapping. The results indicate that our association mapping approach is suitable for fine mapping of the complex traits in rapeseed.

A high-density SNP map for accurate mapping of seed fibre QTL in Brassica napus L

A high density genetic linkage map for the complex allotetraploid crop species Brassica napus (oilseed rape) was constructed in a late-generation recombinant inbred line (RIL) population, using genome-wide single nucleotide polymorphism (SNP) markers assayed by the Brassica 60 K Infinium BeadChip Array. The linkage map contains 9164 SNP markers covering 1832.9 cM. 1232 bins account for 7648 of the markers. A subset of 2795 SNP markers, with an average distance of 0.66 cM between adjacent markers, was applied for QTL mapping of seed colour and the cell wall fiber components acid detergent lignin (ADL), cellulose and hemicellulose. After phenotypic analyses across four different environments a total of 11 QTL were detected for seed colour and fiber traits. The high-density map considerably improved QTL resolution compared to the previous low-density maps. A previously identified major QTL with very high effects on seed colour and ADL was pinpointed to a narrow genome interval on chromosome A09, while a minor QTL explaining 8.1% to 14.1% of variation for ADL was detected on chromosome C05. Five and three QTL accounting for 4.7% to 21.9% and 7.3% to 16.9% of the phenotypic variation for cellulose and hemicellulose, respectively, were also detected. To our knowledge this is the first description of QTL for seed cellulose and hemicellulose in B. napus, representing interesting new targets for improving oil content. The high density SNP genetic map enables navigation from interesting B. napus QTL to Brassica genome sequences, giving useful new information for understanding the genetics of key seed quality traits in rapeseed.

Cross-validation in association mapping and its relevance for the estimation of QTL parameters of complex traits

Association mapping has become a widely applied genomic approach to identify quantitative trait loci (QTL) and dissect the genetic architecture of complex traits. However, approaches to assess the quality of the obtained QTL results are lacking. We therefore evaluated the potential of cross-validation in association mapping based on a large sugar beet data set. Our results show that the proportion of the population that should be used as estimation and validation sets, respectively, depends on the size of the mapping population. Generally, a fivefold cross-validation, that is, 20% of the lines as independent validation set, appears appropriate for commonly used population sizes. The predictive power for the proportion of genotypic variance explained by QTL was overestimated by on average 38% indicating a strong bias in the estimated QTL effects. The cross-validated predictive power ranged between 4 and 50%, which are more realistic estimates of this parameter for complex traits. In addition, QTL frequency distributions can be used to assess the precision of QTL position estimates and the robustness of the detected QTL. In summary, cross-validation can be a valuable tool to assess the quality of QTL parameters in association mapping.

Construction of genetic linkage map and mapping of QTL for seed color in Brassica rapa

A genetic linkage map of Brassica rapa L. was constructed using recombinant inbred lines (RILs) derived from a cross between yellow-seeded cultivar Sampad and a yellowish brown seeded inbred line 3-0026.027. The RILs were evaluated for seed color under three conditions: field plot, greenhouse, and controlled growth chambers. Variation for seed color in the RILs ranged from yellow, like yellow sarson, to dark brown/black even though neither parent had shown brown/black colored seeds. One major QTL (SCA9-2) and one minor QTL (SCA9-1) on linkage group (LG) A9 and two minor QTL (SCA3-1, SCA5-1) on LG A3 and LG A5, respectively, were detected. These collectively explained about 67% of the total phenotypic variance. SCA9-2 mapped in the middle of LG A9, explained about 55% phenotypic variance, and consistently expressed in all environments. The second QTL on LG A9 was ~70 cM away from SCA9-2, suggesting that independent assortment of these QTLs is possible. A digenic epistatic interaction was found between the two main effect QTL on LG A9; and the epistasis × environment interaction was nonsignificant, suggesting stability of the interaction across the environments. The QTL effect on LG A9 was validated using simple sequence repeat (SSR) markers from the two QTL regions of this LG on a B(1)S(1) population (F(1) backcrossed to Sampad followed by self-pollination) segregating for brown and yellow seed color, and on their self-pollinated progenies (B(1)S(2)). The SSR markers from the QTL region SCA9-2 showed a stronger linkage association with seed color as compared with the marker from SCA9-1. This suggests that the QTL SCA9-2 is the major determinant of seed color in the A genome of B. rapa.

A knockout mutation in the lignin biosynthesis gene CCR1 explains a major QTL for acid detergent lignin content in Brassica napus seeds

Seed coat phenolic compounds represent important antinutritive fibre components that cause a considerable reduction in value of seed meals from oilseed rape (Brassica napus). The nutritionally most important fibre compound is acid detergent lignin (ADL), to which a significant contribution is made by phenylpropanoid-derived lignin precursors. In this study, we used bulked-segregant analysis in a population of recombinant inbred lines (RILs) from a cross of the Chinese oilseed rape lines GH06 (yellow seed, low ADL) and P174 (black seed, high ADL) to identify markers with tight linkage to a major quantitative trait locus (QTL) for seed ADL content. Fine mapping of the QTL was performed in a backcross population comprising 872 BC(1)F(2) plants from a cross of an F(7) RIL from the above-mentioned population, which was heterozygous for this major QTL and P174. A 3:1 phenotypic segregation for seed ADL content indicated that a single, dominant, major locus causes a substantial reduction in ADL. This locus was successively narrowed to 0.75 cM using in silico markers derived from a homologous Brassica rapa sequence contig spanning the QTL. Subsequently, we located a B. rapa orthologue of the key lignin biosynthesis gene CINNAMOYL CO-A REDUCTASE 1 (CCR1) only 600 kbp (0.75 cM) upstream of the nearest linked marker. Sequencing of PCR amplicons, covering the full-length coding sequences of Bna.CCR1 homologues, revealed a locus in P174 whose sequence corresponds to the Brassica oleracea wild-type allele from chromosome C8. In GH06, however, this allele is replaced by a homologue derived from chromosome A9 that contains a loss-of-function frameshift mutation in exon 1. Genetic and physical map data infer that this loss-of-function allele has replaced a functional Bna.CCR1 locus on chromosome C8 in GH06 by homoeologous non-reciprocal translocation.