Genome sequencing and assembly of Indian golden silkmoth, Antheraea assamensis Helfer (Saturniidae, Lepidoptera).

Muga silkworm (Antheraea assamensis), one of the economically important wild silkmoths, is unique among saturniid silkmoths. It is confined to the North-eastern part of India. Muga silk has the highest value among the other silks. Unlike other silkmoths, A. assamensis has a low chromosome number (n = 15), and ZZ/ZO sex chromosome system. Here, we report the first high-quality draft genome of A. assamensis, assembled by employing the Illumina and PacBio sequencing platforms. The assembled genome of A. assamensis is 501.18 Mb long, with 2697 scaffolds and an N50 of 683.23 Kb. The genome encompasses 18,385 protein-coding genes, 86.29% of which were functionally annotated. Phylogenetic analysis of A. assamensis revealed its divergence from other Antheraea species approximately 28.7 million years ago. Moreover, an investigation into detoxification-related gene families, CYP450, GST, and ABC-transporter, revealed a significant expansion in A. assamensis as compared to the Bombyx mori. This expansion is comparable to Spodoptera litura, suggesting adaptive responses linked to the polyphagous behavior observed in these insects. This study provides valuable insights into the molecular basis of evolutionary divergence and adaptations in muga silkmoth. The genome assembly reported in this study will significantly help in the functional genomics studies on A. assamensis and other Antheraea species along with comparative genomics analyses of Bombycoidea insects.

Intrinsically Disordered Kiwellin Protein-Like Effectors Target Plant Chloroplasts and are Extensively Present in Rust Fungi.

The effector proteins produced by plant pathogens are one of the essential components of host-pathogen interaction. Despite being important, most of the effector proteins remain unexplored due to the diversity in their primary sequence generated by the high selection pressure of the host immune system. However to maintain the primary function in the infection process, these effectors may tend to maintain their native protein fold to perform the corresponding biological function. In the present study, unannotated candidate secretory effector proteins of sixteen major plant fungal pathogens were analyzed to find the conserved known protein folds using homology, ab initio, and Alpha Fold/Rosetta Fold protein dimensional (3D) structure approaches. Several unannotated candidate effector proteins were found to match various known conserved protein families potentially involved in host defense manipulation in different plant pathogens. Surprisingly a large number of plant Kiwellin proteins fold like secretory proteins (> 100) were found in studied rust fungal pathogens. Many of them were predicted as potential effector proteins. Furthermore, template independent modelling using Alpha Fold/Rosetta Fold analysis and structural comparison of these candidates also predicted them to match with plant Kiwellin proteins. We also found plant Kiwellin matching proteins outside rusts including several non-pathogenic fungi suggesting the broad function of these proteins. One of the highest confidently modeled Kiwellin matching candidates effectors, Pstr_13960 (97.8%), from the Indian P. striiformis race Yr9 was characterized using overexpression, localization, and deletion studies in Nicotiana benthamiana. The Pstr_13960 suppressed the BAX-induced cell death and localized in the chloroplast. Furthermore, the expression of the Kiwellin matching region (Pst_13960_kiwi) alone suppressed the BAX-induced cell death in N. benthamiana despite the change of location to the cytoplasm and nucleus, suggesting the novel function of the Kiwellin core fold in rust fungi. Molecular docking showed that Pstr_13960 can interact with plant Chorismate mutases (CMs) using three loops conserved in plant and rust Kiwellins. Further analysis of Pstr_13960 showed to contain Intrinsically disordered regions (IDRs) in place of the N-terminal ß1/ß2 region found in plant Kiwellins suggesting the evolution of rust Kiwellins-like effectors (KLEs). Overall, this study reports the presence of a Kiwellin protein-like fold containing a novel effector protein family in rust fungi depicting a classical example of the evolution of effectors at the structure level as Kiwellin effectors show very low significant similarity to plant Kiwellin at the sequence level.

