First Report of Leptosphaerulina australis Isolated from Soil in Korea

Research Article

Abstract

The fungal strain KNU16-004 was isolated from a field soil sample collected in Seoul. The isolate was identified as Leptosphaerulina australis based on morphological characterization and phylogenetic analysis using the internal transcribed spacer (ITS), large subunit (LSU) rDNA regions, and β-tubulin (Tub2). This is the first report of Leptosphaerulina australis in Korea.

Keyword



INTRODUCTION

The genus Leptosphaerulina was first introduced by McAlpine in 1902, by designating Leptosphaerulina australis as the type species [1]. In 1961, Graham and Luttrell [2] described three existing species and three new species of Leptosphaerulina found on forage plants; the pathogenic Leptosphaerulina arachidicola, L. briosiana, and L. trifolii and the saprobic L. australis, L. argentinensis, and L. americana. However, Irwin and Davis [3] reported that L. argentinensis is pathogenic to Stylosanthes guianensis. Leptosphaerulina species have ostiolate, papillate ascomata, which are immersed to erumpent, and bitunicate asci that are distinctly saccate. The ascospores are oblong to cylindrical, generally muriform, and mostly hyaline, but become brown at maturity [2, 4–6]. The type species of the genus, L. australis, was found on apricot leaves (Prunus armeniaca L.) and has also been reported to be associated with Dolichos Poa, Lolium and Vitis spp. [1, 2]. In the United States, L. australis has been shown to colonize many species from the following genera: Agrostis Festuca, Ligustrum, Lolium, Panicum, Poa, Rosa, and Trifolium [2, 7, 8].

In the present study, a fungal strain, designated KNU16-004, which was isolated from field soil sample, was identified as L. australis based on morphological characterization and phylogenetic analysis using a concatenated dataset of the partial ribosomal large subunit (LSU), internal transcribed spacer (ITS) region, and β-tubulin (Tub2) gene sequences.

MATERIALS and METHODS

Sampling and fungi isolation

