Heterologous Expression of the Phytochelatin Synthase CaPCS2 from Chlamydomonas acidophila and Its Effect on Different Stress Factors in Escherichia coli.

Phytochelatins (PCs) are cysteine-rich small peptides, enzymatically synthesized from reduced glutathione (GSH) by cytosolic enzyme phytochelatin synthase (PCS). The open reading frame (ORF) of the phytochelatin synthase CaPCS2 gene from the microalgae Chlamydomonas acidophila was heterologously expressed in Escherichia coli strain DH5α, to analyze its role in protection against various abiotic agents that cause cellular stress. The transformed E. coli strain showed increased tolerance to exposure to different heavy metals (HMs) and arsenic (As), as well as to acidic pH and exposure to UVB, salt, or perchlorate. In addition to metal detoxification activity, new functions have also been reported for PCS and PCs. According to the results obtained in this work, the heterologous expression of CaPCS2 in E. coli provides protection against oxidative stress produced by metals and exposure to different ROS-inducing agents. However, the function of this PCS is not related to HM bioaccumulation.

Rare and widespread: integrating Bayesian MCMC approaches, Sanger sequencing and Hyb-Seq phylogenomics to reconstruct the origin of the enigmatic Rand Flora genus Camptoloma.

PREMISE: Genera that are widespread, with geographically discontinuous distributions and represented by few species, are intriguing. Is their achieved disjunct distribution recent or ancient in origin? Why are they species-poor? The Rand Flora is a continental-scale pattern in which closely related species appear codistributed in isolated regions over the continental margins of Africa. Genus Camptoloma (Scrophulariaceae) is the most notable example, comprising three species isolated from each other on the northwest, eastern, and southwest Africa. METHODS: We employed Sanger sequencing of nuclear and plastid markers, together with genomic target sequencing of 2190 low-copy nuclear genes, to infer interspecies relationships and the position of Camptoloma within Scrophulariaceae by using supermatrix and multispecies-coalescent approaches. Lineage divergence times and ancestral ranges were inferred with Bayesian Markov chain Monte Carlo (MCMC) approaches. The population history was estimated with phylogeographic coalescent methods. RESULTS: Camptoloma rotundifolium, restricted to Southern Africa, was shown to be a sister species to the disjunct clade formed by C. canariense, endemic to the Canary Islands, and C. lyperiiflorum, distributed in the Horn of Africa-Southern Arabia. Camptoloma was inferred to be sister to the mostly South African tribes Teedieae and Buddlejeae. Stem divergence was dated in the Late Miocene, while the origin of the extant disjunction was inferred as Early Pliocene. CONCLUSIONS: The current disjunct distribution of Camptoloma across Africa was likely the result of fragmentation and extinction and/or population bottlenecking events associated with historical aridification cycles during the Neogene; the pattern of species divergence, from south to north, is consistent with the "climatic refugia" Rand Flora hypothesis.

Gene-rich UV sex chromosomes harbor conserved regulators of sexual development.

Nonrecombining sex chromosomes, like the mammalian Y, often lose genes and accumulate transposable elements, a process termed degeneration. The correlation between suppressed recombination and degeneration is clear in animal XY systems, but the absence of recombination is confounded with other asymmetries between the X and Y. In contrast, UV sex chromosomes, like those found in bryophytes, experience symmetrical population genetic conditions. Here, we generate nearly gapless female and male chromosome-scale reference genomes of the moss Ceratodon purpureus to test for degeneration in the bryophyte UV sex chromosomes. We show that the moss sex chromosomes evolved over 300 million years ago and expanded via two chromosomal fusions. Although the sex chromosomes exhibit weaker purifying selection than autosomes, we find that suppressed recombination alone is insufficient to drive degeneration. Instead, the U and V sex chromosomes harbor thousands of broadly expressed genes, including numerous key regulators of sexual development across land plants.

