Genome-wide characterization of the NBLRR gene family provides evolutionary and functional insights into blast resistance in pearl millet (Cenchrus americanus (L.) Morrone).

MAIN CONCLUSION: The investigation is the first report on genome-wide identification and characterization of NBLRR genes in pearl millet. We have shown the role of gene loss and purifying selection in the divergence of NBLRRs in Poaceae lineage and candidate CaNBLRR genes for resistance to Magnaporthe grisea infection. Plants have evolved multiple integral mechanisms to counteract the pathogens' infection, among which plant immunity through NBLRR (nucleotide-binding site, leucine-rich repeat) genes is at the forefront. The genome-wide mining in pearl millet (Cenchrus americanus (L.) Morrone) revealed 146 CaNBLRRs. The variation in the branch length of NBLRRs showed the dynamic nature of NBLRRs in response to evolving pathogen races. The orthology of NBLRRs showed a predominance of many-to-one orthologs, indicating the divergence of NBLRRs in the pearl millet lineage mainly through gene loss events followed by gene gain through single-copy duplications. Further, the purifying selection (Ka/Ks < 1) shaped the expansion of NBLRRs within the lineage of pear millet and other members of Poaceae. Presence of cis-acting elements, viz. TCA element, G-box, MYB, SARE, ABRE and conserved motifs annotated with P-loop, kinase 2, RNBS-A, RNBS-D, GLPL, MHD, Rx-CC and LRR suggests their putative role in disease resistance and stress regulation. The qRT-PCR analysis in pearl millet lines showing contrasting responses to Magnaporthe grisea infection identified CaNBLRR20, CaNBLRR33, CaNBLRR46 CaNBLRR51, CaNBLRR78 and CaNBLRR146 as putative candidates. Molecular docking showed the involvement of three and two amino acid residues of LRR domains forming hydrogen bonds (histidine, arginine and threonine) and salt bridges (arginine and lysine) with effectors. Whereas 14 and 20 amino acid residues of CaNBLRR78 and CaNBLRR20 showed hydrophobic interactions with 11 and 9 amino acid residues of effectors, Mg.00g064570.m01 and Mg.00g006570.m01, respectively. The present investigation gives a comprehensive overview of CaNBLRRs and paves the foundation for their utility in pearl millet resistance breeding through understanding of host-pathogen interactions.

Comprehensive genomic analysis of Bacillus subtilis and Bacillus paralicheniformis associated with the pearl millet panicle reveals their antimicrobial potential against important plant pathogens.

BACKGROUND: Plant microbiome confers versatile functional roles to enhance survival fitness as well as productivity. In the present study two pearl millet panicle microbiome member species Bacillus subtilis PBs 12 and Bacillus paralicheniformis PBl 36 found to have beneficial traits including plant growth promotion and broad-spectrum antifungal activity towards taxonomically diverse plant pathogens. Understanding the genomes will assist in devising a bioformulation for crop protection while exploiting their beneficial functional roles. RESULTS: Two potential firmicute species were isolated from pearl millet panicles. Morphological, biochemical, and molecular characterization revealed their identities as Bacillus subtilis PBs 12 and Bacillus paralicheniformis PBl 36. The seed priming assays revealed the ability of both species to enhance plant growth promotion and seedling vigour index. Invitro assays with PBs 12 and PBl 36 showed the antibiosis effect against taxonomically diverse plant pathogens (Magnaporthe grisea; Sclerotium rolfsii; Fusarium solani; Alternaria alternata; Ganoderma sp.) of crops and multipurpose tree species. The whole genome sequence analysis was performed to unveil the genetic potential of these bacteria for plant protection. The complete genomes of PBs 12 and PBl 36 consist of a single circular chromosome with a size of 4.02 and 4.33 Mb and 4,171 and 4,606 genes, with a G + C content of 43.68 and 45.83%, respectively. Comparative Average Nucleotide Identity (ANI) analysis revealed a close similarity of PBs 12 and PBl 36 with other beneficial strains of B. subtilis and B. paralicheniformis and found distant from B. altitudinis, B. amyloliquefaciens, and B. thuringiensis. Functional annotation revealed a majority of pathway classes of PBs 12 (30) and PBl 36 (29) involved in the biosynthesis of secondary metabolites, polyketides, and non-ribosomal peptides, followed by xenobiotic biodegradation and metabolism (21). Furthermore, 14 genomic regions of PBs 12 and 15 of PBl 36 associated with the synthesis of RiPP (Ribosomally synthesized and post-translationally modified peptides), terpenes, cyclic dipeptides (CDPs), type III polyketide synthases (T3PKSs), sactipeptides, lanthipeptides, siderophores, NRPS (Non-Ribosomal Peptide Synthetase), NRP-metallophone, etc. It was discovered that these areas contain between 25,458 and 33,000 secondary metabolite-coding MiBiG clusters which code for a wide range of products, such as antibiotics. The PCR-based screening for the presence of antimicrobial peptide (cyclic lipopeptide) genes in PBs 12 and 36 confirmed their broad-spectrum antifungal potential with the presence of spoVG, bacA, and srfAA AMP genes, which encode antimicrobial compounds such as subtilin, bacylisin, and surfactin. CONCLUSION: The combined in vitro studies and genome analysis highlighted the antifungal potential of pearl millet panicle-associated Bacillus subtilis PBs12 and Bacillus paralicheniformis PBl36. The genetic ability to synthesize several antimicrobial compounds indicated the industrial value of PBs 12 and PBl 36, which shed light on further studies to establish their action as a biostimulant for crop protection.

