Variability of United States isolates of Macrophomina phaseolina based on simple sequence repeats and cross genus transferability to related genera within botryosphaeriaceae

Molecular interactions between the soilborne pathogenic fungus Macrophomina phaseolina and its host plants

Mentioned for the first time in an article 1971, the occurrence of the term "Macrophomina phaseolina" has experienced a steep increase in the scientific literature over the past 15 years. Concurrently, incidences of M. phaseolina-caused crop diseases have been getting more frequent. The high levels of diversity and plasticity observed for M. phasolina genomes along with a rich equipment of plant cell wall degrading enzymes, secondary metabolites and putative virulence effectors as well as the unusual longevity of microsclerotia, their asexual reproduction structures, make this pathogen very difficult to control and crop protection against it very challenging. During the past years several studies have emerged reporting on host defense measures against M. phaseolina, as well as mechanisms of pathogenicity employed by this fungal pathogen. While most of these studies have been performed in crop systems, such as soybean or sesame, recently interactions of M. phaseolina with the model plant Arabidopsis thaliana have been described. Collectively, results from various studies are hinting at a complex infection cycle of M. phaseolina, which exhibits an early biotrophic phase and switches to necrotrophy at later time points during the infection process. Consequently, responses of the hosts are complex and seem coordinated by multiple defense-associated phytohormones. However, at this point no robust and strong host defense mechanism against M. phaseolina has been described.

Comparative analysis of nine Tilletia indica genomes for the development of novel microsatellite markers for genetic diversity and population structure analysis

Karnal bunt (KB; Tilletia indica) is the prime quarantine concern for quality wheat production throughout the world. The most effective approach to dealing with this biotic stress is to breed KB-resistant wheat varieties, which warrants a better understanding of T. indica genome architecture. In India, the North Western Plain Zone is the prime hot spot for KB disease, but only limited efforts have been made to decipher T. indica diversity at the genomic level. Microsatellites offer a powerful and robust typing system for the characterization and genetic diversity assessment of plant pathogens. At present, inadequate information is available with respect to the development of genome-derived markers for revealing genetic variability in T. indica populations. In current research, nine complete genome sequences of T. indica (PSWKBGH_1, PSWKBGH_2, PSWKBGD_1_3, RAKB_UP_1, TiK_1, Tik, DAOMC236408, DAOMC236414, and DAOMC236416) that exist in the public domain were explored to know the dynamic distribution of microsatellites. Comparative genome analysis revealed a high level of relative abundance and relative density of microsatellites in the PSWKBGH_1 genome in contrast to other genomes. No significant correlation between microsatellite distribution for GC content and genome size was established. All the genomes showed the dominance of tri-nucleotide motifs, followed by mono-, di-, tetra-, hexa-, and penta-nucleotide motifs. Out of 50 tested markers, 36 showed successful amplification in T. indica isolates and produced 52 different alleles. A PCR assay along with analysis of the polymorphic information content (PIC) revealed 10 markers as neutral and polymorphic loci (PIC 0.37). The identified polymorphic SSR loci grouped a geographically distinct T. indica population of 50 isolates representing seven Indian regions (Jammu, Himachal Pradesh, Punjab, Haryana, Uttarakhand, Uttar Pradesh, and Rajasthan) into four distinct clusters. The results of the analysis of molecular variance identified 94% genetic variation within the population and 6% among the population. Structure analysis also confirmed the existence of four genetically diverse groups containing admixtures of T. indica isolates across populations. In nutshell, the current study was successful in identifying novel, neutral and polymorphic microsatellite markers that will be valuable in offering deep insight into the evolutionary relationship and dynamics of the T. indica population for devising effective KB management strategies in wheat.

Population genomic analysis reveals geographic structure and climatic diversification for Macrophomina phaseolina isolated from soybean and dry bean across the United States, Puerto Rico, and Colombia

