NRPS-like ATRR in Plant-Parasitic Nematodes Involved in Glycine Betaine Metabolism to Promote Parasitism

Plant-parasitic nematodes (PPNs) are among the most serious phytopathogens and cause widespread and serious damage in major crops. In this study, using a genome mining method, we identified nonribosomal peptide synthetase (NRPS)-like enzymes in genomes of plant-parasitic nematodes, which are conserved with two consecutive reducing domains at the N-terminus (A-T-R1-R2) and homologous to fungal NRPS-like ATRR. We experimentally investigated the roles of the NRPS-like enzyme (MiATRR) in nematode (Meloidogyne incognita) parasitism. Heterologous expression of Miatrr in Saccharomyces cerevisiae can overcome the growth inhibition caused by high concentrations of glycine betaine. RT-qPCR detection shows that Miatrr is significantly upregulated at the early parasitic life stage (J2s in plants) of M. incognita. Host-derived Miatrr RNA interference (RNAi) in Arabidopsis thaliana can significantly decrease the number of galls and egg masses of M. incognita, as well as retard development and reduce the body size of the nematode. Although exogenous glycine betaine and choline have no obvious impact on the survival of free-living M. incognita J2s (pre-parasitic J2s), they impact the performance of the nematode in planta, especially in Miatrr-RNAi plants. Following application of exogenous glycine betaine and choline in the rhizosphere soil of A. thaliana, the numbers of galls and egg masses were obviously reduced by glycine betaine but increased by choline. Based on the knowledge about the function of fungal NRPS-like ATRR and the roles of glycine betaine in host plants and nematodes, we suggest that MiATRR is involved in nematode-plant interaction by acting as a glycine betaine reductase, converting glycine betaine to choline. This may be a universal strategy in plant-parasitic nematodes utilizing NRPS-like ATRR to promote their parasitism on host plants.

Nematicidal glycosylated resorcylic acid lactones from the fungus Pochonia chlamydosporia PC-170 and their key biosynthetic genes

Chemical study of the nematicidal biocontrol fungus Pochonia chlamydosporia PC-170 led to discovery of six resorcylic acid lactones (RALs), including three nematicidal glycosylated RALs, monocillin VI glycoside (1), colletogloeolactone A (2) and monocillin II glycoside (3), and three antibacterial non-glycosylated RALs, monocillin VI (4), monocillin IV (5) and monocillin II (6). The planar structure of the new compound monocillin VI glycoside (1) was elucidated using HRESIMS and NMR data, and its monosaccharide configuration was further determined through sugar hydrolysis experiment and GC-MS analysis method. Furthermore, their two biosynthetic-related PKS genes, pchE and pchI, were identified through the gene knockout experiment. The glycosylated RALs 1-3 exhibited nematicidal activity against Meloidogyne incognita, with LC50 values of 94, 152 and 64 μg/mL, respectively, and thus had great potential in the development of new nematicidal natural products to control M. incognita in the future.

The root-knot nematode effector Mi2G02 hijacks a host plant trihelix transcription factor to promote nematode parasitism

Root-knot nematodes (RKNs) cause huge agricultural losses every year. They secrete a repertoire of effectors to facilitate parasitism through the induction of plant-derived giant feeding cells, which serve as their sole source of nutrients. However, the mode of action of these effectors and their targeted host proteins remain largely unknown. In this study, we investigated the role of the effector Mi2G02 in Meloidogyne incognita parasitism. Host-derived Mi2G02 RNA interference in Arabidopsis thaliana affected giant cell development, whereas ectopic expression of Mi2G02 promoted root growth and increased plant susceptibility to M. incognita. We used various combinations of approaches to study the specific interactions between Mi2G02 and A. thaliana GT-3a, a trihelix transcription factor. GT-3a knockout in A. thaliana affected feeding-site development, resulting in production of fewer egg masses, whereas GT-3a overexpression in A. thaliana increased susceptibility to M. incognita and also root growth. Moreover, we demonstrated that Mi2G02 plays a role in maintaining GT-3a protein stabilization by inhibiting the 26S proteasome-dependent pathway, leading to suppression of TOZ and RAD23C expression and thus promoting nematode parasitism. This work enhances our understanding of how a pathogen effector manipulates the role and regulation of a transcription factor by interfering with a proteolysis pathway to reprogram gene expression for development of nematode feeding cells.

Genome-wide transcriptome profiling reveals molecular response pathways of Trichoderma harzianum in response to salt stress

Trichoderma harzianum exhibits a strong biological control effect on many important plant pathogens, such as Fusarium oxysporum, Botrytis cinerea, and Meloidogyne. However, its biocontrol effectiveness is weakened or reduced under salt stress. The aim of this study was to investigate the molecular response of T. harzianum to salt stress at the whole-genome level. Here, we present a 44.47 Mb near-complete genome assembly of the T. harzianum qt40003 strain for the first time, which was assembled de novo with 7.59 Gb Nanopore sequencing long reads (~170-fold) and 5.2 Gb Illumina short reads (~116-fold). The assembled qt40003 genome contains 12 contigs, with a contig N50 of 4.81 Mb, in which four of the 12 contigs were entirely reconstructed in a single chromosome from telomere to telomere. The qt40003 genome contains 4.27 Mb of repeat sequences and 12,238 protein-coding genes with a BUSCO completeness of 97.5%, indicating the high accuracy and completeness of our gene annotations. Genome-wide transcriptomic analysis was used to investigate gene expression changes related to salt stress in qt40003 at 0, 2% (T2), and 4% (T4) sodium chloride concentrations. A total of 2,937 and 3,527 differentially expressed genes (DEGs) were obtained under T2 and T4 conditions, respectively. GO enrichment analysis showed that the T2-treatment DEGs were highly enriched in detoxification (p < 0.001), while the T4 DEGs were mainly enriched in cell components, mostly in cellular detoxification, cell surface, and cell wall. KEGG metabolic pathway analysis showed that 91 and 173 DEGs were significantly enriched in the T2 and T4 treatments, respectively (p < 0.01), mainly in the glutathione metabolism pathway. We further experimentally analyzed the differentially expressed glutathione transferase genes in the glutathione metabolic pathway, most of which were downregulated (13/15). In addition, we screened 13 genes related to active oxygen clearance, including six upregulated and seven downregulated genes, alongside five fungal hydrophobic proteins, of which two genes were highly expressed. Our study provides high-quality genome information for the use of T. harzianum for biological control and offers significant insights into the molecular responses of T. harzianum under salt-stress conditions.

Parasitism of Hirsutella rhossiliensis on Different Nematodes and Its Endophytism Promoting Plant Growth and Resistance against Root-Knot Nematodes

The endoparasitic fungus Hirsutella rhossiliensis is an important biocontrol agent of cyst nematodes in nature. To determine the potential parasitism of the fungus on a non-natural host, the pinewood nematode (Bursaphelenchus xylophilus) living in pine trees and the endophytic ability of the fungus on plants, in this paper, we first constructed and utilized a green fluorescent protein (GFP)-tagged H. rhossiliensis HR02 transformant to observe the fungal infection process on B. xylophilus and its colonization on Arabidopsis roots. Then, we compared the fungal parasitism on three species of nematodes with different lifestyles, and we found that the fungal parasitism is correlated with nematode species and stages. The parasitic effect of H. rhossiliensis on adults of B. xylophilus is similar to that on second-stage juveniles (J2) of the root-knot nematode Meloidogyne incognita after 24 h of inoculation, although the virulence of the fungus to second-stage juveniles of M. incognita is stronger than that to those of B. xylophilus and Caenorhabditis elegans. Moreover, the endophytism of H. rhossiliensis was confirmed. By applying an appropriate concentration of H. rhossiliensis conidial suspension (5 × 106 spores/mL) in rhizosphere soil, it was found that the endophytic fungus can promote A. thaliana growth and reproduction, as well as improve host resistance against M. incognita. Our results provide a deeper understanding of the fungus H. rhossiliensis as a promising biocontrol agent against plant-parasitic nematodes.

Variation and stability of rhizosphere bacterial communities of Cucumis crops in association with root-knot nematodes infestation

Introduction

Root-knot nematodes (RKN) disease is a devastating disease in Cucumis crops production. Existing studies have shown that resistant and susceptible crops are enriched with different rhizosphere microorganisms, and microorganisms enriched in resistant crops can antagonize pathogenic bacteria. However, the characteristics of rhizosphere microbial communities of Cucumis crops after RKN infestation remain largely unknown.

Methods

