Breeding signature of combining ability improvement revealed by a genomic variation map from recurrent selection population in Brassica napus

Fine mapping of BnDM1-the gene regulating indeterminate inflorescence in Brassica napus

KEY MESSAGE

A candidate gene Bndm1 related to determinate inflorescence was mapped to a 128-kb interval on C02 in Brassica napus. Brassica napus plants with determinate inflorescence exhibit improved traits in field production, such as lower plant height, improved lodging resistance, and consistent maturity. Compared to plants with indeterminate inflorescence, such features are favorable for mechanized harvesting techniques. Here, using a natural mutant 6138 with determinate inflorescence, it is demonstrated that determinate inflorescence reduces plant height significantly without affecting thousand-grain weight and yield per plant. Determinacy was regulated by a single recessive gene, Bndm1. Using a combination of SNP arrays and map-based cloning, we mapped the locus of determinacy to a 128-kb region on C02. Based on sequence comparisons and the reported functions of candidate genes in this region, we predicted BnaC02.knu (a homolog of KNU in Arabidopsis) as a possible candidate gene of Bndm1 for controlling determinate inflorescence. We found a 623-bp deletion in a region upstream of the KNU promoter in the mutant. This deletion led to the significant overexpression of BnaC02.knu in the mutant compared to that in the ZS11 line. The correlation between this deletion and determinate inflorescence was examined in natural populations. The results indicated that the deletion affected the normal transcription of BnaC02.knu in the plants with determinate inflorescence and played an important role in maintaining flower development. This study presents as a new material for optimizing plant architecture and breeding novel canola varieties suitable for mechanized production. Moreover, our findings provide a theoretical basis for analyzing the molecular mechanisms underlying the formation of determinate inflorescence in B. napus.

Improvement of Resistance to Clubroot Disease in the Ogura CMS Restorer Line R2163 of Brassica napus

Oilseed rape (Brassica napus) has significant heterosis and Ogura CMS is a major way to use it. Ogura CMS has the advantages of complete and stable male sterility and easy-to-breed maintainers. Therefore, to breed better restorers has become an important goal for this system. Incidentally, clubroot is a soil-borne disease that is difficult to control by fungicidal chemicals, and it has been the main disease of oilseed rape in recent years in China, severely restricting the development of the oilseed rape industry. At present, the most effective method for controlling clubroot disease is to cultivate resistant varieties. One Ogura CMS restorer line (R2163) has shown much better combining ability, but lacks the clubroot disease resistance. This study was carried out to improve R2163 through marker-assisted backcross breeding (MABB). The resistant locus PbBa8.1 was introduced into the restorer R2163, and we then selected R2163R with clubroot disease resistance. Using the new restorer R2163R as the male parent and the sterile lines 116A and Z11A as the female parent, the improved, new resistant hybrids Kenyouza 741R and Huayouza 706R performed well, providing strong resistance and good agronomic traits. This work advances the utilization of heterosis and breeding for clubroot disease resistance in B. napus.

BnaCRCs with domestication preference positively correlate with the seed-setting rate of canola

Canola (Brassica napus) is an important oil crop worldwide. The seed-setting rate (SS) is a critical factor in determining its yield, and the development of pistils affects pollination and seed sets. However, research on seed-setting defects has been limited owing to difficulties in the identification of phenotypes, mutations, and complex genetic mechanisms. In this study, we found a stigma defect (sd) mutant in B. napus, which had no nectary. The SS of sd mutants in the field was approximately 93.4% lower than that of the wild type. Scanning and transmission electron microscopy imaging of sd mutants showed a low density of stigma papillary cells and stigma papillary cell vacuoles that disappeared 16 h after flowering. Genetic analysis of segregated populations showed that two recessive nuclear genes are responsible for the mutant phenotype of sd. Based on re-sequencing and map-based cloning, we reduced the candidate sites on ChrA07 (BnaSSA07) and ChrC06 (BnaSSC06) to 30 and 67 kb, including six and eight predicted genes, respectively. Gene analyses showed that a pair of CRABS CLAW (CRC) homeologous genes at BnaSSA07 and BnaSSC06 were associated with the development of carpel and nectary. BnaSSA07.CRC and BnaSSC06.CRC candidate genes were found to be expressed in flower organs only, with significant differences in their expression in the pistils of the near-isogenic lines. DNA sequencing showed transposon insertions in the upstream region and intron of the candidate gene BnaSSA07.crc. We also found that BnaSSC06.crc exists widely in the natural population and we give possible reasons for its widespread existence.

