Abstract
Apple cankers are extremely destructive diseases threatening the global apple industry through direct and indirect losses. The population structure of the pathogens is of paramount significance for the development of efficient management strategies. Therefore, phenotypic, pathogenic, and genetic diversity of Diplodia seriata causing black rot canker of apple was investigated in this study. All the isolates were included for investigating the in vitro mycelial growth, conidial dimensions, and pathogenic variability on two-year-old potted apple seedlings. The ISSR approach was used to investigate the molecular diversity of D. seriata. Mycelial growth rates were found to vary significantly amongst the isolates; however, there were no major variations seen between the different geographical groupings of isolates. Pathogenicity tests revealed variations in the size of cankers among the isolates indicating the presence of virulence variability. The isolates were segregated into three virulence groups based on canker length. The Bayesian analyses of ISSR data divided the isolates into two genetic clusters. The genetic clustering of the isolates revealed no relationship with geographical origin of the isolates. Furthermore, no direct relationship of genetic clustering was observed with morphological or pathogenic variability. The ISSR primers revealed very high level of variability in D. seriata; however, no distinct populations of the pathogen existed which is an indication of high level of gene flow between the diverse geographical populations. According to our knowledge, this is the first thorough investigation on the diversity of D. seriata associated with apple black rot canker in India.
Similar content being viewed by others
Introduction
Apple (Malus domestica) is the oldest temperate fruit which is very popular around the world and assumes an excellent global market. It has varied size, shape, color, taste, texture, aroma, and nutritive value fulfilling the needs and desires of global population. Apple is cultivated in temperate and subtropical areas and ranks fourth globally in production after citrus, grapes, and banana1.
Diplodia seriata has been reported on diverse plant genera due to its cosmopolitan and plurivorous nature. The fungus has been reported from 35 plant species2. Diplodia seriata causes canker, dieback, fruit rot and leaf spots on economically important forest and horticultural species3. A few examples of the diseases and first reports caused by Diplodia seriata are branch canker and Pinus pinea decline4, twig and branch canker in mulberry5, spring bud mortality, leaf chlorosis, fruit rot and trunk dieback on grapevines 6,7, first report of fruit rot of loquat8, first repot of canker and dieback on hawthorn9 and first report of dieback disease on cedars10. Nevertheless, Systematic Mycology and Microbiology Laboratory Fungal Database (ARS, USDA) indicates that the host range for D. seriata now is significantly wider, with over 200 species and 130 genera of herbaceous and woody plants functioning as hosts11. It has been reported to produce different symptoms like cankers, dieback of shoots, fruit rot, and leaf spots on a wide category of horticultural crops12. Diplodia seriata (syn. Botryosphaeria obtusa) is an important pathogen of apples inciting diseases like frog-eye leaf spot, black rot, dieback, and cankers13. The occurrence of black rot on apple fruit was first reported by Peek from New York in 1879, whereas its leaf spot phase was reported by Alwood in 1892 and its destructive canker phase (Black rot canker) was reported by Paddock in 189914. Mundkur and Kheswala15 reported the smoky blight (Sphaeropsis malorum Berk.—synonym of D. seriata) from imported apple plants in Mysore, India. Infection occurs through wounds or by direct penetration through natural openings16,17,18 followed by colonization of the healthy tissues.
The smoky or black rot canker develops on branches, limbs and trunks as small, sunken, reddish-brown lesions which rapidly enlarge lengthwise, become elliptical, develop alternate rings, turn smoky and partially or completely girdle the affected part of the plant. While superficial roughening of the bark is typical, in certain instances the bark is destroyed up to the wood and totally shattered19,20. However, Singh21 recorded the development of cracks along the lesion margins; thus, temporarily restricting canker elongation. In advanced stages the bark becomes loose and gets separated from the wood. The initial symptoms are identified by brown discoloration of the bark on branches and twigs followed by leaf chlorosis. The bark cracks as the disease spreads further into the timber, revealing the reddish-brown discolored wood below. Internal browning of the woody tissues is seen on splitting the infected twigs along the axis. The spurs, twigs and branches above the cankered portion are killed without developing any canker symptoms. On the bark of diseased or dead branches and twigs, many black pycnidia of the fungus may be seen22.
Pathogenicity of D. seriata has been evaluated on shoots, fruits, and leaves of apple in some countries23,24. The experimental data on the pathogenic nature of Diplodia spp. on woody plants is controversial25. The studies on the pathogenicity of different D. seriata strains on various hosts and in different regions have shown inconsistent results26. In England and New Zealand, it is considered as secondary and weak pathogen of apples27, while reported as highly pathogenic on apples in the USA14,28,29. Similarly, contrary reports on the pathogenic nature of D. seriata on grapevines have been reported from different parts of the world12. Although different symptoms like dark vascular discoloration, extended necrotic lesions on shoots and formation of cankers have been linked to pathogenicity of Botryosphaeriaceae species30, pathogenicity and virulence of the canker causing fungi is normally estimated by measuring their capability to develop brown necrotic lesions beneath the bark and the lesion size, respectively12. Larignon31 studied the pathogenic variability in D. seriata and divided the isolates into four virulence groups based on variations in lesion size. The variations in pathogenicity level of D. seriata isolates on grapevines has been reported in earlier studies6,30,32. Savocchia33 also recorded pathogenic variability in D. seriata isolates based on intraspecific variations within this fungus.
For genetic diversity analysis, inter-simple sequence repeat (ISSR) or random amplified microsatellite (RAMS) technique34,53 have been proved effective35. The relative participation of telomorphic and anamorphic forms contributes to genetic diversity of a species. For Botryosphaeriaceae species, the telomorphs are rarely found in nature6 and the reproduction is mainly asexual under natural conditions30,36. However, exchange of genetic material can take place between compatible genotypes by parasexual recombination37. In Spain the genetic diversity of D. seriata isolates derived from grapevine was studied by Martín38 using AFLP markers but no data on pathogenic diversity was provided. Elena12 studied phenotypic, pathogenic, and genetic diversity of D. seriata on grapevine and reported high level of diversity within populations but no relationships among the genetic grouping and the virulence of different isolates was observed. Investigations on genetic variance between the populations of Lasiodiplodia theobromae revealed very low gene flow and genotypic diversity between the geographic populations, and increased gene flow and genotypic diversity among the pathogen populations resident on different host plants39. Genetic and virulence variability in the isolates of Neofusicoccum parvum was reported by Baskarathevan40. Nevertheless, no association was found between genetic clustering and isolate virulence. Genetic diversity existed between 66 isolates of Neofusicoccum luteum belonging to different pathogenic groups; however, the pathotype–genotype relationships were not found41. The genetic diversity studies on Togninia minima reported no significant differences in virulence between the genetic groups42. To the best of our knowledge, no comprehensive studies on the combined genetic, phenotypic, and pathogenic characterization on the isolates of D. seriata associated with black rot canker of apple have been undertaken.
The main purpose of the present investigation was to acquire more satisfactory comprehension of the phenotypic, pathogenic, and genetic diversity of D. seriata responsible for black rot canker of apple trees in India. The pathogenicity assays conducted on young potted apple seedlings were used to estimate the pathogenic variability, whereas mycelial growth rates and conidial dimensions were used for phenotypic characterization. The final aim of our investigations was to verify the potential existence of distinct populations of D. seriata in different geographical regions based on combined morphological, pathogenic, and genetic studies, and the expansion of our understanding regarding D. seriata which may help us in understanding the pathogenic nature of D. seriata in apple and will have an impact on the efficient disease management strategies.
Results
Morphological variability
Mycelial growth
The radial mycelial growth of isolates of D. seriata was assessed by growing the isolates on PDA. All the isolates were involved in this investigation (Fig. 1). The isolates showed significant differences in mycelial growth rates (P < 0·001). The mean mycelial growth rate for all the isolates with 95% confidence interval was 2.66 ± 0.87 cm, with minimum and maximum mycelial growth rates as 2.13 cm and 3.18 cm, respectively. The variation in radial mycelial growth of the isolates had no relation with the grouping of isolates based on geographical origin as well as with virulence groups. The statistically significant differences in growth rates were not observed among the geographical groups of isolates (Table 1).
Conidial characteristics
All the isolates of D. seriata under experimental conditions produced conidia on oatmeal agar medium. The mean conidial size of the isolates grouped based on geographical origin was not found to present considerably different values in length, width, and length:width ratio between the groups. Among the geographical groups, the absolute differences calculated between the group means were 0.87 µm between group 1 and group 2 (group 1, 22.91 ± 0.61 µm; group 2, 22.04 ± 1.22 µm), 0.23 µm between group 1 and group 3 (group 3, 23.14 ± 0.44 µm), 0.81 µm between group 1 and group 4 (group 4, 23.72 ± 0.65 µm), 1.10 µm between group 2 and 3, 1.68 µm between group 2 and 4, 0.58 µm between group 3 and 4 in length. Similarly the absolute differences in means of width among the geographical groups were as; 0.14 µm for group 1 and 2 (group 1, 10.58 ± 0.31 µm; group 2, 10.72 ± 0.44 µm), 0.25 µm between group 1 and 3 (group 3, 10.33 ± 0.36 µm), 0.50 µm between group 1 and 4 (group 4, 10.08 ± 0.15 µm), 0.39 µm between group 2 and 3, 0.64 µm between group 2 and 4, 0.25 µm between group 3 and 4, respectively. The mean L:W ratio of the geographical groups was 2.18 ± 0.08 for group 1, 2.07 ± 0.12 for group 2, 2.26 ± 0.07 for group 3 and 2.36 ± 0.07 for group 4 respectively (Table 2). The mean length: width ratio for all the isolates was 2.22 ± 0.04. For all the isolates, the min., maximum and mean ± 95% confidence interval values were in the following range: length (14.13) − 22.95 ± 0.77 μm − (29·24) × width (8.91) − 10.43 ± 0.32 μm − (13·75).
Pathogenic variability
All the isolates of D. seriata were tested for pathogenicity on potted apple seedlings. Each of the isolates produced smoky canker symptoms characterized by water soaked, yellow to brown discoloration at the site of incision. The lesion showed upward and downward extension, girdling and rupturing of bark. Characteristic concentric ring development was also observed in the lesions (Fig. 2). During the investigation, significant differences in canker length between the isolates were observed in statistical analysis. All the isolates produced cankers considerably larger than control (Fig. 3) with mean lesion length varying from 7.08 cm (isolate 12) to 17.43 cm (isolate 41) and an overall mean lesion length of 10.35 cm. In each geographical group of isolates the mean lesion length was considerably longer than that of control but no considerable differences were observed amongst the groups (Table 1). The mycelial growth rates recorded under in vitro conditions and the length of canker lesions recorded during pathogenicity tests were not found to have any relationships. The overall average positive re-isolation from the inoculated twigs was 2.94 out of 3 twigs, i.e., 98% re-isolation frequency; thus, proving that all the isolates were pathogenic in nature.
The isolates exhibited a considerable variability in incubation periods during pathogenicity tests. Based on incubation period the isolate were divided into three groups. The first group comprised of isolates for which the mean incubation period was 31.93 ± 2.60 days to develop canker symptoms on the inoculated plants. A total of 37.50% isolates belonged to this group. The second group comprised of 35.41% isolates which exhibited an incubation period of 44.59 ± 2.61 days, while as the third group comprising of 27.08% isolates was found to show the longest incubation period of 59.79 ± 2.49 days (Table 3).
The isolates also varied in their virulence demonstration during pathogenicity trials. Based on average lesion size, the isolates were divided into three virulence groups viz. Group I-least virulent which produced a canker lesion of < 100 mm length; Group II-moderately virulent which produced a lesion size of 100–150 mm length and Group III-highly virulent which produced a canker lesion of > 150 mm length (Table 4). Cankers show more extension longitudinally than radial expansion, for that reason canker length was used as a standard for comparing canker sizes. Group I accommodated 47.91 per cent of isolates (23 isolates) and was categorized as least virulent group. Group II accommodated 41.66 per cent isolates (20 isolates) and was considered as moderately virulent. Group III accommodated only 10.41 per cent of isolates (5 isolates) and was considered as highly virulent group (Table 4).
Genetic diversity
Initially 17 ISSR primers were screened for standardization using 10 randomly selected isolates of the pathogen. Among these 12 ISSR primers failed to produce polymorphic, reproducible, and scorable bands and hence did not yield any information. The remaining 5 primers resulted in polymorphic and reproducible bands, hence, were adopted for studying genetic diversity between the isolates of D. seriata. DNA fingerprinting of D. seriata isolates derived from diverse geographical locations clearly showed distinct banding pattern. The amplification products produced 65 bands in total ranging from 300 to 2000 bp, out of which 63 (96.92%) were polymorphic (Table 5).
Under similar experimental conditions, all the primers amplified the DNA from ten randomly chosen isolates while testing for reproducibility of the results. Out of 50 total combinations (10 isolates × 5 primers) band patterns identical to band arrangements in ISSR characterization were detected in 21 and the proportion of coincident bands exceeding 90% were observed in 35 combinations. Among all the combinations the coincident bands were estimated to be present with an overall percentage of 94%.
Genetic variation analysis
Four distinct clustering techniques were used to determine the most probable number of communities that the samples obtained would represent, assign individuals to these populations, and estimate the degrees of variance. Except for one isolate, i.e. DS38, which split as an independent lineage, the cluster analysis based on the unweighted Neighbor-joining approach yielded a dendrogram that grouped all the isolates representing four geographical populations into two primary clusters. The cluster first comprised of 26 isolates while as cluster second comprised of 21 isolates. The assignment of isolates into these two clusters was independent of their geographical background as both the clusters contained individuals from all the four geographical populations. Both the major clusters were further divided into smaller sub-clusters which were also independent of the geographical origin of the isolates. This demonstrates that genetic diversity exists both inside and among population isolates (Fig. 4).
The average number of observed alleles (Na) and effective alleles (Ne) across all the loci varied between 1.385 and 1.662 and 1.345 and 1.393, respectively, across the four distinct regional populations (PL, SH, BR, and GN). The highest number of observed alleles (1.662) were found in the population SH belonging to district Shopian of the valley, while as lowest number of observed alleles (1.345) were found in the population GN belonging to district Ganderbal, with mean value of 1.504. Again, the highest effective allele number (1.393) was found in the population SH while as the lowest number of effective alleles (1.345) were found in the population PL belonging to district Pulwama, with an average of 1.372. Similarly, the values of gene diversity (h) and unbiased gene diversity (uh) ranged from 0.218 to 0.242 and 0.238 to 0.264, respectively. The highest values of h (0.242) were recorded for the population SH pertaining to geographic locations of district Shopian and lowest values of h (0.218) were recorded for the population PL pertaining to geographic locations of district Pulwama, with an average of 0.228. Similarly, the highest values of uh (0.264) were recorded for the population SH belonging to geographic locations of district Shopian and lowest values (0.238) were recorded for the population PL, with an average of 0.248. Based on the Shannon’s information index (I), the highest diversity was recorded in the population SH (I = 0.377) and lowest diversity was recorded in the population PL and GN (I = 0.341), with an average of 0.350. The percentage of polymorphic loci among the populations ranged from 66.15% in GN to 81.54% in SH with an average of 72.31% (Table 6).
The average genetic distance and identity of Nei among D. seriata populations was 0.018–0.023 and 0.977–0.982, respectively (Table 7). The BR population of Baramullah, the SH population of Shopian, and the GN population of Ganderbal were found to have the highest genetic differences. The SH populations of Shopian and the PL populations of Pulwama were found to have the least genetic distance. Likewise, the SH population of Shopian and the PL population of Pulwama demonstrated the highest levels of genetic similarity. Ganderbal's GN and Baramulla's BR populations had the fewest genetic characteristics. The findings of Nei's genetic distance analysis showed that groups that were geographically close to one another had the lowest genetic distances, whereas populations that originated in various locations had the greatest genetic distances.
The geographically defined populations had considerable genetic diversity, according to the analysis of molecular variance (AMOVA). Although just 18% of the overall genetic diversity was found across populations from various geographical locations, most of the genetic variation (82%) was found within populations (among individual isolates) (Table 8; Fig. 5).
It was noted that there were statistically significant differences between the populations (P = 0.01). Pairwise PhiPT values estimated the degree of genetic divergence between the four populations. The average degree of genetic similarity was calculated to be 0.178 (average PhiPT value) (Table 9). The populations PL and GN showed the largest genetic difference (0.230), whereas the populations GN and BR showed the lowest genetic differentiation (0.076) (Table 9). The pairwise PhiPT values between the populations of South Kashmir (SH and PL; PhiPT = 0.103) were much lower than those between the populations of South and North Kashmir (PL-BR, SH-GN, and SH-BR), indicating that there was less genetic difference between these groups. The Mantel test confirmed any potential association between the geographic and genetic distances between the corresponding groups. The test's results revealed no connection between genetic and geographic distances (r = 0.591; p = 0.210) (Fig. 6).
The analysis performed with the help of structure software divided the 48 D. seriata isolates into two genetically distinct clusters. The log likelihood values (Ln P(D)) increased from K = 1 to K = 2 and was maximum at K = 2 (ln P(D) = 2145·68). The standard deviation showed a slow increase from K = 1 to K = 2 and thereafter it showed a rapid increase from K = 3 to K = 10. This showed that the population structure of the D. seriata pathogen population is best explained by K = 2 clusters. This is further strengthened by the value of ∆K which was highest at K = 2 (∆K = 101·39) and after that it decreased at a very fast rate. The highest value of ∆K that determines the best fit of K was derived by employing STRUCTURE HARVESTER (Fig. 7).
The clusters formed in the structure analysis are shown in the figure below. Most of the individual isolates were find out to possess 80% of their heredity from one of the two distinct populations (Fig. 8), while as 17 isolates were found to have mixed ancestry with less than 80% posterior probability to any of the two clusters. The cluster in red consisted of 14 D. seriata isolates while as the cluster in green consisted of 17 isolates. The distribution of isolates into different clusters was irrespective of their geographical origin, hence both the clusters contained isolates from all the four predetermined geographical populations.
The principal coordinate analysis (PCoA) was carried out with Genalex software on the whole data set to verify the consistency of genetic differentiation clustering performed by different clustering methods. The two-dimensional PCoA plot revealed that 15.45% and 9.20% of the total genetic variation is represented by the coordinate 1 and coordinate 2, respectively. It is evident from the PCoA plot that except two individuals as outliers and two small groups were formed, and all the other individuals formed a common group, however the individuals in all the groups were intermixed (Fig. 9).
This shows that the PCoA clustering was independent of the geographical origin of the isolates which is in accordance with other clustering analysis methodologies. The heatmap constructed with R software also divided the population of all the isolates into two major clusters which are further subdivided into sub-clusters. It is clear from the heatmap that the genetic variation is spread among the major clusters and sub-clusters irrespective of the geographical origin of the individual isolates (Fig. 10).
Discussion
In the present investigation, the phenotypic traits of the fungal isolates collected from apple black rot canker were in complete accord with those of the genus Dilpodia43. Depending on the color and septation of their conidia, the species of Diplodia may be roughly classified into two distinct groups based on morphological descriptions44. Most Diplodia species, including D. malorum, D. mutila, D. corticola, D. africana, and D. rosulata, have aseptate, thick-walled conidia that may stay hyaline for a long time before becoming brown43,45,46,47,48. Nevertheless, in other species, including D. alatafructa, D. pinea, D. scrobiculata, D. seriata, and D. pseudoseriata, the conidia become colored before discharge from the pycnidia and mostly stay non-septate43,44. All of the isolates morphological characteristics in this investigation were consistent with those previously described for D. seriata1,13. Additionally, no significant differences were identified across regional or genetic groups.
Each of the isolates produced characteristic black rot canker symptoms on potted apple trees and the isolates were recovered from 98% of the inoculated apple twigs, hence proving the isolates to be potential pathogens of apple. Botryosphaeriaceae fungi being cosmopolitan in nature13 have been reported to be common pathogens of a vast diversity of woody hosts including apple49,50,51. Pathogenicity of D. seriata has been evaluated on shoots, fruits, and leaves of apple in some countries23,24. The experimental data on the pathogenic nature of Diplodia spp. on woody plants is contentious. For instance, D. seriata is regarded as a secondary and weak pathogen of apples in England and New Zealand27, while it is reported as highly pathogenic on USA apples where it is responsible for apple cankers, leaf spots and fruit rot14,28,29. Similarly, D. seriata, being highly virulent on grapevines has been considered as a primary and virulent pathogen31,33 but other workers have observed it to be saprophyte or weak pathogen on grapevines6,36,43. These contradictions in pathogenicity may be due to differences in virulence among the isolates. However, in the present investigation all the isolates of D. seriata were found to be highly pathogenic to apple as reflected by the length of lesions produced on the inoculated apple seedlings which were much larger than the lesions in control.
A significant variation was observed in isolates of D. seriata with respect to lesion length on artificially inoculated apple trees of cultivar Red Delicious. The variation in the length of canker lesions caused by different isolates under pathogenicity tests confirms the existence of virulence variability among the isolates of D. seriata. The isolates were categorized into 3 virulence groups depending on average canker length. The isolates showed no relevance between the virulence and geographical origin of isolates which was later supported by genetic grouping as well. Based on length of lesions produced in pathogenicity tests on grapevine canes, the D. seriata isolates were divided into four virulence groups31. Variation in D. seriata isolates with reference to wood colonization on grapevines has also been reported6. Similar findings on the variability in virulence of D. seriata have been recorded by other workers as well30,32. While discussing the contrasting reports regarding the pathogenicity of D. seriata, Savocchia et al.33 concluded that these differences may be because of the presence of intraspecific variations within this fungus. Our findings are well supported by the observations of Elena et al.12 who reported “no relationship between the genetic groups and the virulence of isolates of D. seriata on grapevine”. Similar conclusions were drawn by other workers e.g., Billones-Baaijens et al.41 and Baskarathevan et al.40 who observed no relationship between genotypic and pathogenic variability in the isolates Neofusicoccum luteum and N. parvum, respectively. Again, no correlation was found to exist between the radial mycelial growth and the virulence variability of isolates because the isolates which showed maximum radial mycelial growth did not always produce longest cankers as it should have been. This could be attributed to the reasons like production of phytotoxic secondary metabolites by D. seriata as reported by Martos et al.52, or presence of some other mechanism besides fungus mycelium growth inside the host which may be influencing the virulence of the pathogen12.
The ISSR technique proved effective for ascertaining the genetic diversity between the isolates of D. seriata having diverse geographical origins. Five ISSR primers recorded a high degree of polymorphism in the isolates of D. seriata, hence proving that ISSR markers are very relevant in verifying the genetic diversity in this fungus. The ISSR markers have been found relevant for investigating the genetic variability in D. seriata and other fungal plant pathogens by various other workers12,54,56,57,58,59.
For investigating the population structure of D. seriata isolates Bayesian clustering of ISSR markers proved very useful and provided the basis to investigate the genetic structure. The analyses performed with STRUCTURE divided the isolates into two main clusters however; it was observed that some of the isolates did not constantly belong to these clusters. Further, the isolates originating from diverse geographical regions did not fall in different genetic groups, indicating that the isolates belong to a common ancestral background irrespective of their geographical origin. Elena et al.12, reported high genetic diversity in D. seriata on grapevine from Spain using ISSR markers and recorded that there was no relation between the clustering groups and the geographic origin of isolates. While performing ITS sequencing of D. seriata isolates Phillips et al.13 concluded that no correlation existed between the clustering structure and the geographical origin of the isolates.
Clustering of D. seriata isolates at regional level revealed that most of the variability was found within groups than between the groups. This indicates that the genetic diversity between the geographic populations is less but within populations diversity is high. Padder et al.56 reported high genetic diversity in Venturia inaequalis populations on apple in Kashmir valley using ISSR markers and found greater variability within populations than among populations. Dinolfo et al.57 studied intraspecific diversity within the isolates Fusarium poae derived from Argentina and England but there was no definite relation between genetic diversity and geographic origin, and diversity was high within rather than between the Argentinean and English populations. Fan et al.55 analyzed genetic variation in Chinese and Californian isolates of Monilinia fructicola using ISSR markers and found that the diversity was high within population rather than between populations. Many other researchers have also reported high genetic diversity in fungal plant pathogens using ISSR markers58. Although only two clusters were derived in the STRUCTURE analysis of the data generated from the 48 isolates of D. seriata, the ISSR technique revealed a high level of polymorphism (96%). Similar results were recorded by Elena et al.12 who reported 88% polymorphism in 83 D. seriata isolates from Spain. Employing AFLP markers Martin38 recorded 86% polymorphism in the isolates of D. seriata derived from Spain. With UP-PCR technique Billones-Baaijens et al.41 recorded 93% polymorphism in the isolates of N. luteum. Using the same technique, Baskarathevan et al.40 also recorded 93% polymorphism in 60 isolates of N. parvum. Many factors such as gene flow, recombination (sexual or parasexual), mutation and adaptation of the pathogen to diverse host species may be responsible for the presence of high genetic diversity in D. seriata. Although Hesler recorded the formation of perfect state (Botryosphaira obtusa) for the smoky canker pathogen on apple twigs but other researchers from India and other countries failed to report any sexual state in field22. Thus, further investigations are required to be undertaken both under natural field conditions and at genetic levels to verify whether some hidden telomorphic state is responsible for the high genetic diversity in this fungal species. Estimates of pairwise FST were low in geographical populations. This indicates that there is frequent gene flow between populations probably due to movement of infected planting material.
Conclusions
The morphological differences shown by the isolates in phenotypic characterization did not show any correlation with the genetic clusters. Moreover, the pathogenicity tests revealed D. seriata as highly pathogenic on apple and indicated variability in virulence, but the pathogenic variability was not correlated to either the genetic variability or the geographical origin of the isolates. To the best of our knowledge, this combined phenotypic, pathogenic, and genetic variability of D. seriata isolates causing black rot canker, dieback, and fruit rot of apple in India is studied for the first time and it will help in understanding the pathogenic nature of the D. seriata populations and development of more effective management strategies in future.
Materials and methods
Fungal isolates
The disease samples collected from different parts of the apple trees during the years 2016 and 2017 were cleaned of by using distilled water. The bark at the diseased portions was carefully removed and the area was surface sterilized with 70% ethanol. Wood sections (5 mm2) were taken along the margins of canker lesions by using sterilized blade and the tissue bits were sterilized for 60 s in 1% sodium hypochlorite solution. The bits were washed thrice with sterilized distilled water to eliminate the remains of chemical followed by drying on blotting paper and inoculated on acidified PDA. The inoculated plates were placed under incubation at 25 ± 1 °C in the darkness. All the isolates of D. seriata collected from different geographical areas were purified by single spore technique60. D. seriata was identified based on morphological characters as per the descriptions of fungus given by Phillips13. In all, 48 isolates (Table 10) of D. seriata having different geographical origins were maintained for further studies and the stock cultures were preserved under refrigeration at 4 °C inside the refrigerator.
Pathogenicity tests
Pathogenic variability was studied by performing pathogenicity tests of all the 48 isolates on healthy two-year-old potted apple seedlings (cv. Red Delicious grafted on Quince root stock) following the procedure adopted by Spiers61. The selected plants were first treated with copper oxychloride 50 WP (Blitox 50) @ 0.3% as a spray to eliminate the existence of any pathogens which may be present on the host surface. The pots were kept inside a polyhouse and timely irrigated to prevent the young apple seedlings from water stress and intermittently sprayed with distilled sterilized water to maintain high humidity inside the chambers. The inoculations were done on dormant apple twigs. Each isolate was inoculated on three twigs and the tests were performed as per the following procedure.
On the selected twigs, after surface sterilizing the designated area with 70% ethanol, a 5 mm horizontal cut was given cutting the bark of the twig. Starting from the mid-point of the first cut, another vertical cut of 5 mm was made. The bark of this T-shaped cut was loosened with the help of a sterilized scalpel. A mycelial disc (4 mm) was put beneath the bark of incision and a cotton swab moistened with sterilized water was placed on the inoculated cut followed by wrapping with adhesive paper tape. The control was inoculated with plain media (PDA) disc and the incisions were covered as described earlier. The pathogenicity tests were regularly tracked for canker formation. Observations in terms of symptom development and lesion size by different isolates were recorded. The pathogen was re-isolated from the inoculated twigs to fulfill Koch's postulates. For re-isolations 1 cm pieces of wood after removing bark were cut at the two ends of canker lesions and sterilized in ethanol (70%) for 4 min, blotted dry and inoculated on Petri plates containing acidified PDA. The plates were kept under incubation at 25 ± 1 °C and observed for morphological characters of D. seriata.
Phenotypic characterization
Radial mycelial growth
All the 48 isolates of D. seriata from the four geographical locations were included in this study. Each Petri plate containing PDA was inoculated in the center with 5 mm culture bit of the actively growing isolate of D. seriata followed by incubation of the plates at 25 ± 1 °C. Each isolate was grown on four PDA plates. The observations on radial mycelial growth were recorded after 72 h by taking average of four radial measurements perpendicular to each other. The test was conducted thrice for the purpose of accuracy.
Conidial size and shape
Conidial characteristics such as shape, size, color and septation of all the 48 isolates were assessed on oatmeal agar medium. The cultures were kept on a bench in the laboratory for 4 weeks at about 20–25 °C for inducing pycnidia formation under diffused sunlight. The spore measurements of each isolate were recorded under light microscope [HF (bright field, filter 1); Carl Zeiss] at 100× magnification by preparing semi-permanent slides in lactophenol under 10 random microscopic fields after calibrating the microscope following standard procedures64. For each isolate the conidial dimensions (length and width) of 50 conidia were recorded, and the length: width ratio was estimated.
Genomic DNA extraction and PCR amplifications
For raising enough mycelium for DNA extraction, all the isolates of D. seriata were grown on PDA broth in Erlenmeyer flasks. Mycelium was harvested by filtering the contents of the flask through a double layered filter paper, dried and stored at − 80 °C. Adopting modified CTAB (Cetyl trimethyl ammonium bromide) method62, total genomic DNA was extracted from all the isolates. The DNA was kept at 4° C for amplification and stored at − 20° C for future use.
A set of 17 ISSR primers namely (CAG)5, (ACA)5, BDB(ACA)5, (AG)5, (AT)8B, (CA)5, (GA)5YC, (TC)8RA, (CCA)5, (CGA)5, HBH(AG)7, (CT)5, (TG)5, (AG)5YC, (AG)5YA, (TG)8RC and BBD(AAC)5 were chosen for conducting preliminary screening against 10 randomly selected isolates to study the polymorphism. Five primers which were consistent in polymorphism during preliminary screening were selected for ISSR characterization of all the 48 isolates of D. seriata. PCR amplifications were performed in 0.2 ml PCR tubes in a T-Gradient Whatman Biometra thermal cycler with 50–60 ng of genomic DNA per isolate. The final volume of PCR mix was adjusted to 25 μl per reaction with RNAse-free water.
The reaction mixture in PCR tubes were given short spin or vortexed in microfuge and placed in 96-well thermal cycler. PCR amplifications were conducted in thermal cycler programmed for initial step at 94 °C for five minutes followed by 35 denaturation cycles at 94 °C for one minute, annealing at appropriate temperature (Table 11) for one minute, elongation at 72 °C for two minutes and a final step of extension for 10 min. To 25 μl of the amplified PCR product, 5 μl of 6X loading dye was added to have final concentration of 1X.
The quality of PCR products was evaluated on 2% agarose gel. TAE (0.5X) was used to prepare the gel and ethidium bromide at 0.5 μg/μl was added to it. The gel was loaded with amplified products of all the 48 isolates of D. seriata and run at 5 V/cm and visualized under UV light. The photograph of resultant banding pattern was taken using Alfa Imager gel documentation system. The size of ISSR alleles was determined by their relative positions with respect to DNA ladder (100 base pair). The amplification reactions of all the primers were repeated thrice for a single DNA extraction of each of the isolates. Further, for band scoring only clear and reproducible bands were taken into consideration.
ISSR characterization and analysis of genetic variability
The banding pattern of each primer was scored as 1 (band present) and 0 (band absent) and the pattern was converted into a binary matrix with the help of MS Excel. The binary data matrix was used to work out TNB (total number of bands), NPB (number of polymorphic bands), and PPB (percentage of polymorphic bands). The Jaccard’s coefficient of similarity was calculated in the PAST (v 4.0) software by making pairwise comparisons between the individual isolates. Furthermore, a dendrogram was constructed with the help of same software by adopting the UPGMA method. The cophenetic correlation coefficient was worked out by the PAST program for measuring the goodness of fit between the cophenetic matrix and the distance matrix69. For each marker the alleles were calculated in all the isolates by numbering the alleles as 1, 2, 3 and so on, starting from the allele possessing maximum molecular weight. The bands were visually scored as 0 (absent) and 1 (present). The binary data generated by ISSR primers was analyzed by cluster analysis utilizing the program NTSYS-pc version 2.0263. Based on shared bands among the isolates, similarity was estimated from the data generated by all the 5 ISSR primers using Jaccard’s coefficient. The NTSYS software was utilized for determining the Nei and Li65 coefficients to ascertain the similarity between the isolates of D. seriata. The investigation of genetic diversity among the isolates of D. seriata was done by analysis of the band arrangements using a Bayesian clustering approach executed in STRUCTURE v.2.2.367 considering a population admixture model wherein within population correlations of allele frequencies are worked out. The runs for each K value (1 to 10) were set for 20 independent replications having a burn‐in period of 100,000 accompanied by 500,000 iterations of the Markov Chain Monte Carlo (MCMC). The posterior probability of the data (ln P(D)) was estimated utilizing the ad hoc ΔK statistic, based on rate of change between successive K values, and this estimate was later employed for inferring the uppermost level of structure in the dataset67. The highest values of Ln P(D), which represents the posterior probability of the data for a certain K, and ΔK were used to determine K (the rate of change in the log probability of data between successive values of K)76. The sample populations were divided into genetic groupings using bar graphs that represent individual assignment probability. By using Structure Harvester v6.0 to generate the K values, the best-fit population cluster number (K) was verified76. Employing ARLEQUIN v. 3.568, AMOVA was used to find out the differences in assigning the isolates to different clusters which were obtained after Bayesian analyses.
Further, the pairwise population differentiation (PhiPT) was also worked out by utilizing Genalex 6.503 package69,70,71. Microsatellite data set were analyzed for basic genetic diversity indices such as, N: Sample size; Ppl: percentage of polymorphic loci; Na: observed number of alleles; Ne: effective number of alleles; I: Shannon’s information index; h: diversity; uh: Unbiased diversity using population genetic software GenALex 6.50371.
The principal coordinate analysis (PCoA) was conducted by employing Genalex software to generate a coordinate graph of the most significant axes to verify the clustering pattern. Distribution of the genetic variation into two components viz. within populations and among the populations was calculated by performing the analysis of molecular variance (AMOVA) in Genalex at 999 permutations.
Statistical analysis
To estimate the population genetic structure and examine the distribution of genetic diversity, the isolates of D. seriata were primarily grouped into populations based on geographical origin. The pathogenicity experiments were conducted following randomized complete block design under poly house conditions. Groups based on geographic origin and those generated by the STRUCTURE analysis were considered as independent variables. ANOVA (One‐way) and posterior Tukey's Honestly Significant Difference tests were performed on dependent variables to find out the differences between the mean values. The molecular data generated during the study was analyzed with NTSYS- 2.02e; Arlequin v. 3.5; GenAlex 6.503; PAST 4.0 and STRUCTURE v. 2.2.3 as described earlier. The remaining data was analyzed by SPSS (ver. 18) and SAS (ver. 9.2) software.
Data availability
All data are within the manuscript.
References
O’Rourke, D. World production, trade, consumption and economic outlook for apples. in Apples: Botany, Production and Uses 15–29 (CABI, 2003). https://doi.org/10.1079/9780851995922.0015.
Punithalingam, E. & Waller, J. M. Botryosphaeria obtusa. Descriptions of pathogenic fungi and bacteria, No. 394. Commonw. Mycol. Inst., Kew, Surrey, England (1973).
Farr, D. F., Rossman, A. Y. 2020 – Fungal Databases, U.S. National Fungus Collections, ARS, USDA. (2020) https://nt.ars-grin.gov/fungaldatabases/.
Hlaiem, S. et al. – Characterization and pathogenicity of Diplodia seriata causing branch canker on Pinus pinea in Tunisia. Plant Pathol. Quar. 12(1), 124–132. https://doi.org/10.5943/ppq/12/1/10) (2022).
Arzanlou, M. & Dokhanchi, H. Morphological and molecular characterization of Diplodia seriata, the causal agent of canker and twig dieback disease on mulberry in Iran. Arch. Phytopathol. Plant Protect. 46(6), 682–694. https://doi.org/10.1080/03235408.2012.749703 (2013).
van Niekerk, J. M., Crous, P. W., Groenewald, J. Z. & Ewald, P. H. Fourie & Halleen, F. DNA phylogeny, morphology and pathogenicity of Botryosphaeria species on grapevines. Mycologia 96, 781–798. https://doi.org/10.2307/3762112 (2004).
Urbez-Torres, J. The status of Botryosphaeriaceae species infecting grapevines. Phytopathol. Mediterr. 50(4), 5–45. https://doi.org/10.14601/Phytopathol_Mediterr-9316 (2011).
Abbas, M. F. & Naz, F. First report of Diplodia seriata causing fruit rot of loquat in Pakistan. J. Plant Pathol. 100, 325. https://doi.org/10.1007/s42161-018-0048-4 (2018).
Kurt, Ş et al. First report of Diplodia seriata causing canker and dieback on hawthorn trees in Turkey. J. Plant Pathol. 104, 1543. https://doi.org/10.1007/s42161-022-01169-0 (2022).
Jiao, Z. J., Kan, Y. C. & Huang, S. L. First report of dieback disease on cedars caused by Diplodia seriata in China. Plant Dis. 98(9), 1279. https://doi.org/10.1094/PDIS-03-14-0309-PDN (2014).
Farr, D. F. Fungal databases (2010). http://nt.ars-grin.gov/fungaldatabases/
Elena, G., Garcia-Figueres, F., Reigada, S. & Luque, J. Intraspecific variation in Diplodia seriata isolates occurring on grapevines in Spain. Plant Pathol. 64, 680–689. https://doi.org/10.1111/ppa.12296 (2015).
Phillips, A. J. L., Crous, P. W. & Alves, A. Diplodia seriata, the anamorph of Botryosphaeria obtusa. Fungal Divers. 25, 141–155 (2007).
Brown, E. A. & Britton, K. O. Botryosphaeria diseases of apple and peach in the southeastern United States. Plant Dis. 70, 480–484. https://doi.org/10.1094/PD-70-480 (1986).
Mundkur, B. B. & Kheswala, K. F. A canker of apple trees in Mysore. Indian J. Agric. Sci. 13, 397–398 (1943).
Brown, E. A. & Hendrix, F. F. Pathogenicity and histopathology of Botryosphaeria dothidea on apple stems. Phytopathology 71, 375–379. https://doi.org/10.1094/Phyto-71-375 (1981).
Britton, K. O. Infection of Peach Buds by Botryosphaeria obtusa. Plant Dis. 73, 65. https://doi.org/10.1094/PD-73-0065 (1989).
Kim, Y. K., Xiao, C. L. & Rogers, J. D. Influence of culture media and environmental factors on mycelial growth and pycnidial production of Sphaeropsis pyriputrescens. Mycologia 97, 25–32. https://doi.org/10.1080/15572536.2006.11832835 (2005).
Sutton, T. B. Production and dispersal of ascospores and conidia by Physalospora obtusa and Botryosphaeria dothidea in apple orchards. Phytopathology 71, 584. https://doi.org/10.1094/Phyto-71-584 (1981).
Ellis, M. Black Rot and Frog-Eye Leaf Spot of Apple. (2008). https://ohioline.osu.edu/factsheet/plpath-fru-07
Singh, R. S. Diseases of Fruit Crops (Science Publishers, Inc., 2000). https://www.cabdirect.org/cabdirect/abstract/20013037056
Khan, N. A. Status and Etiology of Canker Diseases of Apple in Kashmir. Ph. D. Thesis, (Division of Plant Pathology, SKUAST-Kashmir, India, 2010) 118p.
Foster, H. H. Studies of the pathogenicity of Physalospora obtusa. Phytopathology 27, 803–823 (1937).
Borovinova, M. Susceptibility of apple cultivars to Botryosphaeria obtusa. Not. Bot. Horti Agrobot. Cluj-Napoca 34, 100–105. https://doi.org/10.15835/nbha341298 (2006).
Phillips, A. J. L., Lopes, J., Abdollahzadeh, J., Bobev, S. & Alves, A. Resolving the Diplodia complex on apple and other Rosaceae hosts. Persoonia - Mol. Phylogeny Evol. Fungi 29, 29–38. https://doi.org/10.3767/003158512X658899 (2012).
Abdollahzadeh, J. Diplodia bulgarica, as a new pathogen and potential threat to the apple industry in Iran. Phytopathol. Mediterr. 54, 128–132. https://doi.org/10.14601/Phytopathol_Mediterr-14686 (2015).
Laundon, G. F. Botryosphaeria obtusa, B. Stevensii and Olthia spiraceae in New Zealand. Trans. Br. Mycol. Soc. 61, 369–374. https://doi.org/10.1016/S0007-1536(73)80158-5 (1973).
Stevens, N. E. Two apple black rot fungi in the United States. Mycologia 25, 536. https://doi.org/10.2307/3753998 (1933).
Brown-Rytlewski, D. E. & McManus, P. S. Virulence of Botryosphaeria dothidea and Botryosphaeria obtusa on Apple and Management of Stem Cankers with Fungicides. Plant Dis. 84, 1031–1037 (2000).
Urbez-Torres, J. R. & Gubler, W. D. Pathogenicity of Botryosphaeriaceae species isolated from grapevine cankers in California. Plant Dis. 93, 584–592. https://doi.org/10.1094/PDIS-93-6-0584 (2009).
Larignon, P., Fulchic, R., Cere, L. & Dubos, B. Observation on black dead arm in French vineyards. Phytopathol. Mediterr. 40, 336–342. https://doi.org/10.14601/Phytopathol_Mediterr-1629 (2001).
Luque, J., Martos, S., Aroca, A., Raposo, R. & Garcia-Figueres, F. Symptoms and fungi associated with declining mature grapevine plants in northeast Spain. J. Plant Pathol. 91, 381–390. https://doi.org/10.4454/JPP.V91I2.968 (2009).
Savocchia, S., Steel, C. C., Stodart, B. J. & Somers, A. Pathogenicity of Botryosphaeria species isolated from declining grapevines in sub tropical regions of Eastern Australia. Vitis-geilweilerhof 46, 27 (2007).
Zietkiewicz, E., Rafalski, A. & Labuda, D. Genome fingerprinting by simple sequence repeat (SSR)-anchored polymerase chain reaction amplification. Genomics 20, 176–183. https://doi.org/10.1006/geno.1994.1151 (1994).
Hantula, J., Dusabenyagasani, M. & Hamelin, R. C. Random amplified microsatellites (RAMS) ? a novel method for characterizing genetic variation within fungi. For. Pathol. 26, 159–166. https://doi.org/10.1111/j.1439-0329.1996.tb00720.x (1996).
Phillips, A. J. L. Excoriose, cane blight and related diseases of grapevines: A taxonomic review of the pathogens. Phytopathol. Mediterr. https://doi.org/10.1400/57866 (2000).
Leslie, J. F. Fungal vegetative compatibility. Annu. Rev. Phytopathol. 31, 127–150. https://doi.org/10.1146/annurev.py.31.090193.001015 (1993).
Martín, L., Sáenz de Miera, L. E. & Martín, M. T. Applicability of AFLP fingerprinting markers to molecular discrimination of the three main fungal species associated with grapevine decline in Spain. J. Phytopathol. 161, 689–695. https://doi.org/10.1111/jph.12118 (2013).
Mohali, S., Burgess, T. I. & Wingfield, M. J. Diversity and host association of the tropical tree endophyte Lasiodiplodia theobromae revealed using simple sequence repeat markers. For. Pathol. 35, 385–396. https://doi.org/10.1111/j.1439-0329.2005.00418.x (2005).
Baskarathevan, J., Jaspers, M. V., Jones, E. E., Cruickshank, R. H. & Ridgway, H. J. Genetic and pathogenic diversity of Neofusicoccum parvum in New Zealand vineyards. Fungal Biol. 116, 276–288. https://doi.org/10.1016/j.funbio.2011.11.010 (2012).
Billones-Baaijens, R., Jones, E. E., Ridgway, H. J. & Jaspers, M. V. Virulence affected by assay parameters during grapevine pathogenicity studies with Botryosphaeriaceae nursery isolates. Plant Pathol. 62, 1214–1225. https://doi.org/10.1111/ppa.12051 (2013).
Gramaje, D., Armengol, J. & Ridgway, H. J. Genetic and virulence diversity, and mating type distribution of Togninia minima causing grapevine trunk diseases in Spain. Eur. J. Plant Pathol. 135, 727–743. https://doi.org/10.1007/s10658-012-0110-636 (2013).
Phillips, A., Alves, A., Correia, A. & Luque, J. Two new species of Botryosphaeria with brown, 1-septate ascospores and Dothiorella anamorphs. Mycologia 97, 513–529. https://doi.org/10.3852/mycologia.97.2.513 (2005).
Arzanlou, M. & Dokhanchi, H. Morphological and molecular characterization of Diplodia seriata the causal agent of canker and twig dieback disease on mulberry in Iran. Arch. Phytopathol. Plant Prot. 46, 682–694. https://doi.org/10.1080/03235408.2012.749703 (2013).
Alves, A., Correia, A., Luque, J. & Phillips, A. Botryosphaeria corticola, sp. Nov. on Quercus species, with notes and description of Botryosphaeria stevensii and its anamorph, Diplodia mutila. Mycologia 96, 598–613. https://doi.org/10.2307/3762177 (2004).
Alves, A., Correia, A. & Phillips, A. J. L. Multi-gene genealogies and morphological data support Diplodia cupressi sp. Nov., previously recognized as D. pinea f. sp. cupressi, as a distinct species. Fungal Divers. 23, 1–15 (2006).
Damm, U., Crous, P. W. & Fourie, P. H. Botryosphaeriaceae as potential pathogens of Prunus species in South Africa, with descriptions of Diplodia africana and Lasiodiplodia plurivora sp. nov. Mycologia 99, 664–680 (2007).
Lazzizera, C., Frisullo, S., Alves, A., Lopes, J. & Phillips, A. J. L. Phylogeny and morphology of Diplodia species on olives in Southern Italy and description of Diplodia olivarum sp. nov. Fungal Divers. https://www.researchgate.net/publication/230583054 (2008).
Barr, M. E. Prodromus to class Loculoascomycetes. (Lubrecht & Cramer, Limited: Amherst, MA, USA, 2004). ISBN 0-934454-51-5.
Slippers, B. & Wingfield, M. J. Botryosphaeriaceae as endophytes and latent pathogens of woody plants: Diversity, ecology and impact. Fungal Biol. Rev. 21, 90–106. https://doi.org/10.1016/j.fbr.2007.06.002 (2007).
Phillips, A. J. L. et al. The Botryosphaeriaceae: Genera and species known from culture. Stud. Mycol. 76, 51–167. https://doi.org/10.3114/sim0021 (2013).
Martos, S. et al. Production of phytotoxic metabolites by five species of Botryosphaeriaceae causing decline on grapevines, with special interest in the species Neofusicoccum luteum and N. parvum. Eur. J. Plant Pathol. 121, 451–461. https://doi.org/10.1007/s10658-007-9263-0 (2008).
Guo, L. D., Xu, L., Zheng, W. H. & Hyde, K. D. Genetic variation of Alternaria alternata, an endophytic fungus isolated from Pinus tabulaeformis as determined by random amplified microsatellites (RAMS). Fungal Divers. 16, 53–65 (2004).
Bayraktar, H. & Dolar, F. S. Genetic diversity of wilt and root rot pathogens of chickpea, as assessed by RAPD and ISSR. Turkish J. Agric. For. 33, 1–10. https://doi.org/10.3906/tar-0709-37 (2009).
Fan, J., Guo, L., Xu, J., Luo, Y. & Michailides, T. J. Genetic diversity of populations of Monilinia fructicola (Fungi, Ascomycota, Helotiales) from China. J. Eukaryot. Microbiol. 57, 206–212. https://doi.org/10.1111/j.1550-7408.2009.00467.x (2010).
Padder, B. A. et al. Virulence and molecular diversity of Venturia inaequalis in commercial apple growing regions in Kashmir. J. Phytopathol. 161, 271–279. https://doi.org/10.1111/jph.12061 (2013).
Dinolfo, M. I. et al. ISSR markers detect high genetic variation among Fusarium poae isolates from Argentina and England. Eur. J. Plant Pathol. 127, 483–491. https://doi.org/10.1007/s10658-010-9613-1 (2010).
Sudisha, J., Kumar, S. A., Thakur, R. P., Rao, V. P. & Shetty, H. S. Molecular characterization of Sclerospora graminicola, the incitant of pearl millet downy mildew using ISSR markers. J. Phytopathol. 157, 748–755. https://doi.org/10.1111/j.1439-0434.2009.01547.x (2009).
Villarino, M., Larena, I., Martinez, F., Melgarejo, P. & De Cal, A. Analysis of genetic diversity in Monilinia fructicola from the Ebro Valley in Spain using ISSR and RAPD markers. Eur. J. Plant Pathol. 132, 511–524. https://doi.org/10.1007/2Fs10658-011-9895-y (2012).
Tuite, J. Plant Pathologist Methods, Fungi and Bacteria. (Burgess Publishing, Minieapolis, USA, 1969). https://doi.org/10.4236/ajps.2018.913187
Spiers, A. G. Botryosphaeria dothidea infection of Salix species in New Zealand. Plant Dis. Rep. 61, 664–667. https://biotanz.landcareresearch.co.nz/references/1cb0e1bd-36b9-11d5-9548-00d0592d548c (1977).
Murray, M. G. & Thompson, W. F. Rapid isolation of high molecular weight plant DNA. Nucl. Acids Res. 8, 4321–4326. https://doi.org/10.1093/nar/8.19.4321 (1980).
Rohlf, F. J. NTSYS: numerical taxonomy and multivariate analysis system version 2.02. Applied Biostatistics Inc., Setauket, NY (1998).https://doi.org/10.4236/psych.2010.12018
Aneja, K. R. Experiments in Microbiology, Plant Pathology, Tissue Culture and Mushroom Production Technology. 3rd edn. 10–15 (New Age Intern. Pvt. Ltd. Pub., New Delhi, 2001).
Nei, M. & Li, W. H. Mathematical model for studying genetic variation in terms of restriction endonucleases. Proc. Natl. Acad. Sci. 76, 5269–5273. https://doi.org/10.1073/pnas.76.10.5269 (1979).
Falush, D., Stephens, M. & Pritchard, J. K. Inference of population structure using multilocus genotype data: Linked loci and correlated allele frequencies. Genetics 164, 1567–1587. https://doi.org/10.1093/genetics/164.4.1567 (2003).
Evanno, G., Regnaut, S. & Goudet, J. Detecting the number of clusters of individuals using the software STRUCTURE: A simulation study. Mol. Ecol. 14, 2611–2620. https://doi.org/10.1111/j.1365-294X.2005.02553.x (2005).
Excoffier, L., Laval, G. & Schneider, S. Arlequin (version 3.0): An integrated software package for population genetics data analysis. Evol. Bioinforma. 1, 117693430500100. https://doi.org/10.1177/117693430500100003 (2005).
Yeh, F. & Boyle, T. Population genetic analysis of co-dominant and dominant markers and quantitative traits. Belgian J. Bot. 129, 57 (1997).
Pritchard, J. K., Stephens, M. & Donnelly, P. Inference of population structure using multilocus genotype data. Genetics 155, 945–959. https://doi.org/10.1093/genetics/155.2.945 (2000).
Peakall, R. & Smouse, P. E. GenAlEx 6.5: Genetic analysis in Excel. Population genetic software for teaching and research—an update. Bioinformatics 28, 2537–2539. https://doi.org/10.1093/bioinformatics/bts460 (2012).
Slippers, B. et al. Phylogenetic lineages in the Botryosphaeriales: A systematic and evolutionary framework. Stud. Mycol. 76, 31–49. https://doi.org/10.3114/sim0020 (2013).
Dissanayake, A. Botryosphaeriaceae: Current status of genera and species. Mycosphere 7, 1001–1073. https://doi.org/10.5943/mycosphere/si/1b/13 (2016).
Crous, P. W. et al. Phylogenetic lineages in the Botryosphaeriaceae. Stud. Mycol. 55, 235–253. https://doi.org/10.3114/sim.55.1.235 (2006).
Epstein, L., Sukhwinder, K. & VanderGheynst, J. S. Botryosphaeria-related dieback and control investigated in noncoastal California grapevines. Calif. Agric. 62, 161–166. https://doi.org/10.3733/ca.v062n04p161 (2008).
Earl, D. A. & VonHoldt, B. M. STRUCTURE HARVESTER: A website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv. Genet. Resour. 4, 359–361. https://doi.org/10.1007/s12686-011-9548-7 (2012).
Acknowledgements
The Article Processing Charge is financed by Wroclaw University of Environmental and Life Science. The authors extend their appreciation to the Researchers supporting project number (RSP2023R306), King Saud University, Riyadh, Saudi Arabia.
Funding
The Article Processing Charge is financed by Wroclaw University of Environmental and Life Science. The authors extend their appreciation to the Researchers supporting project number (RSP2023R306), King Saud University, Riyadh, Saudi Arabia.
Author information
Authors and Affiliations
Contributions
Conceptualization, A.H.B and M.A.A.; methodology, A.H.B., F.M.A. and P.S.; software, A.H.B.; B.A.P. and U.F.; validation, A.H.B.; M.D.S., B.A.P., E.A.A.D., and M.A.A.; formal analysis, A.H.B. and P.S.; investigation, A.H.B. and M.D.S.; resources, M.D.S., B.A.P. and E.A.A.S.; data curation, A.H.B. Z.A.S and U.F.; writing—original draft preparation, A.H.B, Z.A.S, E.A.D, P.S., M.A.A.A. and E.A.A.S.; writing—review and editing, A.H.B., E.A.D, M.S.N., M.A.A., F.M.A. and P.S.; visualization, Z.A.S, E.A.D and M.S.N.; supervision, M.D.S. and B.A.P.; project administration, A.H.B. and funding acquisition, M.A.A., F.M.A. and P.S. All authors have read and agreed to the published version of the manuscript.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Bhat, A.H., Shah, M.D., Padder, B.A. et al. Morphological, pathogenic and genetic diversity in Diplodia seriata associated with black rot canker of apple in India. Sci Rep 13, 15682 (2023). https://doi.org/10.1038/s41598-023-42534-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-023-42534-y
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.