Comparative genomics and association analysis identifies virulence genes of Cercospora sojina in soybean

Genome Sequencing and Comparative Genomic Analysis of Attenuated Strain Gibellulopsis nigrescens GnVn.1 Causing Mild Wilt in Sunflower

Gibellulopsis nigrescens, previously classified in the Verticillium genus until 2007, is an attenuated pathogen known to provide cross-protection against Verticillium wilt in various crops. To investigate the potential mechanisms underlying its reduced virulence, we conducted genome sequencing, annotation, and a comparative genome analysis of G. nigrescens GnVn.1 (GnVn.1), an attenuated strain isolated from sunflower. The genome sequencing and annotation results revealed that the GnVn.1 genome consists of 22 contigs, with a total size of 31.79 Mb. We predicted 10,876 genes, resulting in a gene density of 342 genes per Mb. The pathogenicity gene prediction results indicated 1733 high-confidence pathogenicity factors (HCPFs), 895 carbohydrate-active enzymes (CAZys), and 359 effectors. Moreover, we predicted 40 secondary metabolite clusters (SMCs). The comparative genome analysis indicated that GnVn.1 contains more CAZys, SMCs, predicted effectors, and HCPF genes than Verticillium dahliae (VdLs.17) and Verticillium alfalfae (VaMas.102). The core-pan analysis results showed that GnVn.1 had more specific HCPFs, effectors, CAZys, and secreted protein (SP) genes, and lost many critical pathogenic genes compared to VdLs.17 and VaMs.102. Our results indicate that the GnVn.1 genome harbors more pathogenicity-related genes than the VdLs.17 and VaMs.102 genomes. These abundant genes may play critical roles in regulating virulence. The loss of critical pathogenic genes causes weak virulence and confers biocontrol strategies to GnVn.1.

Population Structure and Mating Type Distribution of Cercospora sojina from Soybeans in Indiana, United States

Frogeye leaf spot on soybeans is traditionally considered as a southern disease in the United States but its impact in North Central USA has been rising in recent years. In this study, we investigated the population structure and mating type distribution in the C. sojina population from Indiana, USA. Based on 27 single nucleotide polymorphism markers, 49 multi-locus genotypes (MLGs) were identified in 234 isolates collected from 29 counties in Indiana in 2020. Bayesian analysis grouped the 49 MLGs into three clusters. This grouping was supported by principal coordinate analysis and, in large part, by the unweighted pair group method with arithmetic mean and minimal spanning tree. Only one mating-type idiomorph was found in each isolate and in each MLG. The MAT1-1 idiomorph was found in 22 MLGs and the MAT1-2 idiomorph was found in 27 MLGs. Based on clone-corrected data, the distribution of mating-type idiomorphs did not deviate significantly from 1:1 ratio in Indiana as a whole and in 22 out of 24 counties where two or more MLGs were found. Thirty MLGs contained QoI-resistant isolates and 22 MLGs contained QoI-sensitive isolates, with three MLGs containing both types of isolates. MLG1, the most common MLG with 90 isolates, contained mostly QoI-resistant isolates. Interestingly, MLG1 was also the dominant genotype in the Tennessee population collected in 2015, suggesting that MLG1 has been a dominant genotype in a wider region for many years. Based on the standard index of association (r¯d), the Indiana population as a whole was in significant linkage disequilibrium. However, in five out of 16 counties where three or more MLGs were found, the null hypothesis of linkage equilibrium was not rejected. Tests of linkage disequilibrium between locus pairs showed that 33.3% of locus pairs on the same contigs were in significant disequilibrium and 17.7% of locus pairs on different contigs were in significant disequilibrium. The possibility of a cryptic sexual stage was discussed.

Pan-genome survey of Septoria pistaciarum, causal agent of Septoria leaf spot of pistachios, across three Aegean sub-regions of Greece