In this study, we compared the changes in rhizosphere bacterial communities between highly RKN-resistant Cucumis metuliferus (cm3) and highly RKN-susceptible Cucumis sativus (cuc) after RKN infection through a pot experiment.

Results

The results showed that the strongest response of rhizosphere bacterial communities of Cucumis crops to RKN infestation occurred during early growth, as evidenced by changes in species diversity and community composition. However, the more stable structure of the rhizosphere bacterial community in cm3 was reflected in less changes in species diversity and community composition after RKN infestation, forming a more complex and positively co-occurrence network than cuc. Moreover, we observed that both cm3 and cuc recruited bacteria after RKN infestation, but the bacteria enriched in cm3 were more abundant including beneficial bacteria Acidobacteria, Nocardioidaceae and Sphingomonadales. In addition, the cuc was enriched with beneficial bacteria Actinobacteria, Bacilli and Cyanobacteria. We also found that more antagonistic bacteria than cuc were screened in cm3 after RKN infestation and most of them were Pseudomonas (Proteobacteria, Pseudomonadaceae), and Proteobacteria were also enriched in cm3 after RKN infestation. We hypothesized that the cooperation between Pseudomonas and the beneficial bacteria in cm3 could inhibit the infestation of RKN.

Discussion

Thus, our results provide valuable insights into the role of rhizosphere bacterial communities on RKN diseases of Cucumis crops, and further studies are needed to clarify the bacterial communities that suppress RKN in Cucumis crops rhizosphere.

Establishment of a CRISPR/Cas9-Mediated Efficient Knockout System of Trichoderma hamatum T21 and Pigment Synthesis PKS Gene Knockout

Trichoderma hamatum is a filamentous fungus that serves as a biological control agent for multiple phytopathogens and as an important resource promising for fungicides. However, the lack of adequate knockout technologies has hindered gene function and biocontrol mechanism research of this species. This study obtained a genome assembly of T. hamatum T21, with a 41.4 Mb genome sequence comprising 8170 genes. Based on genomic information, we established a CRISPR/Cas9 system with dual sgRNAs targets and dual screening markers. CRISPR/Cas9 plasmid and donor DNA recombinant plasmid were constructed for disruption of the Thpyr4 and Thpks1 genes. The result indicates the consistency between phenotypic characterization and molecular identification of the knockout strains. The knockout efficiencies of Thpyr4 and Thpks1 were 100% and 89.1%, respectively. Moreover, sequencing revealed fragment deletions between dual sgRNA target sites or GFP gene insertions presented in knockout strains. The situations were caused by different DNA repair mechanisms, nonhomologous end joining (NHEJ), and homologous recombination (HR). Overall, we have successfully constructed an efficient and convenient CRISPR/Cas9 system in T. hamatum for the first time, which has important scientific significance and application value for studies on functional genomics of Trichoderma and other filamentous fungi.

Chromosome-scale genome assembly-assisted identification of Mi-9 gene in Solanum arcanum accession LA2157, conferring heat-stable resistance to Meloidogyne incognita

Root-knot nematodes (RKNs) are infamous plant pathogens in tomato production, causing considerable losses in agriculture worldwide. Mi-1 is the only commercially available RKN-resistance gene; however, the resistance is inactivated when the soil temperature is over 28 °C. Mi-9 in wild tomato (Solanum arcanum LA2157) has stable resistance to RKNs under high temperature but has not been cloned and applied. In this study, a chromosome-scale genome assembly of S. arcanum LA2157 was constructed through Nanopore and Hi-C sequencing. Based on molecular markers of Mi-9 and comparative genomic analysis, the localization region and candidate Mi-9 genes cluster consisting of seven nucleotide-binding sites and leucine-rich repeat (NBS-LRR) genes were located. Transcriptional expression profiles confirmed that five of the seven candidate genes were expressed in root tissue. Moreover, virus-induced gene silencing of the Sarc_034200 gene resulted in increased susceptibility of S. arcanum LA2157 to Meloidogyne incognita, and genetic transformation of the Sarc_034200 gene in susceptible Solanum pimpinellifolium conferred significant resistance to M. incognita at 25 °C and 30 °C and showed hypersensitive responses at nematode infection sites. This suggested that Sarc_034200 is the Mi-9 gene. In summary, we cloned, confirmed and applied the heat-stable RKN-resistance gene Mi-9, which is of great significance to tomato breeding for nematode resistance.

Comprehensive analysis of the WRKY gene family in Cucumis metuliferus and their expression profile in response to an early stage of root knot nematode infection

Root-knot nematode (RKN) is a major factor that limits the growth and productivity of important Cucumis crops, such as cucumber and melon, which lack RKN-resistance genes in their genome. Cucumis metuliferus is a wild Cucumis species that displays a high degree of RKN-resistance. WRKY transcription factors were involved in plant response to biotic stresses. However, little is known on the function of WRKY genes in response to RKN infection in Cucumis crops. In this study, Cucumis metuliferus 60 WRKY genes (CmWRKY) were identified in the C. metuliferus genome, and their conserved domains were classified into three main groups based on multiple sequence alignment and phylogenetic analysis. Synteny analysis indicated that the WRKY genes were highly conserved in Cucumis crops. Transcriptome data from of C. metuliferus roots inoculated with RKN revealed that 16 CmWRKY genes showed differential expression, of which 13 genes were upregulated and three genes were downregulated, indicating that these CmWRKY genes are important to C. metuliferus response to RKN infection. Two differentially expression CmWRKY genes (CmWRKY10 and CmWRKY28) were selected for further functional analysis. Both CmWRKY genes were localized in nucleus, indicating they may play roles in transcriptional regulation. This study provides a foundation for further research on the function of CmWRKY genes in RKN stress resistance and elucidation of the regulatory mechanism.

Methylorubrum rhodesianum M520 as a biocontrol agent against Meloidogyne incognita (Tylenchida: Heteroderidae) J2s infecting cucumber roots

AIMS

Root-knot nematodes (RKNs) are plant pathogens that cause huge economic losses worldwide. The biological management of RKNs may be a sustainable alternative to chemical control methods. Here, the biocontrol potential of Methylorubrum rhodesianum M520 against the RKN Meloidogyne incognita was investigated to theoretically support its application as a biocontrol agent in field production.

METHODS AND RESULTS

In-vitro assays showed 91.9% mortality of M. incognita second-stage juveniles in the presence of strain M520 and that the hatching rate of M. incognita eggs was 21.7% lower than that of eggs treated with sterile water. In pot experiments, the M520 treatment caused 70.8% reduction in root-knots and increased plant shoot length and stem and root fresh weights, compared to control plant values. In split-root experiments, cucumber roots treated with M520 showed 25.6% decrease in root gall number, compared to that in control roots.

CONCLUSION

M520 has multiple mechanisms against RKNs and might be used as a biocontrol agent against M. incognita in cucumber, laying a foundation for further studying M520 biocontrol against RKNs.

URA3 as a Selectable Marker for Disruption and Functional Assessment of PacC Gene in the Entomopathogenic Fungus Isaria javanica

An effective selection marker is necessary for genetic engineering and functional genomics research in the post-genomic era. Isaria javanica is an important entomopathogenic fungus with a broad host range and prospective biocontrol potentials. Given that no antibiotic marker is available currently in this fungus, developing an effective selection marker is necessary. In this study, by applying overlap PCR and split-marker deletion strategy, combining PEG-mediated protoplasm transformation method, the uridine auxotrophy gene (ura3) in the I. javanica genome was knocked out. Then, using this transformation system, the pH response transcription factor gene (IjpacC) was disrupted successfully. Loss of IjpacC gene results in an obvious decrease in conidial production, but little impact on mycelial growth. The virulence of the ΔIjpacC mutant on caterpillars is similar to that of the wild-type strain. RT-qPCR detection shows that expression level of an acidic-expressed S53 gene (IF1G_06234) in ΔIjpacC mutant is more significantly upregulated than in the wild-type strain during the fungal infection on caterpillars. Our results indicate that a markerless transformation system based upon complementation of uridine auxotrophy is successfully developed in I. javanica, which is useful for exploring gene function and for genetic engineering to enhance biological control potential of the fungus.

The Genome of Fusarium oxysporum f. sp. phaseoli Provides Insight into the Evolution of Genomes and Effectors of Fusarium oxysporum Species

