Comparative genomic analysis of monosporidial and monoteliosporic cultures for unraveling the complexity of molecular pathogenesis of Tilletia indica pathogen of wheat

Genome and transcriptome based comparative analysis of Tilletia indica to decipher the causal genes for pathogenicity of Karnal bunt in wheat

Tilletia indica Mitra causes Karnal bunt (KB) in wheat by pathogenic dikaryophase. The present study is the first to provide the draft genomes of the dikaryon (PSWKBGD-3) and its two monosporidial lines (PSWKBGH-1 and 2) using Illumina and PacBio reads, their annotation and the comparative analyses among the three genomes by extracting polymorphic SSR markers. The trancriptome from infected wheat grains of the susceptible wheat cultivar WL711 at 24 h, 48h, and 7d after inoculation of PSWKBGH-1, 2 and PSWKBGD-3 were also isolated. Further, two transcriptome analyses were performed utilizing T. indica transcriptome to extract dikaryon genes responsible for pathogenesis, and wheat transcriptome to extract wheat genes affected by dikaryon involved in plant-pathogen interaction during progression of KB in wheat. A total of 54, 529, and 87 genes at 24hai, 48hai, and 7dai, respectively were upregulated in dikaryon stage while 21, 35, and 134 genes of T. indica at 24hai, 48hai, and 7dai, respectively, were activated only in dikaryon stage. While, a total of 23, 17, and 52 wheat genes at 24hai, 48hai, and 7dai, respectively were upregulated due to the presence of dikaryon stage only. The results obtained during this study have been compiled in a web resource called TiGeR ( http://backlin.cabgrid.res.in/tiger/ ), which is the first genomic resource for T. indica cataloguing genes, genomic and polymorphic SSRs of the three T. indica lines, wheat and T. indica DEGs as well as wheat genes affected by T. indica dikaryon along with the pathogenecity related proteins of T. indica dikaryon during incidence of KB at different time points. The present study would be helpful to understand the role of dikaryon in plant-pathogen interaction during progression of KB, which would be helpful to manage KB in wheat, and to develop KB-resistant wheat varieties.

Draft genome sequencing of Tilletia caries inciting common bunt of wheat provides pathogenicity-related genes

Common bunt of wheat caused by Tilletia caries is an important disease worldwide. The T. caries TC1_MSG genome was sequenced using the Illumina HiSeq 2500 and Nanopore ONT platforms. The Nanopore library was prepared using the ligation sequencing kit SQK-LSK110 to generate approximately 24 GB for sequencing. The assembly size of 38.18 Mb was generated with a GC content of 56.10%. The whole genome shotgun project was deposited at DDBJ/ENA/GenBank under the accession number JALUTQ000000000. Forty-six contigs were obtained with N50 of 1,798,756 bp. In total, 10,698 genes were predicted in the assembled genome. Out of 10,698 genes, 10,255 genes were predicted significantly in the genome. The repeat sequences made up approximately 1.57% of the genome. Molecular function, cellular components, and biological processes for predicted genes were mapped into the genome. In addition, repeat elements in the genome were assessed. In all, 0.89% of retroelements were observed, followed by long terminal repeat elements (0.86%) in the genome. In simple sequence repeat (SSR) analysis, 8,582 SSRs were found in the genome assembly. The trinucleotide SSR type (3,703) was the most abundant. Few putative secretory signal peptides and pathogenicity-related genes were predicted. The genomic information of T. caries will be valuable in understanding the pathogenesis mechanism as well as developing new methods for the management of the common bunt disease of wheat.

Comparative genomics reveals low levels of inter- and intraspecies diversity in the causal agents of dwarf and common bunt of wheat and hint at conspecificity of Tilletia caries and T. laevis

Tilletia caries and T. laevis, which are the causal agents of common bunt, as well as T. controversa, which causes dwarf bunt of wheat, threaten especially organic wheat farming. The three closely related fungal species differ in their teliospore morphology and partially in their physiology and infection biology. The gene content as well as intraspecies variation in these species and the genetic basis of their separation is unknown. We sequenced the genome of four T. caries, five T. controversa, and two T. laevis and extended this dataset with five publicly available ones. The genomes of the three species displayed microsynteny with up to 94.3% pairwise aligned regions excluding repetitive regions. The majority of functionally characterized genes involved in pathogenicity, life cycle, and infection of corn smut, Ustilago maydis, were found to be absent or poorly conserved in the draft genomes and the biosynthetic pathway for trimethylamine in Tilletia spp. could be different from bacteria. Overall, 75% of the identified protein-coding genes comprising 84% of the total predicted carbohydrate utilizing enzymes, 72.5% putatively secreted proteins, and 47.4% of effector-like proteins were conserved and shared across all 16 isolates. We predicted nine highly identical secondary metabolite biosynthesis gene clusters comprising in total 62 genes in all species and none were species-specific. Less than 0.1% of the protein-coding genes were species-specific and their function remained mostly unknown. Tilletia controversa had the highest intraspecies genetic variation, followed by T. caries and the lowest in T. laevis. Although the genomes of the three species are very similar, employing 241 single copy genes T. controversa was phylogenetically distinct from T. caries and T. laevis, however these two could not be resolved as individual monophyletic groups. This was in line with the genome-wide number of single nucleotide polymorphisms and small insertions and deletions. Despite the conspicuously different teliospore ornamentation of T. caries and T. laevis, a high degree of genomic identity and scarcity of species-specific genes indicate that the two species could be conspecific.

Centenary of Soil and Air Borne Wheat Karnal Bunt Disease Research: A Review

Karnal bunt (KB) of wheat (Triticum aestivum L.), known as partial bunt has its origin in Karnal, India and is caused by Tilletia indica (Ti). Its incidence had grown drastically since late 1960s from northwestern India to northern India in early 1970s. It is a seed, air and soil borne pathogen mainly affecting common wheat, durum wheat, triticale and other related species. The seeds become inedible, inviable and infertile with the precedence of trimethylamine secreted by teliospores in the infected seeds. Initially the causal pathogen was named Tilletia indica but was later renamed Neovossia indica. The black powdered smelly spores remain viable for years in soil, wheat straw and farmyard manure as primary sources of inoculum. The losses reported were as high as 40% in India and also the cumulative reduction of national farm income in USA was USD 5.3 billion due to KB. The present review utilizes information from literature of the past 100 years, since 1909, to provide a comprehensive and updated understanding of KB, its causal pathogen, biology, epidemiology, pathogenesis, etc. Next generation sequencing (NGS) is gaining popularity in revolutionizing KB genomics for understanding and improving agronomic traits like yield, disease tolerance and disease resistance. Genetic resistance is the best way to manage KB, which may be achieved through detection of genes/quantitative trait loci (QTLs). The genome-wide association studies can be applied to reveal the association mapping panel for understanding and obtaining the KB resistance locus on the wheat genome, which can be crossed with elite wheat cultivars globally for a diverse wheat breeding program. The review discusses the current NGS-based genomic studies, assembly, annotations, resistant QTLs, GWAS, technology landscape of diagnostics and management of KB. The compiled exhaustive information can be beneficial to the wheat breeders for better understanding of incidence of disease in endeavor of quality production of the crop.

Refinement of Draft Genome Assemblies of Pigeonpea (Cajanus cajan)

Genome assembly of short reads from large plant genomes remains a challenge in computational biology despite major developments in next generation sequencing. Of late several draft assemblies have been reported in sequenced plant genomes. The reported draft genome assemblies of Cajanus cajan have different levels of genome completeness, a large number of repeats, gaps, and segmental duplications. Draft assemblies with portions of genome missing are shorter than the referenced original genome. These assemblies come with low map accuracy affecting further functional annotation and the prediction of gene components as desired by crop researchers. Genome coverage, i.e., the number of sequenced raw reads mapped onto a certain location of the genome is an important quality indicator of completeness and assembly quality in draft assemblies. The present work aimed to improve the coverage in reported de novo sequenced draft genomes (GCA_000340665.1 and GCA_000230855.2) of pigeonpea, a legume widely cultivated in India. The two recently sequenced assemblies, A1 and A2 comprised 72% and 75% of the estimated coverage of the genome, respectively. We employed an assembly reconciliation approach to compare the draft assemblies and merge them, filling the gaps by employing an algorithm size sorting mate-pair library to generate a high quality and near complete assembly with enhanced contiguity. The majority of gaps present within scaffolds were filled with right-sized mate-pair reads. The improved assembly reduced the number of gaps than those reported in draft assemblies resulting in an improved genome coverage of 82.4%. Map accuracy of the improved assembly was evaluated using various quality metrics and for the presence of specific trait-related functional genes. Employed pair-end and mate-pair local libraries helped us to reduce gaps, repeats, and other sequence errors resulting in lengthier scaffolds compared to the two draft assemblies. We reported the prediction of putative host resistance genes against Fusarium wilt disease by their performance and evaluated them both in wet laboratory and field phenotypic conditions.

Refinement of Draft Genome Assemblies of Pigeonpea (Cajanus cajan)

Genome assembly of short reads from large plant genomes remains a challenge in computational biology despite major developments in next generation sequencing. Of late several draft assemblies have been reported in sequenced plant genomes. The reported draft genome assemblies of Cajanus cajan have different levels of genome completeness, a large number of repeats, gaps, and segmental duplications. Draft assemblies with portions of genome missing are shorter than the referenced original genome. These assemblies come with low map accuracy affecting further functional annotation and the prediction of gene components as desired by crop researchers. Genome coverage, i.e., the number of sequenced raw reads mapped onto a certain location of the genome is an important quality indicator of completeness and assembly quality in draft assemblies. The present work aimed to improve the coverage in reported de novo sequenced draft genomes (GCA_000340665.1 and GCA_000230855.2) of pigeonpea, a legume widely cultivated in India. The two recently sequenced assemblies, A1 and A2 comprised 72% and 75% of the estimated coverage of the genome, respectively. We employed an assembly reconciliation approach to compare the draft assemblies and merge them, filling the gaps by employing an algorithm size sorting mate-pair library to generate a high quality and near complete assembly with enhanced contiguity. The majority of gaps present within scaffolds were filled with right-sized mate-pair reads. The improved assembly reduced the number of gaps than those reported in draft assemblies resulting in an improved genome coverage of 82.4%. Map accuracy of the improved assembly was evaluated using various quality metrics and for the presence of specific trait-related functional genes. Employed pair-end and mate-pair local libraries helped us to reduce gaps, repeats, and other sequence errors resulting in lengthier scaffolds compared to the two draft assemblies. We reported the prediction of putative host resistance genes against Fusarium wilt disease by their performance and evaluated them both in wet laboratory and field phenotypic conditions.