Molecular insights on differential responses of antioxidant and pathogenesis-related genes in pigeonpea [ Cajanus cajan (L.) Millsp.] to Fusarium wilt ( Fusarium udum (L.) CURRENT STATUS: POSTED

Background: Pigeonpea [Cajanus cajan (L.) Millsp.] is a major dietary protein source to a large vegetarian population in the world. However, Fusarium wilt (FW) caused by Fusarium udum Butler (FU) is one of the challenging factors to pigeonpea production. FW resistance in pigeonpea is restricted to specific agro-climatic zones. Therefore, exploring a superior pigeonpea genotype from landraces or local cultivars cultivated in various rural regions and validating its resistance response at biological and molecular level could be an alternative to crop improvement program. The present study was carried out to identify superior genotypes of pigeonpea through selection of innate resistance to FW using different biological and molecular approaches. Result: Five distinct genotypes of pigeonpea along with well-known one resistant (ICP2894) and susceptible (ICP2376) control were selected on the basis of percentage of FW disease incidence from three different states of India after an extensive survey work. Among them, Richa had shown very less FW incidence (10.0%) under genotype evaluation study and was further selected for experimental validation at molecular level for its innate resistance to FW. Molecular characterization and expression profiling of biotic stress responsive genes representing antioxidant enzymes (APX and SOD) and Pathogenesis related proteins (CHS and β-1,3-glucanase)families were carried out. Under FU challenged conditions in Richa, expression patterns of both AO genes found to be similar, while expression of both PR genes was gradually increased with inoculation time. Putative 3D protein structure of CHS and SOD were generated. Molecular docking with fungal effect or protein cellobio-hydrolase-c (CBH-c) also revealed that their role in blocking the hydrolytic activity of CBH-c during pathogen invasion. Conclusion: The obtained results of expression study and in-silico analysis showed that the higher level of induction of PR and AO genes and strong interaction of their putative proteins with fungal CBH-cprotein established their antifungal activity conferring early plant defense responses to restrict the pathogen growth in resistant genotype Richa. Our study displaying a strong and combinatorial approach involving biological assay, molecular experiment and in-silico analysis to identify a superior genotype of pigeonpea resistant to FW across a major biogeographic region. investigate expression through quantitative PCR (qPCR) present report envisages an elaborative biological and molecular study to investigate the genetics of tolerance mechanism of some well-known defense responsive and antioxidant genes in pigeonpea against FW.


Background
The world population is continuously increasing and aggravate by 1 billion in 2050 to an estimated number of 9.6 billion from the present population of 7.2 billion (UNPAN, 2010) resulting in more requirement for food grains. With a long term planning for sustainable genetic improvement of major staple crops like rice, wheat and maize, it is necessary, at the same, time to gain a breakthrough in uninterrupted production of proteinaceous foods to strike a balance for reducing global hunger and malnutrition. Legumes, preferably pigeonpea, can be considered as a good alternative to easily available protein resource and offers a handful quantity of food proteins in Indian and African subcontinent with less cultivation care and low inputs. Pigeonpea belongs to family Fabaceae and is contributing by fulfilling the protein demand of 20% of global population and also serves as a major source of vitamins. In India, pigeonpea stands second after chickpea as the most essential food legume crop [1]. But now days, a major restraint in pigeonpea production is caused by biotic stresses wherein Fusarium wilt (FW) is considered as the most devastating disease followed by Sterility mosaic disease (SMD) and Phytophthora blight respectively. The former accounts for 30-100% yield losses in pigeonpea genotypes vulnerable to this disease [2].
Assessment of the effect of any disease on yield of a crop with the help of a complete survey is prerequisite for preparing rational disease management program. The system biology approach encompasses that the impact of biotic stress especially pathogen attack in combination with other variable abiotic stress factors leads to the molecular fine tuning of the degree of resistance, tolerance or susceptibility on a specific host plant [3]. On the other hand, sources of resistance identified in one region do not necessarily perform with the same degree of resistance in other regions, thereby indicating pathogenic variability in the fungus [4]. Several new sources of resistance have been reported against FW but there are still a lot of probabilities of getting improved and potential germplasms or indigenous genotypes at regional/rural level by searching, collecting and evaluating such genotypes for better resistance by following standard inoculation method [5]. The present investigation of germplasm evaluation against FW also revealed the importance of studying epidemiology of wilt disease. The genetics of FW resistance and or tolerance in pigeonpea is still in an ambiguous mean, and many more genes are supposed to be involved from a single dominant gene to two complementary genes and there might even be involvement of multiple genetic factors [6]. The in-depth knowledge on genetics of resistance mechanisms in host plants and genomic insights of wilt causing pathogen Fusarium udum are equally crucial for the genesis of effective strategies for efficient transfer and stable function of such resistant genes into disease susceptible cultivars. A number of resistance sources of germplasm against FW were identified by screening various pigeonpea genotypes over the last two decades at national and international level. Some of these promising genotypes were effectively used in crop improvement programs [5][6][7][8].
The fabrication of defense molecules such as phytoalexins, phenylpropanoids; pathogenesis-related proteins (PRs) viz.chalcone synthases (CHS), chitinases, and β-1,3-glucanases have very imperative roles because of their direct attacking activity on fungal and insect structural components. Beside these, a large number of biotic stresses related cellular enzymes such as superoxide dismutase (SOD), catalases, peroxidases are over-expressed to get rid of the pathogenic entry [9]. Two PR protein associated enzymes, namely chitinase and β-1,3-glucanase were successfully extracted and purified by gel filtration from the fungal pathogen induced chickpea plant and their higher activity was further confirmed by inhibiting the growth of Fusarium oxysporum f. sp. Ciceris (Foc) and other phytopathogenic fungi [10]. The defense network response in chickpea plant against Foc was also established by highlighting the studies on confocal microscopy to trace pathogen invasion, various enzymatic assays to detect structural changes in plant cellular membrane and validation of defense responsive transcripts by qRT-PCR approach [11]. Microarray based genome wide transcriptional study in F. oxysporum f. sp. lycopersici (Fol) and Tomato mosaic virus (ToMV) induced Tomato plant revealed a strong homeostatic defense response where chitinase and other PR gene families displayed high-level positive regulation in terms of fold changes [12]. Molecular cloning and characterization of β-1,3-glucanase gene in wheat inoculated with the stripe rust pathogen demonstrated a high level of expression pattern of this gene under qPCR analysis [13].
Considering the economic importance of pigeonpea and severity of FW disease, the present study was conducted to explore the existence of host resistance in case of highly wilt tolerant pigeonpea genotypes. A preliminary screening of pigeonpea germplasms collected from different states of India against a highly virulent FU strain was performed to identify resistant/tolerant genotypes. Attempts were made to identify the defense-related genes in tolerant genotype by implying molecular and insilico approaches and thereby investigate their expression pattern through quantitative PCR (qPCR) during pathogen infection. The present report envisages an elaborative biological and molecular study to investigate the genetics of tolerance mechanism of some well-known defense responsive and antioxidant genes in pigeonpea against FW.

Survey and epidemiological study
During the extensive survey in the period between October and January 2014-2015, 25 different villages/Tehsils from four major pigeonpea growing state of India viz. TN, MP, AP and WB were successfully covered to carry out epidemiological study and searching of new pigeonpea sources tolerant to FW from existing germplasms indigenously grown in the surveyed areas. The maximum number of areas was covered in MP because there is no such survey report found in this state, which is contributing around 10% of total pigeonpea production in India (Table 1). As per DI %, the lowest DI % (4.5%) was observed in Perayambetti gate village of Dharmapuri district in TN, whereas the highest DI % (26.0%) was observed in Anjad village in Barwani district of MP, followed by Sanawad village of Khargaon district of same state ( Figure 2). As stated earlier, the surveyed areas were maximally covered in MP where three districts (Dhar, Khargon and Barwani) displayed the highest level of epidemiological variation in terms of DI % ( Table 1). As per the overall disease incidence is concerned, all the four states had a wide range of variations in terms of disease epidemiology. The seed samples of pigeonpea were exclusively collected from Warrangal and Medak regions of AP where pigeonpea is cultivated abundantly. Some local genotypes cultivated in the high and less disease prone areas of MP had natural characteristics of showing extreme tolerance and susceptibility to FW respectively. Maximum types of seed genotypes were collected from major pigeonpea growing regions of MP viz. Dharampuri, Manawar and Kukshi villages of Dhar district. The pigeonpea seeds of a total number of eight unique genotypes were collected from three different states viz. MP, AP and WB ( Table 2). The selective collection of seeds of different genotypes was truly accomplished based on epidemiological data where genotypes were mostly cultivated with a higher and lower level of disease severity.

Stress treatment and screening for tolerance
Out of seven genotypes used for resistance screening, only two genotypes, ICP8863 and Richa were found to be promising showing very low wilt incidence. In known FW resistant genotype ICP8863, only 3.33% disease incidence was recorded. In genotype Richa, the disease incidence was 10% and was thus put into highly resistant category, whereas another genotype, Desi Nimar came under moderately tolerant genotype because of 23.3% DI. All the remaining four genotypes, Parwati, Desi Tur, WB-20/105 and ICP2376 were found as highly susceptible to FW (Table 2).

In-silico study
For in-silico analysis of biotic stress defence-responsive genes, a total of twelve diversified genes were initially screened as all of them are directly or indirectly involved in defense response in plants where some are falling under a specific class of functional gene family. Prior to multiple sequence alignment, all the selected sequences of a single gene were graphically aligned to understand the level of divergence and convergence within ortholog sequences. In case of CHS, the CDS region covered from 58-1227 nucleotides with a length of 1170 bases covering complete CDS sequence. In a similar way, full length CDS coverage in maximum orthologs of a given gene was sharply observed in case of isoflavone reductase (IFR), methionine sulphoxide reductase B2 (MSR), glycosyltransferase (GTF), chitinase, β-1,3-glucanase and osmotin genes. In another case, protease and ascorbate peroxidase (APX) genes have the full length CDS for single gene sequence, but major orthologs did not cover this region of target gene due to their short size (Table 3). Glutathione peroxidase (GPX) gene was found to have a lowest match in the CDS region (120 ntd.) from full length CDS of 648 bases. Further, superoxide dismutase (SOD) gene consisted the lowest CDS (649 ntd.) but only conserved 500 bases (1-500 ntd.) in ortholog sequence. Beside the graphical alignment of all gene accessions, multiple alignments were also performed using all mentioned orthologs of each gene to revalidate the conservation in CDS region. For above twelve gene sequences, the maximum conservation in terms of length was visually observed in the case of nine genes except osmotin, protease and IFR. The phylogenetic tree was constructed for selected nine gene sequences by taking only the full-length CDS region of each gene and partial CDS where the full-length sequence was not available at all. From the cladogram analysis of all selected phylogenetic trees, maximum number of accessions grouped with a target sequence in a single cluster were observed in the case of APX, SOD (as an antioxidant enzyme family) and CHS and β-1,3-glucanase (PR protein family) with a defined biological and molecular defense response against biotic and abiotic stresses. Therefore, these four genes were finally selected for further molecular experiments.

Reverse transcription polymerase chain reaction (RT-PCR)
After several optimizations of gradient scale in PCR for each transcript, the best annealing temperatures found to be reproducible for successful amplification were 52.6 o C for CHS, 60.0 o C for APX, 57.0 o C for β1,3 glucan and SOD. The desirable amplicons with expected size (CHS≈ 560bp, SOD≈ 400bp, APX≈ 350bp and β1,3 glucan≈ 450bp) were obtained respectively followed by analyses on 1.0% agarose/EtBr gel to ensure the successful molecular identification of these four defense responsive genes (Supplementary Figure 1A-D). The identification of these four genes was performed truly upto RT-PCR level as all the amplicons were derived from cDNA rather than genomic DNA showing that these are specific, intensified and unique bands of desired length. Therefore, our study ensured the presence of gene transcripts at mRNA level and was sufficient to authenticate their expression pattern by using real time PCR approach.

Validation of FW stress responsive genes by qRT-PCR in pigeonpea
In this study, expression profiling of four defense-related genes in susceptible (ICP 2376) and resistant (Richa) genotypes of pigeonpea during FW development was determined through qPCR. In the comparative analysis of antioxidant gene expression in pigeonpea during FU infection, it was observed that both the genes, APX and SOD were up regulated in inoculated plants irrespective to genotypes. In overall expression study, it was observed that the expression of the SOD gene was relatively higher than APX gene in later stage of infection in susceptible genotype ICP 2376. Although the expression levels were varied in different time points at 7DAI and 15DAI, the expression was differentially regulated by a combination of positive (ICP 2376) and negative induction (Richa) ( Figure   3A). It can be hypothesized that the degree of antioxidant activity needs to be enhanced with the time of infection and it delayed to make a balance in hypersensitivity response in case of susceptible pigeonpea genotypes to Fusarium. In case of genotype Richa, the difference in APX gene induction (F.C-2.42 to 1.73) in two different TP could be due to activation of anti-oxidative response immediately after the fungal attack and thereby decreased gradually as the fungus was restricted by other defense mechanisms of that genotype. After immediate enzymatic activation, SOD gene product might have accumulated in sufficient amounts to necessitate the switching off of the gene.
The effect of FW pathogenicity stress on SOD gene expression in both genotypes was quite interesting which has shown a variable range of positive induction with a gradual increase with inoculation time in case of Richa whereas the trend for ICP2376 cv. was reverse, showing a slight decrease from 2.80 to 2.41 as TP extended up to 15 days ( Figure 3A). Thus, it was confirmed that both genes of same functional family were induced immediately after the fungal attack, but the expression needed to extend as time of infection increased to sustain the susceptible genotype with maximum defense response whereas the resistant one can recover from the infection after some time and the APX gene came to its native state.
In case of comparative study of expression pattern between CHS and β-1,3glucanase genes, these PR genes had a mix response in terms of induction in both susceptible and the new genotype with unknown HR functions. Initially, after increment of inoculation time, the CHS gene in Richa was increasing its expression up to 5.03 fold from an initial level of 2.21 fold compared to the nearly stagnant expression in ICP 2376 (FC 1.79-2.09) confirming that the said enzyme gets involved in stimulating the initiation of defense response at ground level ( Figure 3B). On the other hand, as the inoculation time increased from 7DAI to 15DAI, the relative expression of β-1, 3 glucanase gene in ICP2376 is more or less constant with an apparent fold change (FC) of 2.40 to 2.89 ( Figure 3B). In case of Richa, a continuous change in the inoculation time further accelerated gene activity by gaining induction up to 15.61 fold from FC-7.47. The overall expression pattern of this gene in both genotypes had shown zigzag type model and it could be possibly explained by the fact that most gene expression is dynamic in nature which is not expected to be up-regulated in all points of factor.

In-silico analysis and protein modelling
To validate the transcript abundance in Fusarium inoculated pigeonpea genotype (Richa), the PCR purified products of partially amplified CHS (537bp) and SOD (361bp) genes were sequenced. The  Table 1) and both proteins were classified as stable protein based on its instability index (CHS: 28.02 and SOD: 5.66). SOD is considered as highly stable structure due to lowest instability index. Further Ramachandran plots were studied where phi(φ)/psi(ψ) angles of the amino acid residues were 96.89% and 95.19% for CHS and SOD respectively (Figure 4d, 4e) thereby proving the stability of both proteins. The QMEAN value was analyzed for the quality of 3D structure. The QMEAN scores for CHS and SOD are 0.84 and 0.85 respectively. The reliable 3-dimension structure lies between the range of 0.5-1. The structure of fungal CBH-C was also checked for stability. The instability score (39.12) and QMEAN value (0.61) confirms its reliability. The docking was done using FRODOCK (Figure 5a, 5b) and simulation was carried out using online server GRAMMX.
Clashes/contact score was calculated for both the interaction models (Figure 5c, 5d). The maximum correlation score was obtained for CHS i.e, 189394.84 and for SOD it was 1.59. CHS docking with CBHc has the clash score of 20 and 138 contacts whereas SOD interacts with CBH-c with 20 clashes and 179 contacts. CBH-C bind more strongly with CHS as compared to SOD. In relation to the attack of FU, CBH gene present in fungal cell can be considered as a potent factor for the synthesis of cellobiohydrolase and this cellular enzyme is secreted by fungi to degrade host plant cellulose through the discharge of cellobiose for the deployment of their pathogenic activities [27]. The present interaction study of both plant genes with FU-CBH gene could be correlated with the qPCR analysis where higher rate in fold change expression of CHS (7.47 at 7DAI to 15.61 at 15DAI) than SOD gene (2.41 at 7DAI to 2.61 at 15DAI) was observed. This can be concluded that the CHS gene transcript is more hyperactive than SOD with a tendency to release more number of functional proteins through translation process after pathogen inoculation and it could further functionally block the activity of CBH protein by making a strong binding affinity.

