Search

Abstract Background Adelphocoris suturalis (Hemiptera: Miridae) is a notorious agricultural pest, which causes serious economic losses to a diverse range of agricultural crops around the world. The poor understanding of its genomic characteristics has seriously hindered the establishment of sustainable and environment-friendly agricultural pest management through biotechnology and biological insecticides. Results Here, we report a chromosome-level assembled genome of A. suturalis by integrating Illumina short reads, PacBio, 10x Chromium, and Hi-C mapping technologies. The resulting 1.29 Gb assembly contains twelve chromosomal pseudomolecules with an N50 of 1.4 and 120.6 Mb for the contigs and scaffolds, respectively, and carries 20,010 protein-coding genes. The considerable size of the A. suturalis genome is predominantly attributed to a high amount of retrotransposons, especially long interspersed nuclear elements (LINEs). Transcriptomic and phylogenetic analyses suggest that A. suturalis-specific candidate effectors, and expansion and expression of gene families associated with omnivory, insecticide resistance and reproductive characteristics, such as digestion, detoxification, chemosensory receptors and long-distance migration likely contribute to its strong environmental adaptability and ability to damage crops. Additionally, 19 highly credible effector candidates were identified and transiently overexpressed in Nicotiana benthamiana for functional assays and potential targeting for insect resistance genetic engineering. Conclusions The high-quality genome of A. suturalis provides an important genomic landscape for further investigations into the mechanisms of omnivory, insecticide resistance and survival adaptation, and for the development of integrated management strategies.

Background: Adelphocoris suturalis (Hemiptera: Miridae) is a notorious agricultural pest, which causes serious economic losses to a diverse range of agricultural crops around the world. The poor understanding of its genomic characteristics has seriously hindered the establishment of sustainable and environment-friendly agricultural pest management through biotechnology and biological insecticides. Results: Here, we report a chromosome-level assembled genome of A. suturalis by integrating Illumina short reads, PacBio, 10x Chromium, and Hi-C mapping technologies. The resulting 1.29 Gb assembly contains twelve chromosomal pseudomolecules with an N50 of 1.4 and 120.6 Mb for the contigs and scaffolds, respectively, and carries 20,010 protein-coding genes. The considerable size of the A. suturalis genome is predominantly attributed to a high amount of retrotransposons, especially long interspersed nuclear elements (LINEs). Transcriptomic and phylogenetic analyses suggest that A. suturalis-specific candidate effectors, and expansion and expression of gene families associated with omnivory, insecticide resistance and reproductive characteristics, such as digestion, detoxification, chemosensory receptors and long-distance migration likely contribute to its strong environmental adaptability and ability to damage crops. Additionally, 19 highly credible effector candidates were identified and transiently overexpressed in Nicotiana benthamiana for functional assays and potential targeting for insect resistance genetic engineering. Conclusions: The high-quality genome of A. suturalis provides an important genomic landscape for further investigations into the mechanisms of omnivory, insecticide resistance and survival adaptation, and for the development of integrated management strategies. [ABSTRACT FROM AUTHOR]

Sporisorium scitamineum is a biotrophic fungus responsible for sugarcane smut disease. To investigate the key genes involved in S. scitamineum infection, we conducted RNA sequencing of sugarcane sprouts inoculated with S. scitamineum teliospores. A weighted gene co-expression network analysis (WGCNA) showed that two co-expressed gene modules, MEdarkturquoise and MEpurple—containing 66 and 208 genes, respectively—were associated with S. scitamineum infection. The genes in these two modules were further studied using Gene Ontology (GO) enrichment analysis, pathogen-host interaction (PHI) database BLASTp, and small secreted cysteine-rich proteins (SCRPs) prediction. The top ten hub genes in each module were identified using the Cytohubba plugin. The GO enrichment analysis found that endoplasmic reticulum-related and catabolism-related genes were expressed during S. scitamineum infection. A total of 83 genes had homologs in the PHI database, 62 of which correlated with pathogen virulence. A total of 21 proteins had the characteristics of small secreted cysteine-rich proteins (SCRPs), a common source of fungal effectors. The top ten hub genes in each module were identified, and seven were annotated as Mig1-Mig1 protein, glycosyl hydrolase, beta-N-acetylglucosaminidase, secreted chorismate mutase, collagen, mRNA export factor, and pleckstrin homology domain protein, while the remaining three were unknown. Two SCRPs—SPSC_06609 and SPSC_04676—and three proteins—SPSC_01958, SPSC_02155, and SPSC_00940—identified in the PHI database were also among the top ten hub genes in the MEdarkturquoise and MEpurple modules, suggesting that they may play important roles in S. scitamineum infection. A S. scitamineum infection model was postulated based on current findings. These findings help to deepen the current understanding of early events in S. scitamineum infection. [ABSTRACT FROM AUTHOR]

