Slow algae, fast fungi: exceptionally high nucleotide substitution rate differences between lichenized fungi Omphalina and their symbiotic green algae Coccomyxa

Omphalina basidiolichens are obligate mutualistic associations of a fungus of the genus Omphalina (the exhabitant) and a unicellular green alga of the genus Coccomyxa (the inhabitant). It has been suggested that symbiotic inhabitants have a lower rate of genetic change compared to exhabitants because the latter are more exposed to abiotic environmental variation and competition from other organisms. In order to test this hypothesis we compared substitution rates in the nuclear ribosomal internal transcribed spacer region (ITS1, 5.8S, ITS2) among fungal species with rates among their respective algal symbionts. To ensure valid comparisons, only taxon pairs (12) with a common evolutionary history were used. On average, substitution rates in the ITS1 portion of Omphalina pairs were 27.5 times higher than rates in the corresponding pairs of Coccomyxa since divergence from their respective ancestor at the base of the Omphalina/Coccomyxa lineage. Substitution rates in the 5.8S and the ITS2 portions were 2.4 and 18.0 times higher, respectively. The highest rate difference (43.0) was found in the ITS1 region. These are, to our knowledge, the highest differences of substitution rates reported for symbiotic organisms. We conclude that the Omphalina model system conforms to the proposed hypothesis of lower substitution rates in the inhabitant, but that the mode of transmission of the inhabitant (vertical versus horizontal) could be a prevailing factor in the regulation of unequal rates of nucleotide substitution between co-evolving symbionts. Our phylogenetic study of Coccomyxa revealed three main lineages within this genus, corresponding to free-living Coccomyxa, individuals isolated from basidiolichens Omphalina and Coccomyxa isolated from ascolichens belonging to the Peltigerales.     

Multiple origins of the symbioses in Paramecium bursaria

Many organisms have symbioses with photosynthetic algae as typified by corals, clams, lichens, and some protozoa. Paramecium bursaria contains green algal symbionts and this unicellular ciliate is a textbook example used for microscopic observation in junior high school science projects. We have determined molecular phylogenies for the green algal symbionts. The symbiotic algae are the main constituent of the Paramecium cytoplasm, and we have recognized a total of four species, of which two were newly discovered in the present study. One should be regarded genetically as Chlorella vulgaris, and it belongs phylogenetically to the Chlorella clade (Chlorellaceae, Trebouxiophyceae) as well as "American" and "European" groups, which we previously introduced. Their genetic dissimilarities are 0.50-0.83% in 18S rDNA comparisons, but those of the internal transcribed spacer 2 (ITS2) reach an unambiguous level (22.6-26.6%). These dissimilarities suggest that they are equivalent to discrete species derived from multiple origins as paramecian symbionts. Another newcomer was clearly separated from the Chlorellaceae, and this alga clustered with Coccomyxa spp. in ITS2 analyses. These symbiotic relations indicate multiple origins of symbionts.     

Phylogenetic position of endosymbiotic green algae in Paramecium bursaria Ehrenberg from Japan

Endosymbiotic green algae of Japanese Paramecium bursaria were phylogenetically analyzed based on DNA sequences from the ribosomal DNA operon (18S rDNA, ITS1, 5.8S rDNA, and ITS2). Phylogenetic trees constructed using 18S rDNA sequences showed that the symbionts belong to the Chlorella sensu stricto (Trebouxiophyceae) group. They are genetically closer to the C. vulgaris Beijerinck group than to C. kessleri Fott et Nováková as proposed previously. Branching order in C. vulgaris group was unresolved in 18S rDNA trees. Compared heterogeneities of 18S rDNA, ITS1, 5.8S r, and ITS2 among symbionts and two Chlorella species, indicated that the ITS2 region (and probably also ITS1) is better able to resolve phylogenetic problems in such closely related taxa. All six symbiotic sequences obtained here (approximately 4000-bp sequences of 18S rDNA, ITS1, 5.8S rDNA, and ITS2) were completely identical in each, strongly suggesting a common origin.     

Phylogeny of lichen- and non-lichen-forming omphalinoid mushrooms and the utility of testing for combinability among multiple data sets

