Microsatellite markers used for genome-wide association mapping of partial resistance to Sclerotinia sclerotiorum in a world collection of Brassica napus

Development of simple sequence repeat markers for sugarcane from data mining of expressed sequence tags

Sugarcane (Saccharum spp. hybrids) is a worldwide acclaimed important agricultural crop used primarily for sugar production and biofuel. Sugarcane's genetic complexity, aneuploidy, and extreme heterozygosity make it a challenging crop in developing improved varieties. The molecular breeding programs promise to develop nutritionally improved varieties for both direct consumption and commercial application. Therefore, to address these challenges, the development of simple sequence repeats (SSRs) has been proven to be a powerful molecular tool in sugarcane. This study involved the collection of 285216 expressed sequence tags (ESTs) from sugarcane, resulting in 23666 unigenes, including 4547 contigs. Our analysis identified 4120 unigenes containing a total of 4960 SSRs, with the most abundant repeat types being monomeric (44.33%), dimeric (13.10%), and trimeric (39.68%). We further chose 173 primers to analyze the banding pattern in 10 sugarcane accessions by PAGE analysis. Additionally, functional annotation analysis showed that 71.07%, 53.6%, and 10.3% unigenes were annotated by Uniport, GO, and KEGG, respectively. GO annotations and KEGG pathways were distributed across three functional categories: molecular (46.46%), cellular (33.94%), and biological pathways (19.6%). The cluster analysis indicated the formation of four distinct clusters among selected sugarcane accessions, with maximum genetic distance observed among the varieties. We believe that these EST-SSR markers will serve as valuable references for future genetic characterization, species identification, and breeding efforts in sugarcane.

Genetic and molecular analysis of stem rot (Sclerotinia sclerotiorum) resistance in Brassica napus (canola type)

Identifying the molecular and genetic basis of resistance to Sclerotinia stem rot (Sclerotinia sclerotiorum) is critical for developing long-term and cost-effective management of this disease in rapeseed/canola (Brassica napus). Current cultural or chemical management options provide, at best, only partial and/or sporadic control. Towards this, a B. napus breeding population (Mystic x Rainbow), including the parents, F1, F2, BC1P1 and BC1P2, was utilized in a field study to determine the inheritance pattern of Sclerotinia stem rot resistance (based on stem lesion length, SLL). Broad sense heritability was 0.58 for SLL and 0.44 for days to flowering (DTF). There was a significant negative correlation between SLL and stem diameter (SD) (r = -0.39) and between SLL and DTF (r = -0.28), suggesting co-selection of SD and DTF traits, along with SLL, should assist in improving overall resistance. Non-additive genetic variance was evident for SLL, DTF, and SD. In a genome wide association study (GWAS), a significant quantitative trait locus (QTL) was identified for SLL. Several putative candidate marker trait associations (MTA) were located within this QTL region. Overall, this study has provided valuable new understanding of inheritance of resistance to S. sclerotiorum, and has identified QTL, MTAs and transgressive segregants with high-level resistances. Together, these will foster more rapid selection for multiple traits associated with Sclerotinia stem rot resistance, by enabling breeders to make critical choices towards selecting/developing cultivars with enhanced resistance to this devastating pathogen.

Fungal effectors versus defense-related genes of B. juncea and the status of resistant transgenics against fungal pathogens

Oilseed brassica has become instrumental in securing global food and nutritional security. B. juncea, colloquially known as Indian mustard, is cultivated across tropics and subtropics including Indian subcontinent. The production of Indian mustard is severely hampered by fungal pathogens which necessitates human interventions. Chemicals are often resorted to as they are quick and effective, but due to their economic and ecological unsustainability, there is a need to explore their alternatives. The B. juncea-fungal pathosystem is quite diverse as it covers broad-host range necrotrophs (Sclerotinia sclerotiorum), narrow-host range necrotrophs (Alternaria brassicae and A. brassicicola) and biotrophic oomycetes (Albugo candida and Hyaloperonospora brassica). Plants ward off fungal pathogens through two-step resistance mechanism; PTI which involves recognition of elicitors and ETI where the resistance gene (R gene) interacts with the fungal effectors. The hormonal signalling is also found to play a vital role in defense as the JA/ET pathway is initiated at the time of necrotroph infection and SA pathway is induced when the biotrophs attack plants. The review discuss the prevalence of fungal pathogens of Indian mustard and the studies conducted on effectoromics. It covers both pathogenicity conferring genes and host-specific toxins (HSTs) that can be used for a variety of purposes such as identifying cognate R genes, understanding pathogenicity and virulence mechanisms, and establishing the phylogeny of fungal pathogens. It further encompasses the studies on identifying resistant sources and characterisation of R genes/quantitative trait loci and defense-related genes identified in Brassicaceae and unrelated species which, upon introgression or overexpression, confer resistance. Finally, the studies conducted on developing resistant transgenics in Brassicaceae have been covered in which chitinase and glucanase genes are mostly used. The knowledge gained from this review can further be used for imparting resistance against major fungal pathogens.

QTL mapping and transcriptome analysis identify novel QTLs and candidate genes in Brassica villosa for quantitative resistance against Sclerotinia sclerotiorum

