Characterizing genetic diversity of Sclerotium rolfsii isolates by biomapping of mycelial compatibility groupings and multilocus sequence analysis.

Genetic diversity in Sclerotium rolfsii is useful for understanding its population structure, identifying different mycelial compatibility groups (MCGs), and developing targeted strategies for disease management in affected crops. In our study, a comprehensive genetic analysis was conducted on 50 isolates of S. rolfsii, collected from various geographic regions and host plants. Two specific genes, TEF1α and RPB2, were utilized to assess the genetic diversity and relationships among these isolates. Notably, out of 1225 pairings examined, only 154 exhibited a compatible reaction, while the majority displayed antagonistic reactions, resulting in the formation of a barrier zone. The isolates were grouped into 10 distinct MCGs. These MCGs were further characterized using genetic sequencing. TEF1α sequences distinguished the isolates into 17 distinct clusters, and RPB2 sequences classified them into 20 clusters. Some MCGs shared identical gene sequences within each gene, while others exhibited unique sequences. Intriguingly, when both TEF1α and RPB2 sequences were combined, all 10 MCGs were effectively differentiated, even those that appeared identical with single-gene analysis. This combined approach provided a comprehensive understanding of the genetic diversity and relationships among the S. rolfsii isolates, allowing for precise discrimination between different MCGs. The results shed light on the population structure and genetic variability within this plant pathogenic fungus, providing valuable insights for disease management and control strategies. This study highlights the significance of comprehending the varied virulence characteristics within S. rolfsii isolates, categorizing them into specific virulence groups based on disease severity index (DSI) values. The association with MCGs provides additional insights into the genetic underpinnings of virulence in this pathogen. Furthermore, the identification of geographical patterns in virulence implies the influence of region-specific factors, with potential implications for disease control and crop protection strategies.Please confirm if the author names are presented accurately and in the correct sequence (given name, middle name/initial, family name). Author 1 Given name: [G. M. Sandeep] Last name [Kumar]. Author 2 Given name: [Praveen Kumar] Last name [Singh]. Also, kindly confirm the details in the metadata are correct.I confirm that the given names are accurate and presented in the correct sequence.

Microbiome of Pukzing Cave in India shows high antimicrobial activity against plant and animal pathogens.

Pukzing cave, the largest cave of Mizoram, India was explored for bacterial diversity. Culture dependent method revealed 235 bacterial isolates using three different treatments. Identity of the microbial species was confirmed by 16S rDNA sequencing. The highest bacterial population was recovered from heat treatment (n = 97;41.2%) followed by normal (n = 79;33.6%) and cold treatment (n = 59;25.1%) indicating dominance of moderate thermophiles. Antimicrobial potential of isolates showed 20.4% isolates having antimicrobial ability against tested pathogens. Amplicon sequencing of PKSI, PKSII and NRP specific genes revealed presence of AMP genes in the microbial population. Six microbial pathogens were selected for screening as they are well known for different disease cause organism in various fields such as agriculture and human health. Cave environment harbors unique microbial flora and hypervariable region V4 is more informative. Higher activity of AMP assay against these microbes indicates that cave microbial communities could be potential source of future genomic resources.

Genome-Wide Association Study and Post-genome-Wide Association Study Analysis for Spike Fertility and Yield Related Traits in Bread Wheat.

Spike fertility and associated traits are key factors in deciding the grain yield potential of wheat. Genome-wide association study (GWAS) interwoven with advanced post-GWAS analysis such as a genotype-phenotype network (geno-pheno network) for spike fertility, grain yield, and associated traits allow to identify of novel genomic regions and represents attractive targets for future marker-assisted wheat improvement programs. In this study, GWAS was performed on 200 diverse wheat genotypes using Breeders' 35K Axiom array that led to the identification of 255 significant marker-trait associations (MTAs) (-log10P ≥ 3) for 15 metric traits phenotyped over three consecutive years. MTAs detected on chromosomes 3A, 3D, 5B, and 6A were most promising for spike fertility, grain yield, and associated traits. Furthermore, the geno-pheno network prioritised 11 significant MTAs that can be utilised as a minimal marker system for improving spike fertility and yield traits. In total, 119 MTAs were linked to 81 candidate genes encoding different types of functional proteins involved in various key pathways that affect the studied traits either way. Twenty-two novel loci were identified in present GWAS, twelve of which overlapped by candidate genes. These results were further validated by the gene expression analysis, Knetminer, and protein modelling. MTAs identified from this study hold promise for improving yield and related traits in wheat for continued genetic gain and in rapidly evolving artificial intelligence (AI) tools to apply in the breeding program.

