Nematodes alter the taxonomic and functional profiles of benthic bacterial communities: A metatranscriptomic approach.

Marine sediments cover 70% of the Earth's surface, and harbour diverse bacterial communities critical for marine biogeochemical processes, which affect climate change, biodiversity and ecosystem functioning. Nematodes, the most abundant and species-rich metazoan organisms in marine sediments, in turn, affect benthic bacterial communities and bacterial-mediated ecological processes, but the underlying mechanisms by which they affect biogeochemical cycles remain poorly understood. Here, we demonstrate using a metatranscriptomic approach that nematodes alter the taxonomic and functional profiles of benthic bacterial communities. We found particularly strong stimulation of nitrogen-fixing and methane-oxidizing bacteria in the presence of nematodes, as well as increased functional activity associated with methane metabolism and degradation of various carbon compounds. This study provides empirical evidence that the presence of nematodes results in taxonomic and functional shifts in active bacterial communities, indicating that nematodes may play an important role in benthic ecosystem processes.

Using environmental DNA metabarcoding to monitor fish communities in small rivers and large brooks: Insights on the spatial scale of information.

Monitoring fish communities is central to the evaluation of ecological health of rivers. Both presence/absence of fish species and their relative quantity in local fish assemblages are crucial parameters to measure. Fish communities in lotic systems are traditionally monitored via electrofishing, characterized by a known limited efficiency and high survey costs. Analysis of environmental DNA could serve as a non-destructive alternative for detection and quantification of lotic fish communities, but this approach still requires further insights in practical sampling schemes incorporating transport and dilution of the eDNA particles; optimization of predictive power and quality assurance of the molecular detection method. Via a controlled cage experiment, we aim to extend the knowledge on streamreach of eDNA in small rivers and large brooks, as laid out in the European Water Framework Directive's water typology. Using a high and low source biomass in two river transects of a species-poor river characterized by contrasting river discharge rates, we found strong and significant correlations between the eDNA relative species abundances and the relative biomass per species in the cage community. Despite a decreasing correlation over distance, the underlying community composition remained stable from 25 to 300 m, or up to 1 km downstream of the eDNA source, depending on the river discharge rate. Such decrease in similarity between relative source biomass and the corresponding eDNA-based community profile with increasing distance downstream from the source, might be attributed to variation in species-specific eDNA persistence. Our findings offer crucial insights on eDNA behaviour and characterization of riverine fish communities. We conclude that water sampled from a relatively small river offers an adequate eDNA snapshot of the total fish community in the 300-1000 m upstream transect. The potential application for other river systems is further discussed.

Monitoring of spatiotemporal occupancy patterns of fish and amphibian species in a lentic aquatic system using environmental DNA.

To effectively monitor, manage and protect aquatic species and understand their interactions, knowledge of their spatiotemporal distribution is needed. In this study, we used a fine-scale spatiotemporal water sampling design, followed by environmental DNA (eDNA) 12S metabarcoding, to investigate occupancy patterns of a natural community of fish and amphibian species in a lentic system. In the same system, we experimentally estimated the spatial and temporal dispersion of eDNA by placing a community of different fish and amphibian species in cages at one side of the pond, creating a controlled point of eDNA emission. Analyses of this cage community revealed a sharp spatial decline in detection rates and relative eDNA quantities at a distance of 5-10 m from the source, depending on the species and its abundance. In addition, none of the caged species could be detected 1 week after removal from the system. This indicates high eDNA decay rates and limited spatial eDNA dispersal, facilitating high local resolution for monitoring spatial occupancy patterns of aquatic species. Remarkably, for seven of the nine cage species, the presence of a single individual could be detected by pooling water of subsamples taken across the whole water body, illustrating the high sensitivity of the eDNA sampling and detection method applied. Finally, our work demonstrated that a fine-scale sampling design in combination with eDNA metabarcoding can cover total biodiversity very precisely and allows the construction of consistent spatiotemporal patterns of relative abundance and local distribution of free-living fish and amphibian species in a lentic ecosystem.

Daring to be differential: metabarcoding analysis of soil and plant-related microbial communities using amplicon sequence variants and operational taxonomical units.