Abstract Background The pine wood nematode (PWN), Bursaphelenchus xylophilus, is a devastating pathogen of many Pinus species in China. The aim of this study was to understand the interactive molecular mechanism of PWN and its host by comparing differentially expressed genes and candidate effectors from three transcriptomes of B. xylophilus at different infection stages. Results In total, 62, 69 and 46 candidate effectors were identified in three transcriptomes (2.5 h postinfection, 6, 12 and 24 h postinoculation and 6 and 15 d postinfection, respectively). In addition to uncharacterized pioneers, other candidate effectors were involved in the degradation of host tissues, suppression of host defenses, targeting plant signaling pathways, feeding and detoxification, which helped B. xylophilus survive successfully in the host. Seven candidate effectors were identified in both our study and the B. xylophilus transcriptome at 2.5 h postinfection, and one candidate effector was identified in all three transcriptomes. These common candidate effectors were upregulated at infection stages, and one of them suppressed pathogen-associated molecular pattern (PAMP) PsXEG1-triggered cell death in Nicotiana benthamiana. Conclusions The results indicated that B. xylophilus secreted various candidate effectors, and some of them continued to function throughout all infection stages. These various candidate effectors were important to B. xylophilus infection and survival, and they functioned in different ways (such as breaking down host cell walls, suppressing host defenses, promoting feeding efficiency, promoting detoxification and playing virulence functions). The present results provide valuable resources for in-depth research on the pathogenesis of B. xylophilus from the perspective of effectors.

Fusarium blight of wheat is usually caused by Fusarium graminearum, and the pathogenic fungi will secrete effectors into the host plant tissue to affect its normal physiological process, so as to make it pathogenic. The CFEM (Common in Fungal Extracellular Membrane) protein domain is unique to fungi, but it is not found in all fungi. The CFEM protein contained in F. graminearum may be closely related to pathogenicity. In this study, 23 FgCFEM proteins were identified from the F. graminearum genome. Then, features of these proteins, such as signal peptide, subcellular localization, and transmembrane domains, etc., were analyzed and candidate effectors were screened out. Sequence alignment results revealed that each FgCFEM protein contains one CFEM domain. The amino acids of the CFEM domain are highly conserved and contain eight spaced cysteines, with the exception that FgCFEM8, 9, and 15 lack two cysteines and three cysteines were missed in FgCFEM18 and FgCFEM22. A recently identified CFEM_DR motif was detected in 11 FgCFEMs, and importantly we identified two new conserved motifs containing about 29 and 18 amino acids (CFEM_WR and CFEM_KF), respectively, in some of FgCFEM proteins. Transcriptome analysis of the genes encoding CFEM proteins indicated that all the CFEM-containing genes were expressed during wheat infection, with seven and six genes significantly up- and down-regulated, respectively, compared with in planta and in vitro. Based on the above analysis, FgCFEM11 and FgCFEM23 were predicted to be F. graminearum effectors. This study provides the basis for future functional analyses of CFEM proteins in F. graminearum. [ABSTRACT FROM AUTHOR]

Abstract Background Gymnosporangium spp. are fungal plant pathogens causing rust disease and most of them are known to infect two different host plants (heteroecious) with four spore stages (demicyclic). In the present study, we sequenced the transcriptome of G. japonicum teliospores on its host plant Juniperus chinensis and we performed comparison to the transcriptomes of G. yamadae and G. asiaticum at the same life stage, that happens in the same host but on different organs. Results Functional annotation for the three Gymnosporangium species showed the expression of a conserved genetic program with the top abundant cellular categories corresponding to energy, translation and signal transduction processes, indicating that this life stage is particularly active. Moreover, the survey of predicted secretomes in the three Gymnosporangium transcriptomes revealed shared and specific genes encoding carbohydrate active enzymes and secreted proteins of unknown function that could represent candidate pathogenesis effectors. A transcript encoding a hemicellulase of the glycoside hydrolase 26 family, previously identified in other rust fungi, was particularly highly expressed suggesting a general role in rust fungi. The comparison between the transcriptomes of the three Gymnosporangium spp. and selected Pucciniales species in different taxonomical families allowed to identify lineage-specific protein families that may relate to the biology of teliospores in rust fungi. Among clustered gene families, 205, 200 and 152 proteins were specifically identified in G. japonicum, G. yamadae and G. asiaticum, respectively, including candidate effectors expressed in teliospores. Conclusions This comprehensive comparative transcriptomics study of three Gymnosporangium spp. identified gene functions and metabolic pathways particularly expressed in teliospores, a stage of the life cycle that is mostly overlooked in rust fungi. Secreted protein encoding transcripts expressed in teliospores may reveal new candidate effectors related to pathogenesis. Although this spore stage is not involved in host plant infection but in the production of basidiospores infecting plants in the Amygdaloideae, we speculate that candidate effectors may be expressed as early as the teliospore stage for preparing further infection by basidiospores.