Discussions
Fusarium wilt is a vital ailment of pigeonpea and is realized to be the worst problem for cultivation of this crop in India [28]. The annual economic losses due to wilt in pigeonpea have been assessed at US $ 36 million in India and the subsequent yield loss can also be increased up to 100% in susceptible cultivars [29]. So, it is necessary to find a way out in sustainable manner to prevent the disease in this crop. Among the different control measures, early diagnosis and control in the light of advancement of molecular technologies are quite promising. The importance of wide level of survey is not only related to epidemiology of a disease to understand, describe, compare and predict epidemic, but also make a chance to explore the most promising highly tolerant genotype/s against a particular disease from surveyed areas. From 1980 to till date, a number of surveys were carried out in different pigeonpea growing countries. The present survey-based study reported regional diversity of more than 25 pigeonpea-cultivated areas representing the major agro-climatic regions of India under pigeonpea production, which were found to be suitable for growing local genotypes/cultivars due to their potential nature of tolerance/susceptibility mechanism to FW. Taking into consideration of the most variable geographical region, MP represented a great diversification in disease incidence range from resistant to highly susceptible areas. As far as lowest disease epidemiology is concerned, TN was found to be the lowest disease prone state with no class as highly infested area and very highly infested area. Besides India, this disease was reported from East Africa and Malawi, where yield losses were reported to have crossed above 50% in some places. Yield losses in pigeonpea due to FW was also reported form countries like Bangladesh, Indonesia, Grenada, Myanmar, Mauritius, Nevis, Nepal, Venezuela, Trinidad and Tobago [30] This pathogen was also recorded for causing severe disease in the Southern Zambezia province under South Africa [31].
The study of exploration of pigeonpea genotypes to FW displayed an overall contrasting response for all selected genotypes. The experimental reference variety ICP-8863 displayed a negligible or near to zero level of symptoms (D.I-3.33%) in one out of 30 plants as it is already reported for its best resistance mechanism against FW independent to environmental and geographical conditions. It is finally concluded that Richa was one of the most highly tolerant genotypes (%D. I-6.67) after ICP-8863 with a higher degree of tolerance potential and that could be due to specific morphological modification in the xylem vessel of the root system or targeted activation of some defense related proteins and transcription factors during fungal stress or due to both reasons. On other hand, Parwati was authenticated as one of the most susceptible among all experimental genotypes and it could be due to its weak defense system. Although ICP 8863 and ICP2376 were confirmed for their extreme resistance and susceptible characteristics, respectively in our experiment and were earlier validated so many times for these features [28,32]. Therefore, these three genotypes viz. ICP 8863, ICP2376 and Richa were picked for further molecular biology studies. In overall evaluation data, all the four classes of disease scales were fitted for selected genotype population. A number of efforts have been taken to find out the most potential genotypes or cultivars resistant against FW over time and pathogenic strains. First time observations on wilt resistance in pigeonpea was reported by Butler in 1908 [33]. Apart from this, 950 genotypes screened had shown zero resistance to FW with less than 10% wilt incidence in 19 genotypes [29], 16 out of 31 pigeonpea cultivars screened were resistant to wilt with highest degree of resistance (2.15 %) in BWR 369 cultivar [34]; six cultivars were reported to be resistant to FW among 216 late maturing pigeonpea germplasm evaluated [7]; resistance screening in 976 genotypes, germplasm and breeding lines FW using wilt sick plot confirmed 92 genotypes resistant after a rigorous two years testing at Patancheru, India [35]. Another study on evaluation of new elite pigeonpea germplasm using wilt sick plots demonstrated a consistent rate of high level of resistance (DI < 20.0 %) in ICEAP-00040 genotype to FW reported in Kenya, Malawi and Tanzania [31].
For the identification, quantification and expression profiling of highly active plant defense responsive genes, proteins under enzyme group and in most cases, transcription factors and high throughput molecular techniques have been developed and customized day by day for their successful implication. In our current study, understanding of spatio-temporal expression of major defense responsive promising genes during pathogenesis and the study of different patterns of symptomatology based tolerance in various plant genotypes specific to various agro-climatic regions of India can provide the vital knowledge for crop improvement with well-developed tolerance mechanism. Higher plants directly defend a number of biotic and abiotic stresses with their broad range defensive mechanism, especially the production and accumulation of pathogenesis related proteins in plants in response to invading pathogen [36]. A wheat β-1,3-glucanase gene (TaGlu) induced by Puccinia striiformis (Pst) was cloned, characterized and qPCR analyses confirmed that the transcription of TaGlu was induced in both compatible and incompatible interactions [13]. In our present molecular biology study, the modified RNA extraction method is equally useful to extract stressed sample as well as fungal RNA with maximum reproducibility to trace their gene's functions and behaviour under various stresses imposed during pathogen attack [16]. By using isolated RNA from modified method, RT-PCR based detection of four genes belonging to AO enzyme and PR protein groups was performed and their successful amplifications made an inference that the genes are actively present in plant system during FW stress. The corresponding mRNAs were synthesized at the nucleus and circulated over cytoplasm followed by their translation to the consequent proteins or enzymes, which are directly or indirectly involved in defense response to significantly enhance their tolerance to pathogen attack. To precisely monitor the expression level of a gene at various stages, qPCR-based transcript profiling approach is highly adorable and the present comparative analysis of expression of AO genes with PR protein genes has revealed a variable mode of expression between susceptible and biologically validated resistant cultivars. It was also indicated that no significant differential induction has been located in case of SOD for both genotypes and CHS for ICP2376. It can be concluded that these genes did not have any significant expression variation under FW stress by considering the varietal difference. Higher magnitude of transcript expression at the last TP confirms that the expression was specific to infection stage rather than gene and genotype. This may help the plant to cope with the deleterious effects of pathogen attack at early stage critical for later growth and development [37]. Comparatively, the imperative variations were observed in the case of APX and β-1, 3 glucanase genes. An attractive finding of the present study is the differential expression of

