Flip‐flop organization in the chloroplast genome of Capsosiphon fulvescens (Ulvophyceae,Chlorophyta)

To better understand organelle genome evolution of the ulvophycean green alga Capsosiphon fulvescens, we sequenced and characterized its complete chloroplast genome. The circular chloroplast genome was 111,561 bp in length with 31.3% GC content that contained 108 genes including 77 protein‐coding genes, two copies of rRNA operons, and 27 tRNAs. In this analysis, we found the two types of isoform, called heteroplasmy, were likely caused by a flip‐flop organization. The flip‐flop mechanism may have caused structural variation and gene conversion in the chloroplast genome of C. fulvescens. In a phylogenetic analysis based on all available ulvophycean chloroplast genome data, including a new C. fulvescens genome, we found three major conflicting signals for C. fulvescens and its sister taxon Pseudoneochloris marina within 70 individual genes: (i) monophyly with Ulotrichales, (ii) monophyly with Ulvales, and (iii) monophyly with the clade of Ulotrichales and Ulvales. Although the 70‐gene concatenated phylogeny supported monophyly with Ulvales for both species, these complex phylogenetic signals of individual genes need further investigations using a data‐rich approach (i.e., organelle genome data) from broader taxon sampling.     

Nelumbonaceae:Systematic position and species diversification revealed by the complete chloroplast genome

<正>Nelumbonaceae is a morphologically unique family of angiosperms and was traditionally placed in Nymphaeales;more recently,it was placed in Proteales based on molecular data,or in an order of its own,Nelumbonales. To determine the systematic position of the family and to date the divergence time of the family and the divergence time of its two intercontinentally disjunct species,we sequenced the entire chloroplast genome of Nelumbo lutea and most of the chloroplast genes of,N.nucifera.We carried out phylogenetic and molecular dating analyses of the two species and representatives of 47 other plant families,representing the major lineages of angiosperms, using 83 plastid genes.The N.lutea genome was 163 510 bp long,with a total of 130 coding genes and an overall GC content of 38%.No significant structural differences among the genomes of N.lutea,Nymphaea alba, and Platanus occidentalis were observed.The phylogenetic relationships based on the 83 plastid genes revealed a close relationship between Nelumbonaceae and Platanaceae.The divergence times were estimated to be 109 Ma between the two families and 1.5 Ma between the two Nelumbo species.The estimated time was only slightly longer than the age of known Nelumbo fossils,suggesting morphological stasis within Nelumbonaceae.We conclude that Nelumbonaceae holds a position in or close to Proteales.We further conclude that the two species of Nelumbo diverged recently from a common ancestor and do not represent ancient relicts on different continents.     

The complete chloroplast genome of Gracilariopsis lemaneiformis (Rhodophyta) gives new insight into the evolution of family Gracilariaceae

The complete chloroplast genome of Gracilariopsis lemaneiformis was recovered from a Next Generation Sequencing data set. Without quadripartite structure, this chloroplast genome (183,013 bp, 27.40% GC content) contains 202 protein‐coding genes, 34 tRNA genes, 3 rRNA genes, and 1 tmRNA gene. Synteny analysis showed plasmid incorporation regions in chloroplast genomes of three species of family Gracilariaceae and in Grateloupia taiwanensis of family Halymeniaceae. Combined with reported red algal plasmid sequences in nuclear and mitochondrial genomes, we postulated that red algal plasmids may have played an important role in ancient horizontal gene transfer among nuclear, chloroplast, and mitochondrial genomes. Substitution rate analysis showed that purifying selective forces maintaining stability of protein‐coding genes of nine red algal chloroplast genomes over long periods must be strong and that the forces acting on gene groups and single genes of nine red algal chloroplast genomes were similar and consistent. The divergence of Gp. lemaneiformis occurred ~447.98 million years ago (Mya), close to the divergence time of genus Pyropia and Porphyra (443.62 Mya).     

The distribution of the thermally tolerant symbiont lineage (Symbiodinium clade D) in corals from Hawaii: correlations with host and the history of ocean thermal stress