BACKGROUND: Microorganisms are not only indispensable to ecosystem functioning, they are also keystones for emerging technologies. In the last 15 years, the number of studies on environmental microbial communities has increased exponentially due to advances in sequencing technologies, but the large amount of data generated remains difficult to analyze and interpret. Recently, metabarcoding analysis has shifted from clustering reads using Operational Taxonomical Units (OTUs) to Amplicon Sequence Variants (ASVs). Differences between these methods can seriously affect the biological interpretation of metabarcoding data, especially in ecosystems with high microbial diversity, as the methods are benchmarked based on low diversity datasets. RESULTS: In this work we have thoroughly examined the differences in community diversity, structure, and complexity between the OTU and ASV methods. We have examined culture-based mock and simulated datasets as well as soil- and plant-associated bacterial and fungal environmental communities. Four key findings were revealed. First, analysis of microbial datasets at family level guaranteed both consistency and adequate coverage when using either method. Second, the performance of both methods used are related to community diversity and sample sequencing depth. Third, differences in the method used affected sample diversity and number of detected differentially abundant families upon treatment; this may lead researchers to draw different biological conclusions. Fourth, the observed differences can mostly be attributed to low abundant (relative abundance < 0.1%) families, thus extra care is recommended when studying rare species using metabarcoding. The ASV method used outperformed the adopted OTU method concerning community diversity, especially for fungus-related sequences, but only when the sequencing depth was sufficient to capture the community complexity. CONCLUSIONS: Investigation of metabarcoding data should be done with care. Correct biological interpretation depends on several factors, including in-depth sequencing of the samples, choice of the most appropriate filtering strategy for the specific research goal, and use of family level for data clustering.

Complex patterns in tolerance and resistance to pests and diseases underpin the domestication of tomato.

A frequent hypothesis explaining the high susceptibility of many crops to pests and diseases is that, in the process of domestication, crops have lost defensive genes and traits against pests and diseases. Ecological theory predicts trade-offs whereby resistance and tolerance go at the cost of each other. We used wild relatives, early domesticated varieties, traditional local landraces and cultivars of tomato (Solanum lycopersicum) to test whether resistance and tolerance trade-offs were phylogenetically structured or varied according to degree of domestication. We exposed tomato genotypes to the aphid Macrosiphum euphorbiae, the cotton leafworm Spodoptera littoralis, the root knot nematode Meloidogyne incognita and two common insect-transmitted plant viruses, and reconstructed their phylogenetic relationships using Genotyping-by-Sequencing. We found differences in the performance and effect of pest and diseases but such differences were not related with domestication degree nor genetic relatedness, which probably underlie a complex genetic basis for resistance and indicate that resistance traits appeared at different stages and in unrelated genetic lineages. Still, wild and early domesticated accessions showed greater resistance to aphids and tolerance to caterpillars, nematodes and diseases than modern cultivars. Our findings help to understand how domestication affects plant-pest interactions and underline the importance of tolerance in crop breeding.

Whole-Genome Deep Sequencing Reveals Host-Driven in-planta Evolution of Columnea Latent Viroid (CLVd) Quasi-Species Populations.

Columnea latent viroid (CLVd) is one of the most serious tomato diseases. In general, viroids have high mutation rates. This generates a population of variants (so-called quasi-species) that co-exist in their host and exhibit a huge level of genetic diversity. To study the population of CLVd in individual host plants, we used amplicon sequencing using specific CLVd primers linked with a sample-specific index sequence to amplify libraries. An infectious clone of a CLVd isolate Chaipayon-1 was inoculated on different solanaceous host plants. Six replicates of the amplicon sequencing results showed very high reproducibility. On average, we obtained 133,449 CLVd reads per PCR-replicate and 79 to 561 viroid sequence variants, depending on the plant species. We identified 19 major variants (>1.0% mean relative abundance) in which a total of 16 single-nucleotide polymorphisms (SNPs) and two single nucleotide insertions were observed. All major variants contained a combination of 4 to 6 SNPs. Secondary structure prediction clustered all major variants into a tomato/bolo maka group with four loops (I, II, IV and V), and a chili pepper group with four loops (I, III, IV and V) at the terminal right domain, compared to the CLVd Chaipayon-1 which consists of five loops (I, II, III, IV and V).

