Incorporating pleiotropic quantitative trait loci in dissection of complex traits: seed yield in rapeseed as an example

Key message

A comprehensive linkage atlas for seed yield in rapeseed.

Abstract

Most agronomic traits of interest for crop improvement (including seed yield) are highly complex quantitative traits controlled by numerous genetic loci, which brings challenges for comprehensively capturing associated markers/genes. We propose that multiple trait interactions underlie complex traits such as seed yield, and that considering these component traits and their interactions can dissect individual quantitative trait loci (QTL) effects more effectively and improve yield predictions. Using a segregating rapeseed (Brassica napus) population, we analyzed a large set of trait data generated in 19 independent experiments to investigate correlations between seed yield and other complex traits, and further identified QTL in this population with a SNP-based genetic bin map. A total of 1904 consensus QTL accounting for 22 traits, including 80 QTL directly affecting seed yield, were anchored to the B. napus reference sequence. Through trait association analysis and QTL meta-analysis, we identified a total of 525 indivisible QTL that either directly or indirectly contributed to seed yield, of which 295 QTL were detected across multiple environments. A majority (81.5%) of the 525 QTL were pleiotropic. By considering associations between traits, we identified 25 yield-related QTL previously ignored due to contrasting genetic effects, as well as 31 QTL with minor complementary effects. Implementation of the 525 QTL in genomic prediction models improved seed yield prediction accuracy. Dissecting the genetic and phenotypic interrelationships underlying complex quantitative traits using this method will provide valuable insights for genomics-based crop improvement.

     

Identification of candidate genes of QTLs for seed weight in <Emphasis Type="Italic">Brassica napus</Emphasis> through comparative mapping among <Emphasis Type="Italic">Arabidopsis</Emphasis> and <Emphasis Type="Italic">Brassica</Emphasis>species

Background

Map-based cloning of quantitative trait loci (QTLs) in polyploidy crop species remains a challenge due to the complexity of their genome structures. QTLs for seed weight in B. napus have been identified, but information on candidate genes for identified QTLs of this important trait is still rare.

Results

In this study, a whole genome genetic linkage map for B. napus was constructed using simple sequence repeat (SSR) markers that covered a genetic distance of 2,126.4 cM with an average distance of 5.36 cM between markers. A procedure was developed to establish colinearity of SSR loci on B. napus with its two progenitor diploid species B. rapa and B. oleracea through extensive bioinformatics analysis. With the aid of B. rapa and B. oleracea genome sequences, the 421 homologous colinear loci deduced from the SSR loci of B. napus were shown to correspond to 398 homologous loci in Arabidopsis thaliana. Through comparative mapping of Arabidopsis and the three Brassica species, 227 homologous genes for seed size/weight were mapped on the B. napus genetic map, establishing the genetic bases for the important agronomic trait in this amphidiploid species. Furthermore, 12 candidate genes underlying 8 QTLs for seed weight were identified, and a gene-specific marker for BnAP2 was developed through molecular cloning using the seed weight/size gene distribution map in B. napus.

Conclusions

Our study showed that it is feasible to identify candidate genes of QTLs using a SSR-based B. napus genetic map through comparative mapping among Arabidopsis and B. napus and its two progenitor species B. rapa and B. oleracea. Identification of candidate genes for seed weight in amphidiploid B. napus will accelerate the process of isolating the mapped QTLs for this important trait, and this approach may be useful for QTL identification of other traits of agronomic significance.

     

QTL analysis and the development of closely linked markers for days to flowering in spring oilseed rape (<Emphasis Type="Italic">Brassica napus</Emphasis> L.)

