Front. Microbiol. Frontiers in Microbiology Front. Microbiol. 1664-302X Frontiers Media S.A. 10.3389/fmicb.2020.593401 Microbiology Original Research The Microbiome of Leonardo da Vinci’s Drawings: A Bio-Archive of Their History Piñar Guadalupe 1 * Sclocchi Maria Carla 2 Pinzari Flavia 3 Colaizzi Piero 2 Graf Alexandra 4 Sebastiani Maria Letizia 2 Sterflinger Katja 1 1Department of Biotechnology, Institute of Microbiology and Microbial Biotechnology, University of Natural Resources and Life Sciences (BOKU), Vienna, Austria 2Laboratorio di Biologia, Istituto Centrale per la Patologia degli Archivi e del Libro (ICPAL), Rome, Italy 3Institute for Biological Systems (ISB), Council of National Research of Italy (CNR), Monterotondo, Italy 4Applied Life Sciences/Bioengineering/Bioinformatics, FH Campus, Vienna, Austria

Edited by: Massimiliano Marvasi, University of Florence, Italy

Reviewed by: Jay Prakash Verma, Banaras Hindu University, India; Giovanna Cristina Varese, University of Turin, Italy

*Correspondence: Guadalupe Piñar, guadalupe.pinar@boku.ac.at

Present address: Katja Sterflinger, Institute of Natural Sciences and Technology in the Arts, Academy of Fine Arts Vienna, Vienna, Austria

This article was submitted to Microbiotechnology, a section of the journal Frontiers in Microbiology

20 11 2020 2020 11 593401 10 08 2020 26 10 2020 Copyright © 2020 Piñar, Sclocchi, Pinzari, Colaizzi, Graf, Sebastiani and Sterflinger. 2020 Piñar, Sclocchi, Pinzari, Colaizzi, Graf, Sebastiani and Sterflinger

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

Seven emblematic Leonardo da Vinci’s drawings were investigated through third generation sequencing technology (Nanopore). In addition, SEM analyses were carried out to acquire photographic documentation and to infer the nature of the micro-objects removed from the surface of the drawings. The Nanopore generated microbiomes can be used as a “bio-archive” of the drawings, offering a kind of fingerprint for current and future biological comparisons. This information might help to create a biological catalog of the drawings (cataloging), a microbiome-fingerprint for each single analyzed drawing, as a reference dataset for future studies (monitoring) and last but not least a bio-archive of the history of each single object (added value). Results showed a relatively high contamination with human DNA and a surprising dominance of bacteria over fungi. However, it was possible to identify typical bacteria of the human microbiome, which are mere contaminants introduced by handling of the drawings as well as other microorganisms that seem to have been introduced through vectors, such as insects and their droppings, visible through the SEM analyses. All drawings showed very specific bio-archives, but a core microbiome of bacteria and fungi that are repeatedly found in this type of material as true degraders were identified, such as members of the phyla Proteobacteria, Actinobacteria, and Firmicutes among bacteria, and fungi belonging to the classes Sordariomycetes and Eurotiomycetes. In addition, some similarities were observed that could be influenced by their geographical location (Rome or Turin), indicating the influence of this factor and denoting the importance of environmental and storage conditions on the specific microbiomes.

third generation sequencing nanopore technology biological diagnosis paper material microbiome Leonardo da Vinci insect droppings bio-archive