Fusarium oxysporum f. sp. phaseoli, the causal agent of cowpea fusarium wilt, is a serious threat to cowpea production in China. In this study, a sample of cowpea fusarium wilt was identified as Fusarium oxysporum f. sp. phaseoli using the methods of morphological characters and molecular detection. We further reported the first genome assembly for Fusarium oxysporum f. sp. phaseoli, with 53.7 Mb genome sequence comprising 14,694 genes. Comparative genomic analysis among five Fusarium oxysporum genomes showed that four accessory chromosomes in the five Fusarium oxysporum display similar characteristics, with low sequence similarity (55.35%, vs. overall average of 81.76%), low gene density (2.18 genes/10 kb vs. 3.02 genes/Mb) and highly transposable element density (TEs) (15.01/100 kb vs. 4.89/100 kb), indicating that variable accessory chromosomes are the main source of Fusarium oxysporum evolution. We identified a total of 100 Fusarium oxysporum f. sp. phaseoli-specific effectors in the genome and found 13 specific effector genes located in large insertion or deletion regions, suggesting that insertion or deletion events can cause the emergence of species-specific effectors in Fusarium oxysporum. Our genome assembly of Fusarium oxysporum f. sp. phaseoli provides a valuable resource for the study of cowpea fusarium wilt, and the comparative genomic study of Fusarium oxysporum could contribute to the knowledge of genome and effector-associated pathogenicity evolution in Fusarium oxysporum study.

Does Climate Change Increase Crop Water Requirements of Winter Wheat and Summer Maize in the Lower Reaches of the Yellow River Basin?

With increasing water resources stress under climate change, it is of great importance to deeply understand the spatio-temporal variation of crop water requirements and their response to climate change for achieving better water resources management and grain production. However, the quantitative evaluation of climate change impacts on crop water requirements and the identification of determining factors should be further explored to reveal the influencing mechanism and actual effects thoroughly. In this study, the water requirements of winter wheat and summer maize from 1981 to 2019 in the lower reaches of the Yellow River Basin were estimated based on the Penman-Monteith model and crop coefficient method using daily meteorological data. Combined with trends test, sensitivity and contribution analysis, the impacts of different meteorological factors on crop water requirement variation were explored, and the dominant factors were then identified. The results indicated that the temperature increased significantly (a significance level of 0.05 was considered), whereas the sunshine duration, relative humidity and wind speed decreased significantly from 1981 to 2019 in the study area. The total water requirements of winter wheat and summer maize presented a significant decreasing trend (-1.36 mm/a) from 1981 to 2019 with a multi-year average value of 936.7 mm. The crop water requirements of winter wheat was higher than that of summer maize, with multi-year average values of 546.6 mm and 390.1 mm, respectively. In terms of spatial distribution patterns, the crop water requirement in the north was generally higher than that in the south. The water requirements of winter wheat and summer maize were most sensitive to wind speed, and were less sensitive to the minimum temperature and relative humidity. Wind speed was the leading factor of crop water requirement variation with the highest contribution rate of 116.26% among the considered meteorological factors. The results of this study will provide important support for strengthening the capacity to cope with climate change and realizing sustainable utilization of agricultural water resources in the lower reaches of the Yellow River Basin.

Negative regulation of root-knot nematode parasitic behavior by root-derived volatiles of wild relatives of Cucumis metuliferus CM3

Root-knot nematodes (RKN; Meloidogyne spp.) cause a significant decrease in the yield of cucumber crops every year. Cucumis metuliferus is an important wild germplasm that has resistance to RKN in which plant root volatiles are thought to play a role. However, the underlying molecular mechanism is unclear. To investigate it, we used the resistant C. metuliferus line CM3 and the susceptible cucumber line Xintaimici (XTMC). CM3 roots repelled Meloidogyne incognita second-stage larvae (J2s), while the roots of XTMC plants attracted the larvae. CM3 and XTMC were found to contain similar amounts of root volatiles, but many volatiles, including nine hydrocarbons, three alcohols, two aldehydes, two ketones, one ester, and one phenol, were only detected in CM3 roots. It was found that one of these, (methoxymethyl)-benzene, could repel M. incognita, while creosol and (Z)-2-penten-1-ol could attract M. incognita. Interestingly, creosol and (Z)-2-penten-1-ol effectively killed M. incognita at high concentrations. Furthermore, we found that a mixture of CM3 root volatiles increased cucumber resistance to M. incognita. The results provide insights into the interaction between the host and plant-parasitic nematodes in the soil, with some compounds possibly acting as nematode biofumigation, which can be used to manage nematodes.

A Meloidogyne incognita C-type lectin effector targets plant catalases to promote parasitism

Root-knot nematodes, Meloidogyne spp., secrete effectors to modulate plant immune responses and establish a parasitic relationship with host plants. However, the functions and plant targets of C-type lectin (CTL)-like effectors of Meloidogyne incognita remain unknown. Here, we characterized a CTL-like effector of M. incognita, MiCTL1a, and identified its target and role in nematode parasitism. In situ hybridization demonstrated the expression of MiCTL1 in the subventral glands; and in planta, immunolocalization showed its secretion during M. incognita parasitism. Virus-induced gene silencing of the MiCTL1 reduced the infection ability of M. incognita in Nicotiana benthamiana. The ectopic expression in Arabidopsis not only increased susceptibility to M. incognita but also promoted root growth. Yeast two-hybrid and co-immunoprecipitation assays revealed that MiCTL1a interacts with Arabidopsis catalases, which play essential roles in hydrogen peroxide homeostasis. Knockout or overexpression of catalases showed either increased or reduced susceptibility to M. incognita, respectively. Moreover, MiCTL1a not only reduced catalase activity in vitro and in planta but also modulated stress-related gene expressions in Arabidopsis. Our data suggest that MiCTL1a interacts with plant catalases and interferes with catalase activity, allowing M. incognita to establish a parasitic relationship with its host by fine-tuning responses mediated by reactive oxygen species.

Protective Effects of Sodium Para-aminosalicylic Acid on Manganese-Induced Damage in Rat Pancreas

Sodium p-aminosalicylic acid (PAS-Na) has been previously shown to protect the brain from manganese (Mn)-induced toxicity. However, the efficacy of PAS-Na in protecting other organs from Mn toxicity and the mechanisms associated with this protection have yet to be addressed. Therefore, here, we assessed pancreatic damage in response to Mn treatment and the efficacy of PAS-Na in limiting this effect, along with specific mechanisms that mediate PAS-Na's protection. Mn exposure led to increased blood Mn content in dose- and time-dependent manner. Furthermore, subchronic Mn exposure (20 mg/kg for 8 weeks) led to pancreatic damage in a dose-dependent manner. In addition, the elevated Mn levels increased iron and decreased zinc and magnesium content in the pancreas. These effects were noted even 8 weeks after Mn exposure cessation. Mn exposure also affected the levels of amylase, lipase, and inflammatory factors such as tumor necrosis factor (TNF-α) and interleukin-1 β (IL-1β). PAS-Na significantly inhibited the increase in Mn concentration in both blood and pancreas, restored Mn-induced pancreatic damage, reversed the Mn-induced alterations in metal levels, and restored amylase and lipase concentrations. Taken together, we conclude that in rats, PAS-Na shows pharmacological efficacy in protecting the pancreas from Mn-induced damage.

Management of Ralstonia solanacearum in Tomato Using ZnO Nanoparticles Synthesized Through Matricaria chamomilla

Matricaria chamomilla flower extract was used as a biocompatible material for synthesis of zinc oxide nanoparticles (ZnONPs). The synthesized NPs were evaluated for their antibacterial potential in vitro and in vivo against the Gram-negative bacterium Ralstonia solanacearum, which causes devastating bacterial wilt disease in tomato and other crops. Synthesized ZnONPs were further analyzed by UV-visible spectroscopy, Fourier transform infrared spectroscopy, x-ray diffraction, transmission electron microscopy, and scanning electron microscopy with energy-dispersive spectroscopy. The synthesized polydisperse ZnONPs were found to be in the size range of 8.9 to 32.6 nm, and at 18.0 µg ml^-1 exhibited maximum in vitro growth inhibition of the pathogen R. solanacearum. Scanning electron microscopy analysis of affected bacterial cells showed morphological deformation such as disruption of the cell membrane and wall, and the leakage of cell contents. Results of in vivo studies also showed that application of ZnONPs to the artificially inoculated tomato plants with the pathogen R. solanacearum significantly enhanced the plant growth by reducing bacterial soil population and disease severity as compared with the untreated control. Biosynthesized ZnONPs could be an effective approach to control the bacterium R. solanacearum.

High-quality chromosome-level genomes of Cucumis metuliferus and Cucumis melo provide insight into Cucumis genome evolution