Virus Detection by High-Throughput Sequencing of Small RNAs: Large-Scale Performance Testing of Sequence Analysis Strategies.

Recent developments in high-throughput sequencing (HTS), also called next-generation sequencing (NGS), technologies and bioinformatics have drastically changed research on viral pathogens and spurred growing interest in the field of virus diagnostics. However, the reliability of HTS-based virus detection protocols must be evaluated before adopting them for diagnostics. Many different bioinformatics algorithms aimed at detecting viruses in HTS data have been reported but little attention has been paid thus far to their sensitivity and reliability for diagnostic purposes. Therefore, we compared the ability of 21 plant virology laboratories, each employing a different bioinformatics pipeline, to detect 12 plant viruses through a double-blind large-scale performance test using 10 datasets of 21- to 24-nucleotide small RNA (sRNA) sequences from three different infected plants. The sensitivity of virus detection ranged between 35 and 100% among participants, with a marked negative effect when sequence depth decreased. The false-positive detection rate was very low and mainly related to the identification of host genome-integrated viral sequences or misinterpretation of the results. Reproducibility was high (91.6%). This work revealed the key influence of bioinformatics strategies for the sensitive detection of viruses in HTS sRNA datasets and, more specifically (i) the difficulty in detecting viral agents when they are novel or their sRNA abundance is low, (ii) the influence of key parameters at both assembly and annotation steps, (iii) the importance of completeness of reference sequence databases, and (iv) the significant level of scientific expertise needed when interpreting pipeline results. Overall, this work underlines key parameters and proposes recommendations for reliable sRNA-based detection of known and unknown viruses.

Reduced-Particle-Size Wheat Bran Is Efficiently Colonized by a Lactic Acid-Producing Community and Reduces Levels of Enterobacteriaceae in the Cecal Microbiota of Broilers.

In the present study, we investigated whether reducing the particle size of wheat bran affects the colonizing microbial community using batch fermentations with cecal inocula from seven different chickens. We also investigated the effect of in-feed administration of regular wheat bran (WB; 1,690 µm) and wheat bran with reduced particle size (WB280; 280 µm) on the cecal microbial community composition of broilers. During batch fermentation, WB280 was colonized by a lactic acid-producing community (Bifidobacteriaceae and Lactobacillaceae) and by Lachnospiraceae that contain lactic acid-consuming butyric acid-producing species. The relative abundances of the Enterobacteriaceae decreased in the particle-associated communities for both WB and WB280 compared to that of the control. In addition, the community attached to wheat bran was enriched in xylan-degrading bacteria. When administered as a feed additive to broilers, WB280 significantly increased the richness of the cecal microbiota and the abundance of bacteria containing the butyryl-coenzyme A (CoA):acetate CoA-transferase gene, a key gene involved in bacterial butyrate production, while decreasing the abundances of Enterobacteriaceae family members in the ceca. Particle size reduction of wheat bran thus resulted in the colonization of the bran particles by a very specific lactic acid- and butyric acid-producing community and can be used to steer toward beneficial microbial shifts. This can potentially increase the resilience against pathogens and increase animal performance when the reduced-particle-size wheat bran is administered as a feed additive to broilers.IMPORTANCE Prebiotic dietary fibers are known to improve the gastrointestinal health of both humans and animals in many different ways. They can increase the bulking capacity, improve transit times, and, depending on the fiber, even stimulate the growth and activity of resident beneficial bacteria. Wheat bran is a readily available by-product of flour processing and is a highly concentrated source of (in)soluble dietary fiber. The intake of fiber-rich diets has been associated with increased Firmicutes and decreased Proteobacteria numbers. Here, we show that applying only 1% of a relatively simple substrate which was technically modified using relatively simple techniques reduces the concentration of Enterobacteriaceae This could imply that in future intervention studies, one should take the particle size of dietary fibers into account.

