Strong population genomic structure of the toxic dinoflagellate Alexandrium minutum inferred from meta-transcriptome samples.

Despite theoretical expectations, marine microeukaryote population are often highly structured and the mechanisms behind such patterns remain to be elucidated. These organisms display huge census population sizes, yet genotyping usually requires clonal strains originating from single cells, hindering proper population sampling. Estimating allelic frequency directly from population wide samples, without any isolation step, offers an interesting alternative. Here, we validate the use of meta-transcriptome environmental samples to determine the population genetic structure of the dinoflagellate Alexandrium minutum. Strain and meta-transcriptome based results both indicated a strong genetic structure for A. minutum in Western Europe, to the level expected between cryptic species. The presence of numerous private alleles, and even fixed polymorphism, would indicate ancient divergence and absence of gene flow between populations. Single nucleotide polymorphisms (SNPs) displaying strong allele frequency differences were distributed throughout the genome, which might indicate pervasive selection from standing genetic variation (soft selective sweeps). However, a few genomic regions displayed extremely low diversity that could result from the fixation of adaptive de novo mutations (hard selective sweeps) within the populations.

Differential Influence of Life Cycle on Growth and Toxin Production of three Pseudo-nitzschia Species (Bacillariophyceae).

We used a multistrain approach to study the intra- and interspecific variability of the growth rates of three Pseudo-nitzschia species - P. australis, P. fraudulenta, and P. pungens - and of their domoic acid (DA) production. We carried out mating and batch experiments to investigate the respective effects of strain age and cell size, and thus the influence of their life cycle on the physiology of these species. The cell size - life cycle relationship was characteristic of each species. The influence of age and cell size on the intraspecific variability of growth rates suggests that these characteristics should be considered cautiously for the strains used in physiological studies on Pseudo-nitzschia species. The results from all three species do not support the hypothesis of a decrease in DA production with time since isolation from natural populations. In P. australis, the cellular DA content was rather a function of cell size. More particularly, cells at the gametangia stage of their life cycle contained up to six times more DA than smaller or larger cells incapable of sexual reproduction. These findings reveal a link between P. australis life cycle and cell toxicity. This suggest that life cycle dynamics in Pseudo-nitzschia natural populations may influence bloom toxicity.

Contrasting effects of historical contingency on phenotypic and genomic trajectories during a two-step evolution experiment with bacteria.

BACKGROUND: The impact of historical contingency, i.e. the past evolutionary history of a population, on further adaptation is mostly unknown at both the phenotypic and genomic levels. We addressed this question using a two-step evolution experiment. First, replicate populations of Escherichia coli were propagated in four different environmental conditions for 1000 generations. Then, all replicate populations were transferred and propagated for further 1000 generations to a single new environment. RESULTS: Using this two-step experimental evolution strategy, we investigated, at both the phenotypic and genomic levels, whether and how adaptation in the initial historical environments impacted evolutionary trajectories in a new environment. We showed that both the growth rate and fitness of the evolved populations obtained after the second step of evolution were contingent upon past evolutionary history. In contrast however, the genes that were modified during the second step of evolution were independent from the previous history of the populations. CONCLUSIONS: Our work suggests that historical contingency affects phenotypic adaptation to a new environment. This was however not reflected at the genomic level implying complex relationships between environmental factors and the genotype-to-phenotype map.

Evolutionary processes and cellular functions underlying divergence in Alexandrium minutum.

Understanding divergence in the highly dispersive and seemingly homogeneous pelagic environment for organisms living as free drifters in the water column remains a challenge. Here, we analysed the transcriptome-wide mRNA sequences, as well as the morphology of 18 strains of Alexandrium minutum, a dinoflagellate responsible for harmful algal blooms worldwide, to investigate the functional bases of a divergence event. Analysis of the joint site frequency spectrum (JSFS) pointed towards an ancestral divergence in complete isolations followed by a secondary contact resulting in gene flow between the two diverging groups, but heterogeneous across sites. The sites displaying fixed SNPs were associated with a highly restricted gene flow and a strong overrepresentation of nonsynonymous polymorphism, suggesting the importance of selective pressures as drivers of the divergence. The most divergent transcripts were homologs to genes involved in calcium/potassium fluxes across the membrane, calcium transduction signal and saxitoxin production. The implication of these results in terms of ecological divergence and build-up of reproductive isolation is discussed. Dinoflagellates are especially difficult to study in the field at the ecological level due to their small size and the dynamic nature of their natural environment, but also at the genomic level due to their huge and complex genome and the absence of closely related model organism. This study illustrates the possibility to identify the traits of primary importance in ecology and evolution starting from high-throughput sequencing data, even for such organisms.

