Xylella fastidiosa Infection Reshapes Microbial Composition and Network Associations in the Xylem of Almond Trees

Xylella fastidiosa represents a major threat to important crops worldwide including almond, citrus, grapevine, and olives. Nowadays, there are no efficient control measures for X. fastidiosa , and the use of preventive measures and host resistance represent the most practical disease management strategies. Research on vessel-associated microorganisms is gaining special interest as an innate natural defense of plants to cope against infection by xylem-inhabiting pathogens. The objective of this research has been to characterize, by next-generation sequencing (NGS) analysis, the microbial communities residing in the xylem sap of almond trees affected by almond leaf scorch disease (ALSD) in a recent X. fastidiosa outbreak occurring in Alicante province, Spain. We also determined community composition changes and network associations occurring between xylem-inhabiting microbial communities and X. fastidiosa . For that, a total of 91 trees with or without ALSD symptoms were selected from a total of eight representative orchards located in five municipalities within the X. fastidiosa -demarcated area. X. fastidiosa infection in each tree was verified by quantitative polymerase chain reaction (qPCR) analysis, with 54% of the trees being tested X. fastidiosa -positive. Globally, Xylella (27.4%), Sphingomonas (13.9%), and Hymenobacter (12.7%) were the most abundant bacterial genera, whereas Diplodia (30.18%), a member of the family Didymellaceae (10.7%), and Aureobasidium (9.9%) were the most predominant fungal taxa. Furthermore, principal coordinate analysis (PCoA) of Bray–Curtis and weighted UniFrac distances differentiated almond xylem bacterial communities mainly according to X. fastidiosa infection, in contrast to fungal community structure that was not closely related to the presence of the pathogen. Similar results were obtained when X. fastidiosa reads were removed from the bacterial data set although the effect was less pronounced. Co-occurrence network analysis revealed negative associations among four amplicon sequence variants (ASVs) assigned to X. fastidiosa with different bacterial ASVs belonging to 1174-901-12 , Abditibacterium , Sphingomonas , Methylobacterium–Methylorubrum , Modestobacter , Xylophilus , and a non-identified member of the family Solirubrobacteraceae. Determination of the close-fitting associations between xylem-inhabiting microorganisms and X. fastidiosa may help to reveal specific microbial players associated with the suppression o