Fungal Diversity in the Neptune Forest: Comparison of the Mycobiota of Posidonia oceanica, Flabellia petiolata, and Padina pavonica

Fungi are widely distributed in the Oceans, interact with other organisms and play roles that range from pathogenic to mutualistic. The present work focuses on the characterization of the cultivable mycobiota associated with the seagrass Posidonia oceanica (L.) Delile collected off the Elba Island (Italy). We identified 102 taxa (mainly Ascomycota) by the mean of a polyphasic approach. Leaves, rhizomes, roots and matte were characterized by unique mycobiota revealing a “plant-part-specificity.” The comparison with the mycobiota associated with the green alga Flabellia petiolata and the brown alga Padina pavonica underlined a “substrate specificity.” Indeed, despite being part of the same phytocoenosis, these photosynthetic organisms recruit different fungal communities. The mycobiota seems to be necessary for the host’s defense and protection, playing, in this way, remarkable ecological roles. Among the 61 species detected in association with P. oceanica (including two species belonging to the newly introduced genus Paralulworthia), 37 were reported for the first time from the Mediterranean Sea.


INTRODUCTION
Posidonia oceanica (L.) Delile, otherwise known as Neptune grass, is the most important seagrass in the Mediterranean Sea (Personnic et al., 2014). This endemic long-lived plant (Magnoliophyta, kingdom Archaeplastida) relies mainly on vegetative reproduction, and its largest clones can spread over 15 km (Arnaud-Haond et al., 2012). Neptune grass meadows, or rather forests, are a spawning ground, a nursery and a permanent habitat for many species: the presence of over 400 different seaweeds species and several thousands of animal species makes this underwater habitat a unique biodiversity hotspot (Boudouresque, 2004).
Rhizomes, roots and senescent leaf sheaths of P. oceanica form a characteristic peat-like seabed layer (matte) that can be several meters thick and thousands of years old. Matte are extremely resistant to microbiological decay, however, herbivores can feed on them with the aid of symbiotic microorganisms such as fungi or bacteria capable of degrading organic matters (Vohnik et al., 2017). To this respect, it is important to underline that likewise other seagrasses, P. oceanica supports a complex ecosystem that supplies large amounts of nutrients to the organisms living within and in adjacent waters. Fungi play an important role in the functioning of this ecosystem: while parasites and patho ens can negatively affect the wellness of the meadows, mutualists may be essential for their health (Vohnik et al., 2019). Finally, saprobes transform dead plants into nutritive detritus through their degrading activity (Raghukumar, 2017). Overall, Posidonia forests are home to 20-25% of Mediterranean species and represent the climax community of the Mediterranean Sea (Minelli et al., 2008). However, in the aftermath of the anthropic impact, of the presence of exotic species and of the climate change, this climax community is seriously threatened and since 1992 is a priority in Annex I of the Council Directive 92/43/EEC on the conservation of natural habitats and of wild fauna and flora.
Marine fungi are a huge part of the microbial diversity hosted in the sea. Indeed, they have been retrieved worldwide from biotic and abiotic substrates such as algae, sediments, invertebrates, drift wood, etc. (Jones and Pang, 2012;Garzoli et al., 2014Garzoli et al., , 2015Garzoli et al., , 2018Gnavi et al., 2017;Raghukumar, 2017;Bovio et al., 2018). Nevertheless, despite the total fungal diversity has been estimated to range between 10,000 and 12,500 taxa (Jones et al., 2019), only 1281 species have been described to date 1 . These organisms are important primary decomposers and living as mutualists (ecto-and endosymbionts), parasites, pathogens and saprobes, contribute greatly to nutrient cycling and food web (Amend et al., 2019;Grossart et al., 2019).
As part of a thorough investigation aimed to uncover the fungal diversity of the Mediterranean Sea, we focused on the isolation and identification of the cultivable mycobiota associated with P. oceanica, collected off the Elba Island (Tuscany). Once defined, the mycobiota was compared to those associated with the green and brown algae Flabellia petiolata  and Padina pavonica , two important components of the same phytocoenosis, collected during the same sampling and in close proximity. Finally, this paper introduces the novel genus Paralulworthia represented by two new species P. gigaspora and P. padinae.