As an initial step toward developing a model system to study requirements for and consequences of transitions to mutualism, the phylogeny of a group of closely related lichenized and nonlichenized basidiomycetes (Omphalina) was reconstructed. The phylogenetic analyses are based on four data sets representing different regions of the nuclear ribosomal repeat unit (ITS1, 5.8S, ITS2, and 25S) obtained from 30 species of Omphalina and related genera. The resulting phylogenetic trees from each of these four data sets, when analyzed separately, were not identical. Testing for the combinability of these four data sets suggested that they could not be combined in their entirety. The removal of ambiguous alignments and saturated sites was sufficient, after reapplying the combinability test on the pruned data sets, to explain the topological discrepancies. In this process, the first of two complementary tests developed by Rodrigo et al. (1993, N.Z. J. Bot. 31:257-268) to assess whether two data sets are the result of the same phylogenetic history was found to be biased, rejecting the combinability of two data sets even when they are samples of the same phylogenetic history. Combining the four pruned data sets yielded phylogenies that suggest the five lichen-forming species of Omphalina form a monophyletic group. The sister group to this symbiotic clade consists mostly of dark brown Omphalina species intermixed with species from the genera Arrhenia and Phaeothellus. The genera Omphalina and Gerronema are shown to be polyphyletic. The lichen-forming species O. ericetorum and the nonmutualistic species O. velutipes, O. epichysium, and O. sphagnicola are the best candidates for experimental work designed to gain a better understanding of mechanisms involved in symbiotic interactions and the role symbiosis has played in the evolution of fungi.     

Phylogenetic relationships in Termitomyces (Family Agaricaceae) based on the nucleotide sequence of ITS: a first approach to elucidate the evolutionary history of the symbiosis between fungus-growing termites and their fungi

Termitomyces constitutes a very poorly known genus of fungi whose essential characteristic is that all representatives of the genus are cultivated by termites (Macrotermitinae) in their nest and that all the fungi cultivated by termites belong to this genus. For the first time, the phylogenetic relationships of several African Termitomyces species was studied by the sequencing of their internal transcriber spacer region (ITS1--5.8S--ITS2). It appeared that this group is clearly monophyletic and belongs to the Tricholomataceae family. The total homology of the ITS zone of several Termitomyces symbionts of different termite genera indicated that the specific diversity of this group is in fact less important than previously supposed. Finally, the comparison between the Termitomyces phylogenetic tree and the taxonomic tree of Macrotermitinae showed that if for certain genera the hypothesis of termite/fungus coevolution is acceptable, it should not be applied for all symbiosis.     

Genetic variation within <Emphasis Type="Italic">Symbiodinium</Emphasis> clade B from the coral genus <Emphasis Type="Italic">Madracis</Emphasis> in the Caribbean (Netherlands Antilles)

The internal transcribed spacer (ITS) region was sequenced in symbiotic dinoflagellates (zooxanthellae) from five morphospecies in the genus Madracis. The phylogeny of the symbionts is congruent with a companion phylogeny of the coral host. Comparison with known clade B symbiont ITS types reveals that M. mirabilis contains the B13 symbiont and that the other morphospecies contain the B7 symbiont. Madracis formosa also contains a previously undescribed type. The B7 and B13 symbionts appear to be highly specific to morphospecies in the genus Madracis. The host specificity between the B13 symbionts and its coral host may be the result of co-evolution of the coral-symbiont association and suggests that the brooding species, M. mirabilis, is reproductively isolated. Microhabitat differentiation associated with light utilization independent of depth is discussed.     

Evolutionary Rates in Veronica L. (Plantaginaceae): Disentangling the Influence of Life History and Breeding System

The evolutionary rate at which DNA sequences evolve is known to differ between different groups of organisms. However, the reasons for these different rates are seldom known. Among plants, the generation-time hypothesis, which states that organisms that reproduce faster also have more DNA substitutions per time, has gained most popularity. We evaluate the generation-time hypothesis using 131 DNA sequences from the plastid trnLF region and the nuclear ribosomal ITS region of the genus Veronica (Plantaginaceae). We also examine the alternative hypothesis that a higher substitution rate is correlated with selfing breeding system. Selfing is associated with annual life history in many organisms and may thus often be the underlying reason for observed correlations of annual life history with other characters. We provide evidence that annual life history is more likely to be the responsible factor for higher substitution rates in Veronica than a selfing breeding system. Nevertheless, the way in which annual life history may influence substitution rate in detail remains unknown, and some possibilities are discussed.     

