Metagenomic Survey Research Articles

Global metagenomic survey identifies sewage-derived hgcAB+ microorganisms as key contributors to riverine methylmercury production

Methylmercury (MeHg) in aquatic systems poses a serious public health risk through bioaccumulation in the aquatic food web. In recent years, MeHg has been observed to increase to concerning levels globally in rivers near cities; however, the causes of this increase are not well understood. Here, we demonstrate the significant role of sewage contamination by analyzing over 1,300 publicly available metagenomes in urban rivers worldwide, and conducting experiments with water samples across China. We find that sewage contamination significantly increases the abundance of mercury (Hg)-methylating microorganisms in urban rivers globally. This increase is primarily attributed to the high abundance of active Hg-methylating microorganisms in sewage, which migrate to rivers via direct discharge or combined sewer overflows (CSOs), becoming key contributors to elevated riverine MeHg levels. Our findings underscore the importance of effectively eliminating Hg-methylating microorganisms from sewage to mitigate the public health risks associated with MeHg in urban rivers.

Metagenomic survey reveals hydrocarbon biodegradation potential of Canadian high Arctic beaches

BackgroundDecreasing sea ice coverage across the Arctic Ocean due to climate change is expected to increase shipping activity through previously inaccessible shipping routes, including the Northwest Passage (NWP). Changing weather conditions typically encountered in the Arctic will still pose a risk for ships which could lead to an accident and the uncontrolled release of hydrocarbons onto NWP shorelines. We performed a metagenomic survey to characterize the microbial communities of various NWP shorelines and to determine whether there is a metabolic potential for hydrocarbon degradation in these microbiomes.ResultsWe observed taxonomic and functional gene evidence supporting the potential of NWP beach microbes to degrade various types of hydrocarbons. The metagenomic and metagenome-assembled genome (MAG) taxonomy showed that known hydrocarbon-degrading taxa are present in these beaches. Additionally, we detected the presence of biomarker genes of aerobic and anaerobic degradation pathways of alkane and aromatic hydrocarbons along with complete degradation pathways for aerobic alkane degradation. Alkane degradation genes were present in all samples and were also more abundant (33.8 ± 34.5 hits per million genes, HPM) than their aromatic hydrocarbon counterparts (11.7 ± 12.3 HPM). Due to the ubiquity of MAGs from the genus Rhodococcus (23.8% of the MAGs), we compared our MAGs with Rhodococcus genomes from NWP isolates obtained using hydrocarbons as the carbon source to corroborate our results and to develop a pangenome of Arctic Rhodococcus. Our analysis revealed that the biodegradation of alkanes is part of the core pangenome of this genus. We also detected nitrogen and sulfur pathways as additional energy sources and electron donors as well as carbon pathways providing alternative carbon sources. These pathways occur in the absence of hydrocarbons allowing microbes to survive in these nutrient-poor beaches.ConclusionsOur metagenomic analyses detected the genetic potential for hydrocarbon biodegradation in these NWP shoreline microbiomes. Alkane metabolism was the most prevalent type of hydrocarbon degradation observed in these tidal beach ecosystems. Our results indicate that bioremediation could be used as a cleanup strategy, but the addition of adequate amounts of N and P fertilizers, should be considered to help bacteria overcome the oligotrophic nature of NWP shorelines.

Viral metagenomic survey of Caspian seals.

Viral diseases of pinnipeds cause substantial mortality and morbidity and can influence population demography. Viral metagenomic studies can therefore play an important role in pinniped health assessments and disease surveillance relevant to both individual species and in a "One Health" context. This study used a metagenomic approach with high throughput sequencing to make the first assessment of viral diversity in Caspian seals (Pusa caspica), the only marine mammal species endemic to the Caspian Sea. Sequencing libraries from 35 seals sampled 2009-2020 were analysed, finding sequences from the viral families Circoviridae, Parvoviridae, Herpesviridae, Papillomaviridae, Picornaviridae, Caliciviridae, Cruciviridae, Anelloviridae, Smacoviridae, and Orthomyxoviridae, with additional detection of Adenoviridae via PCR. The similarity of viral contigs from Caspian seal to sequences recovered from other pinnipeds ranged from 63.74% (San Miguel sea lion calicivirus) to 78.79% (Seal anellovirus 4). Some findings represent novel viral species, but overall, the viral repertoire of Caspian seals is similar to available viromes from other pinnipeds. Among the sequences recovered were partial contigs for influenza B, representing only the second such molecular identification in marine mammals. This work provides a foundation for further studies of viral communities in Caspian seals, the diversity of viromes in pinnipeds more generally, and contributes data relevant for disease risk assessments in marine mammals.