Ecological and evolutionary dynamics of coexisting lineages during a long-term experiment with Escherichia coli.

Closely related organisms usually occupy similar ecological niches, leading to intense competition and even extinction. Such competition also can promote rapid phenotypic evolution and ecological divergence. This process may end with the stable occupation of distinct niches or, alternatively, may entail repeated bouts of evolution. Here we examine two Escherichia coli lineages, called L and S, that coexisted for more than 30,000 generations after diverging from a common ancestor. Both lineages underwent sustained phenotypic evolution based on global transcription and resource utilization profiles, with L seeming to encroach over time on the catabolic profile of S. Reciprocal invasion experiments with L and S clones from the same or different generations revealed evolutionary changes in their interaction, including an asymmetry that confirmed the encroachment by L on the niche of the S lineage. In general, L and S clones from the same generation showed negative frequency-dependent effects, consistent with stable coexistence. However, L clones could invade S clones from both earlier and later generations, whereas S clones could invade only L clones from earlier generations. In this system, the long-term coexistence of competing lineages evidently depended on successive rounds of evolution, rather than on initial divergence followed by a static equilibrium.

Evolutionary history and genetic parallelism affect correlated responses to evolution.

We investigated the relationship between genomic and phenotypic evolution among replicate populations of Escherichia coli evolved for 1000 generations in four different environments. By resequencing evolved genomes, we identified parallel changes in genes encoding transcription regulators within and between environments. Depending on both the environment and the altered gene, genetic parallelism at the gene level involved mutations that affected identical codons, protein domains or were widely distributed across the gene. Evolved clones were characterized by parallel phenotypic changes in their respective evolution environments but also in the three alternative environments. Phenotypic parallelism was high for clones that evolved in the same environment, even in the absence of genetic parallelism. By contrast, clones that evolved in different environments revealed a higher parallelism in correlated responses when they shared mutated genes. Altogether, this work shows that after an environmental change or the colonization of a new habitat, similar ecological performance might be expected from individuals that share mutated genes or that experienced similar past selective pressures.

Genetic association of toxin production in the dinoflagellate Alexandrium minutum.

Dinoflagellates of the genus Alexandrium are responsible for harmful algal blooms and produce paralytic shellfish toxins (PSTs). Their very large and complex genomes make it challenging to identify the genes responsible for toxin synthesis. A family-based genomic association study was developed to determine the inheritance of toxin production in Alexandrium minutum and identify genomic regions linked to this production. We show that the ability to produce toxins is inheritable in a Mendelian way, while the heritability of the toxin profile is more complex. We developed the first dinoflagellate genetic linkage map. Using this map, several major results were obtained: 1. A genomic region related to the ability to produce toxins was identified. 2. This region does not contain any polymorphic sxt genes, known to be involved in toxin production in cyanobacteria. 3. The sxt genes, known to be present in a single cluster in cyanobacteria, are scattered on different linkage groups in A. minutum. 4. The expression of two sxt genes not assigned to any linkage group, sxtI and sxtG, may be regulated by the genomic region related to the ability to produce toxins. Our results provide new insights into the organization of toxicity-related genes in A. minutum, suggesting a dissociated genetic mechanism for the production of the different analogues and the ability to produce toxins. However, most of the newly identified genes remain unannotated. This study therefore proposes new candidate genes to be further explored to understand how dinoflagellates synthesize their toxins.

Dinophyceae can use exudates as weapons against the parasite Amoebophrya sp. (Syndiniales).

