Integrated Analysis of Gene Expression, SNP, InDel, and CNV Identifies Candidate Avirulence Genes in Australian Isolates of the Wheat Leaf Rust Pathogen Puccinia triticina

Candidate effectors for leaf rust resistance gene Lr28 identified through transcriptome and in-silico analysis

Puccinia spp. causing rust diseases in wheat and other cereals secrete several specialized effector proteins into host cells. Characterization of these proteins and their interaction with host's R proteins could greatly help to limit crop losses due to diseases. Prediction of effector proteins by combining the transcriptome analysis and multiple in-silico approaches is gaining importance in revealing the pathogenic mechanism. The present study involved identification of 13 Puccinia triticina (Pt) coding sequences (CDSs), through transcriptome analysis, that were differentially expressed during wheat-leaf rust interaction; and prediction of their effector like features using different in-silico tools. NCBI-BLAST and pathogen-host interaction BLAST (PHI-BLAST) tools were used to annotate and classify these sequences based on their most closely matched counterpart in both the databases. Homology between CDSs and the annotated sequences in the NCBI database ranged from 79 to 94% and with putative effectors of other plant pathogens in PHI-BLAST from 24.46 to 54.35%. Nine of the 13 CDSs had effector-like features according to EffectorP 3.0 (≥0.546 probability of these sequences to be effector). The qRT-PCR expression analysis revealed that the relative expression of all CDSs in compatible interaction (HD2329) was maximum at 11 days post inoculation (dpi) and that in incompatible interactions (HD2329 + Lr28) was maximum at 3 dpi in seven and 9 dpi in five CDSs. These results suggest that six CDSs (>0.8 effector probability as per EffectorP 3.0) could be considered as putative Pt effectors. The molecular docking and MD simulation analysis of these six CDSs suggested that candidate Lr28 protein binds more strongly to candidate effector c14094_g1_i1 to form more stable complex than the remaining five. Further functional characterization of these six candidate effectors should prove useful for a better understanding of wheat-leaf rust interaction. In turn, this should facilitate effector-based leaf rust resistance breeding in wheat.

Integrated Analysis of Transcriptome and microRNA Profile Reveals the Toxicity of Euphorbia Factors toward Human Colon Adenocarcinoma Cell Line Caco-2

Euphorbia factors, lathyrane-type diterpenoids isolated from the medical herb Euphorbia lathyris L. (Euphorbiaceae), have been associated with intestinal irritation toxicity, but the mechanisms underlying this phenomenon are still unknown. The objective of this study was to evaluate the transcriptome and miRNA profiles of human colon adenocarcinoma Caco-2 cells in response to Euphorbia factors L1 (EFL1) and EFL2. Whole transcriptomes of mRNA and microRNA (miRNA) were obtained using second generation high-throughput sequencing technology in response to 200 μM EFL treatment for 72 h, and the differentially expressed genes and metabolism pathway were enriched. Gene structure changes were analyzed by comparing them with reference genome sequences. After 72 h of treatment, 16 miRNAs and 154 mRNAs were differently expressed between the EFL1 group and the control group, and 47 miRNAs and 1101 mRNAs were differentially expressed between the EFL2 group and the control. Using clusters of orthologous protein enrichment, the sequenced mRNAs were shown to be mainly involved in transcription, post-translational modification, protein turnover, chaperones, signal transduction mechanisms, intracellular trafficking, secretion, vesicular transport, and the cytoskeleton. The differentially expressed mRNA functions and pathways were enriched in transmembrane transport, T cell extravasation, the IL-17 signaling pathway, apoptosis, and the cell cycle. The differentially expressed miRNA EFLs caused changes in the structure of the gene, including alternative splicing, insertion and deletion, and single nucleotide polymorphisms. This study reveals the underlying mechanism responsible for the toxicity of EFLs in intestinal cells based on transcriptome and miRNA profiles of gene expression and structure.

Genotype Groups of the Wheat Leaf Rust Fungus Puccinia triticina in the United States as Determined by Genotyping by Sequencing