Days to flowering (DTF) is an important trait impacting cultivar performance in oilseed rape (Brassica napus L.), but the interaction of all loci controlling this trait in spring-type oilseed rape is not fully understood. We identified quantitative trait loci (QTL) for variation in DTF in a doubled haploid (DH) population from the Qinghai–Tibet Plateau that includes 217 lines derived from a cross between spring-type oilseed rape (B. napus L.) line No. 5246 and line No. 4512, the latter of which is responsive to the effective accumulated temperature (EAT). A linkage map was constructed for the DH population, using 202 SSR and 293 AFLP markers. At least 22 DTF QTL were found in multiple environments. Four major QTL were located on linkage groups A7, C2, C8 and C8. Among these QTL, cqDTFA7a and cqDTFC2a were identified in five environments and individually explained 10.4 and 23.0 % of the trait variation, respectively. cqDTFC8, a major QTL observed in spring environments, and a unique winter environment QTL, qDTFC8-3, were identified; these QTL explained 10.0 and 46.5 % of the phenotypic variation, respectively. Minor QTL (for example, cqDTFC2c) and epistatic interactions seemed evident in this population. Two closely linked SSR markers for cqDTFA7a and cqDTFC8 were developed (G1803 and S034). BnAP1, a B. napus gene with homology to Arabidopsis thaliana that was identified as a cqDTFA7a candidate gene, played a major role in this study. The allelic effects of the major and minor QTL on DTF were further validated in the DH population and in 93 breeding genotypes.     

Mapping of days to flower and seed yield in spring oilseed <Emphasis Type="Italic">Brassica napus</Emphasis> carrying genome content introgressed from <Emphasis Type="Italic">Brassica oleracea</Emphasis>

Earliness of flowering and maturity and high seed yield are important objectives of breeding spring Brassica napus canola. Previously, we have introgressed earliness of flowering from Brassica oleracea into spring B. napus canola through interspecific crossing between these two species. In this paper, we report quantitative trait locus (QTL) mapping of days to flower and seed yield by use of publicly available markers and markers designed based on flowering time genes and a doubled haploid population, derived from crossing of the spring canola parent and an early flowering line developed from a B. napus × B. oleracea cross, tested in nine field trials for over 5 years. Five genomic regions associated with days to flower were identified on C1, C2, C3, and C6 of which the single QTL of C1 was detected in all trials; in all cases, the allele introgressed from B. oleracea reduced the number of days to flower. BLASTn search in the Brassica genomes located the physical position of the QTL markers and identified putative flowering time genes in these regions. In the case of seed yield, ten QTL from eight linkage groups were detected; however, none could be consistently detected in all trials. The QTL region of C1 associated with days to flower did not show significant association with seed yield in more than 80% of the field trials; however, in a single trial, the allele introgressed from B. oleracea exerted a negative effect on seed yield. Thus, the genomic regions and molecular markers identified in this research could potentially be used in breeding for the development of early flowering B. napus canola cultivars without affecting seed yield in a majority of the environments.     

Molecular basis of the high-palmitic acid trait in sunflower seed oil

Sunflower (Helianthus annuus L.) seed oil with high palmitic acid content has enhanced thermo-oxidative stability, which makes it well suited to high-temperature uses. CAS-5 is a sunflower mutant line that accumulates over 25 % palmitic acid in its seed oil, compared to 5–8 % in conventional cultivars. The objective of this study was to investigate the molecular basis of the high-palmitic acid trait in CAS-5 through both candidate gene and QTL mapping approaches. An F₂ population derived from the cross between CAS-5 and the conventional line HA-89 was developed. A 3-ketoacyl-ACP synthase II (KASII) locus on a telomeric region of linkage group (LG) 9 of the sunflower genetic map was found to co-segregate with palmitic acid content in this population. The KASII locus explained the vast majority of the phenotypic variation (98 %) of the trait. Two minor QTL affecting palmitic acid content were also found on the lower half of LG 9 and on LG 17. Additionally, QTL associated with other major fatty acids (stearic, oleic, and linoleic acid) were identified on LG 1, 6, and 10. This result may reflect untapped genetic variation that could exist among sunflower cultivars for genes determining fatty acid composition. In addition to demonstrating the major role of a KASII locus in the accumulation of high levels of palmitic acid in CAS-5 seeds, this study stressed the importance of characterizing genes with minor effects on fatty acid profile in order to establish optimal breeding strategies for modifying fatty acid composition in sunflower seed oil.     

Genetic characterization and fine mapping for multi-inflorescence in <Emphasis Type="Italic">Brassica napus</Emphasis> L.

Key message

A major QTL for multi-inflorescence was mapped to a 27.18-kb region on A05 in Brassica napus by integrating QTL mapping, microarray analysis and whole-genome sequencing.

Abstract

