Rice production in the Anhui province is threatened by fungal diseases. We obtained twenty-five fungal isolates from rice and wild rice leaves showing leaf spot disease collected along the Yangtze River. A phylogenetic analysis based on internal transcribed spacer (ITS), translation elongation factor 1 alpha (TEF1-α), and beta tubulin (TUB2) sequences revealed one isolate (SS-2-JB-1B) grouped with Nigrospora sphaerica, one (QY) with Nigrospora chinensis, twenty-two with Nigrospora oryzae, and one isolate (QY-2) grouped in its own clade, which are related to but clearly different from N. oryzae. Nineteen tested isolates, including sixteen strains from the N. oryzae clade and the three isolates of the other three clades, caused disease on detached rice leaves. The three isolates that did not belong to N. oryzae were also able to cause disease in rice seedlings, suggesting that they were rice pathogens. Isolate QY-2 differed from the other isolates in terms of colony morphology, cell size, and susceptibility to fungicides, indicating that this isolate represents a new species that we named Nigrospora anhuiensis. Our analysis showed that N. sphaerica, N. chinensis, and the new species, N. anhuiensis, can cause rice leaf spot disease in the field. This research provides new knowledge for understanding rice leaf spot disease.

Terpenoids are important compounds associated with the pest and herbivore resistance mechanisms of plants; consequently, it is essential to identify and explore terpene synthase (TPS) genes in maize. In the present study, we identified 31 TPS genes based on a pan-genome of 26 high-quality maize genomes containing 20 core genes (present in all 26 lines), seven dispensable genes (present in 2 to 23 lines), three near-core genes (present in 24 to 25 lines), and one private gene (present in only 1 line). Evaluation of ka/ks values of TPS in 26 varieties revealed that TPS25 was subjected to positive selection in some varieties. Six ZmTPS had ka/ks values less than 1, indicating that they were subjected to purifying selection. In 26 genomes, significant differences were observed in ZmTPS25 expression between genes affected by structural variation (SV) and those not affected by SV. In some varieties, SV altered the conserved structural domains resulting in a considerable number of atypical genes. The analysis of RNA-seq data of maize Ostrinia furnacalis feeding revealed 10 differentially expressed ZmTPS, 9 of which were core genes. However, many atypical genes for these responsive genes were identified in several genomes. These findings provide a novel resource for functional studies of ZmTPS.

Terpene synthase (TPS) catalyzes the synthesis of terpenes and plays an important role in plant defense. This study identified 45 OsTPS genes (32 core genes and 13 variable genes) based on the high-quality rice gene-based pan-genome. This indicates limitations in OsTPS gene studies based on a single reference genome. In the present study, through collinearity between multiple rice genomes, one OsTPS gene absent in the reference (Nipponbare) genome was found and two TPS genes in the reference genome were found to have atypical structures, which would have been ignored in single genome analysis. OsTPS genes were divided into five groups and TPS-b was lost according to the phylogenetic tree. OsTPSs in TPS-c and TPS-g were all core genes indicating these two groups were stable during domestication. In addition, through the analysis of transcriptome data, some structural variations were found to affect the expression of OsTPS genes. Through the Ka/Ks calculation of OsTPS genes, we found that different OsTPS genes were under different selection pressure during domestication; for example, OsTPS22 and OsTPS29 experienced stronger positive selection than the other OsTPS genes. After Chilo suppressalis larvae infesting, 25 differentially expressed OsTPS genes were identified, which are involved in the diterpene phytoalexins precursors biosynthesis and ent-kaurene biosynthesis pathways. Overall, the present study conducted a bioinformatics analysis of OsTPS genes using a high-quality rice pan-genome, which provided a basis for further study of OsTPS genes.