Evolutionary history of the mediterranean Pinus halepensis-brutia species complex using gene-resequencing and transcriptomic approaches.

KEY MESSAGE: Complementary gene-resequencing and transcriptomic approaches reveal contrasted evolutionary histories in a species complex. Pinus halepensis and Pinus brutia are closely related species that can intercross, but occupy different geographical ranges and bioclimates. To study the evolution of this species complex and to provide genomic resources for further research, we produce and analyze two new complementary sets of genetic resources: (i) a set of 172 re-sequenced genomic target loci analyzed in 45 individuals, and (ii) a set of 11 transcriptome assemblies. These two datasets provide insights congruent with previous studies: P. brutia displays high level of genetic diversity and no genetic sub-structure, while P. halepensis shows three main genetic clusters, the western Mediterranean and North African clusters displaying much lower genetic diversity than the eastern Mediterranean cluster, the latter cluster having similar genetic diversity to P. brutia. In addition, these datasets provide new insights on the timing of the species-complex history: the two species would have split at the end of the tertiary, and the changing climatic conditions of the Mediterranean region at the end of the Tertiary-beginning of the Quaternary, together with the distinct species tolerance to harsh climatic conditions would have resulted in different geographic distributions, demographic histories and genetic patterns of the two pines. The multiple glacial-interglacial cycles during the Quaternary would have led to the expansion of P. brutia in the Middle East, while P. halepensis would have been through bottlenecks. The last glaciations, from 0.6 Mya on, would have affected further the Western genetic pool of P. halepensis.

Toxicity, Physiological, and Ultrastructural Effects of Arsenic and Cadmium on the Extremophilic Microalga Chlamydomonas acidophila.

The cytotoxicity of cadmium (Cd), arsenate (As(V)), and arsenite (As(III)) on a strain of Chlamydomonas acidophila, isolated from the Rio Tinto, an acidic environment containing high metal(l)oid concentrations, was analyzed. We used a broad array of methods to produce complementary information: cell viability and reactive oxygen species (ROS) generation measures, ultrastructural observations, transmission electron microscopy energy dispersive x-ray microanalysis (TEM-XEDS), and gene expression. This acidophilic microorganism was affected differently by the tested metal/metalloid: It showed high resistance to arsenic while Cd was the most toxic heavy metal, showing an LC50 = 1.94 µM. Arsenite was almost four-fold more toxic (LC50= 10.91 mM) than arsenate (LC50 = 41.63 mM). Assessment of ROS generation indicated that both arsenic oxidation states generate superoxide anions. Ultrastructural analysis of exposed cells revealed that stigma, chloroplast, nucleus, and mitochondria were the main toxicity targets. Intense vacuolization and accumulation of energy reserves (starch deposits and lipid droplets) were observed after treatments. Electron-dense intracellular nanoparticle-like formation appeared in two cellular locations: inside cytoplasmic vacuoles and entrapped into the capsule, around each cell. The chemical nature (Cd or As) of these intracellular deposits was confirmed by TEM-XEDS. Additionally, they also contained an unexpected high content in phosphorous, which might support an essential role of poly-phosphates in metal resistance.

Bridging the micro- and macroevolutionary levels in phylogenomics: Hyb-Seq solves relationships from populations to species and above.

Reconstructing phylogenetic relationships at the micro- and macroevoutionary levels within the same tree is problematic because of the need to use different data types and analytical frameworks. We test the power of target enrichment to provide phylogenetic resolution based on DNA sequences from above species to within populations, using a large herbarium sampling and Euphorbia balsamifera (Euphorbiaceae) as a case study. Target enrichment with custom probes was combined with genome skimming (Hyb-Seq) to sequence 431 low-copy nuclear genes and partial plastome DNA. We used supermatrix, multispecies-coalescent approaches, and Bayesian dating to estimate phylogenetic relationships and divergence times. Euphorbia balsamifera, with a disjunct Rand Flora-type distribution at opposite sides of Africa, comprises three well-supported subspecies: western Sahelian sepium is sister to eastern African-southern Arabian adenensis and Macaronesian-southwest Moroccan balsamifera. Lineage divergence times support Late Miocene to Pleistocene diversification and climate-driven vicariance to explain the Rand Flora pattern. We show that probes designed using genomic resources from taxa not directly related to the focal group are effective in providing phylogenetic resolution at deep and shallow evolutionary levels. Low capture efficiency in herbarium samples increased the proportion of missing data but did not bias estimation of phylogenetic relationships or branch lengths.