Spatially intimate symbioses, such as those between scleractinian corals and unicellular algae belonging to the genus Symbiodinium, can potentially adapt to changes in the environment by altering the taxonomic composition of their endosymbiont communities. We quantified the spatial relationship between the cumulative frequency of thermal stress anomalies (TSAs) and the taxonomic composition of Symbiodinium in the corals Montipora capitata, Porites lobata, and Porites compressa across the Hawaiian archipelago. Specifically, we investigated whether thermally tolerant clade D Symbiodinium was in greater abundance in corals from sites with high frequencies of TSAs. We recovered 2305 Symbiodinium ITS2 sequences from 242 coral colonies in lagoonal reef habitats at Pearl and Hermes Atoll, French Frigate Shoals, and Kaneohe Bay, Oahu in 2007. Sequences were grouped into 26 operational taxonomic units (OTUs) with 12 OTUs associated with Montipora and 21 with Porites. Both coral genera associated with Symbiodinium in clade C, and these co‐occurred with clade D in M. capitata and clade G in P. lobata. The latter represents the first report of clade G Symbiodinium in P. lobata. In M. capitata (but not Porites spp.), there was a significant correlation between the presence of Symbiodinium in clade D and a thermal history characterized by high cumulative frequency of TSAs. The endogenous community composition of Symbiodinium and an association with clade D symbionts after long‐term thermal disturbance appear strongly dependent on the taxa of the coral host.     

Comparative phylogenetic analyses of Chinese Horsfieldia (Myristicaceae) using complete chloroplast genome sequences

The biologist's ruler for biodiversity is the species; accurate species identification is fundamental to the conservation of endangered species and in-depth biological scientific exploration. However, the delimitation and affinities of Horsfieldia in China has been controversial, owing in part to very low levels of molecular divergence within the family Myristicaceae. Because species boundaries and phylogenetic relationships within Horsfieldia are also unclear, 13 samples were collected across its distribution in China and their genomes were subjected to shotgun sequencing using Illumina platforms. A total of 40 487 994–84 801 416 pair-end clean reads were obtained and, after assembly, the complete chloroplast genome was recovered for all samples. Annotation analysis revealed a total of 112 genes, including 78 protein-coding genes, 30 transfer RNA, and 4 ribosomal RNA genes. Six variable loci (petN-psbM, trnH-psbA, ndhC-trnV, psbJ-psbL, ndhF, and rrn5-rrn23) were identified. Phylogenetic analyses strongly support the presence of four distinct species of Horsfieldia in China. In addition, samples that had been identified previously as Horsfieldia kingii (Hook. f.) Warb. were indistinguishable from those of H. prainii (King) Warb., suggesting that if H. kingii does occur in China, it was not collected in this study. Similarly, the chloroplast genome of one H. hainanensis Merr. sample from Guangxi province was identical to H. tetratepala C. Y. Wu, suggesting that the distribution range of H. hainanensis might be narrower than assumed previously. The phylogenetic relationships between the Chinese Horsfieldia species based on the whole chloroplast genomes was supported strongly, indicating the potential for using entire chloroplast genomes as super-barcodes for further resolution of the phylogeny of the genus Horsfieldia.     

Phylogenetic relationships within Adiaphanida (phylum Platyhelminthes) and the status of the crustacean‐parasitic genus Genostoma

The existence of the platyhelminth clade Adiaphanida—an assemblage comprising the well‐studied order Tricladida as well as two lesser known taxa, Prolecithophora and the obligate parasitic Fecampiida—is among the more surprising results of flatworm molecular systematics. Each of these three clades is itself largely well‐defined from a morphological point of view, although Adiaphanida at large, despite its strong support in molecular phylogenetic analyses, lacks known morphological synapomorphies. However, one taxon, the genus Genostoma, a parasite of the leptostracan crustacean Nebalia, rests uneasily within its current classification within the fecampiid family Genostomatidae; ultrastructural investigations on this taxon have uncovered a spermatogenesis reminiscent of Kalyptorhynchia, and a dorsal syncytium resembling the neodermatan tegument. Here, we provide molecular sequence data (nearly complete 18S and 28S rRNA) from a representative of Genostoma, with which we test hypotheses on the phylogenetic position of this taxon within Platyhelminthes, expanding upon a recently published phylum‐wide analysis, and applying novel alignment algorithms and substitution models. These analyses unequivocally position Genostoma as the sister group of Prolecithophora. However, even in taxon‐rich analyses, support for the position of the root of Adiaphanida is lacking, highlighting the need for new data types to study the phylogeny of this clade. Interestingly, our analyses also do not recover the monophyly of several taxa previously proposed, notably Continenticola within Tricladida and Protomonotresidae within Prolecithophora. In light of this phylogeny and the distinctive morphology (especially, spermatogenesis) of Genostoma, we advocate for a redefinition of the family Genostomatidae, outside of both Fecampiida and Prolecithophora, to encompass the members of this unique genus of parasites. Within Fecampiida, the family Piscinquilinidae fam. nov. is erected to accommodate the vertebrate‐parasitic Piscinquilinus, formerly Genostomatidae.     