Multi-inflorescence is a desirable trait for the genetic improvement of rapeseed (Brassica napus L.). However, the genetic mechanism underlying the multi-inflorescence trait is not well understood. In the present study, a doubled haploid (DH) population derived from a cross between single- and multi-inflorescence lines was investigated for the penetrance of multi-inflorescence across 3 years and genotyped with 257 simple sequence repeat and sequence-related amplified polymorphism loci. A major quantitative trait locus (QTL) for penetrance of multi-inflorescence was mapped to a 9.31-Mb region on chromosome A05, explaining 45.81% of phenotypic variance on average. Subsequently, 13 single-inflorescence and 15 multi-inflorescence DH lines were genotyped with the Brassica microarray, and the QTL interval of multi-inflorescence was narrowed to a 0.74-Mb region with 37 successive single nucleotide polymorphisms between single- and multi-inflorescence groups. A 27.18-kb QTL interval was detected by screening 420 recessive F₂ individuals with genome-specific markers. These results will be valuable for gene cloning and molecular breeding of multi-inflorescence in rapeseed.

     

Genetic mapping of clubroot resistance genes in oilseed rape

Clubroot caused by the obligate biotrophic protist Plasmodiophora brassicae is a major disease of Brassica species. Clubroot resistances introduced from B. oleracea var. ‘Böhmerwaldkohl’ and resistance from B. rapa ECD-04 were genetically mapped in oilseed rape (B. napus L.). A doubled haploid (DH) population of rape seed was developed by crossing a resistant DH-line derived from a resynthesized B. napus with the susceptible cultivar ‘Express’. The DH population was tested in the greenhouse against seven P. brassicae isolates showing low and high virulence toward B. oleracea or/and B. rapa. DH-lines with highest or lowest disease scores were used in a bulked segregant analysis (BSA), and 43 polymorphic AFLPs were identified. A genetic map of the whole genome was constructed using 338 AFLP and 156 anchored SSR markers. Nineteen QTL were detected on chromosomes N02, N03, N08, N13, N15, N16 and N19 giving resistance to seven different isolates. Race-specific effects were observed for all QTL, none of the QTL conferred resistance to all isolates. The phenotypic variance explained by the respective QTL ranged between 10.3 and 67.5%. All QTL could be assigned to both ancestral genomes of B. napus. In contrast to previous reports, a clear differentiation into major QTL from B. rapa and minor QTL from B. oleracea could not be found. Composite interval mapping confirmed the linkage relationships determined by BSA, thus demonstrating that markers for oligogenic traits can be selected by merely testing the distributional extremes of a segregating population.     

Impact of seed protein alleles from three soybean sources on seed composition and agronomic traits

Key message

Evaluation of seed protein alleles in soybean populations showed that an increase in protein concentration is generally associated with a decrease in oil concentration and yield.

Abstract

Soybean [Glycine max (L.) Merrill] meal is one of the most important plant-based protein sources in the world. Developing cultivars high in seed protein concentration and seed yield is a difficult task because the traits have an inverse relationship. Over two decades ago, a protein quantitative trait loci (QTL) was mapped on chromosome (chr) 20, and this QTL has been mapped to the same position in several studies and given the confirmed QTL designation cqSeed protein-003. In addition, the wp allele on chr 2, which confers pink flower color, has also been associated with increased protein concentration. The objective of our study was to evaluate the effect of cqSeed protein-003 and the wp locus on seed composition and agronomic traits in elite soybean backgrounds adapted to the Midwestern USA. Segregating populations of isogenic lines were developed to test the wp allele and the chr 20 high protein QTL alleles from Danbaekkong (PI619083) and Glycine soja PI468916 at cqSeed protein-003. An increase in protein concentration and decrease in yield were generally coupled with the high protein alleles at cqSeed protein-003 across populations, whereas the effects of wp on protein concentration and yield were variable. These results not only demonstrate the difficulty in developing cultivars with increased protein and yield but also provide information for breeding programs seeking to improve seed composition and agronomic traits simultaneously.

     

Genetic dissection of seed weight by QTL analysis and detection of allelic variation in Indian and east European gene pool lines of <Emphasis Type="Italic">Brassica juncea</Emphasis>

Key message

Seed weight QTL identified in different populations were synthesized into consensus QTL which were shown to harbor candidate genes by in silico mapping. Allelic variation inferred would be useful in breeding B. juncea lines with high seed weight.

Abstract

