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.

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.

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.

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.

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.

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.

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.

Determination of the complete genomic sequence of grapevine virus H, a novel vitivirus infecting grapevine.

The present work reports the discovery and the complete genome sequencing of a novel member of the genus Vitivirus in the family Betaflexiviridae (subfamily Trivirinae) from a symptomless grapevine of unknown variety from Portugal. Total RNAs extracted from phloem scrapings were sequenced using Illumina technology. Bioinformatic analysis of the RNA-seq data revealed a mixed infection involving three viruses and two viroids in addition to a novel vitivirus. Completion and analysis of the genome sequence (7446 nt excluding the polyA tail) showed a typical vitivirus genomic organization. Phylogenetic analysis of the various ORFs clearly showed the new virus to belong in the genus Vitivirus, but sequence divergence firmly establishes it as a member of a new species, for which the name "Grapevine virus H" is proposed.

Complete genomic sequence of Raphanus sativus cryptic virus 4 (RsCV4), a novel alphapartitivirus from radish.

The present work reports the discovery and complete genome sequencing of a virus from symptomless radish seedlings, classifiable as a novel member of the genus Alphapartitivirus, family Partitiviridae. Total RNA extracted from germinating seedlings was sequenced using Illumina technology. Bioinformatic analysis of the RNA-seq data revealed two contigs representing the near full-length genomic sequences of two genomic RNAs representing a new virus. Analysis of the genome sequence (excluding the polyA tail, RNA1: 1976 nt and RNA2: 1751 nt, respectively) showed a genomic organization typical of viruses classed within the Partitiviridae, with each genomic RNA encoding a single open reading frame (ORF). Phylogenetic analysis of the RNA dependent RNA polymerase (RNA1 ORF) and of the capsid protein (RNA2 ORF) clearly showed the new virus can be classified within the genus Alphapartitivirus, but sequence divergence establishes it as a new species, for which the name "Raphanus sativus cryptic virus 4" is proposed.

Genome characterization of a divergent isolate of the mycovirus Botrytis virus F from a grapevine metagenome.

As part of a grapevine metagenome study, total RNA extracted from grapevine phloem scrapings was analyzed by Illumina sequencing. For one 420A rootstock sample, reads mapping against a reference database and BLAST annotation of contigs identified the presence of a divergent isolate of Botrytis virus F (BVF). The full genome sequence of this isolate (IVC-5-77) was determined (6,828 nucleotides [nt], excluding the poly(A) tail) and shown to be collinear with that of the BVF reference isolate, with the two open reading frames encoding a replication-associated protein (REP) and a coat protein (CP). The IVC-5-77 isolate, however, is very divergent, showing only 81.3-81.6% nucleotide sequence identity to the two other sequenced BVF isolates. The internal non-coding region was also found to be highly variable between BVF isolates. Analysis of the RNASeq reads demonstrated that close to 20% of them belong to Botrytis cinerea, the putative host of the IVC-5-77 isolate. These results extend our knowledge of the diversity and variability of BVF and demonstrate its detectability, together with that of its B. cinerea host, in total RNA RNASeq data from grapevine phloem scrapings.

Complete genome sequence of lettuce chordovirus 1 isolated from cultivated lettuce in France.

Double-stranded RNAs purified from cultivated (Lactuca sativa) or wild (L. serriola) lettuce from southwest France were analyzed by high-throughput sequencing. For both samples, BLAST annotation revealed contigs with homology to Betaflexiviridae family members. The full genome sequence of the isolate from cultivated lettuce (JG1) was completed (8,536 nucleotides [nt], excluding the poly(A) tail). The sequence of the 3' half of the genome (4,800 nt) of a wild lettuce isolate (P22) was determined and found to share 95.1% nt sequence identity with the JG1 isolate. The JG1 genome contains four open reading frames, encoding a replicase, a movement protein, a capsid protein, and a protein of unknown function, respectively. Based on genome organization and phylogenetic relationships, the lettuce virus is most closely related to the recently described carrot chordoviruses 1 and 2 in the family Betaflexiviridae. Considering the species demarcation criteria in this family, the two lettuce viruses represent isolates of a new chordovirus species for which the name "lettuce chordovirus 1" (LeCV1) is proposed.

Molecular characterization of a novel fusarivirus infecting the plant-pathogenic fungus Neofusicoccum luteum.