Genome-wide identification and characterization of NBLRR genes in finger millet (Eleusine coracana L.) and their expression in response to Magnaporthe grisea infection.

BACKGROUND: The nucleotide binding site leucine rich repeat (NBLRR) genes significantly regulate defences against phytopathogens in plants. The genome-wide identification and analysis of NBLRR genes have been performed in several species. However, the detailed evolution, structure, expression of NBLRRs and functional response to Magnaporthe grisea are unknown in finger millet (Eleusine coracana (L.) Gaertn.). RESULTS: The genome-wide scanning of the finger millet genome resulted in 116 NBLRR (EcNBLRRs1-116) encompassing 64 CC-NB-LRR, 47 NB-LRR and 5 CCR-NB-LRR types. The evolutionary studies among the NBLRRs of five Gramineae species, viz., purple false brome (Brachypodium distachyon (L.) P.Beauv.), finger millet (E. coracana), rice (Oryza sativa L.), sorghum (Sorghum bicolor L. (Moench)) and foxtail millet (Setaria italica (L.) P.Beauv.) showed the evolution of NBLRRs in the ancestral lineage of the target species and subsequent divergence through gene-loss events. The purifying selection (Ka/Ks < 1) shaped the expansions of NBLRRs paralogs in finger millet and orthologs among the target Gramineae species. The promoter sequence analysis showed various stress- and phytohormone-responsive cis-acting elements besides growth and development, indicating their potential role in disease defence and regulatory mechanisms. The expression analysis of 22 EcNBLRRs in the genotypes showing contrasting responses to Magnaporthe grisea infection revealed four and five EcNBLRRs in early and late infection stages, respectively. The six of these nine candidate EcNBLRRs proteins, viz., EcNBLRR21, EcNBLRR26, EcNBLRR30, EcNBLRR45, EcNBLRR55 and EcNBLRR76 showed CC, NB and LRR domains, whereas the EcNBLRR23, EcNBLRR32 and EcNBLRR83 showed NB and LRR somains. CONCLUSION: The identification and expression analysis of EcNBLRRs showed the role of EcNBLRR genes in assigning blast resistance in finger millet. These results pave the foundation for in-depth and targeted functional analysis of EcNBLRRs through genome editing and transgenic approaches.

Enhancing defense against rice blast disease: Unveiling the role of leaf endophytic firmicutes in antifungal antibiosis and induced systemic resistance.