Temporal Dynamics of Bacterial and Fungal Colonization on Plastic Debris in the North Sea.

Despite growing evidence that biofilm formation on plastic debris in the marine environment may be essential for its biodegradation, the underlying processes have yet to be fully understood. Thus, far, bacterial biofilm formation had only been studied after short-term exposure or on floating plastic, yet a prominent share of plastic litter accumulates on the seafloor. In this study, we explored the taxonomic composition of bacterial and fungal communities on polyethylene plastic sheets and dolly ropes during long-term exposure on the seafloor, both at a harbor and an offshore location in the Belgian part of the North Sea. We reconstructed the sequence of events during biofilm formation on plastic in the harbor environment and identified a core bacteriome and subsets of bacterial indicator species for early, intermediate, and late stages of biofilm formation. Additionally, by implementing ITS2 metabarcoding on plastic debris, we identified and characterized for the first time fungal genera on plastic debris. Surprisingly, none of the plastics exposed to offshore conditions displayed the typical signature of a late stage biofilm, suggesting that biofilm formation is severely hampered in the natural environment where most plastic debris accumulates.

Redirection of auxin flow in Arabidopsis thaliana roots after infection by root-knot nematodes.

Plant-parasitic root-knot nematodes induce the formation of giant cells within the plant root, and it has been recognized that auxin accumulates in these feeding sites. Here, we studied the role of the auxin transport system governed by AUX1/LAX3 influx proteins and different PIN efflux proteins during feeding site development in Arabidopsis thaliana roots. Data generated via promoter-reporter line and protein localization analyses evoke a model in which auxin is being imported at the basipetal side of the feeding site by the concerted action of the influx proteins AUX1 and LAX3, and the efflux protein PIN3. Mutants in auxin influx proteins AUX1 and LAX3 bear significantly fewer and smaller galls, revealing that auxin import into the feeding sites is needed for their development and expansion. The feeding site development in auxin export (PIN) mutants was only slightly hampered. Expression of some PINs appears to be suppressed in galls, probably to prevent auxin drainage. Nevertheless, a functional PIN4 gene seems to be a prerequisite for proper nematode development and gall expansion, most likely by removing excessive auxin to stabilize the hormone level in the feeding site. Our data also indicate a role of local auxin peaks in nematode attraction towards the root.

Bacterial Community Profiling of Plastic Litter in the Belgian Part of the North Sea.

Bacterial colonization of marine plastic litter (MPL) is known for over four decades. Still, only a few studies on the plastic colonization process and its influencing factors are reported. In this study, seafloor MPL was sampled at different locations across the Belgian part of the North Sea to study bacterial community structure using 16S metabarcoding. These marine plastic bacterial communities were compared with those of sediment and seawater, and resin pellets sampled on the beach, to investigate the origin and uniqueness of plastic bacterial communities. Plastics display great variation of bacterial community composition, while each showed significant differences from those of sediment and seawater, indicating that plastics represent a distinct environmental niche. Various environmental factors correlate with the diversity of MPL bacterial composition across plastics. In addition, intrinsic plastic-related factors such as pigment content may contribute to the differences in bacterial colonization. Furthermore, the differential abundance of known primary and secondary colonizers across the various plastics may indicate different stages of bacterial colonization, and may confound comparisons of free-floating plastics. Our studies provide insights in the factors that shape plastic bacterial colonization and shed light on the possible role of plastic as transport vehicle for bacteria through the aquatic environment.

Comparative genome analysis of pathogenic and non-pathogenic Clavibacter strains reveals adaptations to their lifestyle.