Cucumis metuliferus (African horned cucumber), a wild relative of Cucumis sativus (cucumber) and Cucumis melo (melon), displays high-level resistance to several important plant pathogens (e.g., root-knot nematodes and several viruses). Here, we report a chromosome-level genome assembly for C. metuliferus, with a 316 Mb genome sequence comprising 29 039 genes. Phylogenetic analysis of related species in family Cucurbitaceae indicated that the divergence time between C. metuliferus and melon was 17.8 million years ago. Comparisons between the C. metuliferus and melon genomes revealed large structural variations (inversions and translocations >1 Mb) in eight chromosomes of these two species. Gene family comparison showed that C. metuliferus has the largest number of resistance-related nucleotide-binding site leucine-rich repeat (NBS-LRR) genes in Cucurbitaceae. The loss of NBS-LRR loci caused by large insertions or deletions (indels) and pseudogenization caused by small indels explained the loss of NBS-LRR genes in Cucurbitaceae. Population structure analysis suggested that C. metuliferus originated in Zimbabwe, then spread to other southern African regions where it likely underwent similar domestic selection as melon. This C. metuliferus reference sequence will accelerate the understanding of the molecular evolution of resistance-related genes and enhance cucurbit crop improvement efforts.

The root-knot nematode effector MiPDI1 targets a stress-associated protein (SAP) to establish disease in Solanaceae and Arabidopsis

Large amounts of effectors are secreted by the oesophageal glands of plant-parasitic nematodes, but their molecular mode of action remains largely unknown. We characterized a Meloidogyne incognita protein disulphide isomerase (PDI)-like effector protein (MiPDI1) that facilitates nematode parasitism. In situ hybridization showed that MiPDI1 was expressed specifically in the subventral glands of M. incognita. It was significantly upregulated during parasitic stages. Immunolocalization demonstrated MiPDI1 secretion in planta during nematode migration and within the feeding cells. Host-induced silencing of the MiPDI1 gene affected the ability of the nematode to infect the host, whereas MiPDI1 expression in Arabidopsis increased susceptibility to M. incognita, providing evidence for a key role of MiPDI1 in M. incognita parasitism. Yeast two-hybrid, bimolecular fluorescence complementation and coimmunoprecipitation assays showed that MiPDI1 interacted with a tomato stress-associated protein (SlSAP12) orthologous to the redox-regulated AtSAP12, which plays an important role in plant responses to abiotic and biotic stresses. SAP12 silencing or knocking out in Nicotiana benthamiana and Arabidopsis increased susceptibility to M. incognita. Our results suggest that MiPDI1 acts as a pathogenicity factor promoting disease by fine-tuning SAP-mediated responses at the interface of redox signalling, defence and stress acclimation in Solanaceae and Arabidopsis.

Bioactive Secondary Metabolites from Trichoderma spp. against Phytopathogenic Fungi

Phytopathogenic fungi, causing significant economic and production losses, are becoming a serious threat to global food security. Due to an increase in fungal resistance and the hazardous effects of chemical fungicides to human and environmental health, scientists are now engaged to explore alternate non-chemical and ecofriendly management strategies. The use of biocontrol agents and their secondary metabolites (SMs) is one of the potential approaches used today. Trichoderma spp. are well known biocontrol agents used globally. Many Trichoderma species are the most prominent producers of SMs with antimicrobial activity against phytopathogenic fungi. Detailed information about these secondary metabolites, when grouped together, enhances the understanding of their efficient utilization and further exploration of new bioactive compounds for the management of plant pathogenic fungi. The current literature provides the information about SMs of Trichoderma spp. in a different context. In this review, we summarize and group different antifungal SMs of Trichoderma spp. against phytopathogenic fungi along with a comprehensive overview of some aspects related to their chemistry and biosynthesis. Moreover, a brief overview of the biosynthesis pathway, action mechanism, and different approaches for the analysis of SMs and the factors affecting the regulation of SMs in Trichoderma is also discussed.

Discovery of a new antifungal lipopeptaibol from Purpureocillium lilacinum using MALDI-TOF-IMS

Fungi are considered to be rich in biologically active natural products for agricultural and medicinal purposes. The discovery and accurate identification of the bioactive fungal natural products is important for their efficient utilization. During the course of our continuing search for the new natural products from the fungal agents, we found the well-known bio-control fungus Purpureocillium lilacinum showed in vitro activity against Botrytis cinerea, an airborne plant pathogenic fungus causing gray mold disease in many vegetables and fruits. The co-culture of two fungi on agar plate showed that P. lilacinum inhibited the growth of B. cinerea which means P. lilacinum has potential to produce some bioactive secondary metabolites against B. cinerea. In this study, we applied matrix-assisted laser desorption ionization-time of flight mass spectrometry imaging mass spectrometry (MALDI-TOF-IMS), as a fast identification tool, for the discovery of a new antifungal lipopeptaibol (leucinostatin Z) from P. lilacinum against B. cinerea. The planar structure of leucinostatin Z was further established by using the LC-HRESI-MS-MS analysis. MALDI-TOF-IMS is becoming a new approach that allows us to observe the bioactive natural products directly on growth media between the colonies of two fungi, which is faster and more effective than the traditional techniques to discover new bioactive compounds in fungi.

MiMIF-2 Effector of Meloidogyne incognita Exhibited Enzyme Activities and Potential Roles in Plant Salicylic Acid Synthesis

Plant-parasitic nematodes secrete a series of effectors to promote parasitism by modulating host immunity, but the detailed molecular mechanism is ambiguous. Animal parasites secrete macrophage migration inhibitory factor (MIF)-like proteins for evasion of host immune systems, in which their biochemical activities play essential roles. Previous research demonstrated that MiMIF-2 effector was secreted by Meloidogyne incognita and modulated host immunity by interacting with annexins. In this study, we show that MiMIF-2 had tautomerase activity and protected nematodes against H₂O₂ damage. MiMIF-2 expression not only decreased the amount of H₂O₂ generation during nematode infection in Arabidopsis, but also suppressed Bax-induced cell death by inhibiting reactive oxygen species burst in Nicotiana benthamiana. Further, RNA-seq transcriptome analysis and RT-qPCR showed that the expression of some heat-shock proteins was down regulated in MiMIF-2 transgenic Arabidopsis. After treatment with flg22, RNA-seq transcriptome analysis indicated that the differentially expressed genes in MiMIF-2 expressing Arabidopsis were pointed to plant hormone signal transduction, compound metabolism and plant defense. RT-qPCR and metabolomic results confirmed that salicylic acid (SA) related marker genes and SA content were significantly decreased. Our results provide a comprehensive understanding of how MiMIF-2 modulates plant immunity and broaden knowledge of the intricate relationship between M. incognita and host plants.

Effects of α-pinene on the pinewood nematode (Bursaphelenchus xylophilus) and its symbiotic bacteria

The pinewood nematode (PWN), Bursaphelenchus xylophilus, is an important plant-parasitic nematode that can cause severe mortality of pine trees. This PWN-induced harm to plants may be closely related to the abundance and diversity of the symbiotic microorganisms of the parasitic nematode. In this study, nematodes were divided into untreated and antibiotic-treated groups. Nematodes were treated by fumigation with different amounts of α-pinene, and the resultant mortality rates were analyzed statistically. Concentrations of symbiotic bacteria were calculated as colony-forming units per nematode. High-throughput sequencing was used to investigate the bacterial community structure. The results showed that the mortality of nematodes increased slightly with an increasing concentration of α-pinene, and nematodes untreated with antibiotics were more sensitive to α-pinene than those treated with antibiotics. The highest abundance of symbiotic bacteria was obtained via medium and low levels of α-pinene, but for which community diversity was the lowest (Shannon and Simpson indexes). The proportion of Pseudomonas spp. in the symbiotic bacteria of nematodes without antibiotics was relatively high (more than 70%), while that of Stenotrophomonas spp. was low (6%–20%). However, the proportion of Stenotrophomonas spp. was larger than that of Pseudomonas spp in the symbiotic bacteria associated with the antibiotic-treated nematodes. Pseudomonas sp. increased after pinene treatment, whereas Stenotrophomonas spp. decreased. These results indicate that although α-pinene has low toxicity to PWNs over a short time period, α-pinene ultimately influences the abundance and community diversity of the symbiotic bacteria of these nematodes; this influence may potentially disturb the development and reproduction of nematodes in the process of infecting pine trees.

Antibacterial Radicicol Analogues from Pochonia chlamydosporia and Their Biosynthetic Gene Cluster