Fungal strain KNU16-004 was isolated from a field soil sample collected in Seoul, Korea (N 37°30'59.5", E 127°07'12.6"). The soil sample was collected from the ground, air-dried, and then stored in a plastic bag at 4°C until analysis. The fungus was isolated as previously described [9]. One gram of the soil sample was suspended in 10 mL of sterile distilled water, and the suspension was then vortexed and serially diluted. One hundred microliters of each dilution were spread onto potato dextrose agar (PDA; Difco, Detroit, MI, USA) plates, which were incubated at 25°C for 3~7 days until fungal colony growth was observed. Single colonies from these plates were purified by transferring them onto fresh plates and incubating the plates under the same conditions. The pure cultures were routinely cultured on PDA agar at 25°C and maintained as a glycerol suspension (20%, w/v) at −70°C.

For morphological analysis, isolate KNU16-004 was cultured on PDA and incubated at 25°C in the dark for 7 days. Following incubation, colony characteristics such as colour, size, and shape were recorded. In addition, growth was examined on malt extract agar (MEA), synthetic low-nutrient agar (SNA), and oatmeal agar (OA) plates. In order to observe ascomata, ascospores, and asci, the colonies were exposed to illumination with the near-UV light on a 12-hr diurnal cycle. After 7 days of incubation, samples were collected daily for one week and observed for the presence of ascospores using a model BX-50 light microscope (Olympus, Tokyo, Japan).

DNA extraction, PCR, and sequencing analysis

For phylogenic analysis, the genomic DNA of KNU16-004 was extracted using the HiGene Genomic DNA Prep Kit (BIOFACT, Daejeon, Korea) following the manufacturer’s instructions. Primers ITS1F and ITS4 were used to amplify the ITS regions including the 5.8S rRNA region [10]. Primers LROR [11] and LR5 [12] were used for LSU amplification and primers Btub2Fd and Btub4Rd [13] for amplification of the partial Tub2 gene region. The amplified PCR products were purified using ExoSAP-IT (Thermo Fisher Scientific, Waltham, MA, USA) and sequenced (SolGent, Daejeon, Korea). Phylogenetic neighbors were identified and similarities with closely related species were calculated using the BLAST search program on the National Center for Biotechnology Information (NCBI) website (http://blast.ncbi.nlm.nih.gov/Blast.cgi). The resulting sequences of KNU16-004 have been deposited in the NCBI GenBank/EMBL/DDBJ databases with accession numbers LC420014 for LSU rDNA, LC420015 for ITS rDNA, and LC420016 for Tub2.

Phylogenetic analysis

The LSU, ITS, and Tub2 gene sequences were aligned using the Clustal X computer program [14]. Phylogenetic trees were constructed using a concatenated dataset of partial LSU, ITS, and Tub2 gene sequences, assuming that analysis using a longer sequence would result in better resolution and reliability. Nine related taxa were obtained from GenBank and edited using the BioEdit program [15]; these strains are listed in Table 1. Multiple sequence alignments were performed using the Clustal X program. Gaps and the 5′ and 3′ ends of the alignments were manually edited in BioEdit [15]. Evolutionary distance matrices were generated as described by Kimura [16] and phylogenetic trees were constructed using the neighbor-joining [17], maximum-likelihood [18], and maximum-parsimony [19] algorithms in the MEGA7 program [20], with bootstrap values calculated based on 1,000 replications.

Table 1. Isolates used in this study and their GenBank accession numbers

http://dam.zipot.com:8080/sites/kjom/images/N0320460402_image/Table_kjom_46_04_02_T1.jpg

CBS, Westerdijk Fungal Biodiversity Institute (formerly CBSKNAW), Utrecht, The Netherlands; T, ex-type strains;

ITS, internal transcribed spacer regions 1 & 2 including the 5.8S nrDNA gene; LSU, 28S large subunit of the nrRNA gene; TUB, β -tubulin.

The newly generated sequences are indicated in bold.

RESULTS and DISCUSSION

Morphology of isolate KNU16-004

At 25°C, isolate KNU16-004 formed flat colonies with sparse aerial mycelium and smooth and lobate margins. After 7 days of incubation, the growth diameter reached 22 mm on PDA, MEA, and OA agar, but only 12 mm on SNA agar (Fig. 1A~1D). Mycelium initially white and then changed to deep black from the center of the colony. The surface of the colony was covered by short, white aerial mycelium, the reverse side of the colony was initially black, and the end of the colony reached the edge of the PDA plate within 21 days of incubation; ascospores and asci were produced on PDA agar after exposure to the above-mentioned illumination and observed under a light microscope. The ascospores were ellipsoid to oblong with a gelatinous sheath, 30.8~34.9 × 8.1~15.8 μm in size, and the asci were bitunicate and thick-walled, 73.5~87.5 × 23.7~37.9 μm in size (Fig. 1F, 1H). As shown in Table 2, these morphological characteristics of isolate KNU16-004 were completely consistent with those previously reported for L. australis [21, 22].

http://dam.zipot.com:8080/sites/kjom/images/N0320460402_image/Figure_kjom_46_04_02_F1.jpg

Fig. 1. Culture and morphological characteristics of Leptosphaerulina australis KNU16-004. A, Colony on oatmeal agar; B, Colony on potato dextrose agar; C, Colony on malt extract agar; D, Colony on synthetic low-nutrient agar; E, G, Colony sporulating on ascomata; F, Ascospores; H, Asci (scale bar = 10 μm).

Table 2. Morphological characteristics of KNU16-004 isolated in this study and comparison with previously reported isolates

http://dam.zipot.com:8080/sites/kjom/images/N0320460402_image/Table_kjom_46_04_02_T2.jpg

PDA, potato dextrose agar.

aSource of description [21, 22].

Molecular phylogeny of isolate KNU16-004

To identify the isolated fungal strain at the species level, the nucleotide sequences of the LSU, ITS, and Tub2 (890 bp, 519 bp, and 345 bp, respectively) were obtained and compared with the GenBank database using BLAST. The LSU sequence of KNU16-004 exhibited high similarities (99.8~100%) with various L. australis strains, the ITS sequence shared 100% identity with various L. australis strains, and the Tub2 sequence was identical (100%) to L. australis strain CBS 939.69. Phylogenetic analysis based on a concatenated alignment of the LSU, ITS, and Tub2 nucleotide sequences performed using the neighbor-joining, maximum-likelihood maximum-parsimony methods revealed that the strain KNU16-004 belongs to the genus Leptosphaerulina and is grouped with two L. australis strains with a high bootstrap value of 100% (Fig. 2).

http://dam.zipot.com:8080/sites/kjom/images/N0320460402_image/Figure_kjom_46_04_02_F2.jpg

Fig. 2. Neighbor-joining phylogenetic tree based on a concatenated alignment of the large subunit, internal transcribed spacer, and β-tubulin sequences. The phylogenic analysis shows the position of Leptosphaerulina australis KNU16-004 among related Leptosphaerulina spp. strains. Bootstrap values (based on 1,000 replications) greater than 50% are shown at the branch points. Filled circles indicate that the corresponding nodes were also recovered in trees generated with the maximum-likelihood and maximum-parsimony algorithms. Open circles indicate that the corresponding nodes were also recovered in the tree generated with the maximum-likelihood algorithm. The tree was rooted using Herpotrichia juniperi CBS 468.64 as an outgroup. Bar, 0.02 substitutions per nucleotide position. CBS, Westerdijk Fungal Biodiversity Institute (formerly CBSKNAW), Utrecht, The Netherlands.

Pathogenicity status of L. australis is uncertain [23]. However, L. australis has been detected in imported agricultural produce, namely on Brassica oleracea var. capitate leaf spots [21], thus, it could be relevant to agricultural product trade. Further studies will be needed to investigate its potential pathogenicity and identify the relative susceptibility of available cultivars. This is the first report of Leptosphaerulina australis in Korea.

Acknowledgements

This research was supported by a grant from the National Institute of Biological Resources (NIBR), funded by the Ministry of Environment (MOE) of the Republic of Korea for the project on survey and discovery of indigenous fungal species, and Brain Pool Program (Grant No. 2018H1D3A2065415) through the National Research Foundation (NRF) funded by the Ministry of Science and ICT, Republic of Korea.

References

1  1. McAlpine D. Fungus diseases of stone fruits in Australia and their treatment. Melbourne: R. S. Brain, Government Printer; 1902. 

2 2. Graham JH, Luttrell ES. Species of Leptosphaerulina on forage plants. Phytopathology 1961;51:680-93.  

3  3. Irwin JA, Davis RD. Taxonomy of some Leptosphaerulina spp. on legumes in eastern Australia. Aust J Bot 1985;33:233-7.  

4 

5 5. Abler SW. Ecology and taxonomy of Leptosphaerulina spp. associated with turfgrasses in the United States [dissertation]. Blacksburg (VA): Virginia Polytechnic Institute and State University; 2003. 

6  6. Zhang Y, Crous PW, Schoch CL, Hyde KD. Pleosporales. Fungal Divers 2012;53:1-221. 

7 7. Couch HB. Diseases of turfgrasses. Malabar: Krieger Publishing Co.; 1995  

8 8. Wehmeyer LE. The development of the ascocarp in Pseudoplea gaeumannii . Mycologia 1955;47:163-76. 

9  9. Park S, Ten L, Lee SY, Back CG, Lee JJ, Lee HB, Jung HY. New recorded species in three genera of the Sordariomycetes in Korea. Mycobiology 2017;45:64-72. 

10 10. Raja HA, Miller AN, Pearce CJ, Oberlies NH. Fungal identification using molecular tools: a primer for the natural products research community. J Nat Prod 2017;80:756-70.  

11 11. Rehner SA, Samuels GJ. Taxonomy and phylogeny of Gliocladium analysed from nuclear large subunit ribosomal DNA sequences. Mycol Res 1994;98:625-34. 

12  12. Vilgalys R, Hester M. Rapid genetic identification and mapping of enzymatically amplified ribosomal DNA from several Cryptococcus species. J Bacteriol 1990;172:4238-46. 

13 13. Woudenberg JH, Aveskamp MM, de Gruyter J, Spiers AG, Crous PW. Multiple Didymella teleomorphs are linked to the Phoma clematidina morphotype. Persoonia 2009;22:56-62. 

14 14. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res 1997;25:4876-82. 

15  15. Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser (Oxf) 1999;41:95-8. 

16  16. Kimura M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol 1980;16:111-20.  

17  17. Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol 1987;4:406-25.  

18 18. Felsenstein J. Evolutionary trees from DNA sequences: a maximum likelihood approach. J Mol Evol 1981;17:368-76. 

19 19. Fitch WM. Toward defining the course of evolution: minimum change for a specific tree topology. Syst Zool 1971;20:406-16.  

20 20. Kumar S, Stecher G, Tamura K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol 2016;33:1870-4.  

21 21. Hyun IH, Heo NY, Chang SY, Heo JY, Mel’nik V. Identification of three fungi newly intercepted from importing plants in Korea. Mycobiology 2005;33:243-4. 

22  22. Crous PW, Summerell BA, Swart L, Denman S, Taylor JE, Bezuidenhout CM, Palm ME, Marincowitz S, Groenewald JZ. Fungal pathogens of Proteaceae. Persoonia 2011;27:20-45. 

23  23. Mitkowski NA, Browning M.