Understanding the Dynamics of Blast Resistance in Rice-Magnaporthe oryzae Interactions.

Rice is a global food grain crop for more than one-third of the human population and a source for food and nutritional security. Rice production is subjected to various stresses; blast disease caused by Magnaporthe oryzae is one of the major biotic stresses that has the potential to destroy total crop under severe conditions. In the present review, we discuss the importance of rice and blast disease in the present and future global context, genomics and molecular biology of blast pathogen and rice, and the molecular interplay between rice-M. oryzae interaction governed by different gene interaction models. We also elaborated in detail on M. oryzae effector and Avr genes, and the role of noncoding RNAs in disease development. Further, rice blast resistance QTLs; resistance (R) genes; and alleles identified, cloned, and characterized are discussed. We also discuss the utilization of QTLs and R genes for blast resistance through conventional breeding and transgenic approaches. Finally, we review the demonstrated examples and potential applications of the latest genome-editing tools in understanding and managing blast disease in rice.

Risk Factors for Contracting Invasive Meningococcal Disease and Related Mortality: A Systematic Literature Review and Meta-analysis.

OBJECTIVES: To describe risk factors (RFs) and quantify their effects in invasive meningococcal disease (IMD) and associated mortality across all age groups based on the available published literature. METHODS: A systematic literature review (SLR) was conducted via MEDLINE® and Embase. Study selection, data extraction, and quality assessment were performed by two independent reviewers. Associations between RFs and outcomes were quantified via a meta-analysis (MA). RESULTS: Seventy-four studies (date range 1950 - 2018) were included in the SLR. Statistically significant RFs for contracting IMD identified from the SLR (within-study) included previous IMD infection and young age (0 - 4 years). MA indicated that significant RFs for contracting IMD (11 studies) were: HIV-positive status, passive smoke exposure, and crowded living space. In the MA for IMD-related mortality risk (11 studies), age 25 - 45 years (vs. 0 - 5 years) and serogroup C (vs. serogroup B) were significantly associated with increased risk. CONCLUSIONS: Previous findings of higher risk for IMD contraction with smoke exposure and crowded living conditions in children/adolescents have been extended by this SLR/MA to all age groups. We provide strong evidence for higher risk of IMD in HIV-positive individuals, and confirm previous findings of higher IMD-related mortality risk in adults aged 25 - 45.

Genetic characterisation of an Iflavirus associated with a vomiting disease in the Indian tropical tasar silkworm, Antheraea mylitta.

Antheraea mylitta, the Tropical tasar silkworm, is frequently affected by a vomiting disease called Virosis by sericulturists although not confirmed being of viral origin. Based on the symptoms and the disease pattern, the causal agent is however suspected to be a virus. The condition involves a series of characteristic and progressive symptoms that generally culminates in the death of the larva. The disease is common in autumn season (Sep- Oct), with widespread distribution causing severe damage to the tasar silk industry. The leads for this study were obtained from a transcript identified in the EST database in a different study, which matched the positive strand of Iflavirus, an RNA virus known to infect insects. In the present study the genome of this novel Iflavirus was characterised and the full length of the genome was found to be 9728 nucleotides long encoding for a single large open reading frame (ORF) with flanking NTR regions at 5' and 3' ends and a natural poly A tail at the 3' end. The ORF encoded structural proteins at the N-terminal end and non-structural proteins at the C-terminal end with a predicted 2967 amino acid long polyprotein. The structural proteins consisted of 4 proteins (VP1-VP4) and the non-structural proteins consisted of helicase, RNA-dependent RNA polymerase and 3C-protease. The virus is found in various tissues (midgut, fatbody, trachea, Malpighian tubules and silk gland) and also has a vertical route of transmission, i.e., from gravid females to the offspring. Based on the available data, the virus is a new member of Iflaviridae for which we propose the name Antheraea mylitta Iflavirus (AmIV).

