Corrigendum: Evaluation of the rbcL marker for metabarcoding of marine diatoms and inference of population structure of selected genera.

[This corrects the article DOI: 10.3389/fmicb.2023.1071379.].

Evaluation of the rbcL marker for metabarcoding of marine diatoms and inference of population structure of selected genera.

Diatoms are one of the most important phytoplankton groups in the world's oceans. There are responsible for up to 40% of the photosynthetic activity in the Ocean, and they play an important role in the silicon and carbon cycles by decoupling carbon from atmospheric interactions through sinking and export. These processes are strongly influenced by the taxonomic composition of diatom assemblages. Traditionally, these have been assessed using microscopy, which in some cases is not reliable or reproducible. Next-generation sequencing enabled us to study diversity in a high-throughput manner and uncover new distribution patterns and diversity. However, phylogenetic markers used for this purpose, such as various 18S rDNA regions, are often insufficient because they cannot distinguish between some taxa. In this work, we demonstrate the performance of the chloroplast-encoded rbcL marker for metabarcoding marine diatoms compared to microscopy and 18S-V9 metabarcoding using a series of monthly samples from the Gulf of Trieste (GoT), northern Adriatic Sea. We demonstrate that rbcL is able to detect more taxa compared to 18S-V9 metabarcoding or microscopy, while the overall structure of the diatom assemblage was comparable to the other two methods with some variations, that were taxon dependent. In total, 6 new genera and 22 new diatom species for the study region were identified. We were able to spot misidentification of genera obtained with microscopy such as Pseudo-nitzschia galaxiae, which was mistaken for Cylindrotheca closterium, as well as genera that were completely overlooked, such as Minidiscus and several genera from the Cymatosiraceae family. Furthermore, on the example of two well-studied genera in the region, namely Chaetoceros and particularly Pseudo-nitzschia, we show how the rbcL method can be used to infer even deeper phylogenetic and ecologically significant differences at the species population level. Despite a very thorough community analysis obtained by rbcL the incompleteness of reference databases was still evident, and we shed light on possible improvements. Our work has further implications for studies dealing with taxa distribution and population structure, as well as carbon and silica flux models and networks.

Amplicon sequence variant-based meiofaunal community composition revealed by DADA2 tool is compatible with species composition.

The present study is aimed at implementing the morphological identification-free amplicon sequence variant (ASV) concept for describing meiofaunal species composition, while strongly indicating reasonable compatibility with the underlying species. A primer pair was constructed and demonstrated to PCR amplify a 470-490 bp 18S barcode from a variety of meiofaunal taxa, high throughput sequenced using the Illumina 300 × 2 bps platform. Sixteen 18S multi-species HTS assemblies were created from meiofaunal samples and merged to one assembly of ~2,150,000 reads. Five quality scores (q = 35, 30, 25, 20, 15) were implemented to filter five 18S barcode assemblies, which served as inputs for the DADA2 software, ending with five reference ASV libraries. Each of these libraries was clustered, applying 3% dissimilarity threshold, revealed an average number of 1.38 ±â€¯0.078 ASVs / cluster. Hence, demonstrating high level of ASV uniqueness. The libraries which were based on q ≤ 25 reached a near-asymptote number of ASVs which together with the low average number of ASVs / cluster, strongly indicated fair representation of the actual number of the underlying species. Hence, the q = 25 library was selected to be used as metabarcoding reference library. It contained 461 ASVs and 342-3% clusters with average number of 1.34 ±â€¯1.036 ASV / cluster and their BLASTN annotation elucidated a variety of expected meiofaunal taxa. The sixteen assemblies of sample-specific paired reads were mapped to this reference library and sample ASV profiles, namely the list of ASVs and their proportional copy numbers were created and clustered.

DNA metabarcoding and morphological analysis - Assessment of zooplankton biodiversity in transitional waters.