Chemical investigation of fungus Pochonia chlamydosporia strain 170, derived from rice fermentation sediment samples, afforded seven radicicol analogues, including two new compounds, monocillin VI (1) and monocillin VII (2), and five known compounds, monocillin II (3), monorden D (4), monocillin IV (5), monocillin V (6), and pochonin M (7). The structures of compounds 1-7 were established primarily by analysis of nuclear magnetic resonance data, and the absolute configurations of the secondary alcohol in compounds 1 and 2 were assigned by the modified Mosher method. All seven compounds have modest antibacterial activities, with a minimal inhibitory concentration (MIC) of 25.6 μg/mL for compounds 1 and 3-7 and 51.2 μg/mL for compound 2, on inhibition of the growth of the plant pathogen Xanthomonas campestris (the positive control ampicillin showed a MIC value of 12.8 μg/mL), indicating that the fungus has the potential to control bacterial disease. The biosynthetic gene cluster and putative biosynthetic pathways of these radicicol analogues in the P. chlamydosporia genome were proposed. These findings increase our knowledge of the chemical potential of P. chlamydosporia and may allow us to better utilize the fungus as a biological control agent.

Genomic Characterization Provides New Insights Into the Biosynthesis of the Secondary Metabolite Huperzine a in the Endophyte Colletotrichum gloeosporioides Cg01

A reliable source of Huperzine A (HupA) meets an urgent need due to its wide use in Alzheimer's disease treatment. In this study, we sequenced and characterized the whole genomes of two HupA-producing endophytes, Penicillium polonicum hy4 and Colletotrichum gloeosporioides Cg01, to clarify the mechanism of HupA biosynthesis. The whole genomes of hy4 and Cg01 were 33.92 and 55.77 Mb, respectively. We compared the differentially expressed genes (DEGs) between the induced group (with added extracts of Huperzia serrata) and a control group. We focused on DEGs with similar expression patterns in hy4 and Cg01. The DEGs identified in GO (Gene ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways were primarily located in carbon and nitrogen metabolism and nucleolus, ribosome, and rRNA processing. Furthermore, we analyzed the gene expression for HupA biosynthesis genes proposed in plants, which include lysine decarboxylase (LDC), copper amine oxidase (CAO), polyketides synthases (PKS), etc. Two LDCs, one CAO, and three PKSs in Cg01 were selected as prime candidates for further validation. We found that single candidate biosynthesis-gene knock-out did not influence the HupA production, while both LDC gene knock-out led to increased HupA production. These results reveal that HupA biosynthesis in endophytes might differ from that proposed in plants, and imply that the HupA-biosynthesis genes in endophytic fungi might co-evolve with the plant machinery rather than being acquired through horizontal gene transfer (HGT). Moreover, we analyzed the function of the differentially expressed epigenetic modification genes. HupA production of the histone acetyltransferase (HAT) deletion mutant ΔCgSAS-2 was not changed, while that of the histone methyltransferase (HMT) and histone deacetylase (HDAC) deletion mutants ΔCgClr4, ΔCgClr3, and ΔCgSir2-6 was reduced. Recovery of HupA-biosynthetic ability can be achieved by retro-complementation, demonstrating that HMT and HDACs associated with histone modification are involved in the regulation of HupA biosynthesis in endophytic fungi. This is the first report on epigenetic modification in high value secondary metabolite- producing endophytes. These findings shed new light on HupA biosynthesis and regulation in HupA-producing endophytes and are crucial for industrial production of HupA from fungi.

Inheritance and QTL mapping of cucumber mosaic virus resistance in cucumber (Cucumis Sativus L.)

The commercial yield of cucurbit crops infected with Cucumber mosaic virus (CMV) severely decreases. Chemical treatments against CMV are not effective; therefore, genetic resistance is considered the primary line of defense. Here, we studied resistance to CMV in cucumber inbred line '02245' using a recombinant inbred line (RIL) population generated from a cross between '65G' and '02245' as susceptible and resistant parents, respectively. Genetic analysis revealed that CMV resistance in cucumber is quantitatively inherited. Analysis of the RIL population revealed that a quantitative trait locus (QTL) was found on chromosome 6; named cmv6.1, this QTL was delimited by SSR9-56 and SSR11-177 and explained 31.7% of the phenotypic variation in 2016 and 28.2% in 2017. The marker SSR11-1, which is close to the locus, was tested on 78 different cucumber accessions and found to have an accuracy of 94% in resistant and moderately resistant lines but only 67% in susceptible lines. The mapped QTL was delimited within a region of 1,624.0 kb, and nine genes related to disease resistance were identified. Cloning and alignment of the genomic sequences of these nine genes between '65G' and '02245' revealed that Csa6M133680 had four single-base substitutions within the coding sequences (CDSs) and two single-base substitutions in its 3'-untranslated region, and the other eight genes showed 100% nucleotide sequence identity in their exons. Expression pattern analyses of Csa6M133680 in '65G' and '02245' revealed that the expression levels of Csa6M133680 significantly differed between '65G' and '02245' at 80 h after inoculation with CMV and that the expression in '02245' was 4.4 times greater than that in '65G'. The above results provide insights into the fine mapping and marker-assisted selection in cucumber breeding for CMV resistance.

A Meloidogyne incognita effector MiISE5 suppresses programmed cell death to promote parasitism in host plant

Root-knot nematodes (RKNs) are highly specialized parasites that interact with their host plants using a range of strategies. The esophageal glands are the main places where nematodes synthesize effector proteins, which play central roles in successful invasion. The Meloidogyne incognita effector MiISE5 is exclusively expressed within the subventral esophageal cells and is upregulated during early parasitic stages. In this study, we show that MiISE5 can be secreted to barley cells through infectious hyphae of Magnaporthe oryzae. Transgenic Arabidopsis plants expressing MiISE5 became significantly more susceptible to M. incognita. Inversely, the tobacco rattle virus (TRV)-mediated silence of MiISE5 decreased nematode parasitism. Moreover, transient expression of MiISE5 suppressed cell death caused by Burkholderia glumae in Nicotiana benthamiana. Based on transcriptome analysis of MiISE5 transgenic sample and the wild-type (WT) sample, we obtained 261 DEGs, and the results of GO and KEGG enrichment analysis indicate that MiISE5 can interfere with various metabolic and signaling pathways, especially the JA signaling pathway, to facilitate nematode parasitism. Results from the present study suggest that MiISE5 plays an important role during the early stages of parasitism and provides evidence to decipher the molecular mechanisms underlying the manipulation of host immune defense responses by M. incognita.

Genome and secretome analysis of Pochonia chlamydosporia provide new insight into egg-parasitic mechanisms

Pochonia chlamydosporia infects eggs and females of economically important plant-parasitic nematodes. The fungal isolates parasitizing different nematodes are genetically distinct. To understand their intraspecific genetic differentiation, parasitic mechanisms, and adaptive evolution, we assembled seven putative chromosomes of P. chlamydosporia strain 170 isolated from root-knot nematode eggs (~44 Mb, including 7.19% of transposable elements) and compared them with the genome of the strain 123 (~41 Mb) isolated from cereal cyst nematode. We focus on secretomes of the fungus, which play important roles in pathogenicity and fungus-host/environment interactions, and identified 1,750 secreted proteins, with a high proportion of carboxypeptidases, subtilisins, and chitinases. We analyzed the phylogenies of these genes and predicted new pathogenic molecules. By comparative transcriptome analysis, we found that secreted proteins involved in responses to nutrient stress are mainly comprised of proteases and glycoside hydrolases. Moreover, 32 secreted proteins undergoing positive selection and 71 duplicated gene pairs encoding secreted proteins are identified. Two duplicated pairs encoding secreted glycosyl hydrolases (GH30), which may be related to fungal endophytic process and lost in many insect-pathogenic fungi but exist in nematophagous fungi, are putatively acquired from bacteria by horizontal gene transfer. The results help understanding genetic origins and evolution of parasitism-related genes.

The Novel Secreted Meloidogyne incognita Effector MiISE6 Targets the Host Nucleus and Facilitates Parasitism in Arabidopsis