Genome-Wide Analysis of Four Pathotypes of Wheat Rust Pathogen (Puccinia graminis) Reveals Structural Variations and Diversifying Selection.

Diseases caused by Puccinia graminis are some of the most devastating diseases of wheat. Extensive genomic understanding of the pathogen has proven helpful not only in understanding host- pathogen interaction but also in finding appropriate control measures. In the present study, whole-genome sequencing of four diverse P. graminis pathotypes was performed to understand the genetic variation and evolution. An average of 63.5 Gb of data per pathotype with about 100× average genomic coverage was achieved with 100-base paired-end sequencing performed with Illumina Hiseq 1000. Genome structural annotations collectively predicted 9273 functional proteins including ~583 extracellular secreted proteins. Approximately 7.4% of the genes showed similarity with the PHI database which is suggestive of their significance in pathogenesis. Genome-wide analysis demonstrated pathotype 117-6 as likely distinct and descended through a different lineage. The 3-6% more SNPs in the regulatory regions and 154 genes under positive selection with their orthologs and under negative selection in the other three pathotypes further supported pathotype 117-6 to be highly diverse in nature. The genomic information generated in the present study could serve as an important source for comparative genomic studies across the genus Puccinia and lead to better rust management in wheat.

Identification and expression analysis of pathogenicity-related genes of Rhizoctonia solani anastomosis groups infecting rice.

Sheath blight disease caused by Rhizoctonia solani Kuhn (teleomorph; Thanatephorus cucumeris) is a major constraint in rice production. Among the different anastomosis groups (AGs) of Rhizoctonia solani, AG1-IA causes sheath blight of rice, which induce necrotic lesions on leaf sheaths of the infected plants. Several reports contradict the host specificity of anastomosis groups in Rhizoctonia solani. There is lack of information on the pathogenicity genes of these Rhizoctonia solani anastomosis groups during sheath blight infection in rice. In the present study, Rhizoctonia solani isolates collected from diverse rice growing regions of India were screened for anastomosis groups and two groups namely, AG1-IA, AG2-2 were identified. Accordingly, comparative studies were made with AG1-IA (GenBank ID: 16,395) and AG2-2 (GenBank ID: 2,318,768) group sequences, which enabled the identification of specific gene clusters (119 in AG1-IA and 604 in AG2-2) belonging to these groups. Pathogen Host Interaction (PHI) blast with these specific gene clusters could further identify genes involved in host pathogen interaction (38 in AG1_IA and 150 in AG2-2), which were shortlisted for qRT-PCR validation based on qcov cutoff values representing different phenotypic categories of PHI blast. Expression analysis-based validation in sheath blight susceptible (Pusa Basmati 1) and resistant (Pusa 1908-13-12-5) rice genotypes showed that most of the genes expressed significantly higher in the susceptible variety Pusa Basmati 1. The genes like inorganic phosphate transporter (AG1_IPT), Bromodomain containing protein (AG1_BrD), Aldehyde dehydrogenase (AG1_AldD), AMP binding domain (AG1_AMP) and Heme peroxidase (AG1_HmPr) were upregulated in the susceptible genotype, PB 1 at 72hpi compared to Pusa 1908-13-12-5. Among these, inorganic phosphate transporter (AG1_IPT), Bromodomain containing protein (AG1_BrD) and Heme peroxidase (AG1_HmPr) were specific to Rhizoctonia solani AG1-IA. Through the present study, we could demonstrate the AG1-IA-specific interactions of Rhizoctonia solani causing sheath blight disease of rice, which is a step forward in understanding the specificity of Rhizoctonia solani with reference to sheath blight disease of rice. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s13205-021-02934-1.

Development and optimization of a TaqMan assay for Nosema bombycis, causative agent of pébrine disease in Bombyx mori silkworm, based on the ß-tubulin gene.