Rice remains the primary staple for more than half of the world's population, yet its cultivation faces numerous challenges, including both biotic and abiotic stresses. One significant obstacle is the prevalence of rice blast disease, which substantially diminishes productivity and increases cultivation costs due to frequent fungicide applications. Consequently, the presence of fungicide residues in rice raises concerns about compliance with international maximum residue limits (MRLs). While host resistance has proven effective, it often remains vulnerable to new variants of the Magnaporthe oryzae pathogen. Therefore, there is a critical need to explore innovative management strategies that can complement or enhance existing methods. An unexplored avenue involves harnessing endophytic bacterial communities. To this end, the present study investigates the potential of eleven endophytic Bacillus spp. in suppressing Pyricularia oryzae, promoting plant growth, and eliciting a defense response through phyllobacterization. The results indicate that the secreted metabolome and volatilome of seven tested isolates demonstrate inhibitory effects against P.oryzae, ranging from a minimum of 40% to a maximum of 70%. Bacillus siamensis L34, B. amyloliquefaciens RA37, B. velezensis L12, and B. subtilis B18 produce antifungal antibiotics targeting P.oryzae. Additionally, B. subtilis S4 and B. subtilis S6 emerge as excellent inducers of systemic resistance against blast disease, as evidenced by elevated activity of biochemical defense enzymes such as peroxidase, polyphenol oxidase, and total phenol content. However, a balance between primary metabolic activity (e.g., chlorophyll content, chlorophyll fluorescence, and photosynthetic rate) and defense activity is observed. Furthermore, specific endophytic Bacillus spp. significantly stimulates defense-related genes, including OsPAD4, OsFMO1, and OsEDS1. These findings underscore the multifaceted potential of endophytic Bacillus in managing blast disease through antibiosis and induced systemic resistance. In conclusion, this study highlights the promising role of endophytic Bacillus spp. as a viable option for blast disease management. Their ability to inhibit the pathogen and induce systemic resistance makes them a valuable addition to the existing strategies. However, it is crucial to consider the trade-off between primary metabolic activity and defense response when implementing these bacteria-based approaches.

Composition and Codon Usage Pattern Results in Divergence of the Zinc Binuclear Cluster (Zn(II)2Cys6) Sequences among Ascomycetes Plant Pathogenic Fungi.

Zinc binuclear cluster proteins (ZBC; Zn(II)2Cys6) are unique to the fungi kingdom and associated with a series of functions, viz., the utilization of macromolecules, stress tolerance, and most importantly, host-pathogen interactions by imparting virulence to the pathogen. Codon usage bias (CUB) is the phenomenon of using synonymous codons in a non-uniform fashion during the translation event, which has arisen because of interactions among evolutionary forces. The Zn(II)2Cys6 coding sequences from nine Ascomycetes plant pathogenic species and model system yeast were analysed for compositional and codon usage bias patterns. The clustering analysis diverged the Ascomycetes fungi into two clusters. The nucleotide compositional and relative synonymous codon usage (RSCU) analysis indicated GC biasness toward Ascomycetes fungi compared with the model system S. cerevisiae, which tends to be AT-rich. Further, plant pathogenic Ascomycetes fungi belonging to cluster-2 showed a higher number of GC-rich high-frequency codons than cluster-1 and was exclusively AT-rich in S. cerevisiae. The current investigation also showed the mutual effect of the two evolutionary forces, viz. natural selection and compositional constraints, on the CUB of Zn(II)2Cys6 genes. The perseverance of GC-rich codons of Zn(II)2Cys6 in Ascomycetes could facilitate the invasion process. The findings of the current investigation show the role of CUB and nucleotide composition in the evolutionary divergence of Ascomycetes plant pathogens and paves the way to target specific codons and sequences to modulate host-pathogen interactions through genome editing and functional genomics tools.

Structured Framework and Genome Analysis of Magnaporthe grisea Inciting Pearl Millet Blast Disease Reveals Versatile Metabolic Pathways, Protein Families, and Virulence Factors.