BACKGROUND: The genus Clavibacter harbors economically important plant pathogens infecting agricultural crops such as potato and tomato. Although the vast majority of Clavibacter strains are pathogenic, there is an increasing number of non-pathogenic isolates reported. Non-pathogenic Clavibacter strains isolated from tomato seeds are particularly problematic because they affect the current detection and identification tests for Clavibacter michiganensis subsp. michiganensis (Cmm), which is regulated with a zero tolerance in tomato seed. Their misidentification as pathogenic Cmm hampers a clear judgment on the seed quality and health. RESULTS: To get more insight in the genetic features linked to the lifestyle of these bacteria, a whole-genome sequence of the tomato seed-borne non-pathogenic Clavibacter LMG 26808 was determined. To gain a better understanding of the molecular determinants of pathogenicity, the genome sequence of LMG 26808 was compared with that of the pathogenic Cmm strain (NCPPB 382). The comparative analysis revealed that LMG 26808 does not contain plasmids pCM1 and pCM2 and also lacks the majority of important virulence factors described so far for pathogenic Cmm. This explains its apparent non-pathogenic nature in tomato plants. Moreover, the genome analysis of LMG 26808 detected sequences from a plasmid originating from a member of Enterobacteriaceae/Klebsiella relative. Genes received that way and coding for antibiotic resistance may provide a competitive advantage for survival of LMG 26808 in its ecological niche. Genetically, LMG 26808 was the most similar to the pathogenic Cmm NCPPB 382 but contained more mobile genetic elements. The genome of this non-pathogenic Clavibacter strain contained also a high number of transporters and regulatory genes. CONCLUSIONS: The genome sequence of the non-pathogenic Clavibacter strain LMG 26808 and the comparative analyses with other pathogenic Clavibacter strains provided a better understanding of the genetic bases of virulence and adaptation mechanisms present in the genus Clavibacter.

Correction: Haegeman et al. Looking beyond Virus Detection in RNA Sequencing Data: Lessons Learned from a Community-Based Effort to Detect Cellular Plant Pathogens and Pests. Plants 2023, 12, 2139.

In the original publication [...].

Transcriptional analysis through RNA sequencing of giant cells induced by Meloidogyne graminicola in rice roots.

One of the reasons for the progressive yield decline observed in aerobic rice production is the rapid build-up of populations of the rice root knot nematode Meloidogyne graminicola. These nematodes induce specialized feeding cells inside root tissue, called giant cells. By injecting effectors in and sipping metabolites out of these cells, they reprogramme normal cell development and deprive the plant of its nutrients. In this research we have studied the transcriptome of giant cells in rice, after isolation of these cells by laser-capture microdissection. The expression profiles revealed a general induction of primary metabolism inside the giant cells. Although the roots were shielded from light induction, we detected a remarkable induction of genes involved in chloroplast biogenesis and tetrapyrrole synthesis. The presence of chloroplast-like structures inside these dark-grown cells was confirmed by confocal microscopy. On the other hand, genes involved in secondary metabolism and more specifically, the majority of defence-related genes were strongly suppressed in the giant cells. In addition, significant induction of transcripts involved in epigenetic processes was detected inside these cells 7 days after infection.

Detection of Fusarium oxysporum f.sp. lactucae race 1 and 4 via race-specific real-time PCR and target enrichment.

Fusarium oxysporum f.sp. lactucae (Fol) causes a vascular disease in lettuce that results in significant yield losses. Race-specific and sensitive real-time PCR assays were developed for Fol races 1 and 4, which are prevalent in Europe. Using genotyping-by-sequencing, unique DNA loci specific to each race were identified and subsequently used for the design of primers and hydrolysis probes. Two assays per race were developed to ensure specificity. The two assays of each race could be run in duplex format, while still giving a sensitivity of 100 fg genomic DNA for all assays. Sample preparation methods were developed for plant tissue, soil, and surfaces, with an extra enrichment step when additional sensitivity was required. By controlling the incubation conditions during the enrichment step, the real-time PCR signal could be matched to the number of spore equivalents in the original sample. When enriching naturally infested soil, down to six conidiospore equivalents L-1 soil could be detected. As enrichment ensures sensitive detection and focuses on living Fol propagules, it facilitates the evaluation of control measures. The developed detection methods for soil and surfaces were applied to samples from commercial lettuce farms and confirmed the prevalence of Fol race 4 in Belgium. Monitoring of soil disinfestation events revealed that despite a dramatic decrease in quantity, the pathogen could still be detected either immediately after sheet steaming or after harvesting the first new crop. The detection method for plant tissue was successfully used to quantify Fol in lettuce inoculated with race 1, race 4 or a combination of both. Under the temperature conditions used, race 4 was more aggressive than race 1, as reflected in larger amounts of DNA of race 4 detected in the roots. These newly developed assays are a promising tool for epidemiological research as well as for the evaluation of control measures.