Muscle transcriptome signature and gene regulatory network analysis in two divergent lines of a hilly bovine species Mithun (Bos frontalis).

A massive bovine, Bos frontalis, also known as Mithun or Gayal, found at higher altitude is very promising meat and milk animal. For candidate gene and marker discovery, RNA-seq data was generated from longissimus dorsi muscle tissues with Illumina-HiSeq. Such markers can be used in future for genetic gain of traits like feed conversion efficiency (FCE) and average daily gain (ADG). Analysis revealed 297differentially expressed genes (DEGs) having 173 up and 124 down-regulated unigenes. Extensive conservation was found in genic region while comparing with Bos taurus. Analysis revealed 57 pathways having 112 enzymes, 72 transcriptional factors and cofactors, 212 miRNAs regulating 71 DEGs, 25,855 SSRs, mithun-specific 104,822 variants and 7288 indels, gene regulatory network (GRN) having 24 hub-genes and transcriptional factors regulating cell proliferation, immune tolerance and myogenesis. This is first report of muscle transcriptome depicting candidate genes with GRN controlling FCE and ADG. Reported putative molecular markers, candidate genes and hub proteins can be valuable genomic resources for association studies in genetic improvement programme.

Characterization of genetic diversity and population structure in wheat using array based SNP markers.

Genetic diversity is crucial for successful adaptation and sustained improvement in crops. India is bestowed with diverse agro-climatic conditions which makes it rich in wheat germplasm adapted to various niches. Germplasm repository consists of local landraces, trait specific genetic stocks including introgressions from wild relatives, exotic collections, released varieties, and improved germplasm. Characterization of genetic diversity is done using morpho-physiological characters as well as by analyzing variations at DNA level. However, there are not many reports on array based high throughput SNP markers having characteristics of genome wide coverage employed in Indian spring wheat germplasm. Amongst wheat SNP arrays, 35K Axiom Wheat Breeder's Array has the highest SNP polymorphism efficiency suitable for genetic mapping and genetic diversity characterization. Therefore, genotyping was done using 35K in 483 wheat genotypes resulting in 14,650 quality filtered SNPs, that were distributed across the B (~ 50%), A (~ 39%), and D (~ 10%) genomes. The total genetic distance coverage was 4477.85 cM with 3.27 SNP/cM and 0.49 cM/SNP as average marker density and average inter-marker distance, respectively. The PIC ranged from 0.09 to 0.38 with an average of 0.29 across genomes. Population structure and Principal Coordinate Analysis resulted in two subpopulations (SP1 and SP2). The analysis of molecular variance revealed the genetic variation of 2% among and 98% within subpopulations indicating high gene flow between SP1 and SP2. The subpopulation SP2 showed high level of genetic diversity based on genetic diversity indices viz. Shannon's information index (I) = 0.648, expected heterozygosity (He) = 0.456 and unbiased expected heterozygosity (uHe) = 0.456. To the best of our knowledge, this study is the first to include the largest set of Indian wheat genotypes studied exclusively for genetic diversity. These findings may serve as a potential source for the identification of uncharacterized QTL/gene using genome wide association studies and marker assisted selection in wheat breeding programs.

Harmonizing technological advances in phenomics and genomics for enhanced salt tolerance in rice from a practical perspective.

Half of the global human population is dependent on rice as a staple food crop and more than 25% increase in rice productivity is required to feed the global population by 2030. With increase in irrigation, global warming and rising sea level, rising salinity has become one of the major challenges to enhance the rice productivity. Since the loss on this account is to the tune of US$12 billion per annum, it necessitates the global attention. In the era of technological advancement, substantial progress has been made on phenomics and genomics data generation but reaping benefit of this in rice salinity variety development in terms of cost, time and precision requires their harmonization. There is hardly any comprehensive holistic review for such combined approach. Present review describes classical salinity phenotyping approaches having morphological, physiological and biochemical components. It also gives a detailed account of invasive and non-invasive approaches of phenomic data generation and utilization. Classical work of rice salinity QLTs mapping in the form of chromosomal atlas has been updated. This review describes how QTLs can be further dissected into QTN by GWAS and transcriptomic approaches. Opportunities and progress made by transgenic, genome editing, metagenomics approaches in combating rice salinity problems are discussed. Major aim of this review is to provide a comprehensive over-view of hitherto progress made in rice salinity tolerance research which is required to understand bridging of phenotype based breeding with molecular breeding. This review is expected to assist rice breeders in their endeavours by fetching greater harmonization of technological advances in phenomics and genomics for better pragmatic approach having practical perspective.

