Genome-wide association analysis and differential expression analysis of resistance to Sclerotinia stem rot in Brassica napus

Natural variation in BnaA07.MKK9 confers resistance to Sclerotinia stem rot in oilseed rape

Sclerotinia stem rot (SSR), caused by the necrotrophic fungus Sclerotinia sclerotiorum, is one of the most devastating diseases for several major oil-producing crops. Despite its impact, the genetic basis of SSR resistance in plants remains poorly understood. Here, through a genome-wide association study, we identify a key gene, BnaA07. MKK9, that encodes a mitogen-activated protein kinase kinase that confers SSR resistance in oilseed rape. Our functional analyses reveal that BnaA07.MKK9 interacts with BnaC03.MPK3 and BnaC03.MPK6 and phosphorylates them at the TEY activation motif, triggering a signaling cascade that initiates biosynthesis of ethylene, camalexin, and indole glucosinolates, and promotes accumulation of H2O2 and the hypersensitive response, ultimately conferring resistance. Furthermore, variations in the coding sequence of BnaA07.MKK9 alter its kinase activity and improve SSR resistance by ~30% in cultivars carrying the advantageous haplotype. These findings enhance our understanding of SSR resistance and may help engineer novel diversity for future breeding of oilseed rape.

Identification and characterization of the plasma membrane H+-ATPase genes in Brassica napus and functional analysis of BnHA9 in salt tolerance

As a primary proton pump, plasma membrane (PM) H+-ATPase plays critical roles in regulating plant growth, development, and stress responses. PM H+-ATPases have been well characterized in many plant species. However, no comprehensive study of PM H+-ATPase genes has been performed in Brassica napus (rapeseed). In this study, we identified 32 PM H+-ATPase genes (BnHAs) in the rapeseed genome, and they were distributed on 16 chromosomes. Phylogenetical and gene duplication analyses showed that the BnHA genes were classified into five subfamilies, and the segmental duplication mainly contributed to the expansion of the rapeseed PM H+-ATPase gene family. The conserved domain and subcellular analyses indicated that BnHAs encoded canonical PM H+-ATPase proteins with 14 highly conserved domains and localized on PM. Cis-acting regulatory element and expression pattern analyses indicated that the expression of BnHAs possessed tissue developmental stage specificity. The 25 upstream open reading frames with the canonical initiation codon ATG were predicted in the 5' untranslated regions of 11 BnHA genes and could be used as potential target sites for improving rapeseed traits. Protein interaction analysis showed that BnBRI1.c associated with BnHA2 and BnHA17, indicating that the conserved activity regulation mechanism of BnHAs may be present in rapeseed. BnHA9 overexpression in Arabidopsis enhanced the salt tolerance of the transgenic plants. Thus, our results lay a foundation for further research exploring the biological functions of PM H+-ATPases in rapeseed.

Phenylpropanoids Following Wounding and Infection of Sweet Sorghum Lines Differing in Responses to Stalk Pathogens

Sweet sorghum (Sorghum bicolor) lines M81-E and Colman were previously shown to differ in responses to Fusarium thapsinum and Macrophomina phaseolina, stalk rot pathogens that can reduce the yields and quality of biomass and extracted sugars. Inoculated tissues were compared for transcriptomic, phenolic metabolite, and enzymatic activity during disease development 3 and 13 days after inoculation (DAI). At 13 DAI, M81-E had shorter mean lesion lengths than Colman when inoculated with either pathogen. Transcripts encoding monolignol biosynthetic and modification enzymes were associated with transcriptional wound (control) responses of both lines at 3 DAI. Monolignol biosynthetic genes were differentially coexpressed with transcriptional activator SbMyb76 in all Colman inoculations, but only following M. phaseolina inoculation in M81-E, suggesting that SbMyb76 is associated with lignin biosynthesis during pathogen responses. In control inoculations, defense-related genes were expressed at higher levels in M81-E than Colman. Line, treatment, and timepoint differences observed in phenolic metabolite and enzyme activities did not account for observed differences in lesions. However, generalized additive models were able to relate metabolites, but not enzyme activities, to lesion length for quantitatively modeling disease progression: in M81-E, but not Colman, sinapic acid levels positively predicted lesion length at 3 DAI when cell wall-bound syringic acid was low, soluble caffeic acid was high, and lactic acid was high, suggesting that sinapic acid may contribute to responses at 3 DAI. These results provide potential gene targets for development of sweet sorghum varieties with increased stalk rot resistance to ensure biomass and sugar quality.

Genome-wide association study of powdery mildew resistance in cultivated soybean from Northeast China

Powdery mildew (PMD), caused by the pathogen Microsphaera diffusa, leads to substantial yield decreases in susceptible soybean under favorable environmental conditions. Effective prevention of soybean PMD damage can be achieved by identifying resistance genes and developing resistant cultivars. In this study, we genotyped 331 soybean germplasm accessions, primarily from Northeast China, using the SoySNP50K BeadChip, and evaluated their resistance to PMD in a greenhouse setting. To identify marker-trait associations while effectively controlling for population structure, we conducted genome-wide association studies utilizing factored spectrally transformed linear mixed models, mixed linear models, efficient mixed-model association eXpedited, and compressed mixed linear models. The results revealed seven single nucleotide polymorphism (SNP) loci strongly associated with PMD resistance in soybean. Among these, one SNP was localized on chromosome (Chr) 14, and six SNPs with low linkage disequilibrium were localized near or in the region of previously mapped genes on Chr 16. In the reference genome of Williams82, we discovered 96 genes within the candidate region, including 17 resistance (R)-like genes, which were identified as potential candidate genes for PMD resistance. In addition, we performed quantitative real-time reverse transcription polymerase chain reaction analysis to evaluate the gene expression levels in highly resistant and susceptible genotypes, focusing on leaf tissues collected at different times after M. diffusa inoculation. Among the examined genes, three R-like genes, including Glyma.16G210800, Glyma.16G212300, and Glyma.16G213900, were identified as strong candidates associated with PMD resistance. This discovery can significantly enhance our understanding of soybean resistance to PMD. Furthermore, the significant SNPs strongly associated with resistance can serve as valuable markers for genetic improvement in breeding M. diffusa-resistant soybean cultivars.

Screening of microRNAs and target genes involved in Sclerotinia sclerotiorum (Lib.) infection in Brassica napus L

BACKGROUND

Rapeseed (Brassica napus L.) is the third largest source of vegetable oil in the world, and Sclerotinia sclerotiorum (Lib.) is a major soil-borne fungal plant pathogen that infects more than 400 plant species, including B. napus. Sclerotinia stem rot caused an annual loss of 10 - 20% in rapeseed yield. Exploring the molecular mechanisms in response to S. sclerotiorum infection in B. napus is beneficial for breeding and cultivation of resistant varieties. To gain a better understanding of the mechanisms regarding B. napus tolerance to Sclerotinia stem rot, we employed a miRNAome sequencing approach and comprehensively investigated global miRNA expression profile among five relatively resistant lines and five susceptible lines of oilseed at 0, 24, and 48 h post-inoculation.

RESULTS

In this study, a total of 40 known and 1105 novel miRNAs were differentially expressed after S. sclerotiorum infection, including miR156, miR6028, miR394, miR390, miR395, miR166, miR171, miR167, miR164, and miR172. Furthermore, 8,523 genes were predicted as targets for these differentially expressed miRNAs. These target genes were mainly associated with disease resistance (R) genes, signal transduction, transcription factors, and hormones. Constitutively expressing miR156b (OX156b) plants strengthened Arabidopsis resistance against S. sclerotiorum accompanied by smaller necrotic lesions, whereas blocking miR156 expression in Arabidopsis (MIM156) led to greater susceptibility to S. sclerotiorum disease, associated with extensive cell death of necrotic lesions.

CONCLUSIONS

This study reveals the distinct difference in miRNA profiling between the relatively resistant lines and susceptible lines of B. napus in response to S. sclerotiorum. The identified differentially expressed miRNAs related to sclerotinia stem rot resistance are involved in regulating resistance to S. sclerotiorum in rapeseed by targeting genes related to R genes, signal transduction, transcription factors, and hormones. miR156 positively modulates the resistance to S. sclerotiorum infection by restricting colonization of S. sclerotiorum mycelia. This study provides a broad view of miRNA expression changes after S. sclerotiorum infection in oilseed and is the first to elucidate the function and mechanism underlying the miR156 response to S. sclerotiorum infection in oilseed rape.