Seed weight is an important yield influencing trait in oilseed Brassicas and is a multigenic trait. Among the oilseed Brassicas, Brassica juncea harbors the maximum phenotypic variation wherein thousand seed weight varies from around 2.0 g to more than 7.0 g. In this study, we have undertaken quantitative trait locus/quantitative trait loci (QTL) analysis of seed weight in B. juncea using four bi-parental doubled-haploid populations. These four populations were derived from six lines (three Indian and three east European lines) with parental phenotypic values for thousand seed weight ranging from 2.0 to 7.6 g in different environments. Multi-environment QTL analysis of the four populations identified a total of 65 QTL ranging from 10 to 25 in each population. Meta-analysis of these component QTL of the four populations identified six ‘consensus’ QTL (C-QTL) in A3, A7, A10 and B3 by merging 33 of the 65 component Tsw QTL from different bi-parental populations. Allelic diversity analysis of these six C-QTL showed that Indian lines, Pusajaikisan and Varuna, hold the most positive allele in all the six C-QTL. In silico mapping of candidate genes with the consensus QTL localized 11 genes known to influence seed weight in Arabidopsis thaliana and also showed conserved crucifer blocks harboring seed weight QTL between the A subgenomes of B. juncea and B. rapa. These findings pave the way for a better understanding of the genetics of seed weight in the oilseed crop B. juncea and reveal the scope available for improvement of seed weight through marker-assisted breeding.

     

Integrated analysis of leaf morphological and color traits in different populations of Chinese cabbage (<Emphasis Type="Italic">Brassica rapa</Emphasis> ssp. <Emphasis Type="Italic">pekinensis</Emphasis>)

Key message

QTLs and candidate gene markers associated with leaf morphological and color traits were identified in two immortalized populations of Brassica rapa, which will provide genetic information for marker-assisted breeding.

Abstract

Brassica rapa is an important leafy vegetable consumed worldwide and morphology is a key character for its breeding. To enhance genetic control, quantitative trait loci (QTLs) for leaf color and plant architecture were identified using two immortalized populations with replications of 2 and 4 years. Overall, 158 and 80 QTLs associated with 23 and 14 traits were detected in the DH and RIL populations, respectively. Among them, 23 common robust-QTLs belonging to 12 traits were detected in common loci over the replications. Through comparative analysis, five crucifer genetic blocks corresponding to morphology trait (R, J&U, F and E) and color trait (F, E) were identified in three major linkage groups (A2, A3 and A7). These might be key conserved genomic regions involved with the respective traits. Through synteny analysis with Arabidopsis, 64 candidate genes involved in chlorophyll biosynthesis, cell proliferation and elongation were co-localized within QTL intervals. Among them, SCO3, ABI3, FLU, HCF153, HEMB1, CAB3 were mapped within QTLs for leaf color; and CYCD3;1, CYCB2;4, AN3, ULT1 and ANT were co-localized in QTL regions for leaf size. These robust QTLs and their candidate genes provide useful information for further research into leaf architecture with crop breeding.

     

QTL analysis of falling number and seed longevity in wheat (<Emphasis Type="Italic">Triticum aestivum</Emphasis> L.)

Pre-harvest sprouting (PHS) and seed longevity (SL) are complex biological processes of major importance for agricultural production. In the present study, a recombinant inbred line (RIL) population derived from a cross between the German winter wheat (Triticum aestivum L.) cultivars History and Rubens was used to identify genetic factors controlling these two physiological seed traits. A falling number (FN) test was employed to evaluate PHS, while SL was measured using a germination test (and the speed of germination) after controlled deterioration. FN of the population was assessed in four environments; SL traits were measured in one environment. Four major quantitative trait loci (QTL) for FN were detected on chromosomes 4D, 5A, 5D, and 7B, whereas for SL traits, a major QTL was found on chromosome 1A. The FN QTL on chromosome 4D that coincided with the position of the dwarfing gene Rht-D1b only had effects in environments that were free of PHS. The remaining three QTL for FN were mostly pronounced under conditions conducive to PHS. The QTL on the long arm of chromosome 7B corresponded to the major gene locus controlling late maturity α-amylase (LMA) in wheat. The severity of the LMA phenotype became truly apparent under sprouting conditions. The position on the long arm of chromosome 1A of the QTL for SL points to a new QTL for this important regenerative seed trait.     