Uncovering Genomic Regions Associated With 36 Agro-Morphological Traits in Indian Spring Wheat Using GWAS.

Wheat genetic improvement by integration of advanced genomic technologies is one way of improving productivity. To facilitate the breeding of economically important traits in wheat, SNP loci and underlying candidate genes associated with the 36 agro-morphological traits were studied in a diverse panel of 404 genotypes. By using Breeders' 35K Axiom array in a comprehensive genome-wide association study covering 4364.79 cM of the wheat genome and applying a compressed mixed linear model, a total of 146 SNPs (-log10 P ≥ 4) were found associated with 23 traits out of 36 traits studied explaining 3.7-47.0% of phenotypic variance. To reveal this a subset of 260 genotypes was characterized phenotypically for six quantitative traits [days to heading (DTH), days to maturity (DTM), plant height (PH), spike length (SL), awn length (Awn_L), and leaf length (Leaf_L)] under five environments. Gene annotations mined ∼38 putative candidate genes which were confirmed using tissue and stage specific gene expression data from RNA Seq. We observed strong co-localized loci for four traits (glume pubescence, SL, PH, and awn color) on chromosome 1B (24.64 cM) annotated five putative candidate genes. This study led to the discovery of hitherto unreported loci for some less explored traits (such as leaf sheath wax, awn attitude, and glume pubescence) besides the refined chromosomal regions of known loci associated with the traits. This study provides valuable information of the genetic loci and their potential genes underlying the traits such as awn characters which are being considered as important contributors toward yield enhancement.

Development of species specific putative miRNA and its target prediction tool in wheat (Triticum aestivum L.).