Macrophomina phaseolina causes charcoal rot, which can significantly reduce yield and seed quality of soybean and dry bean resulting from primarily environmental stressors. Although charcoal rot has been recognized as a warm climate-driven disease of increasing concern under global climate change, knowledge regarding population genetics and climatic variables contributing to the genetic diversity of M. phaseolina is limited. This study conducted genome sequencing for 95 M. phaseolina isolates from soybean and dry bean across the continental United States, Puerto Rico, and Colombia. Inference on the population structure using 76,981 single nucleotide polymorphisms (SNPs) revealed that the isolates exhibited a discrete genetic clustering at the continental level and a continuous genetic differentiation regionally. A majority of isolates from the United States (96%) grouped in a clade with a predominantly clonal genetic structure, while 88% of Puerto Rican and Colombian isolates from dry bean were assigned to a separate clade with higher genetic diversity. A redundancy analysis (RDA) was used to estimate the contributions of climate and spatial structure to genomic variation (11,421 unlinked SNPs). Climate significantly contributed to genomic variation at a continental level with temperature seasonality explaining the most variation while precipitation of warmest quarter explaining the most when spatial structure was accounted for. The loci significantly associated with multivariate climate were found closely to the genes related to fungal stress responses, including transmembrane transport, glycoside hydrolase activity and a heat-shock protein, which may mediate climatic adaptation for M. phaseolina. On the contrary, limited genome-wide differentiation among populations by hosts was observed. These findings highlight the importance of population genetics and identify candidate genes of M. phaseolina that can be used to elucidate the molecular mechanisms that underly climatic adaptation to the changing climate.

Genetic Diversity of Phymatotrichopsis omnivora Based on Mating Type and Microsatellite Markers Reveals Heterothallic Mating System

Phymatotrichopsis omnivora is a member of Pezizomycetes and causes root rot disease on a broad range of dicotyledonous plants. Using recently generated draft genome sequence data from four P. omnivora isolates, we developed simple sequence repeat (SSR) markers and identified both mating type genes (MAT1-1-1 and MAT1-2-1) in this fungus. To understand the genetic diversity of P. omnivora isolates (n = 43) and spore mats (n = 29) collected from four locations (Oklahoma, Texas, Arizona, and Mexico) and four host crops (cotton, alfalfa, peach, and soybean), we applied 24 SSR markers and showed that of the 72 P. omnivora isolates and spore mats tested, 41 were distinct genotypes. Furthermore, the developed SSR markers did not show cross-transferability to other close relatives of P. omnivora in the class Pezizomycetes. A multiplex PCR detecting both mating type idiomorphs and a reference gene (TUB2) was developed to screen P. omnivora isolates. Based on the dataset we tested, P. omnivora is a heterothallic fungus with both mating types present in the United States in a ratio close to 1:1. We tested P. omnivora spore mats obtained from spatially distinct disease rings that developed in a center-pivot alfalfa field and showed that both mating types can be present not only in the same field but also within a single spore mat. This study shows that P. omnivora has the genetic toolkit for generating sexually diverse progeny, providing impetus for future studies that focus on identifying sexual morphs in nature.

A novel Arabidopsis pathosystem reveals cooperation of multiple hormonal response-pathways in host resistance against the global crop destroyer Macrophomina phaseolina

Dubbed as a "global destroyer of crops", the soil-borne fungus Macrophomina phaseolina (Mp) infects more than 500 plant species including many economically important cash crops. Host defenses against infection by this pathogen are poorly understood. We established interactions between Mp and Arabidopsis thaliana (Arabidopsis) as a model system to quantitatively assess host factors affecting the outcome of Mp infections. Using agar plate-based infection assays with different Arabidopsis genotypes, we found signaling mechanisms dependent on the plant hormones ethylene, jasmonic acid and salicylic acid to control host defense against this pathogen. By profiling host transcripts in Mp-infected roots of the wild-type Arabidopsis accession Col-0 and ein2/jar1, an ethylene/jasmonic acid-signaling deficient mutant that exhibits enhanced susceptibility to this pathogen, we identified hundreds of genes potentially contributing to a diverse array of defense responses, which seem coordinated by complex interplay between multiple hormonal response-pathways. Our results establish Mp/Arabidopsis interactions as a useful model pathosystem, allowing for application of the vast genomics-related resources of this versatile model plant to the systematic investigation of previously understudied host defenses against a major crop plant pathogen.

Soybean-Macrophomina phaseolina-Specific Interactions and Identification of a Novel Source of Resistance

Charcoal rot, caused by the fungus Macrophomina phaseolina, is an economically important disease of soybean (Glycine max) worldwide. Objectives of the present research were to (i) study the genetic and pathogenic diversity in a collection of M. phaseolina isolates from Argentina and Paraguay and (ii) develop an improved in vitro phenotyping method to evaluate disease response of soybean genotypes to M. phaseolina isolates. Cluster analysis showed no clear association among simple sequence repeat profiles, year of collection, pathogenicity, and geographical origin of the isolates from Argentina and Paraguay. Subsequently, the response of four soybean genotypes against seven M. phaseolina isolates was evaluated in the field and the results were confirmed using the in vitro assay developed. This assay, which is based on root disease development on soybean seedlings, allowed the detection of a differential level of aggressiveness among the isolates on four soybean genotypes. The results suggest the existence of specific interactions among soybean genotypes and M. phaseolina isolates. In addition, cultivar Munasqa RR showed a superior response against M. phaseolina compared with DT 97-4290 (moderately resistant), thus becoming a novel source of resistance to charcoal rot.