Phylogenetic relationships among Myrceugenia, Blepharocalyx, and Luma (Myrtaceae) based on paired-sites models and the secondary structures of ITS and ETS sequences

To establish relationships among the genera Blepharocalyx, Luma, and Myrceugenia, the phylogeny of tribe Myrteae was reconstructed based on secondary structures of the sequences of ITS (ITS1-5.8S-ITS2) and ETS regions from 93 taxa belonging to 29 genera. The DNA sequences were aligned according to their secondary structure and analysed with Bayesian inference (BI) using base substitution models for RNA, which can differentiate between the regions consisting of base pairs (helices) and unpaired bases (loops). The best-fit models were RNA7C for ITS and RNA7B for ETS, which differ in how they predict double substitutions and symmetry of the frequency of base pairs. The phylogenetic hypothesis was compared with results obtained using classical 4-state models, the results of maximum parsimony, and using sequence alignment without considering the secondary structure. Analyses show low support along the backbone of the consensus trees, those using substitution models of paired sites showing more highly supported derived clades than substitution models of independent sites. All tests were consistent and differed in the positioning of certain groups, but they also showed that Blepharocalyx, Luma, and Myrceugenia arise from three independent lineages that are not closely related. The substitution model for RNA shows that Luma is a monophyletic group and sister to the other genera of Myrteae, except for Myrtus communis. Myrceugenia appears as a monophyletic group, except for M. fernandeziana, which is related to Blepharocalyx, a genus that appears to be biphyletic.     

Population structure between environmentally transmitted vibrios and bobtail squids using nested clade analysis

Squids from the genus Euprymna (Cephalopoda: Sepiolidae) and their symbiotic bacteria Vibrio fischeri form a mutualism in which vibrios inhabit a complex light organ within the squid host. A host-mediated daily expulsion event seeds surrounding seawater with symbiotically capable V. fischeri that environmentally colonize newly hatched axenic Euprymna juveniles. Competition experiments using native and non-native Vibrio have shown that this expulsion/re-colonization phenomenon has led to cospeciation in this system in the Pacific Ocean; however, the genetic architecture of these symbiotic populations has not been determined. Using genetic diversity and nested clade analyses we have examined the variation and history of three allopatric Euprymna squid species (E. scolopes of Hawaii, E. hyllebergi of Thailand, and E. tasmanica from Australia) and their respective Vibrio symbionts. Euprymna populations appear to be very genetically distinct from each other, exhibiting little or no migration over large geographical distances. In contrast, Vibrio symbiont populations contain more diverse haplotypes, suggesting both host presence and unidentified factors facilitating long-distance migration structure in Pacific Vibrio populations. Findings from this study highlight the importance of how interactions between symbiotic organisms can unexpectedly shape population structure in phylogeographical studies.     

Molecular tools for isolate and community studies of Pyrenomycete fungi

The Pyrenomycetes, defined physiologically by the formation of a flask-shaped fruiting body present in the sexual form, are a monophyletic group of fungi that consist of a wide diversity of populations including human and plant pathogens. Based on sequence analysis of 18S ribosomal DNA (rDNA), rDNA regions conserved among the Pyrenomycetes but divergent among other organisms were identified and used to develop selective PCR primers and a highly specific primer set. The primers presented here were used to amplify large portions of the 18S rDNA as well as the entire internal transcribed spacer (ITS) region (ITS 1, 5.8S rDNA, and ITS 2). In addition to database searches, the specificity of the primers was verified by PCR amplification of DNA extracted from pure culture isolates and by sequence analysis of fungal rDNA PCR-amplified from environmental samples. In addition, denaturing gradient gel electrophoresis (DGGE) analyses were performed on closely related Colletotrichum isolates serving as a model pathogenic genus of the Pyrenomycetes. Although both ITS and 18S rDNA DGGE analyses of Colletotrichum were consistent with a phylogeny established from sequence analysis of the ITS region, DGGE analysis of the ITS region was found to be more sensitive than DGGE analysis of the 18S rDNA. This study introduces molecular tools for the study of Pyrenomycete fungi by the development of two specific primers, demonstration of the enhanced sensitivity of ITS-DGGE for typing of closely related isolates and application of these tools to environmental samples.     

Betaproteobacterial symbionts of the ciliate Euplotes: origin and tangled evolutionary path of an obligate microbial association