Conclusion
The evaluation of the surveyed results revealed the feasibility of pre-screening of such genotypes falling under the first group (DI 0-10%) for genetic diversity analysis and their potential association with FW tolerance. Both well-known pigeonpea varieties (ICP8863 and ICP2376) showed their similar mode of sensitivity to FU as reported earlier, where two genotypes from MP showed tolerance ability and two were very susceptible to FU stress. Therefore, Richa is considered here as the best performing highly tolerant genotype to FW based on evaluation of the overall survey and resistant screening results. Under molecular study, our aim was to identify promising genes and robust alleles for fungal disease resistance which can be further utilized for the development of crops challenging to FW infection. Therefore, we performed the molecular detection of four defense-related genes through RT-PCR based approach and thereby quantifying their expression level by using qPCR technology. The β-1, 3-glucanase under PR protein family is highly functional and CHS gene under polyketide synthase super family is moderately functional intolerant genotype rather than susceptible one. In another case, APX and SOD under AO group have also shown positive regulation to trigger and modulate resistance response with less variation in two extreme genotypes. The identification of such genes would be potentially used as molecular markers for screening resistance genes from pigeonpea as well as could be used for developing transgenic and/or improved pigeonpea varieties resistant to FW infection suitable under changing global climate.  (Table-1). To understand the epidemiological intensity of the disease in each village under a particular surveyed district, the grading of the % DI was further set into five categories (Figure-2).