Development of Molecular Methods to Detect Macrophomina phaseolina from Strawberry Plants and Soil

Macrophomina phaseolina is a broad-host-range fungus that shows some degree of host preference on strawberry, and causes symptoms that include crown rot and root rot. Recently, this pathogen has affected strawberry production as fumigation practices have changed, leaving many growers in California and around the world in need of accurate, rapid diagnostic tools for M. phaseolina in soil and infected plants. This study uses next-generation sequencing and comparative genomics to identify a locus that is unique to isolates within a main genotype shared by a majority of isolates that infect strawberry. This locus was used to develop a quantitative single-tube nested TaqMan polymerase chain reaction assay which is able to quantify as little as 2 to 3 microsclerotia/g of soil with 100% genotype specificity. An isothermal assay using recombinase polymerase amplification was developed from the same locus and has been validated on over 200 infected strawberry plants with a diagnostic sensitivity of 93% and a diagnostic specificity of 99%. Together, this work demonstrates the value of using new approaches to identify loci for detection and provides valuable diagnostic tools that can be used to monitor soil and strawberry plant samples for M. phaseolina.

Next Generation Sequencing-Based Molecular Marker Development: A Case Study in Betula Alnoides

Betula alnoides is a fast-growing valuable indigenous tree species with multiple uses in the tropical and warm subtropical regions in South-East Asia and southern China. It has been proved to be tetraploid in most parts of its distribution in China. In the present study, next generation sequencing (NGS) technology was applied to develop numerous SSR markers for B. alnoides, and 64,376 contig sequences of 106,452 clean reads containing 164,357 candidate SSR loci were obtained. Among the derived SSR repeats, mono-nucleotide was the main type (77.05%), followed by di- (10.18%), tetra- (6.12%), tri- (3.56%), penta- (2.14%) and hexa-nucleotide (0.95%). The short nucleotide sequence repeats accounted for 90.79%. Among the 291 repeat motifs, AG/CT (46.33%) and AT/AT (44.15%) were the most common di-nucleotide repeats, while AAT/ATT (48.98%) was the most common tri-nucleotide repeats. A total of 2549 primer sets were designed from the identified putative SSR regions of which 900 were randomly selected for evaluation of amplification successfulness and detection of polymorphism if amplified successfully. Three hundred and ten polymorphic markers were obtained through testing with 24 individuals from B. alnoides natural forest in Jingxi County, Guangxi, China. The number of alleles (NA) of each marker ranged from 2 to 19 with a mean of 5.14. The observed (HO) and expected (HE) heterozygosities varied from 0.04 to 1.00 and 0.04 to 0.92 with their means being 0.64 and 0.57, respectively. Shannon-Wiener diversity index (I) ranged from 0.10 to 2.68 with a mean of 1.12. Cross-species transferability was further examined for 96 pairs of SSR primers randomly selected, and it was found that 48.96⁻84.38% of the primer pairs could successfully amplify each of six related Betula species. The obtained SSR markers can be used to study population genetics and molecular marker assisted breeding, particularly genome-wide association study of these species in the future.

Genetic Architecture of Charcoal Rot (Macrophomina phaseolina) Resistance in Soybean Revealed Using a Diverse Panel

Charcoal rot (CR) disease caused by Macrophomina phaseolina is responsible for significant yield losses in soybean production. Among the methods available for controlling this disease, breeding for resistance is the most promising. Progress in breeding efforts has been slow due to the insufficient information available on the genetic mechanisms related to resistance. Genome-wide association studies (GWAS) enable unraveling the genetic architecture of resistance and identification of causal genes. The aims of this study were to identify new sources of resistance to CR in a collection of 459 diverse plant introductions from the USDA Soybean Germplasm Core Collection using field and greenhouse screenings, and to conduct GWAS to identify candidate genes and associated molecular markers. New sources for CR resistance were identified from both field and greenhouse screening from maturity groups I, II, and III. Five significant single nucleotide polymorphism (SNP) and putative candidate genes related to abiotic and biotic stress responses are reported from the field screening; while greenhouse screening revealed eight loci associated with eight candidate gene families, all associated with functions controlling plant defense response. No overlap of markers or genes was observed between field and greenhouse screenings suggesting a complex molecular mechanism underlying resistance to CR in soybean with varied response to different environments; but our findings provide useful information for advancing breeding for CR resistance as well as the genetic mechanism of resistance.