Meloidogyne incognita is highly specialized parasite that interacts with host plants using a range of strategies. The effectors are synthesized in the esophageal glands and secreted into plant cells through a needle-like stylet during parasitism. In this study, based on RNA-seq and bioinformatics analysis, we predicted 110 putative Meloidogyne incognita effectors that contain nuclear localization signals (NLSs). Combining the Burkholderia glumae-pEDV based screening system with subcellular localization, from 20 randomly selected NLS effector candidates, we identified an effector MiISE6 that can effectively suppress B. glumae-induced cell death in Nicotiana benthamiana, targets to the nuclei of plant cells, and is highly expressed in early parasitic J2 stage. Sequence analysis showed that MiISE6 is a 157-amino acid peptide, with an OGFr_N domain and two NLS motifs. Hybridization in situ verified that MiISE6 is expressed in the subventral esophageal glands. Yeast invertase secretion assay validated the function of the signal peptide harbored in MiISE6. Transgenic Arabidopsis thaliana plants expressing MiISE6 become more susceptible to M. incognita. Inversely, the host-derived RNAi of MiISE6 of the nematode can decrease its parasitism on host. Based on transcriptome analysis of the MiISE6 transgenic Arabidopsis samples and the wild-type samples, we obtained 852 differentially expressed genes (DEGs). Integrating Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses, we found that expression of MiISE6 in Arabidopsis can suppress jasmonate signaling pathway. In addition, the expression of genes related to cell wall modification and the ubiquitination proteasome pathway also have detectable changes in the transgenic plants. Results from the present study suggest that MiISE6 is involved in interaction between nematode-plant, and plays an important role during the early stages of parasitism by interfering multiple signaling pathways of plant. Moreover, we found homologs of MiISE6 in other sedentary nematodes, Meloidogyne hapla and Globodera pallida. Our experimental results provide evidence to decipher the molecular mechanisms underlying the manipulation of host immune defense responses by plant parasitic nematodes, and transcriptome data also provide useful information for further study nematode-plant interactions.

Complete mitogenome of the high ethanol production fungus Fusarium oxysporum Mh2-2

Fusarium spp. are significantly important plant pathogens, and some of them are ethanol-producing strains. During infection and/or ethanol production, Fusarium requires a plenty of energy that is mainly provided by mitochondria. Here we report the first mitogenome from a selected Fusarium oxysporum strain mh2-2 that produces ethanol from glucose and xylose. The size of this mitogenome, 46 kb, is different from the size of any reported Fusarium mitogenome. Our results provide insight into the functions and evolution of mitochondrial genes and genomes.

pBACode: a random-barcode-based high-throughput approach for BAC paired-end sequencing and physical clone mapping

Applications that use Bacterial Artificial Chromosome (BAC) libraries often require paired-end sequences and knowledge of the physical location of each clone in plates. To facilitate obtaining this information in high-throughput, we generated pBACode vectors: a pool of BAC cloning vectors, each with a pair of random barcodes flanking its cloning site. In a pBACode BAC library, the BAC ends and their linked barcodes can be sequenced in bulk. Barcode pairs are determined by sequencing the empty pBACode vectors, which allows BAC ends to be paired according to their barcodes. For physical clone mapping, the barcodes are used as unique markers for their linked genomic sequence. After multi-dimensional pooling of BAC clones, the barcodes are sequenced and deconvoluted to locate each clone. We generated a pBACode library of 94,464 clones for the flounder Paralichthys olivaceus and obtained paired-end sequence from 95.4% of the clones. Incorporating BAC paired-ends into the genome preassembly improved its continuity by over 10-fold. Furthermore, we were able to use the barcodes to map the physical locations of each clone in just 50 pools, with up to 11 808 clones per pool. Our physical clone mapping located 90.2% of BAC clones, enabling targeted characterization of chromosomal rearrangements.

Transcriptome profiling of Cucumis metuliferus infected by Meloidogyne incognita provides new insights into putative defense regulatory network in Cucurbitaceae

Root-knot nematodes (RKN) represent extensive challenges to Cucurbitaceae crops. However, Cucumis metuliferus (Cm) is known to be resistant to Meloidogyne incognita (Mi) infections. Thus, analysis of differentially expressed genes may lead to a comprehensive gene expression profiling of the incompatible Cm-Mi interaction. In this study, the time-course transcriptome of Cm against Mi infection was monitored using RNA-Seq. More than 170000 transcripts were examined in Cm roots, and 2430 genes were subsequently identified as differentially expressed in response to Mi infection. Based on function annotation and orthologs finding, the potential mechanism of transcriptional factor, cytoskeleton, pathogen-related genes and plant hormone were assessed at the transcription level. A comparison of gene expression levels between Mi-infected Cm and cucumber plants revealed that cytoskeleton-related genes are key regulators of Cm resistance to Mi. We herein discuss the dual nature of cytoskeleton-related genes in the susceptibility and resistance of plant hosts to Mi. Our observations provide novel insights into the responses of Cm to Mi at the transcriptome level. The data generated in this study may be useful for elucidating the mechanism underlying resistance to RKNs in cucurbitaceous crops.

Genome-wide analysis of microRNA targeting impacted by SNPs in cucumber genome

Background

microRNAs (miRNAs) are endogenous small RNAs that play important regulatory functions in plant development. Genetic variations in miRNAs sequences or their target-binding sites (microRNA-target interaction sites) can alter miRNA targets in animal and human. Whether these single nucleotide polymorphisms (SNPs) in plant are functional have not yet been determined.

Results

In this study, we constructed leaf, root, and stem-derived small libraries of cucumber (Cucumis sativus) line 9930 (cultivated China-group cucumber) and C. sativus var. hardwickii (wild India group cucumber). A total of 22 conserved miRNA families, nine less-conserved miRNA families, and 49 cucumber-specific miRNAs were identified in both line 9930 and hardwickii. We employed cucumber resequencing data to perform a genome-wide scan for SNPs in cucumber miRNA-target interaction sites, including miRNA mature sequences and miRNA-target binding sites. As a result, we identified a total of 19 SNPs in mature miRNA sequences and 113 SNPs in miRNA-target binding sites with the potential to affect miRNA-target interactions. Furthermore, we experimentally confirmed that these SNPs produced 14 9930-unique targets mRNAs and 15 hardwickii-unique targets mRNA for cucumber miRNAs. This is the first experimental validation of SNPs in miRNA-target interaction sites affecting miRNA-target binding in plants.

Conclusions

Our results indicate that SNPs can alter miRNA function and produce unique miRNA targets in cultivated and wild cucumbers. Therefore, miRNA-related SNPs may have played important in events that led to the agronomic differences between domestic and wild cucumber.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3665-y) contains supplementary material, which is available to authorized users.

The complete mitochondrial genome of huperzine A-producing endophytic fungus Penicillium polonicum

Huperzine A-producing fungus Penicillium polonicum Hy4 (CCTCC No.M2010086) was isolated from Huperzia serrata (Thunb) Trev. The complete mitochondrial genome of P. polonicum is 28 192 bp in length, containing 15 protein-encoding genes, 27 tRNA genes and two rRNA genes. The whole mitogenome is high in AT content (74.40%) and low in GC content (25.60%). The mitochondrial gene order and arrangement of P. polonicum are identical to those of other Penicillium. Phylogenetic analysis based on 14 concatenated protein-encoding genes showed that P. polonicum was close to P. solitum. This study reports the complete mitogenome of P. polonicum for the first time and provides valuable information for further exploration of mitochondrial evolution.

A SIX1 Homolog in Fusarium oxysporum f. sp. conglutinans Is Required for Full Virulence on Cabbage

Fusarium oxysporum is a soil-born fungus that induces wilt and root rot on a variety of plants. F. oxysporum f. sp. conglutinans (Foc) can cause wilt disease on cabbage. This study showed that a homolog of SIX1 protein in the Arabidopsis infecting isolate Fo5176 (Fo5176-SIX1) had four isoforms in the conidia of Foc by proteomic analysis. Thus, we analyzed the roles of protein Foc-SIX1. Gene expression analysis showed that, compared to the expression in mycelia, dramatically altered expression of Foc-SIX1 could be detected after infecting cabbages, and Foc-SIX1 was highly expressed in conidia under axenic culture condition. Furthermore, we knocked out the Foc-SIX1 gene and found that Foc-ΔSIX1 mutants had significantly reduced virulence compared with wild type isolate, and full virulence was restored by complementation of Foc-ΔSIX1 mutants with Foc-SIX1. Thus, we concluded that SIX1 in Foc was required for full virulence on cabbage. We also complemented Foc-ΔSIX1 with SIX1 gene in F. oxysporum f. sp. lycopersici (Fol) and found Foc-ΔSIX1::Fol-SIX1 mutants did not affect the virulence of Foc-ΔSIX1. The results confirmed that Fol-SIX1 was not capable of replacing the role of Foc-SIX1 in Foc on the disease symptom development of cabbage. The roles of Fol-SIX1 on virulence might rely on host specificity.