Parasites in the genus Amoebophrya sp. infest dinoflagellate hosts in marine ecosystems and can be determining factors in the demise of blooms, including toxic red tides. These parasitic protists, however, rarely cause the total collapse of Dinophyceae blooms. Experimental addition of parasite-resistant Dinophyceae (Alexandrium minutum or Scrippsiella donghaienis) or exudates into a well-established host-parasite coculture (Scrippsiella acuminata-Amoebophrya sp.) mitigated parasite success and increased the survival of the sensitive host. This effect was mediated by waterborne molecules without the need for a physical contact. The strength of the parasite defenses varied between dinoflagellate species, and strains of A. minutum and was enhanced with increasing resistant host cell concentrations. The addition of resistant strains or exudates never prevented the parasite transmission entirely. Survival time of Amoebophrya sp. free-living stages (dinospores) decreased in presence of A. minutum but not of S. donghaienis. Parasite progeny drastically decreased with both species. Integrity of the dinospore membrane was altered by A. minutum, providing a first indication on the mode of action of anti-parasitic molecules. These results demonstrate that extracellular defenses can be an effective strategy against parasites that protects not only the resistant cells producing them, but also the surrounding community.

From the sxtA4 Gene to Saxitoxin Production: What Controls the Variability Among Alexandrium minutum and Alexandrium pacificum Strains?

Paralytic shellfish poisoning (PSP) is a human foodborne syndrome caused by the consumption of shellfish that accumulate paralytic shellfish toxins (PSTs, saxitoxin group). In PST-producing dinoflagellates such as Alexandrium spp., toxin synthesis is encoded in the nuclear genome via a gene cluster (sxt). Toxin production is supposedly associated with the presence of a 4th domain in the sxtA gene (sxtA4), one of the core genes of the PST gene cluster. It is postulated that gene expression in dinoflagellates is partially constitutive, with both transcriptional and post-transcriptional processes potentially co-occurring. Therefore, gene structure and expression mode are two important features to explore in order to fully understand toxin production processes in dinoflagellates. In this study, we determined the intracellular toxin contents of twenty European Alexandrium minutum and Alexandrium pacificum strains that we compared with their genome size and sxtA4 gene copy numbers. We observed a significant correlation between the sxtA4 gene copy number and toxin content, as well as a moderate positive correlation between the sxtA4 gene copy number and genome size. The 18 toxic strains had several sxtA4 gene copies (9-187), whereas only one copy was found in the two observed non-toxin producing strains. Exploration of allelic frequencies and expression of sxtA4 mRNA in 11 A. minutum strains showed both a differential expression and specific allelic forms in the non-toxic strains compared with the toxic ones. Also, the toxic strains exhibited a polymorphic sxtA4 mRNA sequence between strains and between gene copies within strains. Finally, our study supported the hypothesis of a genetic determinism of toxin synthesis (i.e., the existence of several genetic isoforms of the sxtA4 gene and their copy numbers), and was also consistent with the hypothesis that constitutive gene expression and moderation by transcriptional and post-transcriptional regulation mechanisms are the cause of the observed variability in the production of toxins by A. minutum.

Epistasis and frequency dependence influence the fitness of an adaptive mutation in a diversifying lineage.

The opportunity for a mutation to invade a population can dramatically vary depending on the context in which this mutation occurs. Such context dependence is difficult to document as it requires the ability to measure how a mutation affects phenotypes and fitness and to manipulate the context in which the mutation occurs. We identified a mutation in a gene encoding a global regulator in one of two ecotypes that diverged from a common ancestor during 1200 generations of experimental evolution. We replaced the ancestral allele by the mutant allele, and vice versa, in several clones isolated during the time course of the evolution experiment, and compared the phenotype and fitness of clones isogenic except for the focal mutation. We show that the fitness and phenotype of the mutation are strongly affected by epistatic interactions between genes in the same genome, as well as by frequency dependent selection resulting from biotic interactions between individuals in the same population. We conclude that amongst the replicate population in which it spread, the mutation we identified is only adaptive when occurring in specific genomes and competing with specific individuals. This study thus demonstrates that the opportunity for an adaptive mutation to spread in an evolutionary lineage can only be understood in the light of its genomic and competitive environments.

Strong constitutive expression divergence among strains but no evidence of differential expression associated with sexual reproduction in Alexandrium minutum.

