Genetic Diversity Studies Based on Morphological Variability, Pathogenicity and Molecular Phylogeny of the Sclerotinia sclerotiorum Population From Indian Mustard (Brassica juncea)

White mold or stem rot disease are ubiquitously distributed throughout the world and the causal organism of this disease Sclerotinia sclerotiorum (Lib.) de Bary, is known to infect over 400 plant species. Sclerotinia stem rot is one of the most devastating fungal diseases and poses a serious threat to the worldwide cultivation of oilseed Brassica including India. S. sclerotiorum pathogen usually infects the stem but in severe cases leaves and pods also affected at different developmental stages that deteriorate not only the oil quality but also causing the seed and oil yield losses up to 90% depending on the severity of the disease infestation. This study investigated the morphological and molecular characterization of pathogenic S. sclerotiorum (Lib) de Bary geographical isolates from oilseed Brassica including Brassica juncea (Indian mustard). The aim of this study was to compare isolates of S. sclerotiorum originated from different agro-climatic conditions and to analyse similarity or differences between them as well as to examine the virulence of this pathogen specifically in Brassica for the first time. The collection of S. sclerotiorum isolates from symptomatic Brassica plants was done and analyzed for morphological features, and molecular characterization. The virulence evaluation test of 65 isolates on four Brassica cultivars has shown 5 of them were highly virulent, 46 were virulent and 14 were moderately virulent. Phylogenetic analysis encompassing all the morphological features, SSR polymorphism, and ITS sequencing has shown the existence of high genetic diversity among the isolates that categorized all the isolates in three evolutionary lineages in the derived dendrogram. Further, genetic variability analysis based on sequences variation in ITS region of all the isolates has shown the existence of either insertions or deletions of the nucleotides in the ITS region has led to the interspecies variability and observed the variation were in a clade-specific manner. Together this analysis observed the existence of higher heterogeneity and genetic variability in S. sclerotiorum isolates collection and indicates the presence of clonal and sexual progenies of the pathogen in the mustard growing regions of India surveyed in this study. With a higher level of genetic variability and diversity among the S. sclerotiorum population needs robust screening approaches to identify the donor parent and utilize them in resistance breeding program for effectively counter the menace of stem rot disease in Brassica.

Genetic Diversity Studies Based on Morphological Variability, Pathogenicity and Molecular Phylogeny of the Sclerotinia sclerotiorum Population From Indian Mustard (Brassica juncea) White mold or stem rot disease are ubiquitously distributed throughout the world and the causal organism of this disease Sclerotinia sclerotiorum (Lib.) de Bary, is known to infect over 400 plant species. Sclerotinia stem rot is one of the most devastating fungal diseases and poses a serious threat to the worldwide cultivation of oilseed Brassica including India. S. sclerotiorum pathogen usually infects the stem but in severe cases leaves and pods also affected at different developmental stages that deteriorate not only the oil quality but also causing the seed and oil yield losses up to 90% depending on the severity of the disease infestation. This study investigated the morphological and molecular characterization of pathogenic S. sclerotiorum (Lib) de Bary geographical isolates from oilseed Brassica including Brassica juncea (Indian mustard). The aim of this study was to compare isolates of S. sclerotiorum originated from different agro-climatic conditions and to analyse similarity or differences between them as well as to examine the virulence of this pathogen specifically in Brassica for the first time. The collection of S. sclerotiorum isolates from symptomatic Brassica plants was done and analyzed for morphological features, and molecular characterization. The virulence evaluation test of 65 isolates on four Brassica cultivars has shown 5 of them were highly virulent, 46 were virulent and 14 were moderately virulent. Phylogenetic analysis encompassing all the morphological features, SSR polymorphism, and ITS sequencing has shown the existence of high genetic diversity among the isolates that categorized all the isolates in three evolutionary lineages in the derived dendrogram. Further, genetic variability analysis based on sequences variation in ITS region of all the isolates has shown the existence of either insertions or deletions of the nucleotides in the ITS region has led to the interspecies variability and observed the variation were in a clade-specific manner. Together this analysis observed the existence of higher heterogeneity and genetic variability in S. sclerotiorum isolates collection and indicates the presence of clonal and sexual progenies of the pathogen in the mustard growing regions of India surveyed in this study. With a