Identification and characterization of CYP71 subclade cytochrome P450 enzymes involved in the biosynthesis of bitterness compounds in Cichorium intybus.

Industrial chicory (Cichorium intybus var. sativum) and witloof (C. intybus var. foliosum) are crops with an important economic value, mainly cultivated for inulin production and as a leafy vegetable, respectively. Both crops are rich in nutritionally relevant specialized metabolites with beneficial effects for human health. However, their bitter taste, caused by the sesquiterpene lactones (SLs) produced in leaves and taproot, limits wider applications in the food industry. Changing the bitterness would thus create new opportunities with a great economic impact. Known genes encoding enzymes involved in the SL biosynthetic pathway are GERMACRENE A SYNTHASE (GAS), GERMACRENE A OXIDASE (GAO), COSTUNOLIDE SYNTHASE (COS) and KAUNIOLIDE SYNTHASE (KLS). In this study, we integrated genome and transcriptome mining to further unravel SL biosynthesis. We found that C. intybus SL biosynthesis is controlled by the phytohormone methyl jasmonate (MeJA). Gene family annotation and MeJA inducibility enabled the pinpointing of candidate genes related with the SL biosynthetic pathway. We specifically focused on members of subclade CYP71 of the cytochrome P450 family. We verified the biochemical activity of 14 C. intybus CYP71 enzymes transiently produced in Nicotiana benthamiana and identified several functional paralogs for each of the GAO, COS and KLS genes, pointing to redundancy in and robustness of the SL biosynthetic pathway. Gene functionality was further analyzed using CRISPR/Cas9 genome editing in C. intybus. Metabolite profiling of mutant C. intybus lines demonstrated a successful reduction in SL metabolite production. Together, this study increases our insights into the C. intybus SL biosynthetic pathway and paves the way for the engineering of C. intybus bitterness.

Industrial chicory genome gives insights into the molecular timetable of anther development and male sterility.

Industrial chicory (Cichorium intybus var. sativum) is a biannual crop mostly cultivated for extraction of inulin, a fructose polymer used as a dietary fiber. F1 hybrid breeding is a promising breeding strategy in chicory but relies on stable male sterile lines to prevent self-pollination. Here, we report the assembly and annotation of a new industrial chicory reference genome. Additionally, we performed RNA-Seq on subsequent stages of flower bud development of a fertile line and two cytoplasmic male sterile (CMS) clones. Comparison of fertile and CMS flower bud transcriptomes combined with morphological microscopic analysis of anthers, provided a molecular understanding of anther development and identified key genes in a range of underlying processes, including tapetum development, sink establishment, pollen wall development and anther dehiscence. We also described the role of phytohormones in the regulation of these processes under normal fertile flower bud development. In parallel, we evaluated which processes are disturbed in CMS clones and could contribute to the male sterile phenotype. Taken together, this study provides a state-of-the-art industrial chicory reference genome, an annotated and curated candidate gene set related to anther development and male sterility as well as a detailed molecular timetable of flower bud development in fertile and CMS lines.

Phylogeography and population structure of the global, wide host-range hybrid pathogen Phytophthora × cambivora.

Invasive, exotic plant pathogens pose a major threat to native and agricultural ecosystems. Phytophthora × cambivora is an invasive, destructive pathogen of forest and fruit trees causing severe damage worldwide to chestnuts (Castanea), apricots, peaches, plums, almonds and cherries (Prunus), apples (Malus), oaks (Quercus), and beech (Fagus). It was one of the first damaging invasive Phytophthora species to be introduced to Europe and North America, although its origin is unknown. We determined its population genetic history in Europe, North and South America, Australia and East Asia (mainly Japan) using genotyping-by-sequencing. Populations in Europe and Australia appear clonal, those in North America are highly clonal yet show some degree of sexual reproduction, and those in East Asia are partially sexual. Two clonal lineages, each of opposite mating type, and a hybrid lineage derived from these two lineages, dominated the populations in Europe and were predominantly found on fagaceous forest hosts (Castanea, Quercus, Fagus). Isolates from fruit trees (Prunus and Malus) belonged to a separate lineage found in Australia, North America, Europe and East Asia, indicating the disease on fruit trees could be caused by a distinct lineage of P. × cambivora, which may potentially be a separate sister species and has likely been moved with live plants. The highest genetic diversity was found in Japan, suggesting that East Asia is the centre of origin of the pathogen. Further surveys in unsampled, temperate regions of East Asia are needed to more precisely identify the location and range of the centre of diversity.