Entomophthovirus: an insect-derived iflavirus that infects a behavior-manipulating fungal pathogen of dipterans.

We report a virus infecting Entomophthora muscae, a behavior-manipulating fungal pathogen of dipterans. The virus, which we name Berkeley Entomophthovirus, is a positive-strand RNA virus in the iflaviridae family of capsid-forming viruses, which are mostly known to infect insects. The viral RNA is expressed at high levels in fungal cells in vitro and during in vivo infections of Drosophila melanogaster, and virus particles can be seen intracellularly in E. muscae. This virus, of which we find two closely related variants in our culture of E. muscae, is also closely related to three different viruses reported from metagenomic surveys, two of which were isolated from wild dipterans, and a third isolated from wild ticks. By analyzing sequencing data from these earlier reports, we find abundant reads aligning to E. muscae specifically in the samples from which viral reads were sequenced. These data establish a wide and perhaps obligate association with E. muscae in the wild, consistent with our laboratory data that E. muscae is the host for these closely related viruses. Because of this, we propose the name Entomophthovirus (EV) for this group of highly related virus variants. As other members of the iflaviridae have been reported to cause behavioral changes in insects, we speculate on the possibility that EV plays a role in the behavioral manipulation of flies infected with E. muscae.

Highly diverse RNA viruses and phage sequences concealed within birds.

The diversity of birds in most parts of the world is very high, and thus, they may carry different types of highly differentiated and unknown viruses. Thanks to advanced sequencing technologies, studies on the diversity of bird-associated viruses have increased over the past few years. In this study, a large-scale viral metagenomics survey was performed on cloacal swabs of 2,990 birds from nine provinces of the Chinese mainland. To detect undescribed RNA viruses in birds, more than 1,800 sequences sharing relatively low (<60%) amino acid sequence identity with the best match in the GenBank database were screened. Potentially novel viruses related to vertebrates have been identified, and several potential recombination signals were found. Additionally, hundreds of RNA viral sequences related to plants, fungi, and insects were detected, including previously unknown viruses. Furthermore, we investigated the novelty, functionality, and classification of the phages examined in this study. These viruses occupied topological positions on the evolutionary trees to a certain extent and might form novel putative families, genera, or species, thus providing information to fill the phylogenetic gaps of related viruses. These findings provided new insights into bird-associated viruses, but the interactions among these viruses remain unknown and require further investigation.IMPORTANCEStudying the diversity of RNA viruses in birds and mammals is crucial due to their potential impact on human health and the global ecosystem. Many RNA viruses, such as influenza and coronaviruses, have been shown to cross the species barrier and cause zoonotic diseases. In this metagenomics study involving 2,990 birds from at least 82 species, we identified over 1,800 RNA sequences with distant relationships to known viruses, some of which are rare in birds. The study highlights the scope and diversity of RNA viruses in birds, providing data to predict disease risks and monitor potential viral threats to wildlife, livestock, and human health. This information can aid in the development of strategies for disease prevention and control.

Known and novel viruses in Belgian honey bees: yearly differences, spatial clustering, and associations with overwintering loss.