INTRODUCTION
Globally India continues to be at a 3rd position after Canada and China in acreage (19.3%) and after China and Canada in production (11.1%) of rapeseed-mustard. In India, among nine edible oilseed crops, the share of rapeseed-mustard is about one-fourth of total area and one-third of total oil production in the country. During 2015-2016, production (6.82 mt) and productivity (1184 kg/ha) was achieved (Anonymous, 2016). Rapeseed-mustard is the major source of income especially for the marginal and small farmers in rainfed areas which are about 25% of the total cultivated area. In spite of its increase in demand for the year the production of oilseed Brassica remains to stagnate over the year and most of the demands are being met through import from outside the India. The main reason behind productivity stagnation in Indian Brassica is its susceptibility and damages caused to the crop by various insect pests and disease infestation in addition to the other yield-limiting factors. Out of thirty diseases known to infest the Brassica crops in India, stem rot has been found one of the most devastating diseases that heavily damages the crops during the flowering stage of development. The stem rot disease which is caused by fungal pathogens, Sclerotinia sclerotiorum (Lib) de Bary, ubiquitously found throughout the world is a polyphagous, soil-borne plant pathogen that infects more than 400 plant species of diverse phylogenetic origin (Boland and Hall, 1994;Saharan and Mehta, 2008;Sharma et al., 2015). In India, during the eighties and nineties, the stem rot (SR) disease in rapeseed-mustard was of a minor importance, because of its seldom appearance over the ground level of the isolated plants after mycelial infection. A widely adopted monocropping practices and cultivation of rapeseed-mustard under irrigated condition has significantly increased the sclerotial population in the soil that has made SR very serious disease of oilseed Brassica crops in states including Rajasthan, Haryana, Punjab, Uttar Pradesh, Bihar, Assam, West Bengal and Madhya Pradesh (Sharma et al., 2015). This fungus has been long considered as prototypical necrotrophs as it begins highly pathogenic phase by releasing oxalic acids and cellulolytic enzymes immediately upon host cuticle penetration followed by mycelial proliferation inside the host cell followed by a saprophytic phase that supports the sclerotia formation (Hegedus and Rimmer, 2005). However, the recent studies decipher the fact of evidence for the occurrence of a brief biotrophic phase just within the apoplastic space next after the establishment of the host-pathogen connection and hence based on these it is more appropriately classified as a hemibiotroph (Kabbage et al., 2015). The information related to the genetic diversity of the pathogen and their effective virulence over the target crop is the foremost requirement for taking the breeding program for development of pathogen resistance in the release of the regionspecific cultivars. Various diversity analysis tools based on the molecular methods like microsatellite haplotype (Aldrich-Wolfe et al., 2015), SSR (simple sequence repeat or microsatellitebased marker; Meinhardt et al., 2002), AFLP (amplified fragment length polymorphism; Cubeta et al., 1997), and SRAP (sequencerelated amplified polymorphism technique; Li et al., 2009) have been used in analyzing the genetic diversity of the pathogen S. sclerotiorum from different host species. Very limited variability was observed in ITS (internal transcribed spacer) sequences in S. sclerotiorum isolates from various host species (Njambere et al., 2008) and thus a universal barcode markers were developed from the nearly conserved nature of the nuclear ribosomal internal transcribed spacer (ITS) region for imparting the individual identity to the fungus up to genus level (Schoch et al., 2012). Furthermore, MCGs is another diversity analysis method based on the mycelial compatibility grouping (MCG) has been used in establishing the kinship among S. sclerotiorum isolates from chickpea (Kull et al., 2004;Li et al., 2008). In addition to it, the diversity based on the morphological appearance of sclerotia, mycelial growth, and ascospores formation have also been reported in analyzing the genetic diversity of S. sclerotiorum isolates in previous studies (Li et al., 2008;Sharma et al., 2013). However, polymorphism and genetic diversity of the S. sclerotiorum isolates at the morphological and DNA sequence level has not been comprehensively studied so far especially for the isolates from Brassica species of India.
Being a polyphagous nature, S. sclerotiorum pathogen is usually infecting not only the majority of the economically important dicotyledonous species but also serve as the major pathogen for several monocotyledonous plant species (Boland and Hall, 1994). The yield loss estimated with the S. sclerotiorum infestation has exceeded hundred million dollars annually because of the lack of resistance cultivar of the crop species and also because of lack of the effective management practices (Tok et al., 2016). In general, the management of S. sclerotiorum borne disease is not much easy because of its widespread existence, irregular incidence, and the long-term survival by producing huge numbers of sclerotia in the soil. Although several control measures like chemical and cultural methods have been devised and adopted for countering the Sclerotinia stem rot menace (Rousseau et al., 2007) none of them were found fully effective in preventing either the process of disease infection or pathogenesis progression after infection. The extent of genetic diversity of the pathogen and their widespread distribution among host species across the growing regions play an important role in determining and devising the control strategies to efficiently control the diseases in the more effective way. Hence, the availability of the genetic variability information related to the target pathogen is the foremost requirement for designing the management means to counter the disease incidence more effectively. In plant-pathogen interactions, development of new pathogenic races, and the breakdown of host resistance are the limiting factors in resistance deployment against plant diseases. The pathogen's life history characteristics and evolutionary potential are major factors leading to the pathogen overcoming host resistance. Therefore, major efforts should be focused not only in understanding the genetic structure of the fungal populations but also to determine how populations will evolve in response to different control strategies (McDonald and Linde, 2002).
In recent past, India observed the frequent incidence of the S. sclerotiorum infestation on cereals and horticultural species that draws the wide attention of the researcher over this fungal pathogen. In pursuance of basic understanding about pathogenicity, diversity and distribution pattern of the S. sclerotiorum were extensively studied on isolates collected from various host species like chickpea (Mandal and Dubey, 2012), vegetable crops (Choudhary and Prasad, 2012), cumin (Prasad et al., 2017), carnation (Kumar et al., 2015), oilseed Brassica (Sharma et al., 2013), and their identity were established on the basis of morphological features and cultural conditions. However, major variation in the growth characteristics has been reported in the growing collection of S. sclerotiorum isolates even as they belong to the same Sclerotinia species. Indeed, as projected with diversity analysis S. sclerotiorum isolates has been reported of possessing the variation in morpho-physiological, biochemical properties, molecular features and pathogenicity in terms of virulence because of the presence of clonal and sexual progenies together even in the same crop and in same region (Atallah et al., 2004;Sexton and Howlett, 2004;Irani et al., 2011). For establishing the console features of the pathogen especially in Brassica growing regions of India, the present study was aimed at determining the genetic diversity within S. sclerotiorum population from various Brassica species based on morphological characteristics, genotyping with simple sequence repeats (SSR) markers and molecular phylogeny by ITS sequence analysis.

