A comprehensive, large-scale analysis of "terroir" cheese and milk microbiota reveals profiles strongly shaped by both geographical and human factors.

An exhaustive analysis was performed on more than 2000 microbiotas from French Protected Designation of Origin (PDO) cheeses, covering most cheese families produced throughout the world. Thanks to a complete and accurate set of associated metadata, we have carried out a deep analysis of the ecological drivers of microbial communities in milk and "terroir" cheeses. We show that bacterial and fungal microbiota from milk differed significantly across dairy species while sharing a core microbiome consisting of four microbial species. By contrast, no microbial species were detected in all ripened cheese samples. Our network analysis suggested that the cheese microbiota was organized into independent network modules. These network modules comprised mainly species with an overall relative abundance lower than 1%, showing that the most abundant species were not those with the most interactions. Species assemblages differed depending on human drivers, dairy species, and geographical area, thus demonstrating the contribution of regional know-how to shaping the cheese microbiota. Finally, an extensive analysis at the milk-to-cheese batch level showed that a high proportion of cheese taxa were derived from milk under the influence of the dairy species and protected designation of origin.

Serial fermentation in milk generates functionally diverse community lineages with different degrees of structure stabilization.

Microbial communities offer considerable potential for tackling environmental challenges by improving the functioning of ecosystems. Top-down community engineering is a promising strategy that could be used to obtain communities of desired function. However, the ecological factors that control the balance between community shaping and propagation are not well understood. Dairy backslopping, which consists of using part of the previous production to inoculate a new one, can be used as a model engineering approach to investigate community dynamics during serial propagations. In this study, 26 raw milk samples were serially propagated 6 times each, giving rise to 26 community lineages. Bacterial community structures were analyzed by metabarcoding, and acidification was recorded by pH monitoring. The results revealed that different types of community lineages could be obtained in terms of taxonomic composition and dynamics. Five lineages reached a repeatable community structure in a few propagation steps, with little variation between the final generations, giving rise to stable acidification kinetics. Moreover, these stabilized communities presented a high variability of structure and diverse acidification properties between community lineages. Besides, the other lineages were characterized by different levels of dynamics leading to parallel or divergent trajectories. The functional properties and dynamics of the communities were mainly related to the relative abundance and the taxonomic composition of lactic acid bacteria within the communities. These findings highlight that short-term schemes of serial fermentation can produce communities with a wide range of dynamics and that the balance between community shaping and propagation is intimately linked to community structure. IMPORTANCE: Microbiome applications require approaches for shaping and propagating microbial communities. Shaping allows the selection of communities with desired taxonomic and functional properties, while propagation allows the production of the biomass required to inoculate the engineered communities in the target ecosystem. In top-down community engineering, where communities are obtained from a pool of mixed microorganisms by acting on environmental variables, a major challenge is to master the balance between shaping and propagation. However, the ecological factors that favor high dynamics of community structure and, conversely, those that favor stability during propagation are not well understood. In this work, short-term dairy backslopping was used to investigate the key role of the taxonomic composition and structure of bacterial communities on their dynamics. The results obtained open up interesting prospects for the biotechnological use of microbiomes, particularly in the field of dairy fermentation, to diversify approaches for injecting microbial biodiversity into cheesemaking processes.

Genomic Characterization of Wild Lactobacillus delbrueckii Strains Reveals Low Diversity but Strong Typicity.

Investigating the diversity of a given species could give clues for the development of autochthonous starter cultures. However, few studies have focused on the intraspecies diversity of Lactobacillus delbrueckii strains, a technologically important lactic acid bacterium for the dairy industry. For this reason, Lactobacillus delbrueckii strains from the Saint-Nectaire Protected Designation of Origin (PDO) area were isolated and characterized. Genetic diversity was determined based on core genome phylogenetic reconstruction and pangenome analysis, while phenotypic assessments encompassed proteolysis and volatile compound production potential. A total of 15 L. delbrueckii ssp. lactis unique new strains were obtained. The genetic analysis and further proteolytic activities measurement revealed low variability among these Saint-Nectaire strains, while substantial genetic variability was observed within the L. delbrueckii ssp. lactis subspecies as a whole. The volatile compound profiles slightly differed among strains, and some strains produced volatile compounds that could be of particular interest for cheese flavor development. While the genetic diversity among Saint-Nectaire strains was relatively modest compared to overall subspecies diversity, their distinct characteristics and pronounced differentiation from publicly available genomes position them as promising candidates for developing autochthonous starter cultures for cheese production.