Melon (Cucumismelo L.) is an important vegetable crop that has been subjected to domestication and improvement. Several varieties of melons with diverse phenotypes have been produced. In this study, we constructed a melon pan-genome based on 297 accessions comprising 168 Mb novel sequences and 4,325 novel genes. Based on the results, there were abundant genetic variations among different melon groups, including 364 unfavorable genes in the IMP_A vs. LDR_A group, 46 favorable genes, and 295 unfavorable genes in the IMP_M vs. LDR_M group. The distribution of 709 resistance gene analogs (RGAs) was also characterized across 297 melon lines, of which 603 were core genes. Further, 106 genes were found to be variable, 55 of which were absent in the reference melon genome. Using gene presence/absence variation (PAV)-based genome-wide association analysis (GWAS), 13 gene PAVs associated with fruit length, fruit shape, and fruit width were identified, four of which were located in pan-genome additional contigs.

Soft rot caused by Dickeya zeae is an important bacterial disease affecting rice and other plants worldwide. In this study, Nanopore and Illumina sequencing platforms were used to sequence the high-quality complete genome of a novel D. zeae strain WH1 (size: 4.68 Mb; depth: 322.37x for Nanopore, 243.51x for Illumina; GC content: 53.59%), which was isolated from healthy rice root surface together with Paenibacillus polymyxa, a potential biocontrol bacterium against D. zeae strain WH1. However, the pure WH1 culture presented severe pathogenicity. Multilocus sequence analysis (MLSA) indicated that strains WH1, EC1, and EC2 isolated from rice were grouped into a clade differentiated from other D. zeae strains. The average nucleotide identity (ANI) and DNA-DNA hybridization (DDH) analyses demonstrated that WH1 was phylogenetically closest to EC2. Furthermore, the pathogenicity determinants and virulence factors of WH1 were mainly analyzed through genomic comparison with complete genomes of other D. zeae strains with high virulence (EC1, EC2, MS1, and MS2). The results revealed that plant cell wall-degrading extracellular enzymes (PCWDEs), flagellar and chemotaxis, and quorum sensing were highly conserved in all analyzed genomes, which were confirmed through phenotypic assays. Besides, WH1 harbored type I, II, III, and VI secretion systems (T1SS, T2SS, T3SS, and T6SS), but lost T4SS and T5SS. Like strains MS1 and MS2 isolated from bananas, WH1 harbored genes encoding both capsule polysaccharide (CPS) and exopolysaccharide (EPS) biosynthesis. The results of pathogenicity assays demonstrated that WH1 produced severe soft rot symptoms on potato tubers, carrots, radishes, and Chinese cabbage. Meanwhile, WH1 also produced phytotoxin(s) to inhibit rice seed germination with an 87% inhibitory rate in laboratory conditions. More importantly, we confirmed that phytotoxin(s) produced by WH1 are different from zeamines produced by EC1. Comparative genomics analyses and phenotypic and pathogenicity assays suggested that WH1 likely evolved through a pathway different from the other D. zeae strains from rice, producing a new type of rice foot rot pathogen. These findings highlight the emergence of a new type of D. zeae strain with high virulence, causing soft rot in rice and other plants.

Spodoptera frugiperda is a common and severely damaging agricultural pest. In-depth analysis of its population genomics and transcriptomics is crucial for providing references for pest control efforts. This study, focused on the extensive variation in the genome size of S. frugiperda, constructed its pan-genome and identified 1.37 Gb of non-reference sequences, highlighting significant genetic variation within the population. Analysis of Long Terminal Repeat (LTR) Presence/Absence Variation (PAV) suggests that LTR alterations may be one of the driving factors for genome size variation. Additionally, population gene PAV analysis revealed that variable genes are enriched in functions like acetyltransferase activity, which might be associated with detoxification, implying diverse selection pressures related to detoxification in different S. frugiperda populations. Moreover, 19 horizontal gene transfer (HGT) acquired genes were identified in the reference genome used in this study, which responded to 16 different treatments. Notably, three HGT-acquired genes (SFR02618, SFR05248, and SFR05249) co-expressed with heat shock protein family and responded under treatments with Avermectin and Cypermethrin. This may indicate their involvement in a detoxification mechanism coordinated with heat shock proteins. These results offering new insights into its genomic evolution and the potential functions of HGT-acquired genes.