Zooplankton biodiversity assessment is a crucial element in monitoring marine ecosystem processes and community responses to environmental alterations. In order to evaluate the suitability of metabarcoding for zooplankton biodiversity assessment and biomonitoring as a fast and more cost-effective method, seasonal zooplankton sampling was carried out in the Venice Lagoon and the nearby coastal area (Northern Adriatic Sea). The molecular analysis showed higher taxa richness compared to the classical morphological method (224 vs. 88 taxa), discriminating better the meroplanktonic component, morphologically identified only up to order level. Both methods revealed a similar spatio-temporal distribution pattern and the sequence abundances and individual counts were significantly correlated for various taxonomic groups. These results indicate that DNA metabarcoding is an efficient tool for biodiversity assessments in ecosystems with high spatial and temporal variability, where high sampling effort is required as well as fast alert systems for non-native species (NIS).

Environmental DNA assessment of airborne plant and fungal seasonal diversity.

Environmental DNA (eDNA) metabarcoding and metagenomics analyses can improve taxonomic resolution in biodiversity studies. Only recently, these techniques have been applied in aerobiology, to target bacteria, fungi and plants in airborne samples. Here, we present a nine-month aerobiological study applying eDNA metabarcoding in which we analyzed simultaneously airborne diversity and variation of fungi and plants across five locations in North and Central Italy. We correlated species composition with the ecological characteristics of the sites and the seasons. The most abundant taxa among all sites and seasons were the fungal genera Cladosporium, Alternaria, and Epicoccum and the plant genera Brassica, Corylus, Cupressus and Linum, the latter being much more variable among sites. PERMANOVA and indicator species analyses showed that the plant diversity from air samples is significantly correlated with seasons, while that of fungi varied according to the interaction between seasons and sites. The results consolidate the performance of a new eDNA metabarcoding pipeline for the simultaneous amplification and analysis of airborne plant and fungal particles. They also highlight the promising complementarity of this approach with more traditional biomonitoring frameworks and routine reports of air quality provided by environmental agencies.

Ecological time series and integrative taxonomy unveil seasonality and diversity of the toxic diatom Pseudo-nitzschia H. Peragallo in the northern Adriatic Sea.

Pseudo-nitzschia H. Peragallo (1900) is a globally distributed genus of pennate diatoms that are important components of phytoplankton communities worldwide. Some members of the genus produce the neurotoxin domoic acid, so regular monitoring is in place. However, the identification of toxic members in routine samplings remains problematic. In this study, the diversity and seasonal occurrence of Pseudo-nitzschia species were investigated in the Gulf of Trieste, a shallow gulf in the northern Adriatic Sea. We used time series data from 2005 to 2018 to describe the seasonal and inter-annual occurrence of the genus in the area and its contribution to the phytoplankton community. On average, the genus accounted for about 15 % of total diatom abundance and peaked in spring and autumn, with occasional outbreaks during summer and large inter-annual fluctuations. Increased water temperature and decreased salinity positively affected the presence of some members of the genus, while strong effects could be masked by an unsuitable definition of the species complexes used for monitoring purposes. Therefore, combining morphological (TEM) and molecular analyses by sequencing the ITS, 28S and rbcL markers, eight species were identified from 83 isolated monoclonal strains: P. calliantha, P. fraudulenta, P. delicatissima, P. galaxiae, P. mannii, P. multistriata, P. pungens and P. subfraudulenta. A genetic comparison between the isolated strains and other strains in the Mediterranean was carried out and rbcL was inspected as a potential barcode marker in respect to our results. This is the first study in the Gulf of Trieste on Pseudo-nitzschia time series from a long-term ecological research (LTER) site coupled with molecular data. We show that meaningful ecological conclusions can be drawn by applying integrative methodology, as opposed to the approach that only considers species complexes. The results of this work will provide guidance for further monitoring efforts as well as research activities, including population genetics and genomics, associated with seasonal distribution and toxicity profiles.

PLANiTS: a curated sequence reference dataset for plant ITS DNA metabarcoding.