Novel QTLs and candidate genes for Sclerotinia-resistance were identified in B. villosa, a wild Brassica species, which represents a new genetic source for improving oilseed rape resistance to SSR. Sclerotinia stem rot (SSR), caused by Sclerotinia sclerotiorum, is one of the most destructive diseases in oilseed rape growing regions. To date, there is no effective genetic resistance against S. sclerotiorum in the B. napus germplasm and knowledge of the molecular plant-fungal interaction is also limited. To identify new resistance resources, we screened a set of wild Brassica species and identified B. villosa (BRA1896) with a high level of Sclerotinia-resistance. Two segregating F₂ populations for Sclerotinia-resistance, generated by interspecific crosses between the resistant B. villosa (BRA1896) and the wild susceptible B. oleracea (BRA1909) were assessed for Sclerotinia-resistance. Genetic mapping using a 15-k Illumina Infinium SNP-array resulted in a high-density genetic map containing 1,118 SNP markers and spanning a total genetic length of 792.2 cM. QTL analysis revealed seven QTLs explaining 3.8% to 16.5% of phenotypic variance. Intriguingly, RNAseq-based transcriptome analysis identified genes and pathways specific to B. villosa, of which a cluster of five genes encoding putative receptor-like kinases (RLKs) and two pathogenesis-related (PR) proteins are co-localized within a QTL on chromosome C07. Furthermore, transcriptomic analysis revealed enhanced ethylene (ET)-activated signaling in the resistant B. villosa, which is associated with a stronger plant immune response, depressed cell death, and enhanced phytoalexin biosynthesis compared to the susceptible B. oleracea. Our data demonstrates that B. villosa represents a novel and unique genetic source for improving oilseed rape resistance against SSR.

Characterization of Simple Sequence Repeat (SSR) Markers Mined in Whole Grape Genomes

SSR (simple sequence repeat) DNA markers are widely used for genotype DNA identification, QTL mapping, and analyzing genetic biodiversity. However, SSRs in grapes are still in their early stages, with a few primer pairs accessible. With the whole-genome sequencing (WGS) of several grape varieties, characterization of grape SSR changed to be necessary not only to genomics but to also help SSR development and utility. Based on this, we identified the whole-genome SSR of nine grape cultivars (‘PN40024’, ‘Cabernet Sauvignon’, ‘Carménère’, ‘Chardonnay’, ‘Merlot’, ‘Riesling’, ‘Zinfandel’, ‘Shine Muscat’, and ‘Muscat Hamburg’) with whole-genome sequences released publicly and found that there are great differences in the distribution of SSR loci in different varieties. According to the difference in genome size, the number of SSRs ranged from 267,385 (Cabernet Sauvignon) to 627,429 (Carménère), the density of the SSR locus in the genome of nine cultivars was generally 1 per Kb. SSR motif distribution characteristic analysis of these grape cultivars showed that the distribution patterns among grape cultivars were conservative, mainly enriched in A/T. However, there are some differences in motif types (especially tetranucleotides, pentanucleotides, and hexanucleotides), quantity, total length, and average length in different varieties, which might be related to the size of the assembled genome or the specificity of variety domestication. The distribution characteristics of SSRs were revealed by whole-genome analysis of simple repeats of grape varieties. In this study, 32 pairs of primers with lower polymorphism have been screened, which provided an important research foundation for the development of molecular markers of grape variety identification and the construction of linkage maps of important agronomic traits for crop improvement.

Genetic diversity and population structure assessment of Western Canadian barley cooperative trials

Studying the population structure and genetic diversity of historical datasets is a proposed use for association analysis. This is particularly important when the dataset contains traits that are time-consuming or costly to measure. A set of 96 elite barley genotypes, developed from eight breeding programs of the Western Canadian Cooperative Trials were used in the current study. Genetic diversity, allelic variation, and linkage disequilibrium (LD) were investigated using 5063 high-quality SNP markers via the Illumina 9K Barley Infinium iSelect SNP assay. The distribution of SNPs markers across the barley genome ranged from 449 markers on chromosome 1H to 1111 markers on chromosome 5H. The average polymorphism information content (PIC) per locus was 0.275 and ranged from 0.094 to 0.375. Bayesian clustering in STRUCTURE and principal coordinate analysis revealed that the populations are differentiated primarily due to the different breeding program origins and ear-row type into five subpopulations. Analysis of molecular variance based on PhiPT values suggested that high values of genetic diversity were observed within populations and accounted for 90% of the total variance. Subpopulation 5 exhibited the most diversity with the highest values of the diversity indices, which represent the breeding program gene pool of AFC, AAFRD, AU, and BARI. With increasing genetic distance, the LD values, expressed as r², declined to below the critical r^{2 =} 0.18 after 3.91 cM, and the same pattern was observed on each chromosome. Our results identified an important pattern of genetic diversity among the Canadian barley panel that was proposed to be representative of target breeding programs and may have important implications for association mapping in the future. This highlight, that efforts to identify novel variability underlying this diversity may present practical breeding opportunities to develop new barley genotypes.

The evolutionary and molecular features of the broad-host-range plant pathogen Sclerotinia sclerotiorum

Sclerotinia sclerotiorum is a pathogenic fungus that infects hundreds of plant species, including many of the world's most important crops. Key features of S. sclerotiorum include its extraordinary host range, preference for dicotyledonous plants, relatively slow evolution, and production of protein effectors that are active in multiple host species. Plant resistance to this pathogen is highly complex, typically involving numerous polymorphisms with infinitesimally small effects, which makes resistance breeding a major challenge. Due to its economic significance, S. sclerotiorum has been subjected to a large amount of molecular and evolutionary research. In this updated pathogen profile, we review the evolutionary and molecular features of S. sclerotiorum and discuss avenues for future research into this important species.

Genetic mapping and genomic prediction of sclerotinia stem rot resistance to rapeseed/canola (Brassica napus L.) at seedling stage