The Polynucleobacter-Euplotes association is an obligatory symbiotic system between a monophyletic group of ciliate species belonging to the genus Euplotes and bacteria of the species Polynucleobacter necessarius (Betaproteobacteria). Both organisms are unable to survive independently. Several studies revealed the existence of free-living populations of Polynucleobacter bacteria which are phylogenetically closely related to the endosymbiotic ones, but never share associations with Euplotes in the natural environment. Hence, following the most parsimonious explanation on the origin of the association, this symbiosis should represent a synapomorphic character for the hosts' clade. Nevertheless, phylogenetic analyses performed on an increased number of strains here presented suggest that Euplotes species, during their evolution, recruited Polynucleobacter bacteria as symbionts more than once. Moreover, in three cases, we observed different bacteria as obligate symbionts. These symbionts are the first characterized representatives of a phylogenetic lineage branching in a basal position with respect to the genus Polynucleobacter. The hypothesis that the original obligate symbionts belonged to this newly discovered clade and that, only subsequently, in most cases they have been replaced by Polynucleobacter bacteria recruited from the environment is proposed and discussed. The evolutionary path of this association seems anyway to have been more complex than so far supposed.     

Life at the Limits: Capacities of Isolated and Cultured Lichen Symbionts to Resist Extreme Environmental Stresses

Lichens are described as a symbiosis formed by a myco- and photobiont, capable of colonizing habitats where their separate symbionts would not be able to survive. Space simulation studies on the separated symbionts of the lichen Xanthoria elegans have been performed to test their capacity to resist the most extreme conditions. The isolated cultured symbiont cells were exposed to different doses of the UV spectrum, and to vacuum. Cultures of both symbionts were analysed by specific vitality tests (LIVE/DEAD-staining detected by Confocal Laser Scanning Microscopy). Growth capacity of symbiont cultures on different media was analysed after exposure to extreme environmental stresses. The data obtained support the hypothesis that the symbiotic state considerably enhances the ability of the respective symbionts to survive exposure to extreme conditions, including the conditions of space simulation. Species such as X. elegans may, therefore, be suitable for use as model organisms in exobiological studies.     

Identification of endophytic fungi from Livistona chinensis based on morphology and rDNA sequences

A survey of the endophytic fungi in fronds of Livistona chinensis was carried out in Hong Kong. The endophyte assemblages identified using morphological characters consisted of 16 named species and 19 'morphospecies', the latter grouped based on cultural morphology and growth rates. Arrangement of taxa into morphospecies does not reflect species phylogeny, and therefore selected morphospecies were further identified based on ribosomal DNA (rDNA) sequence analysis. The 5.8S gene and flanking internal transcribed spacers (ITS1 and ITS2) regions of rDNA from 19 representative morphospecies were amplified by the polymerase chain reaction and sequenced. Phylogenetic analysis based on 5.8S gene sequences showed that these morphospecies were filamentous Ascomycota, belonging in the Loculoascomycetes and Pyrenomycetes. Further identification was conducted by means of sequence comparison and phylogenetic analysis of both the ITS and 5.8S regions. Results showed that MS704 belonged to the genus Diaporthe and its anamorph Phomopsis of the Valsaceae. MS594 was inferred to be Mycosphaerella and its anamorph Cladosporium of the Mycosphaerellaceae. MS339, MS366, MS370, MS395, MS1033, MS1083 and MS1092 were placed in the genus Xylaria of the Xylariaceae. MS194, MS375 and MS1028 were close to the Clypeosphaeriaceae. MS191 and MS316 were closely related to the Pleosporaceae within the Dothideales. The other 5 morphospecies, MS786, MS1043, MS1065, MS1076 and MS1095, probably belong in the Xylariales. The value of using DNA sequence analysis in the identification of endophytes is discussed.     

Evidence for functional divergence in arbuscular mycorrhizal fungi from contrasting climatic origins

A considerable amount of phenotypic, genetic and symbiotic functional variability has been documented in arbuscular mycorrhizal fungi (AMF). However, little is known about whether distinct AMF ecotypes have evolved within their geographic range. We tested the hypothesis that AMF growing at temperatures closer to those prevalent within their origin would benefit their host and grow more than isolates distant from their native conditions. For each of six AMF species, we chose pairs of isolates that originated from distant areas with contrasting climates. Each isolate was grown in association with two grass species of different thermal optima at two temperature settings. Thus, we also tested whether AMF from different climatic origins were dependent on the thermal adaptation of the host plant species or to temperature per se. Although fungal growth was not directly affected by temperature, we found that AMF isolates originating from contrasting climates consistently and differentially altered plant growth. Our results suggest that AMF from contrasting climates have altered symbiotic function, thus linking an abiotic factor to ecotypic differentiation of putatively important symbionts.     