Integrating genome-wide association study with transcriptomic data to predict candidate genes influencing Brassica napus root and biomass-related traits under low phosphorus conditions

BACKGROUND

Rapeseed (Brassica napus L.) is an essential source of edible oil and livestock feed, as well as a promising source of biofuel. Breeding crops with an ideal root system architecture (RSA) for high phosphorus use efficiency (PUE) is an effective way to reduce the use of phosphate fertilizers. However, the genetic mechanisms that underpin PUE in rapeseed remain elusive. To address this, we conducted a genome-wide association study (GWAS) in 327 rapeseed accessions to elucidate the genetic variability of 13 root and biomass traits under low phosphorus (LP; 0.01 mM P +). Furthermore, RNA-sequencing was performed in root among high/low phosphorus efficient groups (HP1/LP1) and high/low phosphorus stress tolerance groups (HP2/LP2) at two-time points under control and P-stress conditions.

RESULTS

Significant variations were observed in all measured traits, with heritabilities ranging from 0.47 to 0.72, and significant correlations were found between most of the traits. There were 39 significant trait-SNP associations and 31 suggestive associations, which integrated into 11 valid quantitative trait loci (QTL) clusters, explaining 4.24-24.43% of the phenotypic variance observed. In total, RNA-seq identified 692, 1076, 648, and 934 differentially expressed genes (DEGs) specific to HP1/LP1 and HP2/LP2 under P-stress and control conditions, respectively, while 761 and 860 DEGs common for HP1/LP1 and HP2/LP2 under both conditions. An integrated approach of GWAS, weighted co-expression network, and differential expression analysis identified 12 genes associated with root growth and development under LP stress. In this study, six genes (BnaA04g23490D, BnaA09g08440D, BnaA09g04320D, BnaA09g04350D, BnaA09g04930D, BnaA09g09290D) that showed differential expression were identified as promising candidate genes for the target traits.

CONCLUSION

11 QTL clusters and 12 candidate genes associated with root and development under LP stress were identified in this study. Our study's phenotypic and genetic information may be exploited for genetic improvement of root traits to increase PUE in rapeseed.

LESION MIMIC MUTANT 1 confers basal resistance to Sclerotinia sclerotiorum in rapeseed via a salicylic acid-dependent pathway

Rapeseed (Brassica napus) is a major edible oilseed crop consumed worldwide. However, its yield is seriously affected by infection from the broad-spectrum non-obligate pathogen Sclerotinia sclerotiorum due to a lack of highly resistant germplasm. Here, we identified a Sclerotinia-resistant and light-dependent lesion mimic mutant from an ethyl methanesulfonate-mutagenized population of the rapeseed inbred Zhongshuang 11 (ZS11) named lesion mimic mutant 1 (lmm1). The phenotype of lmm1 is controlled by a single recessive gene, named LESION MIMIC MUTANT 1 (LMM1), which mapped onto chromosome C04 by bulked segregant analysis within a 2.71-Mb interval. Histochemical analysis indicated that H2O2 strongly accumulated and cell death occurred around the lesion mimic spots. Among 877 differentially expressed genes (DEGs) between ZS11 and lmm1 leaves, 188 DEGs were enriched in the defense response, including 95 DEGs involved in systemic acquired resistance, which is consistent with the higher salicylic acid levels in lmm1. Combining bulked segregant analysis and transcriptome analysis, we identified a significantly up-regulated gene, BnaC4.PR2, which encodes β-1,3-glucanase, as the candidate gene for LMM1. Overexpression of BnaC4.PR2 may induce a reactive oxygen species burst to trigger partial cell death and systemic acquired resistance. Our study provides a new genetic resource for S. sclerotiorum resistance as well as new insights into disease resistance breeding in B. napus.

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.

New insight into the genetic basis of oil content based on noninvasive three-dimensional phenotyping and tissue-specific transcriptome in Brassica napus

BACKGROUND

Increasing seed oil content is the most important breeding goal in Brassica napus, and phenotyping is crucial to dissect its genetic basis in crops. To date, QTL mapping for oil content has been based on whole seeds, and the lipid distribution is far from uniform in different tissues of seeds in B. napus. In this case, the phenotype based on whole seeds was unable to sufficiently reveal the complex genetic characteristics of seed oil content.

RESULTS

Here, the three-dimensional (3D) distribution of lipid was determined for B. napus seeds by magnetic resonance imaging (MRI) and 3D quantitative analysis, and ten novel oil content-related traits were obtained by subdividing the seeds. Based on a high-density genetic linkage map, 35 QTLs were identified for 4 tissues, the outer cotyledon (OC), inner cotyledon (IC), radicle (R) and seed coat (SC), which explained up to 13.76% of the phenotypic variation. Notably, 14 tissue-specific QTLs were reported for the first time, 7 of which were novel. Moreover, haplotype analysis showed that the favorable alleles for different seed tissues exhibited cumulative effects on oil content. Furthermore, tissue-specific transcriptomes revealed that more active energy and pyruvate metabolism influenced carbon flow in the IC, OC and R than in the SC at the early and middle seed development stages, thus affecting the distribution difference in oil content. Combining tissue-specific QTL mapping and transcriptomics, 86 important candidate genes associated with lipid metabolism were identified that underlie 19 unique QTLs, including the fatty acid synthesis rate-limiting enzyme-related gene CAC2, in the QTLs for OC and IC.

CONCLUSIONS

The present study provides further insight into the genetic basis of seed oil content at the tissue-specific level.

Functional and evolutionary study of MLO gene family in the regulation of Sclerotinia stem rot resistance in Brassica napus L

BACKGROUND

Oilseed rape (Brassica napus L.) is known as one of the most important oilseed crops cultivated around the world. However, its production continuously faces a huge challenge of Sclerotinia stem rot (SSR), a destructive disease caused by the fungus Sclerotinia sclerotiorum, resulting in huge yield loss annually. The SSR resistance in B. napus is quantitative and controlled by a set of minor genes. Identification of these genes and pyramiding them into a variety are a major strategy for SSR resistance breeding in B. napus.

RESULTS

Here, we performed a genome-wide association study (GWAS) using a natural population of B. napus consisting of 222 accessions to identify BnaA08g25340D (BnMLO2_2) as a candidate gene that regulates the SSR resistance. BnMLO2_2 was a member of seven homolog genes of Arabidopsis Mildew Locus O 2 (MLO2) and the significantly SNPs were mainly distributed in the promoter of BnMLO2_2, suggesting a role of BnMLO2_2 expression level in the regulation of SSR resistance. We expressed BnMLO2_2 in Arabidopsis and the transgenic plants displayed an enhanced SSR resistance. Transcriptome profiling of different tissues of B. napus revealed that BnMLO2_2 had the most expression level in leaf and silique tissues among all the 7 BnMLO2 members and also expressed higher in the SSR resistant accession than in the susceptible accession. In Arabidopsis, mlo2 plants displayed reduced resistance to SSR, whereas overexpression of MLO2 conferred plants an enhanced SSR resistance. Moreover, a higher expression level of MLO2 showed a stronger SSR resistance in the transgenic plants. The regulation of MLO2 in SSR resistance may be associated with the cell death. Collinearity and phylogenetic analysis revealed a large expansion of MLO family in Brassica crops.

CONCLUSION

Our study revealed an important role of BnMLO2 in the regulation of SSR resistance and provided a new gene candidate for future improvement of SSR resistance in B. napus and also new insights into understanding of MLO family evolution in Brassica crops.

Jasmonic Acid as a Mediator in Plant Response to Necrotrophic Fungi