DNA metabarcoding combines DNA barcoding with high-throughput sequencing to identify different taxa within environmental communities. The ITS has already been proposed and widely used as universal barcode marker for plants, but a comprehensive, updated and accurate reference dataset of plant ITS sequences has not been available so far. Here, we constructed reference datasets of Viridiplantae ITS1, ITS2 and entire ITS sequences including both Chlorophyta and Streptophyta. The sequences were retrieved from NCBI, and the ITS region was extracted. The sequences underwent identity check to remove misidentified records and were clustered at 99% identity to reduce redundancy and computational effort. For this step, we developed a script called 'better clustering for QIIME' (bc4q) to ensure that the representative sequences are chosen according to the composition of the cluster at a different taxonomic level. The three datasets obtained with the bc4q script are PLANiTS1 (100 224 sequences), PLANiTS2 (96 771 sequences) and PLANiTS (97 550 sequences), and all are pre-formatted for QIIME, being this the most used bioinformatic pipeline for metabarcoding analysis. Being curated and updated reference databases, PLANiTS1, PLANiTS2 and PLANiTS are proposed as a reliable, pivotal first step for a general standardization of plant DNA metabarcoding studies. The bc4q script is presented as a new tool useful in each research dealing with sequences clustering. Database URL: https://github.com/apallavicini/bc4q; https://github.com/apallavicini/PLANiTS.

Multi-marker metabarcoding approach to study mesozooplankton at basin scale.

Zooplankton plays a pivotal role in marine ecosystems and the characterisation of its biodiversity still represents a challenge for marine ecologists. In this study, mesozooplankton composition from 46 samples collected in summer along the western Adriatic Sea, was retrieved by DNA metabarcoding analysis. For the first time, the highly variable fragments of the mtDNA COI and the V9 region of 18S rRNA genes were used in a combined matrix to compile an inventory of mesozooplankton at basin scale. The number of sequences retrieved after quality filtering were 824,148 and 223,273 for COI and 18S (V9), respectively. The taxonomical assignment against reference sequences, using 95% (for COI) and 97% (for 18S) similarity thresholds, recovered 234 taxa. NMDS plots and cluster analysis divided coastal from offshore samples and the most representative species of these clusters were distributed according to the dominant surface current pattern of the Adriatic for the summer period. For selected sampling sites, mesozooplankton species were also identified under a stereo microscope providing insights on the strength and weakness of the two approaches. In addition, DNA metabarcoding was shown to be helpful for the monitoring of non-indigenous marine metazoans and spawning areas of commercial fish species. We defined pros and cons of applying this approach at basin scale and the benefits of combining the datasets from two genetic markers.

DNA metabarcoding uncovers fungal diversity of mixed airborne samples in Italy.

Fungal spores and mycelium fragments are particles which become and remain airborne and have been subjects of aerobiological studies. The presence and the abundance of taxa in aerobiological samples can be very variable and impaired by changeable climatic conditions. Because many fungi produce mycotoxins and both their mycelium fragments and spores are potential allergens, monitoring the presence of these taxa is of key importance. So far data on exposure and sensitization to fungal allergens are mainly based on the assessment of few, easily identifiable taxa and focused only on certain environments. The microscopic method used to analyze aerobiological samples and the inconspicuous fungal characters do not allow a in depth taxonomical identification. Here, we present a first assessment of fungal diversity from airborne samples using a DNA metabarcoding analysis. The nuclear ITS2 region was selected as barcode to catch fungal diversity in mixed airborne samples gathered during two weeks in four sites of North-Eastern and Central Italy. We assessed the taxonomic composition and diversity within and among the sampled sites and compared the molecular data with those obtained by traditional microscopy. The molecular analyses provide a tenfold more comprehensive determination of the taxa than the traditional morphological inspections. Our results prove that the metabarcoding analysis is a promising approach to increases quality and sensitivity of the aerobiological monitoring. The laboratory and bioinformatic workflow implemented here is now suitable for routine, high-throughput, regional analyses of airborne fungi.

ITS2 metabarcoding analysis complements lichen mycobiome diversity data.

Lichen thalli harbor complex fungal communities (mycobiomes) of species with divergent trophic and ecological strategies. The complexity and diversity of lichen mycobiomes are still largely unknown, despite surveys combining culture-based methods and high-throughput sequencing (HTS). The results of such surveys are strongly influenced by the barcode locus chosen, its sensitivity in discriminating taxa, and the depth to which public sequence repositories cover the phylogenetic spectrum of fungi. Here, we use HTS of the internal transcribed spacer 2 (ITS2) to assess the taxonomic composition and diversity of a well-characterized, alpine rock lichen community that includes thalli symptomatically infected by lichenicolous fungi as well as asymptomatic thalli. Taxa belonging to the order Chaetothyriales are the major components of the observed lichen mycobiomes. We predict sequences representative of lichenicolous fungi characterized morphologically and assess their asymptomatic presence in lichen thalli. We demonstrated the limitations of metabarcoding in fungi and show how the estimation of species diversity widely differs when ITS1 or ITS2 are used as barcode, and particularly biases the detection of Basidiomycota. The complementary analysis of both ITS1 and ITS2 loci is therefore required to reliably estimate the diversity of lichen mycobiomes.