GWAS detected ninety-eight significant SNPs associated with Sclerotinia sclerotiorum resistance. Six statistical models resulted in medium to high predictive ability, depending on trait, indicating potential of genomic prediction for disease resistance breeding. The lack of complete host resistance and a complex resistance inheritance nature between rapeseed/canola and Sclerotinia sclerotiorum often limits the development of functional molecular markers that enable breeding for sclerotinia stem rot (SSR) resistance. However, genomics-assisted selection has the potential to accelerate the breeding for SSR resistance. Therefore, genome-wide association (GWA) mapping and genomic prediction (GP) were performed using a diverse panel of 337 rapeseed/canola genotypes. Three-week-old seedlings were screened using the petiole inoculation technique (PIT). Days to wilt (DW) up to 2 weeks and lesion phenotypes (LP) at 3, 4, and 7 days post-inoculation (dpi) were recorded. A strong correlation (r = - 0.90) between DW and LP_4dpi implied that a single time point scoring at four days could be used as a proxy trait. GWA analyses using single-locus (SL) and multi-locus (ML) models identified a total of 41, and 208 significantly associated SNPs, respectively. Out of these, ninety-eight SNPs were identified by a combination of the SL model and any of the ML models, at least two ML models, or two traits. These SNPs explained 1.25-12.22% of the phenotypic variance and considered as significant, could be associated with SSR resistance. Eighty-three candidate genes with a function in disease resistance were associated with the significant SNPs. Six GP models resulted in moderate to high (0.42-0.67) predictive ability depending on SSR resistance traits. The resistant genotypes and significant SNPs will serve as valuable resources for future SSR resistance breeding. Our results also highlight the potential of genomic selection to improve rapeseed/canola breeding for SSR resistance.

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.

Sources of genomic diversity in the self-fertile plant pathogen, Sclerotinia sclerotiorum, and consequences for resistance breeding

The ascomycete, Sclerotinia sclerotiorum, has a broad host range and causes yield loss in dicotyledonous crops world wide. Genomic diversity was determined in a population of 127 isolates obtained from individual canola (Brassica napus) fields in western Canada. Genotyping with 39 simple sequence repeat (SSR) markers revealed each isolate was a unique haplotype. Analysis of molecular variance showed 97% was due to isolate and 3% due to geographical location. Testing of mycelium compatibility among 133 isolates identified clones of mutually compatible isolates with 86-95% similar SSR haplotype, whereas incompatible isolates were highly diverse. In the Province of Manitoba, 61% of isolates were compatible forming clones and stings of pairwise compatible isolates not described before. In contrast, only 35% of isolates were compatible in Alberta without forming clones and strings, while 39% were compatible in Saskatchewan with a single clone, but no strings. These difference can be explained by wetter growing seasons and more susceptible crop species in Manitoba favouring frequent mycelium interaction and more life cycles over time, which might also explain similar differences observed in other geographical areas and host crops. Analysis of linkage disequilibrium rejected random recombination, consistent with a self-fertile fungus, restricted outcrossing due to mycelium incompatibility, and only a single annual opportunity for genomic recombination during meiosis in the ascospore stage between non-sister chromatids in the rare event nuclei from different isolates come together. More probable sources of genomic diversity is slippage during DNA replication and point mutation affecting single nucleotides that accumulate and likely increase mycelium incompatibility in a population over time. A phylogenetic tree based on SSR haplotype grouped isolates into 17 sub-populations. Aggressiveness was tested by inoculating one isolate from each sub-population onto B. napus lines with quantitative resistance. Analysis of variance was significant for isolate, line, and isolate by line interaction. These isolates represent the genomic and pathogenic diversity in western Canada, and are suitable for resistance screening in canola breeding programs.

Quantitative Inheritance of Sclerotinia Stem Rot Resistance in Brassica napus and Relationship to Cotyledon and Leaf Resistances

Sclerotinia sclerotiorum is a necrotrophic fungus causing devastating stem rot and associated yield losses of canola/rapeseed (Brassica napus) worldwide, including in Australia. Developing host resistance against Sclerotinia stem rot is critical if this disease in canola/rapeseed is to be successfully managed, as cultural or chemical control options provide only partial or sporadic control. Three B. napus breeding populations, C2, C5 and C6, including the parents, F₁, F₂, BC₁P1, and BC₂P2, were used in a field study with an objective of exploring the inheritance pattern of disease resistance (based on stem lesion length [SLL]) and the genetic relationships of disease with stem diameter (SD) or days to first flowering (DTF), and to compare these new adult plant stem resistances against S. sclerotiorum with those of seedling (cotyledon and leaf) resistances in earlier studies. Heritability (broad sense) of SLL was 0.57 and 0.73 for population C2 at 3 and 5 weeks postinoculation and 0.21 for population C5 at 5 weeks postinoculation. Additive genetic variance was evident within all 3 populations for DTF but not for SD. Narrow-sense heritability for DTF was 0.48 (C2), 0.42 (C5), and 0.32 (C6). SD, DTF, and SLL were all inherited independently, with no significant genetic covariance between traits in bivariate analysis. Genetic variance for SLL in populations C2 and C5 was entirely nonadditive, and there was significant nonadditive genetic covariance of SLL at 3 and 5 weeks postinoculation. Generation means analysis in population C2 supported the conclusion that complex epistatic interactions controlled SLL. Several C2 and C5 progeny showed high adult plant stem resistance, which may be critical in developing enhanced stem resistance in canola/rapeseed. Although population C6 showed no genetic variation for SLL resistance in this study, it showed significant nonadditive genetic variance at the cotyledon and leaf stages in earlier studies. We conclude that host resistance varies across different plant growth stages, and breeding must be targeted for resistance at each growth stage. In populations C2, C5, and C6, resistance to S. sclerotiorum in stem, leaf, and cotyledon was always controlled by nonadditive effects such as complex epistasis or dominance. Overall, our findings in relation to the quantitative inheritance of Sclerotinia stem rot resistance, together with the new high-level resistances identified, will enable breeders to select/develop genotypes with enhanced resistances to S. sclerotiorum.