Genetic diversity of Nostoc microsymbionts from Gunnera tinctoria revealed by PCR-STRR fingerprinting

The cyanobacteria belonging to the genus Nostoc fix atmospheric nitrogen, both as free-living organisms and in symbiotic associations with a wide range of hosts, including bryophytes, gymnosperms (cycads), the small water fern Azolla (Pteridophyte), the angiosperm genus Gunnera, and fungi (lichens). The Gunnera–Nostoc symbiosis is the only one that involves a flowering plant. In Chile, 12 species of Gunnera have been described with a broad distribution in the temperate region. We examined the genetic diversity of Nostoc symbionts from three populations of Gunnera tinctoria from Abtao, Chiloé Island, southern Chile, and microsymbionts from other two species of Gunnera from southern Chile, using PCR amplification of STRR (short tandemly repeated repetitive) sequences of the Nostoc infected tissue. To our knowledge, this is the first report of PCR fingerprinting obtained directly from symbiotic tissue of Gunnera. Genetic analyses revealed that Nostoc symbionts exhibit important genetic diversity among host plants, both within and between Gunnera populations. It was also found that only one Nostoc strain, or closely related strains, established symbiosis with an individual plant host.     

CONSPECIFICITY AND INDO-PACIFIC DISTRIBUTION OF SYMBIODINIUM GENOTYPES (DINOPHYCEAE) FROM GIANT CLAMS

We previously reported the occurrence of genetically‐diverse symbiotic dinoflagellates (zooxanthellae) within and between 7 giant clam species (Tridacnidae) from the Philippines based on the algal isolates' allozyme and random amplified polymorphic DNA (RAPD) patterns. We also reported that these isolates all belong to clade A of the Symbiodinium phylogeny with identical 18S rDNA sequences. Here we extend the genetic characterization of Symbiodinium isolates from giant clams and propose that they are conspecific. We used the combined DNA sequences of the internal transcribed spacer (ITS)1, 5.8S rDNA, and ITS2 regions (rDNA‐ITS region) because the ITS1 and ITS2 regions evolve faster than 18S rDNA and have been shown to be useful in distinguishing strains of other dinoflagellates. DGGE of the most variable segment of the rDNA‐ITS region, ITS1, from clonal representatives of clades A, B, and C showed minimal intragenomic variation. The rDNA‐ITS region shows similar phylogenetic relationships between Symbiodinium isolates from symbiotic bivalves and some cnidarians as does 18S rDNA, and that there are not many different clade A species or strains among cultured zooxanthellae (CZ) from giant clams. The CZ from giant clams had virtually identical sequences, with only a single nucleotide difference in the ITS2 region separating two groups of isolates. These data suggest that there is one CZ species and perhaps two CZ strains, each CZ strain containing individuals that have diverse allozyme and RAPD genotypes. The CZ isolated from giant clams from different areas in the Philippines (21 isolates, 7 clam species), the Australian Great Barrier Reef (1 isolate, 1 clam species), Palau (8 isolates, 7 clam species), and Okinawa, Japan (1 isolate, 1 clam species) shared the same rDNA‐ITS sequences. Furthermore, analysis of fresh isolates from giant clams collected from these geographical areas shows that these bivalves also host indistinguishable clade C symbionts. These data demonstrate that conspecific Symbiodinium genotypes, particularly clade A symbionts, are distributed in giant clams throughout the Indo‐Pacific.     

Geographic mosaic of symbiont selectivity in a genus of epiphytic cyanolichens

