1
|
Nallathambi P, Umamaheswari C, Reddy B, Aarthy B, Javed M, Ravikumar P, Watpade S, Kashyap PL, Boopalakrishnan G, Kumar S, Sharma A, Kumar A. Deciphering the Genomic Landscape and Virulence Mechanisms of the Wheat Powdery Mildew Pathogen Blumeria graminis f. sp. tritici Wtn1: Insights from Integrated Genome Assembly and Conidial Transcriptomics. J Fungi (Basel) 2024; 10:267. [PMID: 38667938 PMCID: PMC11051031 DOI: 10.3390/jof10040267] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2023] [Revised: 03/16/2024] [Accepted: 03/19/2024] [Indexed: 04/28/2024] Open
Abstract
A high-quality genome sequence from an Indian isolate of Blumeria graminis f. sp. tritici Wtn1, a persistent threat in wheat farming, was obtained using a hybrid method. The assembly of over 9.24 million DNA-sequence reads resulted in 93 contigs, totaling a 140.61 Mb genome size, potentially encoding 8480 genes. Notably, more than 73.80% of the genome, spanning approximately 102.14 Mb, comprises retro-elements, LTR elements, and P elements, influencing evolution and adaptation significantly. The phylogenomic analysis placed B. graminis f. sp. tritici Wtn1 in a distinct monocot-infecting clade. A total of 583 tRNA anticodon sequences were identified from the whole genome of the native virulent strain B. graminis f. sp. tritici, which comprises distinct genome features with high counts of tRNA anticodons for leucine (70), cysteine (61), alanine (58), and arginine (45), with only two stop codons (Opal and Ochre) present and the absence of the Amber stop codon. Comparative InterProScan analysis unveiled "shared and unique" proteins in B. graminis f. sp. tritici Wtn1. Identified were 7707 protein-encoding genes, annotated to different categories such as 805 effectors, 156 CAZymes, 6102 orthologous proteins, and 3180 distinct protein families (PFAMs). Among the effectors, genes like Avra10, Avrk1, Bcg-7, BEC1005, CSEP0105, CSEP0162, BEC1016, BEC1040, and HopI1 closely linked to pathogenesis and virulence were recognized. Transcriptome analysis highlighted abundant proteins associated with RNA processing and modification, post-translational modification, protein turnover, chaperones, and signal transduction. Examining the Environmental Information Processing Pathways in B. graminis f. sp. tritici Wtn1 revealed 393 genes across 33 signal transduction pathways. The key pathways included yeast MAPK signaling (53 genes), mTOR signaling (38 genes), PI3K-Akt signaling (23 genes), and AMPK signaling (21 genes). Additionally, pathways like FoxO, Phosphatidylinositol, the two-component system, and Ras signaling showed significant gene representation, each with 15-16 genes, key SNPs, and Indels in specific chromosomes highlighting their relevance to environmental responses and pathotype evolution. The SNP and InDel analysis resulted in about 3.56 million variants, including 3.45 million SNPs, 5050 insertions, and 5651 deletions within the whole genome of B. graminis f. sp. tritici Wtn1. These comprehensive genome and transcriptome datasets serve as crucial resources for understanding the pathogenicity, virulence effectors, retro-elements, and evolutionary origins of B. graminis f. sp. tritici Wtn1, aiding in developing robust strategies for the effective management of wheat powdery mildew.
Collapse
Affiliation(s)
- Perumal Nallathambi
- ICAR-Indian Agricultural Research Institute, Regional Station, Wellington 643231, Tamil Nadu, India; (P.N.); (C.U.); (B.A.); (P.R.)
| | - Chandrasekaran Umamaheswari
- ICAR-Indian Agricultural Research Institute, Regional Station, Wellington 643231, Tamil Nadu, India; (P.N.); (C.U.); (B.A.); (P.R.)
| | - Bhaskar Reddy
- ICAR-Indian Agricultural Research Institute, Pusa Campus, New Delhi 110012, Delhi, India; (M.J.); (G.B.)
| | - Balakrishnan Aarthy
- ICAR-Indian Agricultural Research Institute, Regional Station, Wellington 643231, Tamil Nadu, India; (P.N.); (C.U.); (B.A.); (P.R.)
| | - Mohammed Javed
- ICAR-Indian Agricultural Research Institute, Pusa Campus, New Delhi 110012, Delhi, India; (M.J.); (G.B.)
| | - Priya Ravikumar
- ICAR-Indian Agricultural Research Institute, Regional Station, Wellington 643231, Tamil Nadu, India; (P.N.); (C.U.); (B.A.); (P.R.)
| | - Santosh Watpade
- ICAR-Indian Agricultural Research Institute, Regional Station, Shimla 171004, Himachal Pradesh, India;
| | - Prem Lal Kashyap
- ICAR-Indian Institute of Wheat and Barley Research, Karnal 132001, Haryana, India; (P.L.K.); (S.K.); (A.S.)
| | | | - Sudheer Kumar
- ICAR-Indian Institute of Wheat and Barley Research, Karnal 132001, Haryana, India; (P.L.K.); (S.K.); (A.S.)
| | - Anju Sharma
- ICAR-Indian Institute of Wheat and Barley Research, Karnal 132001, Haryana, India; (P.L.K.); (S.K.); (A.S.)
| | - Aundy Kumar
- ICAR-Indian Agricultural Research Institute, Pusa Campus, New Delhi 110012, Delhi, India; (M.J.); (G.B.)
| |
Collapse
|
2
|
The First Telomere-to-Telomere Chromosome-Level Genome Assembly of Stagonospora tainanensis Causing Sugarcane Leaf Blight. J Fungi (Basel) 2022; 8:jof8101088. [PMID: 36294653 PMCID: PMC9605480 DOI: 10.3390/jof8101088] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2022] [Revised: 10/08/2022] [Accepted: 10/13/2022] [Indexed: 11/17/2022] Open
Abstract
The sexual morph Leptosphaeria taiwanensis Yen and Chi and its asexual morph Stagonospora tainanensis W. H. Hsieh is an important necrotrophic fungal phytopathogen, which causes sugarcane leaf blight, resulting in loss of cane tonnage and sucrose in susceptible sugarcane varieties. Decoding the genome and understanding of the basis of virulence is vitally important for devising effective disease control strategies. Here, we present a 38.25-Mb high-quality genome assembly of S. tainanensis strain StFZ01, denovo assembled with 10.19 Gb Nanopore sequencing long reads (~267×) and 3.82 Gb Illumina short reads (~100×). The genome assembly consists of 12 contigs with N50 of 2.86 Mb of which 5 belong to the telomere to telomere (T2T) chromosome. It contains 13.20% repeat sequences, 12,543 proteins, and 12,206 protein-coding genes with the BUSCO completeness 99.18% at fungi (n = 758) and 99.87% at ascomycota (n = 1706), indicating the high accuracy and completeness of our gene annotations. The virulence analysis in silico revealed the presence of 2379 PHIs, 599 CAZys, 248 membrane transport proteins, 191 cytochrome P450 enzymes, 609 putative secreted proteins, and 333 effectors in the StFZ01 genome. The genomic resources presented here will not only be helpful for development of specific molecular marker and diagnosis technique, population genetics, molecular taxonomy, and disease managements, it can also provide a significant precise genomic reference for investigating the ascomycetous genome, the necrotrophic lifestyle, and pathogenicity in the future.
Collapse
|