Genome-wide association mapping and genomic prediction for adult stage sclerotinia stem rot resistance in Brassica napus (L) under field environments

Sclerotinia stem rot (SSR) is a fungal disease of rapeseed/canola that causes significant seed yield losses and reduces its oil content and quality. In the present study, the reaction of 187 diverse canola genotypes to SSR was characterized at full flowering stage using the agar plug to stem inoculation method in four environments. Genome-wide association study (GWAS) using three different algorithms identified 133 significant SNPs corresponding with 123 loci for disease traits like stem lesion length (LL), lesion width (LW), and plant mortality at 14 (PM_14D) and 21 (PM_21D) days. The explained phenotypic variation of these SNPs ranged from 3.6 to 12.1%. Nineteen significant SNPs were detected in two or more environments, disease traits with at least two GWAS algorithms. The strong correlations observed between LL and other three disease traits evaluated, suggest they could be used as proxies for SSR resistance phenotyping. Sixty-nine candidate genes associated with disease resistance mechanisms were identified. Genomic prediction (GP) analysis with all the four traits employing genome-wide markers resulted in 0.41-0.64 predictive ability depending on the model specifications. The highest predictive ability for PM_21D with three models was about 0.64. From our study, the identified resistant genotypes and stable significant SNP markers will serve as a valuable resource for future SSR resistance breeding. Our study also suggests that genomic selection holds promise for accelerating canola breeding progress by enabling breeders to select SSR resistance genotypes at the early stage by reducing the need to phenotype large numbers of genotypes.

Modeling first order additive × additive epistasis improves accuracy of genomic prediction for sclerotinia stem rot resistance in canola

The fungus Sclerotinia sclerotiorum infects hundreds of plant species including many crops. Resistance to this pathogen in canola (Brassica napus L. subsp. napus) is controlled by numerous quantitative trait loci (QTL). For such polygenic traits, genomic prediction may be useful for breeding as it can capture many QTL at once while also considering nonadditive genetic effects. Here, we test application of common regression models to genomic prediction of S. sclerotiorum resistance in canola in a diverse panel of 218 plants genotyped at 24,634 loci. Disease resistance was scored by infection with an aggressive isolate and monitoring over 3 wk. We found that including first-order additive × additive epistasis in linear mixed models (LMMs) improved accuracy of breeding value estimation between 3 and 40%, depending on method of assessment, and correlation between phenotypes and predicted total genetic values by 14%. Bayesian models performed similarly to or worse than genomic relationship matrix-based models for estimating breeding values or overall phenotypes from genetic values. Bayesian ridge regression, which is most similar to the genomic relationship matrix-based approach in the amount of shrinkage it applies to marker effects, was the most accurate of this family of models. This confirms several studies indicating the highly polygenic nature of sclerotinia stem rot resistance. Overall, our results highlight the use of simple epistasis terms for prediction of breeding values and total genetic values for a complex disease resistance phenotype in canola.

Genome-wide simple sequence repeats (SSR) markers discovered from whole-genome sequence comparisons of multiple spinach accessions

The availability of well-assembled genome sequences and reduced sequencing costs have enabled the resequencing of many additional accessions in several crops, thus facilitating the rapid discovery and development of simple sequence repeat (SSR) markers. Although the genome sequence of inbred spinach line Sp75 is available, previous efforts have resulted in a limited number of useful SSR markers. Identification of additional polymorphic SSR markers will support genetics and breeding research in spinach. This study aimed to use the available genomic resources to mine and catalog a large number of polymorphic SSR markers. A search for SSR loci on six chromosome sequences of spinach line Sp75 using GMATA identified a total of 42,155 loci with repeat motifs of two to six nucleotides in the Sp75 reference genome. Whole-genome sequences (30x) of additional 21 accessions were aligned against the chromosome sequences of the reference genome and in silico genotyped using the HipSTR program by comparing and counting repeat numbers variation across the SSR loci among the accessions. The HipSTR program generated SSR genotype data were filtered for monomorphic and high missing loci, and a final set of the 5986 polymorphic SSR loci were identified. The polymorphic SSR loci were present at a density of 12.9 SSRs/Mb and were physically mapped. Out of 36 randomly selected SSR loci for validation, two failed to amplify, while the remaining were all polymorphic in a set of 48 spinach accessions from 34 countries. Genetic diversity analysis performed using the SSRs allele score data on the 48 spinach accessions showed three main population groups. This strategy to mine and develop polymorphic SSR markers by a comparative analysis of the genome sequences of multiple accessions and computational genotyping of the candidate SSR loci eliminates the need for laborious experimental screening. Our approach increased the efficiency of discovering a large set of novel polymorphic SSR markers, as demonstrated in this report.

Genome-wide simple sequence repeats (SSR) markers discovered from whole-genome sequence comparisons of multiple spinach accessions