"Pébrine" is a devastating disease of Bombyx mori silkworms that is highly contagious and can completely destroy an entire crop of silkworms and is thus a serious threat for the viability and profitability of sericulture. The disease is most commonly attributed to microsporidians of the genus Nosema, which are obligate intracellular parasites that are transmitted through spores. Nosema infections in silkworms are diagnosed primarily through light microscopy, which is labour intensive and less reliable, sensitive, and specific than PCR-based techniques. Here, we present the development and optimization of a new TaqMan based assay targeting the ß-tubulin gene in the pébrine disease causing agent Nosema bombycis in silkworms. The assay displayed excellent quantification linearity over multiple orders of magnitude of target amounts and a limit of detection (LOD) of 6.9 × 102 copies of target per reaction. The method is highly specific to N. bombycis with no cross-reactivity to other Nosema species commonly infecting wild silkworms. This specificity was due to three nucleotides in the probe-binding region unique to N. bombycis. The assay demonstrated a high reliability with a Coefficient of variation (CV) <5% for both intra-assay and inter-assay variability. The assay was used to trace experimental N. bombycis infection of silkworm larvae, in the fat body, midgut and ovary tissues, through pupation and metamorphosis to the emerging female moth, and her larval off-spring, confirming the vertical transmission of N. bombycis in silkworms. The TaqMan assay revealed a gradual increase in infection levels in the post-infection samples. The assay is reliable and simple to implement and can be a suitable complement to microscopy for routine diagnostics and surveillance in silkworm egg production centres with appropriate infrastructure.

Deciphering the role of microRNAs during Pi54 gene mediated Magnaporthe oryzae resistance response in rice.

The broad-spectrum resistance gene Pi54 confers resistance to multiple isolates of Magnaporthe oryzae in rice. In order to decipher the molecular mechanism underlying the Pi54 mediated resistance in rice line Taipei309 Pi54 (carrying Pi54), miRNAome study was performed at 24 h post-inoculation (hpi) with M. oryzae. A total of 222 known miRNAs representing 101 miRNA families were found in this study. Of these, 29 and 24 miRNAs were respectively up- and down-regulated in the resistant Taipei309 Pi54 . Defence response (DR) genes, like, NBSGO35, and OsWAK129b, and genes related to transcription factors were up-regulated in Taipei309 Pi54 line. The vast array of miRNA candidates identified here are miR159c, miR167c, miR2100, miR2118o, miR2118l, miR319a, miR393, miR395l, miR397a, miR397b, miR398, miR439g, miR531b, miR812f, and miR815c, and they manifest their role in balancing the interplay between various DR genes during Pi54 mediated resistance. We also validated miRNA/target gene pairs involved in hormone signalling, and cross-talk among hormone pathways regulating the rice immunity. This study suggests that the Pi54 gene mediated blast resistance is influenced by several microRNAs through PTI and ETI components in the rice line Taipei309 Pi54 , leading to incompatible host-pathogen interaction.

Identification and functional prediction of long non-coding RNAs of rice (Oryza sativa L.) at reproductive stage under salinity stress.

Salinity adversely affects the yield and growth of rice (Oryza sativa L.) plants severely, particularly at reproductive stage. Long non-coding RNAs (lncRNAs) are key regulators of diverse molecular and cellular processes in plants. Till now, no systematic study has been reported for regulatory roles of lncRNAs in rice under salinity at reproductive stage. In this study, total 80 RNA-seq data of Horkuch (salt-tolerant) and IR-29 (salt-sensitive) genotypes of rice were used and found 1626 and 2208 transcripts as putative high confidence lncRNAs, among which 1529 and 2103 were found to be novel putative lncRNAs in root and leaf tissue respectively. In Horkuch and IR-29, 14 and 16 lncRNAs were differentially expressed in root tissue while 18 and 63 lncRNAs were differentially expressed in leaf tissue. Interaction analysis among the lncRNAs, miRNAs and corresponding mRNAs indicated that these modules are involved in different biochemical pathways e.g. phenyl propanoid pathway during salinity stress in rice. Interestingly, two differentially expressed lncRNAs such as TCONS_00008914 and TCONS_00008749 were found as putative target mimics of known rice miRNAs. This study indicates that lncRNAs are involved in salinity adaptation of rice at reproductive stage through certain biochemical pathways.