香京julia种子在线播放

    1. <form id=HxFbUHhlv><nobr id=HxFbUHhlv></nobr></form>
      <address id=HxFbUHhlv><nobr id=HxFbUHhlv><nobr id=HxFbUHhlv></nobr></nobr></address>

      Introduction

      Our cultural heritage is an asset that we must protect and monitor in order to pass it on to our future generations. One of the greatest pleasures is its contemplation in various permanent or temporary exhibitions. Most valuable graphic art pieces are displayed or stored in a particular museum, library or exhibition site under optimal conditions for their conservation and safety. However, it is common for these objects to be loaned out for certain thematic exhibitions in other museums or exhibition halls around the world. In this context, it is necessary to devise a method to monitor the risks of microbial contamination that may occur during packaging, transport and/or exhibition at the other sites and to allow a comparison before and after these trips. This monitoring is carried out by restorers and experts from the various institutions that provide and receive the object in question, but also by experts from renowned third-party institutions, which apply conventional technologies such as visual inspection, microscopy, but also introducing molecular technologies in order to achieve a detailed assessment of the risk of biological contamination posed by this type of transfer.

      Conservators have long since started using new technologies such as miniaturized environmental condition recorders that keep track of how an object is stored during a move or an exhibition, away from the usual storage environment (Valentini et al., 2018). In addition to temperature and humidity monitoring systems, modern methods of classification and identification have also been developed, such as smart labels and barcodes technologies, which make it possible to quickly catalog objects and are already widely used in libraries. For example, the Radio-Frequency Identification systems (Steinberg et al., 2014) leverage radio waves to transmit data from chips to the readers (Alwadi et al., 2017). These systems help to cover traditional work processes such as check-in, check-out, anti-theft control and management of inventories. Even more innovative are the labels based on synthetic DNA (Demidov, 2020). These are patented systems that make it possible to make the art pieces not forgeable. Artists can authenticate their work by labeling it with a univocal synthetic bioengineered DNA sequence that provides an encrypted link between the artwork and a secure database containing the definitive information about the artwork. DNA data can be retrieved by sequencing or using other technologies like fluorescent probes (Banerjee et al., 2016). However, any object already contains DNA, from the people who have handled the object since its creation and from all the organisms and microorganisms that have come into contact with it, also through the deposition of dust and dirt. This natural bio-archive changes with time and can give certain information on some particular events, such as the presence on the materials of deteriorating microorganisms harmful for conservation, but also the permanence of the object in particular geographical locations, as indicated by the molecular traces left by the pollen of endemic plants (Piñar et al., 2019).

      Molecular techniques, including DNA-based fingerprinting techniques in addition to cloning and Sanger sequencing, have been largely applied to investigate and monitor the biological colonization on art objects since more than two decades (Schabereiter-Gurtner et al., 2001). This molecular strategy, together with imaging and chemical analyses, has delivered relevant information about the microbial communities associated with different materials and the microbiological risk they represented on their surface (Sterflinger and Piñar, 2013). Nevertheless, the use of first-generation sequencing technology delivers limited data and is time consuming as well as not free of bias. This strategy has become obsolete and has been replaced by new high-throughput molecular technologies in which massive DNA sequencing is possible, offering a much broader view of the real communities that colonize our cultural heritage. Next Generation Sequencing (NGS) technologies have evolved incorporating revolutionary improvements to address the complexities of genomes and metagenomes at an unprecedented speed (Goodwin et al., 2016), being now well established in the field of cultural heritage. These analyses can be performed either by the so-called “shotgun metagenomic approach,” by sequencing the entire DNA library representing all the molecular components of a given sample (Teasdale et al., 2017; Piñar et al., 2020b), or by focusing on specific conserved sequences such as ribosomal RNA genes (Marvasi et al., 2019). The latter approach called “target amplification approach” has some advantages, such as the reduction of data complexity, the possibility to assign more sequences to specific organisms and finally, facilitates some semi-quantitative analysis. This approach is extremely useful in the field of cultural heritage, as DNA with a poor quality or low concentration can be amplified by using degenerate primers and PCR. Many of the studies using the latter strategy have been collected and summarized in the review of Marvasi et al. (2019). These methods are powerful, but as in any PCR-based method, primers and/or exponential amplification introduce bias. In addition, the simultaneous quantitative assessment of all three domains of life is impossible. Furthermore, the amplification of target sequences (mainly ribosomal sequences) do not allow us to draw conclusions on any functional genes and therefore the microbial functional traits remain uncertain. Nevertheless, the further technological development of molecular methods is overcoming many technical biases. Third generation sequencing technologies are emerging and offering some advantages over the limitations suffered by NGS platforms, such as the short-read length that needs to be assembled with the help of bioinformatic tools into the original length template, and the PCR bias introduced by clonal amplification for the detection of base incorporation signal.

      Third generation sequencing technologies offer an intelligent solution to the limitations mentioned above (Schadt et al., 2010). Instead of sequencing a clonally amplified template, a single DNA template is sequenced, leading to a reduced use of reagents and chemicals, simplified protocols, and miniaturization of the entire technical process and equipment. One example is the Nanopore sequencing technology, which involves the use of protein nanopores embedded in a polymer membrane. The sequencing does not necessarily require any intervening PCR amplification or chemical labeling step. Furthermore, this technology offers a pocket-sized sensing device, the MinION, which offers the advantages of a portable device. To date, the Nanopore sequencing technology has been scarcely applied in the field of cultural heritage, and only four studies have reported on the advantages of applying this state-of-the-art technology on valuable art objects. Three of these studies (Šoltys et al., 2019; Grottoli et al., 2020; Kisová et al., 2020) reported on the advantages of sequencing long DNA fragments when targeting ribosomal regions. However, all three studies were based on the target amplification approach that avoided inferring the relative proportion of the different phylogenetic groups. The fourth study (Piñar et al., 2020a) applied for the first time the Nanopore sequencing technology together with a whole genome amplification (WGA) protocol to perform a rapid diagnosis of the biological infection on objects of art, unifying all the advantages that this new technology can offer for metagenomic analyses starting from very low DNA concentrations and reflecting the real proportions of all domains of life.

      Although molecular techniques are not completely free from bias (Sterflinger et al., 2018a), the establishment of next- and third-generation sequencing technologies in the field of cultural heritage is providing a powerful DNA-based approach to monitor the preservation status as well as to foresee the risk of deterioration of valuable objects (Piñar et al., 2020a). The data generated can provide interesting information that can contribute to surprising insights into the objects being investigated, such as the selection of materials at the time of manufacture and the conditions of their storage or about possible displacements and geographical origins, but also uncover important information about the object’s history of use (Teasdale et al., 2017; Piñar et al., 2019, 2020b). This information may help to understand many open questions in a variety of fields, as archeology, history, restoration, philology, criminology but also, they can contribute to giving an historical benefit to the investigated objects in form of a bio-archive (Piñar et al., 2019, 2020b).

      In the present study we apply the Nanopore sequencing technology in combination with a whole genome amplification (WGA) protocol to survey some of Leonardo da Vinci’s most emblematic drawings. This strategy was applied as the only diagnostic method or in combination with SEM analyses, in order to explore the nature of contaminants deposited on the surfaces of the drawings. The sampling strategy, based on a very delicate micro-aspiration, was aimed primarily at collecting dust particles, microbial cells and other debris from small areas of surfaces. The aim was to create a bio-archive of the dust and other collectable material deposited on the drawings, in order to monitor their conservation status and thus be able to use them for possible comparisons in the future. The sampling was not used to assess the percentage of living propagules, however, the molecular strategy adopted could theoretically indicate the presence of fungal or bacterial cells with an anomalous abundance, compared to typical indoor dust composition. Furthermore, the information contained in this bio-archive is an added value for the drawings, that can in turn help to understand some issues related to the history of storage of the investigated drawings. This study represents one more example of how the Nanopore technology, in this case in combination with microscopy techniques, can be a practical tool for the rapid biological diagnosis and monitoring in priceless objects.

      Materials and Methods Objects and Sampling

      Seven emblematic drawings from Leonardo da Vinci were selected for the analysis of their microbiomes (Figure 1). Five of these drawings are currently housed in the Royal Library of Turin, located on the ground floor of the Royal Palace in Turin (Italy): “Autoritratto” (sample L2, 33.3 × 21.3 cm), “Nudi per la battaglia di Anghiari” (sample L3, 25.4 × 19.7 cm), “Studi delle gambe anteriori di un cavallo” (sample L4, 15.4 × 20.5 cm), “Studi di insetti” (sample L5, 12.9 × 11.8 cm) and “Studi di gambe virilirecto, “Figura presso il fuocoverso (sample L6, 14.0 × 15.7 cm). The last two investigated drawings are stored at the Corsinian Library in Rome (Italy): “Uomo della Bitta” (sample L7, 22 × 14.6 cm) and “Studio di panneggio per una figura inginocchiata” (sample L8, 25.7 × 19.0 cm).

      The Leonardo da Vinci’s drawings investigated. The drawings currently housed in the Royal Library of Turin are: “Autoritratto” (sample L2, 33.3 × 21.3 cm), “Nudi per la battaglia di Anghiari” (sample L3, 25.4 × 19.7 cm), “Studi delle gambe anteriori di un cavallo” (sample L4, 15.4 × 20.5 cm), “Studi di insetti” (sample L5, 12.9 × 11.8 cm) and “Studi di gambe virili” recto, “Figura presso il fuoco” verso (sample L6, 14.0 × 15.7 cm). The drawings stored at the Corsinian Library in Rome are: “Uomo della Bitta” (sample L7, 22 × 14.6 cm) and “Studio di panneggio per una figura inginocchiata” (sample L8, 25.7 × 19.0 cm).

      Conservators of the ICPAL (Istituto Centrale per la Patologia degli Archivi e del Libro) collected samples from the recto and the verso of the selected drawings. A micro-sampling method for fragile surfaces has been specifically developed based on a protocol tested mainly on paintings and described by de Carvalho et al. (2019). The system is based on a micro-aspiration aimed primarily at collecting dust particles, and microbial cells from small areas of paper or parchment surfaces (Supplementary Figure 1). The suction force applied was modulable, and the filter membranes used to trap the particulates were sterile and made of a material suitable for subsequent use in both molecular and SEM analysis. The micro-aspiration system was also used to sample the microscopic objects detached from the drawings with needles under an optical microscope, as in the case of single fibers or particles of a few hundred microns. The suction system was an autoclavable polypropylene filter holder, 13 mm diameter (Swinnex filter holder, Millipore) with sterile membranes (Millipore cellulose nitrate membrane filter Ø = 13 mm, Pore size = 0.45 μm) connected to an oil-free diaphragm vacuum pump (DA14 Charles Austen MB2238. Appleton Woods Ltd., United Kingdom).

      The total number of membranes used for sampling varied between the drawings, ranging from 10 to 36 membranes per drawing (depending on the size of the drawing). Each single used membrane was immediately introduced in a single sterile Eppendorf tube and appropriately labeled for the storage of the samples until DNA extraction analyses. The same kind of membranes was used for SEM observation of the particulate sampled from the drawings and its chemical analysis.

      DNA Extraction

      The DNA was extracted using the FastDNA SPIN Kit for soil (MP Biomedicals, Illkrich, France) as recommended by the manufacturers. DNA extraction was performed directly from the membranes by grouping a maximum of four membranes per reaction tube. DNA extracted from each single reaction tube was pooled per sample to obtain a single microbiome from each drawing. The DNA concentrations were assessed by using the Qubit 2.0 fluorometer (Invitrogen Corporation), with the Qubit dsDNA HS Assay Kit.

      Whole Genome Amplification (WGA), Template Preparation and Sequencing

      Whole genome amplification (WGA) and template preparation was performed following the “Premium whole genome amplification protocol” available in the Oxford Nanopore community. For WGA, all reactions were executed in a BioRad C1000 Thermal Cycler using the REPLI-g Midi Kit (Qiagen), which uses the innovative Multiple Displacement Amplification (MDA) technology. The libraries were performed using the Ligation Sequencing kit 1D SQK-LSK109 and the Flow cell Priming kit EXP-FLP001 (Oxford Nanopore Technologies). All steps performed for library preparation were performed as described by Piñar et al. (2020a). A quality control of flow cells (SpotOn Flow cell Mk I R9 Version, FLO-MIN 106D) was performed prior starting the sequencing. The MinKNOWTM software was used to check the number of active pores in the flow cells (Table 1) and for performing the sequencing runs. All seven runs were performed for 48 h.

      Details of the sequencing runs performed on seven drawings of Leonardo da’Vinci using the MinIon (Nanopore sequencing technology).

      Drawing Sample (ID run) Active pores/flow cell Total readsa Analyzed readsb Classified readsc Total yield QSd Average length (Kb)
      Autoritratto L2 (193982) 1518 1,775,691 1,543,693 166,129 7.7 Gb 8.77 4.34
      Nudi per la battaglia di Anghiari L3 (193283) 1639 1,354,462 1,168,117 123,072 5.0 Gb 9.26 3.71
      Studi delle gambe anteriori di un cavallo L4 (193072) 1706 2,796,008 2,408,671 1,009,200 12.4 Gb 9.08 4.44
      Studi di insetti L5 (193067) 1558 2,857,065 2,461,722 317,761 13.9 Gb 9.09 4.86
      Studi di gambe virili”/“Figura presso il fuoco” L6 (192071) 1570 3,168,743 2,796,829 414,016 16.8 Gb 9.51 5.31
      Uomo della Bitta L7 (202641) 1512 1,710,990 1,514,255 65,343 6.8 Gb 10.94 3.97
      Studio di panneggio per una figura inginocchiata L8 (203053) 1523 2,732,413 2,376,913 183,811 14.2 Gb 10.38 5.19
      aTotal reads: Total reads generated in the sequencing run. bAnalyzed reads: Total number of reads processed by WIMP. cClassified reads: Number of reads classified phylogenetically. dQS, Quality Score.
      Data Analyses

      Resulting fast5 data files were basecalled using the Nanopore GPU basecalling with GUPPY 3.0.3 on UBUNTU 16.04 (Nanopore Community Platform). Once the Fastq files were generated, the data were compared with databases using one of the available pipelines for data analyses of the Nanopore Community Platform, following the steps recommended by manufacturers. The selected workflow chosen was “What’s in my pot” (WIMP), which is an EPI2ME workflow for taxonomic classification of basecalled sequences (reads) generated by Nanopore sequencing. WIMP initially filters FASTQ files with a mean q-score below a minimum threshold (defaults to 7).

      Relative abundances, taxonomic clustering, correspondence analysis and alpha diversity analysis were performed using the R packages taxonomic, phyloseq and microbiome for the R version 4.0.0. (McMurdie and Holmes, 2013). Before microbial taxonomic clustering was performed human reads where removed and relative abundances were calculated.

      Nucleotide Sequence Accession Number

      The sequences obtained by the metagenomic analysis are deposited and publicly available under the BioProject ID: PRJNA655381

      Optical Microscopy

      Features of interest in the drawings were documented with a Leica MZ16 stereoscopic microscope fitted with low-temperature optic fiber lights. The system was equipped with a digital camera connected with a computer and software (Leica Application Suite, LAS, Leica Microsystems GmbH Wetzlar, Germany) that allowed the acquisition of images at different magnifications.

      Scanning Electron Microscopy and Microanalysis

      The micro-objects removed from the drawings, fibers and membranes used to sample the drawings were analyzed using a variable pressure SEM instrument (EVO50, Carl-Zeiss Electron Microscopy Group) equipped with a detector for backscattered electron diffraction (BSD). After SEM imaging using Variable Pressure mode, some of the samples were covered with gold for further analysis in High Vacuum mode.

      The samples were placed on a sulfur-free carbon adhesive (Spectrotabs AGG3358, Agar Scientific) glued onto an aluminum stub (1.25 cm diameter, Agar Scientific). Chemical analysis was performed by means of EDS (INCA 250, Oxford Instruments). The SEM was fitted with a tungsten filament and operated at 20 keV, with an average working distance of 12.5 mm, with a chamber pressure between 30 and 150 Pa, chosen to keep unaltered micro-structures. The electron beam current was adjusted (∼400 pA) to generate sufficient X-rays to allow reliable identification of EDS peaks, minimize spectral artifacts and achieve detector dead times < 20%. The beam spot size was 5–30 nm in diameter. Acquisition live time for single analysis was set at 50 s. We performed EDS analyses of several types. Some were point analyses, and some consisted of the approximate entire surface of fibers and particles. Conventional ZAF correction (for atomic number Z, absorption and fluorescence) was applied, integrated into the Oxford INCA 250 microanalysis software used.

      Statistical Analysis of EDS Chemical Data

      The EDS measurements were analyzed using statistical tests to evaluate the relationships between variables and the significance of the differences between samples. One-way analysis of the variance (ANOVA) was applied and the significance of the differences was tested at 95% confidence. The ANOVA model used was “unbalanced” because the number of observations within each category was not the same. ANOVA was followed by a post-hoc analysis using Tukey’s t-test to determine the critical value for significance (McHugh, 2011). The measure of linear dependence (correlation) between couples of analyzed elements from the EDS datasets was obtained using the Pearson’s coefficient ”r. This coefficient has a value between +1 and −1, where +1 is total positive linear correlation, 0 is no linear correlation, and −1 is total negative linear correlation. The p-values that were computed for each coefficient allowed testing the null hypothesis that the coefficients were not significantly different from 0 (Reimann et al., 2008). ANOVA, PCA and Pearson’s correlation coefficient analyses were performed using XLSTAT 2017 software (Addinsoft, Paris, France).

      Results DNA Yield and Sequencing Analyses

      Seven outstanding drawings by Leonardo da Vinci were investigated (Figure 1) by using non-invasive sampling with suction membranes (Supplementary Figure 1), due to the strict policies and requirements for obtaining samples of the invaluable objects under investigation. The DNA was directly extracted from the membranes collected from each drawing and pooled, as explained in the materials and methods section. Results showed an expected low concentration of DNA, ranging from 5.4 to 20 pg/μl and for this reason, a protocol including a WGA, available in the Oxford Nanopore community was chosen (Piñar et al., 2020a), which enables to prepare DNA libraries with as little as 10 pg total DNA as input. Seven independent sequencing runs were carried out, each in an independent Nanopore flow cell loaded with the DNA library prepared from each of the drawings. Table 1 presents the details of the sequencing runs performed with the Nanopore MinION device, showing an average of 2,342,196 total reads generated per sequencing run. The reads analyzed and subsequently classified by the WIMP workflow showed an average of 2,038,600 and 325,619, respectively, indicating that about 16% of the analyzed reads were phylogenetically affiliated. In addition, the total yield obtained in the seven sequencing runs showed an average of 11 Gb, the average length of the fragments sequenced was 4.5 Kb, and the quality score (QS) showed an average of 9.6.

      Microbiomes Generated From Each Drawing

      The results showed that the proportion of the analyzed reads that could finally be phylogenetically assigned using the Nanopore WIMP workflow was different for each drawing (Table 1). However, there was an overall high biodiversity in the microbiome of the drawings investigated, revealing differences but also intriguing similarities between some of them. Drawings L7 and L8 showed the highest biodiversity, followed by drawings L5, L3, L6, and L2, all with similar alpha diversity, while L4 showed the lowest diversity of all (Figure 2). The percentage of human reads varied between 5.35% (L2) and 74.64% (L4) of the classified reads with a median of 11.36%. Therefore, we can assume that the high number of human reads recovered from the surface of drawing L4 masked the detection of reads related to microorganisms in this drawing. Figure 3 shows the distribution of the different superkingdoms among the investigated drawings. The relative abundance of bacteria ranged from 19% (with lowest proportion in drawing L4) to 79% (with highest proportion in drawing L6). Eukaryotes ranged from 20 to 80%, being the highest proportion present in drawing L4, once more mainly due to the massive presence of human DNA in this drawing. Viruses and Archaea represented less than 1% of the total microbiome in all investigated drawings.

      Alpha diversity of drawings calculated from the rarefied classified abundance data using the R package phyloseq.

      Relative abundance of eukaryotes, bacteria, archaea, and viruses in the microbiomes of all seven investigated drawings.

      Concerning bacteria (Figure 4A), the phylum Proteobacteria dominated in all samples (42–94% of bacteria), followed by members of Actinobacteria (3–28% of bacteria) and Firmicutes (1–21% of bacteria). The phylum Cyanobacteria showed to be represented only in the drawing L3 (6% of bacteria) while the phylum Planctomycetes was represented only in drawings L7 and L8 (2 and 1% of bacteria, respectively). The Eukaryotes (Figure 5) were mainly represented by the phylum Ascomycota (3–56% of eukaryotes) with the exception of the microbiome of drawing L4, which showed a massive dominance of the phylum Chordata (over 96% of the eukaryotes). This last phylum showed to be well represented in all drawings due to the presence of human DNA. Members of the phyla Basidiomycota, Mucoromycota, and Microsporidia were identified in all the drawings with relatively low proportions (0.5–5% of eukaryotes), in addition to some other phyla, each representing less than 0.5% of the microbiomes (marked as “others” in Figure 5).

      Relative abundance of the bacterial communities of all seven investigated drawings: (A) at the phylum level and (B) at the class level.

      Relative abundance of the eukaryotic communities of all seven investigated drawings at the phylum level.

      Each drawing owned a very specific microbiome, showing an independent molecular profile or “biological pedigree.” Supplementary Figure 2 shows all genera that represent more than 0.5% of the total microbiome in each of the drawings (cutoff 0.5%). These biological pedigrees are all summarized in Figure 6 in form of a heatmap, showing the relative abundance (%) of each genus in each of the drawings.

      Heatmap reporting the relative abundance (%) of genera in each of the investigated drawings; the genera shown represent greater than 0.5% abundance in at least one sample.

      Bacterial Communities

      The bacterial communities dominated over the eukaryotic communities in all the drawings. The most dominant phylum showed to be the Proteobacteria with dominance of the Gammaproteobacteria class in drawings 2–4 (Figure 4B). Inside this class, members of the order Pseudomonadales were detected, such as the genus Acinetobacter, present in drawings L2, L3, and L7, the genus Moraxella, which was highly abundant in L2 (representing over 51% of total reads) but also relative abundant in L3–L5 (1–5.25% total reads), the genus Psychrobacter being representative solely in L5 and the genus Pseudomonas, which was relative abundant in drawings L7 and L8 (Figure 6). Members of the Enterobacteriales were detected in drawings L2–L6 but absent in drawings L7–L8. The genus Salmonella was detected in L2 and L5 while Escherichia was identified in L4–L6. The genera Proteus, Xenorhabdus and Buchnera were present exclusively in L3 (Figure 6). Interestingly, the genus Alcanivorax (order Oceanospirillales) was present in drawings L7 and L8 and Stenotrophomonas (order Xanthomonadales) was only present in L8 with a high relative abundance (over 10% of total reads, Figure 6).

      In contrast, members of the Betaproteobacteria class were shown to be more abundant in drawings L5–L8 (Figure 4B) with Ralstonia as the dominant genus in all drawings (∼3–40% of total reads) (Figure 6). The genera Burkholderia and Cupriavidus were significative only in drawings L7 and L8 and Paraburkholderia exclusively in drawing L7. Members of the Alphaproteobacteria class were identified in all drawing but in relative lower proportions (Figure 4B). The genus Caulobacter showed to be exclusively present and with a relative high abundance (over 5% of total reads) in drawing L5. The genus Rhizobium was identified in L5 and L8 while Bradyrhizobium only in L8. The genus Sphingopyxis was shown to be present in L7 and L8, while Paracoccus was only identified in L7 and Pannonibacter only in L8 (Figure 6).

      The Actinobacteria were the second most represented phylum within the bacteria, with all members identified belonging to the Actinobacteria class (Figure 4B). Within the order Corynebacteriales, the genus Mycobacterium was present in drawings L3–L7 while the genus Nocardia was only identified in L7 and L8. Members of the order Micrococcales were identified in almost all the drawings, being the genus Intrasporangium only present in L2 and L7 and the genus Cryobacterium widely distributed among L3–L6. In addition, the genus Streptomyces (order Streptomycetales) was well represented in drawings L7 and L8 but also in lower proportion in L5 while Cutibacterium (order Propionibacteriales) was identified in L3 and L5 (Figure 6).

      The phylum Firmicutes was mainly represented by members of the Bacilli class (Figure 4B), with genera of the order Bacillales, such as Bacillus present exclusively in L3 and Staphylococcus, present in a relative high abundance in L3 (over 8% of total reads) but also well represented in L5 and L6 (4 and 3.15% of total reads, respectively) and in lower proportion in L8 (Figure 6). Members of the order Lactobacillales showed to be widely represented in almost all drawings, with the genera Streptococcus (in drawings L3–L6), Lactobacillus (in drawings L5 and L6), Lactococcus (in L4, L6, and L7) and Leuconostoc (solely in L6). Members of the Clostridia class were detected (Figure 4B), but each of the identified genera contributed less than 0.5% of the total reads.

      The same occurred with the Flavobacteria and Cytophagia classes within the phylum Bacteroidetes, no single genus accounted for more than 0.5% of the total reads (Figure 4B). The phylum Planctomycetes showed to be represented only in drawings L7 and L8 with members of the Planctomycetia class. Inside this class, only the genus Isosphaera accounted for more than 0.5% of the total reads in L7. Finally, the phylum Cyanobacteria was only significantly represented in L3 with the genera Anabaena and Oscillatoria, both accounting for more than 0.5% of total reads in this drawing (Figure 6).

      Eukaryotic Communities

      First of all, it is important to note that an average of 34% of the eukaryotic sequences that could be identified in the drawings showed to be affiliated with the Chordata phylum (Figure 5), being all related to human sequences, with the exception of drawing L4, where these sequences represented more than 96% of the eukaryotic sequences. The remaining eukaryotic sequences were affiliated with fungi distributed in different phyla, i.e., Ascomycota (3–56% of eukaryotes) and in lower proportions Basidiomycota, Mucoromycota and Microsporidia (0.5–5% of eukaryotes), in addition to some other phyla, each representing less than 0.5% of the microbiomes.

      Only a few taxa of the phylum Ascomycota represented more than 0.5% of the total reads in the microbiomes of the drawings (Figure 6). Thus, members of the Eurotiomycetes class proved to be widely distributed among the drawings, with species of the genus Aspergillus identified in L3, L5–L8, while species of the genera Penicillium (in L3) and Penicilliopsis (L7) were representative in only one of the drawings. Members of the Sordariomycetes class of the genus Colletotrichum showed to be representative in L2, L7, and L8. The genus Trichoderma was representative only in the last two drawings (L7 and L8) while the genus Grosmannia was solely identified in L5 in a relatively high proportion (12.6% of total reads) (Figure 6). Within the Dothideomycetes class, species of the genus Alternaria were representative in drawings L3 and L5 and of the genus Bipolaris in L7 and L8. Finally, members of the Saccharomycetes class showed to be widespread in drawings L2, L3, L5, and L6, with species of the genera Ascoidea and Candida (Figure 6).

      Canonical Correspondence Analyses

      Canonical correspondence analysis was used as an ordination method, providing a spatial representation of the relationship of taxa and samples. In the first CA axis of the bacterial community (Figure 7A), the drawing L2 showed a strong separation from the other drawings, indicating that the bacterial community in this drawing differs significantly from that detected in the others, mainly due to the high impact of the genus Moraxella on this sample. Looking at the fungal composition (Figure 7B), the drawing L5 clustered distinctly from the others in the first CA axis, confirming the high impact of the fungus Grosmannia on the fungal community. In the second CA axis, the separation between samples from Turin and from Rome can be seen within fungal as well as bacterial species, indicating a geographical influence on the microbiome of the drawings and denoting the importance of environmental and storage conditions on the specific microbiomes.

      Canonical correspondence analysis plot of the first and second CA axis for bacterial (A) and fungal (B) abundance data. Dots represent samples (sample name is indicated within the dots) and colors represent the city where the drawings are located—Turin (blue) and Rom (red). Percentages indicate how much of the variability in the data is explained by the respective axis.

      SEM Analyses

      For obvious reasons of relevance, fragility and uniqueness of the objects studied, the microscopic analysis of the drawings was mainly aimed at acquiring photographic documentation using non-invasive techniques, using above all a stereoscopic microscope with fiber optic illumination. The nature of micro-objects that the restorers were able to remove from the surface without causing damage could then be explored further by SEM-EDS. These were mostly particles superimposed on the fibers without any apparent historical value. In addition, some of the membranes used to remove dust from the surface of the designs undergoing molecular investigation were analyzed with SEM-EDS mainly to assess the effectiveness of the sampling method applied.

      On many of the drawings analyzed, roundish waxy brown incrustations were specifically observed, often surrounded by greasy looking halos on the underlying fibers. These recurrent “objects” have been documented and in some cases removed and chemically analyzed. Figure 8 shows the incrustations on L2 and L4 drawings, where they were particularly abundant, before their removal and their corresponding SEM images. The chemical analysis of three waxy incrustations removed from each of the two drawings is shown in Table 2, which shows the average values of a minimum of 10 EDS observations for each particle and the analysis of the statistical significance of the differences found between them. The EDS elemental analysis is not quantitative but rather qualitative. The values obtained are therefore mostly indicative. However, the comparison between the chemical composition of samples of similar nature and size, analyzed in the same way, may have a statistical value. Similarly, the correlation, and therefore the stoichiometric ratio, between pairs of chemical elements in each sample is significant because the technique measures the chemical composition of a small volume, in the order of nanometers, which therefore reflects the concentration of the elements and their relative percentage in that volume. In light of this, it is interesting what has emerged from the composition of the objects, which differ between the two drawings, and which contain magnesium phosphate in one case (L2 drawing) and potassium sulfate in the other (L4 drawing) (Supplementary Figure 3).

      Insects’ droppings. Waxy brown incrustations documented on the fibers of drawing L4 (Studi delle gambe anteriori di un cavallo) and drawing L2 (Autoritratto). These “objects” have been identified as insects’ droppings. They have been documented and, where possible, removed and analyzed by SEM-EDS. (a) an incrustation that was defacing the verso of drawing L4; (a1,a2) SEM images of the waxy material. The images were obtained in variable pressure, with a backscattered detector on non-metallized samples. (a1) shows the appearance of the surface of the object, covered with fungal mycelium, the detail in the square frame corresponds to (a3). The arrows in (a2) show round-shaped crystals embedded in the amorphous matrix. (b,b1,c,c1,d,d1) Correspond to incrustations that were present on the drawing L2. These were made both of an amorphous matrix, as in (b1), but also by spherulites, that consist of spherical crystals made of calcium and magnesium phosphate. (e,e1) Show an incrustation of the L4 drawing where in the amorphous matrix were documented columnar and bipyramidal crystals (arrow).

      EDS analysis of three insects’ droppings removed and analyzed from drawing L4 (“Studi delle gambe anteriori di un cavallo”) and drawing L2 (“Autoritratto”).

      C N O Na Mg P S K Ca Fe
      L4.1 75.41a 4.66bc 14.10b 0.70a 0.24c 0.65b 2.45a 1.32b 0.41b 0.06b
      L4.2 67.68b 8.58bc 15.89b 0.51b 0.48c 1.25b 2.50a 2.10a 0.83b 0.18ab
      L4.3 56.03c 23.96a 15.48b 0.38c 0.28c 0.62b 1.69b 1.12b 0.41b 0.04b
      L2.1 57.35c 0.00c 31.23a 0.08e 1.13a 3.12a 0.86c 0.13c 5.63a 0.48a
      L2.2 59.29c 0.00c 30.70a 0.20d 0.90b 3.31a 0.58c 0.00c 5.02a 0.00b
      L2.3 54.49c 12.63b 26.44a 0.02e 0.43c 1.05b 0.51c 0.16c 3.89a 0.38ab
      Pr > F (Model) <0.0001 <0.0001 <0.0001 <0.0001 <0.0001 <0.0001 <0.0001 <0.0001 <0.0001 0.02
      Significant Yes Yes Yes Yes Yes Yes Yes Yes Yes Yes
      Values are the means from a minimum of n = 15 replicates. In each column, different letter indicates statistically significant differences (P < 0.05), one-way ANOVA.

      SEM imaging and EDS analysis of the membranes used for sampling the airborne material deposited on the drawings between the fibers and destined to outline their genetic “pedigree,” showed the effectiveness of the suction system which, without contact with the drawing, made it possible to recover significant micro-objects including fungal spores. Figure 9 shows SEM image of a membrane used to sample the L3 drawing: it is possible to distinguish different objects trapped in the polymeric granular structure of the membrane. The backscattered electron detector (BSD) reconstructs a sample image based on the Z of the chemical elements. Therefore, the particles with a higher atomic number appear brighter, and mineral-based shapes could be easily distinguished from the organic ones.

      SEM images of the surface of a membrane used to sample the surface of the drawing L3 (Nudi per la battaglia di Anghiari). The image was obtained in variable pressure, with a backscattered detector on non-metallized samples. The arrow with the asterisk indicates the appearance of the membrane polymer itself. The other arrows indicate fungal spores (Aspergillus or Penicillium). The EDS spectra corresponding to the labeled particles are reported in the plot below. In this drawing several mineral particles were containing lead, silicon or calcium as the main element.

      Discussion Art Meets Third Generation Sequencing Technologies

      In this study we investigated for the first time some of Leonardo da Vinci’s most emblematic drawings using third generation sequencing technology. One of the most limiting steps found when investigating these invaluable objects, namely the sampling, was carried out using a filtering suction system with small nitrocellulose membranes, which allowed a precise and representative sampling of the total surface of the drawings to obtain all the microbiome on their surfaces.

      Furthermore, the tiny nitrocellulose membranes offered the advantage of being easily pooled and used directly for DNA extraction. The total DNA extracted from each drawing was subjected to a metagenomic approach to amplify the entire genomes using the MinION Nanopore sequencing technology. This molecular strategy was already tested by our working group on cultural heritage objects and was shown to be easy for the operator to handle, fast and of a relatively low cost, in comparison with other NGS technologies (Piñar et al., 2020a). The choice of this molecular strategy, including a whole genome amplification protocol (WGA), allowed amplifying all phylogenetic groups at once, enabling to infer the real metagenomes and the relative proportions among the different phylogenetic taxa: eukaryotes, prokaryotes and viruses (Piñar et al., 2020a). An additional advantage of using the Nanopore technology was the length of the sequences generated, which were 4.5 Kb in average, meaning 10-fold longer than with conventional NGS methodologies, as Ion Torrent or Illumina, which yield an average of 400 b length (Kulski, 2016).

      Bacteria vs. Fungi in the Microbiomes of the Drawings

      A somewhat surprising result was the higher proportion of bacteria relative to fungi in all the drawings, as fungi were thought to be always dominant in the microbial communities that colonized cultural heritage objects made or supported on paper (Pinzari et al., 2006). Thus, our results could be explained by (1) the sampling strategy used, which was addressed at sampling primarily the dust deposited on the surfaces, and (2) the molecular approach employed in this study.

      Concerning the first point, the sampling of visibly biodeteriorated paper often lead to the identification of fungi actively growing and using either cellulose or other compounds. In the case of the Leonardo’s drawings no visible biodeterioration, apart from the foxing stains, could be detected. The high proportion of bacteria in relation to fungi could be an index of the type of deposition and the quality of the storage environment. Bacteria in dust are particularly associated with finer dust particles (Yang J. et al., 2020), which could remain in suspension for longer and even penetrate into the containers where the drawings are stored, and which would then act as a selective filter for certain types of dust. Nevertheless, we cannot rule out that the abundance of bacteria in the drawings investigated in this study may also be due to the source of contamination to which they have been subjected in the recent past, namely human contamination due to restoration works, and therefore with the human skin microbiome, as well as contamination with insect droppings and their gut bacteria.

      Regarding the second point, the technologies employed so far could have biased the results in favor of fungi, either through culture-dependent techniques, due to their rapid germination in the appropriate culture medium, or even through molecular techniques. The molecular techniques used so far in paper-based objects, whether of first- (Michaelsen et al., 2006, 2010) or next-generation sequencing (Kraková et al., 2017), have been mainly based on “target amplification” sequencing, which prevents inferring the real proportions of the different target groups (bacteria and fungi). Taking these possible biases into account, this has led to a low number of studies focusing on bacteria colonizing paper materials (Jurado et al., 2010; Michaelsen et al., 2010; Principi et al., 2011; Kraková et al., 2012, 2017). In contrast, there is an extensive literature focused on the colonization of paper by fungi (Sterflinger and Pinzari, 2012). Many fungal species have been described as spoilers of paper-based materials and many of them have shown cellulolytic activities, as some examples species of the genera Chaetomium, Stachybotrys, Cladosporium, Aspergillus, Penicillium. In this context, one of the main research topics in paper and archival material has been related to the pigments produced by paper colonizing fungi, causing foxing, a deterioration of paper that occurs when brownish and reddish spots are produced by the fungi (Arai, 2000). However, so far, no extensive studies have been conducted to elucidate the total microbiome (bacteria vs. fungi) of these valuable objects. The new technology used in this study has been able to give a new focus to the actual proportions of bacteria/fungi that inhabit paper objects.

      Microbiome Biodiversity in the Drawings

      The microbiomes derived from the drawings showed a high biodiversity, with the exception of drawing L4, whose microbiome was heavily contaminated with human sequences, this contamination being perhaps due to its greater handling. In fact, the L4 drawing had massive traces of modern adhesive that required a prolonged conservation work, which was carried out with the aid of a magnification apparatus and various manipulations to which the other drawings were not subjected. This may have resulted in greater removal of surface dust, but at the same time in greater human contamination. Another element to consider in this drawing is the very different nature of its support. The drawing was analyzed with Infrared False Color Photography and Fiber Optics Reflectance Spectroscopy and results showed that the paper was treated with layers of lead white [basic lead carbonate 2PbCO3 × Pb(OH)2] mixed with indigo (Ruggiero, 2020). Moreover, it appeared that Leonardo followed the recipe described by Cennini and Powell Herringham (1899), which consisted in a preparatory layer made with powdered calcinated chicken bones, white lead, indigo from the plant Indigofera tinctoria mixed with animal gelatin and distributed with a brush on the paper support. The design was finally made with a silver tip. This peculiarity may have interfered both with sampling and with the preservation of the microbiome and the environmental DNA in general. Only recent DNA may have remained, i.e., that of the restorers who manipulated it. In this context, Tarsitani et al. (2014) already noted that the period of restoration was probably the critical moment when the discoloration and deterioration of other emblematic drawings by Leonardo da Vinci began, namely in the Atlantic Codex.

      Bacteria

      The results derived from bacteria showed a dominance of Gammaproteobacteria of the order Pseudomonadales. The genus Moraxella was shown to be most abundant at L2, but was also well represented in the L3–L5 drawings, with almost all sequences identified as M. osloensis, a saprophytic species in human skin and mucous membranes (Alkhatib et al., 2017). Its high abundance in these drawings can be explained by contamination due to human handling. Nevertheless, this species was also identified in a survey investigating bulls of indulgence from the fifteenth to sixteenth centuries (Jurado et al., 2010). The genus Acinetobacter was present in drawings L2, L3, and L7, being most of the sequences related to A. baumannii. This species is frequently isolated from hospital environments, although it has been found in environmental samples (Antunes et al., 2014). Nevertheless, sequences related to Acinetobacter spp. were previously identified in a cloning and sequencing survey performed on a deteriorated Italian manuscript dated from the thirteenth century (Michaelsen et al., 2010). The genus Pseudomonas was relative abundant in drawings L7 and L8, with the dominant species P. parafulva, and the genus Psychrobacter was representative solely in L5 with most of the sequences affiliating with Psychrobacter sp. P11G5, which is a psychrophilic species exhibiting laccase activity (Moghadam et al., 2016). Bacterial laccases (Chauhan et al., 2017) are produced either extracellularly or intracellularly and their activity is in a wide range of temperature and pH. They are different from fungal laccases that are multicopper enzymes belonging to the group of blue oxidases (Janusz et al., 2020).

      Members of the Enterobacteriales were detected in drawings L2–L6 but interestingly absent in drawings L7–L8. The genus Salmonella was detected in L2 and L5, being Salmonella enterica subsp. enterica serovars identified in both drawings (Vila Nova et al., 2019). In contrast, Escherichia was identified in L4–L6, with the species E. coli, which is a common inhabitant of the gut microbiota. The genus Proteus, with the species P. mirabilis, which frequently causes urinary tract infections in humans (Pearson et al., 2008), was present exclusively in L3. Interestingly, two genera identified only in this last drawing, namely Buchnera and Xenorhabdus, showed the presence of species that are obligate endosymbionts of aphids, such as B. aphidicola (Shigenobu and Wilson, 2011) and nematodes, such as X. bovienii (Murfin et al., 2015). Interestingly, X. bovienii is a bacterium which has been isolated from several nematode Steinernema species. Steinernema affects many insects and also flies (Archana et al., 2017). Therefore, it is plausible to relate the presence of these species with the presence of insect’s droppings in this drawing.

      The genus Methylovulum was detected only in drawing L3 and all identified sequences affiliated with M. psychrotolerans, a cold-adapted methanotroph from terrestrial environments (Oshkin et al., 2016). Interestingly, the genus Alcanivorax was present in drawings L7 and L8, with the dominant species A. xenomutans, which has capabilities for alkane degradation and heavy-metal resistance (Fu et al., 2018). Finally, the genus Stenotrophomonas, with the species S. rhizophila and S. maltophilia, was only present in L8 with a high relative abundance. It is important to remark that species of this genus have been identified as colonizers and potential spoilers of oil paintings on canvas support (López-Miras et al., 2013), tempera paintings (Zhgun et al., 2020) as well as ancient paper manuscripts (Michaelsen et al., 2010) and book collections (Kraková et al., 2017; Okpalanozie et al., 2018). In addition, Stenotrophomonas sp. has been identified in relation to the biodeterioration of other cultural heritage materials, such as parchments (Piñar et al., 2015a) and historical stone objects (Alakomi et al., 2006; Ettenauer et al., 2011; Rosado et al., 2020). The reason for the perseverance of Stenotrophomonas in this type of materials may be due to their extracellular enzyme activities such as proteases and chitinases (Kobayashi et al., 2002) and their capability to produce extracellular polysaccharide, enabling their adhesion to different materials and providing biofilm formation on their surfaces (Alakomi et al., 2006).

      The Betaproteobacteria were present in all drawings with the dominant genus Ralstonia, with species affiliating to R. solanacearum, R. pickettii, R. insidiosa, and R. mannitolilytica. Ralstonia species, especially R. pickettii and R. solanacearum are prevalent colonizers of printing paper machines (Väisänen et al., 1998) and have been detected as dominant species inhabiting damaged oil paints on canvas (Capodicasa et al., 2010; Piñar et al., 2020a). Furthermore, R. solanacearum is a phytopathogenic bacterium and may exert an active deterioration on paper materials due to the secretion of extracellular proteins, such as cellulolytic enzymes needed for virulence in plants (Liu et al., 2005). Other Betaproteobacteria showed to be relevant only in drawings L7 and L8, as the genera Burkholderia and Cupriavidus. The species B. pseudomallei and B. multivorans were dominant in L7 and Burkholderia sp. MSMB0852 and B. mesoacidophila in L8. Species of this genus were the most frequently isolated from paper products and printing paper machines (Väisänen et al., 1998). Cupriavidus sp. are generally heavy metal tolerant bacteria with the ability to degrade a variety of aromatic hydrocarbon compounds (Wang et al., 2015). The species Cupriavidus sp. USMAA2-4, C. metallidurans and C. gilardii were identified in these two drawings.

      The Alphaproteobacteria class showed to be well represented in drawing L5, with a relative high abundance of species of the genus Caulobacter, with the dominance of C. segnis. Caulobacter species exert a role in decomposition, including the degradation of cellulose, lignin and polyaromatic hydrocarbons and have the capacity to grow on cellulose (Song et al., 2013). This fact may explain the presence of Caulobacter species in this paper-supported drawing and highlights the risk that this bacterium poses to their integrity. This class showed to be relevant also in drawings L7 and L8. The genus Sphingopyxis dominated in both drawings. Sphingopyxis isolates have been reported to degrade aromatic compounds in various habitats (Yang F. et al., 2020). In addition, the genus Paracoccus was only relevant in drawing L7 with the dominance of P. yeei, which has been isolated from diverse natural environments, including human skin (Lim et al., 2018). In contrast, the genus Pannonibacter was only relevant in L8, being all reads affiliated with P. phragmitetus, which is a halotolerant polycyclic aromatic hydrocarbon (PAH)-degrading bacterium (Wang et al., 2016), but it has been also described as a human pathogen (Wang et al., 2017).

      The Actinobacteria were the second most represented phylum within the bacteria. The genus Cryobacterium, with C. articum, showed relatively high abundance in drawings L3–L6. C. articum is a psychrotolerant bacterium isolated from extreme cold environments (Bajerski et al., 2011). This species was previously identified in two oil paintings on canvas (Piñar et al., 2020a). The genus Mycobacterium, with the dominant species M. haemophilum, was present in all the above-mentioned drawings in addition of L7. M. haemophilum is an emerging pathogen associated most commonly with skin infections (Tufariello et al., 2015). The genus Intrasporangium was present in L2 and L7 with the dominant species I. calvum, firstly isolated from the indoor environment (Kalakoutskii et al., 1967). The genus Nocardia was identified in drawings L7 and L8 with the dominant species N. farcinica. This species causes human infections acquired through the respiratory tract or skin (Torres et al., 2000). The genus Streptomyces was also represented in drawings L7 and L8 and in lower proportion in L5. Streptomycetes are ubiquitous in nature and many species are potent producers of lignocellulolytic enzymes including cellulases, hemicellulases and ligninolytic enzymes (Saini et al., 2015). The most dominant species in L7 and L8 was Streptomyces sp. CNQ-509, while Streptomyces albus dominated in drawing L5. The last species is well-known for its capability to produce xylanases (Saini et al., 2015), which are enzymes responsible for the hydrolysis of β-1,4 bonds in plant xylan, the main component of hemicellulose. The detection of species of this genus in these drawings show a latent risk of deterioration due to the potential cellulolytic activity of these strains. Finally, the genus Cutibacterium (Propionibacterium) was identified in L3 and L5, with the species C. acnes in both drawings. This species is considered as human microbial marker, however, its detection should be considered with caution, as it is a potential pathogen (Fitz-Gibbon et al., 2013).

      The phylum Firmicutes was mainly represented by members of the order Bacillales, with the dominance of members of the genus Staphylococcus, present in L3, L5, L6, and L8 (Figure 6). The dominant species in all mentioned drawings was S. aureus, but other species were also identified in high proportions, such as S. equorum, S. capitis, S. pasteuri, and S. epidermidis, among others. The genus Bacillus was present exclusively in L3, with the dominant species B. mycoides and B. thuringiensis. Species belonging to Bacillus and Staphylococcus have been previously identified as the main bacterial species colonizing paper archival materials (Michaelsen et al., 2010; Kraková et al., 2012, 2017; Karakasidou et al., 2018; Okpalanozie et al., 2018).

      Members of the order Lactobacillales showed to be widely represented in almost all drawings, with the genera Streptococcus (in drawings L3–L6), being the emerging pathogen S. pseudopneumoniae dominant in L3–L5. This species is associated with lower-respiratory-tract infections and a high prevalence of antimicrobial resistance (Garriss et al., 2019). In contrast, S. thermophilus showed to be the dominant species in L6, which is a fermentative lactic acid bacterium (Cui et al., 2016). This species shows acidifying capacity, proteolytic activity, fast growth and production of exopolysaccharide (EPS), with all these capabilities posing a risk to the integrity of the paper support of the drawings. The genus Lactobacillus was identified in drawings L5 and L6, with L. salivarius and L. gallinarum as dominant species in L5 and L6, respectively, which are species frequently isolated from the digestive tract of mammals (Neville and O’Toole, 2010) and chicken. The genus Lactococcus was identified in L4, L6, and L7 with L. piscium as dominant species in L4 and L. lactis in L6 and L7. L. lactis has been isolated from a variety of environmental sources but is most commonly associated with fresh or fermented plant material or with milk products (Kelly et al., 2010). L. piscium is a psychrotrophic lactic acid bacterium and is known to be one of the predominant species within spoilage microbial communities in cold-stored packaged foods (Andreevskaya et al., 2015). Finally, members of the Leuconostoc genus were solely identified in L6. Three species showed to be dominant in this drawing, namely L. mesenteroides subsp. dextranicum, L. kimchii, and L. citreum. All of them are commonly used for souring vegetables, producing fermented foods (Dols et al., 1997). However, it is important to note that Leuconostoc species have been detected in association with the fruit fly (Drosophila melanogaster), whose diet may contain communities of fermentation microbes, including a high percentage of Leuconostoc in the fly microbiome (Chandler et al., 2011). Also, the housefly (Musca domestica L.) lives in close association with Leuconostoc (de Jonge et al., 2020), which was detected mainly in adult flies and associated to the buccal parts of the insects. It is conceivable that the discovery of this bacterium is associated with the presence on the drawings of droppings of insects, mainly related to synanthropic diptera (Jones et al., 1999) and dating back to the time when the drawings were not protected and preserved in an appropriate manner (see section “SEM-EDS”).

      Regarding the members of other phyla relatively less abundant in the microbiomes of the drawings, such as Bacteroidetes, Planctomycetes or Cyanobacteria (Figure 4A), it should be mentioned that most of their genera represented less than 0.5% of the total reads. Thus, only the genus Isosphaera, with the species I. pallida (Planctomycetes phylum) accounted for more than 0.5% of the total reads in L7. Finally, the phylum Cyanobacteria was only significantly represented in L3 with species of the genera Anabaena (A. cylindrica) and Oscillatoria (O. nigro-viridis), both accounting for more than 0.5% of the total reads in this drawing. Airborne dispersal of microalgae occur at very low concentration in the atmosphere and their occurrence is little known due to technical limitations in investigating microalgae from air samples (Tesson et al., 2016). However, recent studies showed that microalgae can actually be air transported and form a consistent component of airborne biological particulate, depending on the environmental conditions and the local deposition rates (Genitsaris et al., 2011). Furthermore, both the O. nigro-viridis and the A. cylindrica are cyanobacteria whose tissue can be mineralized with calcium carbonate (Benzerara et al., 2014). Characteristic stromatolite and travertine formations are produced by species precipitating calcium carbonate in their sheaths and within colonial mucilage (Komárek et al., 2003). Some fillers in paper manufacture are based on carbonates from natural sources, such as fossil sedimentary rocks (i.e., biomicrite) (Gerteiser and Laufmann, 1989), which can sometimes be recognized in the form of microfossils among the paper fibers (Bicchieri et al., 2019). It is therefore possible to speculate that the presence of cyanobacteria is due to the use of carbonate-based filler material and natural rocks, perhaps from recent sedimentation, in the manufacture -if not of the Leonardo drawing- perhaps of the cartoons supporting the drawing itself. Both in the case of airborne cyanobacterial propagules and, even more, of biological traces in the sizing minerals, the data confirm the sensitivity of the molecular method used.

      It is important to note that it is difficult to compare the data generated in this study with those published in the literature, due to the different technologies applied. As mentioned above, there are not many studies focused on the detection of bacteria on paper-materials and even fewer conducting culture-independent approaches. However, the data generated in this study through third-generation sequencing analyses are somewhat consistent with previous studies where first- and second-generation sequencing was performed. The core microbiome found in this study includes mainly bacteria belonging to the phyla Proteobacteria, Actinobacteria, and Firmicutes, which correlates well with the few previously existing molecular studies. These three phyla were previously detected by different cloning approaches performed on photographic material, including cellulose supports (Bučková et al., 2014; Puškárová et al., 2016) as well as on an ancient Italian manuscript (Michaelsen et al., 2010). Principi et al. (2011) also showed the presence of members of the phyla Proteobacteria and Firmicutes in the pages of Leonardo da Vinci’s Atlantic Codex. Later, Kraková et al., 2017 confirmed that the core microbiome of various archival materials, investigated with next generation sequencing technology (Illumina MiSeq), consisted mainly of members of the three phyla mentioned above.

      To summarize our data, the bacterial species detected in this study can be roughly grouped according to their origin on the paper-support of the drawings. A first group is formed by bacteria that are mainly contaminants introduced by human manipulation, mainly species of the genera Staphylococcus, Streptococcus, and Cutibacterium (Propionibacterium), but also species of the genus Moraxella, Paracoccus, Mycobacterium, and Nocardia could have been introduced by a mishandling of the drawings. This result is reinforced by a large amount of human DNA found in all the drawings. A second group is formed by bacteria that may have been deposited in the drawings through vectors, such as insects. These bacteria belong to the genera Salmonella, Escherichia, Buchnera, Xenorhabdus, and Leuconostoc. This result is supported by the finding of insect droppings on the drawings. The third group consists of bacteria that are actually potential spoilers of the paper-support due to their capabilities of producing extracellular enzymes, and they carry a latent risk to the integrity of the drawings. They are members of the genera Ralstonia, Stenotrophomonas, Burkholderia, Caulobacter, Streptomyces, and Bacillus. Some of them could have been introduced into the paper during the production process and others have been deposited on its surface later through the air. However, it is important to note that the bacterial community was influenced by the geographical location of the drawings, Rome or Turin, denoting an environmental influence.

      Fungi

      Among the detected fungi, only some members of the Ascomycota phylum represented more than 0.5% of total reads in the microbiomes of the drawings, among them, members of Sordariomycetes, Eurotiomycetes, Dothideomycetes, and Saccharomycetes. Within the Sordariomycetes, the genus Colletotrichum showed to be present in L2, L7–L8. This genus includes several plant pathogens (Cannon et al., 2012). The species C. higginsianum was dominant in L7–L8 while C. orchidophilum was dominant in L2. Interestingly, the genus Grosmannia was only present at L5 in a very high proportion, with G. clavigera being the sole species identified. This species is the main fungal pathogen associated with the mountain pine beetle and causes blue stain in wood (Wang et al., 2014). The fungus is linked to American conifers. However, it is important to note that individuals of Pseudotsuga menziesii are used as a perennial ornamental conifer species in the northern regions of Italy. It is presumed that the presence of this fungus could be due to a contamination event related to its presence in Rome at the conservation institute where the drawings were analyzed and sampled. In the historical garden of the institute there are in fact ornamental evergreen conifer species. The fungus can also live as a saprotroph in its anamorphic state on bark and cellulosic debris (Leptographium form). A further hypothesis is that the drawing came into contact with furniture or woody material, possibly raw, contaminated by the fungus. The DNA in G. clavigera mycelium remains stable for more than 10 years also in heat-treated wood samples (Wong et al., 2020). Finally, the genus Trichoderma was present in L7 and L8, with the species T. reesei as the dominant species in both drawings. Trichoderma reesei is one of the main producers of cellulolytic and xylanolytic enzymes, being the most studied fungus involved with lignocellulosic degradation (Bischof et al., 2016). The detection of this fungus supposes a great latent risk of deterioration for these paper-supported drawings if the environmental conditions become suitable for its germination.

      Members of the Eurotiomycetes were present in drawings L3, L5–L8, being Aspergillus widespread in all these drawings, while Penicillium showed to dominate in L3 and Penicilliopsis in L7. The species A. glaucus, A. clavatus, and A. niger were dominant in all those drawings. Aspergillus species are well known for the production of cellulases and their main role in the biodeterioration of paper (Das et al., 1997). The most dominant species found in the drawings, namely A. glaucus, is known to have a wide environmental distribution due to its physiological hardiness under extreme conditions and to its ability to grow on substrates with a low water activity (Hubka et al., 2013). A. glaucus has been detected in cultural assets made of wooden (Sterflinger et al., 2018b) and canvas support (Piñar et al., 2020a), but its isolation and undistinguishable morphological identification were difficult in the past (Hubka et al., 2013). Penicillium rubens showed to be the dominant species in L3. This species is a common fungus of the indoor environment and well known for its resilience under low water activities (van Laarhoven et al., 2017). Finally, Penicilliopsis zonata was the dominant members of Eurotiomycetes in L7. It is important to note that in the indoor environments, the relative ratio and abundances of fungal species belonging to Aspergillus, Eurotium and Penicillium genera are used along with indicator species to evaluate the presence of hidden moisture damage in buildings (ISO/DIS 16000-17:2006, 2006; Baudisch et al., 2009). Most indoor contamination studies are based on the composition of the dust in suspension. The fungal species in the dust deposited on surfaces are less studied and there are less significant statistics. However, the fact that the drawings, even those preserved in the same environment, presented a different mycoflora is an interesting fact as it allows to formulate hypotheses on the conservation conditions to which the objects have been exposed. The predominance of xerophilic fungal species such as Penicillium rubens and Aspergillus glaucus, indicates that the drawings have not been preserved in humid environments.

      Within the Dothideomycetes only members of two genera accounted for more than 0.5% of total reads, namely Alternaria and Bipolaris. Alternaria, with the only detected species A. alternata was present in L3 and L5. This fungus is a well-known paper-spoiler (Das et al., 1997; Michaelsen et al., 2006; Mesquita et al., 2009) and also considered as an allergen, imposing a severe health threat on librarians, conservators, and scholars who might be in contact with infected paper-supported materials (Karakasidou et al., 2018). It must be mentioned, however, that the Alternaria is a fungus that needs water available in the materials and that it typically develops after flooding and wetting events in building materials, including wood and paper. Therefore, the presence of hydrophilic fungal species on the drawings, in the absence of any visible manifestation of water damage, is due to their common presence in the air as airborne contaminants, also coming from the external environment, from the decomposing plant material typically present in parks and gardens. Bipolaris sorokiniana was the main species identified in L7 and L8. B. sorokiniana (anamorph of Cochliobolus sativus) is a phytopathogenic fungus that produce cell wall-degrading enzymes, such as cellulases, which aid invasion into host cells (Aich et al., 2017). As far as we know, B. sorokiniana has not been considered as a paper degrader but, in fact, the fungus is provided with cellulases that can degrade paper fibers and thus, this species poses a latent risk of deterioration for the drawings in which it has been detected.

      Members of the Saccharomycetes, with Ascoidea and Candida species were detected in L2–L3 and L5–L6 but curiously, not in drawings L7–L8. A. rubescens, C. orthopsilosis, and C. dubliniensis were the dominant species detected in all four mentioned drawings. Ascoidea species appear primarily to be disseminated by insect vectors (Batra, 1963) and have not been identified in paper materials so far. In contrast, Candida spp. have been previously identified in archival documents and photographic material, but no evidence regarding implication in paper deterioration has been provided (Principi et al., 2011; Borrego et al., 2015; Karakasidou et al., 2018). The detection and prevalence of species of the genus Candida in repositories and archives may be an indication of their ubiquity in this environment, where they can be introduced from the outside environment or contaminated materials after production. Moreover, certain species of Candida are typical human microbiota, which might suggest anthropogenic pollution by handling the collections (Borrego et al., 2015, 2017).

      In summary, most of the fungi detected in this study have been previously identified in other studies focused on paper-supported materials and include species of the genera Trichoderma, Aspergillus, Penicillium, Penicilliopsis, Alternaria, and Candida. On the contrary, as far as we know, species of other genera, such as Colletotrichum and Bipolaris, which are mainly plant pathogens, have not been described as common in paper materials. However, they possess extracellular enzymes that may be directly responsible for paper deterioration. Finally, the genera Grosmannia and Ascoidea include typical species associated with insects and could have been introduced to the surface of the drawings through these vectors. As observed with the bacterial communities, it was also possible to see a correlation of the fungal communities with their geographical location, highlighting an environmental and geographical influence on the microbiomes of the drawings.

      Monitoring of Drawing L2 “<italic>Autoritratto</italic>”

      A previous study of the L2 drawing performed in 2012 showed that both the recto and verso had been affected by a fungal attack by the species Eurotium halophilicum, a fungus capable of developing under conditions of extremely low substrate water content (<65%). In fact, the fungus is considered as xerophilic and needs hypertonified culture media by adding 15% NaCL for its normal growth, which can also be replaced by other osmoactive organic substances such as sugars (Imai, 1964). In association with the fungal structures, abundant calcium oxalate crystals, resulting from the metabolic activity of the fungus itself, were observed (Pinzari et al., 2014; Piñar et al., 2015b). In 2012 the self-portrait did not show the presence of active infections by microorganisms but only traces of previous attacks and general dust contamination. The membranes and swabs used to sample the dust present on the self-portrait were used to extract the DNA of the microorganisms and to carry out both an evaluation of the type of contamination of the self-portrait (microbiological diversity and abundance) and a comparison, using molecular techniques, of the effectiveness of surface sampling. Details of the study are given in Piñar et al. (2015b), where the authors showed that small membranes are particularly effective for non-invasive sampling from paper and for direct use in DNA extraction protocols. Several fungal species were found to be significantly associated with both the recto and the verso of the drawing at that time, such as the Phialosimplex, Penicillium, and Acremonium species. From the molecular analysis, however, it was not possible to obtain gene sequences correlated to the fungus E. halophilicum, which was instead diagnosed by SEM analysis, thanks to some unmistakable morphological characters. A possible explanation lies in the treatments that the portrait has undergone over the years. In 1987 Dr. Fausta Gallo, director of the biology laboratory of the former “Istituto Centrale per la Patologia del Libro (Rome)” examined the tracks produced between 1972 and 1987 by the thermo-hydrographers placed in the drawer where the self-portrait was kept. In the spring-summer periods, the portrait was subjected to a Relative Humidity of 70–80% with peaks of 90%. The alarm was such that Dr. Fausta Gallo wrote the same year (1987) that the damage on the portrait was evolving and that disinfection with ethylene oxide was necessary (Pinzari et al., 2014). Ethylene oxide is a disinfectant gas, widely used in the past to sterilize book material, which works as an alkylating agent of nucleic acids and proteins. The chemical action of this compound makes the DNA and RNA of microorganisms (and materials such as parchment) no longer amplifiable (Michaelsen et al., 2013). There is no certainty in the documentation available at ICPAL that the self-portrait has actually undergone disinfection with ethylene oxide. However, if so, any trace of DNA present on the object and prior to treatment would be undetectable with current analysis techniques. In this case, the infection by the fungus E. halophilicum, evident in the images obtained with SEM but not detected either by molecular techniques based on cloning in 2012, and currently by nanopore technology, would be dated back to before 1987. The design has not been treated in any way since 2012. The new series of molecular analyses, carried out by means of Nanopore sequencing, a much more advanced technique than cloning, should have detected the presence of the fungus, which instead did not emerge. Other typical air contaminating species were detected, despite their much lower density than that documented for the E. halophilicum. Therefore, the absence of the fungal species among those found with Nanopore sequencing supports the hypothesis that the design was actually treated with ethylene oxide in 1987.

      SEM-EDX

      The incrustations analyzed so far were found to be insect droppings (Figure 8). This evidence is due partly to the presence of spherical particles of biogenic nature composed of magnesium and calcium phosphate recognized as “spherulites” which are biogenic stones that form in the gut of insects or in their excretion products (Gower, 2008). These spherical micro-objects are used in archeology as markers of the presence or use of animals’ dung (Lancelotti and Madella, 2012). The other encrustations, those without spherulites, appeared to be instead rich in nitrogen compounds, and potassium sulfate, that are both catabolites of insects’ metabolism (Cochran, 1975). The small monoclinic and prismatic crystals characterizing these objects were scarcely distinguishable in the matrix when visualized with BSD detector, which means that they were not presenting a high atomic number, and could presumably be uric salts (uric acid: C5H4N4O3). Nitrogen could make cellulose potentially more degradable since nitrogen-based compounds are substances normally absent in paper media, but which microorganisms need to grow (Pinzari, 2011). The insects’ droppings found on the L4 drawing were colonized by fungal mycelium. Leonardo’s drawings have been kept for many years locked in vaults, protected by cardboard boxes and other protective systems. To imagine that an insect can rest on a drawing to leave its droppings is unrealistic. Therefore, the waxy incrustations removed and analyzed could date back for many years. That is, when the drawings were exposed to the air, so even before they were hung in glass frames, as some appear to have been around the 1930s and 1950s (Misiti, 2014). The species of insects that cause damage to books and libraries belong to different families, united by having representatives who have adapted to living in man-made environments and have spread throughout the world thanks to commercial activities. According to their biological characteristics, they carry out their life cycle within the materials they feed on or in the surrounding environment (Florian, 1997). In the case of the traces found on the drawings, rather than the classic insects harmful to paper and parchment, which typically cause lacunas and tracks such as cockroaches, woodworms and silverfishes, it seems that the droppings are more characteristic of species of diptera such as domestic flies and related species. House flies (Musca domestica) are synanthropic and commonly found on decaying matter but also feces and food. They have a complex prokaryotic microbiota in their gut, that they can release on materials through droppings and regurgitation of food. Furthermore, the microbiome in domestic flies can be relatively consistent across geographic locations and habitat (Park et al., 2019).

      Finally, we would like to make a general discussion about the conservation over the centuries of these unique and very fragile objects, in terms of the need to remove insect excrements, dust and the remains of fungi and bacteria, since they are all substances rich in nitrogen, and chemically reactive compounds, capable of oxidizing and degrading the paper fibers. Microorganisms that are still viable also pose other risks, as they can proliferate and cause much more serious damage under certain environmental conditions. However, the manipulation of the drawings could also lead to possible additional contamination, as demonstrated in this study by the presence of bacteria from the human microbiome and the abundant human DNA. However, thanks to the “dirt” found on these heritage objects it has been possible to study their history and to draw a distinctive profile that forms part of the drawings themselves and describes them. The information obtained in this study comes from samples that would be eliminated in a routine conservative treatment. Incrustations, amorphous particles, pieces of glue and remains of supporting cardboard are materials that have been removed by restorers up to now. However, in light of what can be discovered from their meticulous study, the process of conserving certain artistic objects probably needs to be reconsidered because part of their history is written in their “dirt”.

      Conclusion

      This study shows for the first time the complete microbiomes of some of Leonardo da Vinci’s most emblematic drawings. The strategy applied for their research, namely the Nanopore sequencing technology, is one of the most advanced so far to investigate all the phylogenetic groups included in each of the microbiomes. In addition, the complementation with SEM analysis allowed the visualization of insect droppings in some of the drawings, which are directly correlated with many of the microorganisms found on their surfaces.

      The results show a surprising dominance of bacteria over fungi. This fact is partly due to the sampling strategy used, aimed at sampling mainly the dust deposited on the surface of the drawings, but also to the source of contamination to which they have been subjected in the recent past, resulting in a high proportion of bacteria typical of the human microbiome, introduced by intensive handling of the drawings during restoration works, in addition to other microorganisms that have obviously been introduced as contaminants through vectors, such as insects and their excrements.

      It is important to clarify that, unlike other case studies that point to visible paper damage caused by fungi, often also associated with bacteria, no clear damage apart from foxing was observed in this study. The sensitivity of the Nanopore sequencing method, which provides information on the abundance of the different taxa in each sample, could have an advantage over other methods for monitoring an ongoing deterioration process, because it allows the assessment of quantitative relationships between taxa and also the deduction of variations in these relationships due to detrimental situations.

      In general, the biological information contained in the investigated drawings showed very specific microbiomes, which can be used as a bio-archive of the history the objects, providing a kind of fingerprint for current and future comparisons. However, some similarities were observed that could be influenced by the geographical location of the drawings. Some taxa were only detected in the drawings stored in the Royal Library of Turin, while other taxa were only detected in the drawings stored in the Corsinian Library in Rome.

      The results show a high proportion of human DNA, especially in drawing L4, which indicates an inadequate manipulation of this drawing, most probably during the previous restoration processes. This result indicates that restoration techniques should be continually reviewed based on scientific evidence and that monitoring studies in the field of conservation and restoration should be encouraged.

      Finally, in order to maximize the value of the samples that can be obtained from unique artworks, it is critical to optimize the analytical processes through protocols that allow the conservation and reuse of the samples in various analyses. Analytical techniques based on the study of environmental DNA will enable the maximum exploitation of natural bio-archives represented by dust, biofilm or other traces of past events present in the works. In perspective, it is possible to imagine a pipeline of measures aimed at obtaining the maximum information from samples taken from unique objects, such as Leonardo’s drawings, which allow very little sampling. Measuring the vitality and activity of individual microbial cells present in materials is probably the next challenge. There are techniques that have already been consolidated in other fields of biology, such as cytofluorimetric measures and single-cell RNA sequencing, which examines the sequence information of individual cells with optimized NGS technologies, which would allow a better understanding of the function of an individual microorganism in the context of its microenvironment. In the light of the new research techniques and their resolving power, some considerations on the impact that cleaning and consolidation techniques may have on the information content of ancient materials are now really urgent. Likewise, new methods and protocols are needed, both for the documentation and for the recovery of these natural “bio-archives” that accompany all the unique objects of our artistic and cultural heritage, without weighing down the already arduous work of conservators and restorers.

      Data Availability Statement

      The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/, Bioproject PRJNA655381.

      Author Contributions

      GP conceived the design of the study, conducted molecular analyses, processed sequencing data, and wrote the manuscript. MSc and PC conducted sampling and participated in SEM analyses. FP performed and analyzed SEM and SEM-EDS and contributed to the writing of the manuscript. AG analyzed the sequencing data and contributed to the writing of the manuscript. MSe and KS conceived and contributed to the design of the study and supervised the project. All authors approved the final version for submission.

      Conflict of Interest

      The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

      Funding. This study was funded by the “Istituto Centrale per la Patologia degli Archivi e del Libro,” Ministero per i Beni e le Attività Culturali e per il Turismo, Rome, Italy.

      We thank the ICPAL staff and the restorers who helped to sample the drawings.

      Supplementary Material

      The Supplementary Material for this article can be found online at: /articles/10.3389/fmicb.2020.593401/full#supplementary-material

      References Aich S. Singh R. K. Kundu P. Pandey S. P. Datta S. (2017). Genome-wide characterization of cellulases from the hemi-biotrophic plant pathogen, Bipolaris sorokiniana, reveals the presence of a highly stable GH7 endoglucanase. Biotechnol. Biofuels 10:135. 10.1186/s13068-017-0822-0 28559926 Alakomi H. L. Paananen A. Suihko M. L. Helander I. M. Saarela M. (2006). Weakening effect of cell permeabilizers on Gram-negative bacteria causing biodeterioration. Appl. Environ. Microbiol. 72 46954703. 10.1128/AEM.00142-06 16820461 Alkhatib N. J. Younis M. H. Alobaidi A. S. Shaath N. M. (2017). An unusual osteomyelitis caused by Moraxella osloensis: a case report. Int. J. Surg. Case Rep. 41 146149. 10.1016/j.ijscr.2017.10.022 29078156 Alwadi A. Kilby J. Gawanmeh A. (2017). Tracking and automating a library system using radio frequency identification technology. Int. J. Smart Sens. Intell. Syst. 10 425450. 10.21307/ijssis-2017-219 Andreevskaya M. Johansson P. Laine P. Smolander O. P. Sonck M. Rahkila R. (2015). Genome sequence and transcriptome analysis of meat-spoilage-associated lactic acid bacterium Lactococcus piscium MKFS47. Appl. Environ. Microbiol. 81 38003811. 10.1128/AEM.00320-15 25819958 Antunes L. C. S. Visca P. Towner K. J. (2014). Acinetobacter baumannii: evolution of a global pathogen. Pathog. Dis. 71 292301. 10.1111/2049-632X.12125 24376225 Arai H. (2000). Foxing caused by fungi: twenty-five years of study. Int. Biodeter. Biodegr. 46 181188. 10.1016/S0964-8305(00)00063-9 Archana M. D’Souza P. E. Patil J. (2017). Efficacy of entomopathogenic nematodes (Rhabditida: Steinernematidae and Heterorhabditidae) on developmental stages of house fly. Musca domestica. J. Parasit. Dis. 41 782794. 10.1007/s12639-017-0889-3 28848279 Bajerski F. Ganzert L. Mangelsdorf K. Lipski A. Wagner D. (2011). Cryobacterium arcticum sp. nov., a psychrotolerant bacterium from an Arctic soil. Int. J. Syst. Evol. Microbiol. 61 18491853. 10.1099/ijs.0.027128-0 20817836 Banerjee A. Pons T. Lequeux N. Dubertret B. (2016). Quantum dots-DNA bioconjugates: synthesis to applications. Interface Focus 6:20160064. 10.1098/rsfs.2016.0064 27920898 Batra L. R. (1963). Life history of species of Ascoidea. Amer. J. Bot. 50:624. Baudisch C. Assadian O. Kramer A. (2009). Concentration of the genera Aspergillus, Eurotium and Penicillium in 63-μm house dust fraction as a method to predict hidden moisture damage in homes. BMC Public Health 9:247. 10.1186/1471-2458-9-247 19615082 Benzerara K. Skouri-Panet F. Li J. Férard C. Gugger M. Laurent T. (2014). Intracellular Ca-carbonate biomineralization is widespread in cyanobacteria. Proc. Natl. Acad. Sci. U.S.A. 111 1093310938. 10.1073/pnas.1403510111 25009182 Bicchieri M. Biocca P. Colaizzi P. Pinzari F. (2019). Microscopic observations of paper and parchment: the archaeology of small objects. Herit. Sci. 7:47. 10.1186/s40494-019-0291-9 Bischof R. H. Ramoni J. Seiboth B. (2016). Cellulases and beyond: the first 70 years of the enzyme producer Trichoderma reesei. Microb. Cell. Fact. 15 113. 10.1186/s12934-016-0507-6 27287427 Borrego S. Molina A. Santana A. (2015). Mold on stored photographs and maps: a case study. Top. Photogr. Preserv. 16 109120. Borrego S. Molina A. Santana A. (2017). Fungi in archive repositories environments and the deterioration of the graphics documents. EC Microbiol. 11 205226. Bučková M. Puškárová A. Sclocchi M. C. Bicchieri M. Colaizzi P. Pinzari F. (2014). Co-occurrence of bacteria and fungi and spatial partitioning during photographic materials biodeterioration. Polym. Degrad. Stab. 108 111. 10.1016/j.polymdegradstab.2014.05.025 Cannon P. F. Damm U. Johnston P. R. Weir B. S. (2012). Colletotrichum - current status and future directions. Stud. Mycol. 73 181213. 10.3114/sim0014 23136460 Capodicasa S. Fedi S. Porcelli A. M. Zannoni D. (2010). The microbial community dwelling on a biodeteriorated 16th century painting. Int. Biodeter. Biodegr. 64 727733. 10.1016/j.ibiod.2010.08.006 Cennini C. Powell Herringham C. J. (1899). The Book of the Art of Cennino Cennini: A Contemporary Practical Treatise On Quattrocento Painting. London: G. Allen & Unwin, ltd. Chandler J. A. Lang J. M. Bhatnagar S. Eisen J. A. Kopp A. (2011). Bacterial communities of diverse Drosophila species: ecological context of a host-microbe model system. PLoS Genet. 7:e1002272. 10.1371/journal.pgen.1002272 21966276 Chauhan P. S. Goradia B. Saxena A. (2017). Bacterial laccase: recent update on production, properties and industrial applications. 3 Biotech 7:323. 10.1007/s13205-017-0955-7 28955620 Cochran D. G. (1975). “Excretion in Insects,” in Insect Biochemistry and Function, eds Candy D. J. Kilby B. A. (Boston, MA: Springer). Cui Y. Xu T. Qu X. Hu T. Jiang X. Zhao C. (2016). New insights into various production characteristics of Streptococcus thermophilus strains. Int. J. Mol. Sci. 17:1701. 10.3390/ijms17101701 27754312 Das M. K. L. Prasad J. S. Ahmad S. K. (1997). Endoglucanase production by paper-degrading mycoflora. Letters Appl. Microbiol. 25 313315. 10.1046/j.1472-765X.1997.00217.x 9418063 de Carvalho H. P. Sequeira S. O. Pinho D. Trovão J. Fernandes da Costa R. M. Egas C. (2019). Combining an innovative non-invasive sampling method and high-throughput sequencing to characterize fungal communities on a canvas painting. Int. Biodeter. Biodegr. 145:104816. 10.1016/j.ibiod.2019.104816 de Jonge N. Michaelsen T. Y. Ejbye-Ernst R. Jensen A. Nielsen M. E. Bahrndorff S. (2020). Housefly (Musca domestica L.) associated microbiota across different life stages. Sci. Rep. 10:7842. 10.1038/s41598-020-64704-y Demidov V. V. (2020). “Hiding and storing messages and data in DNA,” in DNA Beyond Genes, ed. Demidov V. V.. Chap. 2, (Switzerland AG: Springer, Cham. Springer Nature), 723. 10.1007/978-3-030-36434-2 Dols M. Chraibi W. Remaud-Simeon M. Lindley N. D. Monsan P. F. (1997). Growth and energetics of Leuconostoc mesenteroides NRRL B-1299 during metabolism of various sugars and their consequences for dextransucrase production. Appl. Environ. Microbiol. 63 21592165. 10.1128/AEM.63.6.2159-2165.1997 9172334 Ettenauer J. Piñar G. Sterflinger K. Gonzalez-Muñoz M. T. Jroundi F. (2011). Molecular monitoring of the microbial dynamics occurring on historical limestone buildings during and after the in situ application of different bio-consolidation treatments. Sci. Total Environ. 409 53375352. 10.1016/j.scitotenv.2011.08.063 21944202 Fitz-Gibbon S. Tomida S. Chiu B. H. Nguyen L. Du C. Liu M. (2013). Propionibacterium acnes strain populations in the human skin microbiome associated with acne. J. Invest. Dermatol. 133 21522160. 10.1038/jid.2013.21 23337890 Florian M. L. (1997). Heritage Eaters. London: James & James, 164. Fu X. Lai Q. Dong C. Wang W. Shao Z. (2018). Complete genome sequence of Alcanivorax xenomutans P40, an alkane-degrading bacterium isolated from deep seawater. Mar. Genomics. 38 14. 10.1016/j.margen.2017.05.010 Garriss G. Nannapaneni P. Simões A. S. Browall S. Subramanian K. Sá-Leão R. (2019). Genomic characterization of the emerging pathogen Streptococcus pseudopneumoniae. mBio 10:e001286-19. 10.1128/mBio.01286-19 31239383 Genitsaris S. Kormas K. A. Moustaka-Gouni M. (2011). Airborne algae and cyanobacteria: occurrence and related health effects. Front. Biosci. 3:772787. 10.2741/e285 21196350 Gerteiser N. Laufmann M. (1989). Effects of natural ground calcium carbonate on wood free papermaking process and paper quality. Appita 42 295300. Goodwin S. McPherson J. McCombie W. (2016). Coming of age: ten years of next-generation sequencing technologies. Nat. Rev. Genet. 17 333351. 10.1038/nrg.2016.49 27184599 Gower L. B. (2008). Biomimetic model systems for investigating the amorphous precursor pathway and its role in biomineralization. Chem. Rev. 108 45514627. 10.1021/cr800443h 19006398 Grottoli A. Beccaccioli M. Zoppi E. Fratini R. S. Schifano E. Santarelli M. L. (2020). Nanopore sequencing and bioinformatics for rapidly identifying cultural heritage spoilage microorganisms. Front. Mater. 7:14. 10.3389/fmats.2020.00014 Hubka V. Kolarík M. Kubátová A. Peterson S. W. (2013). Taxonomic revision of Eurotium and transfer of species to Aspergillus. Mycologia 105 912937. 10.3852/12-151 Imai M. (1964). On a Tonophilic Fungus Isolated from an Old Agar Slant. Shokubutsugaku Zasshi 77 4348. ISO/DIS 16000-17:2006 (2006). Indoor air – Part 17: Detection and enumeration of moulds – Culture-Based Method., ICS: 13.040.20. Berlin: Beuth-Verlag. Janusz G. Pawlik A. Świderska-Burek U. Polak J. Sulej J. Jarosz-Wilkołazka A. (2020). Laccase properties, physiological functions, and evolution. Int. J. Mol. Sci. 21:966. 10.3390/ijms21030966 32024019 Jones C. Isard S. Cortinas M. (1999). Dispersal of synanthropic Diptera: lessons from the past and technology for the future. Ann. Entomol. Soc. Am. 92 829839. 10.1093/aesa/92.6.829 Jurado V. Porca E. Pastrana M. P. Cuezva S. Fernandez-Cortes A. Saiz-Jimenez C. (2010). Microbiological study of bulls of indulgence of the 15th–16th centuries. Sci. Total Environ. 408 37113715. 10.1016/j.scitotenv.2010.05.008 20595060 Kalakoutskii L. V. Kirillova I. P. Krassilnikov N. A. (1967). A new genus of the Actinomycetales - Intrasporangium gen. nov. J. Gen. Microbiol. 48 7985. Karakasidou K. Nikolouli K. Amoutzias G. D. Pournou A. Manassis C. Tsiamis G. (2018). Microbial diversity in biodeteriorated Greek historical documents dating back to the 19th and 20th century: a case study. MicrobiologyOpen. 7:e00596. 10.1002/mbo3.596 29484839 Kelly W. J. Ward L. J. H. Leahy S. C. (2010). Chromosomal diversity in Lactococcus lactis and the origin of dairy starter cultures. ıGenome Biol. Evol. 2 729744. 10.1093/gbe/evq056 20847124 Kisová Z. Planý M. Pavlović J. Bučková M. Puškárová A. Kraková L. (2020). Biodeteriogens characterization and molecular analyses of diverse funeral accessories from XVII Century. Appl. Sci. 10:5451. 10.3390/app10165451 Kobayashi D. Y. Reedy R. M. Bick J. A. Oudemans P. V. (2002). Characterization of a chitinase gene from Stenotrophomonas maltophilia strain 34S1 and its involvement in biological control. Appl. Environ. Microbiol. 68 10471054. 10.1128/AEM.68.3.1047-1054.2002 11872449 Komárek J. Kling H. Komárková J. (2003). “Filamentous cyanobacteria,” in Aquatic Ecology, Freshwater Algae of North America, eds Wehr J. D. Sheath R. G. (Cambridge, MA: Academic Press), 117196. Kraková L. Chovanová K. Selim S. A. Šimonovičová A. Puškarová A. Maková A. (2012). A multiphasic approach for investigation of the microbial diversity and its biodegradative abilities in historical paper and parchment documents. Int. Biodeter. Biodegr. 70 117125. 10.1016/j.ibiod.2012.01.011 Kraková L. Šoltys K. Otlewska A. Pietrzak K. Purkrtová S. Savická D. (2017). Comparison of methods for identification of microbial communities in book collections: culture-dependent (sequencing and MALDI-TOF MS) and culture-independent (Illumina MiSeq). Int. Biodeter. Biodegr. 131 5159. 10.1016/j.ibiod.2017.02.015 Kulski J. K. (2016). “Next-generation sequencing - an overview of the history, tools, and “Omic” applications,” in Next Generation Sequencing - Advances, Applications and Challenges, ed. Kulski J. (Norderstedt: Books on Demand), 10.5772/61964 Lancelotti C. Madella M. (2012). The ‘invisible’ product: developing markers for identifying dung in archaeological contexts. J. Archaeol. Sci. 39 953963. 10.1016/j.jas.2011.11.007 Lim J. Y. Hwang I. Ganzorig M. Pokhriyal S. Singh R. Lee K. (2018). Complete genome sequence of Paracoccus yeei TT13, isolated from human skin. Genome Announc. 6 e1514e1517. 10.1128/genomeA.01514-17 29348361 Liu H. Zhang S. Schell M. A. Denny T. P. (2005). Pyramiding unmarked deletions in Ralstonia solanacearum shows that secreted proteins in addition to plant cell-wall-degrading enzymes contribute to virulence. Mol. Plant. Microbe Interact. 18 12961305. 10.1094/MPMI-18-1296 16478049 López-Miras M. M. Martín-Sánchez I. Yebra-Rodríguez A. Romero-Noguera J. Bolívar-Galiano F. Ettenauer J. (2013). Contribution of the microbial communities detected on an oil painting on canvas to its biodeterioration. PLoS One 8:e80198. 10.1371/journal.pone.0080198 24312203 Marvasi M. Cavalieri D. Mastromei G. Casaccia A. Perito B. (2019). Omics technologies for an in-depth investigation of biodeterioration of cultural heritage. Int. Biodeter. Biodegr. 144:104736. 10.1016/j.ibiod.2019.104736 McHugh M. L. (2011). Multiple comparison analysis testing in ANOVA. Biochem. Med. 21 203209. 10.11613/bm.2011.029 22420233 McMurdie P. J. Holmes S. (2013). Phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One 8:e61217. 10.1371/journal.pone.0061217 23630581 Mesquita N. Portugal A. Videira S. Rodrıguez-Echeverri A. A. Bandeira A. M. L. Santos M. J. A. (2009). Fungal diversity in ancient documents: a case study on the Archive of the University of Coimbra. Int. Biodeter. Biodegr. 63 626629. 10.1016/j.ibiod.2009.03.010 Michaelsen A. Piñar G. Pinzari F. (2010). Molecular and microscopical investigation of the microflora inhabiting a deteriorated Italian manuscript dated from the thirteenth century. Microb. Ecol. 60 6980. 10.1007/s00248-010-9667-9 20449583 Michaelsen A. Pinzari F. Barbabietola N. Piñar G. (2013). Monitoring the effects of different conservation treatments on paper-infecting fungi. Int. Biodeter. Biodegr. 84 333341. 10.1016/j.ibiod.2012.08.005 24092956 Michaelsen A. Pinzari F. Ripka K. Lubitz W. Piñar G. (2006). Application of molecular techniques for the identification of fungal communities colonizing paper material. Int. Biodeter. Biodegr. 58 133141. 10.1016/j.ibiod.2006.06.019 Misiti M. C. (2014). Atti del Seminario Internazionale I disegni di Leonardo. Diagnostica Conservazione Tutela, ed. Misiti M. C. (Roma: ICRCPAL-Sillabe). Moghadam M. S. Albersmeier A. Winkler A. Cimmino L. Rise K. Hohmann-Marriott M. F. (2016). Isolation and genome sequencing of four Arctic marine Psychrobacter strains exhibiting multicopper oxidase activity. BMC Genomics 17:117. 10.1186/s12864-016-2445-4 26879123 Murfin K. E. Lee M. M. Klassen J. L. McDonald B. R. Larget B. Forst S. (2015). Xenorhabdus bovienii strain diversity impacts coevolution and symbiotic maintenance with Steinernema spp. nematode hosts. mBio 6:e00076-15. 10.1128/mBio.00076-15 26045536 Neville B. A. O’Toole P. W. (2010). Probiotic properties of Lactobacillus salivarius and closely related Lactobacillus species. Future Microbiol. 5 759774. 10.2217/fmb.10.35 20441548 Okpalanozie O. E. Adebusoye S. A. Troiano F. Cattò C. Ilori M. O. Cappitelli F. (2018). Assessment of indoor air environment of a Nigerian museum library and its biodeteriorated books using culture-dependent and –independent techniques. Int. Biodeter. Biodegr. 132 139149. 10.1016/j.ibiod.2018.03.003 Oshkin I. Y. Belova S. E. Danilova O. V. Miroshnikov K. K. Rijpstra W. I. C. Sinninghe Damsté J. S. (2016). Methylovulum psychrotolerans sp. nov., a cold-adapted methanotroph from low-temperature terrestrial environments, and emended description of the genus Methylovulum. Int. J. Syst. Evol. Microbiol. 66 24172423. 10.1099/ijsem.0.001046 27031985 Park R. Dzialo M. C. Spaepen S. Nsabimana D. Gielens K. Devriese H. (2019). Microbial communities of the house fly Musca domestica vary with geographical location and habitat. Microbiome 7:147. 10.1186/s40168-019-0748-9 31699144 Pearson M. M. Sebaihia M. Churcher C. Quail M. A. Seshasayee A. S. Luscombe N. M. (2008). Complete genome sequence of uropathogenic Proteus mirabilis, a master of both adherence and motility. J. Bacteriol. 190 40274037. 10.1128/JB.01981-07 18375554 Piñar G. Poyntner C. Lopandic K. Tafer H. Sterflinger K. (2020a). Rapid diagnosis of biological colonization in cultural artefacts using the MinION Nanopore sequencing technology. Int. Biodeter. Biodegr. 148:104908. 10.1016/j.ibiod.2020.104908 Piñar G. Tafer H. Schreiner M. Miklas H. Sterflinger K. (2020b). Decoding the biological information contained in two ancient Slavonic parchment codices: an added historical value. Environ. Microbiol. 22 32183233. 10.1111/1462-2920.15064 Piñar G. Poyntner C. Tafer H. Sterflinger K. (2019). A time travel story: metagenomic analyses decipher the unknown geographical shift and the storage history of possibly smuggled antique marble statues. Ann. Microbiol. 69 10011021. 10.1007/s13213-019-1446-3 Piñar G. Sterflinger K. Pinzari F. (2015a). Unmasking the measles-like parchment discoloration: molecular and micro-analytical approach. Environ. Microbiol. 17 427443. 10.1111/1462-2920.12471 24684276 Piñar G. Tafer H. Sterflinger K. Pinzari F. (2015b). Amid the possible causes of a very famous foxing: molecular and microscopic insight into Leonardo da Vinci’s selfportrait. Environ. Microbiol. Rep. 7 849859. 10.1111/1758-2229.12313 26111623 Pinzari F. (2011). “Microbial ecology of indoor environments. The ecological and applied aspects of microbial contamination in archives, libraries and conservation environments,” in Sick Building Syndrome in Public Buildings and Workplaces, Chap. 9, ed Abdul-Wahab S.-A. (Berlin: Springer). Pinzari F. Pasquariello G. De Mico A. (2006). Biodeterioration of paper: a SEM study of fungal spoilage reproduced under controlled conditions. Macromol. Symp. 238 5766. 10.1002/masy.200650609 Pinzari F. Sclocchi M. C. Colaizzi P. Valenti P. (2014). “Indagini microbiologiche e microscopiche,” in Atti del Seminario Internazionale I disegni di Leonardo. Diagnostica Conservazione Tutela, eds Misiti M. C. Misiti M. C. (Roma: ICRCPAL-Sillabe), 5665. Principi P. Villa F. Sorlini C. Cappitelli F. (2011). Molecular studies of microbial community structure on stained pages of Leonardo da Vinci’s Atlantic Codex. Microb. Ecol. 61 214222. 10.1007/s00248-010-9741-3 20811884 Puškárová A. Bučková M. Habalová B. Kraková L. Maková A. Pangallo D. (2016). Microbial communities affecting albumen photography heritage: a methodological survey. Sci Rep. 6:20810. 10.1038/srep20810 26864429 Reimann C. Filzmoser P. Garret R. G. Dutter R. (2008). Statistical Data Analysis Explained: Applied Environmental Statistics with R. Correlation. Hoboken, NJ: John Wiley & Sons Ltd, 181191. 10.1002/9780470987605.ch11 Rosado T. Dias L. Lança M. Nogueira C. Santos R. Martins M. R. (2020). Assessment of microbiota present on a Portuguese historical stone convent using high-throughput sequencing approaches. MicrobiologyOpen 9 10671084. 10.1002/mbo3.1030 32352657 Ruggiero D. (2020). “Spettroscopia e riflettografia sulle mediazioni grafiche nei disegni di Leonardo: la tecnologia ci racconta il passato con la voce del futuro,” in Leonardo e i suoi segreti, Publisher: Gangemi, 5373 (in press). Saini A. Aggarwal N. K. Sharma A. Yadav A. (2015). Actinomycetes: a source of lignocellulolytic enzymes. Enzyme Res. 2015:279381. 10.1155/2015/279381 26793393 Schabereiter-Gurtner C. Piñar G. Lubitz W. Rölleke S. (2001). An advanced strategy to identify bacterial communities on art objects. J. Microbiol. Methods 45 7787. 10.1016/s0167-7012(01)00227-5 Schadt E. E. Turner S. Kasarskis A. (2010). A window into third-generation sequencing. Hum. Mol. Genet. 19 227240. 10.1093/hmg/ddq416 20858600 Shigenobu S. Wilson A. C. (2011). Genomic revelations of a mutualism: the pea aphid and its obligate bacterial symbiont. Cell. Mol. Life Sci. 68 12971309. 10.1007/s00018-011-0645-2 21390549 Šoltys K. Planý M. Biocca P. Vianello V. Bučková M. Puškárová A. (2019). Lead soaps formation and biodiversity in a XVIII Century wax seal coloured with minium. Environ. Microbiol. 22 15171534. 10.1111/1462-2920.14735 31291047 Song N. Cai H. Y. Yan Z. S. Jiang H. L. (2013). Cellulose degradation by one mesophilic strain Caulobacter sp. FMC1 under both aerobic and anaerobic conditions. Bioresour. Technol. 131 281287. 10.1016/j.biortech.2013.01.003 23357088 Steinberg M. D. Žura I. Murkovic Steinberg I. (2014). Wireless smart tag with on-board conductometric chemical sensor. Sens. Actuators B Chem. 196 208214. 10.1016/j.snb.2014.02.012 Sterflinger K. Little B. Piñar G. Pinzari F. de los Rios A. Gu J. D. (2018a). Future directions and challenges in biodeterioration research on historic materials and cultural properties. Int. Biodeter. Biodegr. 129 1012. 10.1016/j.ibiod.2017.12.007 Sterflinger K. Voitl C. Lopandic K. Piñar G. Tafer H. (2018b). Big sound and extreme fungi—xerophilic, halotolerant Aspergilli and Penicillia with low optimal temperature as invaders of historic pipe organs. Life 8:22. 10.3390/life8020022 29903995 Sterflinger K. Piñar G. (2013). Microbial deterioration of cultural heritage and works of art-tilting at windmills? Appl. Microbiol. Biotechnol. 97 96379646. 10.1007/s00253-013-5283-1 24100684 Sterflinger K. Pinzari F. (2012). The revenge of time: fungal deterioration of cultural heritage with particular reference to books, paper and parchment. Environ. Microbiol. 14 559566. 10.1111/j.1462-2920.2011.02584.x 22004478 Tarsitani G. Moroni C. Cappitelli F. Pasquariello G. Maggi O. (2014). Microbiological analysis of surfaces of Leonardo da Vinci’s Atlantic Codex: Biodeterioration risk. Int. J. Microbiol. 2014:214364. 10.1155/2014/214364 25574171 Teasdale M. D. Fiddyment S. Vnouček J. Mattiangeli V. Speller C. Binois A. (2017). The York Gospels: a 1000-year biological palimpsest. R. Soc. Open Sci. 4:170988. 10.1098/rsos.170988 29134095 Tesson S. Skjøth C. A. Šantl-Temkiv T. Löndahl J. (2016). Airborne microalgae: insights, opportunities, and challenges. Appl. Environ. Microbiol. 82 19781991. 10.1128/AEM.03333-15 26801574 Torres O. Domingo P. Pericas R. Boiron P. Montiel J. A. Vázquez G. (2000). Infection caused by Nocardia farcinica: case report and review. Eur. J. Clin. Microbiol. Infect. Dis. EJCMID 19 205212. 10.1007/s100960050460 10795594 Tufariello J. M. Kerantzas C. A. Vilchčze C. Calder R. B. Nordberg E. K. Fischer J. A. (2015). The complete genome sequence of the emerging pathogen Mycobacterium haemophilum explains its unique culture requirements. mBio 6 e01313e1315. 10.1128/mBio.01313-15 26578674 Väisänen O. M. Weber A. Bennasar A. Rainey F. A. Busse H. J. Salkinoja-Salonen M. S. (1998). Microbial communities of printing paper machines. J. Appl. Microbiol. 84 10691084. 10.1046/j.1365-2672.1998.00447.x 9717292 Valentini F. Calcaterra A. Antonaroli S. Talamo M. (2018). Smart portable devices suitable for cultural heritage: a review. Sensors 18:2434. 10.3390/s18082434 30050013 van Laarhoven K. A. Peeters L. Bekker M. Huinink H. P. Adan O. (2017). Conidia of Penicillium rubens formed at low water activities can attract more water. MicrobiologyOpen 6:e00526. 10.1002/mbo3.526 28872258 Vila Nova M. Durimel K. La K. Felten A. Bessičres P. Mistouet M. Y. (2019). Genetic and metabolic signatures of Salmonella enterica subsp. enterica associated with animal sources at the pangenomic scale. BMC Genomics 20:814. 10.1186/s12864-019-6188-x 31694533 Wang M. Zhang X. Jiang T. Hu S. Yi Z. Zhou Y. (2017). Liver abscess caused by Pannonibacter phragmitetus: case report and literature review. Front. Med. 4:48. 10.3389/fmed.2017.00048 28487855 Wang X. Chen M. Xiao J. Hao L. Crowley D. E. Zhang Z. (2015). Genome sequence analysis of the naphthenic acid degrading and metal resistant bacterium Cupriavidus gilardii CR3. PLoS One 10:e0132881. 10.1371/journal.pone.0132881 26301592 Wang X. Jin D. Zhou L. Zhang Z. (2016). Draft genome sequence of Pannonibacter phragmitetus strain CGMCC9175, a halotolerant polycyclic aromatic hydrocarbon-degrading bacterium. Genome Announc. 4 e1675e1615. 10.1128/genomeA.01675-15 26823598 Wang Y. Lim L. Madilao L. Lah L. Bohlmann J. Breuil C. (2014). Gene discovery for enzymes involved in limonene modification or utilization by the mountain pine beetle-associated pathogen Grosmannia clavigera. Appl. Environ. Microbiol. 80 45664576. 10.1128/AEM.00670-14 24837377 Wong B. Leal I. Feau N. Dale A. Uzunovic A. Hamelin R. C. (2020). Molecular assays to detect the presence and viability of Phytophthora ramorum and Grosmannia clavigera. PLoS One 15:e0221742. 10.1371/journal.pone.0221742 Yang J. Kim E. K. Park H. J. McDowell A. Kim Y. K. (2020). The impact of bacteria-derived ultrafine dust particles on pulmonary diseases. Exp. Mol. Med. 52 338347. 10.1038/s12276-019-0367-3 32203101 Yang F. Feng H. Massey I. Y. Huang F. Guo J. Zhang X. (2020). Genome-wide analysis reveals genetic potential for aromatic compounds biodegradation of Sphingopyxis. Biomed Res. Int. 2020:5849123. 10.1155/2020/5849123 32596333 Zhgun A. Avdanina D. Shumikhin K. Simonenko N. Lyubavskaya E. Volkov I. (2020). Detection of potential biodeterioration risks for tempera painting in 16th century exhibits from State Tretyakov Gallery. PLoS One 15:e0230591. 10.1371/journal.pone.0230591 32240187
      ‘Oh, my dear Thomas, you haven’t heard the terrible news then?’ she said. ‘I thought you would be sure to have seen it placarded somewhere. Alice went straight to her room, and I haven’t seen her since, though I repeatedly knocked at the door, which she has locked on the inside, and I’m sure it’s most unnatural of her not to let her own mother comfort her. It all happened in a moment: I have always said those great motor-cars shouldn’t be allowed to career about the streets, especially when they are all paved with cobbles as they are at Easton Haven, which are{331} so slippery when it’s wet. He slipped, and it went over him in a moment.’ My thanks were few and awkward, for there still hung to the missive a basting thread, and it was as warm as a nestling bird. I bent low--everybody was emotional in those days--kissed the fragrant thing, thrust it into my bosom, and blushed worse than Camille. "What, the Corner House victim? Is that really a fact?" "My dear child, I don't look upon it in that light at all. The child gave our picturesque friend a certain distinction--'My husband is dead, and this is my only child,' and all that sort of thing. It pays in society." leave them on the steps of a foundling asylum in order to insure [See larger version] Interoffice guff says you're planning definite moves on your own, J. O., and against some opposition. Is the Colonel so poor or so grasping—or what? Albert could not speak, for he felt as if his brains and teeth were rattling about inside his head. The rest of[Pg 188] the family hunched together by the door, the boys gaping idiotically, the girls in tears. "Now you're married." The host was called in, and unlocked a drawer in which they were deposited. The galleyman, with visible reluctance, arrayed himself in the garments, and he was observed to shudder more than once during the investiture of the dead man's apparel. HoME香京julia种子在线播放 ENTER NUMBET 0016www.huhu1.com.cn
      www.epepiy.com.cn
      www.jgycos.com.cn
      www.nxchain.com.cn
      qukk.com.cn
      ngdctl.com.cn
      nyriff.com.cn
      www.tdpgkz.com.cn
      wnchain.com.cn
      xfcztd.com.cn
      处女被大鸡巴操 强奸乱伦小说图片 俄罗斯美女爱爱图 调教强奸学生 亚洲女的穴 夜来香图片大全 美女性强奸电影 手机版色中阁 男性人体艺术素描图 16p成人 欧美性爱360 电影区 亚洲电影 欧美电影 经典三级 偷拍自拍 动漫电影 乱伦电影 变态另类 全部电 类似狠狠鲁的网站 黑吊操白逼图片 韩国黄片种子下载 操逼逼逼逼逼 人妻 小说 p 偷拍10幼女自慰 极品淫水很多 黄色做i爱 日本女人人体电影快播看 大福国小 我爱肏屄美女 mmcrwcom 欧美多人性交图片 肥臀乱伦老头舔阴帝 d09a4343000019c5 西欧人体艺术b xxoo激情短片 未成年人的 插泰国人夭图片 第770弾み1 24p 日本美女性 交动态 eee色播 yantasythunder 操无毛少女屄 亚洲图片你懂的女人 鸡巴插姨娘 特级黄 色大片播 左耳影音先锋 冢本友希全集 日本人体艺术绿色 我爱被舔逼 内射 幼 美阴图 喷水妹子高潮迭起 和后妈 操逼 美女吞鸡巴 鸭个自慰 中国女裸名单 操逼肥臀出水换妻 色站裸体义术 中国行上的漏毛美女叫什么 亚洲妹性交图 欧美美女人裸体人艺照 成人色妹妹直播 WWW_JXCT_COM r日本女人性淫乱 大胆人艺体艺图片 女同接吻av 碰碰哥免费自拍打炮 艳舞写真duppid1 88电影街拍视频 日本自拍做爱qvod 实拍美女性爱组图 少女高清av 浙江真实乱伦迅雷 台湾luanlunxiaoshuo 洛克王国宠物排行榜 皇瑟电影yy频道大全 红孩儿连连看 阴毛摄影 大胆美女写真人体艺术摄影 和风骚三个媳妇在家做爱 性爱办公室高清 18p2p木耳 大波撸影音 大鸡巴插嫩穴小说 一剧不超两个黑人 阿姨诱惑我快播 幼香阁千叶县小学生 少女妇女被狗强奸 曰人体妹妹 十二岁性感幼女 超级乱伦qvod 97爱蜜桃ccc336 日本淫妇阴液 av海量资源999 凤凰影视成仁 辰溪四中艳照门照片 先锋模特裸体展示影片 成人片免费看 自拍百度云 肥白老妇女 女爱人体图片 妈妈一女穴 星野美夏 日本少女dachidu 妹子私处人体图片 yinmindahuitang 舔无毛逼影片快播 田莹疑的裸体照片 三级电影影音先锋02222 妻子被外国老头操 观月雏乃泥鳅 韩国成人偷拍自拍图片 强奸5一9岁幼女小说 汤姆影院av图片 妹妹人艺体图 美女大驱 和女友做爱图片自拍p 绫川まどか在线先锋 那么嫩的逼很少见了 小女孩做爱 处女好逼连连看图图 性感美女在家做爱 近距离抽插骚逼逼 黑屌肏金毛屄 日韩av美少女 看喝尿尿小姐日逼色色色网图片 欧美肛交新视频 美女吃逼逼 av30线上免费 伊人在线三级经典 新视觉影院t6090影院 最新淫色电影网址 天龙影院远古手机版 搞老太影院 插进美女的大屁股里 私人影院加盟费用 www258dd 求一部电影里面有一个二猛哥 深肛交 日本萌妹子人体艺术写真图片 插入屄眼 美女的木奶 中文字幕黄色网址影视先锋 九号女神裸 和骚人妻偷情 和潘晓婷做爱 国模大尺度蜜桃 欧美大逼50p 西西人体成人 李宗瑞继母做爱原图物处理 nianhuawang 男鸡巴的视屏 � 97免费色伦电影 好色网成人 大姨子先锋 淫荡巨乳美女教师妈妈 性nuexiaoshuo WWW36YYYCOM 长春继续给力进屋就操小女儿套干破内射对白淫荡 农夫激情社区 日韩无码bt 欧美美女手掰嫩穴图片 日本援交偷拍自拍 入侵者日本在线播放 亚洲白虎偷拍自拍 常州高见泽日屄 寂寞少妇自卫视频 人体露逼图片 多毛外国老太 变态乱轮手机在线 淫荡妈妈和儿子操逼 伦理片大奶少女 看片神器最新登入地址sqvheqi345com账号群 麻美学姐无头 圣诞老人射小妞和强奸小妞动话片 亚洲AV女老师 先锋影音欧美成人资源 33344iucoom zV天堂电影网 宾馆美女打炮视频 色五月丁香五月magnet 嫂子淫乱小说 张歆艺的老公 吃奶男人视频在线播放 欧美色图男女乱伦 avtt2014ccvom 性插色欲香影院 青青草撸死你青青草 99热久久第一时间 激情套图卡通动漫 幼女裸聊做爱口交 日本女人被强奸乱伦 草榴社区快播 2kkk正在播放兽骑 啊不要人家小穴都湿了 www猎奇影视 A片www245vvcomwwwchnrwhmhzcn 搜索宜春院av wwwsee78co 逼奶鸡巴插 好吊日AV在线视频19gancom 熟女伦乱图片小说 日本免费av无码片在线开苞 鲁大妈撸到爆 裸聊官网 德国熟女xxx 新不夜城论坛首页手机 女虐男网址 男女做爱视频华为网盘 激情午夜天亚洲色图 内裤哥mangent 吉沢明歩制服丝袜WWWHHH710COM 屌逼在线试看 人体艺体阿娇艳照 推荐一个可以免费看片的网站如果被QQ拦截请复制链接在其它浏览器打开xxxyyy5comintr2a2cb551573a2b2e 欧美360精品粉红鲍鱼 教师调教第一页 聚美屋精品图 中韩淫乱群交 俄罗斯撸撸片 把鸡巴插进小姨子的阴道 干干AV成人网 aolasoohpnbcn www84ytom 高清大量潮喷www27dyycom 宝贝开心成人 freefronvideos人母 嫩穴成人网gggg29com 逼着舅妈给我口交肛交彩漫画 欧美色色aV88wwwgangguanscom 老太太操逼自拍视频 777亚洲手机在线播放 有没有夫妻3p小说 色列漫画淫女 午间色站导航 欧美成人处女色大图 童颜巨乳亚洲综合 桃色性欲草 色眯眯射逼 无码中文字幕塞外青楼这是一个 狂日美女老师人妻 爱碰网官网 亚洲图片雅蠛蝶 快播35怎么搜片 2000XXXX电影 新谷露性家庭影院 深深候dvd播放 幼齿用英语怎么说 不雅伦理无需播放器 国外淫荡图片 国外网站幼幼嫩网址 成年人就去色色视频快播 我鲁日日鲁老老老我爱 caoshaonvbi 人体艺术avav 性感性色导航 韩国黄色哥来嫖网站 成人网站美逼 淫荡熟妇自拍 欧美色惰图片 北京空姐透明照 狼堡免费av视频 www776eom 亚洲无码av欧美天堂网男人天堂 欧美激情爆操 a片kk266co 色尼姑成人极速在线视频 国语家庭系列 蒋雯雯 越南伦理 色CC伦理影院手机版 99jbbcom 大鸡巴舅妈 国产偷拍自拍淫荡对话视频 少妇春梦射精 开心激动网 自拍偷牌成人 色桃隐 撸狗网性交视频 淫荡的三位老师 伦理电影wwwqiuxia6commqiuxia6com 怡春院分站 丝袜超短裙露脸迅雷下载 色制服电影院 97超碰好吊色男人 yy6080理论在线宅男日韩福利大全 大嫂丝袜 500人群交手机在线 5sav 偷拍熟女吧 口述我和妹妹的欲望 50p电脑版 wwwavtttcon 3p3com 伦理无码片在线看 欧美成人电影图片岛国性爱伦理电影 先锋影音AV成人欧美 我爱好色 淫电影网 WWW19MMCOM 玛丽罗斯3d同人动画h在线看 动漫女孩裸体 超级丝袜美腿乱伦 1919gogo欣赏 大色逼淫色 www就是撸 激情文学网好骚 A级黄片免费 xedd5com 国内的b是黑的 快播美国成年人片黄 av高跟丝袜视频 上原保奈美巨乳女教师在线观看 校园春色都市激情fefegancom 偷窥自拍XXOO 搜索看马操美女 人本女优视频 日日吧淫淫 人妻巨乳影院 美国女子性爱学校 大肥屁股重口味 啪啪啪啊啊啊不要 操碰 japanfreevideoshome国产 亚州淫荡老熟女人体 伦奸毛片免费在线看 天天影视se 樱桃做爱视频 亚卅av在线视频 x奸小说下载 亚洲色图图片在线 217av天堂网 东方在线撸撸-百度 幼幼丝袜集 灰姑娘的姐姐 青青草在线视频观看对华 86papa路con 亚洲1AV 综合图片2区亚洲 美国美女大逼电影 010插插av成人网站 www色comwww821kxwcom 播乐子成人网免费视频在线观看 大炮撸在线影院 ,www4KkKcom 野花鲁最近30部 wwwCC213wapwww2233ww2download 三客优最新地址 母亲让儿子爽的无码视频 全国黄色片子 欧美色图美国十次 超碰在线直播 性感妖娆操 亚洲肉感熟女色图 a片A毛片管看视频 8vaa褋芯屑 333kk 川岛和津实视频 在线母子乱伦对白 妹妹肥逼五月 亚洲美女自拍 老婆在我面前小说 韩国空姐堪比情趣内衣 干小姐综合 淫妻色五月 添骚穴 WM62COM 23456影视播放器 成人午夜剧场 尼姑福利网 AV区亚洲AV欧美AV512qucomwwwc5508com 经典欧美骚妇 震动棒露出 日韩丝袜美臀巨乳在线 av无限吧看 就去干少妇 色艺无间正面是哪集 校园春色我和老师做爱 漫画夜色 天海丽白色吊带 黄色淫荡性虐小说 午夜高清播放器 文20岁女性荫道口图片 热国产热无码热有码 2015小明发布看看算你色 百度云播影视 美女肏屄屄乱轮小说 家族舔阴AV影片 邪恶在线av有码 父女之交 关于处女破处的三级片 极品护士91在线 欧美虐待女人视频的网站 享受老太太的丝袜 aaazhibuo 8dfvodcom成人 真实自拍足交 群交男女猛插逼 妓女爱爱动态 lin35com是什么网站 abp159 亚洲色图偷拍自拍乱伦熟女抠逼自慰 朝国三级篇 淫三国幻想 免费的av小电影网站 日本阿v视频免费按摩师 av750c0m 黄色片操一下 巨乳少女车震在线观看 操逼 免费 囗述情感一乱伦岳母和女婿 WWW_FAMITSU_COM 偷拍中国少妇在公车被操视频 花也真衣论理电影 大鸡鸡插p洞 新片欧美十八岁美少 进击的巨人神thunderftp 西方美女15p 深圳哪里易找到老女人玩视频 在线成人有声小说 365rrr 女尿图片 我和淫荡的小姨做爱 � 做爱技术体照 淫妇性爱 大学生私拍b 第四射狠狠射小说 色中色成人av社区 和小姨子乱伦肛交 wwwppp62com 俄罗斯巨乳人体艺术 骚逼阿娇 汤芳人体图片大胆 大胆人体艺术bb私处 性感大胸骚货 哪个网站幼女的片多 日本美女本子把 色 五月天 婷婷 快播 美女 美穴艺术 色百合电影导航 大鸡巴用力 孙悟空操美少女战士 狠狠撸美女手掰穴图片 古代女子与兽类交 沙耶香套图 激情成人网区 暴风影音av播放 动漫女孩怎么插第3个 mmmpp44 黑木麻衣无码ed2k 淫荡学姐少妇 乱伦操少女屄 高中性爱故事 骚妹妹爱爱图网 韩国模特剪长发 大鸡巴把我逼日了 中国张柏芝做爱片中国张柏芝做爱片中国张柏芝做爱片中国张柏芝做爱片中国张柏芝做爱片 大胆女人下体艺术图片 789sss 影音先锋在线国内情侣野外性事自拍普通话对白 群撸图库 闪现君打阿乐 ady 小说 插入表妹嫩穴小说 推荐成人资源 网络播放器 成人台 149大胆人体艺术 大屌图片 骚美女成人av 春暖花开春色性吧 女亭婷五月 我上了同桌的姐姐 恋夜秀场主播自慰视频 yzppp 屄茎 操屄女图 美女鲍鱼大特写 淫乱的日本人妻山口玲子 偷拍射精图 性感美女人体艺木图片 种马小说完本 免费电影院 骑士福利导航导航网站 骚老婆足交 国产性爱一级电影 欧美免费成人花花性都 欧美大肥妞性爱视频 家庭乱伦网站快播 偷拍自拍国产毛片 金发美女也用大吊来开包 缔D杏那 yentiyishu人体艺术ytys WWWUUKKMCOM 女人露奶 � 苍井空露逼 老荡妇高跟丝袜足交 偷偷和女友的朋友做爱迅雷 做爱七十二尺 朱丹人体合成 麻腾由纪妃 帅哥撸播种子图 鸡巴插逼动态图片 羙国十次啦中文 WWW137AVCOM 神斗片欧美版华语 有气质女人人休艺术 由美老师放屁电影 欧美女人肉肏图片 白虎种子快播 国产自拍90后女孩 美女在床上疯狂嫩b 饭岛爱最后之作 幼幼强奸摸奶 色97成人动漫 两性性爱打鸡巴插逼 新视觉影院4080青苹果影院 嗯好爽插死我了 阴口艺术照 李宗瑞电影qvod38 爆操舅母 亚洲色图七七影院 被大鸡巴操菊花 怡红院肿么了 成人极品影院删除 欧美性爱大图色图强奸乱 欧美女子与狗随便性交 苍井空的bt种子无码 熟女乱伦长篇小说 大色虫 兽交幼女影音先锋播放 44aad be0ca93900121f9b 先锋天耗ばさ无码 欧毛毛女三级黄色片图 干女人黑木耳照 日本美女少妇嫩逼人体艺术 sesechangchang 色屄屄网 久久撸app下载 色图色噜 美女鸡巴大奶 好吊日在线视频在线观看 透明丝袜脚偷拍自拍 中山怡红院菜单 wcwwwcom下载 骑嫂子 亚洲大色妣 成人故事365ahnet 丝袜家庭教mp4 幼交肛交 妹妹撸撸大妈 日本毛爽 caoprom超碰在email 关于中国古代偷窥的黄片 第一会所老熟女下载 wwwhuangsecome 狼人干综合新地址HD播放 变态儿子强奸乱伦图 强奸电影名字 2wwwer37com 日本毛片基地一亚洲AVmzddcxcn 暗黑圣经仙桃影院 37tpcocn 持月真由xfplay 好吊日在线视频三级网 我爱背入李丽珍 电影师傅床戏在线观看 96插妹妹sexsex88com 豪放家庭在线播放 桃花宝典极夜著豆瓜网 安卓系统播放神器 美美网丝袜诱惑 人人干全免费视频xulawyercn av无插件一本道 全国色五月 操逼电影小说网 good在线wwwyuyuelvcom www18avmmd 撸波波影视无插件 伊人幼女成人电影 会看射的图片 小明插看看 全裸美女扒开粉嫩b 国人自拍性交网站 萝莉白丝足交本子 七草ちとせ巨乳视频 摇摇晃晃的成人电影 兰桂坊成社人区小说www68kqcom 舔阴论坛 久撸客一撸客色国内外成人激情在线 明星门 欧美大胆嫩肉穴爽大片 www牛逼插 性吧星云 少妇性奴的屁眼 人体艺术大胆mscbaidu1imgcn 最新久久色色成人版 l女同在线 小泽玛利亚高潮图片搜索 女性裸b图 肛交bt种子 最热门有声小说 人间添春色 春色猜谜字 樱井莉亚钢管舞视频 小泽玛利亚直美6p 能用的h网 还能看的h网 bl动漫h网 开心五月激 东京热401 男色女色第四色酒色网 怎么下载黄色小说 黄色小说小栽 和谐图城 乐乐影院 色哥导航 特色导航 依依社区 爱窝窝在线 色狼谷成人 91porn 包要你射电影 色色3A丝袜 丝袜妹妹淫网 爱色导航(荐) 好男人激情影院 坏哥哥 第七色 色久久 人格分裂 急先锋 撸撸射中文网 第一会所综合社区 91影院老师机 东方成人激情 怼莪影院吹潮 老鸭窝伊人无码不卡无码一本道 av女柳晶电影 91天生爱风流作品 深爱激情小说私房婷婷网 擼奶av 567pao 里番3d一家人野外 上原在线电影 水岛津实透明丝袜 1314酒色 网旧网俺也去 0855影院 在线无码私人影院 搜索 国产自拍 神马dy888午夜伦理达达兔 农民工黄晓婷 日韩裸体黑丝御姐 屈臣氏的燕窝面膜怎么样つぼみ晶エリーの早漏チ○ポ强化合宿 老熟女人性视频 影音先锋 三上悠亚ol 妹妹影院福利片 hhhhhhhhsxo 午夜天堂热的国产 强奸剧场 全裸香蕉视频无码 亚欧伦理视频 秋霞为什么给封了 日本在线视频空天使 日韩成人aⅴ在线 日本日屌日屄导航视频 在线福利视频 日本推油无码av magnet 在线免费视频 樱井梨吮东 日本一本道在线无码DVD 日本性感诱惑美女做爱阴道流水视频 日本一级av 汤姆avtom在线视频 台湾佬中文娱乐线20 阿v播播下载 橙色影院 奴隶少女护士cg视频 汤姆在线影院无码 偷拍宾馆 业面紧急生级访问 色和尚有线 厕所偷拍一族 av女l 公交色狼优酷视频 裸体视频AV 人与兽肉肉网 董美香ol 花井美纱链接 magnet 西瓜影音 亚洲 自拍 日韩女优欧美激情偷拍自拍 亚洲成年人免费视频 荷兰免费成人电影 深喉呕吐XXⅩX 操石榴在线视频 天天色成人免费视频 314hu四虎 涩久免费视频在线观看 成人电影迅雷下载 能看见整个奶子的香蕉影院 水菜丽百度影音 gwaz079百度云 噜死你们资源站 主播走光视频合集迅雷下载 thumbzilla jappen 精品Av 古川伊织star598在线 假面女皇vip在线视频播放 国产自拍迷情校园 啪啪啪公寓漫画 日本阿AV 黄色手机电影 欧美在线Av影院 华裔电击女神91在线 亚洲欧美专区 1日本1000部免费视频 开放90后 波多野结衣 东方 影院av 页面升级紧急访问每天正常更新 4438Xchengeren 老炮色 a k福利电影 色欲影视色天天视频 高老庄aV 259LUXU-683 magnet 手机在线电影 国产区 欧美激情人人操网 国产 偷拍 直播 日韩 国内外激情在线视频网给 站长统计一本道人妻 光棍影院被封 紫竹铃取汁 ftp 狂插空姐嫩 xfplay 丈夫面前 穿靴子伪街 XXOO视频在线免费 大香蕉道久在线播放 电棒漏电嗨过头 充气娃能看下毛和洞吗 夫妻牲交 福利云点墦 yukun瑟妃 疯狂交换女友 国产自拍26页 腐女资源 百度云 日本DVD高清无码视频 偷拍,自拍AV伦理电影 A片小视频福利站。 大奶肥婆自拍偷拍图片 交配伊甸园 超碰在线视频自拍偷拍国产 小热巴91大神 rctd 045 类似于A片 超美大奶大学生美女直播被男友操 男友问 你的衣服怎么脱掉的 亚洲女与黑人群交视频一 在线黄涩 木内美保步兵番号 鸡巴插入欧美美女的b舒服 激情在线国产自拍日韩欧美 国语福利小视频在线观看 作爱小视颍 潮喷合集丝袜无码mp4 做爱的无码高清视频 牛牛精品 伊aⅤ在线观看 savk12 哥哥搞在线播放 在线电一本道影 一级谍片 250pp亚洲情艺中心,88 欧美一本道九色在线一 wwwseavbacom色av吧 cos美女在线 欧美17,18ⅹⅹⅹ视频 自拍嫩逼 小电影在线观看网站 筱田优 贼 水电工 5358x视频 日本69式视频有码 b雪福利导航 韩国女主播19tvclub在线 操逼清晰视频 丝袜美女国产视频网址导航 水菜丽颜射房间 台湾妹中文娱乐网 风吟岛视频 口交 伦理 日本熟妇色五十路免费视频 A级片互舔 川村真矢Av在线观看 亚洲日韩av 色和尚国产自拍 sea8 mp4 aV天堂2018手机在线 免费版国产偷拍a在线播放 狠狠 婷婷 丁香 小视频福利在线观看平台 思妍白衣小仙女被邻居强上 萝莉自拍有水 4484新视觉 永久发布页 977成人影视在线观看 小清新影院在线观 小鸟酱后丝后入百度云 旋风魅影四级 香蕉影院小黄片免费看 性爱直播磁力链接 小骚逼第一色影院 性交流的视频 小雪小视频bd 小视频TV禁看视频 迷奸AV在线看 nba直播 任你在干线 汤姆影院在线视频国产 624u在线播放 成人 一级a做爰片就在线看狐狸视频 小香蕉AV视频 www182、com 腿模简小育 学生做爱视频 秘密搜查官 快播 成人福利网午夜 一级黄色夫妻录像片 直接看的gav久久播放器 国产自拍400首页 sm老爹影院 谁知道隔壁老王网址在线 综合网 123西瓜影音 米奇丁香 人人澡人人漠大学生 色久悠 夜色视频你今天寂寞了吗? 菲菲影视城美国 被抄的影院 变态另类 欧美 成人 国产偷拍自拍在线小说 不用下载安装就能看的吃男人鸡巴视频 插屄视频 大贯杏里播放 wwwhhh50 233若菜奈央 伦理片天海翼秘密搜查官 大香蕉在线万色屋视频 那种漫画小说你懂的 祥仔电影合集一区 那里可以看澳门皇冠酒店a片 色自啪 亚洲aV电影天堂 谷露影院ar toupaizaixian sexbj。com 毕业生 zaixian mianfei 朝桐光视频 成人短视频在线直接观看 陈美霖 沈阳音乐学院 导航女 www26yjjcom 1大尺度视频 开平虐女视频 菅野雪松协和影视在线视频 华人play在线视频bbb 鸡吧操屄视频 多啪啪免费视频 悠草影院 金兰策划网 (969) 橘佑金短视频 国内一极刺激自拍片 日本制服番号大全magnet 成人动漫母系 电脑怎么清理内存 黄色福利1000 dy88午夜 偷拍中学生洗澡磁力链接 花椒相机福利美女视频 站长推荐磁力下载 mp4 三洞轮流插视频 玉兔miki热舞视频 夜生活小视频 爆乳人妖小视频 国内网红主播自拍福利迅雷下载 不用app的裸裸体美女操逼视频 变态SM影片在线观看 草溜影院元气吧 - 百度 - 百度 波推全套视频 国产双飞集合ftp 日本在线AV网 笔国毛片 神马影院女主播是我的邻居 影音资源 激情乱伦电影 799pao 亚洲第一色第一影院 av视频大香蕉 老梁故事汇希斯莱杰 水中人体磁力链接 下载 大香蕉黄片免费看 济南谭崔 避开屏蔽的岛a片 草破福利 要看大鸡巴操小骚逼的人的视频 黑丝少妇影音先锋 欧美巨乳熟女磁力链接 美国黄网站色大全 伦蕉在线久播 极品女厕沟 激情五月bd韩国电影 混血美女自摸和男友激情啪啪自拍诱人呻吟福利视频 人人摸人人妻做人人看 44kknn 娸娸原网 伊人欧美 恋夜影院视频列表安卓青青 57k影院 如果电话亭 avi 插爆骚女精品自拍 青青草在线免费视频1769TV 令人惹火的邻家美眉 影音先锋 真人妹子被捅动态图 男人女人做完爱视频15 表姐合租两人共处一室晚上她竟爬上了我的床 性爱教学视频 北条麻妃bd在线播放版 国产老师和师生 magnet wwwcctv1024 女神自慰 ftp 女同性恋做激情视频 欧美大胆露阴视频 欧美无码影视 好女色在线观看 后入肥臀18p 百度影视屏福利 厕所超碰视频 强奸mp magnet 欧美妹aⅴ免费线上看 2016年妞干网视频 5手机在线福利 超在线最视频 800av:cOm magnet 欧美性爱免播放器在线播放 91大款肥汤的性感美乳90后邻家美眉趴着窗台后入啪啪 秋霞日本毛片网站 cheng ren 在线视频 上原亚衣肛门无码解禁影音先锋 美脚家庭教师在线播放 尤酷伦理片 熟女性生活视频在线观看 欧美av在线播放喷潮 194avav 凤凰AV成人 - 百度 kbb9999 AV片AV在线AV无码 爱爱视频高清免费观看 黄色男女操b视频 观看 18AV清纯视频在线播放平台 成人性爱视频久久操 女性真人生殖系统双性人视频 下身插入b射精视频 明星潜规测视频 mp4 免賛a片直播绪 国内 自己 偷拍 在线 国内真实偷拍 手机在线 国产主播户外勾在线 三桥杏奈高清无码迅雷下载 2五福电影院凸凹频频 男主拿鱼打女主,高宝宝 色哥午夜影院 川村まや痴汉 草溜影院费全过程免费 淫小弟影院在线视频 laohantuiche 啪啪啪喷潮XXOO视频 青娱乐成人国产 蓝沢润 一本道 亚洲青涩中文欧美 神马影院线理论 米娅卡莉法的av 在线福利65535 欧美粉色在线 欧美性受群交视频1在线播放 极品喷奶熟妇在线播放 变态另类无码福利影院92 天津小姐被偷拍 磁力下载 台湾三级电髟全部 丝袜美腿偷拍自拍 偷拍女生性行为图 妻子的乱伦 白虎少妇 肏婶骚屄 外国大妈会阴照片 美少女操屄图片 妹妹自慰11p 操老熟女的b 361美女人体 360电影院樱桃 爱色妹妹亚洲色图 性交卖淫姿势高清图片一级 欧美一黑对二白 大色网无毛一线天 射小妹网站 寂寞穴 西西人体模特苍井空 操的大白逼吧 骚穴让我操 拉好友干女朋友3p