Molecular Identification and Genetic Characterization of Macrophomina phaseolina Strains Causing Pathogenicity on Sunflower and Chickpea

Macrophomina phaseolina is the most devastating pathogen which causes charcoal rot and root rot diseases in various economically important crops. Three strains M. phaseolina 1156, M. phaseolina 1160, and M. phaseolina PCMC/F1 were tested for their virulence on sunflower (Helianthus annuus L.) and chickpea (Cicer arietinum L.). The strains showed high virulence on both hosts with a disease score of 2 on chickpea and sunflower. The strains also increased the hydrogen per oxide (H2O2) content by 1.4- to 1.6-fold in root as well as shoot of chickpea and sunflower. A significant increase in antioxidant enzymes was observed in fungal infected plants which indicated prevalence of oxidative stress during pathogen propagation. The M. phaseolina strains also produced hydrolytic enzymes such as lipase, amylase, and protease with solubilization zone of 5-43 mm, 5-45 mm, and 12-35 mm, respectively. The M. phaseolina strains were identified by 18S rRNA and analyzed for genetic diversity by using random amplified polymorphic DNA (RAPD) markers. The findings based on RAPD markers and 18S rRNA sequence analysis clearly indicate genetic variation among the strains collected from different hosts. The genetically diverse strains were found to be pathogenic to sunflower and chickpea.

Comparative analysis of microsatellites in five different antagonistic Trichoderma species for diversity assessment

Microsatellites provide an ideal molecular markers system to screen, characterize and evaluate genetic diversity of several fungal species. Currently, there is very limited information on the genetic diversity of antagonistic Trichoderma species as determined using a range of molecular markers. In this study, expressed and whole genome sequences available in public database were used to investigate the occurrence, relative abundance and relative density of SSRs in five different antagonistic Trichoderma species: Trichoderma atroviride, T. harzianum, T. reesei, T. virens and T. asperellum. Fifteen SSRs loci were used to evaluate genetic diversity of twenty isolates of Trichoderma spp. from different geographical regions of India. Results indicated that relative abundance and relative density of SSRs were higher in T. asperellum followed by T. reesei and T. atroviride. Tri-nucleotide repeats (80.2%) were invariably the most abundant in all species. The abundance and relative density of SSRs were not influenced by the genome sizes and GC content. Out of eighteen primer sets, only 15 primer pairs showed successful amplification in all the test species. A total of 24 alleles were detected and five loci were highly informative with polymorphism information content values greater than 0.40, these markers provide useful information on genetic diversity and population genetic structure, which, in turn, can exploit for establishing conservation strategy for antagonistic Trichoderma isolates.

Identification and characterization of microsatellite from Alternaria brassicicola to assess cross-species transferability and utility as a diagnostic marker

Alternaria blight caused by Alternaria brassicicola (Schwein.) Wiltshire and A. brassicae (Berk.) Sacc., is one of the most important disease of rapeseed-mustard, characterized by the formation of spots on leaves, stem, and siliquae with premature defoliation and stunting of growth. These two species are very difficult to differentiate based on disease symptoms or spore morphology. Therefore, the aim of present investigation was to identify and characterize transferable microsatellite loci from A. brassicicola to A. brassicae for the development of diagnostic marker. A total of 8,457 microsatellites were identified from transcript sequences of A. brassicicola. The average density of microsatellites was one microsatellite per 1.94 kb of transcript sequence screened. The most frequent repeat was tri-nucleotide (74.03 %), whereas penta-nucleotide (1.14 %) was least frequent. Among amino acids, arginine (13.11 %) showed maximum abundance followed by lysine (10.11 %). A total of 32 alleles were obtained across the 31 microsatellite loci for the ten isolates of A. brassicicola. In cross-species amplifications, 5 of the 31 markers amplified the corresponding microsatellite regions in twenty isolates of A. brassicae and showed monomorphic banding pattern. Microsatellite locus ABS28 was highly specific for A. brassicicola, as no amplification was observed from twenty-nine other closely related taxa. Primer set, ABS28F/ABS28R, amplified a specific amplicon of 380 bp from all A. brassicicola isolates. Standard curves were generated for A. brassicicola isolate using SYBR Green I fluorescent dye for detection of amplification in real-time PCR assay. The lowest detection limit of assay was 0.01 ng. Thus, the primer set can be used as diagnostic marker to discriminate and diagnose A. brassicicola from synchronously occurring fungus, A. brassicae associated with rapeseed and mustard.

Characterization of expressed sequence tag-derived simple sequence repeat markers for Aspergillus flavus: emphasis on variability of isolates from the southern United States