Collection of seed materials
For exploring new pigeonpea genotypes/landraces resistant to FW and grown locally, during disease survey a parallel field survey was undertaken for a collection of seeds from apparently resistant genotypes/landraces in such areas (villages) where Fusarium wilt disease was observed naturally (Table-2). Seeds from a total of eight different pigeonpea genotypes/ landraces were collected.
Among them six were exclusively the landraces or local cultivars cultivated by the farmer communities year after year. These landraces were not even exposed and used by the scientific communities or plant breeders to develop superior variety or to conserve their germplasm in regional or national institutes with a specific accession number. After getting permission and confirmation about their local name from the respective farmer communities, the seeds of a specific genotype were exclusively collected from a single plant which was selected based on growth, health and yield.
The crop was identified informally by the main author and respective farmers', as the crop is well known and widely cultivated species.

Plant growth conditions
The apparently bold and healthy seeds of all collected pigeonpea cultivars/genotypes were surface sterilized using 0.1 % HgCl 2 followed by three times washing with sterile water. The treated seeds were then sown in polythene bags filled with sterilized river sand and watered properly. The polythene bags were then kept cautiously in glasshouse maintained at 28-30 ±2°C with natural photoperiod conditions. The 7 days old seedlings were then used for pathogenicity test and other molecular studies.