Comparative transcriptomic analysis reveals the molecular mechanism underlying seedling biomass heterosis in Brassica napus

BACKGROUND

Heterosis is an important biological phenomenon in which the hybrids exceed the parents in many traits. However, the molecular mechanism underlying seedling heterosis remains unclear.

RESULTS

In the present study, we analyzed the leaf transcriptomes of strong hybrids (AM, HM) and weak hybrids (CM, HW) and their parents (A, C, H, M, and W) at two periods. Phenotypically, hybrids had obvious biomass heterosis at the seedling stage, with statistically significant differences between the strong and weak hybrids. The transcriptomic analysis demonstrated that the number of differentially expressed genes (DEGs) between parents was the highest. Further analysis showed that most DEGs were biased toward parental expression. The biological processes of the two periods were significantly enriched in the plant hormone signal transduction and photosynthetic pathways. In the plant hormone signaling pathway, DEG expression was high in hybrids, with expression differences between strong and weak hybrids. In addition, DEGs related to cell size were identified. Similar changes were observed during photosynthesis. The enhanced leaf area of hybrids generated an increase in photosynthetic products, which was consistent with the phenotype of the biomass. Weighted gene co-expression network analysis of different hybrids and parents revealed that hub genes in vigorous hybrid were mainly enriched in the plant hormone signal transduction and regulation of plant hormones.

CONCLUSION

Plant hormone signaling and photosynthesis pathways, as well as differential expression of plant cell size-related genes, jointly regulate the dynamic changes between strong and weak hybrids and the generation of seedling-stage heterosis. This study may elucidate the molecular mechanism underlying early biomass heterosis and help enhance canola yield.

Genomic Selection in Sugarcane: Current Status and Future Prospects

Sugarcane is a C4 and agro-industry-based crop with a high potential for biomass production. It serves as raw material for the production of sugar, ethanol, and electricity. Modern sugarcane varieties are derived from the interspecific and intergeneric hybridization between Saccharum officinarum, Saccharum spontaneum, and other wild relatives. Sugarcane breeding programmes are broadly categorized into germplasm collection and characterization, pre-breeding and genetic base-broadening, and varietal development programmes. The varietal identification through the classic breeding programme requires a minimum of 12-14 years. The precise phenotyping in sugarcane is extremely tedious due to the high propensity of lodging and suckering owing to the influence of environmental factors and crop management practices. This kind of phenotyping requires data from both plant crop and ratoon experiments conducted over locations and seasons. In this review, we explored the feasibility of genomic selection schemes for various breeding programmes in sugarcane. The genetic diversity analysis using genome-wide markers helps in the formation of core set germplasm representing the total genomic diversity present in the Saccharum gene bank. The genome-wide association studies and genomic prediction in the Saccharum gene bank are helpful to identify the complete genomic resources for cane yield, commercial cane sugar, tolerances to biotic and abiotic stresses, and other agronomic traits. The implementation of genomic selection in pre-breeding, genetic base-broadening programmes assist in precise introgression of specific genes and recurrent selection schemes enhance the higher frequency of favorable alleles in the population with a considerable reduction in breeding cycles and population size. The integration of environmental covariates and genomic prediction in multi-environment trials assists in the prediction of varietal performance for different agro-climatic zones. This review also directed its focus on enhancing the genetic gain over time, cost, and resource allocation at various stages of breeding programmes.

Overdominance at the Gene Expression Level Plays a Critical Role in the Hybrid Root Growth of Brassica napus

Despite heterosis contributing to genetic improvements in crops, root growth heterosis in rapeseed plants is poorly understood at the molecular level. The current study was performed to discover key differentially expressed genes (DEGs) related to heterosis in two hybrids with contrasting root growth performance (FO; high hybrid and FV; low hybrid) based on analysis of the root heterosis effect. Based on comparative transcriptomic analysis, we believe that the overdominance at the gene expression level plays a critical role in hybrid roots’ early biomass heterosis. Our findings imply that a considerable increase in up-regulation of gene expression underpins heterosis. In the FO hybrid, high expression of DEGs overdominant in the starch/sucrose and galactose metabolic pathways revealed a link between hybrid vigor and root growth. DEGs linked to auxin, cytokinin, brassinosteroids, ethylene, and abscisic acid were also specified, showing that these hormones may enhance mechanisms of root growth and the development in the FO hybrid. Moreover, transcription factors such as MYB, ERF, bHLH, NAC, bZIP, and WRKY are thought to control downstream genes involved in root growth. Overall, this is the first study to provide a better understanding related to the regulation of the molecular mechanism of heterosis, which assists in rapeseed growth and yield improvement.