The availability of well-assembled genome sequences and reduced sequencing costs have enabled the resequencing of many additional accessions in several crops, thus facilitating the rapid discovery and development of simple sequence repeat (SSR) markers. Although the genome sequence of inbred spinach line Sp75 is available, previous efforts have resulted in a limited number of useful SSR markers. Identification of additional polymorphic SSR markers will support genetics and breeding research in spinach. This study aimed to use the available genomic resources to mine and catalog a large number of polymorphic SSR markers. A search for SSR loci on six chromosome sequences of spinach line Sp75 using GMATA identified a total of 42,155 loci with repeat motifs of two to six nucleotides in the Sp75 reference genome. Whole-genome sequences (30x) of additional 21 accessions were aligned against the chromosome sequences of the reference genome and in silico genotyped using the HipSTR program by comparing and counting repeat numbers variation across the SSR loci among the accessions. The HipSTR program generated SSR genotype data were filtered for monomorphic and high missing loci, and a final set of the 5986 polymorphic SSR loci were identified. The polymorphic SSR loci were present at a density of 12.9 SSRs/Mb and were physically mapped. Out of 36 randomly selected SSR loci for validation, two failed to amplify, while the remaining were all polymorphic in a set of 48 spinach accessions from 34 countries. Genetic diversity analysis performed using the SSRs allele score data on the 48 spinach accessions showed three main population groups. This strategy to mine and develop polymorphic SSR markers by a comparative analysis of the genome sequences of multiple accessions and computational genotyping of the candidate SSR loci eliminates the need for laborious experimental screening. Our approach increased the efficiency of discovering a large set of novel polymorphic SSR markers, as demonstrated in this report.

Sclerotinia Stem Rot Resistance in Rapeseed: Recent Progress and Future Prospects

Sclerotinia stem rot (SSR) of rapeseed (Brassica napus), caused by the soil-borne fungus Sclerotinia sclerotiorum, is one of the main diseases seriously affecting the yield and oil quality of infected rapeseed crops. The complexity of the inheritance of resistance and of the interaction mechanisms between rapeseed and S. sclerotiorum limits resistance gene identification and molecular breeding. In this review, the latest progress of research into resistance to SSR in B. napus is summarized from the following three directions: the pathogenesis mechanisms of S. sclerotiorum, the resistance mechanisms of B. napus toward S. sclerotiorum, and rapeseed breeding for resistance to SSR. This review aims to provide a theoretical basis and useful reference for analyzing the mechanism of the interaction between B. napus and S. sclerotiorum, searching for gene loci associated with the resistance response, and for achieving disease-resistance genetic manipulation and molecular design breeding in rapeseed.

Current Status of the Disease-Resistant Gene(s)/QTLs, and Strategies for Improvement in Brassica juncea

Brassica juncea is a major oilseed crop in tropical and subtropical countries, especially in south-east Asia like India, China, Bangladesh, and Pakistan. The widespread cultivation of genetically similar varieties tends to attract fungal pathogens which cause heavy yield losses in the absence of resistant sources. The conventional disease management techniques are often expensive, have limited efficacy, and cause additional harm to the environment. A substantial approach is to identify and use of resistance sources within the Brassica hosts and other non-hosts to ensure sustainable oilseed crop production. In the present review, we discuss six major fungal pathogens of B. juncea: Sclerotinia stem rot (Sclerotinia sclerotiorum), Alternaria blight (Alternaria brassicae), White rust (Albugo candida), Downy mildew (Hyaloperonospora parasitica), Powdery mildew (Erysiphe cruciferarum), and Blackleg (Leptoshaeria maculans). From discussing studies on pathogen prevalence in B. juncea, the review then focuses on highlighting the resistance sources and quantitative trait loci/gene identified so far from Brassicaceae and non-filial sources against these fungal pathogens. The problems in the identification of resistance sources for B. juncea concerning genome complexity in host subpopulation and pathotypes were addressed. Emphasis has been laid on more elaborate and coordinated research to identify and deploy R genes, robust techniques, and research materials. Examples of fully characterized genes conferring resistance have been discussed that can be transformed into B. juncea using advanced genomics tools. Lastly, effective strategies for B. juncea improvement through introgression of novel R genes, development of pre-breeding resistant lines, characterization of pathotypes, and defense-related secondary metabolites have been provided suggesting the plan for the development of resistant B. juncea.

Tackling Control of a Cosmopolitan Phytopathogen: Sclerotinia

Phytopathogenic members of the Sclerotinia genus cause widespread disease across a broad range of economically important crops. In particular, Sclerotinia sclerotiorum is considered one of the most destructive and cosmopolitan of plant pathogens. Here, were review the epidemiology of the pathogen, its economic impact on agricultural production, and measures employed toward control of disease. We review the broad approaches required to tackle Sclerotinia diseases and include cultural practices, crop genetic resistance, chemical fungicides, and biological controls. We highlight the benefits and drawbacks of each approach along with recent advances within these controls and future strategies.

Genotype-Specific Antioxidant Responses and Assessment of Resistance Against Sclerotinia sclerotiorum Causing Sclerotinia Rot in Indian Mustard

Productivity of Indian mustard, an important oilseed crop of India, is affected by several pathogens. Among them, the hemibiotroph Sclerotinia sclerotiorum, which causes sclerotinia rot disease, is the most devastating fungal pathogen causing up to 90% yield losses. The availability of host resistance is the only efficient approach to control and understand the host-pathogen interaction. Therefore, the present investigation was carried out using six Indian mustard genotypes with contrasting behavior towards sclerotinia rot to study the antioxidant resistance mechanism against S. sclerotiorum. The plants at post-flowering stage were inoculated with five-day-old pure culture of S. sclerotiorum using artificial stem inoculation method. Disease evaluation revealed significant genotypic differences for mean lesion length among the tested genotypes, where genotype DRMR 2035 was found highly resistant, while genotypes RH 1569 and RH 1633 were found highly susceptible. The resistant genotypes had more phenolics and higher activities of peroxidase, catalase and polyphenol oxidase which provide them more efficient and strong antioxidant systems as compared with susceptible genotypes. Studies of antioxidative mechanisms validate the results of disease responses.