De novo assembly of English yew (Taxus baccata) transcriptome and its applications for intra- and inter-specific analyses.

KEY MESSAGE: We provide novel genomic resources for Taxus baccata in the form of a reference transcriptome, SSR and SNP markers, and orthologous single-copy genes, useful for phylogenomic and population genomic applications. English yew (T. baccata) is the only European representative of the Taxaceae family, a conifer group originated in the Jurassic period. The wide extent of environmental heterogeneity within the species' range, together with its long presence in Europe, make English yew an ideal species to investigate adaptive evolution in conifers. To enlarge the genomic resources available for this species, we used Illumina short read sequencing followed by de novo assembly to build the transcriptome of English yew. In addition to a fully annotated transcriptome as well as large sets of new potential SSR and SNP markers for T. baccata, we provide a data set of orthologous single-copy genes across three Taxus species using Picea sitchensis as outgroup, and discuss ortholog uses and limitations for phylogenomic and population genomic applications.

Basis of genetic adaptation to heavy metal stress in the acidophilic green alga Chlamydomonas acidophila.

To better understand heavy metal tolerance in Chlamydomonas acidophila, an extremophilic green alga, we assembled its transcriptome and measured transcriptomic expression before and after Cd exposure in this and the neutrophilic model microalga Chlamydomonas reinhardtii. Genes possibly related to heavy metal tolerance and detoxification were identified and analyzed as potential key innovations that enable this species to live in an extremely acid habitat with high levels of heavy metals. In addition we provide a data set of single orthologous genes from eight green algal species as a valuable resource for comparative studies including eukaryotic extremophiles. Our results based on differential gene expression, detection of unique genes and analyses of codon usage all indicate that there are important genetic differences in C. acidophila compared to C. reinhardtii. Several efflux family proteins were identified as candidate key genes for adaptation to acid environments. This study suggests for the first time that exposure to cadmium strongly increases transposon expression in green algae, and that oil biosynthesis genes are induced in Chlamydomonas under heavy metal stress. Finally, the comparison of the transcriptomes of several acidophilic and non-acidophilic algae showed that the Chlamydomonas genus is polyphyletic and that acidophilic algae have distinctive aminoacid usage patterns.

Altitudinal gradients, biogeographic history and microhabitat adaptation affect fine-scale spatial genetic structure in African and Neotropical populations of an ancient tropical tree species.

The analysis of fine-scale spatial genetic structure (FSGS) within populations can provide insights into eco-evolutionary processes. Restricted dispersal and locally occurring genetic drift are the primary causes for FSGS at equilibrium, as described in the isolation by distance (IBD) model. Beyond IBD expectations, spatial, environmental or historical factors can affect FSGS. We examined FSGS in seven African and Neotropical populations of the late-successional rain forest tree Symphonia globulifera L. f. (Clusiaceae) to discriminate the influence of drift-dispersal vs. landscape/ecological features and historical processes on FSGS. We used spatial principal component analysis and Bayesian clustering to assess spatial genetic heterogeneity at SSRs and examined its association with plastid DNA and habitat features. African populations (from Cameroon and São Tomé) displayed a stronger FSGS than Neotropical populations at both marker types (mean Sp = 0.025 vs. Sp = 0.008 at SSRs) and had a stronger spatial genetic heterogeneity. All three African populations occurred in pronounced altitudinal gradients, possibly restricting animal-mediated seed dispersal. Cyto-nuclear disequilibria in Cameroonian populations also suggested a legacy of biogeographic history to explain these genetic patterns. Conversely, Neotropical populations exhibited a weaker FSGS, which may reflect more efficient wide-ranging seed dispersal by Neotropical bats and other dispersers. The population from French Guiana displayed an association of plastid haplotypes with two morphotypes characterized by differential habitat preferences. Our results highlight the importance of the microenvironment for eco-evolutionary processes within persistent tropical tree populations.