Sporisorium scitamineum is a biotrophic fungus responsible for sugarcane smut disease. To investigate the key genes involved in S. scitamineum infection, we conducted RNA sequencing of sugarcane sprouts inoculated with S. scitamineum teliospores. A weighted gene co-expression network analysis (WGCNA) showed that two co-expressed gene modules, MEdarkturquoise and MEpurple—containing 66 and 208 genes, respectively—were associated with S. scitamineum infection. The genes in these two modules were further studied using Gene Ontology (GO) enrichment analysis, pathogen-host interaction (PHI) database BLASTp, and small secreted cysteine-rich proteins (SCRPs) prediction. The top ten hub genes in each module were identified using the Cytohubba plugin. The GO enrichment analysis found that endoplasmic reticulum-related and catabolism-related genes were expressed during S. scitamineum infection. A total of 83 genes had homologs in the PHI database, 62 of which correlated with pathogen virulence. A total of 21 proteins had the characteristics of small secreted cysteine-rich proteins (SCRPs), a common source of fungal effectors. The top ten hub genes in each module were identified, and seven were annotated as Mig1-Mig1 protein, glycosyl hydrolase, beta-N-acetylglucosaminidase, secreted chorismate mutase, collagen, mRNA export factor, and pleckstrin homology domain protein, while the remaining three were unknown. Two SCRPs—SPSC_06609 and SPSC_04676—and three proteins—SPSC_01958, SPSC_02155, and SPSC_00940—identified in the PHI database were also among the top ten hub genes in the MEdarkturquoise and MEpurple modules, suggesting that they may play important roles in S. scitamineum infection. A S. scitamineum infection model was postulated based on current findings. These findings help to deepen the current understanding of early events in S. scitamineum infection.

Background: The pine wood nematode (PWN), Bursaphelenchus xylophilus, is a devastating pathogen of many Pinus species in China. The aim of this study was to understand the interactive molecular mechanism of PWN and its host by comparing differentially expressed genes and candidate effectors from three transcriptomes of B. xylophilus at different infection stages. Results: In total, 62, 69 and 46 candidate effectors were identified in three transcriptomes (2.5 h postinfection, 6, 12 and 24 h postinoculation and 6 and 15 d postinfection, respectively). In addition to uncharacterized pioneers, other candidate effectors were involved in the degradation of host tissues, suppression of host defenses, targeting plant signaling pathways, feeding and detoxification, which helped B. xylophilus survive successfully in the host. Seven candidate effectors were identified in both our study and the B. xylophilus transcriptome at 2.5 h postinfection, and one candidate effector was identified in all three transcriptomes. These common candidate effectors were upregulated at infection stages, and one of them suppressed pathogen-associated molecular pattern (PAMP) PsXEG1-triggered cell death in Nicotiana benthamiana. Conclusions: The results indicated that B. xylophilus secreted various candidate effectors, and some of them continued to function throughout all infection stages. These various candidate effectors were important to B. xylophilus infection and survival, and they functioned in different ways (such as breaking down host cell walls, suppressing host defenses, promoting feeding efficiency, promoting detoxification and playing virulence functions). The present results provide valuable resources for in-depth research on the pathogenesis of B. xylophilus from the perspective of effectors. [ABSTRACT FROM AUTHOR]