Sclerotinia sclerotiorum Isolates
Brassica growing areas in 10 states of India (Rajasthan, Haryana, Punjab, Delhi, U.P., Bihar, Uttarakhand, Himachal Pradesh, Jammu & Kashmir, and Jharkhand) were surveyed and stem rot disease infected plants were collected from 65 different locations ( Table 1). The sclerotia obtained from the stem rot infested Indian mustard plant samples were first washed with sterile water than surface sterilized with 70% ethanol for 2 min and again washed two times with sterile water. The drained sclerotia over pre-sterilized filter papers were placed on nutrient media (Potato Dextrose Agar; PDA) plates supplemented with 50 µg/ml tetracycline antibiotics to prevent the growth of bacterial contamination. The samples were wrapped in a brown envelope and kept for incubation at 20 • ± 2 • C in dark for 4-5 days. After the development of the white fluffy mass of mycelial growth of S. sclerotiorum, the mycelial plaques were used to subculture the isolates on PDA slants and pure culture of them was stored at 4 • C for future use. Morphological identification of the isolates was based on cultural characteristics of the S. sclerotiorum and morphology of the mycelial mat and sclerotia formation traits.

Morphological Characterization of S. sclerotiorum Isolates
Freshly grown 3-4 days old cultures of Sclerotinia isolates were used for analyzing the morphological features like the radial spread of mycelial growth (mm) at 72 and 96 h after inoculation. Whereas, the number of sclerotia produced from each inoculum per Petri plates (90 mm), length and diameter of the sclerotia (mm) and weight of individual sclerotia (mg) was measured 10-15 days after inoculation.

Pathogenicity of Isolates
An experiment to study pathogenicity and pathogenic variability among 65 geographical isolates of S. sclerotiorum were conducted during 2016-2017 post-rainy season at ICAR-Directorate of Rapeseed-Mustard Research, Bharatpur, India. Seven different oilseed Brassica species were selected for this study. The crop was destroyed after the experiment by cutting and collecting the debris followed by autoclaving before disposing of them.

Inoculation
The inoculum was mass multiplied in the laboratory on autoclaved sorghum grains in glass jars and incorporated into the soil prior to sowing. Further, the plants were inoculated at 60 days after sowing on the stem with the pathogen growing on agar blocks. Stem inoculation procedure was followed as described by Buchwaldt et al. (2003). A single 5 mm mycelial bit cut from S. sclerotiorum colony of 4-5 days old culture growing on potato dextrose agar was used to inoculate each plant. The mycelial bit along with cotton swab soaked in sterilized distilled water was placed on a small piece of parafilm (5-7 cm). Mycelial bit touching the stem at 15 cm height was then secured by wrapping the parafilm strip around the stem. Wet cotton swab maintained high humidity during the infection period (Sharma et al., 2013).

Disease
Disease incidence was assessed by recording the size of stem lesion length (cm) and disease severity (%) 15-21 days after inoculation. This has been shown to be an ideal time to demonstrate the host response to the pathogen (Li et al., 2007). The disease parameters recorded were statistically verified by data analysis using Duncan's test (P < 0.05) and analysis of variance (ANOVA).

DNA Extraction
Sclerotinia sclerotiorum isolates were grown on PDA medium at 20 • ± 2 • C in BOD incubator for 5 days. Mycelial mat was harvested by scraping with a sterile spatula and ∼2 g of the mycelial mat was used for DNA extraction using modified CTAB method. The fine powder of the grounded mycelial mat in liquid nitrogen was transferred in preheated 10 ml DNA extraction  The severity of disease assessed on Brassica cultivars was measured in lesion length (cm). Disease index for the strains listed, accessed on length of the lesion developed on stem (in cm). Lesion size 0.0 cm: Non-virulent; 0.1-3.0 cm: less virulent; 3.1-10.0: Moderately virulent; 10.1-20.0 cm: virulent; >20.0 cm: highly virulent. NRCDR-02 stands for National Research Centre Disease Resistance-02, YS stands for yellow sarson buffer (1% Cetyltrimethylammonium bromide; 250 mM NaCl; 100 mM Tris-HCl, pH 8.0; 100 mM EDTA, pH 8.0; along with 4 mM Spermidine) and incubated at 65 • C with intermittent mixing at regular interval for 60 min. Denatured proteins were removed by extracting once with an equal volume of Trissaturated phenol: chloroform: isoamyl alcohol (25:24:1, v/v/v), followed by repeated extractions with an equal volume of Trissaturated chloroform: isoamyl alcohol (24:1, v/v). The aqueous phase obtained after centrifugation was added with 0.6 volume of chilled isopropanol and stored overnight at −20 • C for DNA precipitation. DNA was pelleted, dried and resuspended in 500 µl of Tris-EDTA buffer [10 mM Tris-HCl, pH 7.5; 1 mM EDTA (ethylenediaminetetraacetic acid)]. The purity and concentration of the DNA were determined through NanoDrop (Thermo Scientific, USA) and stored in aliquots at −20 • C.