Development of genomic tools in a widespread tropical tree, Symphonia globulifera L.f.: a new low-coverage draft genome, SNP and SSR markers.

Population genetic studies in tropical plants are often challenging because of limited information on taxonomy, phylogenetic relationships and distribution ranges, scarce genomic information and logistic challenges in sampling. We describe a strategy to develop robust and widely applicable genetic markers based on a modest development of genomic resources in the ancient tropical tree species Symphonia globulifera L.f. (Clusiaceae), a keystone species in African and Neotropical rainforests. We provide the first low-coverage (11X) fragmented draft genome sequenced on an individual from Cameroon, covering 1.027 Gbp or 67.5% of the estimated genome size. Annotation of 565 scaffolds (7.57 Mbp) resulted in the prediction of 1046 putative genes (231 of them containing a complete open reading frame) and 1523 exact simple sequence repeats (SSRs, microsatellites). Aligning a published transcriptome of a French Guiana population against this draft genome produced 923 high-quality single nucleotide polymorphisms. We also preselected genic SSRs in silico that were conserved and polymorphic across a wide geographical range, thus reducing marker development tests on rare DNA samples. Of 23 SSRs tested, 19 amplified and 18 were successfully genotyped in four S. globulifera populations from South America (Brazil and French Guiana) and Africa (Cameroon and São Tomé island, FST  = 0.34). Most loci showed only population-specific deviations from Hardy-Weinberg proportions, pointing to local population effects (e.g. null alleles). The described genomic resources are valuable for evolutionary studies in Symphonia and for comparative studies in plants. The methods are especially interesting for widespread tropical or endangered taxa with limited DNA availability.

Horizontal Gene Transfer of Phytochelatin Synthases from Bacteria to Extremophilic Green Algae.

Transcriptomic sequencing together with bioinformatic analyses and an automated annotation process led us to identify novel phytochelatin synthase (PCS) genes from two extremophilic green algae (Chlamydomonas acidophila and Dunaliella acidophila). These genes are of intermediate length compared to known PCS genes from eukaryotes and PCS-like genes from prokaryotes. A detailed phylogenetic analysis gives new insight into the complicated evolutionary history of PCS genes and provides evidence for multiple horizontal gene transfer events from bacteria to eukaryotes within the gene family. A separate subgroup containing PCS-like genes within the PCS gene family is not supported since the PCS genes are monophyletic only when the PCS-like genes are included. The presence and functionality of the novel genes in the organisms were verified by genomic sequencing and qRT-PCR. Furthermore, the novel PCS gene in Chlamydomonas acidophila showed very strong induction by cadmium. Cloning and expression of the gene in Escherichia coli clearly improves its cadmium resistance. The gene in Dunaliella was not induced, most likely due to gene duplication.

Comparative Transcriptomic Analysis of the Response of Dunaliella acidophila (Chlorophyta) to Short-Term Cadmium and Chronic Natural Metal-Rich Water Exposures.