The western honey bee (Apis mellifera) is a highly effective pollinator of flowering plants, including many crops, which gives honey bees an outstanding importance both ecologically and economically. Alarmingly high annual loss rates of managed honey bee colonies are a growing concern for beekeepers and scientists and have prompted a significant research effort toward bee health. Several detrimental factors have been identified, such as varroa mite infestation and disease from various bacterial and viral agents, but annual differences are often not elucidated. In this study, we utilize the viral metagenomic survey of the EPILOBEE project, a European research program for bee health, to elaborate on the most abundant bee viruses of Flanders. We complement the existing metagenomic data with absolute viral loads and their spatial and temporal distributions. Furthermore, we identify Apis orthomyxovirus 1 as a potentially emerging pathogen, as we find evidence for its active replication honey bees.

An integrated metagenomic, metabolomic and transcriptomic survey of Populus across genotypes and environments

Bridging molecular information to ecosystem-level processes would provide the capacity to understand system vulnerability and, potentially, a means for assessing ecosystem health. Here, we present an integrated dataset containing environmental and metagenomic information from plant-associated microbial communities, plant transcriptomics, plant and soil metabolomics, and soil chemistry and activity characterization measurements derived from the model tree species Populus trichocarpa. Soil, rhizosphere, root endosphere, and leaf samples were collected from 27 different P. trichocarpa genotypes grown in two different environments leading to an integrated dataset of 318 metagenomes, 98 plant transcriptomes, and 314 metabolomic profiles that are supported by diverse soil measurements. This expansive dataset will provide insights into causal linkages that relate genomic features and molecular level events to system-level properties and their environmental influences.

Effects of operational parameters on bacterial communities in Hong Kong and global wastewater treatment plants.

Wastewater treatment plants (WWTPs) are an indispensable component of modern cities, as they can remove pollutants in wastewater to prevent anthropogenic activities. Activated sludge (AS) is a fundamental wastewater treatment process and it harbors a highly complex microbial community that forms the main components and contains functional groups. Unveiling "who is there" is a long-term goal of the research on AS microbiology. High-throughput sequencing provides insights into the inventory diversity of microbial communities to an unprecedented level of detail. At present, the analysis of communities in WWTPs usually comes from a specific WWTP and lacks comparisons and verification among different WWTPs. The wide-scale and long-term sampling project and research in this study could help us evaluate the AS community more accurately to find the similarities and different results for different WWTPs in Hong Kong and other regions of the world.

Metagenomic survey of viral diversity obtained from feces of piglets with diarrhea

Pigs are natural host to various zoonotic pathogens including viruses. In this study, we analyzed the viral communities in the feces of 89 piglets with diarrhea under one month old which were collected from six farms in Jiangsu Province of the Eastern China, using the unbiased virus metagenomic method. A total of 89 libraries were constructed, and 46937894 unique sequence reads were generated by Illumina sequencing. Overall, the family Picornaviridae accounted for the majority of the total reads of putative mammalian viruses. Ten novel virus genomes from different family members were discovered, including Parvoviridae (n = 2), Picobirnaviridae (n = 4) and CRESS DNA viruses (n = 4). A large number of phages were identified, which mainly belonged to the order Caudovirales and the family Microviridae. Moreover, some identified viruses were closely related to viruses found in non-porcine hosts, highlighting the potential for cross-species virus dissemination. This study increased our understanding of the fecal virus communities of diarrhea piglets and provided valuable information for virus monitoring and preventing.

Inactive hydrothermal vent microbial communities are important contributors to deep ocean primary productivity.

Active hydrothermal vents are oases for productivity in the deep ocean, but the flow of dissolved substrates that fuel such abundant life ultimately ceases, leaving behind inactive mineral deposits. The rates of microbial activity on these deposits are largely unconstrained. Here we show primary production occurs on inactive hydrothermal deposits and quantify its contribution to new organic carbon production in the deep ocean. Measured incorporation of 14C-bicarbonate shows that microbial communities on inactive deposits fix inorganic carbon at rates comparable to those on actively venting deposits. Single-cell uptake experiments and nanoscale secondary ion mass spectrometry showed chemoautotrophs comprise a large fraction (>30%) of the active microbial cells. Metagenomic and lipidomic surveys of inactive deposits further revealed that the microbial communities are dominated by Alphaproteobacteria and Gammaproteobacteria using the Calvin-Benson-Bassham pathway for carbon fixation. These findings establish inactive vent deposits as important sites for microbial activity and organic carbon production on the seafloor.