Jasmonic acid (JA) and its derivatives, all named jasmonates, are the simplest phytohormones which regulate multifarious plant physiological processes including development, growth and defense responses to various abiotic and biotic stress factors. Moreover, jasmonate plays an important mediator’s role during plant interactions with necrotrophic oomycetes and fungi. Over the last 20 years of research on physiology and genetics of plant JA-dependent responses to pathogens and herbivorous insects, beginning from the discovery of the JA co-receptor CORONATINE INSENSITIVE1 (COI1), research has speeded up in gathering new knowledge on the complexity of plant innate immunity signaling. It has been observed that biosynthesis and accumulation of jasmonates are induced specifically in plants resistant to necrotrophic fungi (and also hemibiotrophs) such as mostly investigated model ones, i.e., Botrytis cinerea, Alternaria brassicicola or Sclerotinia sclerotiorum. However, it has to be emphasized that the activation of JA-dependent responses takes place also during susceptible interactions of plants with necrotrophic fungi. Nevertheless, many steps of JA function and signaling in plant resistance and susceptibility to necrotrophs still remain obscure. The purpose of this review is to highlight and summarize the main findings on selected steps of JA biosynthesis, perception and regulation in the context of plant defense responses to necrotrophic 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.

Comparative transcriptome profiling reveals the importance of GmSWEET15 in soybean susceptibility to Sclerotinia sclerotiorum

Soybean sclerotinia stem rot (SSR) is a disease caused by Sclerotinia sclerotiorum that causes incalculable losses in soybean yield each year. Considering the lack of effective resistance resources and the elusive resistance mechanisms, we are urged to develop resistance genes and explore their molecular mechanisms. Here, we found that loss of GmSWEET15 enhanced the resistance to S. sclerotiorum, and we explored the molecular mechanisms by which gmsweet15 mutant exhibit enhanced resistance to S. sclerotiorum by comparing transcriptome. At the early stage of inoculation, the wild type (WT) showed moderate defense response, whereas gmsweet15 mutant exhibited more extensive and intense transcription reprogramming. The gmsweet15 mutant enriched more biological processes, including the secretory pathway and tetrapyrrole metabolism, and it showed stronger changes in defense response, protein ubiquitination, MAPK signaling pathway-plant, plant-pathogen interaction, phenylpropanoid biosynthesis, and photosynthesis. The more intense and abundant transcriptional reprogramming of gmsweet15 mutant may explain how it effectively delayed colonization by S. sclerotiorum. In addition, we identified common and specific differentially expressed genes between WT and gmsweet15 mutant after inoculation with S. sclerotiorum, and gene sets and genes related to gmsweet15_24 h were identified through Gene Set Enrichment Analysis. Moreover, we constructed the protein-protein interaction network and gene co-expression networks and identified several groups of regulatory networks of gmsweet15 mutant in response to S. sclerotiorum, which will be helpful for the discovery of candidate functional genes. Taken together, our results elucidate molecular mechanisms of delayed colonization by S. sclerotiorum after loss of GmSWEET15 in soybean, and we propose novel resources for improving resistance to SSR.

Genetic dissection of fatty acid components in the Chinese peanut (Arachis hypogaea L.) mini-core collection under multi-environments

Peanut (Arachis hypogaea L.) is an important source of edible oil and protein for human nutrition. The quality of peanut seed oil is mainly determined by the composition of fatty acids, especially the contents of oleic acid and linoleic acid. Improving the composition of fatty acids in the seed oil is one of the main objectives for peanut breeding globally. To uncover the genetic basis of fatty acids and broaden the genetic variation in future peanut breeding programs, this study used genome-wide association studies (GWAS) to identify loci associated with target traits and developed diagnostic marker. The contents of eight fatty acid components of the Chinese peanut mini-core collection were measured under four environments. Using the phenotypic information and over one hundred thousand single nucleotide polymorphisms (SNPs), GWAS were conducted to investigate the genetics basis of fatty acids under multi-environments. Overall, 75 SNPs were identified significant trait associations with fatty acid components. Nineteen associations were repeatedly identified in multiple environments, and 13 loci were co-associated with two or three traits. Three stable major associated loci were identified, including two loci for oleic acid and linoleic acid on chromosome A09 [mean phenotypic variation explained (PVE): 38.5%, 10.35%] and one for stearic acid on B06 (mean PVE: 23%). According to functional annotations, 21 putative candidate genes related to fatty acid biosynthesis were found underlying the three associations. The allelic effect of SNP A09-114690064 showed that the base variation was highly correlated with the phenotypic variation of oleic acid and linoleic acid contents, and a cost-effective Kompetitive allele-Specific PCR (KASP) diagnostic marker was developed. Furthermore, the SNP A09-114690064 was found to change the cis-element CAAT (-) in the promoter of ahFAD2A to YACT (+), leading dozens of times higher expression level. The enhancer-like activity of ahFAD2A promoter was identified that was valuable for enriching the regulation mechanism of ahFAD2A. This study improved our understanding on the genetic architecture of fatty acid components in peanut, and the new effective diagnostic marker would be useful for marker-assisted selection of high-oleic peanut breeding.

Transcription factor WRKY28 curbs WRKY33-mediated resistance to Sclerotinia sclerotiorum in Brassica napus

Sclerotinia sclerotiorum causes substantial damage and loss of yield in oilseed rape (Brassica napus). The molecular mechanisms of oilseed rape defense against Sclerotinia remain elusive. In this study, we found that in the early stages of B. napus infection a conserved mitogen-activated protein kinase (MAPK) cascade mediated by BnaA03.MKK5-BnaA06.MPK3/BnaC03.MPK3 module phosphorylates the substrate BnWRKY33, enhancing its transcriptional activity. The activated BnWRKY33 binds to its own promoter and triggers a transcriptional burst of BnWRKY33, thus helping plants effectively resist the pathogenic fungi by enhancing the expression of phytoalexin synthesis-related genes. The expression of BnWRKY33 is fine-tuned during defense. Ongoing Sclerotinia infection induces BnaA03.WRKY28 and BnaA09.VQ12 expression. BnaA09.VQ12 interacts physically with BnaA03.WRKY28 to form a protein complex, causing BnaA03.WRKY28 to outcompete BnWRKY33 and bind to the BnWRKY33 promoter. BnaA03.WRKY28 induction suppresses BnWRKY33 expression in the later stages of infection but promotes branch formation in the leaf axils by regulating the expression of branching-related genes such as BnBRC1. BnaA03.WRKY28 participates in the trade-off between defense and growth. These findings suggest that oilseed rape plants may modulate defense-response strength and develop alternative reproduction and survival strategies in the face of lethal pathogens.

Quantitative trait locus mapping and improved resistance to sclerotinia stem rot in a backbone parent of rapeseed (Brassica napus L.)

There are three main challenges to improving sclerotinia stem rot (SSR) resistance in rapeseed (Brassica napus L.). First, breeding materials such as the backbone parents have not been extensively investigated, making the findings of previous studies difficult to directly implement. Second, SSR resistance and flowering time (FT) loci are typically linked; thus, use of these loci requires sacrifice of the rapeseed growth period. Third, the SSR resistance loci in susceptible materials are often neglected, thereby reducing the richness of resistant resources. This study was conducted to investigate the stem resistance, disease index, and FT of a doubled haploid population consisting of 151 lines constructed from the backbone parent 19514A and conventional rapeseed cultivar ZY50 within multiple environments. Quantitative trait locus (QTL) mapping revealed 13 stem resistance QTLs, 9 disease index QTLs, and 20 FT QTLs. QTL meta-analysis showed that uqA04, uqC03.1, and uqC03.2 were repeatable SSR resistance QTLs derived from different parents but not affected by the FT. Based on these three QTLs, we proposed a strategy for improving the SSR resistance of 19514A and ZY50. This study improves the understanding of the resistance to rapeseed SSR and genetic basis of FT and demonstrates that SSR resistance QTLs can be mined from parents with a minimal resistance level difference, thereby supporting the application of backbone parents in related research and resistance improvement.

Further insight into decreases in seed glucosinolate content based on QTL mapping and RNA-seq in Brassica napus L