Identification of chloroplast genome loci suitable for high‐resolution phylogeographic studies of Colocasia esculenta (L.) Schott (Araceae) and closely related taxa

Recently, we reported the chloroplast genome‐wide association of oligonucleotide repeats, indels and nucleotide substitutions in aroid chloroplast genomes. We hypothesized that the distribution of oligonucleotide repeat sequences in a single representative genome can be used to identify mutational hotspots and loci suitable for population genetic, phylogenetic and phylogeographic studies. Using information on the location of oligonucleotide repeats in the chloroplast genome of taro (Colocasia esculenta), we designed 30 primer pairs to amplify and sequence polymorphic loci. The primers have been tested in a range of intra‐specific to intergeneric comparisons, including ten taro samples (Colocasia esculenta) from diverse geographical locations, four other Colocasia species (C. affinis, C. fallax, C. formosana, C. gigantea) and three other aroid genera (represented by Remusatia vivipara, Alocasia brisbanensis and Amorphophallus konjac). Multiple sequence alignments for the intra‐specific comparison revealed nucleotide substitutions (point mutations) at all 30 loci and microsatellite polymorphisms at 14 loci. The primer pairs reported here reveal levels of genetic variation suitable for high‐resolution phylogeographic and evolutionary studies of taro and other closely related aroids. Our results confirm that information on repeat distribution can be used to identify loci suitable for such studies, and we expect that this approach can be used in other plant groups.     

The chloroplast genome of Phacus orbicularis (Euglenophyceae): an initial datum point for the phacaceae

The Euglenophyceae chloroplast was acquired when a heterotrophic euglenoid engulfed a green alga and subsequently retained the algal chloroplast, in a process known as secondary endosymbiosis. Since this event, Euglenophyceae have diverged widely and their chloroplast genomes (cpGenomes) have as well. Changes to the cpGenome include extensive gene rearrangement and the proliferation of introns, the analyses of which have proven to be useful in examining cpGenome changes throughout the Euglenophyceae. The Euglenales fall into two families, Euglenaceae and Phacaceae. Euglenaceae contains eight genera and at least one cpGenome has been published for each genus. Phacaceae, on the other hand, contains three genera, none of which have had a representative chloroplast genome sequenced. Members of this family have many small disk‐shaped chloroplasts that lack pyrenoids. We sequenced and annotated the cpGenome of Phacus orbicularis in order to fill in the large gap in our understanding of Euglenophyceae cpGenome evolution, especially in regard to intron number and gene order. We compared this cpGenome to those of species from both the Euglenaceae and Eutreptiales of the Euglenophyceae phylogenetic tree. The cpGenome showed characteristics that were more derived than that of the basal species Eutreptia viridis, with extensive gene rearrangements and nearly three times as many introns. In contrast, it contained fewer introns than all but one of the previously reported Euglenaceae cpGenomes, had a smaller estimated genome size, and shared greater synteny with two main branches of that family.     

Corals in the hottest reefs in the world exhibit symbiont fidelity not flexibility

Reef‐building corals are at risk of extinction from ocean warming. While some corals can enhance their thermal limits by associating with dinoflagellate photosymbionts of superior stress tolerance, the extent to which symbiont communities will reorganize under increased warming pressure remains unclear. Here we show that corals in the hottest reefs in the world in the Persian Gulf maintain associations with the same symbionts across 1.5 years despite extreme seasonal warming and acute heat stress (≥35°C). Persian Gulf corals predominantly associated with Cladocopium (clade C) and most also hosted Symbiodinium (clade A) and/or Durusdinium (clade D). This is in contrast to the neighbouring and milder Oman Sea, where corals associated with Durusdinium and only a minority hosted background levels of Cladocopium. During acute heat stress, the higher prevalence of Symbiodinium and Durusdinium in bleached versus nonbleached Persian Gulf corals indicates that genotypes of these background genera did not confer bleaching resistance. Within symbiont genera, the majority of ITS2 rDNA type profiles were unique to their respective coral species, confirming the existence of host‐specific symbiont lineages. Notably, further differentiation among Persian Gulf sites demonstrates that symbiont populations are either isolated or specialized over tens to hundreds of kilometres. Thermal tolerance across coral species was associated with the prevalence of a single ITS2 intragenomic sequence variant (C3gulf), definitive of the Cladocopium thermophilum group. The abundance of C3gulf was highest in bleaching‐resistant corals and at warmer sites, potentially indicating a specific symbiont genotype (or set of genotypes) that may play a role in thermal tolerance that warrants further investigation. Together, our findings indicate that co‐evolution of host–Symbiodiniaceae partnerships favours fidelity rather than flexibility in extreme environments and under future warming.     