Breeding programs of five-needle pines have documented both major gene resistance (MGR) and quantitative disease resistance (QDR) to Cronartium ribicola (Cri), a non-native, invasive fungal pathogen causing white pine blister rust (WPBR). WPBR is one of the most deadly forest diseases in North America. However, Cri virulent pathotypes have evolved and can successfully infect and kill trees carrying resistance (R) genes, including vcr2 that overcomes MGR conferred by the western white pine (WWP, Pinus monticola) R gene (Cr2). In the absence of a reference genome, the present study generated a vcr2 reference transcriptome, consisting of about 20,000 transcripts with 1,014 being predicted to encode secreted proteins (SPs). Comparative profiling of transcriptomes and secretomes revealed vcr2 was significantly enriched for several gene ontology (GO) terms relating to oxidation-reduction processes and detoxification, suggesting that multiple molecular mechanisms contribute to pathogenicity of the vcr2 pathotype for its overcoming Cr2. RNA-seq-based bulked segregant analysis (BSR-Seq) revealed genome-wide DNA variations, including about 65,617 single nucleotide polymorphism (SNP) loci in 7,749 polymorphic genes shared by vcr2 and avirulent (Avcr2) pathotypes. An examination of the distribution of minor allele frequency (MAF) uncovered a high level of genomic divergence between vcr2 and Avcr2 pathotypes. By integration of extreme-phenotypic genome-wide association (XP-GWAS) analysis and allele frequency directional difference (AFDD) mapping, we identified a set of vcr2-associated SNPs within functional genes, involved in fungal virulence and other molecular functions. These included six SPs that were top candidate effectors with putative activities of reticuline oxidase, proteins with common in several fungal extracellular membrane (CFEM) domain or ferritin-like domain, polysaccharide lyase, rds1p-like stress responsive protein, and two Cri-specific proteins without annotation. Candidate effectors and vcr2-associated genes provide valuable resources for further deciphering molecular mechanisms of virulence and pathogenicity by functional analysis and the subsequent development of diagnostic tools for monitoring the virulence landscape in the WPBR pathosystems.

Breeding programs of five-needle pines have documented both major gene resistance (MGR) and quantitative disease resistance (QDR) to Cronartium ribicola (Cri), a non-native, invasive fungal pathogen causing white pine blister rust (WPBR). WPBR is one of the most deadly forest diseases in North America. However, Cri virulent pathotypes have evolved and can successfully infect and kill trees carrying resistance (R) genes, including vcr2 that overcomes MGR conferred by the western white pine (WWP, Pinus monticola) R gene (Cr2). In the absence of a reference genome, the present study generated a vcr2 reference transcriptome, consisting of about 20,000 transcripts with 1,014 being predicted to encode secreted proteins (SPs). Comparative profiling of transcriptomes and secretomes revealed vcr2 was significantly enriched for several gene ontology (GO) terms relating to oxidation-reduction processes and detoxification, suggesting that multiple molecular mechanisms contribute to pathogenicity of the vcr2 pathotype for its overcoming Cr2. RNA-seq-based bulked segregant analysis (BSR-Seq) revealed genome-wide DNA variations, including about 65,617 single nucleotide polymorphism (SNP) loci in 7,749 polymorphic genes shared by vcr2 and avirulent (Avcr2) pathotypes. An examination of the distribution of minor allele frequency (MAF) uncovered a high level of genomic divergence between vcr2 and Avcr2 pathotypes. By integration of extreme-phenotypic genome-wide association (XP-GWAS) analysis and allele frequency directional difference (AFDD) mapping, we identified a set of vcr2 -associated SNPs within functional genes, involved in fungal virulence and other molecular functions. These included six SPs that were top candidate effectors with putative activities of reticuline oxidase, proteins with common in several fungal extracellular membrane (CFEM) domain or ferritin-like domain, polysaccharide lyase, rds1p-like stress responsive protein, and two Cri-specific proteins without annotation. Candidate effectors and vcr2 -associated genes provide valuable resources for further deciphering molecular mechanisms of virulence and pathogenicity by functional analysis and the subsequent development of diagnostic tools for monitoring the virulence landscape in the WPBR pathosystems. [ABSTRACT FROM AUTHOR]