An update on the arsenal: mining resistance genes for disease management of Brassica crops in the genomic era

Brassica species include many economically important crops that provide nutrition and health-promoting substances to humans worldwide. However, as with all crops, their production is constantly threatened by emerging viral, bacterial, and fungal diseases, whose incidence has increased in recent years. Traditional methods of control are often costly, present limited effectiveness, and cause environmental damage; instead, the ideal approach is to mine and utilize the resistance genes of the Brassica crop hosts themselves. Fortunately, the development of genomics, molecular genetics, and biological techniques enables us to rapidly discover and apply resistance (R) genes. Herein, the R genes identified in Brassica crops are summarized, including their mapping and cloning, possible molecular mechanisms, and application in resistance breeding. Future perspectives concerning how to accurately discover additional R gene resources and efficiently utilize these genes in the genomic era are also discussed.

ddRAD sequencing-based identification of inter-genepool SNPs and association analysis in Brassica juncea

BACKGROUND

Narrow genetic base, complex allo-tetraploid genome and presence of repetitive elements have led the discovery of single nucleotide polymorphisms (SNPs) in Brassica juncea (AABB; 2n = 4x = 36) at a slower pace. Double digest RAD (ddRAD) - a genome complexity reduction technique followed by NGS was used to generate a total of 23 million paired-end reads from three genotypes each of Indian (Pusa Tarak, RSPR-01 and Urvashi) and Exotic (Donskaja IV, Zem 1 and EC287711) genepools.

RESULTS

Sequence data analysis led to the identification of 10,399 SNPs in six genotypes at a read depth of 10x coverage among the genotypes of two genepools. A total of 44 hyper-variable regions (nucleotide variation hotspots) were also found in the genome, of which 93% were found to be a part of coding genes/regions. The functionality of the identified SNPs was estimated by genotyping a subset of SNPs on MassARRAY® platform among a diverse set of B. juncea genotypes. SNP genotyping-based genetic diversity and population studies placed the genotypes into two distinct clusters based mostly on the place of origin. The genotypes were also characterized for six morphological traits, analysis of which revealed a significant difference in the mean values between Indian and Exotic genepools for six traits. The association analysis for six traits identified a total of 45 significant marker-trait associations on 11 chromosomes of A- and B- group of progenitor genomes.

CONCLUSIONS

Despite narrow diversity, the ddRAD sequencing was able to identify large number of nucleotide polymorphisms between the two genepools. Association analysis led to the identification of common SNPs/genomic regions associated between flowering and maturity traits, thereby underscoring the possible role of common chromosomal regions-harboring genes controlling flowering and maturity in Brassica juncea.

Unraveling candidate genes underlying biomass digestibility in elephant grass (Cenchrus purpureus)

BACKGROUND

Elephant grass [Cenchrus purpureus (Schumach.) Morrone] is used for bioenergy and animal feed. In order to identify candidate genes that could be exploited for marker-assisted selection in elephant grass, this study aimed to investigate changes in predictive accuracy using genomic relationship information and simple sequence repeats for eight traits (height, green biomass, dry biomass, acid and neutral detergent fiber, lignin content, biomass digestibility, and dry matter concentration) linked to bioenergetics and animal feeding.

RESULTS

We used single-step, genome-based best linear unbiased prediction and genome association methods to investigate changes in predictive accuracy and find candidate genes using genomic relationship information. Genetic variability (p < 0.05) was detected for most of the traits evaluated. In general, the overall means for the traits varied widely over the cuttings, which was corroborated by a significant genotype by cutting interaction. Knowing the genomic relationships increased the predictive accuracy of the biomass quality traits. We found that one marker (M28_161) was significantly associated with high values of biomass digestibility. The marker had moderate linkage disequilibrium with another marker (M35_202) that, in general, was detected in genotypes with low values of biomass digestibility. In silico analysis revealed that both markers have orthologous regions in other C4 grasses such as Setaria viridis, Panicum hallii, and Panicum virgatum, and these regions are located close to candidate genes involved in the biosynthesis of cell wall molecules (xyloglucan and lignin), which support their association with biomass digestibility.

CONCLUSIONS

The markers and candidate genes identified here are useful for breeding programs aimed at changing biomass digestibility in elephant grass. These markers can be used in marker-assisted selection to grow elephant grass cultivars for different uses, e.g., bioenergy production, bio-based products, co-products, bioactive compounds, and animal feed.

Genome-wide association study reveals genomic regions controlling root and shoot traits at late growth stages in wheat

BACKGROUND AND AIMS

Root system morphology is important for sustainable agriculture, but the genetic basis of root traits and their relationship to shoot traits remain to be elucidated. The aim of the present study was to dissect the genetic basis of root traits at late growth stages and its implications on shoot traits in wheat.

METHODS

Among 323 wheat accessions, we investigated phenotypic differences in root traits at booting and mid-grain fill stages in PVC tubes, shoot traits including plant height (PH), canopy temperature (CT) and grain yield per plant (YPP) in a field experiment, and performed a genome-wide association study with a Wheat 660K SNP Array.

KEY RESULTS