Magnaporthe grisea (T.T. Herbert) M.E. Barr is a major fungal phytopathogen that causes blast disease in cereals, resulting in economic losses worldwide. An in-depth understanding of the basis of virulence and ecological adaptation of M. grisea is vital for devising effective disease management strategies. Here, we aimed to determine the genomic basis of the pathogenicity and underlying biochemical pathways in Magnaporthe using the genome sequence of a pearl millet-infecting M. grisea PMg_Dl generated by dual NGS techniques, Illumina NextSeq 500 and PacBio RS II. The short and long nucleotide reads could be draft assembled in 341 contigs and showed a genome size of 47.89 Mb with the N50 value of 765.4 Kb. Magnaporthe grisea PMg_Dl showed an average nucleotide identity (ANI) of 86% and 98% with M. oryzae and Pyricularia pennisetigena, respectively. The gene-calling method revealed a total of 10,218 genes and 10,184 protein-coding sequences in the genome of PMg_Dl. InterProScan of predicted protein showed a distinct 3637 protein families and 695 superfamilies in the PMg_Dl genome. In silico virulence analysis revealed the presence of 51VFs and 539 CAZymes in the genome. The genomic regions for the biosynthesis of cellulolytic endo-glucanase and beta-glucosidase, as well as pectinolytic endo-polygalacturonase, pectin-esterase, and pectate-lyases (pectinolytic) were detected. Signaling pathways modulated by MAPK, PI3K-Akt, AMPK, and mTOR were also deciphered. Multicopy sequences suggestive of transposable elements such as Type LTR, LTR/Copia, LTR/Gypsy, DNA/TcMar-Fot1, and Type LINE were recorded. The genomic resource presented here will be of use in the development of molecular marker and diagnosis, population genetics, disease management, and molecular taxonomy, and also provide a genomic reference for ascomycetous genome investigations in the future.

Deciphering core phyllomicrobiome assemblage on rice genotypes grown in contrasting agroclimatic zones: implications for phyllomicrobiome engineering against blast disease.

BACKGROUND: With its adapted microbial diversity, the phyllosphere contributes microbial metagenome to the plant holobiont and modulates a host of ecological functions. Phyllosphere microbiome (hereafter termed phyllomicrobiome) structure and the consequent ecological functions are vulnerable to a host of biotic (Genotypes) and abiotic factors (Environment) which is further compounded by agronomic transactions. However, the ecological forces driving the phyllomicrobiome assemblage and functions are among the most understudied aspects of plant biology. Despite the reports on the occurrence of diverse prokaryotic phyla such as Proteobacteria, Firmicutes, Bacteroides, and Actinobacteria in phyllosphere habitat, the functional characterization leading to their utilization for agricultural sustainability is not yet explored. Currently, the metabarcoding by Next-Generation-Sequencing (mNGS) technique is a widely practised strategy for microbiome investigations. However, the validation of mNGS annotations by culturomics methods is not integrated with the microbiome exploration program. In the present study, we combined the mNGS with culturomics to decipher the core functional phyllomicrobiome of rice genotypes varying for blast disease resistance planted in two agroclimatic zones in India. There is a growing consensus among the various stakeholder of rice farming for an ecofriendly method of disease management. Here, we proposed phyllomicrobiome assisted rice blast management as a novel strategy for rice farming in the future. RESULTS: The tropical "Island Zone" displayed marginally more bacterial diversity than that of the temperate 'Mountain Zone' on the phyllosphere. Principal coordinate analysis indicated converging phyllomicrobiome profiles on rice genotypes sharing the same agroclimatic zone. Interestingly, the rice genotype grown in the contrasting zones displayed divergent phyllomicrobiomes suggestive of the role of environment on phyllomicrobiome assembly. The predominance of phyla such as Proteobacteria, Actinobacteria, and Firmicutes was observed in the phyllosphere irrespective of the genotypes and climatic zones. The core-microbiome analysis revealed an association of Acidovorax, Arthrobacter, Bacillus, Clavibacter, Clostridium, Cronobacter, Curtobacterium, Deinococcus, Erwinia, Exiguobacterium, Hymenobacter, Kineococcus, Klebsiella, Methylobacterium, Methylocella, Microbacterium, Nocardioides, Pantoea, Pedobacter, Pseudomonas, Salmonella, Serratia, Sphingomonas and Streptomyces on phyllosphere. The linear discriminant analysis (LDA) effect size (LEfSe) method revealed distinct bacterial genera in blast-resistant and susceptible genotypes, as well as mountain and island climate zones. SparCC based network analysis of phyllomicrobiome showed complex intra-microbial cooperative or competitive interactions on the rice genotypes. The culturomic validation of mNGS data confirmed the occurrence of Acinetobacter, Aureimonas, Curtobacterium, Enterobacter, Exiguobacterium, Microbacterium, Pantoea, Pseudomonas, and Sphingomonas in the phyllosphere. Strikingly, the contrasting agroclimatic zones showed genetically identical bacterial isolates suggestive of vertical microbiome transmission. The core-phyllobacterial communities showed secreted and volatile compound mediated antifungal activity on M. oryzae. Upon phyllobacterization (a term coined for spraying bacterial cells on the phyllosphere), Acinetobacter, Aureimonas, Pantoea, and Pseudomonas conferred immunocompetence against blast disease. Transcriptional analysis revealed activation of defense genes such as OsPR1.1, OsNPR1, OsPDF2.2, OsFMO, OsPAD4, OsCEBiP, and OsCERK1 in phyllobacterized rice seedlings. CONCLUSIONS: PCoA indicated the key role of agro-climatic zones to drive phyllomicrobiome assembly on the rice genotypes. The mNGS and culturomic methods showed Acinetobacter, Aureimonas, Curtobacterium, Enterobacter, Exiguobacterium, Microbacterium, Pantoea, Pseudomonas, and Sphingomonas as core phyllomicrobiome of rice. Genetically identical Pantoea intercepted on the phyllosphere from the well-separated agroclimatic zones is suggestive of vertical transmission of phyllomicrobiome. The phyllobacterization showed potential for blast disease suppression by direct antibiosis and defense elicitation. Identification of functional core-bacterial communities on the phyllosphere and their co-occurrence dynamics presents an opportunity to devise novel strategies for rice blast management through phyllomicrobiome reengineering in the future.