Molecular and Biological Characterization of Novel and Known Family Secoviridae Members Infecting Lettuce.

High-throughput sequencing of two lettuces showing virus-like symptoms in France provided evidence of infection by members of the family Secoviridae. One plant (JG1) had a complex mixed infection that involved, among others, a novel waikavirus (lettuce waikavirus 1) and two isolates of a sequivirus related to lettuce mottle virus (LeMoV). The second lettuce plant (JG2) was singly infected by LeMoV. Complete genomic sequences were obtained for all four isolates and, in addition, near complete genome sequences were obtained for other LeMoV or LeMoV-related isolates (from French cultivated and wild lettuces and from a Brazilian cultivated lettuce) and for two isolates of another family Asteraceae-infecting sequivirus, dandelion yellow mosaic virus (DaYMV). Analysis of these genomic sequences allows the proposal of tentative genome organization for the various viruses and clarification of their phylogenetic relationships. Sequence and host range comparisons point to significant differences between the two sequivirus isolates identified in the JG1 plant and LeMoV isolates from France and Brazil, suggesting they belong to a novel species for which the name lettuce star mosaic virus is proposed.

Dietary Live Yeast Supplementation Influence on Cow's Milk, Teat and Bedding Microbiota in a Grass-Diet Dairy System.

The supplementation of animal feed with microbial additives remains questioning for the traditional or quality label raw milk cheeses with regard to microbial transfer to milk. We evaluated the effect of dietary administration of live yeast on performance and microbiota of raw milk, teat skin, and bedding material of dairy cows. Two balanced groups of cows (21 primiparous 114 ± 24 DIM, 18 multiparous 115 ± 33 DIM) received either a concentrate supplemented with Saccharomyces cerevisiae CNCM I-1077 (1 × 1010 CFU/d) during four months (LY group) or no live yeast (C group). The microbiota in individual milk samples, teat skins, and bedding material were analysed using culture dependent techniques and high-throughput amplicon sequencing. The live yeast supplementation showed a numerical increase on body weight over the experiment and there was a tendency for higher milk yield for LY group. A sequence with 100% identity to that of the live yeast was sporadically found in fungal amplicon datasets of teat skin and bedding material but never detected in milk samples. The bedding material and teat skin from LY group presented a higher abundance of Pichia kudriavzevii reaching 53% (p < 0.05) and 10% (p < 0.05) respectively. A significant proportion of bacterial and fungal ASVs shared between the teat skin and the milk of the corresponding individual was highlighted.

Genetic and technological diversity of Streptococcus thermophilus isolated from the Saint-Nectaire PDO cheese-producing area.

Streptococcus thermophilus is of major importance for cheese manufacturing to ensure rapid acidification; however, studies indicate that intensive use of commercial strains leads to the loss of typical characteristics of the products. To strengthen the link between the product and its geographical area and improve the sensory qualities of cheeses, cheese-producing protected designations of origin (PDO) are increasingly interested in the development of specific autochthonous starter cultures. The present study is therefore investigating the genetic and functional diversity of S. thermophilus strains isolated from a local cheese-producing PDO area. Putative S. thermophilus isolates were isolated and identified from milk collected in the Saint-Nectaire cheese-producing PDO area and from commercial starters. Whole genomes of isolates were sequenced, and a comparative analysis based on their pan-genome was carried out. Important functional properties were studied, including acidifying and proteolytic activities. Twenty-two isolates representative of the diversity of the geographical area and four commercial strains were selected for comparison. The resulting phylogenetic trees do not correspond to the geographical distribution of isolates. The clustering based on the pan-genome analysis indicates that isolates are divided into five distinct groups. A Kyoto Encyclopedia of Genes and Genomes (KEGG) functional annotation of the accessory genes indicates that the accessory gene contents of isolates are involved in different functional categories. High variability in acidifying activities and less diversity in proteolytic activities were also observed. These results indicate that high genetic and functional variabilities of the species S. thermophilus may arise from a small (1,800 km2) geographical area and may be exploited to meet demand for use as autochthonous starters.

A new potyvirus from hedge mustard (Sisymbrium officinale (L.) Scop.) sheds light on the evolutionary history of turnip mosaic virus.