Genome-wide dissection of hybridization for fiber quality- and yield-related traits in upland cotton

An evaluation of combining ability can facilitate the selection of suitable parents and superior F₁ hybrids for hybrid cotton breeding, although the molecular genetic basis of combining ability has not been fully characterized. In the present study, 282 female parents were crossed with four male parents in accordance with the North Carolina II mating scheme to generate 1128 hybrids. The parental lines were genotyped based on restriction site-associated DNA sequencing and 306 814 filtered single nucleotide polymorphisms were used for genome-wide association analysis involving the phenotypes, general combining ability (GCA) values, and specific combining ability values of eight fiber quality- and yield-related traits. The main results were: (i) all parents could be clustered into five subgroups based on population structure analyses and the GCA performance of the female parents had significant differences between subgroups; (ii) 20 accessions with a top 5% GCA value for more than one trait were identified as elite parents for hybrid cotton breeding; (iii) 120 significant single nucleotide polymorphisms, clustered into 66 quantitative trait loci, such as the previously reported Gh_A07G1769 and GhHOX3 genes, were found to be significantly associated with GCA; and (iv) identified quantitative trait loci for GCA had a cumulative effect on GCA of the accessions. Overall, our results suggest that pyramiding the favorable loci for GCA may improve the efficiency of hybrid cotton breeding.

Recommendations for Choosing the Genotyping Method and Best Practices for Quality Control in Crop Genome-Wide Association Studies

High-throughput genotyping boosts genome-wide association studies (GWAS) in crop species, leading to the identification of single-nucleotide polymorphisms (SNPs) associated with economically important traits. Choosing a cost-effective genotyping method for crop GWAS requires careful examination of several aspects, namely, the purpose and the scale of the study, crop-specific genomic features, and technical and economic matters associated with each genotyping option. Once genotypic data have been obtained, quality control (QC) procedures must be applied to avoid bias and false signals in genotype–phenotype association tests. QC for human GWAS has been extensively reviewed; however, QC for crop GWAS may require different actions, depending on the GWAS population type. Here, we review most popular genotyping methods based on next-generation sequencing (NGS) and array hybridization, and report observations that should guide the investigator in the choice of the genotyping method for crop GWAS. We provide recommendations to perform QC in crop species, and deliver an overview of bioinformatics tools that can be used to accomplish all needed tasks. Overall, this work aims to provide guidelines to harmonize those procedures leading to SNP datasets ready for crop GWAS.

Identification and fine mapping of a major locus controlling branching in Brassica napus

A candidate branching-controlling gene for qDBA09 was identified after delimiting a Brassica napus recessive locus within a 270-kb interval on chromosome A09. Although branching is an important trait associated with the adaptation and yield potential of rapeseed (Brassica napus), the genetic mechanisms underlining branching in this crop remain poorly understood. In this study, we characterized a naturally occurring rapeseed mutant, db1, which showed an ultrahigh branching density phenotype. By combining bulked segregant analysis (BSA) and the Brassica 60K SNP BeadChip Array, we identified two major quantitative trait loci (QTLs), qDBA09 and qDBC06, which were subsequently confirmed using the traditional QTL-mapping approach. Analysis of 208 individuals from a BC₁F₃ population indicated that the qDBA09 locus is a single Mendelian factor and that the dense branching phenotype is controlled by a single recessive gene. Furthermore, QTL analysis confirmed that qDBA09 explained between 9.5 and 70.5% of the variation in branching-related traits. Using 7785 individuals from the BC₁F₃ population, we mapped qDBA09 to a DNA fragment of approximately 270 kb in length that contained 27 predicted genes, three of which were identified as potentially involved in the control of the dense branching trait. Based on the reported function of these genes, together with sequence comparisons and co-segregation analysis, we identified a potential candidate gene for the qDBA09 locus. The present findings lay the foundations for further in-depth research on the branching mechanisms of B. napus.

Potential of the C Genome of the Different Variants of Brassica oleracea for Heterosis in Spring B. napus Canola