Sexual reproduction remains poorly characterized in dinoflagellates. This is especially the case at the molecular level. Here crossing experiments were performed among strains of the toxic dinoflagellate Alexandrium minutum belonging to two genetically divergent groups. Gene expression was compared between sexually compatible and incompatible crosses at the time of gamete fusion and resting cyst (~zygote) formation. Not a single transcript was identified as differentially expressed between compatible and incompatible crosses at these two crucial time points of the dinoflagellate life cycle. However, several thousands of transcripts displayed constitutive expression differences between strains. This was especially the case between the strains belonging to the genetically divergent groups. A few hundreds of transcripts were also identified as differentially expressed between strains belonging to opposite mating types. Some of these transcripts displayed homology with the SxtA protein, known to be involved in saxitoxin production in cyanobacteria, as well as with proteins potentially involved in mating in fungi.

Alexandrium pacificum and Alexandrium minutum: Harmful or environmentally friendly?

Alexandrium minutum and Alexandrium pacificum are representatives of the dinoflagellate genus that regularly proliferate on the French coasts and other global coastlines. These harmful species may threaten shellfish harvest and human health due to their ability to synthesize neurotoxic alkaloids of the saxitoxin group. However, some dinoflagellates such as A. minutum, and as reported here A. pacificum as well, may also have a beneficial impact on the environment by producing dimethylsulfoniopropionate-DMSP, the precursor of dimethylsulfur-DMS and sulfate aerosols involved in climate balance. However, environmental conditions might influence Alexandrium physiology towards the production of harmful or environmentally friendly compounds. After assessing the influence of two salinity regimes (33 and 38) relative to each species origin (Atlantic French coast and Mediterranean Lagoon respectively), it appears that DMSP and toxin content was variable between the three experimented strains and that higher salinity disadvantages toxin production and tends to favor the production of the osmolytes DMSP and glycine betaine. Hence, this key metabolite production is strain and species-dependent and is influenced by environmental conditions of salinity which in turn, can diversely affect the environment. Widespread coastal blooms of A. minutum and A. pacificum, although being a risk for seafood contamination with toxins, are also a DMSP and DMS source that potentially contribute to the ecosystem structuration and climate. Regarding recent advances in DMSP biosynthesis pathway, 3 dsyB homologs were found in A. minutum but no homolog of the diatom sequence TpMMT.

Species specific gene expression dynamics during harmful algal blooms.

Harmful algal blooms are caused by specific members of microbial communities. Understanding the dynamics of these events requires comparing the strategies developed by the problematic species to cope with environmental fluctuations to the ones developed by the other members of the community. During three consecutive years, the meta-transcriptome of micro-eukaryote communities was sequenced during blooms of the toxic dinoflagellate Alexandrium minutum. The dataset was analyzed to investigate species specific gene expression dynamics. Major shifts in gene expression were explained by the succession of different species within the community. Although expression patterns were strongly correlated with fluctuation of the abiotic environment, and more specifically with nutrient concentration, transcripts specifically involved in nutrient uptake and metabolism did not display extensive changes in gene expression. Compared to the other members of the community, A. minutum displayed a very specific expression pattern, with lower expression of photosynthesis transcripts and central metabolism genes (TCA cycle, glucose metabolism, glycolysis…) and contrasting expression pattern of ion transporters across environmental conditions. These results suggest the importance of mixotrophy, cell motility and cell-to-cell interactions during A. minutum blooms.

Metabolic changes associated with adaptive diversification in Escherichia coli.

During a 1000-generation evolution experiment, two types of morphologically and kinetically distinct bacteria repeatedly diverged from a common ancestor in a fully sympatric seasonal environment containing glucose and acetate. To investigate the metabolic modifications associated with this adaptive diversification, we compared transcription profiles of the two derived types and the common ancestor. Both derived types share a suite of common metabolic changes that may represent adaptation to the environment preceding the diversification event. These include improved translation efficiency, glucose uptake capacity via the mal/lamB genes, upregulation of various transporters during stationary phase, and likely the disruption of the rbs operon. The diversification event is associated with the overexpression of genes involved in the TCA cycle, glyoxylate shunt, acetate consumption, and anaerobic respiration in one type and in acetate excretion in the other. These results reveal that competition for both carbon and oxygen have likely played an important role in the adaptation of Escherichia coli during this adaptive diversification event, where one derived type mainly consumes glucose at a fast rate when oxygen is not limiting, and the other derived type consumes glucose and acetate at a slower rate, even when oxygen is limiting.