Species-Specific PCR Assay
Sclerotinia-specific primers as described by Freeman et al. (2002) ( Table 2) were used to carry out PCR amplification in all the isolates. The PCR reaction of 15 µL contained 1.5 µL of 10x PCR buffer with 15 mM MgCl 2, 50 µM dNTPs, 10 µM forward and reverse primers, 1 U Taq DNA polymerase (Bangalore Genei, India) and 50 ng of template DNA was performed in a thermal cycler (Biometra, ILS, USA) and the program made with initial denaturation at 95 • C for 4 min followed by 35 cycles of denaturation at 94 • C for 45 s, annealing at 57 • C for 45 s and  White et al., 1990 ITS5 GGAAGTAAAAGTCGTAACAAGG extension at 72 • C for 1 min, and a final extension at 72 • C for 10 min. The PCR amplicons were resolved on 1% agarose gel along with 1 kb standard DNA ladder and visualized in a UV transilluminator. Gels were photographed using Alpha image (Alpha innotech, USA) Gel Doc system.

Genetic Diversity
Simple Sequence Repeat (SSR) Marker Analysis The pathogenicity tested 65 S. sclerotiorum isolates were subjected to SSR fingerprinting. A total of 25 SSR primers (Sirjusingh and Kohn, 2001; Table 3) were synthesized from IDT, Agrigenome, Bengaluru, India and used for amplification of microsatellite loci. PCR reactions were performed using a BioRad C1000 thermocycler (BioRad, USA). The PCR reaction was set up in a 25 µL reaction volume with 2.5 µL of 10x Taq buffer (with 15 mM MgCl 2 ), 0.25 mM dNTPs (10 mM dNTPs mix), 10 nM of forward and reverse primers, 1 U of Taq polymerase and 200 ng template DNA. The reproducibility of the amplification was confirmed by repeated PCR with a similar set of ingredients and condition. For each experiment, negative controls were taken with sterile water in place of the template. The thermal cycler program consisted of an initial denaturation for 4 min at 95 • C, followed by 35 cycles of denaturation at 94 • C for 45 s, optimized annealing temperature for 45 s, and extension at 72 • C for 1 min with a final extension of 72 • C for 10 min. The PCR amplicons were resolved on 3.5% agarose gel in 1x TAE buffer by electrophoresis at 60 V cm −1 .

SSR Profiling and Data Analysis:
The standard binomial matrices were followed in SSR profiling and the presence and absence of the allele at a particular locus was scored as 1 and 0, respectively. Though Jaccard's coefficient pairwise distance were calculated and the resultant distance matrices were further used in NTSYSpc2 software version 2.0 (Rohlf, 2000) for clustering the isolates under investigation based on the Unweighted Pair Group method using arithmetic means (UPGMA) method. Allele frequency and Polymorphism information content (PIC) value for each marker was also calculated based on the formula: Hn = 1-Σpi2, where pi is the allele frequency of the ith allele (Nei, 1973).

Internal Transcribed Spacer (ITS) Amplification
The universal primers pair ITS4 and ITS5 (White et al., 1990; Table 2) were used for amplifying the ITS region of the rDNA.
The PCR was carried out in a total volume of 50 µL with 1xPCR buffer (with 15 mM MgCl 2 ), 250 µM each dNTPs, 1 µL of each primer (10 nmol), 1.5 U Taq DNA polymerase, and 200 ng of genomic DNA. The thermal cycler program consisted of an initial denaturation at 95 • C for 4 min followed by cycling conditions included denaturation at 94 • C for 45 s, annealing at 49 • C for 45 s, and elongation at 72 • C for 1 min (35 cycles), followed by a final extension at 72 • C for 10 min. All the PCR reactions were set up along with a negative control (without template DNA). The amplification product was analyzed on 1% agarose gels by running a standard DNA molecular weight 1 kb marker parallel to it.

rDNA Sequence Analysis
The PCR amplified internal transcribed (ITS) rDNA region from each of the isolates was purified using Favorgen PCR purification kit (FAVORGEN Biotech Corp, Taiwan), as per the manufacturer's instructions. The purified ITS amplicons were sequenced by Agrigenome, Bengaluru, India and then the sequences were analyzed in GenBank (http://www.ncbi.nlm.nih. gov/) by using the Mega BLAST sequence analysis tool. The full length ITS sequences were used as queries for establishing kinship with the published sequences and search result with maximum homology and highest score were marked for further analysis. All the ITS sequences obtained from the isolates used in the present study were submitted to GenBank. ClustalW program of the Bioedit sequence alignment editor (Thompson et al., 1994;Hall, 1999) was used to align all the obtained sequences. The resulting multiple-alignment file was used for phylogenetic analyses which were performed using Molecular Evolutionary Genetics Analysis (MEGA 7.0) with Neighbour-Joining method (Kumar et al., 2016).