The QTL hotspots determining seed glucosinolate content instead of only four HAG1 loci and elucidation of a potential regulatory model for rapeseed SGC variation. Glucosinolates (GSLs) are amino acid-derived, sulfur-rich secondary metabolites that function as biopesticides and flavor compounds, but the high seed glucosinolate content (SGC) reduces seed quality for rapeseed meal. To dissect the genetic mechanism and further reduce SGC in rapeseed, QTL mapping was performed using an updated high-density genetic map based on a doubled haploid (DH) population derived from two parents that showed significant differences in SGC. In 15 environments, a total of 162 significant QTLs were identified for SGC and then integrated into 59 consensus QTLs, of which 32 were novel QTLs. Four QTL hotspot regions (QTL-HRs) for SGC variation were discovered on chromosomes A09, C02, C07 and C09, including seven major QTLs that have previously been reported and four novel major QTLs in addition to HAG1 loci. SGC was largely determined by superimposition of advantage allele in the four QTL-HRs. Important candidate genes directly related to GSL pathways were identified underlying the four QTL-HRs, including BnaC09.MYB28, BnaA09.APK1, BnaC09.SUR1 and BnaC02.GTR2a. Related differentially expressed candidates identified in the minor but environment stable QTLs indicated that sulfur assimilation plays an important rather than dominant role in SGC variation. A potential regulatory model for rapeseed SGC variation constructed by combining candidate GSL gene identification and differentially expressed gene analysis based on RNA-seq contributed to a better understanding of the GSL accumulation mechanism. This study provides insights to further understand the genetic regulatory mechanism of GSLs, as well as the potential loci and a new route to further diminish the SGC in rapeseed.

Comparative transcriptome analysis of compatible and incompatible Brassica napus-Xanthomonas campestris interactions

Black rot caused by the vascular pathogenic bacterium Xanthomonas campestris pv. campestris (Xcc) is widespread in Brassicaceae plants and an infectious disease that causes large yield losses in oil seed rape (Brassica napus L.). Improvement of resistance through breeding is a crucial strategy to prevent black rot disease in B. napus, but presently hampered by insufficient understanding of Xcc-Brassica interactions. This study compares two EMS-mutagenized B. napus lines that show contrasting resistance levels to their susceptible progenitor. Patterns of differential gene expression between these B. napus lines were evaluated at three time points post inoculation by comparative RNA-seq analysis. In line with the observed disease phenotypes, the susceptible line ZS9mXccS-1 displayed a steady amount of differentially expressed genes (DEGs) at different time points of infection, whereas the resistant line ZS9mXccR-1 displayed a gradual increase in DEGs throughout the course of infection. Weighted gene co-expression network analysis (WGCNA) pinpointed multiple defense-related hub genes with potential central roles in immunity, including the cell surface receptor genes CRK11 and BIR1, and the associated downstream regulatory genes WRKY11 and PBL30. KEGG analysis of DEGs belonging to two distinct co-expression modules revealed enriched pathways associated with defense, including Ca²⁺-signaling, receptor-mediated immunity, and phytohormone balance. Taken together, our comparative transcriptome analysis provides new avenues to unravel the mechanisms underlying black rot resistance in B. napus.

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.

Analysis of Tissue-Specific Defense Responses to Sclerotinia sclerotiorum in Brassica napus

Sclerotinia stem rot (SSR) caused by Sclerotinia sclerotiorum (S. sclerotiorum) is the main disease threat of oilseed rape (Brassica napus), resulting in huge economic losses every year. SSR resistance manifests as quantitative disease resistance (QDR), and no gene with complete SSR resistance has been cloned or reported so far. Transcriptome analysis has revealed a large number of defense-related genes and response processes. However, the similarities and differences in the defense responses of different tissues are rarely reported. In this study, we analyzed the similarities and differences of different tissues in response to S. sclerotiorum at 24 h post inoculation (hpi) by using the published transcriptome data for respective leaf and stem inoculation. At 24 hpi, large differences in gene expression exist in leaf and stem, and there are more differentially expressed genes and larger expression differences in leaf. The leaf is more sensitive to S. sclerotiorum and shows a stronger response than stem. Different defense responses appear in the leaf and stem, and the biosynthesis of lignin, callose, lectin, chitinase, PGIP, and PR protein is activated in leaf. In the stem, lipid metabolism-mediated defense responses are obviously enhanced. For the common defense responses in both leaf and stem, the chain reactions resulting from signal transduction and biological process take the primary responsibility. This research will be beneficial to exploit the potential of different tissues in plant defense and find higher resistance levels of genotypic variability in different environments. Our results are significant in the identification of resistance genes and analysis of defense mechanisms.

Genome-Wide Identification and Expression Analysis of nsLTP Gene Family in Rapeseed (Brassica napus) Reveals Their Critical Roles in Biotic and Abiotic Stress Responses

Non-specific lipid transfer proteins (nsLTPs) are small cysteine-rich basic proteins which play essential roles in plant growth, development and abiotic/biotic stress response. However, there is limited information about the nsLTP gene (BnLTP) family in rapeseed (Brassica napus). In this study, 283 BnLTP genes were identified in rapeseed, which were distributed randomly in 19 chromosomes of rapeseed. Phylogenetic analysis showed that BnLTP proteins were divided into seven groups. Exon/intron structure and MEME motifs both remained highly conserved in each BnLTP group. Segmental duplication and hybridization of rapeseed’s two sub-genomes mainly contributed to the expansion of the BnLTP gene family. Various potential cis-elements that respond to plant growth, development, biotic/abiotic stresses, and phytohormone signals existed in BnLTP gene promoters. Transcriptome analysis showed that BnLTP genes were expressed in various tissues/organs with different levels and were also involved in the response to heat, drought, NaCl, cold, IAA and ABA stresses, as well as the treatment of fungal pathogens (Sclerotinia sclerotiorum and Leptosphaeria maculans). The qRT-PCR assay validated the results of RNA-seq expression analysis of two top Sclerotinia-responsive BnLTP genes, BnLTP129 and BnLTP161. Moreover, batches of BnLTPs might be regulated by BnTT1 and BnbZIP67 to play roles in the development, metabolism or adaptability of the seed coat and embryo in rapeseed. This work provides an important basis for further functional study of the BnLTP genes in rapeseed quality improvement and stress resistance.

Deciphering the Genetic Basis of Root and Biomass Traits in Rapeseed (Brassica napus L.) through the Integration of GWAS and RNA-Seq under Nitrogen Stress

An excellent root system is responsible for crops with high nitrogen-use efficiency (NUE). The current study evaluated the natural variations in 13 root- and biomass-related traits under a low nitrogen (LN) treatment in a rapeseed association panel. The studied traits exhibited significant phenotypic differences with heritabilities ranging from 0.53 to 0.66, and most of the traits showed significant correlations with each other. The genome-wide association study (GWAS) found 51 significant and 30 suggestive trait–SNP associations that integrated into 14 valid quantitative trait loci (QTL) clusters and explained 5.7–21.2% phenotypic variance. In addition, RNA sequencing was performed at two time points to examine the differential expression of genes (DEGs) between high and low NUE lines. In total, 245, 540, and 399 DEGs were identified as LN stress-specific, high nitrogen (HN) condition-specific, and HNLN common DEGs, respectively. An integrated analysis of GWAS, weighted gene co-expression network, and DEGs revealed 16 genes involved in rapeseed root development under LN stress. Previous studies have reported that the homologs of seven out of sixteen potential genes control root growth and NUE. These findings revealed the genetic basis underlying nitrogen stress and provided worthwhile SNPs/genes information for the genetic improvement of NUE in rapeseed.

Genome-Wide Association Studies of Root-Related Traits in Brassica napus L. under Low-Potassium Conditions

Roots are essential organs for a plant’s ability to absorb water and obtain mineral nutrients, hence they are critical to its development. Plants use root architectural alterations to improve their chances of absorbing nutrients when their supply is low. Nine root traits of a Brassica napus association panel were explored in hydroponic-system studies under low potassium (K) stress to unravel the genetic basis of root growth in rapeseed. The quantitative trait loci (QTL) and candidate genes for root development were discovered using a multilocus genome-wide association study (ML-GWAS). For the nine traits, a total of 453 significant associated single-nucleotide polymorphism (SNP) loci were discovered, which were then integrated into 206 QTL clusters. There were 45 pleiotropic clusters, and qRTA04-4 and qRTC04-7 were linked to TRL, TSA, and TRV at the same time, contributing 5.25–11.48% of the phenotypic variance explained (PVE) to the root traits. Additionally, 1360 annotated genes were discovered by examining genomic regions within 100 kb upstream and downstream of lead SNPs within the 45 loci. Thirty-five genes were identified as possibly regulating root-system development. As per protein–protein interaction analyses, homologs of three genes (BnaC08g29120D, BnaA07g10150D, and BnaC04g45700D) have been shown to influence root growth in earlier investigations. The QTL clusters and candidate genes identified in this work will help us better understand the genetics of root growth traits and could be employed in marker-assisted breeding for rapeseed adaptable to various conditions with low K levels.