The genetic base of Brassica napus canola need to be broadened for exploitation of heterosis at a greater level in the breeding of F1 hybrid canola cultivars. In this study, we evaluated 228 inbred B. napus canola lines derived from six B. napus × B. oleracea interspecific crosses and following two breeding methods (F2- and BC1-derived lines) to understand the effect of the B. oleracea alleles on heterosis for different agronomic and seed quality traits. Test hybrids of the inbreds derived from crosses involving vars. botrytis (cauliflower), alboglabra (Chinese kale) and capitata (cabbage) cv. Badger Shipper, on an average, gave about 10% mid-parent heterosis (MPH), and about 67% of the test hybrids gave higher seed yield than the common B. napus parent indicating that B. oleracea alleles can contribute to heterosis for seed yield in spring B. napus canola hybrids. This was also evident from a positive correlation of the genetic distance of the inbred lines from the common B. napus parent with MPH for seed yield (r = 0.31) as well as with hybrid yield (r = 0.26). Almost no correlation was found between genetic distance and MPH for seed oil and protein content as well as with the performance of the test hybrids for these two traits. The occurrence of positive correlation between seed yield of the inbred lines and test hybrids suggested the importance of the genes exerting additive effect for high seed yield in the hybrids. Very little or almost no heterosis was found for the other agronomic traits as well as for seed oil and protein content. While comparing the two breeding methods, no significant difference was found for seed yield of the test hybrids or the level of MPH; however, the BC1-derived inbred and test hybrid populations flowered and matured earlier and had longer grain-filling period than the F2-derived population. Thus, the results suggested that the B. oleracea gene pool can be used in the breeding of spring B. napus canola to improve seed yield in hybrid cultivars.

Fine mapping of a silique length- and seed weight-related gene in Brassica napus

Using microarray analysis combined with map-based cloning, a major locus positively regulating SL and SW was mapped to a 98.47 kb interval on A09 in rapeseed. In rapeseed, seed yield is closely associated with silique-related traits such as silique length (SL) and seed weight (SW). Previously identified quantitative trait loci (QTLs) revealed that SL and SW are complex traits and many QTLs overlap. However, the genetic characterization of the association between SL and SW is poorly understood. In the present study, a BC₃F₃ near isogenic line developed from a short silique plant and the long silique cultivar 'ZS11' was analyzed to identify the locus related to SL. Map-based cloning indicated that a major locus acting as a single Mendelian factor was mapped to a 98.47 kb region on chromosome A09. BLAST analysis and DNA sequencing showed SNP variations and a fragment replacement in the upstream region of the candidate gene BnaA09g55530D may alter gene expression and influence SL. The results showed that this SL locus may also positively affect SW as well as in the 186 rapeseed accessions identified by the associated markers. Therefore, selecting plants with appropriate SL and developing functional markers for the associated gene could play important roles in the molecular breeding of high-yield rapeseed varieties.

B. K, Satyavathi CT, Kumar S, Singh RB, Yadav RS, Gupta R, Srivastava RK. Discerning combining ability loci for divergent environments using chromosome segment substitution lines (CSSLs) in pearl millet

Pearl millet is an important crop for arid and semi-arid regions of the world. Genomic regions associated with combining ability for yield-related traits under irrigated and drought conditions are useful in heterosis breeding programs. Chromosome segment substitution lines (CSSLs) are excellent genetic resources for precise QTL mapping and identifying naturally occurring favorable alleles. In the present study, testcross hybrid populations of 85 CSSLs were evaluated for 15 grain and stover yield-related traits for summer and wet seasons under irrigated control (CN) and moisture stress (MS) conditions. General combining ability (GCA) and specific combining ability (SCA) effects of all these traits were estimated and significant marker loci linked to GCA and SCA of the traits were identified. Heritability of the traits ranged from 53-94% in CN and 63-94% in MS. A total of 40 significant GCA loci and 36 significant SCA loci were identified for 14 different traits. Five QTLs (flowering time, panicle number and panicle yield linked to Xpsmp716 on LG4, flowering time and grain number per panicle with Xpsmp2076 on LG4) simultaneously controlled both GCA and SCA, demonstrating their unique genetic basis and usefulness for hybrid breeding programs. This study for the first time demonstrated the potential of a set of CSSLs for trait mapping in pearl millet. The novel combining ability loci linked with GCA and SCA values of the traits identified in this study may be useful in pearl millet hybrid and population improvement programs using marker-assisted selection (MAS).

Gene action and heterosis in F1 clonal progenies of cassava for β-Carotene and farmers' preferred traits