Simple sequence repeat (SSR) markers were developed from Aspergillus flavus expressed sequence tag (EST) database to conduct an analysis of genetic relationships of Aspergillus isolates from numerous host species and geographical regions, but primarily from the United States. Twenty-nine primers were designed from 362 tri-nucleotide EST-SSR sequences. Eighteen polymorphic loci were used to genotype 96 Aspergillus species isolates. The number of alleles detected per locus ranged from 2 to 24 with a mean of 8.2 alleles. Haploid diversity ranged from 0.28 to 0.91. Genetic distance matrix was used to perform principal coordinates analysis (PCA) and to generate dendrograms using unweighted pair group method with arithmetic mean (UPGMA). Two principal coordinates explained more than 75 % of the total variation among the isolates. One clade was identified for A. flavus isolates (n = 87) with the other Aspergillus species (n = 7) using PCA, but five distinct clusters were present when the others taxa were excluded from the analysis. Six groups were noted when the EST-SSR data were compared using UPGMA. However, the latter PCA or UPGMA comparison resulted in no direct associations with host species, geographical region or aflatoxin production. Furthermore, there was no direct correlation to visible morphological features such as sclerotial types. The isolates from Mississippi Delta region, which contained the largest percentage of isolates, did not show any unusual clustering except for isolates K32, K55, and 199. Further studies of these three isolates are warranted to evaluate their pathogenicity, aflatoxin production potential, additional gene sequences (e.g., RPB2), and morphological comparisons.

Diversity and specificity of microsatellites within Aspergillus section Fumigati

BACKGROUND

Microsatellites (or short tandem repeats, STRs) are the genetic markers of choice for studying Aspergillus fumigatus molecular epidemiology due to its reproducibility and high discrimination power. However, the specificity of these markers must be investigated in a group of isolates from closely related species. The aim of this work was to test a microsatellite-based PCR multiplex previously designed for A. fumigatus in a set of species belonging to section Fumigati, namely Aspergillus fumigatiaffinis, Aspergillus lentulus, Aspergillus novofumigatus, Aspergillus unilateralis, Aspergillus viridinutans, Neosartorya fischeri, Neosartorya hiratsukae, Neosartorya pseudofischeri and Neosartorya udagawae.

RESULTS

The reference A. fumigatus strain ATCC 46645 was easily genotyped in standard conditions showing a final electrophoretic profile of 8 expected peaks corresponding to each microsatellite locus. Inversely, no peaks were observed for all other species from section Fumigati, with an exception for marker MC6b in A. unilateralis. By screening the genome sequence of Neosartorya fischeri NRRL 181, the results showed that MC3, MC6a and MC7 might be employed for N. fischeri genotyping since these markers present several repeats of each motif. The accumulation of insertions and deletions was frequently observed in the genomic regions surrounding the microsatellites, including those where the A. fumigatus primers are located. The amplification of microsatellite markers in less stringent amplification conditions resulted in a distinct electrophoretic profile for species within section Fumigati.

CONCLUSIONS

Therefore, the microsatellite-based PCR multiplex allow simple identification of A. fumigatus and, with a slight modification of temperature conditions, it also allows discriminating other pathogenic species within section Fumigati, particularly A. fumigatiaffinis, N. fischeri and N. udagawae.

Genetic diversity and effect of temperature and pH on the growth of Macrophomina phaseolina isolates from sunflower fields in Hungary

The effects of temperature and pH on the growth of 45 Hungarian Macrophomina phaseolina isolates from different locations and hosts were compared on the basis of their genetic diversity. One Spanish and two Serbian isolates were also included in the experiment. The most favourable temperature regimes for the development of the isolates ranged between 25 and 35 °C. The optimal pH for the pathogen varied between 4.0 and 6.0, but growth was observed on potato dextrose agar even at pH values of 3.0, 7.0 and 8.0. RAPD analysis with 13 different primer pairs generated 148 unambiguous bands. RFLP analysis involving 8 different restriction endonucleases was performed on a 1550 bp fragment of the rDNA region containing internal transcribed spacers (ITS1, ITS2), the 5.8S rDNA and part of the 25S rDNA. The greatest genetic distance values were obtained for three isolates, two from Hungary and one from Spain, which had similar values, but were quite distinct from all the others. A strong positive correlation was observed between the genetic distances and the growth parameters measured at various temperatures, and between the geographical data and the growth data sets at different pH values, but the correlation was less strong in the latter case. While Hungarian M. phaseolina populations are thought to reproduce clonally, the present results indicate the coexistence of different haplotypes in this area, and besides the geographical dominance of a given haplotype it was found that a closer genetic relationship might exist between spatially distinct haplotypes.