In silico integration of disease resistance QTL, genes and markers with the Brassica juncea physical map

Brassica juncea (AABB), Indian mustard, is a source of disease resistance genes for a wide range of pathogens. The availability of reference genome sequences for B. juncea has made it possible to characterise the genomic structure and distribution of these disease resistance genes. Potentially functional disease resistance genes can be identified by co-localization with genetically mapped disease resistance quantitative trait loci (QTL). Here we identify and characterise disease resistance gene analogs (RGAs), including nucleotide-binding site-leucine-rich repeat (NLR), receptor-like kinase (RLK) and receptor-like protein (RLP) classes, and investigate their association with disease resistance QTL intervals. The molecular genetic marker sequences for four white rust (Albugo candida) disease resistance QTL, six blackleg (Leptosphaeria maculans) disease resistance QTL and BjCHI1, a gene cloned from B. juncea for hypocotyl rot disease, were extracted from previously published studies and used to compare with candidate RGAs. Our results highlight the complications for the identification of functional resistance genes, including the duplicated appearance of genetic markers for several resistance loci, including Ac2(t), AcB1-A4.1, AcB1-A5.1, Rlm6 and PhR2 in both the A and B genomes, due to the presence of homoeologous regions. Furthermore, the white rust loci, Ac2(t) and AcB1-A4.1, mapped to the same position on chromosome A04 and may be different alleles of the same gene. Despite these challenges, a total of nine candidate genomic regions hosting 14 RLPs, 28 NLRs and 115 RLKs were identified. This study facilitates the mapping and cloning of functional resistance genes for applications in crop improvement programs.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01309-5.

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.

Mining of Cloned Disease Resistance Gene Homologs (CDRHs) in Brassica Species and Arabidopsis thaliana

Simple Summary

Developing cultivars with resistance genes (R genes) is an effective strategy to support high yield and quality in Brassica crops. The availability of clone R gene and genomic sequences in Brassica species and Arabidopsis thaliana provide the opportunity to compare genomic regions and survey R genes across genomic databases. In this paper, we aim to identify genes related to cloned genes through sequence identity, providing a repertoire of species-wide related R genes in Brassica crops. The comprehensive list of candidate R genes can be used as a reference for functional analysis.

Abstract

Various diseases severely affect Brassica crops, leading to significant global yield losses and a reduction in crop quality. In this study, we used the complete protein sequences of 49 cloned resistance genes (R genes) that confer resistance to fungal and bacterial diseases known to impact species in the Brassicaceae family. Homology searches were carried out across Brassica napus, B. rapa, B. oleracea, B. nigra, B. juncea, B. carinata and Arabidopsis thaliana genomes. In total, 660 cloned disease R gene homologs (CDRHs) were identified across the seven species, including 431 resistance gene analogs (RGAs) (248 nucleotide binding site-leucine rich repeats (NLRs), 150 receptor-like protein kinases (RLKs) and 33 receptor-like proteins (RLPs)) and 229 non-RGAs. Based on the position and distribution of specific homologs in each of the species, we observed a total of 87 CDRH clusters composed of 36 NLR, 16 RLK and 3 RLP homogeneous clusters and 32 heterogeneous clusters. The CDRHs detected consistently across the seven species are candidates that can be investigated for broad-spectrum resistance, potentially providing resistance to multiple pathogens. The R genes identified in this study provide a novel resource for the future functional analysis and gene cloning of Brassicaceae R genes towards crop improvement.

Multi-omics analysis dissects the genetic architecture of seed coat content in Brassica napus

Background

Brassica napus is an important vegetable oil source worldwide. Seed coat content is a complex quantitative trait that negatively correlates with the seed oil content in B. napus.

Results

Here we provide insights into the genetic basis of natural variation of seed coat content by transcriptome-wide association studies (TWAS) and genome-wide association studies (GWAS) using 382 B. napus accessions. By population transcriptomic analysis, we identify more than 700 genes and four gene modules that are significantly associated with seed coat content. We also characterize three reliable quantitative trait loci (QTLs) controlling seed coat content by GWAS. Combining TWAS and correlation networks of seed coat content-related gene modules, we find that BnaC07.CCR-LIKE (CCRL) and BnaTT8s play key roles in the determination of the trait by modulating lignin biosynthesis. By expression GWAS analysis, we identify a regulatory hotspot on chromosome A09, which is involved in controlling seed coat content through BnaC07.CCRL and BnaTT8s. We then predict the downstream genes regulated by BnaTT8s using multi-omics datasets. We further experimentally validate that BnaCCRL and BnaTT8 positively regulate seed coat content and lignin content. BnaCCRL represents a novel identified gene involved in seed coat development. Furthermore, we also predict the key genes regulating carbon allocation between phenylpropane compounds and oil during seed development in B. napus.

Conclusions

This study helps us to better understand the complex machinery of seed coat development and provides a genetic resource for genetic improvement of seed coat content in B. napus breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02647-5.

Integration of GWAS and transcriptome analyses to identify SNPs and candidate genes for aluminum tolerance in rapeseed (Brassica napus L.)

BACKGROUND

The exchangeable aluminum (Al), released from the acid soils, is another addition to the environmental stress factors in the form of Al toxicity stress. Al stress affects the normal crop development and reduces the overall yield of rapeseed (Brassica napus L.). The response mechanism of plants to Al toxicity is complicated and difficult to understand with few QTL related studies in rapeseed under Al toxicity stress.

RESULT

Using 200,510 SNPs developed by SLAF-seq (specific-locus amplified fragment sequencing) technology, we carried out the genome-wide association analysis (GWAS) in a population of 254 inbred lines of B. napus with large genetic variation and Al-tolerance differences. There were 43 SNPs significantly associated with eight Al-tolerance traits in the seedling stage were detected on 14 chromosomes, and 777 candidate genes were screened at the flanking 100 kb region of these SNPs. Moreover, RNA-seq detected 8291 and 5341 DEGs (the differentially expressed gene) in the Al -tolerant line (ATL) and -sensitive line (ASL), respectively. Based on integration of GWAS and RNA-seq analysis, 64 candidate genes from GWAS analysis differentially expressed at least once in 6 h vs 0 h or 24 h vs 0 h conditions in ATL or ASL. Moreover, four out of sixty-four candidate genes (BnaA03g30320D, BnaA10g11500D, BnaC03g38360D and BnaC06g30030D) were differentially expressed in both 6 h and 24 h compared to 0 h (control) conditions in both lines. The proposed model based on the candidate genes excavated in this study highlighted that Al stress disturb the oxidation-redox balance, causing abnormal synthesis and repair of cell wall and ABA signal transduction, ultimately resulting in inhibition of root elongation.

CONCLUSIONS

The integration of GWAS and transcriptome analysis provide an effective strategy to explore the SNPs and candidate genes, which has a potential to develop molecular markers for breeding Al tolerant rapeseed varieties along with theoretical basis of molecular mechanisms for Al toxicity response of Brassica napus plants.

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.

BnaA03.MKK5-BnaA06.MPK3/BnaC03.MPK3 Module Positively Contributes to Sclerotinia sclerotiorum Resistance in Brassica napus

Brassica napus (oilseed rape) is one of the most important oil crops worldwide, but its growth is seriously threatened by Sclerotinia sclerotiorum. The mechanism of oilseed rape response to this pathogen has rarely been studied. Here, it was identified that BnaA03.MKK5 whose expression was induced by S. sclerotiorum infection was involved in plant immunity. BnaA03.MKK5 overexpression lines exhibited decreased disease symptoms compared to wild-type plants, accompanied by the increased expression of camalexin-biosynthesis-related genes, including BnPAD3 and BnCYP71A13. In addition, two copies of BnMPK3 (BnA06.MPK3 and BnC03.MPK3) were induced by Sclerotinia incubation, and BnaA03.MKK5 interacted with BnaA06.MPK3/BnaC03.MPK3 in yeast. These interactions were confirmed using in vivo co-immunoprecipitation assays. In vitro phosphorylation assays showed that BnaA06.MPK3 and BnaC03.MPK3 were the direct phosphorylation substrates of BnaA03.MKK5. The transgenic oilseed rape plants including BnaA06.MPK3 and BnaC03.MPK3 overexpression lines and BnMPK3 gene editing lines mediated by CRISPR/Cas9 were generated; the results of the genetic transformation of BnaA06.MPK3/BnaC03.MPK3 indicate that BnMPK3 also has a positive role in Sclerotinia resistance. This study provides information about the potential mechanism of B. napus defense against S. Sclerotiorum mediated by a detailed BnaA03.MKK5-BnaA06.MPK3/BnaC03.MPK3 module.