In symbiotic systems, patterns of symbiont diversity and selectivity are crucial for the understanding of fundamental ecological processes such as dispersal and establishment. The lichen genus Nephroma (Peltigerales, Ascomycota) has a nearly cosmopolitan distribution and is thus an attractive model for the study of symbiotic interactions over a wide range of spatial scales. In this study, we analyze the genetic diversity of Nephroma mycobionts and their associated Nostoc photobionts within a global framework. The study is based on Internal Transcribed Spacer (ITS) sequences of fungal symbionts and tRNA^Leu (UAA) intron sequences of cyanobacterial symbionts. The full data set includes 271 Nephroma and 358 Nostoc sequences, with over 150 sequence pairs known to originate from the same lichen thalli. Our results show that all bipartite Nephroma species associate with one group of Nostoc different from Nostoc typically found in tripartite Nephroma species. This conserved association appears to have been inherited from the common ancestor of all extant species. While specific associations between some symbiont genotypes can be observed over vast distances, both symbionts tend to show genetic differentiation over wide geographic scales. Most bipartite Nephroma species share their Nostoc symbionts with one or more other fungal taxa, and no fungal species associates solely with a single Nostoc genotype, supporting the concept of functional lichen guilds. Symbiont selectivity patterns within these lichens are best described as a geographic mosaic, with higher selectivity locally than globally. This may reflect specific habitat preferences of particular symbiont combinations, but also the influence of founder effects.     

Open relationships in the castles of clay: high diversity and low host specificity of Termitomyces fungi associated with fungus-growing termites in Africa

In the African and Asian tropics, termites of the subfamily Macrotermitinae play a major role in the decomposition of dead plant material. Their ecological success lies in the obligate mutualism of the termites with fungi of the genus Termitomyces. Before the advent of molecular studies, the interaction with these fungi was poorly understood. Here, we combined available ITS sequence data from West, Central, and South Africa with data of 39 new samples from East Africa to achieve the most comprehensive view of the diversity and host specificity of Termitomyces symbionts across Africa to date. A high amount of sequence divergence in the ITS sequences was found; 11 different Termitomyces lineages in East Africa and >30 lineages across Africa were identified, and the expected diversity is estimated to be about 41 lineages. The fungal lineages belong to four major clades, each almost exclusively associated with one termite host genus. Analysis of molecular variance revealed that 40% of the ITS sequence variation occurred between host genera, indicating close co-evolution at this level. However, within host genera, fungal lineages and haplotypes were frequently shared among host species and sampling localities, except for fungal symbionts of Odontotermes. Horizontal transmission of fungal symbionts may facilitate the transfer of haplotypes and species among hosts. However, at present, we have little understanding of the maintenance of specificity at the genus level. Possible explanations range from substrate specificity of fungi to an active selection of fungi by termites.     

Diversity of Endosymbiotic Nostoc in Gunnera magellanica (L) from Tierra del Fuego, Chile

Global warming is causing ice retreat in glaciers worldwide, most visibly over the last few decades in some areas of the planet. One of the most affected areas is the region of Tierra del Fuego (southern South America). Vascular plant recolonisation of recently deglaciated areas in this region is initiated by Gunnera magellanica, which forms symbiotic associations with the cyanobacterial genus Nostoc, a trait that likely confers advantages in this colonisation process. This symbiotic association in the genus Gunnera is notable as it represents the only known symbiotic relationship between angiosperms and cyanobacteria. The aim of this work was to study the genetic diversity of the Nostoc symbionts in Gunnera at three different, nested scale levels: specimen, population and region. Three different genomic regions were examined in the study: a fragment of the small subunit ribosomal RNA gene (16S), the RuBisCO large subunit gene coupled with its promoter sequence and a chaperon-like protein (rbcLX) and the ribosomal internal transcribed spacer (ITS) region. The identity of Nostoc as the symbiont was confirmed in all the infected rhizome tissue analysed. Strains isolated in the present study were closely related to strains known to form symbioses with other organisms, such as lichen-forming fungi or bryophytes. We found 12 unique haplotypes in the 16S rRNA (small subunit) region analysis, 19 unique haplotypes in the ITS region analysis and 57 in the RuBisCO proteins region (rbcLX). No genetic variability was found among Nostoc symbionts within a single host plant while Nostoc populations among different host plants within a given sampling site revealed major differences. Noteworthy, interpopulation variation was also shown between recently deglaciated soils and more ancient ones, between eastern and western sites and between northern and southern slopes of Cordillera Darwin. The cell structure of the symbiotic relationship was observed with low-temperature scanning electron microscopy, showing changes in morphology of both cyanobiont cells (differentiate more heterocysts) and plant cells (increased size). Developmental stages of the symbiosis, including cell walls and membranes and EPS matrix states, were also observed.