Mapping of QTL for the seed storage proteins cruciferin and napin in a winter oilseed rape doubled haploid population and their inheritance in relation to other seed traits

Key message

Cruciferin (cru) and napin (nap) were negatively correlated and the cru/nap ratio was closely negative correlated with glucosinolate content indicating a link between the two biosynthetic pathways.

Abstract

Canola-type oilseed rape (Brassica napus L.) is an economically important oilseed crop in temperate zones. Apart from the oil, the canola protein shows potential as a value-added food and nutraceutical ingredient. The two major storage protein groups occurring in oilseed rape are the 2 S napins and 12 S cruciferins. The aim of the present study was to analyse the genetic variation and the inheritance of napin and cruciferin content of the seed protein in the winter oilseed rape doubled haploid population Express 617 × R53 and to determine correlations to other seed traits. Seed samples were obtained from field experiments performed in 2 years at two locations with two replicates in Germany. A previously developed molecular marker map of the DH population was used to map quantitative trait loci (QTL) of the relevant traits. The results indicated highly significant effects of the year and the genotype on napin and cruciferin content as well as on the ratio of cruciferin to napin. Heritabilities were comparatively high with 0.79 for napin and 0.77 for cruciferin. Napin and cruciferin showed a significant negative correlation (?0.36**) and a close negative correlation of the cru/nap ratio to glucosinolate content was observed (?0.81**). Three QTL for napin and two QTL for cruciferin were detected, together explaining 47 and 35 % of the phenotypic variance. A major QTL for glucosinolate content was detected on linkage group N19 whose confidence interval overlapped with QTL for napin and cruciferin content. Results indicate a relationship between seed protein composition and glucosinolate content.     

Novel and major QTL for branch angle detected by using DH population from an exotic introgression in rapeseed (<Emphasis Type="Italic">Brassica napus</Emphasis> L.)

Key message

A high-density SNP map was constructed and several novel QTL for branch angle across six environments in Brassica napus were identified.

Abstract

Branch angle is a major determinant for the ideotype of a plant, while the mechanisms underlying this trait in Brassica napus remain elusive. Herein, we developed one doubled haploid population from a cross involving one Capsella bursa-pastoris derived B. napus intertribal introgression line with the compressed branches and wooden stems, and constructed a high-density SNP map covering the genetic distance of 2242.14 cM, with an average marker interval of 0.73 cM. After phenotypic measurements across six environments, the inclusive composite interval mapping algorithm was conducted to analyze the QTL associated with branch angle. In single-environment analysis, a total of 17 QTL were detected and mainly distributed on chromosomes A01, A03, A09 and C03. Of these, three major QTL, qBA.A03-2, qBA.C03-3 and qBA.C03-4 were steadily expressed, each explaining more than 10% of the phenotypic variation in at least two environments. Compared with other results on rapeseed branch angle, these major QTL were newly detected. In QTL by environment interactions (QEI) mapping, 10 QTL were identified, and the QTL average effect and QEI effect were estimated. Of these, 7 QTL were detected in both single-environment analysis and QEI mapping. Based on the physical positions of SNPs and the functional annotation of the Arabidopsis thaliana genome, 27 genes within the QTL regions were selected as candidate genes, including early auxin-responsive genes, small auxin-up RNA, auxin/indoleacetic acid and gretchenhagen-3. These results may pave the way for deciphering the genetic control of branch angle in B. napus.

     

Inheritance pattern and genetic correlations among growth and wood quality traits in Para rubber tree (<Emphasis Type="Italic">Hevea brasiliensis</Emphasis>) and implications for breeding