Heavy metals are toxic compounds known to cause multiple and severe cellular damage. However, acidophilic extremophiles are able to cope with very high concentrations of heavy metals. This study investigated the stress response under natural environmental heavy metal concentrations in an acidophilic Dunaliella acidophila. We employed Illumina sequencing for a de novo transcriptome assembly and to identify changes in response to high cadmium concentrations and natural metal-rich water. The photosynthetic performance was also estimated by pulse amplitude-modulated (PAM) fluorescence. Transcriptomic analysis highlights a number of processes mainly related to a high constitutive expression of genes involved in oxidative stress and response to reactive oxygen species (ROS), even in the absence of heavy metals. Photosynthetic activity seems to be unaltered under short-term exposition to Cd and chronic exposure to natural metal-rich water, probably due to an increase in the synthesis of structural photosynthetic components preserving their functional integrity. An overrepresentation of Gene Ontology (GO) terms related to metabolic activities, transcription, and proteosomal catabolic process was observed when D. acidophila grew under chronic exposure to natural metal-rich water. GO terms involved in carbohydrate metabolic process, reticulum endoplasmic and Golgi bodies, were also specifically overrepresented in natural metal-rich water library suggesting an endoplasmic reticulum stress response.

Solar Radiation Stress in Natural Acidophilic Biofilms of Euglena mutabilis Revealed by Metatranscriptomics and PAM Fluorometry.

The daily photosynthetic performance of a natural biofilm of the extreme acidophilic Euglena mutabilis from Río Tinto (SW, Spain) under full solar radiation was analyzed by means of pulse amplitude-modulated (PAM) fluorescence measurements and metatrascriptomic analysis. Natural E. mutabilis biofilms undergo large-scale transcriptomic reprogramming during midday due to a dynamic photoinhibition and solar radiation stress. Photoinhibition is due to UV radiation and not to light intensity, as revealed by PAM fluorometry analysis. In order to minimize the negative effects of solar radiation, our data supports the presence of a circadian rhythm in this euglenophyte that increases their opportunity to survive. Differential gene expression throughout the day (at 12:00, 20:00 and night) was monitored by massive Illumina parallel sequencing of metatranscriptomic libraries. The transcription pattern was altered in genes involved in Photosystem II stability and repair, UV damaged DNA repair, non-photochemical quenching and oxidative stress, supporting the photoinhibition detected by PAM fluorometry at midday.

Evolution of the tRNALeu (UAA) Intron and Congruence of Genetic Markers in Lichen-Symbiotic Nostoc.

The group I intron interrupting the tRNALeu UAA gene (trnL) is present in most cyanobacterial genomes as well as in the plastids of many eukaryotic algae and all green plants. In lichen symbiotic Nostoc, the P6b stem-loop of trnL intron always involves one of two different repeat motifs, either Class I or Class II, both with unresolved evolutionary histories. Here we attempt to resolve the complex evolution of the two different trnL P6b region types. Our analysis indicates that the Class II repeat motif most likely appeared first and that independent and unidirectional shifts to the Class I motif have since taken place repeatedly. In addition, we compare our results with those obtained with other genetic markers and find strong evidence of recombination in the 16S rRNA gene, a marker widely used in phylogenetic studies on Bacteria. The congruence of the different genetic markers is successfully evaluated with the recently published software Saguaro, which has not previously been utilized in comparable studies.

Transcriptional response to copper excess and identification of genes involved in heavy metal tolerance in the extremophilic microalga Chlamydomonas acidophila.

High concentrations of heavy metals are typical of acidic environments. Therefore, studies on acidophilic organisms in their natural environments improve our understanding on the evolution of heavy metal tolerance and detoxification in plants. Here we sequenced the transcriptome of the extremophilic microalga Chlamydomonas acidophila cultivated in control conditions and with 500 µM of copper for 24 h. High-throughput 454 sequencing was followed by de novo transcriptome assembly. The reference transcriptome was annotated and genes related to heavy metal tolerance and abiotic stress were identified. Analyses of differentially expressed transcripts were used to detect genes involved in metabolic pathways related to abiotic stress tolerance, focusing on effects caused by increased levels of copper. Both transcriptomic data and observations from PAM fluorometry analysis suggested that the photosynthetic activity of C. acidophila is not adversely affected by addition of high amounts of copper. Up-regulated transcripts include several transcripts related to photosynthesis and carbohydrate metabolism, transcripts coding for general stress response, and a transcript annotated as homologous to the oil-body-associated protein HOGP coding gene. The first de novo assembly of C. acidophila significantly increases transcriptomic data available on extremophiles and green algae and thus provides an important reference for further molecular genetic studies. The differences between differentially expressed transcripts detected in our study suggest that the response to heavy metal exposure in C. acidophila is different from other studied green algae.