QTL Identification for Stem Fiber, Strength and Rot Resistance in a DH Population from an Alien Introgression of Brassica napus

Stem fiber, stem strength and stem-rot resistance are important agronomic traits in Brassica napus. To understand the molecular mechanism that controls the stem-related traits, we investigated the stem lignin (ADL), cellulose (Cel), hemicellulose (Hem) content, S/G monolignol ratio (SG), stem breaking force (BF), breaking strength (F) and Sclerotinia sclerotiorum resistance (SSR). Each trait was significantly positively or negatively correlated with more than three of the other six traits. QTL mapping for ADL, Cel, Hem, SG, BF, F and SSR were performed using a doubled haploid population derived from an intertribal B. napus introgression line 'Y689' crossed with B. napus cv. 'Westar'. A total of 67 additive QTL were identified and integrated into 55 consensus QTL by meta-analysis. Among the 55 consensus QTL, 23 (41.8%) QTL were co-located and were integrated into 11 unique QTL. The QTL by environment (Q × E) interactions were analyzed and 22 combined QTL were identified. In addition, candidate genes within the QTL intervals were proposed based on the known function of Arabidopsis orthologs. These results provided valuable information for improving lodging resistance, S. sclerotiorum resistance and mechanized harvesting of B. napus.

Genome-wide association of single nucleotide polymorphism loci and candidate genes for frogeye leaf spot (Cercospora sojina) resistance in soybean

BACKGROUND

Frogeye leaf spot (FLS) is a destructive fungal disease that affects soybean production. The most economical and effective strategy to control FLS is the use of resistant cultivars. However, the use of a limited number of resistant loci in FLS management will be countered by the emergence of new high-virulence Cercospora sojina races. Therefore, we identified quantitative trait loci (QTL) that control resistance to FLS and identified novel resistant genes using a genome-wide association study (GWAS) on 234 Chinese soybean cultivars.

RESULTS

A total of 30,890 single nucleotide polymorphism (SNP) markers were used to estimate linkage disequilibrium (LD) and population structure. The GWAS results showed four loci (p < 0.0001) distributed over chromosomes (Chr.) 5 and 20, that are significantly associated with FLS resistance. No previous studies have reported resistance loci in these regions. Subsequently, 45 genes in the two resistance-related haplotype blocks were annotated. Among them, Glyma20g31630 encoding pyruvate dehydrogenase (PDH), Glyma05g28980, which encodes mitogen-activated protein kinase 7 (MPK7), and Glyma20g31510, Glyma20g31520 encoding calcium-dependent protein kinase 4 (CDPK4) in the haplotype blocks deserves special attention.

CONCLUSIONS

This study showed that GWAS can be employed as an effective strategy for identifying disease resistance traits in soybean and narrowing SNPs and candidate genes. The prediction of candidate genes in the haplotype blocks identified by disease resistance loci can provide a useful reference to study systemic disease resistance.

Biologia Futura: progress and future perspectives of long non-coding RNAs in forest trees

Forest trees are affected by climate change, anthropogenic pressure, as well as abiotic and biotic stresses. Conventional tree breeding has so far been limited to enhance overall productivity, and our understanding of the genetic basis of quantitative traits is still inadequate. Quantum leaps in next-generation sequencing technologies and bioinformatics have permitted the exploration and identification of various non-coding regions of the genome other than protein coding genes. These genomic regions produce various types of non-coding RNAs and regulate myriads of biological functions at epigenetic, transcriptional and translational levels. Recently, long non-coding RNAs (lncRNAs) which act as molecular switch have been identified to be pivotal molecules in forest trees. This review focuses on progress made in regulatory mechanisms in various developmental phases like wood formation, adventitious rooting and flowering and stress responses. It was predicted that complex regulatory interactions among lncRNA, miRNA and gene exist. LncRNAs can function as a sponge for miRNAs, reducing the suppressive effect of miRNAs on target mRNAs and perhaps adding a new layer of regulatory interactions among non-coding RNA classes in trees. Furthermore, network analysis revealed the interactions of lncRNA and genes during the expression of several important genes. The insights generated about lncRNAs in forest trees would enable improvement of economically important traits including the devastating abiotic and biotic stresses. In addition, solid understanding on the wide range of regulatory functions of lncRNAs on traits influencing biomass productivity and adaptation would aid the applications of biotechnology in genetic improvement of forest trees.

Genetic Dissection of Mature Root Characteristics by Genome-Wide Association Studies in Rapeseed (Brassica napus L.)

Roots are complicated quantitative characteristics that play an essential role in absorbing water and nutrients. To uncover the genetic variations for root-related traits in rapeseed, twelve mature root traits of a Brassica napus association panel were investigated in the field within three environments. All traits showed significant phenotypic variation among genotypes, with heritabilities ranging from 55.18% to 79.68%. Genome-wide association studies (GWAS) using 20,131 SNPs discovered 172 marker-trait associations, including 103 significant SNPs (-log10 (p) > 4.30) that explained 5.24-20.31% of the phenotypic variance. With the linkage disequilibrium r² > 0.2, these significant associations were binned into 40 quantitative trait loci (QTL) clusters. Among them, 14 important QTL clusters were discovered in two environments and/or with phenotypic contributions greater than 10%. By analyzing the genomic regions within 100 kb upstream and downstream of the peak SNPs within the 14 loci, 334 annotated genes were found. Among these, 32 genes were potentially associated with root development according to their expression analysis. Furthermore, the protein interaction network using the 334 annotated genes gave nine genes involved in a substantial number of interactions, including a key gene associated with root development, BnaC09g36350D. This research provides the groundwork for deciphering B. napus' genetic variations and improving its root system architecture.

SNP- and Haplotype-Based GWAS of Flowering-Related Traits in Brassica napus

Traits related to flowering time are the most promising agronomic traits that directly impact the seed yield and oil quality of rapeseed (Brassica napus L.). Developing early flowering and maturity rapeseed varieties is an important breeding objective in B. napus. Many studies have reported on days to flowering, but few have reported on budding, bolting, and the interval between bolting and DTF. Therefore, elucidating the genetic architecture of QTLs and genes regulating flowering time, we presented an integrated investigation on SNP and haplotype-based genome-wide association study of 373 diverse B. napus germplasm, which were genotyped by the 60K SNP array and were phenotyped in the four environments. The results showed that a total of 15 and 37 QTLs were detected from SNP and haplotype-based GWAS, respectively. Among them, seven QTL clusters were identified by haplotype-based GWAS. Moreover, three and eight environmentally stable QTLs were detected by SNP-GWAS and haplotype-based GWAS, respectively. By integrating the above two approaches and by co-localizing the four traits, ten (10) genomic regions were under selection on chromosomes A03, A07, A08, A10, C06, C07, and C08. Interestingly, the genomic regions FT.A07.1, FT.A08, FT.C06, and FT.C07 were identified as novel. In these ten regions, a total of 197 genes controlling FT were detected, of which 14 highly expressed DEGs were orthologous to 13 Arabidopsis thaliana genes after integration with transcriptome results. In a nutshell, the above results uncovered the genetic architecture of important agronomic traits related to flowering time and provided a basis for multiple molecular marker-trait associations in B. napus.

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.

Sclerotinia sclerotiorum Thioredoxin1 (SsTrx1) is required for pathogenicity and oxidative stress tolerance