Hybrid de novo genome-reassembly reveals new insights on pathways and pathogenicity determinants in rice blast pathogen Magnaporthe oryzae RMg_Dl.

Blast disease incited by Magnaporthe oryzae is a major threat to sustain rice production in all rice growing nations. The pathogen is widely distributed in all rice paddies and displays rapid aerial transmissions, and seed-borne latent infection. In order to understand the genetic variability, host specificity, and molecular basis of the pathogenicity-associated traits, the whole genome of rice infecting Magnaporthe oryzae (Strain RMg_Dl) was sequenced using the Illumina and PacBio (RSII compatible) platforms. The high-throughput hybrid assembly of short and long reads resulted in a total of 375 scaffolds with a genome size of 42.43 Mb. Furthermore, comparative genome analysis revealed 99% average nucleotide identity (ANI) with other oryzae genomes and 83% against M. grisea, and 73% against M. poe genomes. The gene calling identified 10,553 genes with 10,539 protein-coding sequences. Among the detected transposable elements, the LTR/Gypsy and Type LINE showed high occurrence. The InterProScan of predicted protein sequences revealed that 97% protein family (PFAM), 98% superfamily, and 95% CDD were shared among RMg_Dl and reference 70-15 genome, respectively. Additionally, 550 CAZymes with high GH family content/distribution and cell wall degrading enzymes (CWDE) such endoglucanase, beta-glucosidase, and pectate lyase were also deciphered in RMg_Dl. The prevalence of virulence factors determination revealed that 51 different VFs were found in the genome. The biochemical pathway such as starch and sucrose metabolism, mTOR signaling, cAMP signaling, MAPK signaling pathways related genes were identified in the genome. The 49,065 SNPs, 3267 insertions and 3611 deletions were detected, and majority of these varinats were located on downstream and upstream region. Taken together, the generated information will be useful to develop a specific marker for diagnosis, pathogen surveillance and tracking, molecular taxonomy, and species delineation which ultimately leads to device improved management strategies for blast disease.

First Draft Genome Sequence of a Pearl Millet Blast Pathogen, Magnaporthe grisea Strain PMg_Dl, Obtained Using PacBio Single-Molecule Real-Time and Illumina NextSeq 500 Sequencing.

The first draft genome sequence of the pearl millet blast pathogen Magnaporthe grisea PMg_Dl from India is presented. The genome information of M. grisea will be useful to understand the Magnaporthe speciation, genetic diversity, environmental adaptation, and pathogenic and host range determinants.

Genome Sequence of a Unique Magnaporthe oryzae RMg-Dl Isolate from India That Causes Blast Disease in Diverse Cereal Crops, Obtained Using PacBio Single-Molecule and Illumina HiSeq2500 Sequencing.

The whole-genome assembly of a unique rice isolate from India, Magnaporthe oryzae RMg-Dl that causes blast disease in diverse cereal crops is presented. Analysis of the 34.82 Mb genome sequence will aid in better understanding the genetic determinants of host range, host jump, survival, pathogenicity, and virulence factors of M. oryzae.