Sampling
Plants of P. oceanica were collected from two sampling sites in the coastal waters of the Elba Island (Livorno, Italy): (1) Ghiaie (UTM WGS84 42 • 49 04 N, 10 • 19 20 E) and (2) Margidore (UTM WGS84 42 • 45 29 N, 10 • 18 24 E). Twelve plants (6 per sampling site) with the surrounding fluctuating senescent leaves (matte) were harvested at depth ranging from 5 to 15 m below sea level (bsl). Plants were maintained in sterile dark containers at 4 • C and processed within 36 h from sampling. Specific permissions to operate in the protected area of "Le Ghiaie" (Ghiaie site) and in the freely accessible Margidore site, were obtained by the port authority of Portoferraio (Livorno, Italy).

Fungal Isolation
Each plant was divided into four parts: leaves, rhizomes, roots and matte (Figure 1). In order to remove unrefined sediments, 1 www.marinefungi.org each sample was individually sonicated (30" each time) and serially washed (three times) in sterilized artificial SeaWater (SW, 3.4% w/v Sea Salt mix -Sigma-Aldrich, St. Louis, MO, United States -in distilled water). Individual samples were homogenized in 20 mL of sterile filtered seawater by means of Ultra-Turrax (IKA, Staufen, Germany). One mL of homogenate was plated onto Petri dishes (15 cm Ø) containing the following culture media: Corn Meal Agar SeaWater (CMASW; 17g CMA -Sigma-Aldrich, St. Louis, MO, United States -in 1 L of filtered SW) or Agar Posidonia (AP; 20 g of P. oceanica heated in 100 mL seawater -60 • C, 30 min -were filtered prior to the addition of 18 g agar and SW up to 1L). Each medium was supplemented with antibiotics (Gentamicin 80 mg L −1 , Piperacillin and Tazobactam 100 mg L −1 -Sigma-Aldrich, St. Louis, United States).
Three replicates per medium and per sample were performed. A total of 288 plates were incubated at 15 • C (spring average temperature of submerged meadows between 5 and 15 m bsl) for 15 days. Afterward, the temperatures was raised to 25 • C for 45 days. This procedure allowed the isolation of both psychrophilic and mesophilic fungi. Colony forming units per gram of dry weight of each plant part (CFU g −1 dw) were recorded. Isolated fungi were maintained in pure culture.

Fungal Identification
Fungi were identified by combining morpho-physiological, molecular and phylogenetic studies. Following preliminary determination of genera according to macroscopic and microscopic features (Domsch et al., 1980;von Arx, 1981;Kiffer and Morelet, 1997), fungal isolates were transferred to the media recommended by the authors of selected genus monographs for species identification. Molecular identification was performed by amplifying and sequencing the appropriate genetic marker as previously described (Poli et al., 2016;Garzoli et al., 2018). PCR products were purified and sequenced at Macrogen Europe Laboratory (Madrid, Spain). The resulting ABI chromatograms were visually inspected, trimmed and assembled to obtain consensus sequences using Sequencer 5.0 (GeneCodes Corporation, Ann Arbor, MI, United States) 2 . The newly generated sequences were compared by BLASTn analyses (default settings) to those available in public nucleotide databases provided by the NCBI (Bethesda MD, United States) and by the Westerdijk Fungal Biodiversity Institute (Utrecht, The Netherlands). Similarity values equal or higher than 98% (e-value > e-100) were considered credible.
Representative strains of each species isolated in pure culture during this work are preserved at Mycotheca Universitatis Taurinensis  Paralulworthia gigaspora and Paralulworthia posidoniae: Growth Condition, Morphological, and Molecular Study Following pre-growth, MUT 435, MUT 672, MUT 5413, and MUT 5261 were inoculated in triplicate onto Petri dishes (9 cm Ø) containing MEASW and incubated at 21 • C. The fungal growth, as well as macro-and microscopic features, was monitored at 4, 7, 14, 21 and 28 days from the inoculum. In an attempt to induce sporulation, sterile specimens of Quercus ruber cork and Pinus pinaster wood (species autochthonous to the Mediterranean area) were placed on 3 weeks old fungal colonies (Panebianco et al., 2002) which were further incubated for 4 weeks at 21 • C. Colonized wood specimens were transferred into 50 mL tubes containing 20 mL sterile sea water and incubated at 21 • C for 1 month. Mature reproductive fungal structures were observed and images captured with an optical microscope (Leica DM4500B, Leica Microsystems GmbH, Germany) equipped with a camera (Leica DFC320, Leica Microsystems GmbH, Germany). For a better inspection of ascomata growing inside wood tissues, scraps of colonized wood embedded in melted paraffin were dissected with a cryostat (Leica CM1850, Leica Microsystems GmbH, Germany).
Genomic DNA was extracted as mentioned above and the partial sequences of the internal transcribed spacers including the 5.8S rDNA gene (nrITS), 28S large ribosomal subunit (nrLSU) and 18S small ribosomal subunit (nrSSU) were amplified by PCR using the primer pairs ITS1/ITS4 (White et al., 1990), LR0R/LR7 (Vilgalys and Hester, 1990), NS1/NS4 (White et al., 1990). A dataset consisting of nrSSU, nrITS, and nrLSU, was assembled on the basis of BLASTn results and of the most recent phylogenetic study focused on Lulworthiales (Azevedo et al., 2017). Reference sequences were retrieved from GenBank ( Table 2). Alignments were generated using MUSCLE (default conditions for gap openings and gap extension penalties), implemented in MEGA v. 7.0 (Molecular Evolutionary Genetics Analysis), visually inspected and trimmed by TrimAl v. 1.2 4 to delimit and discard ambiguously aligned regions. No incongruence was observed among single-loci phylogenetic trees and alignments were concatenated into a single data matrix with SequenceMatrix (Vaidya et al., 2011). The best evolutionary model under the Akaike Information Criterion (AIC) was determined with jModelTest 2 (Darriba et al., 2012). Phylogenetic inference was estimated using Maximum Likehood (ML) and Bayesian Inference (BI) criteria. The ML analysis was produced using RAxML v. 8.1.2 (Stamatakis, 2014) under GTR + I + G evolutionary model and 1,000 bootstrap replicates. BI was performed with MrBayes 3.2.2 (Ronquist et al., 2012) with the same substitution model (GTR + I + G). The alignment was run for 5 million generations with two independent runs each containing four Markov Chains Monte Carlo (MCMC) and sampling every 100 iterations. The first 25% of generated trees were discarded as "burn-in." A consensus tree was generated using the "sumt" function of MrBayes and Bayesian posterior probabilities (BYPP) were calculated. Consensus trees were visualized in FigTree v. 1.4.2 5 . Due to topological similarity of the two resulting trees, only BI analysis with BYPP and MLB values was reported.

Statistical Analysis
Differences among the total fungal loads (CFU g −1 dwt) of each plant portion were inferred by applying the analysis of variance (ANOVA), Bonferroni post hoc test (p < 0.05), using GraphPad Prism 5 for Windows, GraphPad Software (San Diego, CA, United States) 6 .
Significant differences among mycobiota of individual plant parts and organisms were evaluated by applying the PERmutational Multivariate ANalysis Of Variance (PERMANOVA; pseudo-F index; p < 0.05) and visualized by the Canonical Analysis of Principal Coordinates (CAP) or Principal Coordinate Analysis (PCO). The contribution of single species (in percentage) to the diversity observed within and between groups was assessed by SIMilarity PERcentage (SIMPER) analysis. These analyses were performed with the statistical package PRIMER 7 (Plymouth Routines in Multivariate Ecological Research, Albany Auckland, New Zealand).

RESULTS
All plants of P. oceanica were colonized by fungi. The average fungal abundance (CFU g −1 dw) of the four parts ranged from 5.0 × 10 2 CFU g −1 dw to 1.4 × 10 4 CFU g −1 dw (Figure 2). In Ghiaie, rhizomes displayed the highest fungal load, followed by matte, leaves and roots, whereas in Margidore leaves were the richest part, followed by rhizomes and matte evenly, and finally FIGURE 2 | Fungal total load (CFU g -1 dwt) detected in the four parts of P. oceanica (leaves, roots, rhizomes, and matte) in Ghiaie and Margidore. Results are expressed as mean ± SD and analyzed through two-way analysis of variance (ANOVA), Bonferroni post hoc (p < 0.05). Capital letters indicate significant differences within the same sampling site; lowercase letters indicate significant differences between sampling sites. In brackets total number of taxa retrieved.
by roots. Considering the individual plant parts, only roots did not differ significantly between the two sampling sites.
Overall, 287 isolates, representative of 102 taxa were retrieved: 43 from leaves, 70 from matte, 118 from rhizomes, and 56 from roots. The majority of taxa needed a specific medium and incubation temperature: 50 taxa were exclusively isolated from CMASW, 31 from AP and only 21 from both media; 76 taxa grew exclusively at 25 • C and 11 at 15 • C. In comparison to the plants sampled at Ghiaie, the number of taxa retrieved from the plants harvested at Margidore was almost double ( Table 1).
Sixty-eight percent of the strains were identified at species level and 12% at genus level. Those fungi that remained identified at higher rank were sterile in axenic culture and did not accomplished a univocal molecular result. Most of the taxa belonged to Ascomycota (91) while Basidiomycota (8) and Mucoromycota (3) were scarcely represented ( Table 3).
Penicillium and Cladosporium were the most abundant genera with 17 and 9 species, and were isolated from all sites and plant parts, regardless of the culture conditions used.
Out of the 102 taxa isolated, only Cladosporium cladosporioides and Penicillium commune were detected in all four plant parts; 7 and 18 taxa were isolated from three and two parts, respectively. The rest was exclusive to rhizomes (32), matte (15), leaves (14), and roots (13) (Figure 3).
The combined dataset had an aligned length of 2135 characters, of which 1096 were constant, 352 were parsimonyuninformative and 687 parsimony informative.
The strains under investigation formed a strongly supported and separated clade (BYPP = 1.00; MLB = 100%) within Lulworthiaceae, unique family of the order Lulworthiales ( Figure 5). Since this clade was close to the genus Lulworthia sensu latu, we introduced the novel genus Paralulworthia that includeed the two species Paralulworthia gigaspora and Paralulworthia posidoniae. The development of the sexual morphotype of MUT 435 and MUT 5261 on maritime pine wood confirmed the divergence between the two species.
Description: growing efficiently on wood. Ascomata subglobose to ellipsoidal, single, scattered, immersed in the cortex or erumpent, brown to black at maturity; peridium dark brown composed of several layers of thick-walled cells of textura angularis. Asci containing 8 spores, dehiscent. Sterile elements not observed. Ascospores uni-to 3-septate at maturity, constricted at septa (particularly the central one), ellipsoid to fusiform, pointed at both ends, thin-walled, from light to dark brown at maturity, straight.
Etymology: in reference to the large size of the ascospores.
Asexual morph not observed. Colony description: growing on MEASW from 60 to 90 mm Ø after 28 days at 21 • C, mycelium from feltrose to floccose, from whitish to ochre, reverse from light brown to bordeaux. A caramel colored diffusible pigment and exudates were often present.
Asexual morph not observed. Colony description: growing on MEASW, reaching 42 mm Ø after 28 days at 21 • C, mycelium feltrose with irregular edges, whitish with slight shades of purple in the middle, gray-green at the edges; reverse caramel. A caramel colored diffusible pigment and exudates were often present.

DISCUSSION
Broadening the knowledge of marine fungi associated with a vast range of biotic and abiotic substrates is becoming more and more important. The reason behind this lies in the ecological relevance that these organisms may play in the oceans and in their possible exploitation as a source of potentially novel enzymes and bioactive compounds (Bovio et al., 2019;Grossart et al., 2019;Overy et al., 2019;Youssef et al., 2019).
Nowadays, P. oceanica is rapidly disappearing from the Mediterranean Sea, being seriously threatened by a number of factors 7 . Its mycobiota was partially studied in the past (Cuomo et al., 1985;Panno et al., 2013). Recently Vohnik and collaborators described the endophytes colonizing the roots of this endemic seagrass (Vohnik et al., 2017(Vohnik et al., , 2019. Here, we investigated the cultivable mycobiota of P. oceanica with the intent of comparing it with the fungal communities of two photobionts of the same phytocoenosis.

Mycobiota of P. oceanica
The observed total fungal load ranged from 5 × 10 2 to 1.4 × 10 4 CFU g −1 dw, approximately one order of magnitude higher than what reported by Panno et al. (2013) (up to 1.6 × 10 3 CFU g −1 dw). Since the two sampling occurred in the same period of the year (March), we hypothesize that the climate factor did not affect the total load. The fungal abundance through the plant parts was also slightly different: in the present study, the fungal load followed the distribution rhizomes > leaves > matte > roots, whereas in the previous investigation matte and leaves were swapped (Panno et al., 2013). The use of specific culture media and of two incubation temperatures maximized the number of taxa isolated. Indeed, a lower temperature (15 • C) allowed the growth of psychrotolerant/psychrophilic fungi such as four species of Penicillium (P. canescens, P. glabrum, P. paneum, P. vinaceum) ( Table 2). Most probably, this is one of the reasons why we retrieved a number of taxa higher than what reported in 2013 by Panno (102 vs. 88).
As expected, Ascomycota prevailed, confirming once again the dominance of this phylum in the oceans 8 . In detail, Penicillium and Cladosporium, were widely represented in all plant parts and sampling sites, and were characterized by the highest number of species (17 and 9, respectively). Members of Penicillium and Cladosporium are commonly reported in marine environments worldwide, probably due to the high adaptation to the peculiar chemical-physical conditions (Raghukumar, 2017). On the other side, it is not unreasonable to relate the high abundance of these genera to their extraordinary sporulation rate that may end in an overestimation of the substrate colonization (Imhoff et al., 2011).
Interestingly, 37 species were retrieved for the first time from the Mediterranean Sea and 61 from P. oceanica (Table 2), while only 19 species were previously detected in P. oceanica (Cuomo et al., 1985;Panno et al., 2013;Gnavi et al., 2014). Several species, besides being isolated for the first time from the Mediterranean Sea and from P. oceanica, were never found in any marine habitat before. For instance, this is the case of Alfaria dandenongensis, Aspergillus pseudoglaucus, Cerrena unicolor, Didymella fabae, and Neoroussoella solani. A. pseudoglaucus is an opportunistic human pathogen (Siqueira et al., 2018) as well as N. solani, a causal agent of keratomycosis (Mochizuki et al., 2017), while 7 www.eea.europa.eu/data-and-maps/figures/distribution-of-the-marineangiosperm-posidonia-oceanica-and-zostera-sp-in-the-mediterranean 8 www.marinefungi.org their underwater ecological role is still obscure. On the contrary, we can hypothesize the role of the plant pathogen D. fabae and of the leaf litter saprobe A. dandenongensis 9 . Similarly, it is possible to assume an active role for white rot fungi (i.e., Bjerkandera adusta, C. unicolor and Schizophyllum commune) in transforming recalcitrant matrixes . Finally, entomopathogenic fungi like Cordyceps farinosa, due to their ability to degrade chitin, may be capable of attacking the shellfish carapace.

Distribution Among Plant Parts
The mycobiota of P. oceanica changed significantly among the four parts of the plant, with the exception of leaves vs. matte. This last observation could be ascribable to the similar structure of the substrates: the main components of matte are the decaying and fluctuating leaves (Vohnik et al., 2017). The presence of antagonistic organisms, both micro-and macro-, together with the specific localization of toxic metabolites, such as tannic acid in the leaves (Pergent et al., 2008), may affect the distribution among plant parts. It is not unexpected to retrieve unique mycobiota on differently organized and composed parts of a complex plant, as occurs in terrestrial phanerogams.
The only two species common to the four regions were C. cladosporioides and P. commune: the former is a usual inhabitant of the sea (Garzoli et al., 2015), the latter, although observed for the first time on P. oceanica, was already found associated with F. petiolata . P. chrysogenum dominated the mycobiota of the leaves, whereas the marine genus Corollospora prevailed in matte. With respect to rhizomes and roots, two taxa (Lulworthiales sp. and P. antarcticum) were the main components of the fungal communities, even though their proportions were different in the two plant parts. Vohník et al. (2016) found members of the order Lulworthiales in the rhizomes and roots of P. oceanica. Lulworthiales degrade wood and marsh plants in marine and estuarine environments and seem to be the most common seagrasses decomposers followed by Corollospora maritima. Indeed, C. maritima and Lulworthiales spp., are wellknown cellulases producers and can break down complex lignocellulose compounds, making them available for other organisms nourishment (Raghukumar, 2017). Interestingly, two new species of Lulworthiales were found in rhizomes and described here for the first time.
The retrieval of P. chrysogenum in leaves and, to a lower extent in matte, may be related to defense mechanisms against predators and/or pathogens, as demonstrated in terrestrial plants, where the presence of this fungus stimulates the synthesis of secondary metabolites such as salycilic and jasmonic acids (Chen et al., 2018).
The Mycobiota of P. oceanica Compared With the Fungal Communities of F. petiolata and P. pavonica F. petiolata and P. pavonica are two algae that characterize the Neptune forests. Since their mycobiota were recently described Garzoli et al., 2018), we considered worthwhile to compare the fungal communities associated with these three organisms that grow only few centimeters apart. To this end, one must bear in mind that samplings occurred at the same time and in the same sites (Ghiaie and Margidore): it would be interesting to clarify whether a strict specificity between fungi and photobionts exists.
In terms of number of taxa, P. pavonica displayed the richest mycoflora, followed by P. oceanica and finally by F. petiolata, where the number of taxa was almost halved. As shown in Figure 8, P. oceanica and the two algae had 11 species in common: Acrostalagmus luteoalbus, Arthrinium arundinis, Cladosporium allicinum, C. cladosporioides, C. sphaerospermum, Gibellulopsis nigrescens, Massarina rubi, P. antarcticum, P. brevicompactum, P. roseopurpureum, and S. commune. Besides their broad diffusion in the oceans, these fungi, due to the secretion of secondary metabolites (e.g., antialgal, antifungal, and/or antibacterial) may help in deterring colonization of their hosts by other microbes, and in warding off predators (Raghukumar, 2017). On the other side, the three substrates show specific mycobiota, significantly different from each other (Figure 9), with high and equivalent percentages of dissimilarity (94.68-95.74%). The uniqueness of each community derives from the prevalence of exclusive taxa (Figure 8).
Confirming the "substrate specificity" is also the fact that there is a significant difference between the mycobiota of the three organisms even if the "non-photosynthesizing" parts of P. oceanica were excluded from the analysis. The leaves of P. oceanica are in fact characterized by a mycobiota significantly different from that of the other two photobionts (data not shown).  This divergence in "substrate specificity" may be due to the different composition and morphology of the hosts (e.g., P. pavonica with its hairy surface could entrap a high number of fungal propagules, contrary to the slippery texture of F. petiolata) and/or to the different metabolites produced by seagrasses and algae. Extracts of P. oceanica were proved active against bacteria, dermatophytes, and yeasts (Bernard and Pesando, 1989) as well as against Trypanosoma brucei rhodesiense and Leishmania donovani (Orhan et al., 2006). Likewise, the activity of P. pavonica extracts against bacteria and fungi was comparable to the efficiency showed by antibiotics and antimycotics used as control (Al-Enazi et al., 2018). Antibacterial, antiviral, antimitotic and antifungal properties were also detected in raw extracts of F. petiolata (Ballesteros et al., 1992). However, whether the macro-organisms, their microbiota, or the holobiont derived from their strict interaction is responsible for the production of the active compounds is still unclear. With the exclusion of cosmopolitan and/or widespread fungi, we can hypothesize the shaping of a relation between the microorganisms and their host, as recently demonstrated for three Atlantic sponges (Bovio et al., 2018). The establishment of a holobiont depends on both active (selection) and passive (drift) events. Due to the high connectivity of aquatic environments, differences in microbiota may be dependent on a combination of selection and drift. In most cases, fungi spread by passive dispersion of hyphal fragments and/or propagules that happen to adhere to the host. Their further development is consequent to the compatibility and interdependence between the individual organisms. In other words, the holobiont is an organism with its own features, modeled upon a mutual selection. The unimportance of the sampling site supports this idea. Ghiaie and Margidore are two extremely different sites both geographically (Ligurian Sea and Tyrrhenian Sea, respectively) and for the anthropic influence: while the former is located in a marine protected area, the latter is subjected to an intense anthropic disturbance. Nevertheless, no significant difference was detected in terms of fungal biodiversity.
An often-underappreciated ecological function of holobionts is their ability to restrain and nourish rare species: the host provides a niche favorable for the growth of specific communities different from those of the surrounding environment (Dittami et al., 2019). To this regard, we may define "rare" the two novel species here introduced.

Paralulworthia gigaspora and P. posidoniae
This investigation introduces two new species that belong to a novel genus of Lulworthiaceae, the sole family of the order Lulworthiales (Sordariomycetes). Lulworthiaceae spp. are strictly marine species found on a range of substrates (e.g., submerged wood, seaweed, seafoam, etc.). This group includes the following genera: Cumulospora, Halazoon, Hydea, Kohlmeyerella, Lulwoana, Lulworthia, Lindra, Matsusporium, and Moleospora. The new genus Paralulworthia accommodates the four strains MUT 435, MUT 672, MUT 5413, and MUT 5261. P. gigaspora and P. posidoniae formed a strongly supported cluster segregated from the genus Lulworthia sensu latu. The previously suggested polyphyletic nature of Lulworthia spp. (Azevedo et al., 2017), that discriminate Lulworthia sensu latu from Lulworthia sensu strictu is here confirmed and justifies the introduction of the new genus Paralulworthia. Besides, members of the genus Lulworthia are morphologically characterized by filamentous hyaline ascospores (Kohlmeyer et al., 2000), well distinct from the ellipsoidal brown ascospores of the new genus. Furthermore, beyond the phylogenetic distance, the dimension of the ascomata and ascospores allowed the recognition of two species. P. gigaspora and P. posidoniae were isolated exclusively from rhizomes that are charachterized by a high content of lignin (Kaal et al., 2016). Reproductive structures were detected only following colonization of P. pinaster wood specimens submerged in sea-water. These two observations indicate a possible lignicolous nature of these species.

CONCLUSION
The seagrass P. oceanica is home to a huge variety of fungi. The comparison of our data with those previously obtained on two algae collected in the same phytocoenosis underlines a substrate-specificity. Besides, we demonstrated a plant-part-specificity of the fungal communities associated with the host. Marine phanerogames, likewise their terrestrial counterparts, enroll distinct microbiomes in different regions (i.e., phylloplane, rhizoplane etc.). This selective distribution seems to increment the resistance of plant to pathogens and predators by producing chemo-attractants and anti-microbial compounds. This is even more important for Posidonia meadows that are seriously threatened by climate changes, allochthonous species and anthropic activities. Nonetheless, further studies will need to clarify the mechanisms (e.g., molecular, physiological, etc.) at the basis of the specific interaction between microorganisms and hosts responsible for the shaping of a holobiont.

DATA AVAILABILITY STATEMENT
The sequences generated in this study can be found in GenBank. The accession numbers of the sequences deposited in

AUTHOR CONTRIBUTIONS
AP conducted molecular studies, phylogenesis, data analysis, and manuscript writing. EB and LR carried morphological analysis. VP was responsible for morphological analysis and manuscript writing. GV supervised the work.

FUNDING
This study was funded by University of Torino (ex 60%) and Fondazione CRT (Turin, Italy).