A metagenomic survey of bacterial communities from kurut: The fermented cow milk in Kyrgyzstan.

Kurut is a traditional dry dairy product mostly consumed in Central Asia. In this study, the distribution of the dominant bacteria present in kurut samples (n=84) originated from seven (Chuy, Issyk-Kul, Talas, Naryn, Jalal-Abad, Osh, and Batken) regions in Kyrgyzstan were analyzed with Illumina iSeq100 platform. The dominant phylum detected was Firmicutes followed by Proteobacteria, Actinobacteria, Cyanobacteria/Chloroplast, and Tenericutes. The most abundant family detected was Lactobacillaceae followed by Streptococcaceae, Enterococcaceae, Chloroplast, and Leuconostocaceae. At the genus level, Lactobacillus was the predominant one in samples and Streptococcus, Enterococcus, Lactococcus, and Streptophyta followed this. Further comprehensive characterization analyses in kurut samples may have potential applications both in industrial starter culture developments and also future therapeutic approaches based on potential strains with probiotic properties.

Observations on the Antagonistic Relationships between Fungi, Archaea, and Bacteria in Livingston Island (Maritime Antarctica) with the Use of Amplicon-Based Metagenomics

An amplicon-based metagenomic survey of archaea, fungi, and bacteria was performed on Livingston Island, Maritime Antarctica. In many of the samples, patterns of antagonism between these three superkingdoms were observed in the form of an inversely proportional dependence of the richnesses of the three types of microorganisms. The antagonism was quantified—based on the observed numbers of the total tags and the numbers of the operational taxonomic units (OTUs), and on four alpha diversity parameters—using the Shannon, the Simpson, the Chao1, and the ACE indices. We found that the most discriminative results in the antagonism measuring were obtained when the numbers of the OTUs and the ACE community richness estimator were compared. The antagonism between archaea and fungi was most potent, followed by that of archaea and bacteria. The fungi–bacteria antagonism was slightly detectable. Pearson and Spearman correlation analyses also showed a statistically significant negative correlation between the fungal and archaeal effective tags, while the correlation between archaeal and bacterial diversity was positive. Indications of the order of primary microbial succession in barren ecological niches were also observed, demonstrating that archaea and bacteria are the pioneers, followed by fungi, which would displace archaea over time.

Hidden genomic diversity drives niche partitioning in a cosmopolitan eukaryotic picophytoplankton.

Marine eukaryotic phytoplankton are fundamental to the marine food web, yet the lack of reference genomes or just a single genome representing a taxon has led to an underestimation of their taxonomic, adaptive, and functional diversity. Here, we integrated strain isolation with metagenomic binning to recover genomes from the cosmopolitan picophytoplankton genus Bathycoccus, traditionally considered monospecific. Our recovery and analysis of 37 Bathycoccus genomes delineated their global genomic diversity and established four evolutionary clades (BI, BII, BIII, BIV). Our metagenomic abundance survey revealed well-differentiated ecological niches and distinct biogeographic distributions for each clade, predominantly shaped by temperature, salinity, and nutrient availability. Comparative genomics analyses further revealed clade-specific genomic traits that underpin niche adaptation and contribute to the global prevalence of Bathycoccus. Our findings underscore temperature as a major driver of genome diversification in this genus, with clade divergences coinciding with major paleoclimatic events that influenced their contemporary thermal niches. Moreover, the unique enrichment of C2H2 zinc finger and ankyrin repeat gene families in polar-adapted clades suggests previously unrecognized cold-adaptation mechanisms in marine eukaryotic phytoplankton. Our study offers a comprehensive genomic landscape of this crucial eukaryotic picophytoplankton, providing insights into their microdiversity and adaptive evolution in response to changing environments.