Environmental DNA in subterranean biology: range extension and taxonomic implications for Proteus.

Europe's obligate cave-dwelling amphibian Proteus anguinus inhabits subterranean waters of the north-western Balkan Peninsula. Because only fragments of its habitat are accessible to humans, this endangered salamander's exact distribution has been difficult to establish. Here we introduce a quantitative real time polymerase chain reaction-based environmental DNA (eDNA) approach to detect the presence of Proteus using water samples collected from karst springs, wells or caves. In a survey conducted along the southern limit of its known range, we established a likely presence of Proteus at seven new sites, extending its range to Montenegro. Next, using specific molecular probes to discriminate the rare black morph of Proteus from the closely related white morph, we detected its eDNA at five new sites, thus more than doubling the known number of sites. In one of these we found both black and white Proteus eDNA together. This finding suggests that the two morphs may live in contact with each other in the same body of groundwater and that they may be reproductively isolated species. Our results show that the eDNA approach is suitable and efficient in addressing questions in biogeography, evolution, taxonomy and conservation of the cryptic subterranean fauna.

Subclinical alteration of the cervical-vaginal microbiome in women with idiopathic infertility.

Biomarkers have a wide application in research and clinic, they help to choose the correct treatment for diseases. Recent studies, addressing the vaginal microbiome using next generation sequencing (NGS), reported the involvement of bacterial species in infertility. We compared the vaginal microbiome of idiopathic infertile women with that of healthy, including bacterial vaginosis affected women and non-idiopathic infertile women, to identify bacterial species suitable as biomarkers. Information on microorganisms was obtained from the V3-16S rDNA sequencing of cervical-vaginal fluids of 96 women using the Ion Torrent platform. Data were processed with QIIME and classified against the Vaginal 16S rDNA Reference Database. The analysis revealed a significant beta-diversity variation (p < 0.001) between the four groups included in the study. L. iners, L. crispatus, and L. gasseri distinguished idiopathic infertile women from the other groups. In these women, a microbial profile similar to that observed in bacterial vaginosis women has been detected. Our results suggest that the quantitative assessment and identification of specific microorganisms of the cervical-vaginal microflora could increase the accuracy of available tools for the diagnosis of infertility and improve the adoption of therapeutic protocols.

The Pillars of Hercules as a bathymetric barrier to gene flow promoting isolation in a global deep-sea shark (Centroscymnus coelolepis).

Knowledge of the mechanisms limiting connectivity and gene flow in deep-sea ecosystems is scarce, especially for deep-sea sharks. The Portuguese dogfish (Centroscymnus coelolepis) is a globally distributed and near threatened deep-sea shark. C. coelolepis population structure was studied using 11 nuclear microsatellite markers and a 497-bp fragment from the mtDNA control region. High levels of genetic homogeneity across the Atlantic (Φ(ST) = -0.0091, F(ST) = 0.0024, P > 0.05) were found suggesting one large population unit at this basin. The low levels of genetic divergence between Atlantic and Australia (Φ(ST) = 0.0744, P < 0.01; F(ST) = 0.0015, P > 0.05) further suggested that this species may be able to maintain some degree of genetic connectivity even across ocean basins. In contrast, sharks from the Mediterranean Sea exhibited marked genetic differentiation from all other localities studied (Φ(ST) = 0.3808, F(ST) = 0.1149, P < 0.001). This finding suggests that the shallow depth of the Strait of Gibraltar acts as a barrier to dispersal and that isolation and genetic drift may have had an important role shaping the Mediterranean shark population over time. Analyses of life history traits allowed the direct comparison among regions providing a complete characterization of this shark's populations. Sharks from the Mediterranean had markedly smaller adult body size and size at maturity compared to Atlantic and Pacific individuals. Together, these results suggest the existence of an isolated and unique population of C. coelolepis inhabiting the Mediterranean that most likely became separated from the Atlantic in the late Pleistocene.