Background: Rust fungi constitute the largest group of plant fungal pathogens. However, a paucity of data, including genomic sequences, transcriptome sequences, and associated molecular markers, hinders the development of inhibitory compounds and prevents their analysis from an evolutionary perspective. Gymnosporangium yamadae and G. asiaticum are two closely related rust fungal species, which are ecologically and economically important pathogens that cause apple rust and pear rust, respectively, proved to be devastating to orchards. In this study, we investigated the transcriptomes of these two Gymnosporangium species during the telial stage of their lifecycles. The aim of this study was to understand the evolutionary patterns of these two related fungi and to identify genes that developed by selection. Results: The transcriptomes of G. yamadae and G. asiaticum were generated from a mixture of RNA from three biological replicates of each species. We obtained 49,318 and 54,742 transcripts, with N50 values of 1957 and 1664, for G. yamadae and G. asiaticum, respectively. We also identified a repertoire of candidate effectors and other gene families associated with pathogenicity. A total of 4947 pairs of putative orthologues between the two species were identified. Estimation of the non-synonymous/synonymous substitution rate ratios for these orthologues identified 116 pairs with Ka/Ks values greater than1 that are under positive selection and 170 pairs with Ka/Ks values of 1 that are under neutral selection, whereas the remaining 4661 genes are subjected to purifying selection. We estimate that the divergence time between the two species is approximately 5.2 Mya. Conclusion: This study constitutes a de novo assembly and comparative analysis between the transcriptomes of the two rust species G. yamadae and G. asiaticum. The results identified several orthologous genes, and many expressed genes were identified by annotation. Our analysis of Ka/Ks ratios identified orthologous genes subjected to positive or purifying selection. An evolutionary analysis of these two species provided a relatively precise divergence time. Overall, the information obtained in this study increases the genetic resources available for research on the genetic diversity of the Gymnosporangium genus. [ABSTRACT FROM AUTHOR]

Uromyces fabae is a major pathogen of broad bean, Vicia faba. U. fabae has served as a model among rust fungi to elucidate the development of infection structures, expression and secretion of cell wall degrading enzymes and gene expression. Using U. fabae, enormous progress was made regarding nutrient uptake and metabolism and in the search for secreted proteins and effectors. Here, we present results from a genome survey of U. fabae. Paired end Illumina sequencing provided 53 Gb of data. An assembly gave 59,735 scaffolds with a total length of 216 Mb. K-mer analysis estimated the genome size to be 329 Mb. Of a representative set of 23,153 predicted proteins we could annotate 10,209, and predict 599 secreted proteins. Clustering of the protein set indicates families of highly likely effectors. We also found new homologs of RTP1p, a prototype rust effector. The U. fabae genome will be an important resource for comparative analyses with U. appendiculatus and P. pachyrhizi and provide information regarding the phylogenetic relationship of the genus Uromyces with respect to other rust fungi already sequenced, namely Puccinia graminis f. sp. tritici, P. striiformis f. sp. tritici, Melampsora lini, and Melampsora larici-populina.

Uromyces fabae is a major pathogen of broad bean, Vicia faba. U. fabae has served as a model among rust fungi to elucidate the development of infection structures, expression and secretion of cell wall degrading enzymes and gene expression. Using U. fabae, enormous progress was made regarding nutrient uptake and metabolism and in the search for secreted proteins and effectors. Here, we present results from a genome survey of U. fabae. Paired end Illumina sequencing provided 53 Gb of data. An assembly gave 59,735 scaffolds with a total length of 216 Mb. K-mer analysis estimated the genome size to be 329 Mb. Of a representative set of 23,153 predicted proteins we could annotate 10,209, and predict 599 secreted proteins. Clustering of the protein set indicates families of highly likely effectors. We also found new homologs of RTP1p, a prototype rust effector. The U. fabae genome will be an important resource for comparative analyses with U. appendiculatus and P. pachyrhizi and provide information regarding the phylogenetic relationship of the genus Uromyces with respect to other rust fungi already sequenced, namely Puccinia graminis f. sp. tritici, P. striiformis f. sp. tritici, Melampsora lini, and Melampsora larici-populina. [ABSTRACT FROM AUTHOR]