Inter- and Intra-Specific Transcriptional and Phenotypic Responses of Pseudo-nitzschia under Different Nutrient Conditions.

Untangling the functional basis of divergence between closely related species is a step toward understanding species dynamics within communities at both the evolutionary and ecological scales. We investigated cellular (i.e., growth, domoic acid production, and nutrient consumption) and molecular (transcriptomic analyses) responses to varying nutrient concentrations across several strains belonging to three species of the toxic diatom genus Pseudo-nitzschia. Three main results were obtained. First, strains from the same species displayed similar transcriptomic, but not necessarily cellular, responses to the experimental conditions. It showed the importance of considering intraspecific diversity to investigate functional divergence between species. Second, a major exception to the first finding was a strain recently isolated from the natural environment and displaying contrasting gene expression patterns related to cell motility and domoic acid production. This result illustrated the profound modifications that may occur when transferring a cell from the natural to the in vitro environment and asks for future studies to better understand the influence of culture duration and life cycle on expression patterns. Third, transcriptomic responses were more similar between the two species displaying similar ecology in situ, irrespective of the genetic distance. This was especially true for molecular responses related to TCA cycle, photosynthesis, and nitrogen metabolism. However, transcripts related to phosphate uptake were variable between species. It highlighted the importance of considering both overall genetic distance and ecological divergence to explain functional divergence between species.

Cophylogeny of the anther smut fungi and their caryophyllaceous hosts: prevalence of host shifts and importance of delimiting parasite species for inferring cospeciation.

BACKGROUND: Using phylogenetic approaches, the expectation that parallel cladogenesis should occur between parasites and hosts has been validated in some studies, but most others provided evidence for frequent host shifts. Here we examine the evolutionary history of the association between Microbotryum fungi that cause anther smut disease and their Caryophyllaceous hosts. We investigated the congruence between host and parasite phylogenies, inferred cospeciation events and host shifts, and assessed whether geography or plant ecology could have facilitated the putative host shifts identified. For cophylogeny analyses on microorganisms, parasite strains isolated from different host species are generally considered to represent independent evolutionary lineages, often without checking whether some strains actually belong to the same generalist species. Such an approach may mistake intraspecific nodes for speciation events and thus bias the results of cophylogeny analyses if generalist species are found on closely related hosts. A second aim of this study was therefore to evaluate the impact of species delimitation on the inferences of cospeciation. RESULTS: We inferred a multiple gene phylogeny of anther smut strains from 21 host plants from several geographic origins, complementing a previous study on the delimitation of fungal species and their host specificities. We also inferred a multi-gene phylogeny of their host plants, and the two phylogenies were compared. A significant level of cospeciation was found when each host species was considered to harbour a specific parasite strain, i.e. when generalist parasite species were not recognized as such. This approach overestimated the frequency of cocladogenesis because individual parasite species capable of infecting multiple host species (i.e. generalists) were found on closely related hosts. When generalist parasite species were appropriately delimited and only a single representative of each species was retained, cospeciation events were not more frequent than expected under a random distribution, and many host shifts were inferred.Current geographic distributions of host species seemed to be of little relevance for understanding the putative historical host shifts, because most fungal species had overlapping geographic ranges. We did detect some ecological similarities, including shared pollinators and habitat types, between host species that were diseased by closely related anther smut species. Overall, genetic similarity underlying the host-parasite interactions appeared to have the most important influence on specialization and host-shifts: generalist multi-host parasite species were found on closely related plant species, and related species in the Microbotryum phylogeny were associated with members of the same host clade. CONCLUSION: We showed here that Microbotryum species have evolved through frequent host shifts to moderately distant hosts, and we show further that accurate delimitation of parasite species is essential for interpreting cophylogeny studies.

Speciation in fungi.