Double-stranded RNAs from an isolate of Neofusicoccum luteum collected from grapevines were analyzed by high-throughput sequencing. Contig annotations revealed the presence of a potential novel virus belonging to the newly proposed family Fusariviridae. Completion of the viral genome sequence was performed. The genome is 6,244 nucleotide long, excluding the poly(A) tail and contains two putative open reading frames (ORFs). The first one encodes a large polypeptide of 1,552 amino acids (aa) with conserved RNA-dependent RNA polymerase and helicase domains typical of viral replicases. The second ORF encodes a putative 475-aa-long polypeptide showing weak homology to the corresponding ORF of Macrophomina phaseolina single-stranded RNA virus 1, for which no function is known so far. Phylogenetic analyses indicated that this virus should be considered a novel mycovirus belonging to the proposed family Fusariviridae, for which the name "Neofusicoccum luteum fusarivirus 1" (NlFV1) is proposed.

Determination of the complete genomic sequence of Neofusicoccum luteum mitovirus 1 (NLMV1), a novel mitovirus associated with a phytopathogenic Botryosphaeriaceae.

Neofusicoccum luteum species belongs to the Botryosphaeriaceae family and is involved in grapevine wood decay diseases. The present study reports the discovery and the molecular characterization of a novel mitovirus infecting this fungus. Double-stranded RNAs were purified from cultivated N. luteum and analysed by next generation sequencing. Using contigs showing BlastX homology with the RNA-dependent RNA polymerase (RdRp) gene of various members of the family Narnaviridae, a single contig of approximately 1.2 kb was constructed. The genomic sequence was completed and phylogenetic analyses indicated that this virus represents a new member of the genus Mitovirus, for which the name of "Neofusicoccum luteum mitovirus 1" is proposed. The genome is 2,389 nucleotides long and, based on the fungal mitochondrial genetic code, it encodes a putative protein of 710 amino acids, homologous to the RdRps of members of the Narnaviridae family. The neofusicoccum luteus mitovirus 1 (NLMV1) RdRp contains the six conserved motifs previously reported for mitoviral RdRps. Our findings represent the first evidence that a mycovirus can infect N. luteum, an important pathogenic fungus of grapevine.

Two novel Alphaflexiviridae members revealed by deep sequencing of the Vanilla (Orchidaceae) virome.

The genomes of two novel viruses were assembled from 454 pyrosequencing data obtained from vanilla leaves from La Réunion. Based on genome organization and homologies, one agent was unambiguously classified as a member of the genus Potexvirus and named vanilla virus X (VVX). The second one, vanilla latent virus (VLV), is phylogenetically close to three unclassified members of the family Alphaflexiviridae with similarity to allexiviruses, and despite the presence of an additional 8-kDa open reading frame, we propose to include VLV as a new member of the genus Allexivirus. Both VVX and VLV were mechanically transmitted to vanilla plants, resulting in asymptomatic infections.

Complete genomic sequence of barley (Hordeum vulgare) endornavirus (HvEV) determined by next-generation sequencing.

Endornaviruses are unusual plant-, fungus- and oomycete-infecting viruses with a large, ca 14- to 17-kb linear double-stranded RNA (dsRNA) genome and a persistent lifestyle. The complete genome sequence of an endornavirus from the barley (Hordeum vulgare) Nerz variety was determined from paired Illumina MySeq reads derived from purified dsRNAs. The genome is 14,243 nt long, with 5' and 3' non-coding regions of 207 and 47 nt, respectively. It encodes a single large protein of 4663 amino acids that carries conserved motifs for a methyltransferase, a helicase and an RNA-dependent RNA polymerase. The sequence of Hordeum vulgare endornavirus (HvEV) carries all the hallmarks of a typical member of the genus Endornavirus, with the exception of an UDP-glycosyltransferase motif observed in many, but not all, endornaviral genomes.

Beet western yellows virus infects the carnivorous plant Nepenthes mirabilis.

Although poleroviruses are known to infect a broad range of higher plants, carnivorous plants have not yet been reported as hosts. Here, we describe the first polerovirus naturally infecting the pitcher plant Nepenthes mirabilis. The virus was identified through bioinformatic analysis of NGS transcriptome data. The complete viral genome sequence was assembled from overlapping PCR fragments and shown to share 91.1 % nucleotide sequence identity with the US isolate of beet western yellows virus (BWYV). Further analysis of other N. mirabilis plants revealed the presence of additional BWYV isolates differing by several insertion/deletion mutations in ORF5.

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.

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.

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.

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.