Inoculum preparation and pathogenicity test for germplasm screening
The highly virulent strain of FU (ICFU 109) was initially inoculated in potato dextrose agar (PDA) media from slant culture. After seven days of incubation in PDA plate, the young FU culture with profuse conidia was further inoculated in PDB media. The FU inoculated broth culture was kept in an incubator-shaker maintained at 25°C with a continuous shaking of 120 rpm for 6-7 days. For screening of resistance and or tolerance level in pigeonpea genotypes against FU under controlled conditions, the root dip screening technique was used [4]. Conidial suspension of FU was diluted with sterile water to maintain a final concentration of approx 6 × 10 5 spores/ ml. Seven days old, seedlings of each genotype were carefully uprooted from germination trays and washed under running water.
Root tips (0.5 cm long) of each seedling was cut off and dipped in the diluted inoculum suspension for 1-2 minutes for pathogen entry. A similar procedure was followed for experimental control where injured seedlings were treated with sterile distilled water. The inoculated seedlings were transplanted into six-inch pots filled with sterilized soil, sand and farm yard manure (FYM) in a ratio of 1:1:1. Each pot contained three plants and each genotype was transplanted into ten pots (three seedlings/pot). Based on disease incidence, genotypes were categorized as tolerant/resistant (≤10.0%); moderately tolerant/resistant (>10.0 -20.0 %); susceptible (>20.0-40.0%) and highly susceptible (>40%). At each TP the inoculated seedlings were removed from the soil, washed with sterile water, completely dried with tissue paper and quickly snap frozen in Liquid N 2 followed by storing at -80°C.

