Front. Microbiol. Frontiers in Microbiology Front. Microbiol. 1664-302X Frontiers Media S.A. 10.3389/fmicb.2019.00102 Microbiology Original Research Adaptive Strategies in a Poly-Extreme Environment: Differentiation of Vegetative Cells in Serratia ureilytica and Resistance to Extreme Conditions Filippidou Sevasti 1 2 3 Junier Thomas 4 Wunderlin Tina 1 Kooli Wafa M. 1 Palmieri Ilona 1 Al-Dourobi Andrej 1 Molina Veronica 5 Lienhard Reto 6 Spangenberg Jorge E. 7 Johnson Shannon L. 8 Chain Patrick S. G. 8 Dorador Cristina 9 10 Junier Pilar 1 * 1Laboratory of Microbiology, University of Neuchatel, Neuchatel, Switzerland 2Microbial Ecology Group, Centre for Ecology and Hydrology, Wallingford, United Kingdom 3Space Microbiology Research Group, Radiation Biology Department, Institute of Aerospace Medicine, German Aerospace Center (DLR e.V.), Cologne, Germany 4Vital-IT Group, Swiss Institute of Bioinformatics, Lausanne, Switzerland 5Departamento de Biología, Facultad de Ciencias Naturales y Exactas, Universidad de Playa Ancha, Valparaíso, Chile 6ADMED Microbiologie, La Chaux-de-Fonds, Switzerland 7Institute of Earth Surface Dynamics, University of Lausanne, Lausanne, Switzerland 8Los Alamos National Laboratory, Los Alamos, NM, United States 9Laboratorio de Complejidad Microbiana y Ecología Funcional, Departamento de Biotecnología, Facultad de Ciencias del Mar y Recursos Biológicos, Universidad de Antofagasta, Antofagasta, Chile 10Centre for Biotechnology and Bioengineering, Antofagasta, Chile

Edited by: Mónica Sánchez-Román, VU University Amsterdam, Netherlands

Reviewed by: Haitham Sghaier, National Center for Nuclear Sciences and Technology, Tunisia; Surendra Vikram, University of Pretoria, South Africa; James A. Coker, University of Maryland University College, United States

*Correspondence: Pilar Junier, pilar.junier@unine.ch

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

05 02 2019 2019 10 102 14 06 2018 17 01 2019 Copyright © 2019 Filippidou, Junier, Wunderlin, Kooli, Palmieri, Al-Dourobi, Molina, Lienhard, Spangenberg, Johnson, Chain, Dorador and Junier. 2019 Filippidou, Junier, Wunderlin, Kooli, Palmieri, Al-Dourobi, Molina, Lienhard, Spangenberg, Johnson, Chain, Dorador and Junier

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.

Poly-extreme terrestrial habitats are often used as analogs to extra-terrestrial environments. Understanding the adaptive strategies allowing bacteria to thrive and survive under these conditions could help in our quest for extra-terrestrial planets suitable for life and understanding how life evolved in the harsh early earth conditions. A prime example of such a survival strategy is the modification of vegetative cells into resistant resting structures. These differentiated cells are often observed in response to harsh environmental conditions. The environmental strain (strain Lr5/4) belonging to Serratia ureilytica was isolated from a geothermal spring in Lirima, Atacama Desert, Chile. The Atacama Desert is the driest habitat on Earth and furthermore, due to its high altitude, it is exposed to an increased amount of UV radiation. The geothermal spring from which the strain was isolated is oligotrophic and the temperature of 54°C exceeds mesophilic conditions (15 to 45°C). Although the vegetative cells were tolerant to various environmental insults (desiccation, extreme pH, glycerol), a modified cell type was formed in response to nutrient deprivation, UV radiation and thermal shock. Scanning (SEM) and Transmission Electron Microscopy (TEM) analyses of vegetative cells and the modified cell structures were performed. In SEM, a change toward a circular shape with reduced size was observed. These circular cells possessed what appears as extra coating layers under TEM. The resistance of the modified cells was also investigated, they were resistant to wet heat, UV radiation and desiccation, while vegetative cells did not withstand any of those conditions. A phylogenomic analysis was undertaken to investigate the presence of known genes involved in dormancy in other bacterial clades. Genes related to spore-formation in Myxococcus and Firmicutes were found in S. ureilytica Lr5/4 genome; however, these genes were not enough for a full sporulation pathway that resembles either group. Although, the molecular pathway of cell differentiation in S. ureilytica Lr5/4 is not fully defined, the identified genes may contribute to the modified phenotype in the Serratia genus. Here, we show that a modified cell structure can occur as a response to extremity in a species that was previously not known to deploy this strategy. This strategy may be widely spread in bacteria, but only expressed under poly-extreme environmental conditions.

differentiation dormancy extreme environments harsh conditions Serratia sporulation resistance Schweizerischer Nationalfonds zur Förderung der Wissenschaftlichen Forschung10.13039/501100001711