Background: Plant and animal pathogenic fungus Cochliobolus lunatus cause great economic damages worldwide every year. C. lunatus displays an increased temperature dependent-virulence to a wide range of hosts. Nonetheless, this phenomenon is poorly understood due to lack of insights on the coordinated secretome weaponries produced by C. lunatus under heat-stress conditions on putative hosts. To understand the mechanism better, we dissected the secretome of C. lunatus interacting with potato (Solanum tuberosum L.) leaf at different temperature regimes. Results: C. lunatus produced melanized colonizing hyphae in and on potato leaf, finely modulated the ambient pH as a function of temperature and secreted diverse set of proteins. Using two dimensional gel electrophoresis (2-D) and mass spectrometry (MS) technology, we observed discrete secretomes at 20°C, 28°C and 38°C. A total of 21 differentially expressed peptide spots and 10 unique peptide spots (that did not align on the gels) matched with 28 unique protein models predicted from C. lunatus m118 v.2 genome peptides. Furthermore, C. lunatus secreted peptides via classical and non-classical pathways related to virulence, proteolysis, nucleic acid metabolism, carbohydrate metabolism, heat stress, signal trafficking and some with unidentified catalytic domains. Conclusions: We have identified a set of 5 soluble candidate effectors of unknown function from C. lunatus secretome weaponries against potato crop at different temperature regimes. Our findings demonstrate that C. lunatus has a repertoire of signature secretome which mediates thermo-pathogenicity and share a leucine rich “CL[xxxx]LHM”-motif. Considering the rapidly evolving temperature dependent-virulence and host diversity of C. lunatus, this data will be useful for designing new protection strategies. [ABSTRACT FROM AUTHOR]

Fusarium blight of wheat is usually caused by Fusarium graminearum , and the pathogenic fungi will secrete effectors into the host plant tissue to affect its normal physiological process, so as to make it pathogenic. The CFEM (Common in Fungal Extracellular Membrane) protein domain is unique to fungi, but it is not found in all fungi. The CFEM protein contained in F. graminearum may be closely related to pathogenicity. In this study, 23 FgCFEM proteins were identified from the F. graminearum genome. Then, features of these proteins, such as signal peptide, subcellular localization, and transmembrane domains, etc., were analyzed and candidate effectors were screened out. Sequence alignment results revealed that each FgCFEM protein contains one CFEM domain. The amino acids of the CFEM domain are highly conserved and contain eight spaced cysteines, with the exception that FgCFEM8, 9, and 15 lack two cysteines and three cysteines were missed in FgCFEM18 and FgCFEM22. A recently identified CFEM_DR motif was detected in 11 FgCFEMs, and importantly we identified two new conserved motifs containing about 29 and 18 amino acids (CFEM_WR and CFEM_KF), respectively, in some of FgCFEM proteins. Transcriptome analysis of the genes encoding CFEM proteins indicated that all the CFEM-containing genes were expressed during wheat infection, with seven and six genes significantly up- and down-regulated, respectively, compared with in planta and in vitro. Based on the above analysis, FgCFEM11 and FgCFEM23 were predicted to be F. graminearum effectors. This study provides the basis for future functional analyses of CFEM proteins in F. graminearum .

Background: Rust fungi constitute the largest group of plant fungal pathogens. However, a paucity of data, including genomic sequences, transcriptome sequences, and associated molecular markers, hinders the development of inhibitory compounds and prevents their analysis from an evolutionary perspective. Gymnosporangium yamadae and G. asiaticum are two closely related rust fungal species, which are ecologically and economically important pathogens that cause apple rust and pear rust, respectively, proved to be devastating to orchards. In this study, we investigated the transcriptomes of these two Gymnosporangium species during the telial stage of their lifecycles. The aim of this study was to understand the evolutionary patterns of these two related fungi and to identify genes that developed by selection., Results: The transcriptomes of G. yamadae and G. asiaticum were generated from a mixture of RNA from three biological replicates of each species. We obtained 49,318 and 54,742 transcripts, with N50 values of 1957 and 1664, for G. yamadae and G. asiaticum, respectively. We also identified a repertoire of candidate effectors and other gene families associated with pathogenicity. A total of 4947 pairs of putative orthologues between the two species were identified. Estimation of the non-synonymous/synonymous substitution rate ratios for these orthologues identified 116 pairs with Ka/Ks values greater than1 that are under positive selection and 170 pairs with Ka/Ks values of 1 that are under neutral selection, whereas the remaining 4661 genes are subjected to purifying selection. We estimate that the divergence time between the two species is approximately 5.2 Mya., Conclusion: This study constitutes a de novo assembly and comparative analysis between the transcriptomes of the two rust species G. yamadae and G. asiaticum. The results identified several orthologous genes, and many expressed genes were identified by annotation. Our analysis of Ka/Ks ratios identified orthologous genes subjected to positive or purifying selection. An evolutionary analysis of these two species provided a relatively precise divergence time. Overall, the information obtained in this study increases the genetic resources available for research on the genetic diversity of the Gymnosporangium genus.