A novel potyvirus was identified in symptomatic hedge mustard (Sisymbrium officinale (L.) Scop.) and wild radish (Raphanus raphanistrum L.) in France. The nearly complete genome sequence of hedge mustard mosaic virus (HMMV) was determined, demonstrating that it belongs to a sister species to turnip mosaic virus (TuMV). HMMV readily infected several other members of the family Brassicaceae, including turnip, shepherd's purse (Capsella bursa-pastoris), and arabidopsis. The identification of HMMV as a Brassicaceae-infecting virus closely related to TuMV leads us to question the current scenario of TuMV evolution and suggests a possible alternative one in which transition from a monocot-adapted ancestral lifestyle to a Brassicaceae-adapted one could have occurred earlier than previously recognized.Please check and confirm that the authors and their respective affiliations have been correctly identified and amend if necessary.all OK.

Role of Klebsiella pneumoniae Type VI secretion system (T6SS) in long-term gastrointestinal colonization.

Type VI secretion systems (T6SS), recently described in hypervirulent K. pneumoniae (hvKp) strains, are involved in bacterial warfare but their role in classical clinical strains (cKp) has been little investigated. In silico analysis indicated the presence of T6SS clusters (from zero to four), irrespective of the strains origin or virulence, with a high prevalence in the K. pneumoniae species (98%). In the strain CH1157, two T6SS-apparented pathogenicity islands were detected, T6SS-1 and -2, harboring a phospholipase-encoding gene (tle1) and a potential new effector-encoding gene named tke (Type VI Klebsiella effector). Tle1 expression in Escherichia coli periplasm affected cell membrane permeability. T6SS-1 isogenic mutants colonized the highest gastrointestinal tract of mice less efficiently than their parental strain, at long term. Comparative analysis of faecal 16S sequences indicated that T6SS-1 impaired the microbiota richness and its resilience capacity. Oscillospiraceae family members could be specific competitors for the long-term gut establishment of K. pneumoniae.

Functional changes in Bleu d'Auvergne cheese during ripening.

The authentic characteristics of the famous Bleu d'Auvergne cheese were studied. Many parameters were analysed during the ripening of cheeses. Migrations of Na and Ca ions, associated with a pH gradient, occurred between the rind and the core. At 34 days, this cheese had a high salt content (2.87 %), contributing to 23 % of the recommended sodium intake for adults, but significant calcium (6.14 g/kg) and vitamin B12 (1.14 µg/100 g) levels. Thus, a 40 g serving contributed to 25 % of the population reference intake for Ca and 11 % of the adequate intake for B12. Proteolysis, yeast and mould counts strongly increased. Lactococcus and Streptococcus were predominant and correlated with B2 and B6 levels. Bleu d'Auvergne was characterised by salty taste, blue odour and aroma. This cheese has a noticeable B vitamins concentration, but the level of salt should be reduced to meet the nutritional guidelines, possibly by implementing alternative salting methods.

Molecular characterization of Cordyline virus 1 isolates infecting yam (Dioscorea spp).

Cordyline virus 1 (CoV1) is a velarivirus that has so far only been reported in ornamental Ti plants (Cordyline fruticosa). Using high-throughput sequencing, we identified CoV1 infection in yam accessions from Vanuatu. Using a specific RT-PCR assay, we found that CoV1 is also present and highly prevalent in Dioscorea alata, D. cayenensis, and D. trifida in Guadeloupe. Phylogenetic analysis showed that CoV1 isolates infecting yam in Guadeloupe display a low level of molecular diversity. These data provide insights into the transmission of CoV1 in yam in Guadeloupe.

Integration of Multiomic Data to Characterize the Influence of Milk Fat Composition on Cantal-Type Cheese Microbiota.

A previous study identified differences in rind aspects between Cantal-type cheeses manufactured from the same skimmed milk, supplemented with cream derived either from pasture-raised cows (P) or from cows fed with maize silage (M). Using an integrated analysis of multiomic data, the present study aimed at investigating potential correlations between cream origin and metagenomic, lipidomic and volatolomic profiles of these Cantal cheeses. Fungal and bacterial communities of cheese cores and rinds were characterized using DNA metabarcoding at different ripening times. Lipidome and volatolome were obtained from the previous study at the end of ripening. Rind microbial communities, especially fungal communities, were influenced by cream origin. Among bacteria, Brachybacterium were more abundant in P-derived cheeses than in M-derived cheeses after 90 and 150 days of ripening. Sporendonema casei, a yeast added as a ripening starter during Cantal manufacture, which contributes to rind typical aspect, had a lower relative abundance in P-derived cheeses after 150 days of ripening. Relative abundance of this fungus was highly negatively correlated with concentrations of C18 polyunsaturated fatty acids and to concentrations of particular volatile organic compounds, including 1-pentanol and 3-methyl-2-pentanol. Overall, these results evidenced original interactions between milk fat composition and the development of fungal communities in cheeses.