Sclerotinia sclerotiorum infects host plant tissues by inducing necrosis to source nutrients needed for its establishment. Tissue necrosis results from an enhanced generation of reactive oxygen species (ROS) at the site of infection and apoptosis. Pathogens have evolved ROS scavenging mechanisms to withstand host-induced oxidative damage. However, the genes associated with ROS scavenging pathways are yet to be fully investigated in S. sclerotiorum. We selected the S. sclerotiorum Thioredoxin1 gene (SsTrx1) for our investigations as its expression is significantly induced during S. sclerotiorum infection. RNA interference-induced silencing of SsTrx1 in S. sclerotiorum affected the hyphal growth rate, mycelial morphology, and sclerotial development under in vitro conditions. These outcomes confirmed the involvement of SsTrx1 in promoting pathogenicity and oxidative stress tolerance of S. sclerotiorum. We next constructed an SsTrx1-based host-induced gene silencing (HIGS) vector and mobilized it into Arabidopsis thaliana (HIGS-A) and Nicotiana benthamiana (HIGS-N). The disease resistance analysis revealed significantly reduced pathogenicity and disease progression in the transformed genotypes as compared to the nontransformed and empty vector controls. The relative gene expression of SsTrx1 increased under oxidative stress. Taken together, our results show that normal expression of SsTrx1 is crucial for pathogenicity and oxidative stress tolerance of S. sclerotiorum.

Genomic background selection to reduce the mutation load after random mutagenesis

Random mutagenesis is a standard procedure to increase allelic variation in a crop species, especially in countries where the use of genetically modified crops is limited due to legal constraints. The chemical mutagen EMS is used in many species to induce random mutations throughout the genome with high mutation density. The major drawback for functional analysis is a high background mutation load in a single plant that must be eliminated by subsequent backcrossing, a time and resource-intensive activity. Here, we demonstrate that genomic background selection combined with marker-assisted selection is an efficient way to select individuals with reduced background mutations within a short period. We identified BC1 plants with a significantly higher share of the recurrent parent genome, thus saving one backcross generation. Furthermore, spring rapeseed as the recurrent parent in a backcrossing program could accelerate breeding by reducing the generation cycle. Our study depicts the potential for reducing the background mutation load while accelerating the generation cycle in EMS-induced winter oilseed rape populations by integrating genomic background selection.

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

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

Overexpression of Cinnamoyl-CoA Reductase 2 in Brassica napus Increases Resistance to Sclerotinia sclerotiorum by Affecting Lignin Biosynthesis

Sclerotinia sclerotiorum causes severe yield and economic losses for many crop and vegetable species, especially Brassica napus. To date, no immune B. napus germplasm has been identified, giving rise to a major challenge in the breeding of Sclerotinia resistance. In the present study, we found that, compared with a Sclerotinia-susceptible line (J902), a Sclerotinia-resistant line (J964) exhibited better xylem development and a higher lignin content in the stems, which may limit the invasion and spread of S. sclerotiorum during the early infection period. In addition, genes involved in lignin biosynthesis were induced under S. sclerotiorum infection in both lines, indicating that lignin was deposited proactively in infected tissues. We then overexpressed BnaC.CCR2.b, which encodes the first rate-limiting enzyme (cinnamoyl-CoA reductase) that catalyzes the reaction of lignin-specific pathways, and found that overexpression of BnaC.CCR2.b increased the lignin content in the stems of B. napus by 2.28-2.76% under normal growth conditions. We further evaluated the Sclerotinia resistance of BnaC.CCR2.b overexpression lines at the flower-termination stage and found that the disease lesions on the stems of plants in the T₂ and T₃ generations decreased by 12.2-33.7% and 32.5-37.3% compared to non-transgenic control plants, respectively, at 7days post-inoculation (dpi). The above results indicate that overexpression of BnaC.CCR2.b leads to an increase in lignin content in the stems, which subsequently leads to increased resistance to S. sclerotiorum. Our findings demonstrate that increasing the lignin content in the stems of B. napus is an important strategy for controlling Sclerotinia.

A Plant Stress-Responsive Bioreporter Coupled With Transcriptomic Analysis Allows Rapid Screening for Biocontrols of Necrotrophic Fungal Pathogens

Streptomyces are soil-borne Actinobacteria known to produce a wide range of enzymes, phytohormones, and metabolites including antifungal compounds, making these microbes fitting for use as biocontrol agents in agriculture. In this study, a plant reporter gene construct comprising the biotic stress-responsive glutathione S-transferase promoter GSTF7 linked to a luciferase output (GSTF7:luc) was used to screen a collection of Actinobacteria candidates for manipulation of plant biotic stress responses and their potential as biocontrol agents. We identified a Streptomyces isolate (KB001) as a strong candidate and demonstrated successful protection against two necrotrophic fungal pathogens, Sclerotinia sclerotiorum and Rhizoctonia solani, but not against a bacterial pathogen (Pseudomonas syringe). Treatment of Arabidopsis plants with either KB001 microbial culture or its secreted compounds induced a range of stress and defense response-related genes like pathogenesis-related (PR) and hormone signaling pathways. Global transcriptomic analysis showed that both treatments shared highly induced expression of reactive oxygen species and auxin signaling pathways at 6 and 24 h posttreatment, while some other responses were treatment specific. This study demonstrates that GSTF7 is a suitable marker for the rapid and preliminary screening of beneficial bacteria and selection of candidates with potential for application as biocontrols in agriculture, including the Streptomyces KB001 that was characterized here, and could provide protection against necrotrophic fungal pathogens.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Exploring the gene pool of Brassica napus by genomics-based approaches

De novo allopolyploidization in Brassica provides a very successful model for reconstructing polyploid genomes using progenitor species and relatives to broaden crop gene pools and understand genome evolution after polyploidy, interspecific hybridization and exotic introgression. B. napus (AACC), the major cultivated rapeseed species and the third largest oilseed crop in the world, is a young Brassica species with a limited genetic base resulting from its short history of domestication, cultivation, and intensive selection during breeding for target economic traits. However, the gene pool of B. napus has been significantly enriched in recent decades that has been benefit from worldwide effects by the successful introduction of abundant subgenomic variation and novel genomic variation via intraspecific, interspecific and intergeneric crosses. An important question in this respect is how to utilize such variation to breed crops adapted to the changing global climate. Here, we review the genetic diversity, genome structure, and population-level differentiation of the B. napus gene pool in relation to known exotic introgressions from various species of the Brassicaceae, especially those elucidated by recent genome-sequencing projects. We also summarize progress in gene cloning, trait-marker associations, gene editing, molecular marker-assisted selection and genome-wide prediction, and describe the challenges and opportunities of these techniques as molecular platforms to exploit novel genomic variation and their value in the rapeseed gene pool. Future progress will accelerate the creation and manipulation of genetic diversity with genomic-based improvement, as well as provide novel insights into the neo-domestication of polyploid crops with novel genetic diversity from reconstructed genomes.

Proteome Changes in Stem Tissues of Sunflower Lines Inoculated with Culture Filtrate of Sclerotinia sclerotiorum

BACKGROUND

Sclerotinia sclerotiorum (Lib.) de Bary cause a deleterious disease on sunflower plants. Oxalic acid is the main pathogenicity factor of S. sclerotiorum. Two dimensional gel electrophoresis and mass spectrometry have been used in several studies to investigate molecular changes that occur in the plants in response to S. sclerotiorum infection. Comparing responses of resistant and susceptible lines upon pathogen infection provided novel information regarding defense mechanisms against this necrotrophic pathogen.

OBJECTIVES

The present study reports proteome changes of partially resistant and susceptible sunflower lines under pathogen's culture filtrate treatment, resulting in the characterization of up- and down- regulated proteins.

MATERIAL AND METHODS

Sunflower partially resistant and susceptible lines with two true leaves were exposed to fungus culture filtrate. The stems of treated and untreated plants were sampled at 24, 48 and 72 hours after treatment for two-dimensional electrophoresis. Twenty spots showed more than 1.5-fold change in abundance were subjected to MALDI/TOF-TOF MS for further analysis.

RESULTS

The identified proteins were categorized into several classes including carbohydrate and energy metabolism (25%), cellular metabolic process (15%), stress response (15%), plant cell wall biogenesis (10%), photosynthesis (10%), protein metabolism (10%), unknown function (10%) and redox homeostasis (5%).

CONCLUSIONS

Our proteomic investigation demonstrates an increase in the expression of proteins only in partially resistant line, such as proteins involved in carbohydrate metabolism and plant defense responses (malate dehydrogenase and peroxidase), metabolic process (adenosine kinase), regulating cell redox homeostasis (disulfide isomerase) and lignin biosynthetic process (laccase). Moreover, the expression of pyrroline-5-carboxylate reductase, involved in proline biosynthesis, was significantly changed in both sunflower lines in response to pathogen culture filtrate. Proteins which were only up-regulated in the partially resistant lines might have a significant role in mediating the defense against Sclerotinia and could be considered for enhancing resistance against this devastating pathogen.