Transcriptome analysis of Nosema assamensis infecting muga silkworms (Antheraea assamensis) reveals insights into candidate pathogenicity related genes and molecular pathways required for pathogenesis.

Muga silkworms are often prone to many diseases since, these are non-domesticated and are reared outdoors. Microsporidia, an obligate intracellular pathogen with spore as its active form, causes pebrine disease in these silkworms. The study has attempted to categorise the transcript data of the Nosema obtained from the infected muga silkworm using gene ontology and KEGG pathway studies. A total of 2850 unigene sets were identified out of which 2739 unigenes were placed under biological, cellular as well as molecular function categories based on the gene ontology (GO) terms. 1620 out of these unigenes sets found their orthologous partner in the corresponding Nosema bombycis transcriptome. The unigenes were found to be enriched under organic substance metabolic process, organic cyclic compound binding and intracellular anatomical structure for biological process, molecular function and cellular components respectively. The KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analysis indicated majority of the enzymes were found to be mapped under purine and thiamine metabolic pathways, indicating an increase in the energy metabolism required to establish the infection in the silkworms. The putative virulence genes identified in this study are PTP2, PTP3, SWP12 and SWP26 which were found to be expressed in other Nosema species indigenous to India, indicating a probable conservation of these genes, which are primarily involved in establishing host pathogen interactions. The expression of these genes was in detectable levels in the infected silkworm samples. These genes may be validated further through bioassay in order understand their roles in establishing the infection and propagation of the spores. The identified virulence genes may be further targeted to develop diagnostic tools for identification of the pathogen at early stages of infection.

Effector Biology of Biotrophic Plant Fungal Pathogens: Current Advances and Future Prospects.

The interaction of fungal pathogens with their host requires a novel invading mechanism and the presence of various virulence-associated components responsible for promoting the infection. The small secretory proteins, explicitly known as effector proteins, are one of the prime mechanisms of host manipulation utilized by the pathogen to disarm the host. Several effector proteins are known to translocate from fungus to the plant cell for host manipulation. Many fungal effectors have been identified using genomic, transcriptomic, and bioinformatics approaches. Most of the effector proteins are devoid of any conserved signatures, and their prediction based on sequence homology is very challenging, therefore by combining the sequence consensus based upon machine learning features, multiple tools have also been developed for predicting apoplastic and cytoplasmic effectors. Various post-genomics approaches like transcriptomics of virulent isolates have also been utilized for identifying active consortia of effectors. Significant progress has been made in understanding biotrophic effectors; however, most of it is underway due to their complex interaction with host and complicated recognition and signaling networks. This review discusses advances, and challenges in effector identification and highlighted various features of the potential effector proteins and approaches for understanding their genetics and strategies for regulation.

Smut fungi as a stratagem to characterize rust effectors: opportunities and challenges.