Host range and molecular variability of the sadwavirus dioscorea mosaic associated virus.

Dioscorea mosaic associated virus (DMaV) is a member of the genus Sadwavirus, family Secoviridae, that is associated with mosaic symptoms in Dioscorea rotundata in Brazil. The genome of a DMaV isolate detected in D. trifida in Guadeloupe was sequenced by high-throughput sequencing. Using an RT-PCR-based detection assay, we found that DMaV infects D. alata, D. bulbifera, D. cayenensis-rotundata, D. esculenta, and D. trifida accessions conserved in Guadeloupe and Côte d'Ivoire and displays a very high level of molecular diversity in a relatively small region of the genome targeted by the assay. We also provide evidence that DMaV is also present in D. rotundata in Benin and in D. alata in Nigeria.

In vivo investigation of Lcr35® anti-candidiasis properties in Caenorhabditis elegans reveals the involvement of highly conserved immune pathways.

Lactic acid bacteria, including the microorganisms formerly designated as Lactobacillus, are the major representatives of Live Biotherapeutic Microorganisms (LBM) when used for therapeutic purposes. However, in most cases, the mechanisms of action remain unknown. The antifungal potential of LBM has already been demonstrated using preclinical models (cell cultures, laboratory animals). Understanding their mechanisms of action is strategic for the development of new therapeutics for humans. Here, Caenorhabditis elegans was used as an in vivo model to analyze pro-longevity, anti-aging and anti-candidiasis effects of the LBM Lacticaseibacillus rhamnosus (formerly Lactobacillus rhamnosus) Lcr35®. A high-throughput transcriptomic analysis revealed a specific response of C. elegans depending on whether it is in the presence of the LBM L. rhamnosus Lcr35® (structural response), the yeast Candida albicans (metabolic response) or both (structural and metabolic responses) in a preventive and a curative conditions. Studies on C. elegans mutants demonstrated that the p38 MAPK (sek-1, skn-1) and the insulin-like (daf-2, daf-16) signaling pathways were involved in the extended lifespan provided by L. rhamnosus Lcr35® strain whereas the JNK pathway was not involved (jnk-1). In addition, the anti C. albicans effect of the bacterium requires the daf-16 and sek-1 genes while it is independent of daf-2 and skn-1. Moreover, the anti-aging effect of Lcr35®, linked to the extension of longevity, is not due to protection against oxidative stress (H2O2). Taken together, these results formally show the involvement of the p38 MAP kinase and insulin-like signaling pathways for the longevity extension and anti-Candida albicans properties of Lcr35® with, however, differences in the genes involved. Overall, these findings provide new insight for understanding the mechanisms of action of a probiotic strain with antimicrobial potential.

Lactic Starter Dose Shapes S. aureus and STEC O26:H11 Growth, and Bacterial Community Patterns in Raw Milk Uncooked Pressed Cheeses.

Adding massive amounts of lactic starters to raw milk to manage the sanitary risk in the cheese-making process could be detrimental to microbial diversity. Adjusting the amount of the lactic starter used could be a key to manage these adverse impacts. In uncooked pressed cheeses, we investigated the impacts of varying the doses of a lactic starter (the recommended one, 1×, a 0.1× lower and a 2× higher) on acidification, growth of Staphylococcus aureus SA15 and Shiga-toxin-producing Escherichia coli (STEC) O26:H11 F43368, as well as on the bacterial community patterns. We observed a delayed acidification and an increase in the levels of pathogens with the 0.1× dose. This dose was associated with increased richness and evenness of cheese bacterial community and higher relative abundance of potential opportunistic bacteria or desirable species involved in cheese production. No effect of the increased lactic starter dose was observed. Given that sanitary criteria were paramount to our study, the increase in the pathogen levels observed at the 0.1× dose justified proscribing such a reduction in the tested cheese-making process. Despite this, the effects of adjusting the lactic starter dose on the balance of microbial populations of potential interest for cheese production deserve an in-depth evaluation.

rANOMALY: AmplicoN wOrkflow for Microbial community AnaLYsis.