Good and bad dispositions between archaea and bacteria in the human gut: New insights from metagenomic survey and co-occurrence analysis

Archaea are an understudied component of the human microbiome. In this study, the gut archaeome and bacteriome of 60 healthy adults from different region were analyzed by whole-genome shotgun sequencing. Archaea were ubiquitously found in a wide range of abundances, reaching up to 7.2 %. The dominant archaeal phylum was Methanobacteriota, specifically the family Methanobacteriaceae, encompassing more than 50 % of Archaea in 50 samples. The previously underestimated Thermoplasmatota, mostly composed of Methanomassiliicoccaceae, dominated in 10 subjects (>50 %) and was present in all others except one. Halobacteriota, the sole other archaeal phylum, occurred in negligible concentration, except for two samples (4.6–4.8 %). This finding confirmed that the human gut archaeome is primarily composed of methanogenic organisms and among the known methanogenic pathway: i) hydrogenotrophic reduction of CO2 is the predominant, being the genus Methanobrevibacter and the species Methanobrevibacter smithii the most abundant in the majority of the samples; ii) the second pathway, that involved Methanomassiliicoccales, was the hydrogenotrophic reduction of methyl-compounds; iii) dismutation of acetate or methyl-compounds seemed to be absent. Co-occurrence analysis allowed to unravel correlations between Archaea and Bacteria that shapes the overall structure of the microbial community, allowing to depict a clearer picture of the human gut archaeome.

Comparative single-cell genomics of Atribacterota JS1 in the Japan Trench hadal sedimentary biosphere.

Deep-sea and subseafloor sedimentary environments host heterotrophic microbial communities that contribute to Earth's carbon cycling. However, the potential metabolic functions of individual microorganisms and their biogeographical distributions in hadal ocean sediments remain largely unexplored. In this study, we conducted single-cell genome sequencing on sediment samples collected from six sites (7,445-8,023 m water depth) along an approximately 500 km transect of the Japan Trench during the International Ocean Discovery Program Expedition 386. A total of 1,886 single-cell amplified genomes (SAGs) were obtained, offering comprehensive genetic insights into sedimentary microbial communities in surface sediments (<1 m depth) above the sulfate-methane transition zone along the Japan Trench. Our genome data set included 269 SAGs from Atribacterota JS1, the predominant bacterial clade in these hadal environments. Phylogenetic analysis classified SAGs into nine distinct phylotypes, whereas metagenome-assembled genomes were categorized into only two phylotypes, advancing JS1 diversity coverage through a single cell-based approach. Comparative genomic analysis of JS1 lineages from different habitats revealed frequent detection of genes related to organic carbon utilization, such as extracellular enzymes like clostripain and α-amylase, and ABC transporters of oligopeptide from Japan Trench members. Furthermore, specific JS1 phylotypes exhibited a strong correlation with in situ methane concentrations and contained genes involved in glycine betaine metabolism. These findings suggest that the phylogenomically diverse and novel Atribacterota JS1 is widely distributed in Japan Trench sediment, playing crucial roles in carbon cycling within the hadal sedimentary biosphere.IMPORTANCEThe Japan Trench represents tectonically active hadal environments associated with Pacific plate subduction beneath the northeastern Japan arc. This study, for the first time, documented a large-scale single-cell and metagenomic survey along an approximately 500 km transect of the Japan Trench, obtaining high-quality genomic information on hadal sedimentary microbial communities. Single-cell genomics revealed the predominance of diverse JS1 lineages not recoverable through conventional metagenomic binning. Their metabolic potential includes genes related to the degradation of organic matter, which contributes to methanogenesis in the deeper layers. Our findings enhance understanding of sedimentary microbial communities at water depths exceeding 7,000 m and provide new insights into the ecological role of biogeochemical carbon cycling in the hadal sedimentary biosphere.

Metagenomic survey of antimicrobial resistance (AMR) in Maryland surface waters differentiated by high and low human impact.