Taxonomy and biogeography of Diapensia (Diapensiaceae) based on chloroplast genome data

Diapensia L. is the second largest genus of Diapensiaceae. The taxonomic treatment within Diapensia and relationships within Diapensiaceae have been disputed. Chloroplast genome sequence data have proved to be useful for plant phylogenetic analyses and species delimitation. In this study, we de novo sequenced and assembled 22 chloroplast genomes of 15 species of Diapensiaceae, including all recognized species of Diapensia with multiple samples. A super‐matrix containing a total of 107 genes and 18 taxa was constructed for phylogenetic analyses to resolve phylogenetic relationships among genera of the family and within Diapensia. The resulting phylogenetic tree showed the following strongly supported relationships: (Galax, (Pyxidanthera, (Berneuxia, ((Schizocodon, Diapensia), and Shortia s.s.)))). The dated phylogeny and reconstructed lineage‐through‐time plot for the family indicated rapid diversification in the Neogene and an acceleration of diversification rate after c. 8 Ma. Biogeographic analysis suggested that Diapensia originated in the Northeast Asian mountains approximately 6.06 Ma, followed by northward dispersal to the Arctic and southwestward dispersal to the Himalaya–Hengduan Mountains. Phylogenetic relationships within Diapensia were well resolved. Based on the phylogenetic results, we proposed to reinstate the species status of Diapensia bulleyana Forrest ex Diels, and raised D. purpurea f. albida to the species rank (D. albida [W. E. Evans] J. F. Ye comb. & stat. nov.). The distribution ranges of all species delineated based on the phylogenetic results were revised accordingly based on specimen occurrences. Our study adds new examples for the power of plastid genome data for resolving phylogenetic relationships and clarifying taxonomic disputes among closely related species.     

Comparative DNA sequence analyses of Pyramimonas parkeae (Prasinophyceae) chloroplast genomes

Prasinophytes form a paraphyletic assemblage of early diverging green algae, which have the potential to reveal the traits of the last common ancestor of the main two green lineages: (i) chlorophyte algae and (ii) streptophyte algae. Understanding the genetic composition of prasinophyte algae is fundamental to understanding the diversification and evolutionary processes that may have occurred in both green lineages. In this study, we sequenced the chloroplast genome of Pyramimonas parkeae NIES254 and compared it with that of P. parkeae CCMP726, the only other fully sequenced P. parkeae chloroplast genome. The results revealed that P. parkeae chloroplast genomes are surprisingly variable. The chloroplast genome of NIES254 was larger than that of CCMP726 by 3,204 bp, the NIES254 large single copy was 288 bp longer, the small single copy was 5,088 bp longer, and the IR was 1,086 bp shorter than that of CCMP726. Similarity values of the two strains were almost zero in four large hot spot regions. Finally, the strains differed in copy number for three protein‐coding genes: ycf20, psaC, and ndhE. Phylogenetic analyses using 16S and 18S rDNA and rbcL sequences resolved a clade consisting of these two P. parkeae strains and a clade consisting of these plus other Pyramimonas isolates. These results are consistent with past studies indicating that prasinophyte chloroplast genomes display a higher level of variation than is commonly found among land plants. Consequently, prasinophyte chloroplast genomes may be less useful for inferring the early history of Viridiplantae than has been the case for land plant diversification.     

The endemic Cladophorales (Ulvophyceae) of ancient Lake Baikal represent a monophyletic group of very closely related but morphologically diverse species