Bioinformatic tools for marker gene sequencing data analysis are continuously and rapidly evolving, thus integrating most recent techniques and tools is challenging. We present an R package for data analysis of 16S and ITS amplicons based sequencing. This workflow is based on several R functions and performs automatic treatments from fastq sequence files to diversity and differential analysis with statistical validation. The main purpose of this package is to automate bioinformatic analysis, ensure reproducibility between projects, and to be flexible enough to quickly integrate new bioinformatic tools or statistical methods. rANOMALY is an easy to install and customizable R package, that uses amplicon sequence variants (ASV) level for microbial community characterization. It integrates all assets of the latest bioinformatics methods, such as better sequence tracking, decontamination from control samples, use of multiple reference databases for taxonomic annotation, all main ecological analysis for which we propose advanced statistical tests, and a cross-validated differential analysis by four different methods. Our package produces ready to publish figures, and all of its outputs are made to be integrated in Rmarkdown code to produce automated reports.

Leaf-associated fungal and viral communities of wild plant populations differ between cultivated and natural ecosystems.

Plants are colonized by diverse fungal and viral communities that influence their growth and survival as well as ecosystem functioning. Viruses interact with both plants and the fungi they host. Our understanding of plant-fungi-virus interactions is very limited, especially in wild plants. Combining metagenomic and culturomic approaches, we assessed the richness, diversity, and composition of leaf-associated fungal and viral communities from pools of herbaceous wild plants representative of four sites corresponding to cultivated or natural ecosystems. We identified 161 fungal families and 18 viral families comprising 249 RNA-dependent RNA polymerase-based operational taxonomic units (RdRp OTUs) from leaves. Fungal culturomics captured 12.3% of the fungal diversity recovered with metagenomic approaches and, unexpectedly, retrieved viral OTUs that were almost entirely different from those recovered by leaf metagenomics. Ecosystem management had a significant influence on both leaf mycobiome and virome, with a higher fungal community richness in natural ecosystems and a higher viral family richness in cultivated ecosystems, suggesting that leaf-associated fungal and viral communities are under the influence of different ecological drivers. Both the leaf-associated fungal and viral community compositions showed a strong site-specificity. Further research is needed to confirm these trends and unravel the factors structuring plant-fungi-virus interactions in wild plant populations.

Yam asymptomatic virus 1, a novel virus infecting yams (Dioscorea spp.) with significant prevalence in a germplasm collection.

A novel virus infecting yams (Dioscorea spp.), tentatively named "yam asymptomatic virus 1" (YaV1), was characterized and sequenced from an asymptomatic D. alata plant from Vanuatu. Sequence comparisons and phylogenetic analysis showed that YaV1 is a novel ampelovirus and has the smallest genome among "subgroup 1" members. RT-PCR-based screening of a yam germplasm collection conserved in Guadeloupe showed that YaV1 is prevalent in D. alata, D. bulbifera, D. cayennensis subsp. rotundata, D. esculenta and D. trifida accessions but causes no apparent symptoms. Additional phylogenetic analysis revealed a low variability of YaV1 in Guadeloupe in a limited part of the genome, and suggested the occurrence of plant-to-plant transmission.

Phytovirome Analysis of Wild Plant Populations: Comparison of Double-Stranded RNA and Virion-Associated Nucleic Acid Metagenomic Approaches.

Metagenomic studies have indicated that the diversity of plant viruses was until recently far underestimated. As important components of ecosystems, there is a need to explore the diversity and richness of the viruses associated with plant populations and to understand the drivers shaping their diversity in space and time. Two viral sequence enrichment approaches, double-stranded RNA (dsRNA) and virion-associated nucleic acids (VANA), have been used and compared here for the description of the virome of complex plant pools representative of the most prevalent plant species in unmanaged and cultivated ecosystems. A novel bioinformatics strategy was used to assess viral richness not only at the family level but also by determining operational taxonomic units (OTU) following the clustering of conserved viral domains. A large viral diversity dominated by novel dsRNA viruses was detected in all sites, while a large between-site variability limited the ability to draw a clear conclusion on the impact of cultivation. A trend for a higher diversity of dsRNA viruses was nevertheless detected in unmanaged sites (118 versus 77 unique OTUs). The dsRNA-based approach consistently revealed a broader and more comprehensive diversity for RNA viruses than the VANA approach, whatever the assessment criterion. In addition, dissimilarity analyses indicated both approaches to be largely reproducible but not necessarily convergent. These findings illustrate features of phytoviromes in various ecosystems and a novel strategy for precise virus richness estimation. These results allow us to reason methodological choices in phytovirome studies and likely in other virome studies where RNA viruses are the focal taxa.IMPORTANCE There are today significant knowledge gaps on phytovirus populations and on the drivers impacting them but also on the comparative performance-methodological approaches for their study. We used and compared two viral sequence enrichment approaches, double-stranded RNAs (dsRNA) and virion-associated nucleic acids (VANA), for phytovirome description in complex pools representative of the most prevalent plant species in unmanaged and cultivated ecosystems. Viral richness was assessed by determining operational taxonomic units (OTU) following the clustering of conserved viral domains. There is some limited evidence of an impact of cultivation on viral populations. These results provide data allowing us to reason the methodological choices in virome studies. For researchers primarily interested in RNA viruses, the dsRNA approach is recommended because it consistently provided a more comprehensive description of the analyzed phytoviromes, but it understandably underrepresented DNA viruses and bacteriophages.

