What Do We Know about Botryosphaeriaceae? An Overview of a Worldwide Cured Dataset

Botryosphaeriaceae-related diseases occur worldwide in a wide variety of plant hosts. The number of studies targeting the distribution, diversity, ecology, and pathogenicity of Botryosphaeriaceae species are consistently increasing. However, with the lack of consistency in species delimitation, the name of hosts, and the locations of studies, it is almost impossible to quantify the presence of these species worldwide, or the number of different host–fungus interactions that occur. In this review, we collected and organized Botryosphaeriaceae occurrences in a single cured dataset, allowing us to obtain for the first time a complete perspective on species’ global diversity, dispersion, host association, ecological niches, pathogenicity, communication efficiency of new occurrences, and new host–fungus associations. This dataset is freely available through an interactive and online application. The current release (version 1.0) contains 14,405 cured isolates and 2989 literature references of 12,121 different host–fungus interactions with 1692 different plant species from 149 countries.


Introduction
Species of Botryosphaeriaceae (Botryosphaeriales, Ascomycetes) are distributed worldwide and are known to have different ecological roles. These fungi can act as saprobic, endophytic, or latent pathogens [1,2]. Some members of this family are recognized as aggressive plant pathogens on different types of hosts. From agricultural crops to ornamental and forest hosts, these fungi have no boundaries [3][4][5][6][7][8]. Their wide distribution; ability to persist endophytically; ability to become pathogenic only when their hosts are under stress, causing diseases that eventually may lead to host death; and capability to adapt and colonize new hosts [2,9] makes these organisms into a great challenge for plant pathology in a changing world [10][11][12].
Several species of Botryosphaeriaceae currently accepted have been described in the XIX century, such as, for example, Diplodia mutila (1834), Diplodia seriata (1845), Botryosphaeria dothidea (1863), and Diplodia sapinea (1870), among others. At that time, the description of new species, as well as the taxonomic position of the family Botryosphaeriaceae, was exclusively based on their micromorphological characteristics. Therefore, over the years the taxonomic position of these organisms suffered multiple classifications (for a historical overview, see [1]).
With the use of DNA sequencing and phylogenetic methods, our understanding of the taxonomy and diversity of the Botryosphaeriaceae species has changed profoundly over time. Since 1996, several authors using sequence data have studied the phylogenetic relationships of this family, and currently 20 genera and 280 species have been described [1,[13][14][15][16][17][18][19].
Extensive surveys reporting the diversity, distribution, and pathogenicity of Botryosphaeriaceae species have been performed in several countries, providing valuable information in terms of the frequency and diversity of hosts. These include Australia [20], Algeria [21],

MDRBOT Database and Shiny Interface
A literature review was performed, and both datasets were analyzed to fill missing information. Several extra organization-level features were built to allow different filter functions. A shiny interface was created to allow easy access to both datasets.

The Site
The MDRBOT database was built using R 3.6.0 with a web shiny interface. The website includes a search engine to find the cured Botryosphaeriaceae isolates from the Nucleotide-NCBI database and for the host-fungus interaction dataset, where the user can perform multiple field search and download the output in an excel format file (.xlsx). A worldwide map generator tool for species occurrence is available, where the user can select an input species and generate a world occurrence map. A climate data analysis is included, where the user can select an input species and observe isolates with valid geo-

MDRBOT Database and Shiny Interface
A literature review was performed, and both datasets were analyzed to fill missing information. Several extra organization-level features were built to allow different filter functions. A shiny interface was created to allow easy access to both datasets.

The Site
The MDRBOT database was built using R 3.6.0 with a web shiny interface. The website includes a search engine to find the cured Botryosphaeriaceae isolates from the Nucleotide-NCBI database and for the host-fungus interaction dataset, where the user can perform multiple field search and download the output in an excel format file (.xlsx). Despite all the efforts to characterize this fungal family, it is impossible to evaluate the global diversity of these organisms due to the different levels of sampling effort across countries. For this reason, we compared diversity with countries' sampling efforts using the location of isolates from the Nucleotide-GenBank collection. In this analysis, we consider that data on nucleotide are representative for the sampling effort. We found that 138 countries still have no records of Botryosphaeriaceae species deposited in GenBank (represented in white in Figure 2) (17), and Portugal (23). In Figure 2, we produced a bi-variate world map to evaluate countries' diversity vs. sampling effort. The lack of data for several regions of the globe is clear, predominantly in Africa; Central and Western Europe; North, Central, and Western Asia; Bolivia; and several countries in the Caribbean region. Despite all the efforts to characterize this fungal family, it is impossible to evaluate the global diversity of these organisms due to the different levels of sampling effort across countries. For this reason, we compared diversity with countries' sampling efforts using the location of isolates from the Nucleotide-GenBank collection. In this analysis, we consider that data on nucleotide are representative for the sampling effort. We found that 138 countries still have no records of Botryosphaeriaceae species deposited in GenBank (represented in white in Figure 2) and 66% of the countries have records of less than 51 isolates.  (17), and Portugal (23). In Figure 2, we produced a bi-variate world map to evaluate countries' diversity vs. sampling effort. The lack of data for several regions of the globe is clear, predominantly in Africa; Central and Western Europe; North, Central, and Western Asia; Bolivia; and several countries in the Caribbean region. Bi-variate world map analyzing diversity vs. sampling effort for Botryosphaeriaceae isolates. Data obtained from Nucleotide-GenBank (a total of 14,405 isolates: 12,593 with location and 1812 isolates without location). Countries in both variables were divided in quartiles with an equal probability of occurrence (e.g., for species diversity, the probability of randomly selecting a country is equal for the intervals with 1-3, 4-8, and 9-85 species) to facilitate visualization.

Worldwide Occurrence-From Where to Where?
Evolutionary divergence studies show that Botryosphaeriaceae lineages emerged during the late Cretaceous period, over 66 million years ago, in a period dominated by the expansion of angiosperms occupying environments previously dominated by conifers. It is hypothesized that the evolution of modern Botryosphaeriaceae species was driven by Bi-variate world map analyzing diversity vs. sampling effort for Botryosphaeriaceae isolates. Data obtained from Nucleotide-GenBank (a total of 14,405 isolates: 12,593 with location and 1812 isolates without location). Countries in both variables were divided in quartiles with an equal probability of occurrence (e.g., for species diversity, the probability of randomly selecting a country is equal for the intervals with 1-3, 4-8, and 9-85 species) to facilitate visualization.

Worldwide Occurrence-From Where to Where?
Evolutionary divergence studies show that Botryosphaeriaceae lineages emerged during the late Cretaceous period, over 66 million years ago, in a period dominated by the expansion of angiosperms occupying environments previously dominated by conifers. It is hypothesized that the evolution of modern Botryosphaeriaceae species was driven by the evolution and diversification during the Paleocene epoch of their hosts, which we currently know as modern plants [18,35]. Several authors have investigated the origin of some species by performing population studies at the global scale among different hosts and did not identify an obvious phylogeographic origin [6,[36][37][38][39][40][41].
It is known that the infection and colonization of hosts can occur through natural wounds on leaves, branches, or stems and due to other openings such as lenticels and stomata. The distribution of these organisms is favored by the sticky spores dispersed by wind, rain, and insects [2]. At the intercontinental level, human movement, and the international trade of plants and derivates (timber and non-timber products) without appropriate quarantine measures, has caused the dispersion of these organisms [30]. Within this family, several species are known to have limited distributions, and a few species such as Botryosphaeria dothidea, Diplodia sapinea, Diplodia seriata, Dothiorella sarmentorum, Neofusicoccum parvum, and Lasiodiplodia theobromae are recognized to be globally distributed [1,17,41]. To understand the distribution of the main Botryosphaeriaceae species, we analyzed their spatial distributions over continents and terrestrial ecoregions [42,43]. Within the studied species, so far only Diplodia corticola and Neofusicoccum mangiferae have not been reported in every continent. This observation confirms the ability of the remaining species to spread globally. However, if we take into consideration the different terrestrial ecoregions, we can observe that some species are reported only in certain ecoregions. Based on occurrence data, we suggest a latitudinal shift among different types of climates where some species are clearly concentrated in some types of ecosystems. For example, Do. Sarmentorum, D. corticola, D. mutila, N. australe, D. seriata, and D. sapinea are found only in temperate and Mediterranean ecosystems, while others with a wider range of dispersion. Neofusicoccum parvum in our analysis appears to be the most adapted organism, with it being detected from the north to the south, with the exception of boreal forests and montane grasslands. Boreal forests appear to be the most unlikely place to find Botryosphaeriaceae species; so far, only D. sapinea has been reported in these regions ( Figure 3).

Understanding the Process of Host Jumps-Can We Spot Host Specificity?
Defining a host range and understanding which factors favor future host jumps is essential to study emerging and re-emerging fungal pathogens. Several drivers are often mentioned in the literature, such as international trade, the failure of quarantine and preventive measures, changes in land use or agricultural practices, pathogen evolution and plasticity, mechanisms of genome divergence (e.g., mutation, hybridization, sexual recombination, horizontal gene transfer, and others), host-fungus genotype-by-genotype interactions, poor host health, and climate change, among others [44][45][46][47][48][49]. Comparative genomics and omics studies are slowing unveiling host-fungus interaction mechanisms by dissecting plant defence mechanisms, fungal pathogenic strategies, and nutrient uptake pathways [50][51][52][53][54]. To clarify a momentary host range boundary and spot host specificity, a complete overview of all the mentioned areas is crucial.
Over time, our knowledge regarding the biology and ecology of these pathogens has improved significantly. However, the mechanisms behind host jumps and the worldwide extension of host association patterns across Botryosphaeriaceae genera are still poorly understood. In general, the most studied species lack host specificity and have a proven capability to colonize and cause disease in diverse native and introduced plant hosts [2,7-9, [55][56][57][58]. Even species such as D. sapinea and D. corticola, which have been consistently associated with a certain type of host, have been occasionally found to occur in other unrelated hosts in different regions of world [9,[59][60][61].
To investigate possible host specificity, we compared the number of hosts against the total number of hosts of the genera, the number of isolates, the number of countries where these species were detected, the number of reports found in the literature, and differences between the numbers of known associations with angiosperms and gymnosperms ( Figure 4). Overall, L. theobromae is by far the organism with the largest number of known hosts (666 of 749 hosts reported for the genus Lasiodiplodia), with the largest number of isolates in GenBank (1944), the largest number of country occurrences (97), and the largest number of host-fungus interaction reports (365). For the remaining species, it is possible to observe, even with different research efforts in a number of isolates or literature reports, that some species clearly have the capability to colonize a higher number of hosts than other (e.g., N. parvum with 223 hosts in 50 countries and B. dothidea with 403 hosts in 66 countries), and few have been consistently reported worldwide in a low number of hosts (e.g., D. sapinea was reported in 62 countries but only in 102 hosts, of which 83 are gymnosperms, or D. seriata with 121 hosts in 46 countries and Lasiodiplodia pseudotheobromae with 124 hosts in 44 countries).
However, a higher number of different hosts might be related only to the diversity of species among different plant functional groups and not to their ability to colonize different hosts from different plant groups. For example, the Pinopsida class-also known as conifers-has only 615 living species, while the Liliopsida (monocotyledons) and the Magnoliopsida (dicotyledons) classes have an estimated 77,267 and 246,366 living species, respectively. To analyze host jumps among different functional groups, we created a circular graph with shared hosts-fungus interactions where hosts were divided in different plants groups ( Figure 5). Not surprisingly, L. theobromae showed a capability to colonize hosts in all plant groups. However, L. pseudotheobromae and N. parvum, with a considerably lower number of known hosts, present a similar pattern of shared host-fungus interactions, suggesting their clear capability to increase their host range in a similar way to that of L. theobromae. For other fungi, it is possible to explore range expansions: N. mangiferae and L. iraniensis have recently been described in several members of the Malvids group and are occasionally found in other groups. Diplodia sapinea has been consistently reported in the Pinales species and D. corticola in the Quercus species from the Fabids group, but both are also starting to be found in different plant groups. These types of analyses can be useful for institutions and national authorities to guide their studies and to anticipate future host jumps, especially for highly pathogenic organisms. Notably, few studies have proven host specialization among different Botryosphaeriaceae species based on differences in pathogenicity-related genes (e.g., B. dothidea versus Botryosphaeria kuwatsukai) [62].

How Much We Know about Pathogenicity and Plant Mortality?
Uncovering the complexity of host-pathogen interaction is not a stationary science and depends on multiple variables from the environment and the interaction of both hosts' and pathogens' genome artillery [63,64]. Host-pathogen interactions are similar to a chess game, where for a specific outcome several actions are possible from each organism. The reproducibility of these interactions under controlled conditions might not always be representative of what occurs in nature, and common pathogenicity trials do not fully expose pathogenic and resistance mechanisms [65,66]. Nonetheless, the sequence and annotation of both genomes opens the possibility for multi-omics analyses to provide a more complete overview of these interactions [53,67].
The pathogenicity of Botryosphaeriaceae isolates leading to plant mortality has been demonstrated under controlled conditions, mainly in seedlings, in well-watered conditions or under drought stress [9,68,69]. However, plant mortality in nature is often a combination of multiple biotic and abiotic stresses. On one hand, abiotic factors such as drought or heat stress can disrupt plant physiological performance, allowing the colonization of fungal pathogens and increasing disease susceptibility. On the other hand, colonization by fungal pathogens can reduce tolerance to biotic stress, leading to higher mortality rates [63,[70][71][72]. The interaction effects of drought and fungal pathogen infection on plant mortality are shown in Figure 6.
Several factors can induce mortality, such as carbon starvation when non-structural carbohydrates resources are depleted, affecting normal plant maintenance, growth, and defence mechanisms [73]; the toxic effect of fungal metabolites; and the hydraulic failure of xylem water transport to the leaves due to embolism or phloem transport caused by an impaired xylem water potential [74]. These factors can be directly or indirectly induced by both biotic and abiotic stresses. For example, drought and heat stress can induce stomata closure and simultaneously decrease carbon assimilation, which is essential in maintaining plant defence metabolism and functional sapwood maintenance. Pathogen infection can directly impact carbon assimilation by down-regulating the genes involved in photosynthetic activities or triggering carbon starvation by inducing tree defences and/or inhibiting the expression of genes involved in carbon metabolism and transport. Additionally, when these pathogens colonize plant vascular tissues and vascular necrosis occurs, the whole-plant hydraulic conductance is reduced, increasing the risk of hydraulic failure [74]. Furthermore, some species were proven to have the capability to exhibit cytotoxicity against mammalian cell lines, and again temperature was suggested to modulate the expression of toxic compounds [54,79,80].

Climate Sensitivity, a Hidden Pattern?
Temperature growth studies suggest that, in general, Botryosphaeriaceae species present minimal growth rates at 5 °C or over 35 °C and optimal growth rates at around Genome and transcriptome analyses of some Botryosphaeriaceae species have shown a higher number of pathogenicity-related genes associated with cell wall degradation, nutrient uptake, secondary metabolism, and membrane transport functions, which are important for woody plant infection when compared with other fungi with different lifestyles [62,75]. Additionally, several authors have demonstrated that, during infection, gene families related to carbohydrate catabolism; pectin, starch, and sucrose metabolism; pentose and glucuronate interconversion pathways were induced [76][77][78]. Some of these genes' families were even induced by higher temperatures [54,75].
Furthermore, some species were proven to have the capability to exhibit cytotoxicity against mammalian cell lines, and again temperature was suggested to modulate the expression of toxic compounds [54,79,80].

Climate Sensitivity, a Hidden Pattern?
Temperature growth studies suggest that, in general, Botryosphaeriaceae species present minimal growth rates at 5 • C or over 35 • C and optimal growth rates at around 15-25 • C under controlled conditions [1,17]. So far, no major studies have been performed to characterize the natural bioclimatic envelopes of Botryosphaeriaceae species. Based on the geographical coordinates obtained during this review, we analyzed the climatic variability of these records (Figure 7). In terms of annual mean temperature, it is clear that B. dothidea, D. seriata, and Do. iberica are often collected in places with lower annual mean temperatures when compared with L. theobromae, M. phaseolina, N. parvum, and Ne. dimidiatum. For annual precipitation, this pattern is not so clear, with the exception of L. theobromae, where the majority of records were collected in places with more than 1000 mm of annual precipitation (typical of tropical regions). Temperature and precipitation seasonality were calculated using the standard deviation of the mean monthly values. Regions with larger standard deviations have greater temperature and precipitation variability across the year. These metrics are important in understanding the tolerance range for species in the future. For example, a species that is often present in tropical regions with a tolerance to higher annual mean temperatures but without tolerance to strong temperature seasonality may have problems expanding its distribution range in the future to a region with a strong temperature variability across the year, even if the mean temperature increases (e.g., can L. theobromae, often collected in tropical regions, expand its range to Mediterranean and temperate regions with climate change?) These results are consistent with our ecoregions analysis, suggesting that probably the known host-fungus interactions and species distributions are strongly dependent on climate. We encourage authors to provide precise coordinates for occurrence data to improve the understanding of the Botryosphaeriaceae species distribution.

Global Dispersion-How Far Can They Go? Framing Ecological Niche Requirements for Potential Species Distributions Areas
As mentioned before, the Botryosphaeriaceae species produce spores that are naturally dispersed through wind, rain, insects. However, human-induced activities are responsible for the majority of long-distance dispersion [2,30]. Taking into account the large number of potential hosts and the large quantity of plant material moving worldwide every day, it is virtually impossible to verify and detect efficiently latent pathogens living endophytically in symptomatic or asymptomatic material [30,81,82]. Thus, understanding the ecological niche requirements for potential species distribution areas might be a better solution to predict and prevent future outbreaks [83].
We propose an adaptation of the classical BAM diagram ( Figure 8I) to explain the influence of environmental conditions, biotic interactions, and dispersal in shaping species' geographic distribution [85][86][87]. Following the classical BAM diagram, B represents the geographic regions where the interaction factors with other species are favorable for species occurrence, A represents the geographical regions where the climatic conditions are favorable for maintaining a viable long-term population, and M corresponds to the geographical region accessible for species dispersion. However, for endophytes and latent pathogens, this perspective does not fully represent the ecological dynamics of these organisms. We assume that: (1) endophytic latent pathogens are mainly dispersed by human activities such as movement and trade; (2) the introduction of a species in a new environment is likely to occur, as human movement/trade exists and is favored by a lack of preventive and quarantine measures, therefore M is virtually unlimited; (3) the establishment of these species is affected by climate. Nonetheless, unfavorable conditions might hide the presence of those species in asymptomatic hosts or by resistance structures. Seasonal effects might expand or decrease the growth of these organisms, invalidating viable long-term populations. (4) Optimal conditions for disease expression are mainly occasional climatic events that can affect the susceptibility of the host (i.e., a reduction in precipitation or/and temperature increments could cause drought or heat stress to the host [70,88,89]) or/and trigger the pathogenic behavior of these organisms (i.e., variations in temperature, light intensity, or atmospheric ozone could induce phytopathogenic mechanisms [54,79,80,90,91]), (5) Only when host specificity is demonstrated, whether for a fungal species with a limited ability to colonize and persist endophytically in certain hosts or/and for species with limited ability to infect and express disease symptoms in a certain type of hosts, is it assumed that biotic interactions (B) can shape the geographical distribution ( Figure 8II).

Figure 7.
Annual mean temperature and precipitation and temperature and precipitation seasonality for the main Botryosphaeriaceae species occurrence. Data were thinned by 200 km distance to avoid spatial autocorrelation using the R package ELSA (Entropy-Based Local Indicator of Spatial Association) [84]. Climatic variables were extracted from the CHELSA project. Number of valid coordinates are indicated per species.
As a result, for non-specific endophytes and latent pathogens, like many of the Botryosphaeriaceae species, assuming an imperfect quarantine system worldwide, climate is the main variable constraining the geographical distribution of these organisms. or/and trigger the pathogenic behavior of these organisms (i.e., variations in temperature, light intensity, or atmospheric ozone could induce phytopathogenic mechanisms [54,79,80,90,91]), (5) Only when host specificity is demonstrated, whether for a fungal species with a limited ability to colonize and persist endophytically in certain hosts or/and for species with limited ability to infect and express disease symptoms in a certain type of hosts, is it assumed that biotic interactions (B) can shape the geographical distribution ( Figure 8II).

How Good Are We at Reporting New Occurrences and Host-Fungus Associations?
Despite the increase in new Botryosphaeriaceae-related studies worldwide, there is still a lack of standard databases that are consistently curated and maintained over time. The failure in efficiently reporting new occurrences and host-fungus associations increases the gap between science and society and dilutes the scientific effort to improve preventive and quarantine policy measures. Information regarding this fungal family is often reported in indexed scientific journals and should be followed by public sequence data that allow the scientific community to validate and confirm the taxonomic identification of these organisms. If information only present in the literature might be outdated and difficult to verify, information only presented in GenBank or similar databases is also susceptible to be lost without a proper report. Ideally, a report of a new occurrence or a new host-fungus interaction should be documented in the literature and supported with public genomic data. These data should be well-organized in public databases to allow national institutions such as customs or agricultural/forestry authorities to be prepared for possible new threats.
In Figure 9, we compare the country occurrences and host-fungus interactions by country that are reported only in the literature, in both the literature and GenBank, and only reported in GenBank. We also classified countries by the percentage of information that is not properly reported. We found that only 53.69% of the species occurrence by country is reported in both datasets, and if we take in consideration host-fungus interactions by country, only 23.07% of the current knowledge is well reported. Records prior to the massification of DNA-based methodologies (older than the year 2000) are often poorly reported, but there is also a large number of recent studies that fail to provide consistent information about host and location even when public genomic data are available.

Conclusions and Future Perspectives
Our database represents a single effort to clean and organize all Botryosphaeriaceaerelated occurrences. This database will be continuously maintained and researchers working with these organisms are welcome to submit or update their information. Major updates on the source information are expected to happen every January and new features may be added over time as a result of users' feedback.
With the insights of a worldwide cured dataset, we verify a consistent growing interest of these plant pathogenic fungi when taking in account the number of public records and publications and the citation history of papers over the years [30]. However, we demonstrate that 80% of all isolates with public sequence data are concentrated only in 11 countries, leading to a huge knowledge gap regarding Botryosphaeriaceae's occurrence and diversity worldwide. Additionally, even countries with a higher sampling effort are still far from reaching a plateau for species diversity versus sampling effort, suggesting a high number of undescribed species. The problematic of estimating species numbers is a common issue in other fungal families [92], raising the concern of fully understanding the limit of a species boundary and properly identifying cryptic species on species interactions, as well the recognition of the hybridization phenomenon [36,93,94].
In terms of ecology, we reflect about the climate influence on the distribution and dispersion of these organisms for the first time. Despite the worldwide distribution across all continents, with the exception of Antarctica, it was possible to observe a climatic preference for some species (e.g., D. corticola, Do. sarmentorum, or N. australe) in temperate and Mediterranean regions and other species with a wider tolerance in temperate to tropical regions (e.g., B. dothidea, N. parvum, or M. phaseolina). In terms of dispersion, we propose a new framework to define the ecological niche requirements for most fungal latent endophytes. This framework is essential to improve our understanding of the current species distributions areas and to forecast future disease outbreaks [83]. We assume that human movement and trade are the main dispersion routes of these organisms, with climate variability being the major constrain for the occurrence of new stable populations. Additionally, we highlight that disease expression is mainly due to occasional climatic events that

Conclusions and Future Perspectives
Our database represents a single effort to clean and organize all Botryosphaeriaceaerelated occurrences. This database will be continuously maintained and researchers working with these organisms are welcome to submit or update their information. Major updates on the source information are expected to happen every January and new features may be added over time as a result of users' feedback.
With the insights of a worldwide cured dataset, we verify a consistent growing interest of these plant pathogenic fungi when taking in account the number of public records and publications and the citation history of papers over the years [30]. However, we demonstrate that 80% of all isolates with public sequence data are concentrated only in 11 countries, leading to a huge knowledge gap regarding Botryosphaeriaceae's occurrence and diversity worldwide. Additionally, even countries with a higher sampling effort are still far from reaching a plateau for species diversity versus sampling effort, suggesting a high number of undescribed species. The problematic of estimating species numbers is a common issue in other fungal families [92], raising the concern of fully understanding the limit of a species boundary and properly identifying cryptic species on species interactions, as well the recognition of the hybridization phenomenon [36,93,94].
In terms of ecology, we reflect about the climate influence on the distribution and dispersion of these organisms for the first time. Despite the worldwide distribution across all continents, with the exception of Antarctica, it was possible to observe a climatic preference for some species (e.g., D. corticola, Do. sarmentorum, or N. australe) in temperate and Mediterranean regions and other species with a wider tolerance in temperate to tropical regions (e.g., B. dothidea, N. parvum, or M. phaseolina). In terms of dispersion, we propose a new framework to define the ecological niche requirements for most fungal latent endophytes. This framework is essential to improve our understanding of the current species distributions areas and to forecast future disease outbreaks [83]. We assume that human movement and trade are the main dispersion routes of these organisms, with climate variability being the major constrain for the occurrence of new stable populations.
Additionally, we highlight that disease expression is mainly due to occasional climatic events that can affect the susceptibility of the host. Raising the importance to sample asymptomatic hosts for an early detection of new species occurrence (e.g., Diplodia insularis was reported for the first time in Portugal in an asymptomatic host [9]).
To finalize, we evaluate the consistency of known species occurrence and host associations reports in both indexed scientific journals and public sequence databases. We demonstrate the incapability of our society to efficiently use and aggregate data of these emergent plant-pathogens. More than never, we consider that consistent and open plant pathology databases are fundamental to address the challenge of Botryosphaeriaceaerelated diseases in a changing world.