Inheritance pattern of wood traits viz. specific gravity, fibre dimensions and fibre-derived biometrical indices and their interactions among themselves and with that of growth are reported in Hevea brasiliensis. Girth (h² =???0.02?±?0.44 to h² =?0.35?±?0.24) showed moderate genetic control. Among wood traits, specific gravity (h²?=?0.15?±?0.31 to h² =?0.33?±?0.28) was found to be under moderate genetic control. Fibre traits viz., fibre length (h² =???0.26?±?0.30 to h² =?0.50?±?0.34), fibre diameter (h² =?0.19?±?0.49 to h² =?0.70?±?0.11), fibre lumen diameter (h² =???0.18?±?0.35 to h² =?0.56?±?0.47) and fibre wall thickness (h² =???5.17?±?5.26 to h² =?0.50?±?0.50) were under moderate to strong genetic control. Among fibre-derived indices, flexibility coefficient (h² =?0.48?±?0.21 to h² =?0.89?±?0.29) showed moderate to very strong genetic control. The Runkel ratio (h² =???0.40?±?0.27 to h² =?0.42?±?0.29) and slenderness ratio (h² =???0.36?±?0.29 to h² =?0.43?±?0.28) showed moderate genetic control. Girth showed very strong positive genetic correlation with fibre wall thickness and strong positive correlation with fibre width indicating scope of indirect selection potential for these traits. Wood specific gravity was not correlated with either girth or fibre traits. Hence, it would be possible to concomitantly improve growth and fibre traits without adversely affecting wood specific gravity. Moderate to very high estimates of heritability for fibre traits, girth and specific gravity indicated that considerable genetic gain can be realised for these traits. Implications of the above findings in genetic improvement of wood in Hevea are discussed.     

Development of an integrated linkage map of einkorn wheat and its application for QTL mapping and genome sequence anchoring

Key message

An integrated genetic map was constructed for einkorn wheat A genome and provided valuable information for QTL mapping and genome sequence anchoring.

Abstract

Wheat is one of the most widely grown food grain crops in the world. The construction of a genetic map is a key step to organize biologically or agronomically important traits along the chromosomes. In the present study, an integrated linkage map of einkorn wheat was developed using 109 recombinant inbred lines (RILs) derived from an inter sub-specific cross, KT1-1 (T. monococcum ssp. boeoticum) × KT3-5 (T. monococcum ssp. monococcum). The map contains 926 molecular markers assigned to seven linkage groups, and covers 1,377 cM with an average marker interval of 1.5 cM. A quantitative trait locus (QTL) analysis of five agronomic traits identified 16 stable QTL on all seven chromosomes, except 6A. The total phenotypic variance explained by these stable QTL using multiple regressions varied across environments from 8.8 to 87.1 % for days to heading, 24.4–63.0 % for spike length, 48.2–79.6 % for spikelet number per spike, 13.1–48.1 % for plant architecture, and 12.2–26.5 % for plant height, revealing that much of the RIL phenotypic variation had been genetically dissected. Co-localizations of closely linked QTL for different traits were frequently observed, especially on 3A and 7A. The QTL on 3A, 5A and 7A were closely associated with Eps-A ^m 3, Vrn1 and Vrn3 loci, respectively. Furthermore, this genetic map facilitated the anchoring of 237 T. urartu scaffolds onto seven chromosomes with a physical length of 26.15 Mb. This map and the QTL data provide valuable genetic information to dissect important agronomic and developmental traits in diploid wheat and contribute to the genetic ordering of the genome assembly.

     

Evaluation of QTL alleles from exotic sources for hybrid seed yield in the original and different genetic backgrounds of spring-type <Emphasis Type="Italic">Brassica napus</Emphasis> L.

Previously identified alleles at quantitative trait loci (QTL) for hybrid seed yield were re-evaluated in the same genetic background (in hybrid combination with the same tester) as the original QTL mapping study and also evaluated in a different genetic background (in hybrid combination with two different testers). The QTL were identified from wide crosses of exotic germplasm sources with spring-type Brassica napus L., in which alleles from the exotic germplasm sources increased hybrid seed yield. Results from the re-evaluation of six QTL, in the same genetic background and hybrid combination, indicate that several of the exotic donor QTL alleles did increase hybrid seed yield and could be successfully used for improving the original single-cross hybrid. However, results from the evaluation of seven QTL (including the same six previous QTL) in a new genetic background, in combination with two new testers, indicate that the exotic QTL alleles were often no different or produced significantly lower hybrid seed yield than the spring QTL alleles. In all studies, the QTL were also very sensitive to environmental interactions. Thus, our results indicate that although these exotic sources contain favorable QTL alleles when introgressed into one spring hybrid background, the effects are not predictive of other genetic backgrounds or hybrid combinations. Although QTL affecting hybrid seed yield have been identified, comparisons of multiple QTL alleles are needed to determine the most favorable allele at each locus. Characterization of QTL complementation across testers will be required to predict their effects in multiple hybrid combinations.