Septoria pistaciarum, a causal agent of Septoria leaf spot disease of pistachio, is a fungal pathogen that causes substantial losses in the cultivation, worldwide. This study describes the first pan-genome-based survey of this phytopathogen-comprising a total of 27 isolates, with 9 isolates each from 3 regional units of Greece (Pieria, Larissa and Fthiotida). The reference isolate (SPF8) assembled into a total of 43.1 Mb, with 38.6% contained within AT-rich regions of approximately 37.5% G:C. The genomes of the 27 isolates exhibited on average 42% gene-coding and 20% repetitive regions. The genomes of isolates from the southern Fthiotida region appeared to more diverged from each other than the other regions based on SNP-derived trees, and also contained isolates similar to both the Pieria and Larissa regions. In contrast, isolates of the Pieria and Larissa were less diverse and distinct from one another. Asexual reproduction appeared to be typical, with no MAT1-2 locus detected in any isolate. Genome-based prediction of infection mode indicated hemibiotrophic and saprotrophic adaptations, consistent with its long latent phase. Gene prediction and orthology clustering generated a pan-genome-wide gene set of 21,174 loci. A total of 59 ortholog groups were predicted to contain candidate effector proteins, with 36 (61%) of these either having homologs to known effectors from other species or could be assigned predicted functions from matches to conserved domains. Overall, effector prediction suggests that S. pistaciarum employs a combination of defensive effectors with roles in suppression of host defenses, and offensive effectors with a range of cytotoxic activities. Some effector-like ortholog groups presented as divergent versions of the same protein, suggesting region-specific adaptations may have occurred. These findings provide insights and future research directions in uncovering the pathogenesis and population dynamics of S. pistaciarum toward the efficient management of Septoria leaf spot of pistachio.

Functional Evolution of Pseudofabraea citricarpa as an Adaptation to Temperature Change

Citrus target spot, caused by Pseudofabraea citricarpa, was formerly considered a cold-tolerant fungal disease. However, it has now spread from high-latitude regions to warmer low-latitude regions. Here, we conducted physiological observations on two different strains of the fungus collected from distinct regions, and evaluated their pathogenicity. Interestingly, the CQWZ collected from a low-latitude orchard, exhibited higher temperature tolerance and pathogenicity when compared to the SXCG collected from a high-latitude orchard. To further understand the evolution of temperature tolerance and virulence in these pathogens during the spread process, as well as the mechanisms underlying these differences, we performed genomic comparative analysis. The genome size of CQWZ was determined to be 44,004,669 bp, while the genome size of SXCG was determined to be 45,377,339 bp. Through genomic collinearity analysis, we identified two breakpoints and rearrangements during the evolutionary process of these two strains. Moreover, gene annotation results revealed that the CQWZ possessed 376 annotated genes in the "Xenobiotics biodegradation and metabolism" pathway, which is 79 genes more than the SXCG. The main factor contributing to this difference was the presence of salicylate hydroxylase. We also observed variations in the oxidative stress pathways and core pathogenic genes. The CQWZ exhibited the presence of a heat shock protein (HSP SSB), a catalase (CAT2), and 13 core pathogenic genes, including a LysM effector, in comparison to the SXCG. Furthermore, there were significant disparities in the gene clusters responsible for the production of seven metabolites, such as Fumonisin and Brefeldin. Finally, we identified the regulatory relationship, with the HOG pathway at its core, that potentially contributes to the differences in thermotolerance and virulence. As the global climate continues to warm, crop pathogens are increasingly expanding to new territories. Our findings will enhance understanding of the evolution mechanisms of pathogens under climate change.

Efficacy and Profitability of Fungicides for Managing Frogeye Leaf Spot on Soybean in the United States: A 10-Year Quantitative Summary