Deep-rooted accessions had lower CT and higher YPP than those with shallow roots, but no significant relationship was identified between root dry weight and shoot traits. Ninety-three significantly associated loci (SALs) were detected by the mixed linear model, among which three were hub SALs (Co-6A, Co-6B and Co-6D) associated with root depth at both booting and mid-grain fill stages, as well as CT and YPP. Minirhizotron system scanning results suggested that the causal genes in the three SALs may regulate root elongation in the field. The heritable independence between root depth and PH was demonstrated by linkage disequilibrium analysis. The YPP was significantly higher in genotypes which combined favourable marker alleles (FMAs) for root depth and PH, suggesting that a deep root and shorter plant height are suitable traits for pyramiding target alleles by molecular marker-assisted breeding.

CONCLUSIONS

These results uncovered promising genomic regions for functional gene discovery of root traits in the late growth period, enhanced understanding of correlation between root and shoot traits, and will facilitate intensive study on root morphology and breeding through molecular design.

Detection of First Marker Trait Associations for Resistance Against Sclerotinia sclerotiorum in Brassica juncea-Erucastrum cardaminoides Introgression Lines

A set of 96 Brassica juncea-Erucastrum cardaminoides introgression lines (ILs) were developed with genomic regions associated with Sclerotinia stem rot (Sclerotinia sclerotiorum) resistance from a wild Brassicaceous species E. cardaminoides. ILs were assessed for their resistance responses to stem inoculation with S. sclerotiorum, over three crop seasons (season I, 2011/2012; II, 2014/2015; III, 2016-2017). Initially, ILs were genotyped with transferable SSR markers and subsequently through genotyping by sequencing. SSR based association mapping identified six marker loci associated to resistance in both A and B genomes. Subsequent genome-wide association analysis (GWAS) of 84 ILs recognized a large number of SNPs associated to resistance, in chromosomes A03, A06, and B03. Chromosomes A03 and A06 harbored the maximum number of resistance related SNPs. Annotation of linked genomic regions highlighted an array of resistance mechanisms in terms of signal transduction pathways, hypersensitive responses and production of anti-fungal proteins and metabolites. Of major importance was the clustering of SNPs, encoding multiple resistance genes on small regions spanning approximately 885 kb region on chromosome A03 and 74 kb on B03. Five SNPs on chromosome A03 (6,390,210-381) were associated with LRR-RLK (receptor like kinases) genes that encode LRR-protein kinase family proteins. Genetic factors associated with pathogen-associated molecular patterns (PAMPs) and effector-triggered immunity (ETI) were predicted on chromosome A03, exhibiting 11 SNPs (6,274,763-994). These belonged to three R-Genes encoding TIR-NBS-LRR proteins. Marker trait associations (MTAs) identified will facilitate marker assisted introgression of these critical resistances, into new cultivars of B. juncea initially and, subsequently, into other crop Brassica species.

NBS-Encoding Genes in Brassica napus Evolved Rapidly After Allopolyploidization and Co-localize With Known Disease Resistance Loci

Genes containing nucleotide-binding sites (NBS) play an important role in pathogen resistance in plants. However, the evolutionary fate of NBS-encoding genes after formation of allotetraploid Brassica napus (A_nA_nC_nC_n, 2n = 38) is still unknown. We performed a genome-wide comparison of putatively functional NBS-encoding genes in B. napus and its progenitor species Brassica rapa (A_rA_r, 2n = 20) and Brassica oleracea (C_oC_o, 2n = 18), identifying 464, 202, and 146 putatively functional NBS-encoding genes respectively, with genes unevenly distributed in several clusters. The A_n-subgenome of B. napus possessed similar numbers of NBS-encoding genes (191 genes) to the A_r genome of B. rapa (202 genes) and similar clustering patterns. However, the C_n genome of B. napus had many more genes (273) than the B. oleracea C_o genome (146), with different clustering trends. Only 97 NBS-encoding genes (66.4%) in B. oleracea were homologous with NBS-encoding genes in B. napus, while 176 NBS-encoding genes (87.1%) were homologous between B. rapa and B. napus. These results suggest a greater diversification of NBS-encoding genes in the C genome may have occurred after formation of B. napus. Although most NBS-encoding genes in B. napus appeared to derive from the progenitors, the birth and death of several NBS-encoding genes was also putatively mediated by non-homologous recombination. The Ka/Ks values of most homologous pairs between B. napus and the progenitor species were less than 1, suggesting purifying selection during B. napus evolution. The majority of NBS-encoding genes (60% in all species) showed higher expression levels in root tissue (out of root, leaf, stem, seed and flower tissue types). Comparative analysis of NBS-encoding genes with mapped resistance QTL against three major diseases of B. napus (blackleg, clubroot and Sclerotinia stem rot) found 204 NBS-encoding genes in B. napus located within 71 resistance QTL intervals. The majority of NBS-encoding genes were co-located with resistance QTLs against a single disease, while 47 genes were co-located with QTLs against two diseases and 3 genes were co-located with QTLs against all three. Our results revealed significant variation as well as interesting evolutionary trajectories of NBS-encoding genes in the different Brassica subgenomes, while co-localization of NBS-encoding genes and resistance QTL may facilitate resistance breeding in oilseed rape.

Integrated mRNA, sRNA, and degradome sequencing reveal oilseed rape complex responses to Sclerotinia sclerotiorum (Lib.) infection