Double-stranded RNA-mediated interference of dumpy genes in Bursaphelenchus xylophilus by feeding on filamentous fungal transformants

RNA interference (RNAi) is a valuable tool for studying gene function in vivo and provides a functional genomics platform in a wide variety of organisms. The pinewood nematode, Bursaphelenchus xylophilus, is a prominent invasive plant-parasitic nematode and has become a serious worldwide threat to forest ecosystems. Presently, the complete genome sequence of B. xylophilus has been published, and research involving genome-wide functional analyses is likely to increase. In this study, we describe the construction of an effective silencing vector, pDH-RH, which contains a transcriptional unit for a hairpin loop structure. Utilising this vector, double-stranded (ds)RNAs with sequences homologous to the target genes can be expressed in a transformed filamentous fungus via Agrobacterium tumefaciens-mediated transformation technology, and can subsequently induce the knockdown of target gene mRNA expression in B. xylophilus by allowing the nematode to feed on the fungal transformants. Four dumpy genes (Bx-dpy-2, 4, 10 and 11) were used as targets to detect RNAi efficiency. By allowing the nematode to feed on target gene-transformed Fusarium oxysporum strains, target transcripts were knocked down 34-87% compared with those feeding on the wild-type strain as determined by real-time quantitative PCR (RT-qPCR). Morphological RNAi phenotypes were observed, displaying obviously reduced body length; weak dumpy or small (short and thin) body size; or general abnormalities. Moreover, compensatory regulation and non-specific silencing of dpy genes were found in B. xylophilus. Our results indicate that RNAi delivery by feeding in B. xylophilus is a successful technique. This platform may provide a new opportunity for undertaking RNAi-based, genome-wide gene functional studies in vitro in B. xylophilus. Moreover, as B. xylophilus feeds on endophytic fungi when a host has died, RNAi feeding technology will offer the prospect for developing a novel control strategy for the nematode. Furthermore, this platform may also be applicable to other parasitic nematodes that have a facultative, fungivorous habit.

The complete mitochondrial genome of the nematophagous fungus Lecanicillium saksenae

Lecanicillium saksenae CGMCC5329 is a useful biological control agent against plant-parasitic nematodes. The complete mitogenome sequence of L. saksenae is reported for the first time. The mitochondrial genome is 25 919 bp long with 14 typical protein-coding genes, an intronic ORF coding for a putative ribosomal protein (rps3), 2 ribosomal RNA genes and a set of 26 transfer RNA genes. The phylogeny based on 12 protein-coding genes (except the loss of other two genes in Acremonium implicatum) suggests the close phylogenetic relationship between L. saksenae and L. muscarium. Comparative analysis reveals that mitogenome of L. saksenae is 1420 bp larger than L. muscarium, mainly due to the intergenic region between cox2 and trnR. The trnC between cob and cox1 is conserved in the mitogenomes of three nematophagous fungus of Pochonia chlamydosporia, A. implicatum and L. saksenae, but absent in L. muscarium. This study may provide valuable information for further research on mitochondrial evolution of nematophagous fungus and Lecanicillium species.

Enhancement of Palmarumycin C12 and C13 Production by the Endophytic Fungus Berkleasmium sp. Dzf12 in an Aqueous-Organic Solvent System

The endophytic fungus Berkleasmium sp. Dzf12, isolated from Dioscorea zingiberensis, was found to produce palmarumycins C₁₂ and C₁₃ which possess a great variety of biological activities. Seven biocompatible water-immiscible organic solvents including n-dodecane, n-hexadecane, 1-hexadecene, liquid paraffin, dibutyl phthalate, butyl oleate and oleic acid were evaluated to improve palmarumycins C₁₂ and C₁₃ production in suspension culture of Berkleasmium sp. Dzf12. Among the chosen solvents both butyl oleate and liquid paraffin were the most effective to improve palmarumycins C₁₂ and C₁₃ production. The addition of dibutyl phthalate, butyl oleate and oleic acid to the cultures of Berkleasmium sp. Dzf12 significantly enhanced palmarumycin C₁₂ production by adsorbing palmarumycin C₁₂ into the organic phase. When butyl oleate was fed at 5% (v/v) in medium at the beginning of fermentation (day 0), the highest palmarumycin C₁₂ yield (191.6 mg/L) was achieved, about a 34.87-fold increase in comparison with the control (5.3 mg/L). n-Dodecane, 1-hexadecene and liquid paraffin had a great influence on the production of palmarumycin C₁₃. When liquid paraffin was added at 10% (v/v) in medium on day 3 of fermentation, the palmarumycin C₁₃ yield reached a maximum value (134.1 mg/L), which was 4.35-fold that of the control (30.8 mg/L). Application of the aqueous-organic solvent system should be a simple and efficient process strategy for enhancing palmarumycin C₁₂ and C₁₃ production in liquid cultures of the endophytic fungus Berkleasmium sp. Dzf12.

Comparative Proteomics Analyses of Two Races of Fusarium oxysporum f. sp. conglutinans that Differ in Pathogenicity

Fusarium oxysporum is a soil-inhabiting fungus that induces vascular wilt and root rot in a variety of plants. F. oxysporum f. sp. conglutinans (Foc), which comprises two races, can cause wilt disease in cabbage. Compared with race 1 (52557(-TM), R1), race 2 (58385(-TM), R2) exhibits much stronger pathogenicity. Here, we provide the first proteome reference maps for Foc mycelium and conidia and identify 145 proteins with different abundances among the two races. Of these proteins, most of the high-abundance proteins in the R2 mycelium and conidia are involved in carbohydrate, amino acid and ion metabolism, which indicates that these proteins may play important roles in isolate R2's stronger pathogenicity. The expression levels of 20 typical genes demonstrate similarly altered patterns compared to the proteomic analysis. The protein glucanosyltransferase, which is involved in carbohydrate metabolism, was selected for research. We knocked out the corresponding gene (gas1) and found that Foc-∆gas1 significantly reduced growth rate and virulence compared with wild type isolates. These results deepened our understanding of the proteins related to F. oxysporum pathogenicity in cabbage Fusarium wilt and provided new opportunities to control this disease.

Identification of MicroRNAs in Meloidogyne incognita Using Deep Sequencing

MicroRNAs play important regulatory roles in eukaryotic lineages. In this paper, we employed deep sequencing technology to sequence and identify microRNAs in M. incognita genome, which is one of the important plant parasitic nematodes. We identified 102 M. incognita microRNA genes, which can be grouped into 71 nonredundant miRNAs based on mature sequences. Among the 71 miRANs, 27 are known miRNAs and 44 are novel miRNAs. We identified seven miRNA clusters in M. incognita genome. Four of the seven clusters, miR-100/let-7, miR-71-1/miR-2a-1, miR-71-2/miR-2a-2 and miR-279/miR-2b are conserved in other species. We validated the expressions of 5 M. incognita microRNAs, including 3 known microRNAs (miR-71, miR-100b and let-7) and 2 novel microRNAs (NOVEL-1 and NOVEL-2), using RT-PCR. We can detect all 5 microRNAs. The expression levels of four microRNAs obtained using RT-PCR were consistent with those obtained by high-throughput sequencing except for those of let-7. We also examined how M. incognita miRNAs are conserved in four other nematodes species: C. elegans, A. suum, B. malayi and P. pacificus. We found that four microRNAs, miR-100, miR-92, miR-279 and miR-137, exist only in genomes of parasitic nematodes, but do not exist in the genomes of the free living nematode C. elegans. Our research created a unique resource for the research of plant parasitic nematodes. The candidate microRNAs could help elucidate the genomic structure, gene regulation, evolutionary processes, and developmental features of plant parasitic nematodes and nematode-plant interaction.

Qip gene in Fusarium oxysporum is required for normal hyphae morphology and virulence

Ribonucleic acid (RNA)-silencing mechanisms exist in many eukaryotes to regulate a variety of biological processes. The known molecular components are related to Dicers, Argonautes and RNA-dependent RNA polymerases. Previous biochemical studies have also suggested that Qip, with an exonuclease domain, facilitates the conversion of duplex small interfering RNAs into single strands. In our study, the Qip gene in Fusarium oxysporum was disrupted using homologous recombination technology. The deletion of the Qip gene resulted in a decrease in colony growth rates but increased the number of branches. Additionally, the ΔQip mutant had a reduced pathogenicity in cabbage. Our results show Qip gene in F. oxysporum is required for normal hyphae morphology and virulence. The mutant will be useful for elucidating the relationship between the RNA-silencing mechanism and hyphal growth and development in F. oxysporum.