Morphological Identification
All the 65 S. sclerotiorum isolates collected from the major Brassica growing regions of India were observed exhibited morphological characteristics specific to Sclerotinia sclerotiorum.
The mycelia of S. sclerotiorum isolates produced abundantly and appeared white to off-white, fluffy, delicate and generally with a brownish to blackish tinge. Sclerotia produced by the fungus were often dark brown to black and appeared concentrically in some case and at the periphery in others. The sclerotia were formed in abundant at the terminal stage of the growth and vary in its shape and size viz. oval-ellipsoid, straight to curved, 1.8-2.4 × 2.9-6.2 mm in size (Table 4). Usually, the sclerotia were produced solitary in the chain of ring-shaped either at the periphery or in the center and occasionally also formed in pairs randomly. Substantial variability in culture and morphology were exhibited by the S. sclerotiorum isolates during their growth in a controlled environment (Figure 1).

Species-Specific PCR Assay
For ascertaining the originality, the S. sclerotiorum isolates was tested by using the species-specific PCR assay. For this, the species-specific primers designed and verified by Freeman et al. (2002) were used and a single 278 bp DNA fragment amplified by PCR was specific to S. sclerotiorum, which confirms the speciesspecific identity of the isolates ( Figure S10).
The morphological features among S. sclerotiorum collections from Brassica field of different agronomic regions in India (Figure 1), were used in genetic diversity analysis. The dendrogram generated (Jaccard's coefficient) have three distinct clusters at 65% of distance (Figure 2). Cluster I, with the highest degree of genetic variation, comprises 60 isolates from different regions; cluster II has four isolates, collected from distinct geographical locations or provinces, and cluster III contains one isolate, from Haryana. The robustness of the tree was estimated by bootstrap resampling, and the bootstrap values for clusters I, II, and III were very high, 100, 99, and 100%, respectively. Out of three clusters, cluster I was the largest and differentiated into seven distinct sub-clusters: IA-IG. Almost 99% of the isolates from different provinces fell under these sub-clusters. The majority of the S. sclerotiorum isolates from Rajasthan were grouped in the sub-cluster IC to IG (83%) and rest were distributed across the dendrogram. Cluster II had 4 isolates were grouped into two sub-clusters IIA and IIB. All the four isolates one in IIA and three in IIB were from four different states whereas single isolate present in cluster III was from Haryana region. Individual isolates from Jharkhand and J&K distributed in sub-cluster IIB. Interestingly out of 4 isolates in cluster II, 2 of the isolates one each from Rajasthan and Jharkhand was moderately virulent. One notable thing in the dendrogram was the distribution of the highly virulent isolates, one from each states UP, Punjab and Rajasthan had clustered into sub-clad ID whereas another highly virulent isolates one each from UP and Rajasthan falls under IA and IB, respectively.

Pathogenicity
Plants inoculated with S. sclerotiorum developed necrotic and bleached lesions 6-9 days after inoculation, and white cottony mycelial growth appeared on stem surface after successful infestation. Pathogenic variability test conducted on 4 Brassica differentials has shown that all the isolates collected were pathogenic and able to develop stem rot disease after infection. However, the aggressiveness of the isolates was observed varies contrarily on all the host differentials and thus diverse groups were identified for each character by comparing the pathologically quantitative characteristics i.e., stem lesion length and percentage of disease incidence on different host differentials. All the 65 S. sclerotiorum isolates from different geographical regions showed variations in pathogenicity in respect of stem lesion length. After artificial inoculation of S. sclerotiorum, the typical water-soaked lesion was observed on the stems of host differentials. On the basis of longest and shortest stem lesion formed by different geographical isolates of four Brassica species, following 4 groups were made: Group I: (Brassica juncea var. NRCDR-02), Group II: (B. rapa var. toria), Group III: (B. rapa var. yellow sarson), and Group IV: (B. rapa var. brown sarson). Group I, composed of 9 isolates viz., 15,[17][18][19]25,35,36, 65 produced longer stem lesion and among them ESR-36 produced the longest lesion i.e., 49.0 cm, while ESR-26, 40 and 64 isolates produced shortest stem lesion i.e., 1.0 cm (range:1.0-49.0 cm). Similarly, other 6 groups also showed pathogenic variability. Based on the stem lesion formation, out of 7 groups, ESR-36 and ESR-60 isolates were found producing longest and shortest stem lesion on the inoculated stem, respectively.
The results of pathogenicity tests conducted by using the 65 isolates of S. sclerotiorum over the 4 different Brassica differential has shown that the isolates had differences in their pathogenicity, infectivity and disease severity on tested Brassica cultivars and are presented in Table 1. The disease severity was measured in terms of Mean Disease Severity (MDS) by measuring the lesion length appeared on infected stem and that provides the basis for recording the virulence of individual isolate as low (MDS: < 3 cm), medium (MDS: < 3-10 cm) and high (MDS: > 10 cm). Thus, based on the disease severity index of the isolates S. sclerotiorum collection was distributed into three groups viz., highly virulent (5 isolates), virulent (46 isolates), and moderately virulent (14 isolates) based on the variations observed by analyzing the mean disease severity index in the treated mustard variety ( Table 5). The un-inoculated and control plants of test variety inoculated with plain agar media showed no symptoms. The Sclerotinia isolates were further re-isolated and cultured from the infected Brassica plants following Koch's postulates.

SSR Genotyping
The population dynamics study of the pathogen plays a vital role in deciphering the mechanism behind behavioral pattern as well as in understanding the distribution pattern of the pathogen in distinct geographical areas where host species envisioned frequently by the infective pathogen. In this direction, the genetic diversity of S. sclerotiorum isolates were analyzed by simple sequence repeat (SSR) based fingerprinting. All the 25 SSR primers were analyzed after optimizing the melting temperature (Tm) for polymorphism assay on collective set of 65 isolates comprising representatives from the major Brassica growing states of India included in the study: Rajasthan (31), Uttar Pradesh (16), Punjab (8), Haryana (4), Delhi (1), Bihar (1), Madhya Pradesh (1), Uttarakhand (1), Jharkhand (1), and Jammu & Kashmir (1). Among the genome-wide distributed primers, only four primers (9I, 12L, 14N, and 17Q) has generated multiple bands whereas remaining primers were produced a single polymorphic band ( Table 3). Two primers 9I and 17Q were observed produced 2 polymorphic bands of 100-500 bp size range per isolate. Subsequently other two primers 12L and 14N have produced 4 polymorphic bands within a range of 200-1,000 bp ( Figures S1-S9). The differences in genetic diversity among the S. sclerotiorum pathogen population obtained through SSR fingerprinting in this study were found reproducible. Subsequently, UPGMA dendrogram drawn based on presence and absence of the banding patterns showed substantial diversity among the isolates. Based on the topology and similarity indices of the dendrogram, the isolates were grouped into three major clades and each clade was signified by a roman numeral (group I to III; Figure 3).
The detailed analysis of the dendrogram revealed, maximum numbers of isolates grouped in clade I and it was the largest and based on the similarity indices the clade I was further separated into seven divergent sub-clades: IA-IG. Among these sub-clades, IA, IF and IG harbors nearly 60% of the S. sclerotiorum isolates from Rajasthan whereas, remaining 40% were distributed in other clades of the dendrogram. The other two major clades, clade II and III had 2 and 4 isolates and in which 1 isolate in clade II and 2 isolates in clade III were from Rajasthan. Almost all the isolates from UP were clustered in clade I under sub-clades IA, IB, IC, IE, and IG. But, isolates from Punjab has a different trend of distribution across the dendrogram. All the isolates from Haryana were clustered in sub-clades IA and IC. Individual isolates from Delhi and J&K distributed in sub-clade IA, from Bihar, Jharkhand, and MP distributed in IC, and one isolate from Uttarakhand distributed in the major Clade II. The distribution of highly virulent isolates in dendrogram was found interesting were two from UP and one isolates from each Rajasthan and Punjab had clustered into IA sub-clade whereas another highly virulent isolate from Rajasthan was clustered in clade III. The prominent genetic diversity feature in clustering was observed in sub-clade IE, IF, IG and clade II, III those included exclusively virulent or highly virulent isolates and more interestingly these clades had at least one isolate from individual state included in this study.

ITS Sequence Analysis
The fungal-specific ITS4-ITS5 universal primers pair was used in amplifying the internal transcribed spacer (ITS) region from the DNA of all the 65 S. sclerotiorum isolates collected from infected Brassica crop grown in the diverse genetic background. The amplicon lengths and purity was estimated by gel electrophoresis and found about 600 bp in size ( Figure S11). The sequence analysis result of the rRNA data by NCBI BLAST tool observed connoted the morphological variability and supported their individual existence. Further, in NCBI GenBank database the   99-100% closest match was found with S. sclerotiorum. The sequence data information for the ITS rDNA region of all the 65 S. sclerotiorum isolates obtained in this study were submitted in the GenBank database of the NCBI (GenBank Accession numbers MF408233-MF408297; Table 5).
Based on the ITS sequence information a dendrogram was constructed by using the ClustalW and Mega7.0 software with the Neighbor-Joining method (Figure 4). The resulted phylogenetic tree was clustered in total 11 clades and six of them were formed major clades, four clades with two members each and one clade with a single member that clearly indicates the distribution of the isolates across dendrogram irrespective of their properties of pathogenicity or geographic origin. Isolates from Rajasthan and UP were found distributed in nine and six different clades, respectively, whereas isolates from Punjab were found distributed in five different clades. Clade I was found to be the largest and consisted of 21 divergent isolates from Rajasthan, UP, Punjab, Bihar, and Delhi with varied pathogenicity (highly virulent, virulent and moderately virulent). Clade II had virulent and moderately virulent isolates from Rajasthan and UP whereas isolates from Haryana were all virulent and distributed in clades IV, V and VI. Clade VII had only two isolates and both of them were from Rajasthan and virulent. Interestingly, the clade IV and V were also found heterogeneous as of clade I in terms of the isolates of different region and their virulence features. The isolate from Jharkhand was found distinctly in the clade VIII and was moderately virulent. Remarkably clade X had four isolates from the four different states. Two pathogenic isolate each from Rajasthan and Punjab were formed clade III, from UP and Punjab have formed clade IX and from UP and Rajasthan were formed clade X. One moderately pathogenic isolate from Jharkhand formed the single member clade VIII. The highly pathogenic isolates from UP, Rajasthan, and Punjab were distributed in Clades I, V, and VI.
The ITS sequence analysis of all the isolates in the study depicted on an average 95-100% of sequence similarity that entails conspecific nature of the isolates and their host from which the isolates were collected irrespective of different geographic regions.

DISCUSSION
The incidence of stem rot disease in India, especially in oilseed Brassica, has been known to damage the crops since from the mid-1990 but the genetic diversity of the pathogen especially in oilseed Brassica has not been comprehensively studied so far which is primarily required for understanding the mechanism of the pathogenesis and virulence of the pathogen. From the last decade, the stem rot disease in Brassica poses a serious threat to the crop every year in India and causes significant damage to the crop in terms of yield and economic losses to the farmers due to lack of effective management practices and other control measures. From the frequent disease occurring reports and extent of its severity over the oilseed crops, this pathogen has drawn attention and interest of the researchers to analyze the host-pathogen interaction in a broader spectrum as the mustard growing farmers experiencing this problem almost every years during the cropping season. The present investigation thoroughly reporting the existence of S. sclerotiorum pathogen in mustard fields and the genetic diversity among the population distributed throughout the Brassica growing regions of India.
The field-based crop-specific disease survey in different states of the India has shown the high occurrence and well-dispersed existence of the stem rot disease that commonly associated with the scorched lesion on the stem of the infected plants and sclerotia infested soil samples in the severely affected crop field. The regular appearance of the stem rot disease in Brassica indicates it is a recurrent menace for the crop and the lack of promising resistance source to the pathogen in existing cultivated varieties. Although the isolates in the present investigation were collected from Brassica field only but even though they were found genetically diverse in terms of their culture conditions and in morphological appearance as well. This might be due to the divergent geographical location from where the isolates have been collected. The evaluation of the pathogenicity of different Brassica cultivars has revealed that there was differential susceptibility in host species against different isolates of the S. sclerotiorum and thus based on the pathogenicity assay and virulence properties the pathogen population was grouped into three subgroups highly virulent, virulent and moderately virulent isolates. Depending on their virulence all the isolates resulted in the development of lesions of white patches of a different length over the infected stem of the Brassica cultivars that indicates the existence of genetically diverse and variable population of S. sclerotiorum. The pathogenicity assay and virulence variability test have shown that out of 65 S. sclerotiorum isolates 5 of them were highly virulent, 46 were virulent and 14 isolates were falls under moderately virulent category and they disseminated through either contaminated soil or seed across Brassica growing regions in India. The genetic diversity analysis by all the means followed in this study has clearly shown there were not many differences exists in terms of virulence and in groups that comprised from the isolates of different virulence nature belongs to different localities. Our results are in accordance with the observations made in previous studies where aggressiveness of the isolates are being independent of their geographic distribution (Atallah et al., 2004;Xie et al., 2014). Hence, the nature of the virulence of the S. sclerotiorum pathogen possesses a potential threat to the planting of Brassica cultivars susceptible to the stem rot disease.
Sclerotinia sclerotiorum apart from reproducing asexually by myceliogenic germination of sclerotia it is also able to produce sexually through apothecia and ascospores formation. The diversity analysis and genetic architecture of S. sclerotiorum have shown the existence of homogeneity across the isolates with little variation in 18 and 28 s rDNA regions (Kohli et al., 1995). The genetic diversity of the pathogen, severity of the disease and frequency of its occurrence are the important factors that directly correlates with the epidemics of the disease at the regional level and hence their comprehensive studies are the foremost requirement for understanding the behavioral pattern of the pathogens for effectively design the management practices for the control of this disease more efficiently. The polymorphism assay with the genome-wide distributed SSR markers has been proven effective in studying the existence of genetic diversity within the pathogen population. The SSR is the co-dominant molecular markers and it is a PCR-based robust technique that amplifies the putative repeat regions of the DNA and therefore it is able to identify a higher level of polymorphism than PCRbased another less robust markers like RAPD and hence it has been extensively used for analyzing the fungal population. The genetic diversity among S. sclerotiorum isolates based on the morphological features and SSR fingerprinting had observed quite diverse and reveals that some of the isolates grouped together belonging to the same region while others from the same region distributed across the dendrogram. The morphological features of the isolates and their SSR fingerprinting confirms the existence of a highly variable population of S. sclerotiorum even within the same crop and the region. However, among some isolates of S. sclerotiorum those belonging to the same geographic regions has been found having lower genetic diversity and that might be due to the locally adapted cropping system where introduction of newer varieties are usually not preferred over the existing one and hence diversity of the pathogen remained to stagnate due to lack of introduction of the newer genotypes of the pathogen in the area under cultivation.
The result of the present investigation depicting the genetic diversity of S. sclerotiorum population based on the morphological parameters and SSR analysis has been observed in coherence with earlier findings of the studies made on the same pathogen but in different host species. In eggplant, the genetic diversity analysis of S. sclerotiorum population by using the SSR and RAPD molecular markers, Tok et al. (2016) had shown the existence of high level of heterogeneity in the pathogen populations and presence of different isolates of S. sclerotiorum pathogen in the same regions. The presence of high heterogeneity in genetic composition of the S. sclerotiorum isolates in different regions of Turkey has reflected by their distribution across the dendrogram in different groups made from the RAPD and SSR profiling. In a similar study, Barari et al. (2013) have characterized the genetic diversity within the population of S. sclerotiorum isolates from Iran and the variations among them were projected by SSR profiling. The depiction of variability in closely related isolates by SSR markers demonstrated its effectiveness in analyzing the genetic diversity among S. sclerotiorum pathogen. In chickpea, the use of RAPD, ITS-RFLP, ITS sequencing, and mycelial compatibility groupings (MCG) has shown the limited variability among the distinct S. sclerotiorum isolates possess the higher genetic homogeneity and partially correlated with their geographical origin (Mandal and Dubey, 2012). The robustness and reproducibility features of the SSR markers make it more effective and reliable marker for deciphering the genetic relationship among population studies and thus the phylogenetic analysis of S. sclerotiorum based on the SSR genotyping revealed the existence of wider genetic diversity among its population. The present findings are in connotation with the earlier reports that molecular typing by SSR is an efficient technique for establishing the genetic relationship among the diverse population and it is being preferred over other molecular tools for diversity analysis.
Further, the ITS (internal transcribed spacer) sequence based molecular phylogeny was established by using the ITS region as a genetic marker to investigate the genetic association within the S. sclerotiorum isolates. This phylogenetic study has facilitated the revelation of the evolutionary relationship within S. sclerotiorum species collected from a single host Brassica species from different geographical regions of India. Based on the ITS phylogeny all the 65 isolates investigated in the present study were observed grouped into six major and five minor evolutionary lineages of S. sclerotiorum isolates collected from different mustard growing states of India. In all the lineages deviation in virulence was not so evident and this suggests that the isolates from single host species may not differ ominously in their pathogenicity. The least variation in virulence of the pathogen has been observed in the case where isolates from defined geographical regions were taken into the investigation (Alvarez and Molina, 2000;Atallah et al., 2004;Sexton and Howlett, 2004). The differences in virulence can be more evident and conclusive when isolates from the widely distant geographical area will be taken for comparison. Kull et al. (2004) has found the conclusive evidence for the host specialization among S. sclerotiorum isolates.
Although, there is no report of any toxin is usually produced by the S. sclerotiorum pathogen it has been observed the production of oxalic acid (Cessna et al., 2000) and extracellular lytic enzymes production (Riou et al., 1991) invariably by all the isolates during pathogenesis. However, extensive biochemical and physiological studies are required to investigate the production of any natural mycotoxins by the pathogen if any like other fungus and further needs to find the effect of those toxins on the host plants and in herbivores consuming those plants. In this study, we attempted to verify the authenticity of S. sclerotiorum isolates by using the Sclerotinia-specific PCR primers and found all the geographical isolates (65) investigated were genetically diverse.
Development of the Genetic fingerprinting database for distinguishing the pathogen diversity between and within the fungal species is one of the foremost requirement for determining the origin and evaluation of the pathogen. Unequivocally, environmental factors are also influencing the properties of the pathogen and play crucial roles in determining the genetic differentiation within the pathogen population (Weller, 1988). The existence of a mixed population of S. sclerotiorum within a region could be due to the possible human interventions like the movement of infested soil through seeds, packaging materials or agricultural equipment that might carry the pathogen from one region to another. Results of the present study, as well as evidence from the previous studies, indicates the minute changes in physiological properties of the pathogen like changes in virulence due to evolutionary changes happening in due course of time that brings changes among genetically distinct populations of S. sclerotiorum isolates. Such changes may happen because of the varying force of selection over pathogen imposed by the selected set of genes of the host germplasm being tolerant to the invasive pathogen. The genes involved in host defense system and coding for the virulence factors in pathogen have been observed comparatively much more vulnerable for accumulating the mutation under positive selection pressure than the genes those regulating the basic physiological processes in the organism (Matute et al., 2008). As a consequence, in most of the hostpathogen interaction pathogen often adapts themselves to the evolutionary changes arises in response to the variation in the host genotype.

CONCLUSION
The current investigation reports the studies on morphological features, pathogenicity variance over the Brassica cultivars, and genetic diversity based on molecular characteristics of the S. sclerotiorum isolates population collected from 65 locations of the 10 major Brassica growing states of India. This is the first comprehensively studied report on S. sclerotiorum isolates especially from Brassica field of India that indicates the existence of broader genetic variations among the pathogen and this would facilitate the plant breeders to use the information in Brassica breeding program for stem rot disease resistance development. The SSR-based markers utilized in present investigation for diversity analysis can be effectively used for identification and analysis of S. sclerotiorum population from various host species. The present study offers the basis for increasing the basic understanding of host-pathogen interaction and also for deciphering the genes and associated pathways involved in combating the pathogen invasion. The evidence of the present study may provide a broader understanding of the pathogen virulence and that will facilitate the development of the effective disease management strategies based on molecular breeding and other advanced approaches.

AUTHOR CONTRIBUTIONS
PS, NG, and LP conceived and obtained funding from ICAR-EMR project for Sclerotinia work in Indian mustard. PS and VS collected isolates from diseased plants and along with LP done the morphological characterization and pathogenicity evaluation. AS and MR obtained isolates from PS, carried out microsatellite genotyping. NG and RB do the ITS sequencing of S. sclerotiorum isolates and wrote the paper. NG and DM carried out clustering and phylogenetic analyses of microsatellite data and ITS sequences. RB and VS edited the manuscript.