Gene action and heterosis provides information to assist breeder for selecting and generating improved plant recombinants. This study aimed at determining the gene action of selected cassava traits. The F₁ clones exhibited considerable phenotypic variability between families and offsprings. The best F₁ progenies had a higher amount of β-carotene (β-C) of 6.12 mg 100 g⁻¹ against 1.32 mg 100 g⁻¹ of the best parent. This superiority could be attributed to the over-dominance from the recombination of additive gene action and epistasis. The general combining ability (GCA) of parents and specific combining ability (SCA) of combinations were significant for different traits, and indicating the role of additive and non-additive gene action in controlling such traits. The significant GCA for β-C and postharvest physiological deterioration (PPD) indicates the role of additive gene action. The significant SCA for cassava mosaic disease (CMD) and cassava brown streak disease (CBSD) showed a predominance of non-additive gene action. The F₁ progenies from the family Mavoka x Garukunsubire expressed the highest positive heterosis for CMD, dry matter and β-C. The high positive heterosis for β-C and DMC could be linked to transgressive segregation, because one of the parents was poor combiner.

Transcriptome and organellar sequencing highlights the complex origin and diversification of allotetraploid Brassica napus

Brassica napus, an allotetraploid crop, is hypothesized to be a hybrid from unknown varieties of Brassica rapa and Brassica oleracea. Despite the economic importance of B. napus, much is unresolved regarding its phylogenomic relationships, genetic structure, and diversification. Here we conduct a comprehensive study among diverse accessions from 183 B. napus (including rapeseed, rutabaga, and Siberian kale), 112 B. rapa, and 62 B. oleracea and its wild relatives. Using RNA-seq of B. napus accessions, we define the genetic diversity and sub-genome variance of six genetic clusters. Nuclear and organellar phylogenies for B. napus and its progenitors reveal varying patterns of inheritance and post-formation introgression. We discern regions with signatures of selective sweeps and detect 8,187 differentially expressed genes with implications for B. napus diversification. This study highlights the complex origin and evolution of B. napus providing insights that can further facilitate B. napus breeding and germplasm preservation.

Association analysis uncovers the genetic basis of general combining ability of 11 yield-related traits in parents of hybrid rice

Association analysis between constructed single nucleotide polymorphism linkage disequilibrium blocks (SNPLDBs) and general combining ability (GCA) effects is a novel approach to uncover the genetic basis of GCA within the sequence genomes of parents of hybrid rice. Here, we calculated the GCA effect values of 33 parents of hybrid rice and sequenced them to identify genome-wide single nucleotide polymorphisms (SNPs). In total, 64.6 % of the uniquely mapped paired-end short reads revealed a final total of 291 959 SNPs between the 33 parental genomes and the Nipponbare reference genome. The identified SNPs were non-randomly distributed among all chromosomes of rice, whereas one-fourth of the SNPs were situated in the exonic regions with 16 % being non-synonymous. Further, the identified SNPs were merged and optimized for construction of 2612 SNPLDB markers, using linkage disequilibrium information. The single-factor analysis of variance-based association method between the constructed SNPLDB markers and GCA effects values detected 99 significant SNPLDBs for GCA of 11 yield-related traits. The associated SNPLDB markers explained 26.4 % of phenotypic variations with traits, on average. We mined 50 favourable GCA alleles at the associated SNPLDBs regions, distributed across the 33 parental genomes. The parental genomes possessed a small number of favourable GCA alleles for studied traits, with the exception of days to heading and plant height. Our results suggest that the identified GCA alleles could be used to improve the GCA performance of parents of hybrid rice through optimal crossing design. Moreover, favourable GCA alleles should be incorporated in the parental genomes through marker-assisted selection experiments, and the parental lines carrying more alleles could be utilized in breeding as superior parents for developing rice hybrids of desirable characteristics.

Assembly and comparison of two closely related Brassica napus genomes

As an increasing number of plant genome sequences become available, it is clear that gene content varies between individuals, and the challenge arises to predict the gene content of a species. However, genome comparison is often confounded by variation in assembly and annotation. Differentiating between true gene absence and variation in assembly or annotation is essential for the accurate identification of conserved and variable genes in a species. Here, we present the de novo assembly of the B. napus cultivar Tapidor and comparison with an improved assembly of the Brassica napus cultivar Darmor-bzh. Both cultivars were annotated using the same method to allow comparison of gene content. We identified genes unique to each cultivar and differentiate these from artefacts due to variation in the assembly and annotation. We demonstrate that using a common annotation pipeline can result in different gene predictions, even for closely related cultivars, and repeat regions which collapse during assembly impact whole genome comparison. After accounting for differences in assembly and annotation, we demonstrate that the genome of Darmor-bzh contains a greater number of genes than the genome of Tapidor. Our results are the first step towards comparison of the true differences between B. napus genomes and highlight the potential sources of error in future production of a B. napus pangenome.