Analysis of the complete mitochondrial genome of Pochonia chlamydosporia suggests a close relationship to the invertebrate-pathogenic fungi in Hypocreales

Background

The fungus Pochonia chlamydosporia parasitizes nematode eggs and has become one of the most promising biological control agents (BCAs) for plant-parasitic nematodes, which are major agricultural pests that cause tremendous economic losses worldwide. The complete mitochondrial (mt) genome is expected to open new avenues for understanding the phylogenetic relationships and evolution of the invertebrate-pathogenic fungi in Hypocreales.

Results

The complete mitogenome sequence of P. chlamydosporia is 25,615 bp in size, containing the 14 typical protein-coding genes, two ribosomal RNA genes, an intronic ORF coding for a putative ribosomal protein (rps3) and a set of 23 transfer RNA genes (trn) which recognize codons for all amino acids. Sequence similarity studies and syntenic gene analyses show that 87.02% and 58.72% of P. chlamydosporia mitogenome sequences match 90.50% of Metarhizium anisopliae sequences and 61.33% of Lecanicillium muscarium sequences with 92.38% and 86.04% identities, respectively. A phylogenetic tree inferred from 14 mt proteins in Pezizomycotina fungi supports that P. chlamydosporia is most closely related to the entomopathogenic fungus M. anisopliae. The invertebrate-pathogenic fungi in Hypocreales cluster together and clearly separate from a cluster comprising plant-pathogenic fungi (Fusarium spp.) and Hypocrea jecorina. A comparison of mitogenome sizes shows that the length of the intergenic regions or the intronic regions is the major size contributor in most of mitogenomes in Sordariomycetes. Evolutionary analysis shows that rps3 is under positive selection, leading to the display of unique evolutionary characteristics in Hypocreales. Moreover, the variability of trn distribution has a clear impact on gene order in mitogenomes. Gene rearrangement analysis shows that operation of transposition drives the rearrangement events in Pezizomycotina, and most events involve in trn position changes, but no rearrangement was found in Clavicipitaceae.

Conclusions

We present the complete annotated mitogenome sequence of P. chlamydosporia. Based on evolutionary and phylogenetic analyses, we have determined the relationships between the invertebrate-pathogenic fungi in Hypocreales. The invertebrate-pathogenic fungi in Hypocreales referred to in this paper form a monophyletic group sharing a most recent common ancestor. Our rps3 and trn gene order results also establish a foundation for further exploration of the evolutionary trajectory of the fungi in Hypocreales.

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-015-0341-8) contains supplementary material, which is available to authorized users.

The complete mitochondrial genome of the nematophagous fungus Acremonium implicatum

The complete mitochondrial genome of the nematophagous fungus Acremonium implicatum is reported for the first time. The genome is concatenated with 22,367 bp in length, encoding 13 protein-coding genes, 2 ribosomal RNA genes and a set of 17 transfer RNA genes. The synteny analysis reveals that 50.35% of A. implicatum mitochondrial sequences matched to 48.21% of Acremonium chrysogenum mitochondrial sequences with 85.68% identity. Two proteins of cox3 and nad6, as well as seven tRNAs are lost in A. implicatum mitogenome compared to A. chrysogenum mitogenome. The gene orders in A. implicatum and A. chrysogenum mitogenome is different, which is mainly due to the location of nad4 and cox2. In addition, one transposition event related to tRNAs is identified in these two mitogenomes. This study may provide valuable mitochondrial information for research on A. implicatum and facilitate the study of mitochondrial evolution.

Development of polymorphic microsatellites for Meloidogyne incognita, through screening predicted microsatellite loci based on genome sequence

Microsatellites are extensively distributed in the eukaryotic genome, and they are widely used for their high polymorphism and accessibility. The microsatellites in M. incognita, a worldwide agriculture pest, are inadequate for diversity research. A repertoire of 1620 microsatellites appeared appropriate to design primer as markers were identified based on the M. incognita genome. 120 loci were chosen as candidate, from which 88 microsatellites were characterized. Finally, we found 13 polymorphic microsatellites with 2 to 23 alleles in a survey of three nematode populations in China, while other positive loci were monomorphic. These new molecular markers afford to genetic diversity analysis in M. incognita population of poorly investigation. Furthermore, the predicted microsatellites have potential values for other plant parasitic nematodes.

Mapping and analysis of a novel candidate Fusarium wilt resistance gene FOC1 in Brassica oleracea

BACKGROUND

Cabbage Fusarium wilt is a major disease worldwide that can cause severe yield loss in cabbage (Brassica olerecea). Although markers linked to the resistance gene FOC1 have been identified, no candidate gene for it has been determined so far. In this study, we report the fine mapping and analysis of a candidate gene for FOC1 using a double haploid (DH) population with 160 lines and a F2 population of 4000 individuals derived from the same parental lines.

RESULTS

We confirmed that the resistance to Fusarium wilt was controlled by a single dominant gene based on the resistance segregation ratio of the two populations. Using InDel primers designed from whole-genome re-sequencing data for the two parental lines (the resistant inbred-line 99-77 and the highly susceptible line 99-91) and the DH population, we mapped the resistance gene to a 382-kb genomic region on chromosome C06. Using the F2 population, we narrowed the region to an 84-kb interval that harbored ten genes, including four probable resistance genes (R genes): Bol037156, Bol037157, Bol037158 and Bol037161 according to the gene annotations from BRAD, the genomic database for B. oleracea. After correcting the model of the these genes, we re-predicted two R genes in the target region: re-Bol037156 and re-Bol0371578. The latter was excluded after we compared the two genes' sequences between ten resistant materials and ten susceptible materials. For re-Bol037156, we found high identity among the sequences of the resistant lines, while among the susceptible lines, there were two types of InDels (a 1-bp insertion and a 10-bp deletion), each of which caused a frameshift and terminating mutation in the cDNA sequences. Further sequence analysis of the two InDel loci from 80 lines (40 resistant and 40 susceptible) also showed that all 40 R lines had no InDel mutation while 39 out of 40 S lines matched the two types of loci. Thus re-Bol037156 was identified as a likely candidate gene for FOC1 in cabbage.

CONCLUSIONS

This work may lay the foundation for marker-assisted selection as well as for further function analysis of the FOC1 gene.

Analysis of the defence-related mechanism in cucumber seedlings in relation to root colonization by nonpathogenic Fusarium oxysporum CS-20

A defence response can be induced by nonpathogenic Fusarium oxysporum CS-20 in several crops, but the molecular mechanism has not been clearly demonstrated. In the present study, we analysed the defence mechanism of a susceptible cucumber cultivar (Cucumis sativus L. 9930) against a pathogen (F. oxysporum f. sp. cucumerinum) through the root precolonization of CS-20. A challenge inoculation assay indicated that the disease severity index (DSI) was reduced, ranging from 18.83 to 61.67 in comparison with the pathogen control. Root colonization analysis indicated that CS-20 clearly did not appear to influence the growth of cucumber seedlings. Quantitative reverse-transcriptase polymerase chain reaction (qRT-PCR) revealed that CS-20-mediated defence response was activated by PR3, LOX1 and PAL1 and the pathogen-mediated resistance response was regulated by PR1 and PR3. Moreover, both nonpathogenic and pathogenic F. oxysporum were able to upregulate NPR1 expression. In contrast to a pathogen, CS-20 can activate the Ca(2+) /CaM signal transduction pathway, and the gene expression of both CsCam7 and CsCam12 increased significantly. The gene expression analysis indicated that CS-20 strongly enhanced the expression of PR3, LOX1, PAL1, NPR1, CsCam7 and CsCam12 after inoculation. Overall, the defence response induced by CS-20 can be controlled by multiple genes in the cucumber plant.

A genomic variation map provides insights into the genetic basis of cucumber domestication and diversity

Most fruits in our daily diet are the products of domestication and breeding. Here we report a map of genome variation for a major fruit that encompasses ~3.6 million variants, generated by deep resequencing of 115 cucumber lines sampled from 3,342 accessions worldwide. Comparative analysis suggests that fruit crops underwent narrower bottlenecks during domestication than grain crops. We identified 112 putative domestication sweeps; 1 of these regions contains a gene involved in the loss of bitterness in fruits, an essential domestication trait of cucumber. We also investigated the genomic basis of divergence among the cultivated populations and discovered a natural genetic variant in a β-carotene hydroxylase gene that could be used to breed cucumbers with enhanced nutritional value. The genomic history of cucumber evolution uncovered here provides the basis for future genomics-enabled breeding.