Sclerotinia stem rot (SSR), caused by the fungal pathogen Sclerotinia sclerotiorum, is a devastating disease resulting in yield losses and decreases in seed quality in oilseed rape (Brassica napus) worldwide. However, the molecular mechanisms underlying the response of oilseed rape to S. sclerotiorum infection at the transcriptional and post-transcriptional levels are poorly understood. Here, we used an integrated omics approach (transcriptome, sRNAome, and degradome sequencing) on the Illumina platform to compare the RNA expression and post-transcriptional profiles of oilseed rape plants inoculated or not with S. sclerotiorum. In total, 7,065 differentially expressed genes (DEGs) compared with the mock-inoculated control at 48 hours post inoculation were identified. These DEGs were associated with protein kinases, signal transduction, transcription factors, hormones, pathogenesis-related proteins, secondary metabolism, and transport. In the sRNA-Seq analysis, 77 known and 176 novel miRNAs were identified; however, only 10 known and 41 novel miRNAs were differentially expressed between the samples inoculated or not with S. sclerotiorum. Degradome sequencing predicted 80 cleavage sites with 64 miRNAs. Integrated mRNA, sRNA and degradome sequencing analysis reveal oilseed rape complex responses to S. sclerotiorum infection. This study provides a global view of miRNA and mRNA expression profiles in oilseed rape following S. sclerotiorum infection.

Association mapping unveils favorable alleles for grain iron and zinc concentrations in lentil (Lens culinaris subsp. culinaris)

Lentil is a major cool-season grain legume grown in South Asia, West Asia, and North Africa. Populations in developing countries of these regions have micronutrient deficiencies; therefore, breeding programs should focus more on improving the micronutrient content of food. In the present study, a set of 96 diverse germplasm lines were evaluated at three different locations in India to examine the variation in iron (Fe) and zinc (Zn) concentration and identify simple sequence repeat (SSR) markers that associate with the genetic variation. The genetic variation among genotypes of the association mapping (AM) panel was characterized using a genetic distance-based and a general model-based clustering method. The model-based analysis identified six subpopulations, which satisfactorily explained the genetic structure of the AM panel. AM analysis identified three SSRs (PBALC 13, PBALC 206, and GLLC 563) associated with grain Fe concentration explaining 9% to 11% of phenotypic variation and four SSRs (PBALC 353, SSR 317-1, PLC 62, and PBALC 217) were associated with grain Zn concentration explaining 14%, to 21% of phenotypic variation. These identified SSRs exhibited consistent performance across locations. These candidate SSRs can be used in marker-assisted genetic improvement for developing Fe and Zn fortified lentil varieties. Favorable alleles and promising genotypes identified in this study can be utilized for lentil biofortification.

Current Status and Challenges in Identifying Disease Resistance Genes in Brassica napus

Brassica napus is an economically important crop across different continents including temperate and subtropical regions in Europe, Canada, South Asia, China and Australia. Its widespread cultivation also brings setbacks as it plays host to fungal, oomycete and chytrid pathogens that can lead to serious yield loss. For sustainable crop production, identification of resistance (R) genes in B. napus has become of critical importance. In this review, we discuss four key pathogens affecting Brassica crops: Clubroot (Plasmodiophora brassicae), Blackleg (Leptosphaeria maculans and L. biglobosa), Sclerotinia Stem Rot (Sclerotinia sclerotiorum), and Downy Mildew (Hyaloperonospora parasitica). We first review current studies covering prevalence of these pathogens on Brassica crops and highlight the R genes and QTL that have been identified from Brassica species against these pathogens. Insights into the relationships between the pathogen and its Brassica host, the unique host resistance mechanisms and how these affect resistance outcomes is also presented. We discuss challenges in identification and deployment of R genes in B. napus in relation to highly specific genetic interactions between host subpopulations and pathogen pathotypes and emphasize the need for common or shared techniques and research materials or tighter collaboration between researchers to reconcile the inconsistencies in the research outcomes. Using current genomics tools, we provide examples of how characterization and cloning of R genes in B. napus can be carried out more effectively. Lastly, we put forward strategies to breed resistant cultivars through introgressions supported by genomic approaches and suggest prospects that can be implemented in the future for a better, pathogen-resistant B. napus.

Genetic and transcriptomic analyses of lignin- and lodging-related traits in Brassica napus

Candidate genes associated with lignin and lodging traits were identified by combining phenotypic, genotypic, and gene expression data in B. napus. Brassica napus is one of the world's most important oilseed crops, but its yield can be dramatically reduced by lodging, bending, and falling of its vertical stems. Lignin has been shown to contribute to stem mechanical strength. In this study, we found that the syringyl/guaiacyl (S/G) monolignol ratio exhibits a significant negative correlation with disease and lodging resistance. A total of 92 and 50 SNP and SSR loci, respectively, were found to be significantly associated with five traits, breaking force, breaking strength, lodging coefficient, acid detergent lignin content, and the S/G monolignol ratio using GWAS. To identify novel genes involved in lignin biosynthesis, transcriptome sequencing of high- (H) and low (L)-ADL content accessions was performed. The up-regulated genes were mainly involved in glycoside catabolic processes (especially glucosinolate catabolism) and cell wall biogenesis, while down-regulated genes were involved in glucosinolate biosynthesis, indicating that crosstalk exists between glucosinolate metabolic processes and lignin biosynthesis. Integrating this differential expression with the GWAS analysis, we identified four candidate genes regulating lignin, including glycosyl hydrolase (BnaA01g00480D), CYT1 (BnaA04g22820D), and two encoding transcription factors, SHINE1 (ERF family) and DAR6 (LIM family). This study provides insight into the genetic control of lodging and lignin in B. napus.