Lake Baikal, the oldest lake in the world, is home to spectacular biodiversity and extraordinary levels of endemism. While many of the animal species flocks from Lake Baikal are famous examples of evolutionary radiations, the lake also includes a wide diversity of endemic algae that are not well investigated with regards to molecular‐biological taxonomy and phylogeny. The endemic taxa of the green algal order Cladophorales show a range of divergent morphologies that led to their classification in four genera in two families. We sequenced partial large‐ and small‐subunit rDNA as well as the internal transcribed spacer region of 14 of the 16 described endemic taxa to clarify their phylogenetic relationships. One endemic morphospecies, Cladophora kusnetzowii, was shown to be conspecific with the widespread Aegagropila linnaei. All other endemic morphospecies formed a monophyletic group nested within the genus Rhizoclonium (Cladophoraceae), a very surprising result, in stark contrast to their morphological affinities. The Baikal clade represents a species flock of closely related taxa with very low genetic differentiation. Some of the morphospecies were congruent with lineages recovered in the phylogenies, but due to the low phylogenetic signal in the rDNA sequences the relationships within the Baikal clade were not all well resolved. The Baikal clade appears to represent a recent radiation, based on the low molecular divergence within the group, and it is hypothesized that the large morphological variation results from diversification in sympatry from a common ancestor in Lake Baikal.     

The chloroplast genome sequence of Chara vulgaris sheds new light into the closest green algal relatives of land plants

The phylum Streptophyta comprises all land plants and six monophyletic groups of charophycean green algae (Mesostigmatales, Chlorokybales, Klebsormidiales, Zygnematales, Coleochaetales, and Charales). Phylogenetic analyses of four genes encoded in three cellular compartments suggest that the Charales are sister to land plants and that charophycean green algae evolved progressively toward an increasing cellular complexity. To validate this phylogenetic hypothesis and to understand how and when the highly conservative pattern displayed by land plant chloroplast DNAs (cpDNAs) originated in the Streptophyta, we have determined the complete chloroplast genome sequence (184,933 bp) of a representative of the Charales, Chara vulgaris, and compared this genome to those of Mesostigma (Mesostigmatales), Chlorokybus (Chlorokybales), Staurastrum and Zygnema (Zygnematales), Chaetosphaeridium (Coleochaetales), and selected land plants. The phylogenies we inferred from 76 cpDNA-encoded proteins and genes using various methods favor the hypothesis that the Charales diverged before the Coleochaetales and Zygnematales. The Zygnematales were identified as sister to land plants in the best tree topology (T1), whereas Chaetosphaeridium (T2) or a clade uniting the Zygnematales and Chaetosphaeridium (T3) occupied this position in alternative topologies. Chara remained at the same basal position in trees including more land plant taxa and inferred from 56 proteins/genes. Phylogenetic inference from gene order data yielded two most parsimonious trees displaying the T1 and T3 topologies. Analyses of additional structural cpDNA features (gene order, gene content, intron content, and indels in coding regions) provided better support for T1 than for the topology of the above-mentioned four-gene tree. Our structural analyses also revealed that many of the features conserved in land plant cpDNAs were inherited from their green algal ancestors. The intron content data predicted that at least 15 of the 21 land plant group II introns were gained early during the evolution of streptophytes and that a single intron was acquired during the transition from charophycean green algae to land plants. Analyses of genome rearrangements based on inversions predicted no alteration in gene order during the transition from charophycean green algae to land plants.     

Ocean acidification and warming scenarios increase microbioerosion of coral skeletons

Biological mediation of carbonate dissolution represents a fundamental component of the destructive forces acting on coral reef ecosystems. Whereas ocean acidification can increase dissolution of carbonate substrates, the combined impact of ocean acidification and warming on the microbioerosion of coral skeletons remains unknown. Here, we exposed skeletons of the reef‐building corals, Porites cylindrica and Isopora cuneata, to present‐day (Control: 400 μatm – 24 °C) and future pCO₂–temperature scenarios projected for the end of the century (Medium: +230 μatm – +2 °C; High: +610 μatm – +4 °C). Skeletons were also subjected to permanent darkness with initial sodium hypochlorite incubation, and natural light without sodium hypochlorite incubation to isolate the environmental effect of acidic seawater (i.e., Ω_aragonite <1) from the biological effect of photosynthetic microborers. Our results indicated that skeletal dissolution is predominantly driven by photosynthetic microborers, as samples held in the dark did not decalcify. In contrast, dissolution of skeletons exposed to light increased under elevated pCO₂–temperature scenarios, with P. cylindrica experiencing higher dissolution rates per month (89%) than I. cuneata (46%) in the high treatment relative to control. The effects of future pCO₂–temperature scenarios on the structure of endolithic communities were only identified in P. cylindrica and were mostly associated with a higher abundance of the green algae Ostreobium spp. Enhanced skeletal dissolution was also associated with increased endolithic biomass and respiration under elevated pCO₂–temperature scenarios. Our results suggest that future projections of ocean acidification and warming will lead to increased rates of microbioerosion. However, the magnitude of bioerosion responses may depend on the structural properties of coral skeletons, with a range of implications for reef carbonate losses under warmer and more acidic oceans.     