Here, we examine surface waters as a modality to better understand baseline antimicrobial resistance (AMR) across the environment to supplement existing AMR monitoring in pathogens associated with humans, foods, and animals. Data from metagenomic and quasimetagenomic (shotgun sequenced enrichments) are used to describe AMR in Maryland surface waters from high and low human impact classifications.

Oxygen respiration and polysaccharide degradation by a sulfate-reducing acidobacterium

Sulfate-reducing microorganisms represent a globally important link between the sulfur and carbon cycles. Recent metagenomic surveys expanded the diversity of microorganisms putatively involved in sulfate reduction underscoring our incomplete understanding of this functional guild. Here, we use genome-centric metatranscriptomics to study the energy metabolism of Acidobacteriota that carry genes for dissimilation of sulfur compounds in a long-term continuous culture running under alternating anoxic and oxic conditions. Differential gene expression analysis reveals the unique metabolic flexibility of a pectin-degrading acidobacterium to switch from sulfate to oxygen reduction when shifting from anoxic to oxic conditions. The combination of facultative anaerobiosis and polysaccharide degradation expands the metabolic versatility among sulfate-reducing microorganisms. Our results highlight that sulfate reduction and aerobic respiration are not mutually exclusive in the same organism, sulfate reducers can mineralize organic polymers, and anaerobic mineralization of complex organic matter is not necessarily a multi-step process involving different microbial guilds but can be bypassed by a single microbial species.

Identification of a new polyomavirus in distinct human populations and tissues.

Several human polyomaviruses (HPyVs) have been described in the last 15 years. This work aimed to characterize a novel HPyV with cutaneous tropism. Swabs of healthy skin (forehead) of 75 immunocompetent individuals from Argentina were screened for HPyV through sequence amplification techniques. Publicly available metagenomic data sets were also analyzed. A previously unknown polyomavirus sequence was detected in two skin swab samples. A nearly identical sequence was detected in public data sets representing metagenomic surveys of human skin and feces. Further analyses showed that the new polyomavirus diverges from its nearest relative, human polyomavirus 6 (HPyV6), by 17.3%-17.7% (in nucleotides for the large T antigen), which meets criteria for a new species designation in the genus Deltapolyomavirus. The screening also revealed more distant HPyV6 relatives in macaque genital and chimpanzee fecal data sets. Since polyomaviruses are generally thought to cospeciate with mammalian hosts, the high degree of similarity to HPyV6 suggests the new polyomavirus species is human-tropic. Therefore, a novel polyomavirus was identified and characterized from samples of distinct populations and tissues. We suggest the common name human polyomavirus 16 (HPyV16).

Differential effects of wastewater treatment plant effluents on the antibiotic resistomes of diverse river habitats.

Wastewater treatment plants (WWTPs) are key sources of antimicrobial resistance genes (ARGs) that could influence the resistomes of microbial communities in various habitats of the receiving river ecosystem. However, it is currently unknown which habitats are most impacted and whether ARGs, like certain chemical contaminants, could be accumulated or enriched in the river ecosystem. We conducted a systematic metagenomic survey on the antibiotic resistomes of WWTP effluent, four riverine habitats (water, suspended particles, sediment, epilithic biofilm), and freshwater amphipod gut microbiomes. The impact of WWTP effluent on the downstream habitats was assessed in nine Swiss rivers. While there were significant differences in resistomes across habitats, the wastewater resistome was more similar to the resistome of receiving river water than to the resistomes of other habitats, and river water was the habitat most strongly impacted by the WWTPs effluent. The sulfonamide, beta-lactam, and aminoglycoside resistance genes were among the most abundant ARGs in the WWTP effluents, and especially aadA, sul1, and class A beta-lactamase genes showed significantly increased abundance in the river water of downstream compared to upstream locations (p < 0.05). However, this was not the case for the sediment, biofilm, and amphipod gut habitats. Accordingly, evidence for accumulation or enrichment of ARGs through the riverine food web was not identified. Our study suggests that monitoring riverine antimicrobial resistance determinants could be conducted using "co-occurrence" of aadA, sul1, and class A beta-lactamase genes as an indicator of wastewater-related pollution and should focus on the water as the most affected habitat.