The rust pathogens are one of the most complex fungi in the Basidiomycetes. The development of genomic resources for rust and other plant pathogens has opened the opportunities for functional genomics of fungal genes. Despite significant progress in the field of fungal genomics, functional characterization of the genome components has lacked, especially for the rust pathogens. Their obligate nature and lack of standard stable transformation protocol are the primary reasons for rusts to be one of the least explored genera despite its significance. In the recently sequenced rust genomes, a vast catalogue of predicted effectors and pathogenicity genes have been reported. However, most of these candidate genes remained unexplored due to the lack of suitable characterization methods. The heterologous expression of putative effectors in Nicotiana benthamiana and Arabidopsis thaliana has proved to be a rapid screening method for identifying the role of these effectors in virulence. However, no fungal system has been used for the functional validation of these candidate genes. The smuts, from the evolutionary point of view, are closely related to the rust pathogens. Moreover, they have been widely studied and hence could be a suitable model system for expressing rust fungal genes heterologously. The genetic manipulation methods for smuts are also well standardized. Complementation assays can be used for functional validation of the homologous genes present in rust and smut fungal pathogens, while the species-specific proteins can be expressed in the mutant strains of smut pathogens having reduced or no virulence for virulence analysis. We propose that smuts, especially Ustilago maydis, may prove to be a good model system to characterize rust effector proteins in the absence of methods to manipulate the rust genomes directly.

Identification and characterization of Dicer-like genes in leaf rust pathogen (Puccinia triticina) of wheat.

Puccinia triticina (P. triticina) is one of the most devastating fungal pathogens of wheat which causes significant annual yield loss to the crop. Understanding the gene regulatory mechanism of the biotrophic pathogen is one of the important aspects of host-pathogen interaction studies. Dicer-like genes are considered as important mediators of RNAi-based gene regulation. In this study, we report the presence of three Dicer-like genes (Pt-DCL1, Pt-DCL2, Pt-DCL3) in P. triticina genome identified through computational and biological analyses. Quantitative real-time PCR studies revealed an increase in the expression of these genes in germinating spore stages. Heterologous expression combined with mass spectrometry analysis of Pt-DCL2 confirmed the presence of a canonical Dicer-like gene in P. triticina. Phylogenetic analysis of the Pt-DCLs with the Dicer-like proteins from other organisms showed a distinct cluster of rust pathogens from the order Pucciniales. The results indicated a species-specific duplication of Dicer-like genes within the wheat rust pathogens. This study, for the first time, reports the presence of Dicer-dependent RNAi pathway in P. triticina that may play a role in gene regulatory mechanism of the pathogen during its development. Our study serves as a vital source of information for further RNAi-based molecular studies for better understanding and management of the wheat leaf rust disease.

Genome wide characterization revealed MnMLO2 and MnMLO6A as candidate genes involved in powdery mildew susceptibility in mulberry.

Mulberry is a fast growing economically important tree for sericulture industry and contains compounds for preventing and treating several diseases and ailments. The quality and quantity of mulberry leaf available to produce silk fibre and for medicinal purpose is greatly affected by number of foliar diseases, out of which powdery mildew is the major one. Imparting genetic resistance becomes an important approach in disease management in mulberry as spraying of fungicides has harmful effects on silkworm growth and development. Deployment of non-functional susceptible genes such as Mildew resistance Locus O (MLO) against powdery mildew in few crops stimulated to identify and characterize MLO genes in mulberry. In this study, genome wide analysis identified 16 MLO genes in Morus notabilis. Phylogenetic analysis found that MnMLO2, MnMLO6A, MnMLO6B, MnMLO12A and MnMLO12B clustered with functionally characterized MLOs associated with powdery mildew susceptibility in dicot species. Gene expression analysis indicated increased transcript abundance of MnMLO2, MnMLO6A, and MnMLO12A in response to powdery mildew infection. Further, conserved motifs exclusive to functionally characterized MLOs were identified in MnMLO1C, MnMLO2 and MnMLO6A proteins. Combined analysis of the phylogenetic relationship, conserved motif analysis and gene expression in response to infection identified MnMLO2 and MnMLO6A as potential candidate genes involved in powdery mildew susceptibility in mulberry. Identification and deployment of natural and induced mutations in the candidate genes can be useful for mulberry breeding programs to develop powdery mildew resistant varieties.

TeaMiD: a comprehensive database of simple sequence repeat markers of tea.