Mitochondrial genomes of the green macroalga Ulva pertusa (Ulvophyceae,Chlorophyta): novel insights into the evolution of mitogenomes in the Ulvophyceae

To further understand the trends in the evolution of mitochondrial genomes (mitogenomes or mtDNAs) in the Ulvophyceae, the mitogenomes of two separate thalli of Ulva pertusa were sequenced. Two U. pertusa mitogenomes (Up1 and Up2) were 69,333 bp and 64,602 bp in length. These mitogenomes shared two ribosomal RNAs (rRNAs), 28 transfer RNAs (tRNAs), 29 protein‐coding genes, and 12 open reading frames. The 4.7 kb difference in size was attributed to variation in intron content and tandem repeat regions. A total of six introns were present in the smaller U. pertusa mtDNA (Up2), while the larger mtDNA (Up1) had eight. The larger mtDNA had two additional group II introns in two genes (cox1 and cox2) and tandem duplication mutations in noncoding regions. Our results showed the first case of intraspecific variation in chlorophytan mitogenomes and provided further genomic data for the undersampled Ulvophyceae.     

A rare case of plastid protein‐coding gene duplication in the chloroplast genome of Euglena archaeoplastidiata (Euglenophyta)

Gene duplication is an important evolutionary process that allows duplicate functions to diverge, or, in some cases, allows for new functional gains. However, in contrast to the nuclear genome, gene duplications within the chloroplast are extremely rare. Here, we present the chloroplast genome of the photosynthetic protist Euglena archaeoplastidiata. Upon annotation, it was found that the chloroplast genome contained a novel tandem direct duplication that encoded a portion of RuBisCO large subunit (rbcL) followed by a complete copy of ribosomal protein L32 (rpl32), as well as the associated intergenic sequences. Analyses of the duplicated rpl32 were inconclusive regarding selective pressures, although it was found that substitutions in the duplicated region, all non‐synonymous, likely had a neutral functional effect. The duplicated region did not exhibit patterns consistent with previously described mechanisms for tandem direct duplications, and demonstrated an unknown mechanism of duplication. In addition, a comparison of this chloroplast genome to other previously characterized chloroplast genomes from the same family revealed characteristics that indicated E. archaeoplastidiata was probably more closely related to taxa in the genera Monomorphina, Cryptoglena, and Euglenaria than it was to other Euglena taxa. Taken together, the chloroplast genome of E. archaeoplastidiata demonstrated multiple characteristics unique to the euglenoid world, and has justified the longstanding curiosity regarding this enigmatic taxon.     

Ecological Inferences from a deep screening of the Complex Bacterial Consortia associated with the coral,Porites astreoides

The functional role of the bacterial organisms in the reef ecosystem and their contribution to the coral well‐being remain largely unclear. The first step in addressing this gap of knowledge relies on in‐depth characterization of the coral microbial community and its changes in diversity across coral species, space and time. In this study, we focused on the exploration of microbial community assemblages associated with an ecologically important Caribbean scleractinian coral, Porites astreoides, using Illumina high‐throughput sequencing of the V5 fragment of 16S rRNA gene. We collected data from a large set of biological replicates, allowing us to detect patterns of geographical structure and resolve co‐occurrence patterns using network analyses. The taxonomic analysis of the resolved diversity showed consistent and dominant presence of two OTUs affiliated with the order Oceanospirillales, which corroborates a specific pattern of bacterial association emerging for this coral species and for many other corals within the genus Porites. We argue that this specific association might indicate a symbiotic association with the adult coral partner. Furthermore, we identified a highly diverse rare bacterial ‘biosphere’ (725 OTUs) also living along with the dominant bacterial symbionts, but the assemblage of this biosphere is significantly structured along the geographical scale. We further discuss that some of these rare bacterial members show significant association with other members of the community reflecting the complexity of the networked consortia within the coral holobiont.