香京julia种子在线播放

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

      Introduction

      The Biospace corresponds to the physicochemical conditions under which life can be sustained (Cockell, 2015). On Earth, both mesophilic and extremophilic habitats are part of the Biospace, but only extreme habitats are considered as analogs to Young Earth or other planets and exoplanets, making them suitable sites to study the origin and evolution of life, as well as the extraterrestrial Biospace. Extreme habitats are considered challenging ecosystems for life because physicochemical parameters deviate from those used to describe mesophilic conditions sustaining life. However, on Earth, a diverse microbial community inhabit extreme ecosystems. This is supported, for example, by diversity surveys on terrestrial and oceanic geothermal sites (Schubotz et al., 2013), saline environments [salt lakes, seashore evaporations, or estuaries; (Edwards, 1990)], acidic and alkaline environments (Edwards, 1990; Sorokin et al., 2014) and ecosystems that are highly radiated (Sghaier et al., 2007; Cordero et al., 2018). Although based on a biased perception of extremity we tend to consider extreme habitats as more challenging for survival, fluctuations in environmental parameters even in non-extreme habitats, would lead to periods in which active growth is compromised (Sadoff, 1973).

      In order to withstand adversity, life adapts by deploying various survival strategies, one of which is cellular differentiation into a more resistant form. In eukaryotes, cells containing the same genetic information can transform into a myriad of differentiated cell types [for instance in the infamous case of stem cell differentiation to a broad range of specialized cells (indicative reference: Engler et al., 2006)]. In contrast to the complexity and cellular diversity of eukaryotes, prokaryotic cells are perceived as possessing simpler morphologies. However, they do differentiate into various resistant forms, such as endospores (Firmicutes), myxospores (δ-Proteobacteria), akinetes (Cyanobacteria), and exospores (Actinobacteria) (Barton, 2005). These differentiated structures offer the advantage of resisting through adverse times under resting conditions (low metabolic activity), with the possibility to deploy full metabolic activity and cellular division once conditions become favorable.

      Reports of dormancy involving specialized cellular structures outside these four groups are also found in literature. These include the proposed formation of endospores in the purple non-sulfur phototroph α-Proteobacterium Rhodobacter johrii JA192, which produced light-refracting structures visible under phase contrast microscopy and malachite green staining (Girija et al., 2010). Likewise, cyst formation has been shown in Rhodospirillum centenum, a photosynthetic member of the Azospirillum clade with a complex developmental life cycle (Dong and Bauer, 2015). Another example, which has been at the center of recent debate, was the production of phase-bright granules similar to endospores in Mycobacteria (Ghosh et al., 2009). Although data of mRNA expression for some key sporulation-related genes was provided in the case of Mycobacteria (Ghosh et al., 2009), subsequent genomic analysis revealed the absence of the core set of sporulation genes identified for Firmicutes and production of phase-bright structures was irreproducible by other researchers (Traag et al., 2010). This latter example, shows some of the challenges of describing sporulation in a clade previously considered as asporogenic. Formation of endospores has been also suggested for Serratia (Ajithkumar et al., 2003). The genus Serratia belongs to γ-Proteobacteria, to the family Enterobacteriaceae, and some of its members are of clinical importance (Grimont and Grimont, 2015). Production of concentric “spores” in response to heat-shock was reported in Serratia marcescens subsp. sakuensis (Ajithkumar et al., 2003). This strain was isolated from a wastewater treatment tank, alongside several strains belonging to the genus Bacillus, leading the authors to speculate at horizontal gene transfer to explain the origin of sporulation in Serratia. In spite of these reports, the formation of dormant cells outside the four thoroughly studied taxa have always been received with skepticism, and in particular in the case of endospore formation, which is so far the best studied example of a dormant cell in Bacteria. Endospores arise from asymmetrical cell division and a complex morphogenetic differentiation program (Way, 1996). Because of its complexity, endospore formation is supposed to be restricted to Firmicutes (Traag et al., 2010). In the case of other differentiated dormant cells such as exospores, myxospores or cysts, the mechanisms of formation are much less understood (Abel-Santos, 2012), limiting our ability to perform comparative studies in order to assess their prevalence.

      Over the last three decades, the formation of cyst-like resting structures has been proposed for diverse non-spore-formers (Sadasivan and Neyra, 1985; Suzina et al., 2004). These cyst-like cells share the characteristics of refracting light in the phase contrast microscope, survival over long periods of time, and revival under optimal conditions (Muliukin et al., 1997). Although these reports expand our understanding of dormancy as a survival strategy in response to conditions limiting active growth, the current knowledge on environmental triggers for the production of these cyst-like resting cells is mainly limited to nutrient deprivation. Moreover, the resistance limits of these structures and the genetic background of this process are not yet studied.

      In this study, we report the discovery of cyst-like cells in S. ureilytica str. Lr5/4, which was isolated in a geothermal spring in Lirima, Atacama Desert, Chile. We provide evidence that this cyst-like cell form is more resistant than the vegetative cell it derived from. The environmental factors triggering the formation of cyst-like cells, the resistance of these structures and their biochemical properties, as well as the genomic imprints that could be related to cyst production are studied herein. Moreover, whether the phenomenon of resting cells in non-spore forming bacteria is widespread remains unanswered. In order to address this question, phylogenomic and functional analysis of genes related to resting structures was performed. Overall, with this work we attempt to shed light on the survival behavior of cyst-forming mesophilic bacteria in extreme environments and discuss the ecological relevance of these bacteria.

      Materials and Methods Sampling and Isolation

      Lirima is a geothermal site located in the Andean highlands of northern Chile (3.997 m.a.s.l., 19°51.118W, 68°54.402S). The Chilean Altiplano presents multiple co-occurring environmental stressors: the highest UV radiation ever reported on Earth (Cordero et al., 2016); variable daily temperatures (from below 0°C up to 30°C); high altitude (4000 m); and low relative humidity (Cordero et al., 2016). Microbial mats in different geothermal environments (aqueous and non-aqueous, temperature ranging from 30°C to over 90°C) were sampled during a field campaign in April 2011. The microbial mats were collected and stored using sterile material. After sampling, the sample was preserved at 4°C until enrichment. For aerobic enrichment, one gram of homogenized mat was inoculated in nutrient broth (NB) (Biolife, Italy), pH adjusted to 7.2, and incubated at room temperature for 7 days under aerobic conditions. The mismatch between in situ and cultivation temperature was used to favor the selection of dormant forms. 100 μL of the enrichment culture was plated on nutrient agar (NA) (Biolife, Italy) and incubated for 24 h at room temperature (RT) under aerobic conditions. Axenic cultures were isolated after repeated serial dilutions on NA medium at RT for 24 h. Strains were confirmed for purity by microscopic observation. Pure isolates were cryopreserved in 30% (v/v) of glycerol at -80°C. Isolates were screened for phase bright spore-like structures by microscopy after having been selected for dormancy by nutrient deprivation. The isolate Lr5/4 described here was isolated from a double layer microbial mat in a source pond at 54°C (Figure 1), alongside four other strains that were identified as belonging to the genus Bacillus.

      (A) The sampling site of the geothermal pond in Lirima, Chile, from which S. ureilytica Lr5/4 was isolated. (B) A close-up photograph of the microbial mat.

      Physiological Characterization of Lr5/4

      Cell growth was monitored at different temperatures (4, 15, 25, 30, 37, 45, 55, and 60°C) for up to 72 h in nutrient broth medium (Biolife, Italy), measuring optical density at 600 nm (OD600) with a spectrophotometer Genesys 10S UV-Vis Thermo Scientific (Thermo Fisher Scientific Inc., United States). In addition, growth was measured at different pH (2 to 13) and at different NaCl concentrations (1, 2, 3, 3.5, 4, 5, 5.5, 6, 6.5, 8, 9, 10, 11, and 12% w/v) with a microplate reader at 590 nm (UVM 340 Microplate Reader, ASYS, Hiteck, United Kingdom) for up to 72 h. The need for oxygen during growth was verified using the method of thioglycollate medium (Brewer, 1940).

      Phospholipids were extracted from the whole cells using a single-phase (methanol/dichloromethane/phosphate buffer; 2:1:0.8; v/v/v) extraction procedure (Guckert et al., 1985). Fatty acid methyl esters (FAMES) were prepared from the extracted acyl lipid mixture by acid-catalyzed methanolysis (Ichihara and Fukubayashi, 2010). The chemical characterization of the lipids was performed with the Agilent gas chromatograph 6890 coupled to the Agilent 5973 quadrupole mass selective detector (GC-MS; Palo Alto, CA, United States). For the analysis of FAME fractions, the system was equipped with the Agilent free fatty acids phase (FFAP) fused silica capillary column (50 m length, 0.20 mm i.d.) coated with nitroterephthalic acid modified polyethylene glycol stationary phase (film thickness 0.33 μm). A sample aliquot was injected split-less at a temperature of 200°C. Helium was used as carrier gas (1 mL/min flow rate). After an initial period of 2 min at 100°C, the column was heated to 240°C at 5°C/min followed by an isothermal period of 30 min. The MS was operated in the electron impact mode at 70 eV, source temperature of 250°C, emission current of 1 mA and multiple-ion detection with a mass range from 50 to 600 amu. Compound identifications were based on comparison of standards, GC retention time, and mass spectrometric fragmentation patterns.

      Molecular Characterization

      For G+C content analysis and DNA-DNA hybridization, the services of DSMZ, Germany were used. Strain Lr5/4 was tested against S. marcescens (DSMZ 30121) and S. ureilytica (DSMZ 16952). MALDI-TOF mass spectra were acquired on a Bruker Microflex RLF mass spectrometer (Bruker Daltonics, Bremen, Germany). Spectral mass resolutions and signal-to-noise ratios were determined with the software Flex Analysis 3.3.65 (Bruker Daltonics).

      For PCR amplifications, genomic DNA was extracted using the InnuPREP Bacteria DNA Kit (Analytik Jena, Germany), according to the manufacturer’s instructions. DNA was quantified fluorometrically with a Qubit® dsDNA HS Assay Kit with a Qubit® 2.0 Fluorometer (Invitrogen Ltd., Paisley, United Kingdom) according to the manufacturer’s instructions. PCR amplification of the 16S rRNA gene was performed using the primer set GM3F and GM4R (Muyzer et al., 1995). The PCR products were purified with a MultiScreen PCRμ96 Filter Plate (Millipore), according to the manufacturer’s instructions and sequenced using the services of GATC Biotech (Germany). To obtain the full-16S rRNA gene sequence, Lr5/4 PCR products were sequenced in addition with the primers 907r, 520r, and 926f primers (Muyzer et al., 1995; Forney et al., 1997). The 16S rRNA gene sequence was compared in GenBank by using the BLASTn tool (Altschul et al., 1997). After alignment, a series of sequences that had an identity match over 98% were selected, along with 16S rRNA gene sequences from cultured representative strains belonging to different species of the genus Serratia and a 16S rRNA gene sequence of a Geobacillus sp. (serving as an outgroup in order to root the tree). These sequences were used to build-up a cladogram using the online tools of phylogeny.fr website (Dereeper et al., 2008). Multiple sequence alignment has been performed using MUSCLE and the phylogenetic tree has been constructed using PhyML. To determine the confidence values for individual branches, 100 bootstrap replications were calculated for each generated tree.

      Morphological Characterization

      Colonies of strain Lr5/4 were obtained and observed after overnight growth on NA. Gram staining was performed on an overnight solid culture using the Hucker staining method (Gerhardt, 1994). Spores were observed with a phase contrast microscope (Leica DM R, magnification 1000×). Vegetative cells and spores were also observed by Scanning and Transmission Electron Microscopy (SEM and TEM). Vegetative cells were observed from a fresh 24 h culture inoculated in NB. Cyst-like structures were collected from a 3-month old solid culture. Both preparations were fixed in 2.5% glutaraldehyde in a cacodylate buffer (0.1 M; pH7.4) for 2 h at room temperature and then overnight at 4°C. They were then washed by gentle immersion in cacodylate buffer (0.2M; pH7.4) post fixed with 1% OsO4 in the same buffer, and carefully washed with the above buffer. For SEM, the samples were dehydrated in 15 to 100% of ethanol solution and finally fixed on Poly-L-Lysine slides and coated with a 23 nm gold layer in a BaltecSCD005 sputter apparatus. The samples were observed with a Philips XL30 SEM at acceleration voltages of 10–20 kV. For TEM, the samples were dehydrated in 15–100% acetone and embedded in Spurr’s resin. Serial sections were made with a Reichert Ultracut-S microtome, mounted on copper grids, double stained with uranyl acetate and lead citrate, and observed with a Philips CM 100 TEM at 60 or 80 kV.

      Dormant Structure Production and Resistance Assays

      Nutrient deprivation on solid media was one method for the production of light-refracting spherical cysts. Other environmental triggers were also put to test. Overnight liquid cultures in NB were subjected to temperature shocks. These shocks included two rounds of placing the liquid culture in an 80°C water bath for 20 min, then moved to 4°C for 20 min. The culture was then placed at room temperature overnight and observations for light refracting spheres was observed at 24 or 48 h. Overnight solid cultures on NA were subjected to desiccation and UV radiation, in order to verify whether these factors trigger the formation of cyst-like bodies. Overnight solid cultures were placed in a desiccator for 72 h until the solid medium dried to a 1 mm thin layer. Formation of light refracting spheres was examined at 24 and 48 h. Overnight solid cultures were also placed 60 cm from an UVc lamp (approximate wavelength 254 nm). The experiment was held under sterile conditions (BSL II hood), and the lids of the petri dishes were removed so that the bacterial colonies were directly exposed to the radiation. At t = 5, 10, 20, 30, 60, 120, 180, 240, and 720 min, three petri dishes were removed, closed sterilely and kept at room temperature. Observations for light refracting spheres was made at 24 or 48 h. Finally, overnight liquid cultures in NB were autoclaved and subjected to temperature extremes (no heat-shock). The extreme temperatures tested were -80, -20, 4, 60, 80, and 100°C. The overnight liquid cultures were incubated for 24 h at these temperatures, then observed for cyst formation. S. ureilytica str. Lr5/4 was tested along with Bacillus subtilis str. NEU1294, Anoxybacillus sp. str. Lr10/3 (isolated along Lr5/4 from Lirima), S. marcescens DSMZ 30121, S. ureilytica DSMZ 16952, and Escherichia coli str. NEU1006.

      Spores produced by nutrient deprivation and heat-shocks were subjected to resistance tests. A heat resistance test for the spores of strain Lr5/4 was performed as previously described (Ajithkumar et al., 2003), along with S. marcescens (DSMZ 30121), S. ureilytica (DSMZ 16952), Bacillus subtilis (Neu1294), and Anoxybacillus sp. (Lr10/3). Suspensions of these cultures in tryptic soy broth (TSB) medium were heat-shocked at 70 and 75°C for 20 min and then re-cultured under optimal conditions. Growth was macroscopically verified after 24 and 48 h. Cyst-like or spore preparations, as well as vegetative cells from old S. marcescens str. DSMZ 30121, S. ureilytica str. DSMZ 16952, and E. coli str. NEU1006, were also autoclaved and exposed to UV radiation (radiation applied at 60 cm from a UV-C lamp, 200–280 nm) and to temperature extremes (as described above for cyst production). After exposure, biomass from the solid or liquid medium was stroked on solid media and cultured under optimal conditions in order to verify re-growth.

      DPA Measurement

      The presence of dipicolinic acid (DPA) in the spores was assessed on a 3-month old Lr5/4 preparation (20 mg wet weight), according to a previously published method (Brandes Ammann et al., 2011). This preparation contained almost exclusively refracting light structures. Fluorescence was measured within a Perkin-Elmer LS50B fluorometer. The excitation wavelength was set at 272 nm with a slit width of 2.5 nm. Emission was measured at 545 nm (slit width 2.5 nm). The device was set in the phosphorescence mode (equivalent to time-resolved fluorescence). The delay between emission and measurement was set at 50 μs. Measurements were performed every 20 ms. The integration of signal was performed over a duration of 1.2 ms. Values recovered for each measurement corresponded to the mean of the relative fluorescence unit (RFU) values given by the instrument within the 30 s following sample introduction in the device. Finally, to transform RFU units into DPA concentrations, a 10-point standard curve was established using increasing concentrations of DPA from 0.5 up to 10 μM.

      Genome Sequencing

      Genomic DNA was extracted from an overnight culture using the Genomic-tip 20/G Kit (Qiagen GmbH, Germany). Sequencing was performed with PacBio RS II system based on single molecule, real-time (SMRT) technology (Pacific Biosciences, California). Genome assembly was performed using hgap 2.0 in the SmrtPortal. The full closed genome of S. ureilytica strain Lr5/4 presents a unique contig of 5,39 Mbp in size, and a G+C content of 59.2%. Genome annotation utilized an Ergatis based (Hemmerich et al., 2010) workflow with minor manual curation. The circular genome with the additional features was created using DNAplotter (Carver et al., 2009). The genome sequence was submitted to NCBI under Bioproject accession number PRJNA260750.

      Sporulation-Related Genes Profiles

      An extensive literature search was performed in order to compile all spore- and sporulation-specific genes reported for the known spore-forming phyla (Actinobacteria, Cyanobacteria, Firmicutes, and Proteobacteria). Whether a protein encoded by those genes was submitted to UniProt or NCBI was then verified. These proteins, along with their homologs from other spore-forming species were downloaded and multiple sequence alignments were performed using Clustal Omega, with default parameters (Sievers et al., 2011). For these proteins, profiles were built using HMMER version 3.1 with default parameters (Eddy, 2011). An additional extended literature search was performed in order to catalog all reported spore-forming and non-spore-forming species of the four phyla (Actinobacteria, Cyanobacteria, Firmicutes, and Proteobacteria). The genomes of the reported spore-formers were downloaded from NCBI, when available. These genomes were thus screened using the generated profiles and the latter were updated. The updated profiles were used to screen the genome of S. ureilytica Lr5/4.

      Genome Analysis Comparative Genomics

      Serratia ureilytica Lr 4/5 full genome was compared to other genomes of Serratia species using the online comparative genomics tools EDGAR (Blom et al., 2009). A phylogenetic tree based on this comparison was generated. To construct a phylogenetic tree comprising several genomes (NCBI Accession numbers provided in Supplementary Figure 2) the core genes of these genomes were established. Alignments of the core genes were generated using MUSCLE (Edgar, 2004). Non-matching parts of the alignment were masked by GBLOCKS (Castresana, 2000) and subsequently removed. The remaining parts of the alignments were concatenated to one unique alignment used as input for PHYLIP (Felsenstein, 1993). This led to a phylogenetic tree, represented in Newick format.

      Presence of Sporulation Genes

      The genome of S. ureilytica Lr5/4 was screened for orthologs against the four above-mentioned datasets, using tBLASTn, under default parameters and an E-value cutoff of 1e-11. A tBLASTn run on the shuffled protein sequences as a negative control set showed no hit with an E-value lower than 4e-4. Each sporulation-related S. ureilytica Lr 4/5 gene was subjected to the following procedure:

      Finding homologs

      We searched for sequence-level similarities in UniProtKB using BLASTX on the UniProtKB website and default parameters. Then, using the best BLASTX hit, we considered (in this order): (1) links to orthology or phylogenomics databases; (2) links to databases of protein domain families. Links to orthology and phylogenomics databases were not available. Therefore, protein domain data was used instead. In this case, we selected the longest domain that matched our query, and downloaded Pfam (Finn et al., 2014) seed alignments for the domain. This resulted in sequences that were probably homologous to our query gene (since they contained the matching domain), but, in most cases, not necessarily orthologous as a since paralog will contain the domain as well.

      Building alignments

      The query sequence was added to the Pfam seed alignment using Mafft (Katoh et al., 2002), using the option –keep-aligned and xyz so as to keep the original alignment intact and only aligning the query to it.

      Computing and rooting trees

      Using these alignments as inputs, maximum-likelihood phylogenetic trees were computed with PhyML (Le and Gascuel, 2008), using the XYZ model, maximum-likelihood estimates of amino acid and invariant frequencies, 12 substitution rate categories, and best topology search moves (-d aa -f a -v e -c 12 -a e -s BEST). In the general case, it was not possible to determine a suitable outgroup. This is because (i) in most cases the trees contain paralogs, and (ii) the phylogeny of bacteria is still unresolved, which entails that only a eukaryotic or archaeal homolog would be suitable – and these do not always exist. Therefore, we did not attempt to root the trees.

      Closest relatives

      Instead, we used the patristic distance (i.e., distance along tree branches) to infer the closest relative to the genes identified in S. ureilytica Lr 5/4. This differs from a simple genetic distance matrix in that the distances were actually estimated by maximum likelihood over a phylogeny and are therefore arguably more accurate. The distances were extracted with the Newick Utilities (Junier and Zdobnov, 2010).

      Heat maps

      The heat maps were produced with the heatmap() function from the R package.

      Results Isolation of <italic>Serratia ureilytica</italic> Lr5/4

      The characterization of the strain Lr5/4 showed oval cells (0.5 × 1–3 μm) with a negative Gram-staining. Growth tests showed that Lr5/4 is a mesophile (temperature growth range between 10 and 37°C) with an optimum at 25°C. This range is not within the temperature measured at the sampling site (56°C), suggesting that the strain was present in the geothermal spring in a dormant state. Isolate Lr5/4 is also moderately halotolerant (growth up to 8% NaCl), able to develop over a vast pH range (from 3 to 11; optimal 5 to 6), and a facultative anaerobe, which are all characteristics that can be expected from a poly-extremophilic bacterium inhabiting the Chilean Altiplano. Based on the entire 16S rRNA gene sequence, isolate Lr5/4 belongs to the genus Serratia, with more than 97% identity to both Serratia marcescens and Serratia ureilytica. To establish the phylogenetic affiliation of the strain, a series of additional analysis were performed. While, MALDI-TOF characterization suggested that strain Lr5/4 belonged to S. marcescens (Supplementary Figure 1), PFLA analysis (Supplementary Table 1) showed that palmitic acid is the major fatty acid in the three strains. Based on the total fatty acid composition, Lr5/4 is more closely related to S. ureilytica, although the fatty acid profile of S. ureilytica type strain correlates better with S. marcescens than with Lr5/4. DNA/DNA hybridization analysis between Lr5/4 and its closest relatives (S. ureilytica DSM-16952 and S. marcescens DSM-30121) suggests that it is a member of S. ureilytica, with above 99% DNA/DNA relatedness to S. ureilytica DSM-16952. Nonetheless, DNA/DNA relatedness was also found to be above 70% with S. marcescens subsp. marcescens, which is the recommended threshold for delineating bacterial species (18). This observation is intriguing as it has been previously shown that DNA/DNA relatedness between the two reference species (S. ureilytica and S. marcescens subsp. marcescens) is only 43.7% (19). Finally, a phylogenetic tree constructed using all the available full genomes of Serratia isolates was constructed showing the relationship between the species that belong to this genus (Supplementary Figure 2). Considering all biochemical and molecular tests performed, we concluded that strain Lr5/4 represents a novel strain of Serratia ureilytica.

      Morphology of Dormant Cells in <italic>S. ureilytica</italic> Lr5/4

      After nutrient deprivation, S. ureilytica Lr5/4 differentiated into elliptic-spherical cell forms that appeared phase-bright at the phase contrast microscope (Figure 2A). The morphological characterization of dormant cells in strain Lr5/4 was further conducted using electron microscopy. Using scanning electron microscopy, vegetative cells display a bacillus-like shape of approximately 2.0 μm × 0.5 μm in size, while dormant cells corresponded to round spheres of 0.5 μm in diameter (Figure 2B). Using TEM and osmium (OsO4) fixation, the interior of vegetative cells appears stained while dormant cells present a clearer core surrounded by two thin outer layers (Figure 2C and inset). Although the exact composition of these two layers is not yet known, they prevent the penetration of the osmium tetraoxide (OsO4), used for fixation. The latter penetrated easily the cell wall and membrane of the vegetative form, which became dark gray (Figure 2C). OsO4 is used for fixation of spore structures of Firmicutes, Actinomycetes (Lechevalier et al., 1966), and Myxococcus (Müller et al., 2015). An increase in coat thickness hinders the penetration of OsO4, resulting in poor resolution of the inner coat and the spore protoplasm (Stevenson et al., 1972). The same has also been shown during germination experiments. In the latter, the protoplasm of spores is brighter and stains darker during germination as the spore coat is remodeled during outgrow (Hoeniger and Headley, 1968). Finally, experiments on the formation of resting light-refracting structures previously performed (Muliukin et al., 1997), have shown that TEM pictures of cells that resume their growth cycle stain dark, as observed also for S. ureilytica str. Lr5/4 (Figure 2C). This suggests that in spite of the apparent simplicity of the outer layers in the dormant form of strain Lr5/4, they are nonetheless impermeable to OsO4 and prevent staining. Overall the morphological characteristics of the dormancy cells in Serratia suggest that they differ in size and form to the vegetative cells.

      (A) Phase contrast microscopy of a 2-month old culture of Lr5/4. Phase bright objects of different sizes, are surrounded by lysed cells. Bar = 10 μm. (B) Scanning electron microscopy of vegetative (V), lysed (D), and dormant (S) cells. (C) Transmission electron microscopy (TEM) of vegetative (V), lysed (D), and dormant (S) cells. Noticeable black spots inside the lysed cells. Dormant cells are clearly shown to be double-membraned (see also inset).

      Dormancy Triggers

      Although many environmental factors can hinder bacterial survival, nutrient deprivation is a well-studied dormancy triggering factor (Barton, 2005). Accordingly, S. ureilytica Lr5/4 dormant cells developed in response to starvation (after 3 months of growth on agar medium). Other environmental factors may also influence spore-formation, as well as resting cell production in non-spore-formers (Suzina et al., 2006), however, these factors have not been thoroughly tested in previous studies. Therefore, for this study, the roles of other known triggers of dormancy that could be relevant in the geothermal spring from which strain Lr5/4 was isolated (UV radiation, thermal shock, temperature extremes and desiccation) were also measured. Dormant cells of Lr5/4 were produced in response to UV radiation and thermal shock (Supplementary Table 2). In contrast, cells of Lr5/4 were largely tolerant to desiccation. Extreme temperatures (down to -80°C and up to 100°C applied overnight), did not induce dormancy in S. ureilytica Lr5/4. B. subtilis str. NEU1294 and Anoxybacillus sp. str. Lr10/3 responded to desiccation and extreme temperatures by producing endospores. On the contrary, S. marcescens DSMZ 30121, S. ureilytica DSMZ 16952, and E. coli str. NEU1006 did not produce any cyst-like cells. Although cyst-like resting structures have been previously reported for E. coli (Suzina et al., 2004), it appears that the specific environmental triggers that are characteristic of the Atacama Desert did not trigger cyst-production in the E. coli strain tested here. No spore or cyst production was observed in response to autoclaving. Autoclaving killed all cells tested.

      Resistance of Lr5/4 Dormant Cell

      In order to evaluate the resistance of the dormant forms produced by S. ureilytica Lr5/4, cyst-like cells from 90-day old and heat-shocked cultures were challenged with multiple environmental stressors. Dormant forms of Lr5/4 were able to withstand and generate an actively growing culture after wet heat-shock at 75°C, 240 min UV radiation exposure, and incubation at extreme temperatures. This was comparable to the survival of endospores from the endospore-forming bacilli B. subtilis and Anoxybacillus sp. Lr10/3, although the latter resisted to 120 min UV radiation exposure (Supplementary Table 3). However, cells of S. marcescens DSMZ 30121 and S. ureilytica DSMZ 16952, as well as cells of E. coli did not survive these treatments. Autoclaving killed all spores and cells tested. When dormant forms of Lr5/4, B. subtilis and Anoxybacillus sp. Lr10/3 were subjected to multiple stressors applied simultaneously (starvation, desiccation, and UV radiation), only the endospores of the Firmicutes species could be revived, while Lr5/4 did not survive.

      Although for most of the stressors studied, the precise molecular mechanism involved in resistance is unknown, in endospores of Firmicutes the protection of DNA in response to wet heat is conferred by the accumulation of dipicolinic acid (DPA) in the core (Setlow et al., 2006). The presence of DPA (Brandes Ammann et al., 2011) in the dormant cells of Lr5/4 was assessed on dilutions of a 3-month old culture (Supplementary Figure 3). The DPA content per dormant cell was calculated to be between 2.36 and 2.55 × 10-7 μg/dormant cell, a value nearly 14 times lower than the DPA concentration in spores of B. subtilis. DPA was also reported as detectable in the case of S. marcescens subsp. sakuensis [10 mg (g bacteria)-1] (Ajithkumar et al., 2003), but the values cannot be compared to those obtained here as the dormant to vegetative cell ratio was not provided in the case of S. marcescens subsp. sakuensis.

      Genomic Imprints of Dormancy in Lr5/4

      In order to investigate the genomic imprints of dormancy in strain Lr5/4, the complete genome of the strain was sequenced (5.39 Mb in size, 59.2% G+C) and annotated. We next compared the 5,056 annotated genes to genes reported to be involved in sporulation in other clades. The molecular pathways for bacterial cellular differentiation to the formation of a spore are well described in Actinobacteria, Cyanobacteria, Firmicutes, and Proteobacteria. Therefore, a database with all the known genes involved in dormancy in these taxa (Supplementary Table 4) was built and the genome of S. ureilytica Lr5/4 was screened for the presence of those genes. The detection of these genes using a BLAST search is indicative rather than demonstrative of the presence of this genetic determinants and further experiments, including gene expression and functional characterization during sporulation, need to be performed. Therefore, these results should be regarded with caution. Nonetheless, significant homology was detected between 12 Lr5/4 loci and myxospore-formation genes, nine in actinobacterial exospores, and 24 genes involved in spore formation in the Firmicutes (Supplementary Table 5). These genes participate in peptidoglycan synthesis, signaling, and transcription regulation. These genes are dispersed throughout the genome of Lr5/4 (Figure 3). The phylogenetic analysis of the genes that are potentially involved in dormancy show that most of them are likely of proteobacterial origin (Figure 4), suggesting they have not recently been acquired via horizontal gene transfer. A thorough phylogenomic and functional analysis of the dormancy-related genes showed that they contain conserved domains homologous to proteins that, in other species, participate mainly to general stages of cell division and sensing (Table 1), with homologs reported also in non-dormant bacteria.

      Circular representation of the S. ureilytica Lr5/4 genome. Genes related to sporulation are spread throughout the chromosome. Genes indicated in red show homology to genes related to endospore-formation in Firmicutes, in blue to myxospore-formation in delta-Proteobacteria and in green to exospores in Actinobacteria. The inner circle in yellow (A+T%) and purple (G+C%) represents the G+C% content variation throughout the bacterial chromosome.

      Heatmaps showing the phylogenomic neighborhood for four genes related to sporulation found in the genome of S. ureilytica Lr5/4 (highlighted in red). In panels (A–C), the highest similarity for the sporulation genes lies with phylogenetically closely related taxa (indicated in red). This was true for 25 genes. In panel (D), the only exception is shown. Homologs of the lonB gene are more similar between species from distant phyla, than from species of the same phylum. Strain Lr5/4 has a lonB gene that is found to bear closer similarity to its homolog in Myxococcus (green) and Bacillus (purple), rather than to homolog genes in γ-Proteobacteria.

      Presence of the 26 genes detected in the genome of Serratia in multiple phyla.

      A, actinomycetes; C, cyanobacteria; F, firmicutes; P, proteobacteria.
      Discussion

      In this work we report the formation of cyst-like cells in a novel Serratia strain. We show for the first time that S. ureilytica Lr5/4 produces a differentiated cell form that is more resistant than vegetative cells. This is probably due to its denser cellular structure. Furthermore, to date, only the vegetative form has been observed to divide under optimal conditions. When the differentiated cell types are cultured under optimal conditions, only the vegetative cell forms are observed under the microscope. In dormant state, cells contain dipicolinic acid, a compound that protects the genetic material from degradation upon exposure to wet heat (Setlow et al., 2006). For all these reasons, we argue that the differentiated cell type of Serratia Lr5/4 is a resting resistant cellular structure. Our results show that cellular differentiation occurs under starvation, strong temperature fluctuations (shocks) and UV radiation. These parameters characterize the environment of the Chilean Altiplano from which this strain was isolated. Our experiments lead us to the conclusion that the formation of such structures is a likely adaptation to environmental extremes.

      Finding cyst-like cells in S. ureilytica Lr5/4 provides another example in the extending list of resistant cyst-like-producing, non-spore-forming bacteria. This list includes Pseudomonas aeruginosa (Zechman and Casida, 1982), Micrococcus luteus (Mulyukin et al., 1996), Thioalkalivibrio versutus (Loiko et al., 2003), Sinorhizobium meliloti (Loiko et al., 2011), and Arthrobacter globiformis (Demkina et al., 2000), among others (Suzina et al., 2004). Dormancy seems now to be extended to almost every known bacterial phylum, and it has been even suggested that a common “sporulating” ancestor to bacteria might have existed in the past (Tocheva et al., 2016). Indeed, the ability to enter into a dormant and resistant cellular state would have provided a clear advantage to bacteria in an early Earth exposed to frequent meteoritic bombardments and alternating periods of boiling and freezing temperatures (Nicholson et al., 2002; Müller et al., 2014). Although beneficial under fluctuating conditions, untimely dormancy has severe consequences and would be under negative selection under stable conditions (Maughan et al., 2009). Therefore, selective pressure from the environment seems to be an essential condition for this mechanism to be maintained and deployed. In the case of the Chilean Altiplano, poly-extreme conditions could have led S. ureilytica Lr5/4 to evolve or to express this trait. The same could be the case for Archaea, in which the formation of spherical forms in response to a decrease in water activity (Fendrihan et al., 2012), appears to correlate with the survival and re-culturing of halophilic Archaea from salt deposits as old as the Permian (Stan-Lotter et al., 2002; Gruber et al., 2004). In environments favoring dormancy (such as polyextremophilic environments), this survival strategy could thus be found among bacterial and archaeal groups previously considered non-dormant. In contrast, in environments in which dormancy incurs a fitness cost, this trait could have been lost explaining the diversity of non-dormant groups reported nowadays.

      To date, harmful conditions [extreme or autolyzing (Muliukin et al., 1997)] seem to trigger the formation of resting cells in various microorganisms. However, whether the molecular mechanism for the formation of resting structures is similar among bacterial taxa remains to be shown. It has been hypothesized that cellular differentiation in these non-spore-forming organisms follow an alternative molecular pathway to the pathways described in the known spore-forming phyla (Suzina et al., 2004). The genetic analysis of the present study shed light on this unknown molecular mechanism. The results of the comparative analysis of the genome of S. ureilytica Lr5/4 with a database of genes involved in spore formation are striking. With the exception of Cyanobacteria, similarities were found. The homologous sporulation genes found in S. ureilytica Lr5/4 were either transcriptional factors, such as sigG, sigK, or bldH, or are homologous to house-keeping genes related to cellular division, such as spoIIIE, mreB, exoD, or exoF. These observations may lead to the suggestion that specialized transcriptional factors may be used for the alternative transcription of house-keeping genes for cellular differentiation, a rather economic alternative to the maintenance of a costly spore-formation machinery. The existence of spoVR homologs in different organisms appears as an important commonality to spore production. The SpoVR protein is responsible for spore cortex assembly in Firmicutes, but it is also found in myxospore-formation. The spore peptidoglycan in Firmicutes and Myxococcales differs from the one in vegetative cells (Atrih and Foster, 1999; Bui et al., 2009). Therefore, the presence of this gene in a resting cyst producer, like S. ureilytica Lr5/4 also highlights the need for the presence of specialized genes in the molecular pathway of cyst-like production.

      An extensive genetic analysis can only provide candidates to the molecular pathways. Whether these genes are indeed expressed while the bacterial cells are stressed, should be investigated for all cyst-producing organisms using gene expression analysis. A phylogenomic analysis of those genes, however, answers two other fundamental questions: what is the origin of dormancy in strain Lr5/4? And, do these genes correspond to a “dormancy gene core” widespread among unrelated taxa? Two scenarios have been proposed for the presence of dormancy in Serratia. The first scenario proposes dormancy is the result of gene transfer across unrelated bacterial genera (Ajithkumar et al., 2003), a phenomenon that has been reported as commonly occurring among hyperthermophiles in geothermal sites (van Wolferen et al., 2013). This scenario could be possible because genes related to sporulation have been found in genomes of non-sporulating species (Onyenwoke et al., 2004). However, in our opinion, gene transfer is unlikely because the formation of a fully functional spore is an exquisitely complex process that in the case of Firmicutes (the most studied model in terms of genetics and the mechanism of sporulation), requires a minimum of 60 genes, located throughout the bacterial chromosome and acting co-ordinately at different stages (Galperin et al., 2012). The second scenario is that dormancy is an ancient train within this genus, expressed only under particular conditions. The phylogenetic analysis of genes related to dormancy that are found in the genome of S. ureilytica Lr5/4 shows that the phylogenetic neighborhood of these genes is clearly proteobacterial, providing proof of the ancestry of these genes within the clade of Proteobacteria. Our phylogenomic analysis also revealed that these genes are not uniquely found in the genome of S. ureilytica Lr5/4, nor only acquitted among Proteobacteria. On the contrary, these genes are widespread among various sporulating taxa. These genes are required for cellular division, or control at crucial checkpoints. A previous analysis on Firmicutes sporulation-related genes showed that some of these genes may also be found in non-spore formers (Galperin, 2013). These genes are also shown to be involved in cell division and signaling/sensing and can be remnants of a sporulating ancestor (Onyenwoke et al., 2004). With the exception of cotU, which was solely found in Proteobacteria, the other genes detected in the genome of S. ureilytica Lr5/4 are found in two or more sporulating phyla. This observation suggests that sporulating organisms might share common dormancy-related genes. However, although homologs, whether these genes play a role in sporulation or have a function related to dormancy in all phyla needs to be verified in future studies.

      Conclusion

      The widespread phylogenomic distribution of dormancy genes should be seen under the prism of widespread cellular differentiation in bacteria, promoted under harsh conditions. The output of the phylogenomic and functional analysis of dormancy genes highlights that cellular differentiation to a more resistant structure, for instance, is highly related to cellular division. Moreover, it also shows that the choice between division and differentiation depends on cellular sensing. Since many non-spore-formers produce resting structures (Muliukin et al., 1997; Suzina et al., 2004, 2006), the molecular mechanism for such differentiation capability may not involve a complex specialized pathway. As our findings indicate, the production of resistant resting cells from previously known non-spore-forming bacteria, may involve specific cell-division genes. An important step into the understanding of dormancy in S. ureilytica Lr5/4 is to discover all the genes involved in cell differentiation and to reconstruct the molecular pathway, but also study the levels of expression for each gene. Future work should involve an expression analysis under the two morphotypes (vegetative and resting) for all cyst-producing bacteria. However, with this work we support the scenario that dimorphic bacterial cells are a widespread phenomenon, showing that one of the two cell morphologies, the most resistant, is produced after exposure to harsh environmental conditions. The ecological significance of this mechanism is high, since survival seems to be feasible to an increasing number of extreme conditions, but also it seems to be extended to more and more bacterial phyla.

      Author Contributions

      SF performed physiological and bioinformatics analysis and identification and wrote the manuscript. TJ performed bioinformatics analysis. TW participated to the sampling and performed isolation. WK, IP, and AA-D helped for the collection of physiological data. VM participated to the sampling. RL performed the MALDI-TOF analysis. JS performed the lipid analysis. SJ performed the automated genome annotation. PC supervised the automatic genome annotation. CD coordinated access to the sampling site and participated to the in situ collection of data and sampling. PJ designed the study, participated to the data analysis, and wrote the manuscript. All authors contributed and approved the final version of the manuscript.

      Conflict of Interest Statement

      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. We acknowledge funding from the Swiss National Science Foundation projects 31003A_152972, P2NEP3_178561 and from Fondation Pierre Mercier pour la science.

      For the sampling campaign, we would like to thank Dr. Ludovic Russel-Delif and Mrs. Nicole Jeanneret. For help in data analysis, we would like to thank Dr. Nathalie Marinakis.

      Supplementary Material

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

      References Abel-Santos E. (2012). Bacterial Spores: Current Research and Applications. Poole: Horizon Scientific Press. Ajithkumar B. Ajithkumar V. P. Iriye R. Doi Y. Sakai T. (2003). Spore-forming Serratia marcescens subsp. sakuensis subsp. nov., isolated from a domestic wastewater treatment tank. Int. J. Syst. Evol. Microbiol. 53 253258. 10.1099/ijs.0.02158-0 12656181 Altschul S. F. Madden T. L. Schaffer A. A. Zhang J. H. Zhang Z. Miller W. (1997). Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25 33893402. 10.1093/nar/25.17.3389 9254694 Atrih A. Foster S. J. (1999). The role of peptidoglycan structure and structural dynamics during endospore dormancy and germination. Antonie Van Leeuwenhoek 75 299307. 10.1023/A:1001800507443 10510717 Barton L. (2005). Structural and Functional Relationships in Prokaryotes. New York, NY: Springer-Verlag. Blom J. Albaum S. P. Doppmeier D. Pühler A. Vorhölter F.-J. Zakrzewski M. (2009). EDGAR: a software framework for the comparative analysis of prokaryotic genomes. BMC Bioinformatics 10:154. 10.1186/1471-2105-10-154 19457249 Brandes Ammann A. Kölle L. Brandl H. (2011). Detection of bacterial endospores in soil by terbium fluorescence. Int. J. Microbiol. 2011:435281. 10.1155/2011/435281 21754939 Brewer J. H. (1940). Clear liquid mediums for the “aerobic” cultivation of anaerobes. J. Am. Med. Assoc. 115 598600. 10.1001/jama.1940.72810340001009 Bui N. K. Gray J. Schwarz H. Schumann P. Blanot D. Vollmer W. (2009). The peptidoglycan sacculus of Myxococcus xanthus has unusual structural features and is degraded during glycerol-induced Myxospore development. J. Bacteriol. 191 494505. 10.1128/JB.0060808 18996994 Carver T. Thomson N. Bleasby A. Berriman M. Parkhill J. (2009). DNAPlotter: circular and linear interactive genome visualization. Bioinformatics 25 119120. 10.1093/bioinformatics/btn578 18990721 Castresana J. (2000). Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol. Biol. Evol. 17 540552. 10.1093/oxfordjournals.molbev.a026334 10742046 Cockell C. S. (2015). Astrobiology: Understanding Life in the Universe. Hoboken, NJ: John Wiley & Sons. Cordero R. R. Damiani A. Jorquera J. Sepúlveda E. Caballero M. Fernandez S. (2018). Ultraviolet radiation in the Atacama desert. Antonie Van Leeuwenhoek 111 13011313. 10.1007/s10482-018-1075-z 29605897 Cordero R. R. Damiani A. Seckmeyer G. Jorquera J. Caballero M. Rowe P. (2016). The solar spectrum in the atacama desert. Sci. Rep. 6:22457. 10.1038/srep22457 26932150 Demkina E. V. Soina V. S. Registan G. I. E. (2000). Formation of resting forms of Arthrobacter globiformis in autolyzing cell suspensions. Microbiology 69 314318. 10.1007/BF02756740 Dereeper A. Guignon V. Blanc G. Audic S. Buffet S. Chevenet F. (2008). Phylogeny.fr: robust phylogenetic analysis for the non-specialist. Nucleic Acids Res. 36 W465W469. 10.1093/nar/gkn180 18424797 Dong Q. Bauer C. E. (2015). Transcriptome analysis of cyst formation in Rhodospirillum centenum reveals large global changes in expression during cyst development. BMC Genomics 16:68. 10.1186/s12864-015-1250-9 25758168 Eddy S. R. (2011). Accelerated profile HMM searches. PLoS Comput. Biol. 7:e1002195. 10.1371/journal.pcbi.1002195 22039361 Edgar R. C. (2004). MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32 17921797. 10.1093/nar/gkh340 15034147 Edwards C. (1990). Microbiology of Extreme Environments. New York, NY: Mcgraw-Hill. Engler A. J. Sen S. Sweeney H. L. Discher D. E. (2006). Matrix elasticity directs stem cell lineage specification. Cell 126 677689. 10.1016/j.cell.2006.06.044 16923388 Felsenstein J. (1993). PHYLIP (Phylogeny Inference Package), Version 3.5 c. Philadelphia, PA: Joseph Felsenstein. Fendrihan S. Dornmayr-Pfaffenhuemer M. Gerbl F. Holzinger A. Grösbacher M. Briza P. (2012). Spherical particles of halophilic archaea correlate with exposure to low water activity–implications for microbial survival in fluid inclusions of ancient halite. Geobiology 10 424433. 10.1111/j.1472-4669.2012.00337.x 22804926 Finn R. D. Bateman A. Clements J. Coggill P. Eberhardt R. Y. Eddy S. R. (2014). Pfam: the protein families database. Nucleic Acids Res. 42 D222D230. 10.1093/nar/gkt1223 24288371 Forney L. Daae F. Torsvik V. (1997). Distribution of bacterioplankton in meromictic Lake Saelenvannet, as determined by denaturing gradient gel electrophoresis of PCR-amplified gene fragments coding for 16S rRNA. Appl. Environ. Microbiol. 63 33673373. 9292986 Galperin M. Y. (2013). Genome diversity of spore-forming firmicutes. Microbiol Spectr. 1 TBS-0015-2012. 10.1128/microbiolspectrum.TBS-0015-2012 Galperin M. Y. Mekhedov S. L. Puigbo P. Smirnov S. Wolf Y. I. Rigden D. J. (2012). Genomic determinants of sporulation in Bacilli and Clostridia: towards the minimal set of sporulation-specific genes. Environ. Microbiol. 14 28702890. 10.1111/j.1462-2920.2012.02841.x 22882546 Gerhardt P. (1994). Methods for General and Molecular Bacteriology. Washington, DC: American Society for Microbiology. Ghosh J. Larsson P. Singh B. Pettersson B. M. F. Islam N. M. Sarkar S. N. (2009). Sporulation in mycobacteria. Proc. Natl. Acad. Sci. U.S.A. 106 1078110786. 10.1073/pnas.0904104106 19541637 Girija K. R. Sasikala C. Ramana C. V. Spröer C. Takaichi S. Thiel V. (2010). Rhodobacter johrii sp. nov., an endospore-producing cryptic species isolated from semi-arid tropical soils. Int. J. Syst. Evol. Microbiol. 60 20992107. 10.1099/ijs.0.0117180 19854875 Grimont F. Grimont P. A. D. (2015). “Serratia,” in Bergey’s Manual of Systematics of Archaea and Bacteria, eds DeVos P. Chun J. Dedysh S. Hedlund B. Kämpfer P. Raine F. (Hoboken, NJ: John Wiley & Sons, Ltd). Gruber C. Legat A. Pfaffenhuemer M. Radax C. Weidler G. Busse H.-J. (2004). Halobacterium noricense sp. nov., an archaeal isolate from a bore core of an alpine Permian salt deposit, classification of Halobacterium sp. NRC-1 as a strain of H. salinarum and emended description of H. salinarum. Extremophiles 8 431439. 10.1007/s00792-004-0403-6 15290323 Guckert J. B. Antworth C. P. Nichols P. D. White D. C. (1985). Phospholipid, ester-linked fatty acid profiles as reproducible assays for changes in prokaryotic community structure of estuarine sediments. FEMS Microbiol. Lett. 31 147158. 10.1111/j.1574-6968.1985.tb01143.x Hemmerich C. Buechlein A. Podicheti R. Revanna K. V. Dong Q. (2010). An ergatis-based prokaryotic genome annotation web server. Bioinformatics 26 11221124. 10.1093/bioinformatics/btq090 20194626 Hoeniger J. F. M. Headley C. L. (1968). Cytology of spore germination in Clostridium pectinovorum. J. Bacteriol. 96 18351847. 4882027 Ichihara K. Fukubayashi Y. (2010). Preparation of fatty acid methyl esters for gas-liquid chromatography. J. Lipid Res. 51 635640. 10.1194/jlr.D001065 19759389 Junier T. Zdobnov E. M. (2010). The Newick utilities: high-throughput phylogenetic tree processing in the UNIX shell. Bioinformatics 26 16691670. 10.1093/bioinformatics/btq243 20472542 Katoh K. Misawa K. Kuma K. Miyata T. (2002). MAFFT: a novel method for rapid multiple sequence alignment based on fast fourier transform. Nucleic Acids Res. 30 30593066. 10.1093/nar/gkf436 12136088 Le S. Q. Gascuel O. (2008). An improved general amino acid replacement matrix. Mol. Biol. Evol. 25 13071320. 10.1093/molbev/msn067 18367465 Lechevalier H. A. Lechevalier M. P. Holbert P. E. (1966). Electron microscopic observation of the sporangial structure of strains of actinoplanaceae. J. Bacteriol. 92 12281235. 5926747 Loiko N. G. Kryazhevskikh N. A. Suzina N. E. Demkina E. V. Muratova A. Y. Turkovskaya O. V. (2011). Resting forms of Sinorhizobium meliloti. Microbiology 80:472. 10.1134/S0026261711040126 Loiko N. G. Soina V. S. Sorokin D. Y. Mityushina L. L. El’-Registan G. I. (2003). Production of resting forms by the gram-negative Chemolithoautotrophic bacteria Thioalkalivibrio versutus and Thioalkalimicrobium aerophilum. Microbiology 72 285294. 10.1023/A:1024291730779 12901006 Maughan H. Birky C. W. Nicholson W. L. (2009). Transcriptome divergence and the loss of plasticity in Bacillus subtilis after 6,000 generations of evolution under relaxed selection for sporulation. J. Bacteriol. 191 428433. 10.1128/JB.01234-08 18952793 Muliukin A. L. Lusta K. A. Griaznova M. N. Babusenko E. S. Kozlova A. N. Duzha M. V. (1997). [Formation of resting cells in microbial suspensions undergoing autolysis]. Mikrobiologiia 66 4249. 9132531 Müller A. L. de Rezende J. R. Hubert C. R. J. Kjeldsen K. U. Lagkouvardos I. Berry D. (2014). Endospores of thermophilic bacteria as tracers of microbial dispersal by ocean currents. ISME J. 8 11531165. 10.1038/ismej.2013.225 24351936 Müller S. Strack S. N. Ryan S. E. Kearns D. B. Kirby J. R. (2015). Predation by Myxococcus xanthus induces bacillus subtilis to form spore-filled megastructures. Appl. Environ. Microbiol. 81 203210. 10.1128/AEM.0244814 25326308 Mulyukin A. Lusta K. Gryaznova M. Kozlova A. Duzha M. Duda V. (1996). Formation of resting cells by Bacillus cereus and Micrococcus luteus. Microbiology 65 683689. Muyzer G. Teske A. Wirsen C. Jannasch H. (1995). Phylogenetic-relationships of Thiomicrospira species and their identification in deep-sea hydrothermal vent samples by denaturing gradient gel-electrophoresis of 16S rDNA fragments. Arch. Microbiol. 164 165172. 10.1007/BF02529967 7545384 Nicholson W. L. Fajardo-Cavazos P. Rebeil R. Slieman T. A. Riesenman P. J. Law J. F. (2002). Bacterial endospores and their significance in stress resistance. Antonie Van Leeuwenhoek 81 2732. 10.1023/A:1020561122764 Onyenwoke R. U. Brill J. A. Farahi K. Wiegel J. (2004). Sporulation genes in members of the low G+C Gram-type-positive phylogenetic branch (Firmicutes). Arch. Microbiol. 182 182192. 10.1007/s00203-004-0696-y 15340788 Sadasivan L. Neyra C. A. (1985). Flocculation in Azospirillum brasilense and Azospirillum lipoferum: exopolysaccharides and cyst formation. J. Bacteriol. 163 716723. 3894333 Sadoff H. L. (1973). Comparative aspects of morphogenesis in three prokaryotic genera. Annu. Rev. Microbiol. 27 133153. 10.1146/annurev.mi.27.100173.001025 4201689 Schubotz F. Meyer-Dombard D. R. Bradley A. S. Fredricks H. F. Hinrichs K.-U. Shock E. L. (2013). Spatial and temporal variability of biomarkers and microbial diversity reveal metabolic and community flexibility in streamer biofilm communities in the lower geyser basin. Yellowstone National Park. Geobiology 11 549569. 10.1111/gbi.12051 23981055 Setlow B. Atluri S. Kitchel R. Koziol-Dube K. Setlow P. (2006). Role of dipicolinic acid in resistance and stability of spores of Bacillus subtilis with or without DNA-protective alpha/beta-type small acid-soluble proteins. J. Bacteriol. 188 37403747. 10.1128/JB.00212-06 16707666 Sghaier H. Narumi I. Satoh K. Ohba H. Mitomo H. (2007). Problems with the current deinococcal hypothesis: an alternative theory. Theory Biosci. 126 4345. 10.1007/s12064-007-0004-x 18087756 Sievers F. Wilm A. Dineen D. Gibson T. J. Karplus K. Li W. (2011). Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol. Syst. Biol. 7:539. 10.1038/msb.2011.75 21988835 Sorokin D. Y. Berben T. Melton E. D. Overmars L. Vavourakis C. D. Muyzer G. (2014). Microbial diversity and biogeochemical cycling in soda lakes. Extremophiles 18 791809. 10.1007/s00792-014-0670-9 25156418 Stan-Lotter H. Pfaffenhuemer M. Legat A. Busse H.-J. Radax C. Gruber C. (2002). Halococcus dombrowskii sp. nov., an archaeal isolate from a Permian alpine salt deposit. Int. J. Syst. Evol. Microbiol. 52 18071814. 12361290 Stevenson K. E. Vaughn R. H. Crisan E. V. (1972). Fixation of mature spores of Clostridium botulinum. J. Bacteriol. 109 12951297. 4551752 Suzina N. E. Mulyukin A. L. Dmitriev V. V. Nikolaev Y. A. Shorokhova A. P. Bobkova Y. S. (2006). The structural bases of long-term anabiosis in non-spore-forming bacteria. Adv. Space Res. 38 12091219. 10.1016/j.asr.2005.09.020 Suzina N. E. Mulyukin A. L. Kozlova A. N. Shorokhova A. P. Dmitriev V. V. Barinova E. S. (2004). Ultrastructure of resting cells of some non-spore-forming bacteria. Microbiology 73 435447. 10.1023/B:MICI.0000036990.94039.af Tocheva E. I. Ortega D. R. Jensen G. J. (2016). Sporulation, bacterial cell envelopes and the origin of life. Nat. Rev. Microbiol. 14 535542. 10.1038/nrmicro.2016.85 28232669 Traag B. A. Driks A. Stragier P. Bitter W. Broussard G. Hatfull G. (2010). Do mycobacteria produce endospores? Proc. Natl. Acad. Sci. U.S.A. 107 878881. 10.1073/pnas.0911299107 20080769 van Wolferen M. Ajon M. Driessen A. J. M. Albers S.-V. (2013). How hyperthermophiles adapt to change their lives: DNA exchange in extreme conditions. Extremophiles 17 545563. 10.1007/s00792-013-0552-6 23712907 Way J. C. (1996). The mechanism of bacterial asymmetric cell division. Bioessays 18 99101. 10.1002/bies.950180205 8851042 Zechman J. M. Casida L. E. (1982). Death of Pseudomonas aeruginosa in soil. Can. J. Microbiol. 28 788794. 10.1139/m82-120 6816464
      ‘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 0016librbn.com.cn
      judemt.com.cn
      www.eleonline.com.cn
      kokwz.net.cn
      hilegua.com.cn
      fixiapac.com.cn
      www.ruuyue.com.cn
      ssyukd.com.cn
      www.vrfenzi.com.cn
      ufcekj.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