First report of lettuce ring necrosis virus in chili pepper and tomato in Belgium and The Netherlands.

Lettuce ring necrosis virus (LRNV), genus Ophiovirus, was detected by the Netherlands Institute for Vectors, Invasive plants and Plant health (NIVIP) in June and November of 2021 in two samples of chili pepper fruits (Capsicum spp.), both in mixed infection with other viruses. The first sample originated from a production site in Belgium (Sample ID: 40009704) and the second from a production site in the Netherlands (Sample ID: 41115269). One of the fruits of 40009704 showed a light purple circular pattern, while fruits from 41115269 showed colored (ring)spots. The samples were analyzed using Illumina sequencing on a NovaSeq 6000 platform (PE 150) as described previously (Hammond et al., 2021), obtaining 39.9M and 22.8M total reads for 40009704 and 41115269. The corresponding sequence read archives (SRA) were deposited in the NCBI SRA database under BioProject accession number PRJNA917231. From both samples, the nearly complete genome of LRNV (RNA1-4) was obtained and deposited in GenBank (40009704, OQ160823- OQ160826 (7616, 1799, 1502, 1382 nt, mapped reads: 40K, 12K, 114K, 12K , average read coverage (ARC): 0.8K, 0.9K, 11.3K and 1.1K); 41115269, OQ160827- OQ160830 (7616, 1801, 1518, 1389 nt, mapped reads: 112K, 7K, 357K, 55K reads, ARC: 2.2K, 0.6K, 34K and 5.8K)). The shared sequence identities with the Genbank reference sequence of LRNV (NC_006051-NC_006051) were 99.2 and 99.2% (RNA1), 99.1 and 99.1% (RNA2), 98.3 and 98.8% (RNA3), 99.0 and 98.9% (RNA4) for 40009704 and 41115269 respectively. The shared sequence identities between 40009704 and 41115269 were 99.9 (RNA1), 99.0 (RNA2), 99.1 (RNA3) and 99.5% (RNA4). In addition to LRNV, the ophiovirus ranunculus white mottle virus (RWMV) was detected in both samples (OQ160831-OQ160834; OQ160835-OQ160838), while the tobamovirus pepper mild mottle virus (PMMoV) was present in the fruits of 41115269 (OQ160839). Since RWMV has been associated with leaf symptoms in pepper (Gambley et al., 2019; Rivarez et al., 2022) and the colored (ring)spots of 41115269 were very similar to reported symptoms of PMMoV-infected pepper fruits (Martínez-Ochoa et al., 2003), it remains unclear whether LRNV contributed to the observed symptoms. Additionally, LRNV was detected in tomato (Solanum lycopersicum) in Belgium in 2020. In the frame of a metagenomic survey using Virion-Associated Nucleic Acids (VANA)-based protocol (Maclot et al., 2021) on a Nextseq 500 platform (PE 150), partial genome sequences of LRNV were detected in two pools of tomato plants. One pool was made of 44 asymptomatic cultivars from a non-commercial grower (one sample per cultivar) yielding 118K total reads of which 84, 59, 335, and 18 reads mapped on RNA1, 2, 3, and 4, covering 35%, 69%, 100% and 55% of the genome, respectively. The other pool consisted of 15 plants from one cultivar from a production site yielding 3.1M total reads of which 6 and 5 reads mapped on RNA3 and 4, respectively. The detection of LRNV was confirmed for both pooled samples using the real-time RT-PCR method, targeting the CP gene, as described by Maachi et al. (2021). To our knowledge this is the first report of LRNV in pepper anywhere in the world. Additionally, although the disease lettuce ring necrosis in lettuce (Lactuca sativa) has been described in Belgium and the Netherlands before the causal agent was identified (Bos &amp; Huijberts, 1996), this is the first official report of this virus in Belgium and the Netherlands. This publication resulted from pre-publication data sharing of sequences and biological data among plant virologists to provide more context to two independent findings (Hammond et al., 2021).