Frogeye leaf spot (FLS), caused by Cercospora sojina, is an economically important disease of soybean in the United States. Data from 66 uniform fungicide trials (UFTs) conducted from 2012 to 2021 across eight states (Alabama, Arkansas, Illinois, Iowa, Kentucky, Louisiana, Mississippi, and Tennessee) were gathered and analyzed to determine the efficacy and profitability of the following fungicides applied at the beginning pod developmental stage (R3): azoxystrobin + difenoconazole (AZOX + DIFE), difenoconazole + pydiflumetofen (DIFE + PYDI), pyraclostrobin (PYRA), pyraclostrobin + fluxapyroxad + propiconazole (PYRA + FLUX + PROP), tetraconazole (TTRA), thiophanate-methyl (TMET), thiophanate-methyl + tebuconazole (TMET + TEBU), and trifloxystrobin + prothioconazole (TFLX + PROT). A network meta-analytic model was fitted to the log of the means of FLS severity data and to the nontransformed mean yield for each treatment, including the nontreated. The percent reduction in disease severity (%) and the yield response (kg/ha) relative to the nontreated was the lowest for PYRA (11%; 136 kg/ha) and the greatest for DIFE + PYDI (57%; 441 kg/ha). A significant decline in efficacy over time was detected for PYRA (18 percentage points [p.p.]), TTRA (27 p.p.), AZOX + DIFE (18 p.p.), and TMET + TEBU (19 p.p.) by using year as a continuous covariate in the model. Finally, probabilities of breaking even were the greatest (>65%) for the most effective fungicide DIFE + PYDI and the lowest (<55%) for PYRA. Results of this meta-analysis may be useful to support decisions when planning fungicide programs.

Comparative genomics predict specific genes in potential mucorales identification

Mucoralean fungi could cause mucormycosis in humans, particularly in immunodeficient individuals and those with diabetes mellitus or trauma. With plenty of species and genera, their molecular identification and pathogenicity have a large deviation. Reported cases of mucormycosis showed frequent occurrence in Rhizopus species, Mucor species, and Lichtheimia species. We analyzed the whole genome sequences of 25 species of the top 10 Mucorales genera, along with another 22 important pathogenic non-Mucorales species, to dig the target genes for monitoring Mucorales species and identify potential genomic imprints of virulence in them. Mucorales-specific genes have been found in various orthogroups extracted by Python script, while genus-specific genes were annotated covering cellular structure, biochemistry metabolism, molecular processing, and signal transduction. Proteins related to the virulence of Mucorales species varied with distinct significance in copy numbers, in which Orthofinder was conducted. Based on our fresh retrospective analysis of mucormycosis, a comparative genomic analysis of pathogenic Mucorales was conducted in more frequent pathogens. Specific orthologs between Mucorales and non-Mucoralean pathogenic fungi were discussed in detail. Referring to the previously reported virulence proteins, we included more frequent pathogenic Mucorales and compared them in Mucorales species and non-Mucorales species. Besides, more samples are needed to further verify the potential target genes.

Whole Genome Sequencing and Comparative Genomics of Indian Isolates of Wheat Spot Blotch Pathogen Bipolaris sorokiniana Reveals Expansion of Pathogenicity Gene Clusters

Spot blotch is a highly destructive disease in wheat caused by the fungal pathogen Bipolaris sorokiniana (teleomorph, Cochliobolus sativus). It is prevalent in warm and humid areas, including Africa, Asia, Latin America, and the USA. In the present study, twelve isolates of B. sorokiniana were collected from wheat fields in three different geographical locations in India. The pathogenicity of seven sporulating isolates was assessed on 'DDK 1025', a spot blotch-susceptible wheat variety under greenhouse conditions. The isolate 'D2' illustrated the highest virulence, followed by 'SI' and 'BS52'. These three isolates were sequenced using the Illumina HiSeq1000 platform. The estimated genome sizes of the isolates BS52, D2, and SI were 35.19 MB, 39.32 MB, and 32.76 MB, with GC contents of 48.48%, 50.43%, and 49.42%, respectively. The numbers of pathogenicity genes identified in BS52, D2, and SI isolates were 2015, 2476, and 2018, respectively. Notably, the isolate D2 exhibited a relatively larger genome with expanded arsenals of Biosynthetic Gene Clusters (BGCs), CAZymes, secretome, and pathogenicity genes, which could have contributed to its higher virulence among the tested isolates. This study provides the first comparative genome analysis of the Indian isolates of B. sorokiniana using whole genome sequencing.

Genome-Wide Association Analysis and Gene Mining of Resistance to China Race 1 of Frogeye Leaf Spot in Soybean