MicroRNA319a regulates plant resistance to Sclerotinia stem rot

MicroRNA319a (miR319a) controls cell division arrest in plant leaves by inhibiting the expression of TCP (TEOSINTE BRANCHED 1/CYCLOIDEA/PCF) family genes. However, it is unclear whether miR319a influences infection by necrotrophic pathogens and host susceptibility. In this study, we revealed that miR319a affects plant resistance to stem rot disease caused by Sclerotinia sclerotiorum. In Brassica rapa plants infected with S. sclerotiorum, miR319a levels increased while the expression levels of several BraTCP genes significantly decreased compared with those of uninfected plants. Overexpression of BraMIR319a in B. rapa increased the susceptibility of the plants to S. sclerotiorum and aggravated stem rot disease, whereas overexpression of BraTCP4-1 promoted plant resistance. RNA sequencing data revealed a potential relationship between miR319a and pathogen-related WRKY genes. Chromatin immunoprecipitation, electrophoretic mobility shift, and reporter transaction assays showed that BraTCP4-1 could bind to the promoters of WRKY75, WRKY70, and WRKY33 and directly activate these pathogen-related genes. Moreover, the expression levels of WRKY75, WRKY70, and WRKY33 in plants overexpressing BraMIR319a decreased significantly, whereas those of plants overexpressing BraTCP4-1 increased significantly, relative to the wild type. These results suggest that miR319a and its target gene BraTCP4 control stem rot resistance through pathways of WRKY genes.

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.

Genome-wide association study in hexaploid wheat identifies novel genomic regions associated with resistance to root lesion nematode (Pratylenchus thornei)

Root lesion nematode (RLN; Pratylenchus thornei) causes extensive yield losses in wheat worldwide and thus pose serious threat to global food security. Reliance on fumigants (such as methyl bromide) and nematicides for crop protection has been discouraged due to environmental concerns. Hence, alternative environment friendly control measures like finding and deployment of resistance genes against Pratylenchus thornei are of significant importance. In the present study, genome-wide association study (GWAS) was performed using single-locus and multi-locus methods. In total, 143 wheat genotypes collected from pan-Indian wheat cultivation states were used for nematode screening. Genotypic data consisted of  > 7K SNPs with known genetic positions on the high-density consensus map was used for association analysis. Principal component analysis indicated the existence of sub-populations with no major structuring of populations due to the origin. Altogether, 25 significant marker trait associations were detected with - log10 (p value) > 4.0. Three large linkage disequilibrium blocks and the corresponding haplotypes were found to be associated with significant SNPs. In total, 37 candidate genes with nine genes having a putative role in disease resistance (F-box-like domain superfamily, Leucine-rich repeat, cysteine-containing subtype, Cytochrome P450 superfamily, Zinc finger C2H2-type, RING/FYVE/PHD-type, etc.) were identified. Genomic selection was conducted to investigate how well one could predict the phenotype of the nematode count without performing the screening experiments. Prediction value of r = 0.40 to 0.44 was observed when 56 to 70% of the population was used as a training set. This is the first report where GWAS has been conducted to find resistance against root lesion nematode (P. thornei) in Indian wheat germplasm.

The Resistance of Narrow-Leafed Lupin to Diaporthe toxica Is Based on the Rapid Activation of Defense Response Genes

Narrow-leafed lupin (Lupinus angustifolius L.) is a grain legume crop that is advantageous in animal nutrition due to its high protein content; however, livestock grazing on stubble may develop a lupinosis disease that is related to toxins produced by a pathogenic fungus, Diaporthe toxica. Two major unlinked alleles, Phr1 and PhtjR, confer L. angustifolius resistance to this fungus. Besides the introduction of these alleles into modern cultivars, the molecular mechanisms underlying resistance remained unsolved. In this study, resistant and susceptible lines were subjected to differential gene expression profiling in response to D. toxica inoculation, spanning the progress of the infection from the early to latent phases. High-throughput sequencing of stem transcriptome and PCR quantification of selected genes were performed. Gene Ontology term analysis revealed that an early (24 h) response in the resistant germplasm encompassed activation of genes controlling reactive oxygen species and oxylipin biosynthesis, whereas in the susceptible germplasm, it comprised induction of xyloglucan endotransglucosylases/hydrolases. During the first five days of the infection, the number of genes with significantly altered expressions was about 2.6 times higher in resistant lines than in the susceptible line. Global transcriptome reprogramming involving the activation of defense response genes occurred in lines conferring Phr1 and PhtjR resistance alleles about 4–8 days earlier than in the susceptible germplasm.

Revisiting Plant-Microbe Interactions and Microbial Consortia Application for Enhancing Sustainable Agriculture: A Review

The present scenario of agricultural sector is dependent hugely on the use of chemical-based fertilizers and pesticides that impact the nutritional quality, health status, and productivity of the crops. Moreover, continuous release of these chemical inputs causes toxic compounds such as metals to accumulate in the soil and move to the plants with prolonged exposure, which ultimately impact the human health. Hence, it becomes necessary to bring out the alternatives to chemical pesticides/fertilizers for improvement of agricultural outputs. The rhizosphere of plant is an important niche with abundant microorganisms residing in it. They possess the properties of plant growth promotion, disease suppression, removal of toxic compounds, and assimilating nutrients to plants. Utilizing such beneficial microbes for crop productivity presents an efficient way to modulate the crop yield and productivity by maintaining healthy status and quality of the plants through bioformulations. To understand these microbial formulation compositions, it becomes essential to understand the processes going on in the rhizosphere as well as their concrete identification for better utilization of the microbial diversity such as plant growth–promoting bacteria and arbuscular mycorrhizal fungi. Hence, with this background, the present review article highlights the plant microbiome aboveground and belowground, importance of microbial inoculants in various plant species, and their subsequent interactive mechanisms for sustainable agriculture.

Analysis and Identification of QTL for Resistance to Sclerotinia sclerotiorum in Pea (Pisum sativum L.)

White mold caused by Sclerotinia sclerotiorum is an important constraint to field pea (Pisum sativum L.) production worldwide. To transfer white mold resistance into an adapted background, and study the genetics of the disease, two recombinant inbred line (RIL) populations (PRIL17 and PRIL19) were developed by crossing two partially resistant plant introductions with two susceptible pea cultivars. PRIL17 (Lifter × PI240515), and PRIL19 (PI169603 × Medora) were evaluated for resistance to white mold by measuring lesion expansion inhibition (LEI) and nodal transmission inhibition (NTI) at 3, 7, and 14 days post inoculation (dpi) under controlled environmental conditions. Lesion expansion inhibition percentage (LEIP), survival rate (SR), and area under disease progress curves (AUDPC) were also calculated accordingly. Because of a positive correlation between LEI and NTI with height, short and long internode individuals of each population were analyzed separately to avoid any confounding effect of height to pathogen response. A total of 22 short genotypes demonstrated partial resistance based on at least two Porter's resistance criteria. Only two pea genotypes with partial resistance to white mold (PRIL19-18 and PRIL19-124) had both semi-leafless (afila) and short internode traits. Both the RIL populations were genotyped using genotyping by sequencing (GBS). For PRIL17 and PRIL19, genetic maps were constructed from a total of 1,967 and 1,196 single nucleotide polymorphism (SNP) and spanned over 1,494 cM and 1,415 cM representing seven and nine linkage groups, respectively. A consensus map constructed using data from both populations, had 1,486 unique SNPs over 2,461 cM belonging to seven linkage groups. Inclusive composite interval mapping (ICIM) identified thirteen quantitative trait loci (QTL) associated with white mold resistance traits in both populations. Three of them were co-located with height genes (a morphological trait that reduces infection risk and acts as disease avoidance) and the other ten QTL were associated with two forms of physiological resistance (seven for LEI and three for NTI) with LOD and r2 ranging from 3.0 to 28.5 and 5.1 to 64.3, respectively. The development of resistance lines, genetic dissection and identification of markers associated will help accelerate breeding efforts for white mold resistance using molecular breeding approaches.