Survey
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
FEMS Microbiology Ecology, 92, 2016, fiw012 doi: 10.1093/femsec/fiw012 Advance Access Publication Date: 28 January 2016 Research Article RESEARCH ARTICLE Biogeography of thermophilic phototrophic bacteria belonging to Roseiflexus genus Vasil A. Gaisin1,∗ , Denis S. Grouzdev1 , Zorigto B. Namsaraev2,3 , Marina V. Sukhacheva1 , Vladimir M. Gorlenko2 and Boris B. Kuznetsov1 1 Institute of Bioengineering, Research Center of Biotechnology of the Russian Academy of Sciences, Leninsky Ave., 33, bld. 2, Moscow 119071, Russian Federation, 2 Winogradsky Institute of Microbiology, Research Center of Biotechnology of the Russian Academy of Sciences, Leninsky Ave., 33, bld. 2, Moscow 119071, Russian Federation and 3 Department of Biotechnology and Bioenergy, National Research Centre ‘Kurchatov Institute’, Akademika Kurchatova pl., 1, Moscow 123182, Russian Federation ∗ Corresponding author: Institute of Bioengineering, Research Center of Biotechnology of the Russian Academy of Sciences. Leninsky Ave., 33, bld. 2, Moscow 119071, Russian Federation. Tel: +7(499)135-12-40; E-mail: [email protected] One sentence summary: The analysis of 16S rRNA gene sequences and environmental data revealed the existence of the specific phylogeographic pattern for thermophilic phototrophic bacteria belonging to Roseiflexus genus. Editor: Riks Laanbroek ABSTRACT Isolated environments such as hot springs are particularly interesting for studying the microbial biogeography. These environments create an ‘island effect’ leading to genetic divergence. We studied the phylogeographic pattern of thermophilic anoxygenic phototrophic bacteria, belonging to the Roseiflexus genus. The main characteristic of the observed pattern was geographic and geochronologic fidelity to the hot springs within Circum-Pacific and Alpine-HimalayanIndonesian orogenic belts. Mantel test revealed a correlation between genetic divergence and geographic distance among the phylotypes. Cluster analysis revealed a regional differentiation of the global phylogenetic pattern. The phylogeographic pattern is in correlation with geochronologic events during the break up of Pangaea that led to the modern configuration of continents. To our knowledge this is the first geochronological scenario of intercontinental prokaryotic taxon divergence. The existence of the modern phylogeographic pattern contradicts with the existence of the ancient evolutionary history of the Roseiflexus group proposed on the basis of its deep-branching phylogenetic position. These facts indicate that evolutionary rates in Roseiflexus varied over a wide range. Keywords: microbial biogeography; phylogeography; geotype; bacterial speciation; evolutionary rates; Roseiflexus INTRODUCTION The studies of the last 15 years show that microbial distribution is determined not only by environmental factors but also geographical/historical ones (Martiny et al. 2006; Hanson et al. 2012; Lindström and Langenheder 2012). This modern view was summarized by van der Gast: ‘Some things are everywhere and some things are not. Sometimes the environment selects and sometimes it doesn’t’ (van der Gast 2015). Here, the author transformed the Baas-Becking hypothesis: ‘Everything is everywhere, but, the environment selects’ (Baas-Becking 1934) which postulated the exclusive ecological approach in study of microbial diversity and distribution. Data that clearly prove the significance of historical factors were obtained while studying the population structure of Received: 26 October 2015; Accepted: 27 January 2016 C FEMS 2016. All rights reserved. For permissions, please e-mail: [email protected] 1 2 FEMS Microbiology Ecology, 2016, Vol. 92, No. 3 hyperthermophilic archaea belonging to the genera Sulfolobus (Whitaker, Grogan and Taylor 2003) and Pyrococcus (EscobarPáramo, Ghosh and DiRuggiero 2005), as well as after a study of biogeography of thermophilic cyanobacteria belonging to the genera Synechoccocus (Papke et al. 2003) and Mastigocladus (Miller, Castenholz and Pedersen 2007). In all the mentioned cases, specific conditions of hydrothermal springs were responsible for the phenomenon called ‘island effect’ (Hreggvidsson et al. 2012). However, such effect was not revealed for mesophilic cyanobacteria (van Gremberghe et al. 2011); in this case, the higher resolution power of genetic markers is required. Apparently, stenobiont microorganisms, including thermophiles, are affected by the process which is similar to vicariance for animals. The organisms with more restricted ecologic abilities are known to be more endemic (Barberán et al. 2014). Thus, the modern phylogenetic patterns of thermophilic prokaryotes should reflect the historical factors influencing the distribution of genetic lineages. The objects of the current study were bacteriochlorophyll-a containing filamentous thermophilic bacteria, belonging to the Chloroflexi class. These bacteria inhabit thermophilic cyanobacterial mats and in some cases form very thick red-colored layers (Boomer et al. 2000, 2002). Isolated pure cultures of the representatives of this group were described as filamentous anoxygenic phototrophic bacteria (FAPB) of the Roseiflexus genus. At the moment, three pure cultures are available in microbial collections: Roseiflexus castenholzii isolated from Nakabusa hot spring, Japan (Hanada et al. 2002), Roseiflexus sp. RS-1 and RS-2 isolated from Octopus hot spring, Yellowstone, North America (van der Meer et al. 2010). These bacteria are obligate thermophiles growing at 45–60◦ C (Topt. = 50◦ C for R. castenholzii and 55–60◦ C for Roseiflexus sp. RS-1 and RS-2). The presence of Roseiflexus-like phylotypes has been revealed in microbial communities inhabiting hydrotherms of Japan (Hanada et al. 2002; Everroad et al. 2012), Yellowstone (Boomer et al. 2002, 2009; Nübel et al. 2002; Spear et al. 2005), Kamchatka Peninsula (Burgess et al. 2012; Akimov et al. 2013), Thailand (Portillo et al. 2009), Bulgaria (Tomova et al. 2010), Tibet (Lau and Pointing 2009; Lau, Aitchison and Pointing 2009) and Andes (Engel, Johnson and Porter 2013). Recently, a new Roseiflexuslike phylotype with a 3% variance (difference) from the known 16S rRNA gene, was found in the cyanobacterial mat of Alla hot spring (Buryat Republic, Russia) (Gaisin et al. 2015). According to the recently obtained data, the same phylotype is also present in Tsenher hydrotherm (Mongolia). The presence of FAPB belonging to R. castenholzii phylotype in Thermophilny hot spring (Kamchatka) has also been revealed. As a result of the statistical and cluster analysis of the Roseiflexus-like 16S rRNA sequences, along with data on habitats where the phylotypes were observed, we propose the existence of a specific phylogeographic pattern for bacteria belonging to the Roseiflexus genus. MATERIALS AND METHODS Original sequences from Tsenher and Thermophilny hot springs The Roseiflexus-like sequences from Tsenher (Mongolia, Global Positioning System (GPS) coordinates: 47.316667, 101.653611) and Thermophilny (Kamchatka, GPS coordinates: 54.816667, 160.016667) hot springs were presented for the first time in this study. Sequences from the Tsenher hot spring were obtained by analyzing the cloned 16S rRNA gene polymerase chain reaction (PCR) fragments, amplified from the total DNA of the cyanobacterial mat with primers Univ27F/Univ1492R (Lane 1991). Sequence from the Thermophilny hot spring was obtained using Roseiflexus-specific 16S rRNA primers Rof97f/Rof940r. The detailed procedures of DNA extraction and PCR protocols used were described previously (Gaisin et al. 2015). The temperature and pH of the sampling site were 62.2◦ C and 8.4, respectively for the Tsenher hot spring, and 54.0◦ C and 5.9 for the Thermophilny hot spring. The clone sequences from the Tsenher hot spring were deposited in GenBank under accession numbers KT258753– KT258772. The sequence from the Thermophilny hot spring derived by the Roseiflexus-specific primers was deposited in GenBank under accession numbers KT266780. Sequence analysis Roseiflexus-like sequences were retrieved from GenBank using Basic Local Alignment Search Tool on-line service (http://www. ncbi.nlm.nih.gov/blast). Phylogenetic reconstructions were performed using MEGA 5.0 (Research Center of Biotechnology of the Russian Academy of Sciences, Moscow, Russian Federation) (Tamura et al. 2011). Non-metric multidimensional scaling (nMDS) was performed in R using the statistical packages vegan (Oksanen et al. 2007) and ecodist (Goslee and Urban 2007). The three dimensional (3D) plots of nMDS analysis were created with R package rgl (Adler and Murdoch 2012). Statistical analyses Mantel test was used to estimate the correlations between matrices of temperature, water pH, geographic positions of isolation points and genetic distances with a significance testing against 1000 permutations. Mantel test was performed with R using the statistical package ecodist (Goslee and Urban 2007). Vincenty’s formula and GPS coordinates were used to calculate the geographic distances between isolation points. RESULTS Sets of Roseiflexus-like sequences of 16S rRNA gene Sequences with a similarity level higher than 95% to R. castenholzii (CP000804) and Roseiflexus RS-1 (CP000686) were retrieved from GenBank. According to (Tindall et al. 2010), the selected sequences belonged to the representatives of the same genus. The sequences found in GenBank together with the sequences obtained during this study were assembled into the primary dataset of Roseiflexus-like sequences (Table S1, Supporting Information). The majority of the sequences originated from the Yellowstone hot springs (Fig. S1, Supporting Information). For further analysis, the region between the 400th and 1000th nucleotide of the gene was chosen (Fig. S2, Supporting Information) due to its presence in a maximal number of sequences from all habitats where the sequences were found. Thereafter, the following sequences were excluded from the primary set: (i) short sequences not overlapping the above gene region; (ii) sequences from non-thermophilic and unknown locations; (iii) sequences with nucleotide anomalies and unambiguous letters. The final set included 184 sequences (Table S2, Supporting Information), which originated from 26 thermal springs, 11 of which were located in Yellowstone. Also, the final set represented 12 areas of hydrothermal activity. Gaisin et al. 3 sequences. Considering the comparison of the geographic map and the phylogenetic tree, the phylogeography pattern was observed on the global scale (Fig. 2). The result of phylogenetic reconstruction based on phylogenetic dichotomy clustering is similar to the nMDS scheme. The main feature of the resulted phylogenetic tree is that every separate cluster on it included sequences retrieved from a single restricted geographic region. The exception is the fact that sequences obtained from Yellowstone hot springs belonged to two separate distant clusters. One of these clusters included sequences found in hydrotherms of North and South Americas, whereas the second one consisted of sequences found in Eurasia. All Roseiflexus-like sequences were detected in hot springs from regions whose hydrothermal activity is determined by the geological dynamics in the Circum-Pacific and AlpineHimalayan-Indonesian orogenic belts (Collins et al. 2011). Effect of geographic distance on phylogenetic divergence Figure 1. nMDS-derived clusters of Roseiflexus-like sequences. The color of the spheres corresponds to the color of the regional names. Cluster analysis of Roseiflexus-like sequences nMDS was used to analyze the sequence clustering. The results were performed in 3D graphic reconstruction (Fig. S3, Supporting Information). The most informative projection of 3D reconstruction is as shown in Fig. 1. According to the obtained data, Roseiflexus-like sequences formed groups corresponding to the geographic regions from where these sequences were obtained. The sequences from various sources formed separate non-overlapping groups while sequences from Yellowstone springs formed two separate groups. Maximum Likelihood algorithm was used to construct the dendrogram of phylogenetic similarity of 16S rRNA gene Mantel test was used to analyze correlations between genetic and geographic distances, as well as water temperature and pH. The results obtained supported non-random correlation between parameters of genetic and geographic distances with a high probability (Fig. 3), R2 = 0.3446, P < 0.0001. No effect was found for pH and water temperature on clustering; R2 = 0.0026, P = 0.3098 for pH and R2 = 0.0012, P = 0.4926 for water temperature correspondingly (Fig. 3). DISCUSSION The results obtained during our study are in agreement with the ‘Geotype-plus-Boeing’ model proposed by Cohan and Perry (Cohan and Perry 2007). According to the model, the whole group of Roseiflexus bacteria is considered to be a uniform ecotype, that is ‘group of bacteria that are ecologically similar to one another, so Figure 2. Global distribution of Roseiflexus-like sequences. The tree was generated using maximum likelihood method. Bootstrap values were calculated through the comparison of 1000 alternative trees. DTAM-Digital Tectonic Activity Map (Lowman et al. 1999). 4 FEMS Microbiology Ecology, 2016, Vol. 92, No. 3 Figure 3. The results of Mantel test. Relationships between genetic distances and geographic distances (panel a), Relationships between genetic distances and pH (panel b), Relationships between genetic distances and temperatures (panel c). The values on the axis of temperature, pH and genetic distances are given as Euclid distances. similar that genetic diversity within the ecotype is limited by a cohesive force’ (Cohan and Perry 2007). The main features determining the Roseiflexus-ecotype are obligate thermophily and ecological/physiological dependence from cyanobacteria in mats that provide a nutrient feeding (van der Meer et al. 2005; Klatt et al. 2013; Rodionova et al. 2015). The high level of sequence homology found within each phylogenetic cluster supports the hypothesis that natural Roseiflexus populations exist under strong influence of stabilizing selection and these bacteria are stenobionts (Boomer et al. 2002; Klatt et al. 2013; Gaisin et al. 2015). It therefore means that the Roseiflexus population within each hydrothermal region is relatively genetically uniform and does not undergo significant genetic divergence. However, populations from different geographical locations are phylogenetically different (Fig. 2). They form the so-called ‘geotypes’ (Cohan and Perry 2007). Altogether the geotypes form the global phylogeographic pattern of the Roseiflexus ecotype. One of the main features of the Roseiflexus phylogeographic pattern is that all known Roseiflexus-like sequences were retrieved from habitats situated within the areas of geological activity along subduction zones of the Circum-Pacific and Alpine-Himalayan-Indonesian orogenic belts (Collins et al. 2011) (Fig. 2). Microbial communities of Iceland, Middle Atlantic Ridge and African hydrotherms were studied and to our knowledge there are no Roseiflexus-like phylotypes (Skirnisdottir et al. 2000; van der Meer et al. 2008; Tobler and Benning 2011; Tekere et al. 2012; Sahm et al. 2013). Roeselers and co-authors reported about the 16S rRNA gene sequences from hot spring in Greenland that were related to Roseiflexus bacteria (Roeselers et al. 2007). However, additional analysis revealed that these sequences do not belong to the Roseiflexus ecotype, because identity between sequences from Greenland and the Roseiflexus bacteria was 93%. Anomalous geographical distribution was also observed for thermophilic cyanobacteria Synechococcus spp., which was not observed at all in hot springs in Iceland, Alaska and the Azores, even though it would eagerly grow in water from Icelandic hot springs as reported by R.W. Castenholz (Papke et al. 2003). The observed phylogeographic pattern of the Roseiflexuslike bacteria could be explained by the specific physiological requirements of these bacteria. Being obligate thermophiles (Hanada et al. 2002; van der Meer et al. 2010), the Roseiflexus representatives require hydrothermal activity associated mostly with continental rift zones and volcanic activity in plate subduction zones. Moreover, the Roseiflexus does not produce resting forms (spores, cysts, etc.) and thus can- not spread over large distances. Hence, we propose that the Roseiflexus bacteria could migrate from one geothermal area to another adjacent one along active subduction zones in the Circum-Pacific and Alpine-Himalayan-Indonesian orogenic belts. Such scenario is similar to distribution of deep-sea hydrothermal vent biota along the ocean ridges (Vrijenhoek 2010). However, in our case mechanism of transfer is not clear. Most likely, the modern phylogeographic pattern of Roseiflexus-like geotypes reflects ancient historical legacies, just like it has already been mentioned by some authors for extremophilic microorganisms (Takacs-Vesbach et al. 2008; Bahl et al. 2011). The proposal that Roseiflexus geotypes were formed under the influence of historical factors is supported by the region-associated distribution of the Roseiflexus phylotypes (Fig. 2). The regional correlations were also supported by the Mantel test results, whereas no correlations were found with the temperature and pH of spring water (Fig. 3). The revealed correlation of Roseiflexus-like bacteria habitats along with the history of the geological activity allows to propose a timeframe of evolutionary scenario. To our knowledge this is the first geochronological scenario of intercontinental prokaryotic taxon divergence. According to data on the Earth’s tectonic evolution and geochronology, the branching order of Roseiflexuslike phylotypes could be the following. About 200 million years (My) ago (Jurassic), the supercontinent of Pangea started to break up forming the Atlantic Ocean (Neall and Trewick 2008). Simultaneously, the subduction of ocean floor on the borders of Pangea started to form the Circum-Pacific orogenic belt and resulted in an intense volcanic activity (Scotese 1987). After the break up of Pangaea, the North American plate began drifting westward. Around 200 My ago, the Pacific coast turned from a passive margin into an active one with multiple signs of hydrothermal activity that is still active today (Sigloch and Mihalynuk 2013). The formation of South America occurred later, around 110 My ago, during the opening of the South Atlantic Ocean (Orme 2007). The scenario of continental plate formation is also modeled in the PALEOMAP Project (Scotese 2001) (see the beautiful illustration at http://www.scotese.com/Default.htm). Such order of events strongly correlates with the branching order of Roseiflexus-like phylotypes on the corresponding dendrogram; the South American Roseiflexus group is younger than that of North America (Fig. 2). The Eurasian cluster is divided into three subclusters: North Eastern (Kamchatka, Japan and Yellowstone), Central and South Gaisin et al. Asian (Buryatia, Mongolia, Tibet, Southern China and Thailand) and dispersed subcluster which include sequences from Bulgaria and Eastern China. The formation of the North Eastern subcluster could be linked to accretion of the intra-oceanic volcanic arcs with Eurasia and North America. These arcs existed in the central part of the Panthalassa Ocean and were formed in the course of intraoceanic subduction, ∼245–201 My ago (Middle-Upper Triassic) (Van der Meer et al. 2012). After a northward drift, these volcanic arcs accreted to the continental margins of Western North America, North-East Asia and Japan ∼145–100 My ago (EarlyMiddle Cretaceous) (Van der Meer et al. 2012). The Central and South Asian groups could have been formed during a series of events that followed the India-Eurasia collision. After the break up of Pangaea (200 My ago), the Indian plate separated from Africa and Antarctica and started moving northward. In ∼65 My ago, India collided with Eurasia (White and Lister 2012). This resulted in the formation of Alpine-HimalayanIndonesian orogenic belt including uplift of the Himalayas and Tibetan Plateau. The same mechanism may be responsible for the formation of the Baikal rift zone (Mats and Perepelova 2011). These events could lead to the increase of hydrothermal activity along faults formed in the south part of Eurasia (Bulgaria, Tibet, South China and Thailand) and in the Baikal rift zone (Buryatia and Mongolia). It is more difficult to explain the origin of the third group which includes sequences from Bulgaria and Eastern China. Taking into account their deep branching, its formation should have happened a long time ago, while geographical dispersal could have been caused by a recent transportation (the so-called Boeing-effect) or other factors (Cohan and Perry 2007). We cannot exclude a possibility of the fact that the neighbor clustering of Bulgarian and Chinese sequences is an error of the dichotomous Maximum Likelihood analysis. For example, this clustering is not obvious in nMDS analysis (Fig. S3, Supporting Information), but additional data are required. The results of this study contribute to a problem of the evolutionary rates in deep-branching bacteria that is highly disputed in previous years (Dvořák et al. 2015; Plata, Henry and Vitkup 2015; Schopf et al. 2015). Thermophilic FAPB are considered to be one of the most ancient living organisms on Earth (CavalierSmith 2010). According to the modern reconstructions, the average water temperature of Precambrian ocean was about 55–85◦ C (from about 3.5 to 2 Gyr ago) (Knauth 2005; Robert and Chaussidon 2006). The microbial mats were formed by anoxygenic phototrophic bacteria, developed in the photic zone of the Precambrian ocean (Tice and Lowe 2004). Perhaps, thermophilic FAPB were one of these ancient bacteria and had found a way to survive up to date in hydrothermal springs. Earlier results show that Roseiflexus bacteria are one of the most deep-branching phylotypes among modern FAPB (Grouzdev et al. 2015). Contrary to these results, the observed phylogeographic pattern of the Roseiflexus bacteria was probably formed during relatively recent geological events. The continents reached their modern configuration in the Miocene (23–5 My ago), after the break up of Pangaea that started 200 My ago. In this case, we can suggest that the tempo of evolution in the Roseiflexus genus varied over the time from slow long-term to rapid diversification. Also, we did not exclude other scenarios while explaining the observed pattern. An alternative explanation implies a serial convergence of different ancient phylotypes of Roseiflexus through horizontal gene transfer, followed by the divergence of a common phylotype into local populations, thereby erasing the ancient phylogenetic pattern. A similar scenario was sug- 5 gested recently for cyanobacteria (Dvořák et al. 2014). Another explanation is that the ancient phylogenetic diversity of Roseiflexus was devastated by environmental change, whereas the modern diversity of Roseiflexus originated from a survived phylotype. The later scenario was proposed for biogeographic distribution of the local populations of Sulfurihydrogenibium bacteria in the Yellowstone Caldera (Takacs-Vesbach et al. 2008). SUPPLEMENTARY DATA Supplementary data are available at FEMSEC online. ACKNOWLEDGEMENTS The authors are grateful to researchers from the Laboratory of Microbiology at the Institute of General and Experimental Biology, Siberian Branch, RAS, Ulan-Ude, for organizing the expedition to the thermal springs of the Barguzin Valley (Alla) and Tsenher hot spring. They are also grateful to Dr Kovaleva O. from Winogradsky Institute of Microbiology RAS, who kindly provided the sample from Thermophilny (Kamchatka). Special thanks to Kolganova T.V. and Baslerov R.V. from the Molecular Diagnostics Laboratory, Centre for Bioengineering RAS. The work was carried out using the scientific equipment of Core Research Facility of Centre Bioengineering RAS. FUNDING This work was supported by Russian Foundation for Basic Research, grant numbers 15-04-07655\15; Russian Federation President Grant NS 6150.2014.4 and Program of the Presidium of the Russian Academy of Sciences ‘The evolution of the organic world and planetary processes’. Conflict of interest. None declared. REFERENCES Adler D, Murdoch D. Version rgl: 3D visualization device system (OpenGL). R Package 092. 2012, https://cran.r-project.org (10 January 2016, date last accessed). Akimov V, Podosokorskaya O, Shlyapnikov M et al. Dominant phylotypes in the 16S rRNA gene clone libraries from bacterial mats of the uzon caldera (Kamchatka, Russia) hydrothermal springs. Microbiology 2013;82:721–7. Baas-Becking LGM. Geobiologie of inleiding tot de milieukunde. The Hague: WP Van Stockum & Zoon, 1934. Bahl J, Lau MC, Smith GJ et al. Ancient origins determine global biogeography of hot and cold desert cyanobacteria. Nat Commun 2011;2:163. Barberán A, Ramirez KS, Leff JW et al. Why are some microbes more ubiquitous than others? Predicting the habitat breadth of soil bacteria. Ecol Lett 2014;17:794–802. Boomer SM, Lodge DP, Dutton BE et al. Molecular characterization of novel red green nonsulfur bacteria from five distinct hot spring communities in Yellowstone National Park. Appl Environ Microb 2002;68:346–55. Boomer SM, Noll KL, Geesey GG et al. Formation of multilayered photosynthetic biofilms in an alkaline thermal spring in Yellowstone National Park, Wyoming. Appl Environ Microb 2009;75:2464–75. Boomer SM, Pierson BK, Austinhirst R et al. Characterization of novel bacteriochlorophyll-a-containing red filaments from 6 FEMS Microbiology Ecology, 2016, Vol. 92, No. 3 alkaline hot springs in Yellowstone National Park. Arch Microbiol 2000;174:152–61. Burgess EA, Unrine JM, Mills GL et al. Comparative geochemical and microbiological characterization of two thermal pools in the Uzon Caldera, Kamchatka, Russia. Microb Ecol 2012;63:471–89. Cavalier-Smith T. Deep phylogeny, ancestral groups and the four ages of life. Philos T R Soc B 2010;365:111–32. Cohan FM, Perry EB. A systematics for discovering the fundamental units of bacterial diversity. Curr Biol 2007;17: R373–86. Collins WJ, Belousova EA, Kemp AI et al. Two contrasting Phanerozoic orogenic systems revealed by hafnium isotope data. Nat Geosci 2011;4:333–7. Dvořák P, Casamatta DA, Hašler P et al. Putative extremely long evolutionary stasis in bacteria might be explained by serial convergence. P Natl Acad Sci USA 2015;112:E2559. Dvořák P, Casamatta DA, Poulı́čková A et al. Synechococcus: 3 billion years of global dominance. Mol Ecol 2014;23: 5538–51. Engel AS, Johnson LR, Porter ML. Arsenite oxidase gene diversity among Chloroflexi and Proteobacteria from El Tatio Geyser Field, Chile. FEMS Microbiol Ecol 2013;83:745–56. Escobar-Páramo P, Ghosh S, DiRuggiero J. Evidence for genetic drift in the diversification of a geographically isolated population of the hyperthermophilic archaeon Pyrococcus. Mol Biol Evol 2005;22:2297–303. Everroad RC, Otaki H, Matsuura K et al. Diversification of bacterial community composition along a temperature gradient at a thermal spring. Microbes Environ 2012;27:374–81. Gaisin VA, Kalashnikov AM, Sukhacheva MV et al. Filamentous anoxygenic phototrophic bacteria from cyanobacterial mats of Alla hot springs (Barguzin Valley, Russia). Extremophiles 2015;19:1067–76. Gast CJ. Microbial biogeography: the end of the ubiquitous dispersal hypothesis? Environ Microbiol 2015;17:544–6. Goslee SC, Urban DL. The ecodist package for dissimilarity-based analysis of ecological data. J Stat Softw 2007;22:1–19. Grouzdev DS, Kuznetsov BB, Keppen OI et al. Reconstruction of bacteriochlorophyll biosynthesis pathways in the filamentous anoxygenic phototrophic bacterium Oscillochloris trichoides DG-6 and evolution of anoxygenic phototrophs of the order Chloroflexales. Microbiology 2015;161: 120–30. Hanada S, Takaichi S, Matsuura K et al. Roseiflexus castenholzii gen. nov., sp. nov., a thermophilic, filamentous, photosynthetic bacterium that lacks chlorosomes. Int J Syst Evol Micr 2002;52:187–93. Hanson CA, Fuhrman JA, Horner-Devine MC et al. Beyond biogeographic patterns: processes shaping the microbial landscape. Nat Rev Microbiol 2012;10:497–506. Hreggvidsson GO, Petursdottir SK, Björnsdottir SH et al. Microbial speciation in the geothermal ecosystem In: Stan-Lotter H, Fendrihan S (eds). Adaption of Microbial Life to Environmental Extremes. Vienna: Springer, 2012, 37–67. Klatt CG, Liu Z, Ludwig M et al. Temporal metatranscriptomic patterning in phototrophic Chloroflexi inhabiting a microbial mat in a geothermal spring. ISME J 2013;7:1775–89. Knauth LP. Temperature and salinity history of the Precambrian ocean: implications for the course of microbial evolution. Palaeogeogr Palaeocl 2005;219:53–69. Lane D. 16S/23S rRNA sequencing. In: Stackebrandt E, Goodfellow M (eds). Nucleic Acid Techniques in Bacterial Systematics. New York: John Wiley and Sons. 1991, 125–75. Lau MC, Aitchison JC, Pointing SB. Bacterial community composition in thermophilic microbial mats from five hot springs in central Tibet. Extremophiles 2009;13:139–49. Lau MC, Pointing SB. Vertical partitioning and expression of primary metabolic genes in a thermophilic microbial mat. Extremophiles 2009;13:533–40. Lindström ES, Langenheder S. Local and regional factors influencing bacterial community assembly. Environ Microbiol Rep 2012;4:1–9. Lowman P, Yates J, Masuoka P et al. A Digital Tectonic Activity Map of the Earth. J Geosci Ed 1999;47:428–37. Martiny JBH, Bohannan BJ, Brown JH et al. Microbial biogeography: putting microorganisms on the map. Nat Rev Microbiol 2006;4:102–12. Mats VD, Perepelova TI. A new perspective on evolution of the Baikal Rift. Geosci Front 2011;2:349–65. Miller SR, Castenholz RW, Pedersen D. Phylogeography of the thermophilic cyanobacterium Mastigocladus laminosus. Appl Environ Microb 2007;73:4751–9. Neall VE, Trewick SA. The age and origin of the Pacific islands: a geological overview. Philos T R Soc B 2008;363: 3293–308. Nübel U, Bateson MM, Vandieken V et al. Microscopic examination of distribution and phenotypic properties of phylogenetically diverse Chloroflexaceae-related bacteria in hot spring microbial mats. Appl Environ Microb 2002;68:4593–603. Oksanen J, Kindt R, Legendre P et al. The vegan package. Community Ecology Package Version 1.15-1. 2007. http://cran. r-project.org (10 January 2016, date last accessed). Orme AR. The tectonic framework of South America. In: Veblen TT, Young KR, Orme AR (eds). The Physical Geography of South America. Los Angeles: Oxford University Press, 2007, 3–22. Papke RT, Ramsing NB, Bateson MM et al. Geographical isolation in hot spring cyanobacteria. Environ Microbiol 2003;5:650–9. Plata G, Henry CS, Vitkup D. Long-term phenotypic evolution of bacteria. Nature 2015;517:369–72. Portillo M, Sririn V, Kanoksilapatham W et al. Differential microbial communities in hot spring mats from Western Thailand. Extremophiles 2009;13:321–31. Robert F, Chaussidon M. A palaeotemperature curve for the Precambrian oceans based on silicon isotopes in cherts. Nature 2006;443:969–72. Rodionova IA, Li X, Plymale AE et al. Genomic distribution of Bvitamin auxotrophy and uptake transporters in environmental bacteria from the Chloroflexi phylum. Environ Microbiol Rep 2015;7:204–10. Roeselers G, Norris TB, Castenholz RW et al. Diversity of phototrophic bacteria in microbial mats from Arctic hot springs (Greenland). Environ Microbiol 2007;9:26–38. Sahm K, John P, Nacke H et al. High abundance of heterotrophic prokaryotes in hydrothermal springs of the Azores as revealed by a network of 16S rRNA gene-based methods. Extremophiles 2013;17:649–62. Schopf JW, Kudryavtsev AB, Walter MR et al. Sulfur-cycling fossil bacteria from the 1.8-Ga Duck Creek Formation provide promising evidence of evolution’s null hypothesis. P Natl Acad Sci USA 2015;112:2087–92. Scotese CR. Development of the circum-pacific panthallassic ocean during the early paleozoic. In: Monger JWH, Francheteau J (eds). Circum-Pacific Orogenic Belts and Evolution of the Pacific Ocean Basin. Washington, DC: American Geophysical Union, 1987:49–57. Scotese CR. Atlas of earth history, Volume 1, Paleogeography, PALEOMAP Project, Arlington, Texas, 2001. Gaisin et al. Sigloch K, Mihalynuk MG. Intra-oceanic subduction shaped the assembly of Cordilleran North America. Nature 2013;496: 50–6. Skirnisdottir S, Hreggvidsson GO, Hjörleifsdottir S et al. Influence of sulfide and temperature on species composition and community structure of hot spring microbial mats. Appl Environ Microb 2000;66:2835–41. Spear JR, Walker JJ, McCollom TM et al. Hydrogen and bioenergetics in the Yellowstone geothermal ecosystem. P Natl Acad Sci USA 2005;102:2555–60. Takacs-Vesbach C, Mitchell K, Jackson-Weaver O et al. Volcanic calderas delineate biogeographic provinces among Yellowstone thermophiles. Environ Microbiol 2008;10:1681–9. Tamura K, Peterson D, Peterson N et al. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol 2011;28:2731–9. Tekere M, Prinsloo A, Olivier J et al. An evaluation of the bacterial diversity at Tshipise, Mphephu and Sagole hot water springs, Limpopo Province, South Africa. Afr J Microbiol Res 2012;6:4993–5004. Tice MM, Lowe DR. Photosynthetic microbial mats in the 3,416Myr-old ocean. Nature 2004;431:549–52. Tindall BJ, Rosselló-Mora R, Busse H-J et al. Notes on the characterization of prokaryote strains for taxonomic purposes. Int J Syst Evol Micr 2010;60:249–66. Tobler DJ, Benning LG. Bacterial diversity in five Icelandic geothermal waters: temperature and sinter growth rate effects. Extremophiles 2011;15:473–85. Tomova I, Stoilova-Disheva M, Lyutskanova D et al. Phylogenetic analysis of the bacterial community in a geothermal 7 spring, Rupi Basin, Bulgaria. World J Microb Biot 2010;26: 2019–28. Van der Meer D, Torsvik T, Spakman W et al. Intra-Panthalassa Ocean subduction zones revealed by fossil arcs and mantle structure. Nat Geosci 2012;5:215–9. van der Meer MT, Klatt CG, Wood J et al. Cultivation and genomic, nutritional, and lipid biomarker characterization of Roseiflexus strains closely related to predominant in situ populations inhabiting Yellowstone hot spring microbial mats. J Bacteriol 2010;192:3033–42. van der Meer MT, Lammerts L, Skirnisdottir S et al. Distribution and isotopic composition of bacterial lipid biomarkers in microbial mats from a sulfidic Icelandic hot spring. Org Geochem 2008;39:1015–19. van der Meer MT, Schouten S, Bateson MM et al. Diel variations in carbon metabolism by green nonsulfur-like bacteria in alkaline siliceous hot spring microbial mats from Yellowstone National Park. Appl Environ Microb 2005;71: 3978–86. van Gremberghe I, Leliaert F, Mergeay J et al. Lack of phylogeographic structure in the freshwater cyanobacterium Microcystis aeruginosa suggests global dispersal. PLoS One 2011;6: 1–12. Vrijenhoek RC. Genetic diversity and connectivity of deep-sea hydrothermal vent metapopulations. Mol Ecol 2010;19:4391– 11. Whitaker RJ, Grogan DW, Taylor JW. Geographic barriers isolate endemic populations of hyperthermophilic archaea. Science 2003;301:976–8. White L, Lister G. The collision of India with Asia. J Geodyn 2012;56:7–17.