MicroRNA are 20-24 nt, non-coding, single stranded molecule regulating traits and stress response. Tissue and time specific expression limits its detection, thus is major challenge in their discovery. Wheat has limited 119 miRNAs in MiRBase due to limitation of conservation based methodology where old and new miRNA genes gets excluded. This is due to origin of hexaploid wheat by three successive hybridization, older AA, BB and younger DD subgenome. Species specific miRNA prediction (SMIRP concept) based on 152 thermodynamic features of training dataset using support vector machine learning approach has improved prediction accuracy to 97.7%. This has been implemented in TamiRPred ( http://webtom.cabgrid.res.in/tamirpred ). We also report highest number of putative miRNA genes (4464) of wheat from whole genome sequence populated in database developed in PHP and MySQL. TamiRPred has predicted 2092 (>45.10%) additional miRNA which was not predicted by miRLocator. Predicted miRNAs have been validated by miRBase, small RNA libraries, secondary structure, degradome dataset, star miRNA and binding sites in wheat coding region. This tool can accelerate miRNA polymorphism discovery to be used in wheat trait improvement. Since it predicts chromosome-wise miRNA genes with their respective physical location thus can be transferred using linked SSR markers. This prediction approach can be used as model even in other polyploid crops.

Development of model web-server for crop variety identification using throughput SNP genotyping data.

Crop varieties or genotypes of a given species are pivotal for agricultural production and ownership, management and improvement of their germplasm is a great challenge. Its morphological identification requires time, cost and descriptors are often compromised statistically due to phenotypic plasticity. Development of DNA based signature of varieties can overcome these limitations. There is a global need to implement world trade organization (WTO) and intellectual property rights (IPR) guidelines of Plant Breeders Rights (PBR) where DUS (distinctness, uniformity and stability) testing can be supplemented by DNA profile. Universalization and minimization of SNP number without compromising identification accuracy is the major challenge in development of varietal profile by rapid genotype assay. Besides this, there is no server-based approach reducing computational skill with global accessibility of referral phenotypic and genotypic data. We report world's first model web server for crop variety identification using >350 Indian wheat varieties and Axiom 35 K SNP chip data. Standard filtering and linkage disequilibrium approach were used to develop varietal signature in Linux using HTML, Java, PHP and MySQL with provision of QR code generator to facilitate bar-coding. Phylogenetic tree constructed by selected SNPs confirms six major trait based clusters of varieties and their pedigree. Our user friendly server based tool, VISTa (Variety Identification System of Triticum aestivum) ( http://webtom.cabgrid.res.in/vista ) can be used in DUS testing having dispute resolution of sovereignty and access benefit sharing (ABS) issues. This model approach can be used in other crops with pan-global level management of crop germplasm in endeavour of crop productivity.

Author Correction: Transcriptomic signature of drought response in pearl millet (Pennisetum glaucum (L.) and development of web-genomic resources.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has been fixed in the paper.

Transcriptomic signature of drought response in pearl millet (Pennisetum glaucum (L.) and development of web-genomic resources.

Pearl millet, (Pennisetum glaucum L.), an efficient (C4) crop of arid/semi-arid regions is known for hardiness. Crop is valuable for bio-fortification combating malnutrition and diabetes, higher caloric value and wider climatic resilience. Limited studies are done in pot-based experiments for drought response at gene-expression level, but field-based experiment mimicking drought by withdrawal of irrigation is still warranted. We report de novo assembly-based transcriptomic signature of drought response induced by irrigation withdrawal in pearl millet. We found 19983 differentially expressed genes, 7595 transcription factors, gene regulatory network having 45 hub genes controlling drought response. We report 34652 putative markers (4192 simple sequence repeats, 12111 SNPs and 6249 InDels). Study reveals role of purine and tryptophan metabolism in ABA accumulation mediating abiotic response in which MAPK acts as major intracellular signal sensing drought. Results were validated by qPCR of 13 randomly selected genes. We report the first web-based genomic resource ( http://webtom.cabgrid.res.in/pmdtdb/ ) which can be used for candidate genes-based SNP discovery programs and trait-based association studies. Looking at climatic change, nutritional and pharmaceutical importance of this crop, present investigation has immense value in understanding drought response in field condition. This is important in germplasm management and improvement in endeavour of pearl millet productivity.

PolyMorphPredict: A Universal Web-Tool for Rapid Polymorphic Microsatellite Marker Discovery From Whole Genome and Transcriptome Data.

Microsatellites are ubiquitously distributed, polymorphic repeat sequence valuable for association, selection, population structure and identification. They can be mined by genomic library, probe hybridization and sequencing of selected clones. Such approach has many limitations like biased hybridization and selection of larger repeats. In silico mining of polymorphic markers using data of various genotypes can be rapid and economical. Available tools lack in some or other aspects like: targeted user defined primer generation, polymorphism discovery using multiple sequence, size and number limits of input sequence, no option for primer generation and e-PCR evaluation, transferability, lack of complete automation and user-friendliness. They also lack the provision to evaluate published primers in e-PCR mode to generate additional allelic data using re-sequenced data of various genotypes for judicious utilization of previously generated data. We developed the tool (PolyMorphPredict) using Perl, R, Java and launched at Apache which is available at http://webtom.cabgrid.res.in/polypred/. It mines microsatellite loci and computes primers from genome/transcriptome data of any species. It can perform e-PCR using published primers for polymorphism discovery and across species transferability of microsatellite loci. Present tool has been evaluated using five species of different genome size having 21 genotypes. Though server is equipped with genomic data of three species for test run with gel simulation, but can be used for any species. Further, polymorphism predictability has been validated using in silico and in vitro PCR of four rice genotypes. This tool can accelerate the in silico microsatellite polymorphism discovery in re-sequencing projects of any species of plant and animal for their diversity estimation along with variety/breed identification, population structure, MAS, QTL and gene discovery, traceability, parentage testing, fungal diagnostics and genome finishing.

MiSNPDb: a web-based genomic resources of tropical ecology fruit mango (Mangifera indica L.) for phylogeography and varietal differentiation.

Mango is one of the most important fruits of tropical ecological region of the world, well known for its nutritive value, aroma and taste. Its world production is >45MT worth >200 billion US dollars. Genomic resources are required for improvement in productivity and management of mango germplasm. There is no web-based genomic resources available for mango. Hence rapid and cost-effective high throughput putative marker discovery is required to develop such resources. RAD-based marker discovery can cater this urgent need till whole genome sequence of mango becomes available. Using a panel of 84 mango varieties, a total of 28.6 Gb data was generated by ddRAD-Seq approach on Illumina HiSeq 2000 platform. A total of 1.25 million SNPs were discovered. Phylogenetic tree using 749 common SNPs across these varieties revealed three major lineages which was compared with geographical locations. A web genomic resources MiSNPDb, available at http://webtom.cabgrid.res.in/mangosnps/ is based on 3-tier architecture, developed using PHP, MySQL and Javascript. This web genomic resources can be of immense use in the development of high density linkage map, QTL discovery, varietal differentiation, traceability, genome finishing and SNP chip development for future GWAS in genomic selection program. We report here world's first web-based genomic resources for genetic improvement and germplasm management of mango.

Computational deciphering of biotic stress associated genes in tomato (Solanum lycopersicum).

Tomato (Solanum lycopersicum) is one of the major vegetable plant and a model system for fruit development. Its global importance is due to its lycopene pigment which has anti-oxidative and anti-cancerous properties. Though > 1.5 M biotic stress associated ESTs of tomato are available but cumulative analysis to predict genes is warranted. Availability of whole genome de novo assembly can advantageously be used to map them over different chromosome. Further, available 0.14 M catalogued markers can be used to introgress specific desirable genes in varietal improvement program. We report here 57 novel genes associated with biotic stress of tomato along with 50 genes having physical location over different chromosomes. We also report 52 cis-regulating elements and 69 putative miRNAs which are involved in regulation of these biotic stresses associated genes. These putative candidate genes associated with biotic stress can be used in molecular breeding in the endeavor of tomato productivity along with its sustainable germplasm management.

Draft whole genome sequence of groundnut stem rot fungus Athelia rolfsii revealing genetic architect of its pathogenicity and virulence.

Groundnut (Arachis hypogaea L.) is an important oil seed crop having major biotic constraint in production due to stem rot disease caused by fungus, Athelia rolfsii causing 25-80% loss in productivity. As chemical and biological combating strategies of this fungus are not very effective, thus genome sequencing can reveal virulence and pathogenicity related genes for better understanding of the host-parasite interaction. We report draft assembly of Athelia rolfsii genome of ~73 Mb having 8919 contigs. Annotation analysis revealed 16830 genes which are involved in fungicide resistance, virulence and pathogenicity along with putative effector and lethal genes. Secretome analysis revealed CAZY genes representing 1085 enzymatic genes, glycoside hydrolases, carbohydrate esterases, carbohydrate-binding modules, auxillary activities, glycosyl transferases and polysaccharide lyases. Repeat analysis revealed 11171 SSRs, LTR, GYPSY and COPIA elements. Comparative analysis with other existing ascomycotina genome predicted conserved domain family of WD40, CYP450, Pkinase and ABC transporter revealing insight of evolution of pathogenicity and virulence. This study would help in understanding pathogenicity and virulence at molecular level and development of new combating strategies. Such approach is imperative in endeavour of genome based solution in stem rot disease management leading to better productivity of groundnut crop in tropical region of world.

Draft Genome Sequence of Two Monosporidial Lines of the Karnal Bunt Fungus Tilletia indica Mitra (PSWKBGH-1 and PSWKBGH-2).

Karnal bunt disease caused by the fungus Tilletia indica Mitra is a serious concern due to strict quarantines affecting international trade of wheat. We announce here the first draft assembly of two monosporidial lines, PSWKBGH-1 and -2, of this fungus, having approximate sizes of 37.46 and 37.21 Mbp, respectively.

Development of Antimicrobial Peptide Prediction Tool for Aquaculture Industries.

Microbial diseases in fish, plant, animal and human are rising constantly; thus, discovery of their antidote is imperative. The use of antibiotic in aquaculture further compounds the problem by development of resistance and consequent consumer health risk by bio-magnification. Antimicrobial peptides (AMPs) have been highly promising as natural alternative to chemical antibiotics. Though AMPs are molecules of innate immune defense of all advance eukaryotic organisms, fish being heavily dependent on their innate immune defense has been a good source of AMPs with much wider applicability. Machine learning-based prediction method using wet laboratory-validated fish AMP can accelerate the AMP discovery using available fish genomic and proteomic data. Earlier AMP prediction servers are based on multi-phyla/species data, and we report here the world's first AMP prediction server in fishes. It is freely accessible at http://webapp.cabgrid.res.in/fishamp/ . A total of 151 AMPs related to fish collected from various databases and published literature were taken for this study. For model development and prediction, N-terminus residues, C-terminus residues and full sequences were considered. Best models were with kernels polynomial-2, linear and radial basis function with accuracy of 97, 99 and 97 %, respectively. We found that performance of support vector machine-based models is superior to artificial neural network. This in silico approach can drastically reduce the time and cost of AMP discovery. This accelerated discovery of lead AMP molecules having potential wider applications in diverse area like fish and human health as substitute of antibiotics, immunomodulator, antitumor, vaccine adjuvant and inactivator, and also for packaged food can be of much importance for industries.

Low-depth shotgun sequencing resolves complete mitochondrial genome sequence of Labeo rohita.

Labeo rohita, popularly known as rohu, is a widely cultured species in whole Indian subcontinent. In the present study, we used in-silico approach to resolve complete mitochondrial genome of rohu. Low-depth shotgun sequencing using Roche 454 GS FLX (Branford, Connecticut, USA) followed by de novo assembly in CLC Genomics Workbench version 7.0.4 (Aarhus, Denmark) revealed the complete mitogenome of L. rohita to be 16 606 bp long (accession No. KR185963). It comprised of 13 protein-coding genes, 22 tRNAs, 2 rRNAs and 1 putative control region. The gene order and organization are similar to most vertebrates. The mitogenome in the present investigation has 99% similarity with that of previously reported mitogenomes of rohu and this is also evident from the phylogenetic study using maximum-likelihood (ML) tree method. This study was done to determine the feasibility, accuracy and reliability of low-depth sequence data obtained from NGS platform as compared to the Sanger sequencing. Thus, NGS technology has proven to be competent and a rapid in-silico alternative to resolve the complete mitochondrial genome sequence, thereby reducing labors and time.

Assembly and variation analyses of Clarias batrachus mitogenome retrieved from WGS data and its phylogenetic relationship with other catfishes.

Whole genome sequencing (WGS) using next generation sequencing technologies paves the way to sequence the mitochondrial genomes with greater ease and lesser time. Here, we used the WGS data of Clarias batrachus, generated from Roche 454 and Ion Torrent sequencing platforms, to assemble the complete mitogenome using both de novo and reference based approaches. Both the methods yielded almost similar results and the best assembled mitogenome was of 16,510 bp size (GenBank Acc. No. KM259918). The mitogenome annotation resulted in 13 coding genes, 22 tRNA genes, 2 rRNA genes and one control region, and the gene order was found to be identical with other catfishes. Variation analyses between assembled and the reference (GenBank Acc. No. NC_023923) mitogenome revealed 51 variations. The phylogenetic analysis of coding DNA sequences and tRNA supports the monophyly of catfishes. Two SSRs were identified in C. batrachus mitogenome, out of which one was unique to this species. Based on the relative rate of gene evolution, protein coding mitochondrial genes were found to evolve at a much faster pace than the d-loop, which in turn are followed by the rRNAs; the tRNAs showed wide variability in the rate of sequence evolution, and on average evolve the slowest. Among the coding genes, ND2 evolves most rapidly. The variations present in the coding regions of the mitogenome and their comparative analyses with other catfish species may be useful in species conservation and management programs.

Evidence of salicylic acid pathway with EDS1 and PAD4 proteins by molecular dynamics simulation for grape improvement.

Biotic stress is a major cause of heavy loss in grape productivity. In order to develop biotic stress-resistant grape varieties, the key defense genes along with its pathway have to be deciphered. In angiosperm plants, lipase-like protein phytoalexin deficient 4 (PAD4) is well known to be essential for systemic resistance against biotic stress. PAD4 functions together with its interacting partner protein enhanced disease susceptibility 1 (EDS1) to promote salicylic acid (SA)-dependent and SA-independent defense pathway. Existence and structure of key protein of systemic resistance EDS1 and PAD4 are not known in grapes. Before SA pathway studies are taken in grape, molecular evidence of EDS1: PAD4 complex is to be established. To establish this, EDS1 protein sequence was retrieved from NCBI and homologous PAD4 protein was generated using Arabidopsis thaliana as template and conserved domains were confirmed. In this study, computational methods were used to model EDS1 and PAD4 and simulated the interactions of EDS1 and PAD4. Since no structural details of the proteins were available, homology modeling was employed to construct three-dimensional structures. Further, molecular dynamic simulations were performed to study the dynamic behavior of the EDS1 and PAD4. The modeled proteins were validated and subjected to molecular docking analysis. Molecular evidence of stable complex of EDS1:PAD4 in grape supporting SA defense pathway in response to biotic stress is reported in this study. If SA defense pathway genes are explored, then markers of genes involved can play pivotal role in grape variety development especially against biotic stress leading to higher productivity.