Tea is a highly cross-pollinated, woody, perennial tree. High heterozygosity combined with a long gestational period makes conventional breeding a cumbersome process. Therefore, marker-assisted breeding is a better alternative approach when compared with conventional breeding. Considering the large genome size of tea (~3 Gb), information about simple sequence repeat (SSR) is scanty. Thus, we have taken advantage of the recently published tea genomes to identify large numbers of SSR markers in the tea. Besides the genomic sequences, we identified SSRs from the other publicly available sequences such as RNA-seq, GSS, ESTs and organelle genomes (chloroplasts and mitochondrial) and also searched published literature to catalog validated set of tea SSR markers. The complete exercise yielded a total of 935 547 SSRs. Out of the total, 82 SSRs were selected for validation among a diverse set of tea genotypes. Six primers (each with four to six alleles, an average of five alleles per locus) out of the total 27 polymorphic primers were used for a diversity analysis in 36 tea genotypes with mean polymorphic information content of 0.61-0.76. Finally, using all the information generated in this study, we have developed a user-friendly database (TeaMiD; http://indianteagenome.in:8080/teamid/) that hosts SSR from all the six resources including three nuclear genomes of tea and transcriptome sequences of 17 Camellia wild species. Database URL: http://indianteagenome.in:8080/teamid/.

Deciphering signalling network in broad spectrum Near Isogenic Lines of rice resistant to Magnaporthe oryzae.

Disease resistance (R) genes like Pi9, Pita, Pi21, Pi54 are playing important role for broad spectrum blast resistance in rice. Development of near isogenic lines (NILs) using these type of broad spectrum genes and understanding their signalling networks is essential to cope up with highly evolving Magnaporthe oryzae strains for longer duration. Here, transcriptional-level changes were studied in three near-isogenic lines (PB1 + Pi1, PB1 + Pi9 and PB1 + Pi54) of rice resistant to blast infection, to find the loci that are unique to resistant lines developed in the background of Pusa Basmati 1 (PB1). The pathway analysis of loci, unique to resistant NILs compared to susceptible control revealed that plant secondary metabolite synthesis was the common mechanism among all NILs to counter against M. oryzae infection. Comparative transcriptome analysis helped to find out common clusters of co-expressed significant differentially expressed loci (SDEL) in both PB1 + Pi9 and PB1 + Pi54 NILs. SDELs from these clusters were involved in the synthesis and degradation of starch; synthesis and elongation of fatty acids; hydrolysis of phospholipids; synthesis of phenylpropanoid; and metabolism of ethylene and jasmonic acid. Through detailed analysis of loci specific to each resistant NIL, we identified a network of signalling pathways mediated by each blast resistance gene. The study also offers insights into transcriptomic dynamics, points to a set of important candidate genes that serve as module to regulate the changes in resistant NILs. We suggest that pyramiding of the blast resistance gene Pi9 with Pi54 will lead to maximum broad spectrum resistance to M. oryzae.

De novo genome sequencing and secretome analysis of Tilletia indica inciting Karnal bunt of wheat provides pathogenesis-related genes.

Tilletia indica is an internationally quarantined fungal pathogen causing Karnal bunt of wheat. The present study carried out that the whole genome of T. indica was sequenced and identified transposable elements, pathogenicity-related genes using a comparative genomics approach. The T. indica genome assembly size of 33.7 MB was generated using Illumina and Pac Bio platforms with GC content of 55.0%. A total of 1737 scaffolds were obtained with N50 of 58,667 bp. The ab initio gene prediction was performed using Ustilago maydis as the reference species. A total number of 10,113 genes were predicted with an average gene size of 1945 bp out of which functionally annotated genes were 7262. A total number of 3216 protein-coding genes were assigned in different categories. Out of a total number of 1877 transposable elements, gypsy had the highest count (573). Total 5772 simple sequence repeats were identified in the genome assembly, and the most abundant simple sequence repeat type was trinucleotide having 42% of total SSRs. The comparative genome analysis suggested 3751 proteins of T. indica had orthologs in five fungi, whereas 126 proteins were unique to T. indica. Secretome analysis revealed the presence of 1014 secretory proteins and few carbohydrate-active enzymes in the genome. Some putative candidate pathogenicity-related genes were identified in the genome. The whole genome of T. indica will provide a window to understand the pathogenesis mechanism, fungal life cycle, survival of teliospores, and novel strategies for management of Karnal bunt disease of wheat.

Transcriptome sequencing of Himalayan Raspberry (Rubus ellipticus) and development of simple sequence repeat markers.

KEY MESSAGE: Rubus ellipticus is a wild crop having less number of EST-SSR markers. First attempt was made towards the transcriptomics data analysis and generation of EST-SSR markers which were used in genetic diversity and transferability studies. ABSTRACT: Rubus ellipticus is a raspberry with yellow fruits, native to tropical and subtropical India and Asia. Leaves of Rubus ellipticus 'Kumarhatti' collection were utilized for cDNA library construction. More than 15 million sequencing reads were generated using NextSeq 500 Illumina RNA-seq technology. The DNASTAR software was used for de novo assembly from which 7777 unigenes with an average length of 500 bp was obtained. These unigenes were annotated using public databases, including NCBI non-redundant and gene ontology. De novo assembly of R. ellipticus unigenes found the highest similarity to apple than to other members of Rosaceae. This is the first attempt to use the Illumina platform of RNA sequencing and de novo assembly for R. ellipticus without a reference genome. In this study, unigenes were used for SSR marker development. ESTs containing SSR motifs were extracted using an online Microsatellite Identification Tool (MISA). SSR primers were designed from the SSR containing 704 EST sequences using the Websat software. Total 304 EST-SSRs primers were successfully designed, out of which 68 randomly selected primer pairs were custom synthesized and used for validation. Real-time PCR was also performed to analyze the relationship of transcriptional factors with fruit ripening. Out of 68 primer pairs, 61 were found to be informative in R. ellipticus, whereas 65 primer pairs were informative in the five tested genera of Rosaceae, i.e., pear, peach, apple, rose, and strawberry with 95.3% and 93.5% polymorphism, leading to the conclusion that these marker systems are very efficient to carryout diversity and cross transferability study in Rosaceae genera.

Discovery and profiling of small RNAs from Puccinia triticina by deep sequencing and identification of their potential targets in wheat.

Cross-kingdom RNAi is a well-documented phenomenon where sRNAs generated by host and pathogens may govern resistance or susceptible phenotypes during host-pathogen interaction. With the first example of the direct involvement of fungal generated sRNAs in virulence of plant pathogenic fungi Botrytis cinerea and recently from Puccinia striiformis f. sp. tritici, we attempted to identify sRNAs in Puccinia triticina (P. triticina). Four sRNA libraries were prepared and sequenced using Illumina sequencing technology and a total of ~ 1-1.28 million potential sRNAs and two microRNA-like small RNA (mil-RNAs) candidates were identified. Computational prediction of targets using a common set of sRNAs and P. triticina mil-RNAs (pt-mil-RNAs) within P. triticina and wheat revealed the majority of the targets as repetitive elements in P. triticina whereas in wheat, the target genes were identified to be involved in many biological processes including defense-related pathways. We found 9 receptor-like kinases (RLKs) and 14 target genes of each related to reactive oxygen species (ROS) pathway and transcription factors respectively, including significant numbers of target genes from various other categories. Expression analysis of twenty selected sRNAs, targeting host genes pertaining to ROS related, disease resistance, metabolic processes, transporter, apoptotic inhibitor, and transcription factors along with two pt-mil-RNAs by qRT-PCR showed distinct patterns of expression of the sRNAs in urediniospore-specific libraries. In this study, for the first time, we report identification of novel sRNAs identified in P. triticina including two pt-mil-RNAs that may play an important role in biotrophic growth and pathogenicity.