Detection of single nucleotide polymorphisms in virus genomes assembled from high-throughput sequencing data: large-scale performance testing of sequence analysis strategies.

Recent developments in high-throughput sequencing (HTS) technologies and bioinformatics have drastically changed research in virology, especially for virus discovery. Indeed, proper monitoring of the viral population requires information on the different isolates circulating in the studied area. For this purpose, HTS has greatly facilitated the sequencing of new genomes of detected viruses and their comparison. However, bioinformatics analyses allowing reconstruction of genome sequences and detection of single nucleotide polymorphisms (SNPs) can potentially create bias and has not been widely addressed so far. Therefore, more knowledge is required on the limitations of predicting SNPs based on HTS-generated sequence samples. To address this issue, we compared the ability of 14 plant virology laboratories, each employing a different bioinformatics pipeline, to detect 21 variants of pepino mosaic virus (PepMV) in three samples through large-scale performance testing (PT) using three artificially designed datasets. To evaluate the impact of bioinformatics analyses, they were divided into three key steps: reads pre-processing, virus-isolate identification, and variant calling. Each step was evaluated independently through an original, PT design including discussion and validation between participants at each step. Overall, this work underlines key parameters influencing SNPs detection and proposes recommendations for reliable variant calling for plant viruses. The identification of the closest reference, mapping parameters and manual validation of the detection were recognized as the most impactful analysis steps for the success of the SNPs detections. Strategies to improve the prediction of SNPs are also discussed.

Looking beyond Virus Detection in RNA Sequencing Data: Lessons Learned from a Community-Based Effort to Detect Cellular Plant Pathogens and Pests.

High-throughput sequencing (HTS), more specifically RNA sequencing of plant tissues, has become an indispensable tool for plant virologists to detect and identify plant viruses. During the data analysis step, plant virologists typically compare the obtained sequences to reference virus databases. In this way, they are neglecting sequences without homologies to viruses, which usually represent the majority of sequencing reads. We hypothesized that traces of other pathogens might be detected in this unused sequence data. In the present study, our goal was to investigate whether total RNA-seq data, as generated for plant virus detection, is also suitable for the detection of other plant pathogens and pests. As proof of concept, we first analyzed RNA-seq datasets of plant materials with confirmed infections by cellular pathogens in order to check whether these non-viral pathogens could be easily detected in the data. Next, we set up a community effort to re-analyze existing Illumina RNA-seq datasets used for virus detection to check for the potential presence of non-viral pathogens or pests. In total, 101 datasets from 15 participants derived from 51 different plant species were re-analyzed, of which 37 were selected for subsequent in-depth analyses. In 29 of the 37 selected samples (78%), we found convincing traces of non-viral plant pathogens or pests. The organisms most frequently detected in this way were fungi (15/37 datasets), followed by insects (13/37) and mites (9/37). The presence of some of the detected pathogens was confirmed by independent (q)PCRs analyses. After communicating the results, 6 out of the 15 participants indicated that they were unaware of the possible presence of these pathogens in their sample(s). All participants indicated that they would broaden the scope of their bioinformatic analyses in future studies and thus check for the presence of non-viral pathogens. In conclusion, we show that it is possible to detect non-viral pathogens or pests from total RNA-seq datasets, in this case primarily fungi, insects, and mites. With this study, we hope to raise awareness among plant virologists that their data might be useful for fellow plant pathologists in other disciplines (mycology, entomology, bacteriology) as well.