WRKY6 transcription factor restricts arsenate uptake and transposon activation in Arabidopsis.

Stress constantly challenges plant adaptation to the environment. Of all stress types, arsenic was a major threat during the early evolution of plants. The most prevalent chemical form of arsenic is arsenate, whose similarity to phosphate renders it easily incorporated into cells via the phosphate transporters. Here, we found that arsenate stress provokes a notable transposon burst in plants, in coordination with arsenate/phosphate transporter repression, which immediately restricts arsenate uptake. This repression was accompanied by delocalization of the phosphate transporter from the plasma membrane. When arsenate was removed, the system rapidly restored transcriptional expression and membrane localization of the transporter. We identify WRKY6 as an arsenate-responsive transcription factor that mediates arsenate/phosphate transporter gene expression and restricts arsenate-induced transposon activation. Plants therefore have a dual WRKY-dependent signaling mechanism that modulates arsenate uptake and transposon expression, providing a coordinated strategy for arsenate tolerance and transposon gene silencing.

Phylogeny-based comparative methods question the adaptive nature of sporophytic specializations in mosses.

Adaptive evolution has often been proposed to explain correlations between habitats and certain phenotypes. In mosses, a high frequency of species with specialized sporophytic traits in exposed or epiphytic habitats was, already 100 years ago, suggested as due to adaptation. We tested this hypothesis by contrasting phylogenetic and morphological data from two moss families, Neckeraceae and Lembophyllaceae, both of which show parallel shifts to a specialized morphology and to exposed epiphytic or epilithic habitats. Phylogeny-based tests for correlated evolution revealed that evolution of four sporophytic traits is correlated with a habitat shift. For three of them, evolutionary rates of dual character-state changes suggest that habitat shifts appear prior to changes in morphology. This suggests that they could have evolved as adaptations to new habitats. Regarding the fourth correlated trait the specialized morphology had already evolved before the habitat shift. In addition, several other specialized "epiphytic" traits show no correlation with a habitat shift. Besides adaptive diversification, other processes thus also affect the match between phenotype and environment. Several potential factors such as complex genetic and developmental pathways yielding the same phenotypes, differences in strength of selection, or constraints in phenotypic evolution may lead to an inability of phylogeny-based comparative methods to detect potential adaptations.

Reconstruction of structural evolution in the trnL intron P6b loop of symbiotic Nostoc (Cyanobacteria).

In this study we reconstruct the structural evolution of the hyper-variable P6b region of the group I trnLeu intron in a monophyletic group of lichen-symbiotic Nostoc strains and establish it as a useful marker in the phylogenetic analysis of these organisms. The studied cyanobacteria occur as photosynthetic and/or nitrogen-fixing symbionts in lichen species of the diverse Nephroma guild. Phylogenetic analyses and secondary structure reconstructions are used to improve the understanding of the replication mechanisms in the P6b stem-loop and to explain the observed distribution patterns of indels. The variants of the P6b region in the Nostoc clade studied consist of different combinations of five sequence modules. The distribution of indels together with the ancestral character reconstruction performed enables the interpretation of the evolution of each sequence module. Our results indicate that the indel events are usually associated with single nucleotide changes in the P6b region and have occurred several times independently. In spite of their homoplasy, they provide phylogenetic information for closely related taxa. Thus we recognize that features of the P6b region can be used as molecular markers for species identification and phylogenetic studies involving symbiotic Nostoc cyanobacteria.