Wheat leaf rust caused by Puccinia triticina is a widespread disease of wheat in the United States and worldwide. Populations of P. triticina are characterized by virulence phenotypes that change rapidly because of selection by wheat cultivars with leaf rust resistance genes. The objective of this study was to genotype collections of P. triticina from 2011 to 2018 in the United States using restriction site-associated genotyping by sequencing (GBS) to determine if recently identified new virulence phenotypes belong to established genotype groups or to groups previously not detected. A total of 158 isolates were phenotyped for virulence on 20 lines of Thatcher wheat that are isogenic for leaf rust resistance genes and also genotyped for single nucleotide polymorphism. Eight distinct groups of P. triticina genotypes from common wheat were described based on coancestry, nucleotide divergence, and principal coordinate plots. A separate genotype group had isolates with virulence to durum wheat. Isolates within groups had similar virulence phenotypes, and the overall population had high levels of heterozygosity and high levels of linkage disequilibria, which were all indicators of clonality. Two new genotype groups were described, thereby raising the possibility of new introductions of P. triticina; however, genotypes in these groups may have also originated from somatic nuclear exchange and recombination. A genome-wide association study detected 19 single nucleotide polymorphisms that were highly associated with virulence to 11 resistance genes in the Thatcher near-isogenic lines.

Host Adaptation and Virulence in Heteroecious Rust Fungi

Rust fungi (Pucciniales, Basidiomycota) are obligate biotrophic pathogens that cause rust diseases in plants, inflicting severe damage to agricultural crops. Pucciniales possess the most complex life cycles known in fungi. These include an alternation of generations, the development of up to five different sporulating stages, and, for many species, the requirement of infecting two unrelated host plants during different parts of their life cycle, termed heteroecism. These fungi have been extensively studied in the past century through microscopy and inoculation studies, providing precise descriptions of their infection processes, although the molecular mechanisms underlying their unique biology are poorly understood. In this review, we cover recent genomic and life cycle transcriptomic studies in several heteroecious rust species, which provide insights into the genetic tool kits associated with host adaptation and virulence, opening new avenues for unraveling their unique evolution.

A Chromosome-Scale Assembly of the Wheat Leaf Rust Pathogen Puccinia triticina Provides Insights Into Structural Variations and Genetic Relationships With Haplotype Resolution

Despite the global economic importance of the wheat leaf rust pathogen Puccinia triticina (Pt), genomic resources for Pt are limited and chromosome-level assemblies of Pt are lacking. Here, we present a complete haplotype-resolved genome assembly at a chromosome-scale for Pt using the Australian pathotype 64-(6),(7),(10),11 (Pt64; North American race LBBQB) built upon the newly developed technologies of PacBio and Hi-C sequencing. PacBio reads with ∼200-fold coverage (29.8 Gb data) were assembled by Falcon and Falcon-unzip and subsequently scaffolded with Hi-C data using Falcon-phase and Proximo. This approach allowed us to construct 18 chromosome pseudomolecules ranging from 3.5 to 12.3 Mb in size for each haplotype of the dikaryotic genome of Pt64. Each haplotype had a total length of ∼147 Mb, scaffold N 50 of ∼9.4 Mb, and was ∼93% complete for BUSCOs. Each haplotype had ∼29,800 predicted genes, of which ∼2,000 were predicted as secreted proteins (SPs). The investigation of structural variants (SVs) between haplotypes A and B revealed that 10% of the total genome was spanned by SVs, highlighting variations previously undetected by short-read based assemblies. For the first time, the mating type (MAT) genes on each haplotype of Pt64 were identified, which showed that MAT loci a and b are located on two chromosomes (chromosomes 7 and 14), representing a tetrapolar type. Furthermore, the Pt64 assembly enabled haplotype-based evolutionary analyses for 21 Australian Pt isolates, which highlighted the importance of a haplotype resolved reference when inferring genetic relationships using whole genome SNPs. This Pt64 assembly at chromosome-scale with full phase information provides an invaluable resource for genomic and evolutionary research, which will accelerate the understanding of molecular mechanisms underlying Pt-wheat interactions and facilitate the development of durable resistance to leaf rust in wheat and sustainable control of rust disease.