Molecular characterization of defenseresponsive genes e against FW in pigeonpea
In-silico study for targeting stress responsive genes The aim of the present study is to find out some important and promising genes, which have a major biological role directly or indirectly related to defense response against biotic stresses mainly against fungal pathogens. Twelve genes were selected and corresponding nucleotide sequences with their accession numbers were retrieved from NCBI database. The information on plant sources, sequence length (bp) and length of coding sequences (CDS) were collected and tabulated (Table-3). The individual gene sequence was further subjected to blast to NCBI database (http://www.ncbi.nlm.nih.gov/blast) to get best hit ortholog sequences belonging to the Fabaceae crop family where the maximum cut off was fixed upto 80% [14]. All the identified gene sequences from same family were retrieved and subjected to multiple alignments using MULTALIN program to identify the conserved region [15]. The CDS region, which was shared by common conserved region of maximum number of accessions was marked and selected to find out longest possible Open Reading Frame (ORF). The longest stretches of ORF sequence without any intermediate stop codon were subjected to amino acid translation using Expassy translate tools to get the putative protein sequences (http://web.expasy.org/translate/). The translated amino acid sequences were finally used as an input in protein motif finder tool to know whether translated amino acid sequences give the most stable target protein motif or not and for this Conserved Domains Database (CDD) was used (http://www.ncbi.nlm.nih.gov/cdd/). On the basis of obtained functional motif of each protein, the normal and special functional behaviour of said protein was figured out. The specific functions of the proteins were further validated with reported literatures to confirm whether they are actively involved in fungal or any other biotic stress response or not.

Primer designing and constructing RT-PCR profile
Oligonucleotide primer designing for RT-PCR experiment was performed for four defense responsive genes viz. Chalcone synthase and β-1,3-glucanase representing pathogenesis related (PR) and Ascorbate peroxidase and Superoxide dismutase representing antioxidant (AO) gene families which were initially screened through in-silico based gene targeting study as discussed earlier. The online tool of Integrated DNA Technology (http://eu.idtdna.com/primerquest/home/index) was used to design the primers. Total RNA of good quality with adequate quantity from healthy and fungal stress imposed pigeonpea seedlings was isolated after a level of optimization of classical RNA isolation technique [16] The experimental ImProm-II™ Reverse Transcription System (Promega, USA) was used to synthesize first-strand cDNA from isolated RNA. For each reaction, experimental RNA (up to 1μg) was combined with the cDNA primer [Oligo (dT) 15 @ 0.5μg/reaction] in nuclease-free water to make a final volume of 5 μl per RT reaction. All the subsequent reactions of RT-PCR were strictly followed as per manufacturer's instruction. The recommended temperature-time profile was followed after keeping the final 20 μl RT reaction in a thermal cycler. In the next step of PCR reactions, the cDNA mixture was further amplified by PCR (GenAmp, Applied Biosystems, USA) with 25 µl of total reaction volume.
As the calculated Tm values for all selected four primers were in the range between 52 o C to 60 o C, the gradient scale was adjusted for individual gene by keeping five replicate reactions in different gradient blocks. After amplification, all samples were stored at -20 o C until further use. The PCR amplified products were gel electrophoresed (1.2%), visualized with UV Trans-illuminator, and photographed using the Alpha Digi DocTM system.

Expression profiling of FW stress responsive genes by qPCR in pigeonpea
The cDNA were quantified using Nanodrop spectrophotometer (Thermo Scientific NanoDrop 2000C Technologies, Wilmington, USA) and the integrity of cDNA was checked by 1.2% agarose gel electrophoresis and ethidium bromide staining. For determining the expression profile for known PR and AO genes, qPCR study was conducted. The specific oligonucleotide primers for qPCR were designed by IDT tool (Table 4). IF4α-Initiation component 4a gene as most stable housekeeping gene in pigeonpea, the gene was considered as an internal standard for endogenous control [17]. The PCR reaction mixture was made by combining the following chemicals: 2x KAPA SYBR® FAST qPCR Master Mix (Wilmington, USA) @10 μl, PCR primers (Forward and Reverse) @ 10 mM-1 μl each, cDNA template 40 ng/μl-5 μl and PCR grade water-3 μl. The qPCR was carried out by completing 30 cycles using LightCycler ® 480 Real-Time PCR System (Roche, USA) by following the standard cycles as mentioned in instruction manual. All the technical replications of qRT-PCR experiments were conducted in triplicate and repeated twice to consider two biological replications and each time the results did not vary. By utilizing the quantitative value of cDNAs, expression profiling of individual gene was quantified at variable stress conditions by using algorithms truly based on 2 -ΔΔCT method [18].

In silico characterization of stress responsive genes in correlation with fungal infection
Out of four gene transcripts amplified in RT-PCR, two candidate genes namely chalcone synthase (CHS) and superoxide dismutase (SOD) representing PR and AO gene families, respectively were sequenced from both side (Xcelris Genomics, Gujarat, India). The obtained sequences were aligned and purified. To validate, the purified sequences of CHS and SOD were NCBI-BLAST and submitted into GenBank. The transcript sequences of both genes were analyzed at in-silico level. Gene predictions in CHS and SOD sequences were performed by FGENESH (http://linux1.softberry.com/) using Eukaryotic gene finder module [19]. The active sites of the transcripts were predicted by Prosite online server of ExPasy [20] (https://prosite.expasy.org/). The motifs of secondary protein structure were created by PDBsum using v.3.0 of Gail Hutchinson's PROMOTIF program [21]. The created motifs were then used for 3D modelling and structure prediction. CHS and SOD exist in the homodimer state in the cell, therefore, the 3D models were predicted using homology-based modelling by SWISS-MODEL (https://swissmodel.expasy.org/) using default parameters [22]. The 3D structures and predicted models were visualized using UCSF Chimera v.1.13.1 software [23]. The energy minimized models were used for generation of Ramachandran plot of the proteins using PDBsum server (http://www.ebi.ac.uk/thornton-srv/databases/pdbsum). In order to calculate physicochemical properties of the protein sequences of both proteins, CHS and SOD Protparam software was used (http://web.expasy.org/protparam/).The QMEAN server (https://swissmodel.expasy.org/qmean/) was used to find the overall quality of 3D structure [24]. The 3D model for ligand (CBH-C) cellobiohydrolase-c was predicted using SWISS-MODEL and structure analysis was carried out as above mentioned. The protein-protein docking was performed by FRODOCK version 2.0 for Interactive protein-protein docking [25]. Further simulation study was carried out using online server GRAMM-X Protein-Protein Docking Web Server v.     Table 4 List of selected primers used for RT-PCR based identification and qPCR based expression analysis of biotic stress responsive genes.  Figure 1 Schematic representation of the linking pathways between present study and molecular pathogenesis of F. udumconcomitant with ROS regeneration and PR protein induction.