In this review on fungal speciation, we first contrast the issues of species definition and species criteria and show that by distinguishing the two concepts the approaches to studying the speciation can be clarified. We then review recent developments in the understanding of modes of speciation in fungi. Allopatric speciation raises no theoretical problem and numerous fungal examples exist from nature. We explain the theoretical difficulties raised by sympatric speciation, review the most recent models, and provide some natural examples consistent with speciation in sympatry. We describe the nature of prezygotic and postzygotic reproductive isolation in fungi and examine their evolution as functions of temporal and of the geographical distributions. We then review the theory and evidence for roles of cospeciation, host shifts, hybridization, karyotypic rearrangement, and epigenetic mechanisms in fungal speciation. Finally, we review the available data on the genetics of speciation in fungi and address the issue of speciation in asexual species.

Evolution of reproductive isolation within a parasitic fungal species complex.

Despite important advances in the last few years, the evolution of reproductive isolation (RI) remains an unresolved and critical gap in our understanding of speciation processes. In this study, we investigated the evolution of RI among species of the parasitic fungal species complex Microbotryum violaceum, which is responsible for anther smut disease of the Caryophyllaceae. We found no evidence for significant positive assortative mating by M. violaceum even over substantial degrees of genetic divergence, suggesting a lack of prezygotic isolation. In contrast, postzygotic isolation increased with the genetic distance between mating partners when measured as hyphal growth. Total RI, measured as the ability of the pathogen to infect and produce a diploid progeny in the host plant, was significantly and positively correlated with genetic distance, remaining below complete isolation for most of the species pairs. The results of this study, the first one on the time course of speciation in a fungus, are therefore consistent with previous works showing that RI generally evolves gradually with genetic distance, and thus presumably with time. Interestingly, prezygotic RI due to gamete recognition did not increase with genetic distance, in contrast to the pattern found in plants and animals.

Phylogenetic evidence of host-specific cryptic species in the anther smut fungus.

Cryptic structure of species complexes confounds an accurate accounting of biological diversity in natural systems. Also, cryptic sibling species often become specialized to different ecological conditions, for instance, with host specialization by cryptic parasite species. The fungus Microbotryum violaceum causes anther smut disease in plants of Caryophyllaceae, and the degree of specialization and gene flow between strains on different hosts have been controversial in the literature. We conducted molecular phylogenetic analyses on M. violaceum from 23 host species and different geographic origins using three single-copy nuclear genes (beta-tub, gamma-tub, and Ef1alpha). Congruence between the phylogenies identified several lineages that evolved independently for a long time. The lineages had overlapping geographic ranges but were highly specialized on different hosts. These results thus suggest that M. violaceum is a complex of highly specialized sibling species. Two incongruencies between the individual gene phylogenies and one intragene recombination event were detected at basal nodes, suggesting ancient introgression events or speciation events via hybridizations. However, incongruencies and recombination were not detected among terminal branches, indicating that the potentials for cross-infection and experimental hybridization are often not sufficient to suggest that introgressions would likely persist in nature.

Inter and intra-specific growth and domoic acid production in relation to nutrient ratios and concentrations in Pseudo-nitzschia: phosphate an important factor.

The factors responsible for inducing the synthesis of toxins and responses from toxic phytoplankton blooms remain unclear. In this study we compare the influence of genotypic (at both the intra and interspecific levels) and environmental factors (nutrient concentration and ratio) on growth (in terms of cell densities) and domoic acid (DA) production in three Pseudo-nitzschia species: P. australis, P.pungens and P.fradulenta. A strong phosphate effect was detected. More precisely, a low initial concentration in phosphate, even at high initial nitrogen and silicate concentrations, induced the highest DA concentrations and the lowest cell densities in all strains/species studied. In contrast, a low initial concentration of nitrogen and silicate combined, with a higher phosphate concentration resulted in low cell densities, but without high DA production. Inter-species effects were also observed in DA production, where P. australis represented the most toxigenic species of all. Intra-specific variations were only moderate, except for a recently isolated P. australis strain, suggesting the influence of time since isolation on the physiology and DA production of Pseudo-nitzschia species. Overall, the lack of strong interaction between environmental and genotypic factors showed that the various genotypes investigated did not extensively diverge in their ability to respond (in terms of DA production and cell densities) to contrasting nutrient supply.