Characterization of the virome of shallots affected by the shallot mild yellow stripe disease in France.

To elucidate the etiology of a new disease of shallot in France, double-stranded RNAs from asymptomatic and symptomatic shallot plants were analyzed using high-throughput sequencing (HTS). Annotation of contigs, molecular characterization and phylogenetic analyses revealed the presence in symptomatic plants of a virus complex consisting of shallot virus X (ShVX, Allexivirus), shallot latent virus (SLV, Carlavirus) and two novel viruses belonging to the genera Carlavirus and Potyvirus, for which the names of shallot virus S (ShVS) and shallot mild yellow stripe associated virus (SMYSaV), are proposed. Complete or near complete genomic sequences were obtained for all these agents, revealing divergent isolates of ShVX and SLV. Trials to fulfill Koch's postulates were pursued but failed to reproduce the symptoms on inoculated shallots, even though the plants were proved to be infected by the four viruses detected by HTS. Replanting of bulbs from SMYSaV-inoculated shallot plants resulted in infected plants, showing that the virus can perpetuate the infection over seasons. A survey analyzing 351 shallot samples over a four years period strongly suggests an association of SMYSaV with the disease symptoms. An analysis of SMYSaV diversity indicates the existence of two clusters of isolates, one of which is largely predominant in the field over years.

Lavender Decline in France Is Associated with Chronic Infection by Lavender-Specific Strains of "Candidatus Phytoplasma solani".

Lavender decline compromises French lavender production, and preliminary data have suggested the involvement of "Candidatus Phytoplasma solani" in the etiology of the disease. In order to evaluate the epidemiological role of "Ca Phytoplasma solani," a 3-year survey was conducted in southeastern France. "Ca Phytoplasma solani" was detected in 19 to 56% of the declining plants, depending on seasons and cultivars, and its prevalence was correlated with symptom severity. Autumn was more favorable than spring for phytoplasma detection, and "Ca Phytoplasma solani" incidence was higher in Lavandula angustifolia than in Lavandula intermedia hybrids. Detection of the phytoplasma fluctuated over months, supporting the chronicity of infection. Three "Ca Phytoplasma solani" secY genotypes, S17, S16, and S14, were the most prevalent in lavender fields and were also detected in nurseries, whereas strains detected in surrounding bindweed and wild carrots were mostly of the S1 and S4 genotypes. This suggests that lavender is the main pathogen reservoir of the epidemic. Adults and nymphs of the planthopper vector Hyalesthes obsoletus were commonly captured in lavender fields and were shown to harbor mainly the prevalent phytoplasma genotypes detected in lavenders. The "Ca Phytoplasma solani" genotype S17 was transmitted to Catharanthus roseus periwinkle by naturally infected H. obsoletus Finally, the inventory of the bacterial community of declining lavenders that tested negative for "Ca Phytoplasma solani" by 16S rRNA deep sequencing ruled out the involvement of other phloem-limited bacterial pathogens.IMPORTANCE The etiology and main pathways for the spread of lavender decline, an infectious disease affecting French lavender production since the 1960s, have remained unclear, hampering the development of efficient control strategies. An extensive survey of lavender fields led to the conclusion that "Candidatus Phytoplasma solani" was chronically infecting declining lavenders and was associated with large infectious populations of Hyalesthes obsoletus planthoppers living on the crop itself. Lavender appeared to be the main reservoir host for lavender-specific phytoplasma strains, an unusual feature for this phytoplasma, which usually propagates from reservoir weeds to various economically important crops. These results point out the necessity to protect young lavender fields from the initial phytoplasma inoculum coming from surrounding lavender fields or from infected nurseries and to promote agricultural practices that reduce the development of H. obsoletus vector populations.