Soybean frogeye leaf spot (FLS) is a worldwide fungal disease. Its higher occurrence frequency and wider distribution range always led to severe yield losses of soybean, therefore, breeding new cultivars with FLS resistance has been an important breeding goal for soybean breeders. In this study, an association panel of 183 representative soybean accessions was used to evaluate their resistance to FLS race 1, and to identify quantitative trait nucleotides (QTNs) and candidate genes based on genome-wide association study (GWAS) and high-throughput single-nucleotide polymorphisms (SNPs). A total of 23,156 high-quality SNPs were developed using the specific locus-amplified fragment sequencing (SLAF-seq) approach. Finally, 13 novel association signals associated with FLS race 1 resistance were identified by the compressed mixed linear model (CMLM). In addition, 119 candidate genes were found within the 200-kb flanking genomic region of these 13 peak SNPs. Based on the gene-based association analysis, haplotype analysis, expression pattern analysis, and virus-induced gene silencing (VIGS) systems, four genes (Glyma.05G121100, Glyma.17G228300, Glyma.19G006900, and Glyma.19G008700) were preliminarily proved to play an important role in the soybean resistance to FLS race 1.

Genome-wide association of single nucleotide polymorphism loci and candidate genes for frogeye leaf spot (Cercospora sojina) resistance in soybean

BACKGROUND

Frogeye leaf spot (FLS) is a destructive fungal disease that affects soybean production. The most economical and effective strategy to control FLS is the use of resistant cultivars. However, the use of a limited number of resistant loci in FLS management will be countered by the emergence of new high-virulence Cercospora sojina races. Therefore, we identified quantitative trait loci (QTL) that control resistance to FLS and identified novel resistant genes using a genome-wide association study (GWAS) on 234 Chinese soybean cultivars.

RESULTS

A total of 30,890 single nucleotide polymorphism (SNP) markers were used to estimate linkage disequilibrium (LD) and population structure. The GWAS results showed four loci (p < 0.0001) distributed over chromosomes (Chr.) 5 and 20, that are significantly associated with FLS resistance. No previous studies have reported resistance loci in these regions. Subsequently, 45 genes in the two resistance-related haplotype blocks were annotated. Among them, Glyma20g31630 encoding pyruvate dehydrogenase (PDH), Glyma05g28980, which encodes mitogen-activated protein kinase 7 (MPK7), and Glyma20g31510, Glyma20g31520 encoding calcium-dependent protein kinase 4 (CDPK4) in the haplotype blocks deserves special attention.

CONCLUSIONS

This study showed that GWAS can be employed as an effective strategy for identifying disease resistance traits in soybean and narrowing SNPs and candidate genes. The prediction of candidate genes in the haplotype blocks identified by disease resistance loci can provide a useful reference to study systemic disease resistance.

High-quality genome assembly of the soybean fungal pathogen Cercospora kikuchii

Plant diseases caused by the Cercospora genus of ascomycete fungi are a major concern for commercial agricultural practices. Several Cercospora species can affect soybeans, such as Cercospora kikuchii which causes soybean leaf blight. Speciation in Cercospora on soybean has not been adequately studied. Some cryptic groups of Cercospora also cause diseases on soybean. Moreover, it has been known C. kikuchii population genetic structure is different between countries. Consequently, further genomic information could help to elucidate the covert differentiation of Cercospora diseases in soybean. Here, we report for the first time, a chromosome-level genome assembly for C. kikuchii. The genome assembly of 9 contigs was 34.44 Mb and the N50 was 4.19 Mb. Based on ab initio gene prediction, several candidates for pathogenicity-related genes, including 242 genes for putative effectors, 55 secondary metabolite gene clusters, and 399 carbohydrate-active enzyme genes were identified. The genome sequence and the features described in this study provide a solid foundation for comparative and evolutionary genomic analysis for Cercospora species that cause soybean diseases worldwide.

Draft Genome Sequences of Three Fusarium oxysporum f. sp. niveum Isolates Used in Designing Markers for Race Differentiation

Here, we report the draft genome sequences of three Fusarium oxysporum f. sp. niveum isolates that were used to design markers for molecular race differentiation. The isolates were collected from watermelon fields in Georgia (USA) and were determined to be different races of F. oxysporum f. sp. niveum using a traditional bioassay.

Here, we report the draft genome sequences of three Fusarium oxysporum f. sp. niveum isolates that were used to design markers for molecular race differentiation. The isolates were collected from watermelon fields in Georgia (USA) and were determined to be different races of F. oxysporum f. sp. niveum using a traditional bioassay.