Whole transcriptome RNA-Seq analysis reveals extensive cell type-specific compartmentalization in Volvox carteri

Different transcription of novel, functional long non-coding RNA genes by UV-B in green algae, Volvox carteri

Long non-coding RNAs (lncRNAs) are identified as important regulatory molecules related to diverse biological processes. In recent years, benefiting from the rapid development of high-throughput sequencing technology, RNA-seq, and analysis methods, more lncRNAs have been identified and discovered in various plant and algal species. However, so far, only limited studies related to algal lncRNAs are available. Volvox carteri f. nagariensis is the best multicellular model organism to study in developmental and evolutionary biology; therefore, studying and increasing information about this species is important. This study identified lncRNAs in the multicellular green algae Volvox carteri and 1457 lncRNAs were reported, using RNA-seq data and with the help of bioinformatics tools and software. This study investigated the effect of low-dose UV-B radiation on changes in the expression profile of lncRNAs in gonidial and somatic cells. The differential expression of lncRNAs was analyzed between the treatment (UV-B) and the control (WL) groups in gonidial and somatic cells. A total of 37 and 26 lncRNAs with significant differential expression in gonidial and somatic cells, respectively, were reported. Co-expression analysis between the lncRNAs and their neighbor protein-coding genes (in the interval of ± 10 Kb) was accomplished. In gonidial cells, 184 genes with a positive correlation and 13 genes with a negative correlation (greater than 0.95), and in somatic cells, 174 genes with a positive correlation, and 18 genes with a negative correlation were detected. Functional analysis of neighboring coding genes was also performed based on gene ontology. The results of the current work may help gain deeper insight into the regulation of gene expression in the studied model organism, Volvox carteri.

Cell Type-Specific Promoters of Volvox carteri for Molecular Cell Biology Studies

The multicellular green alga Volvox carteri has emerged as a valuable model organism for investigating various aspects of multicellularity and cellular differentiation, photoreception and phototaxis, cell division, biogenesis of the extracellular matrix and morphogenetic movements. While a range of molecular tools and bioinformatics resources have been made available for exploring these topics, the establishment of cell type-specific promoters in V. carteri has not been achieved so far. Therefore, here, we conducted a thorough screening of transcriptome data from RNA sequencing analyses of V. carteri in order to identify potential cell type-specific promoters. Eventually, we chose two putative strong and cell type-specific promoters, with one exhibiting specific expression in reproductive cells (gonidia), the PCY1 promoter, and the other in somatic cells, the PFP promoter. After cloning both promoter regions, they were introduced upstream of a luciferase reporter gene. By using particle bombardment, the DNA constructs were stably integrated into the genome of V. carteri. The results of the expression analyses, which were conducted at both the transcript and protein levels, demonstrated that the two promoters drive cell type-specific expression in their respective target cell types. Transformants with considerably diverse expression levels of the chimeric genes were identifiable. In conclusion, the screening and analysis of transcriptome data from RNA sequencing allowed for the identification of potential cell type-specific promoters in V. carteri. Reporter gene constructs demonstrated the actual usability of two promoters. The investigated PCY1 and PFP promoters were proven to be potent molecular tools for genetic engineering in V. carteri.

The Genetics of Fitness Reorganization during the Transition to Multicellularity: The Volvocine regA-like Family as a Model

The evolutionary transition from single-celled to multicellular individuality requires organismal fitness to shift from the cell level to a cell group. This reorganization of fitness occurs by re-allocating the two components of fitness, survival and reproduction, between two specialized cell types in the multicellular group: soma and germ, respectively. How does the genetic basis for such fitness reorganization evolve? One possible mechanism is the co-option of life history genes present in the unicellular ancestors of a multicellular lineage. For instance, single-celled organisms must regulate their investment in survival and reproduction in response to environmental changes, particularly decreasing reproduction to ensure survival under stress. Such stress response life history genes can provide the genetic basis for the evolution of cellular differentiation in multicellular lineages. The regA-like gene family in the volvocine green algal lineage provides an excellent model system to study how this co-option can occur. We discuss the origin and evolution of the volvocine regA-like gene family, including regA-the gene that controls somatic cell development in the model organism Volvox carteri. We hypothesize that the co-option of life history trade-off genes is a general mechanism involved in the transition to multicellular individuality, making volvocine algae and the regA-like family a useful template for similar investigations in other lineages.

Cell Type-Specific Pherophorins of Volvox carteri Reveal Interplay of Both Cell Types in ECM Biosynthesis

The spheroidal green algae Volvox carteri serves as a model system to investigate the formation of a complex, multifunctional extracellular matrix (ECM) in a relatively simple, multicellular organism with cell differentiation. The V. carteri ECM is mainly composed of hydroxyproline-rich glycoproteins (HRGPs) and there are diverse region-specific, anatomically distinct structures in the ECM. One large protein family with importance for ECM biosynthesis stands out: the pherophorins. The few pherophorins previously extracted from the ECM and characterized, were specifically expressed by somatic cells. However, the localization and function of most pherophorins is unknown. Here, we provide a phylogenetic analysis of 153 pherophorins of V. carteri and its unicellular relative Chlamydomonas reinhardtii. Our analysis of cell type-specific mRNA expression of pherophorins in V. carteri revealed that, contrary to previous assumptions, only about half (52%) of the 102 investigated pherophorin-related genes show stronger expression in somatic cells, whereas about one-third (34%) of the genes show significant higher expression in reproductive cells (gonidia). We fused two pherophorin genes that are expressed by different cell types to yfp, stably expressed them in Volvox and studied the tagged proteins by live-cell imaging. In contrast to earlier biochemical approaches, this genetic approach also allows the in vivo analysis of non-extractable, covalently cross-linked ECM proteins. We demonstrate that the soma-specific pherophorin SSG185 is localized in the outermost ECM structures of the spheroid, the boundary zone and at the flagellar hillocks. SSG185:YFP is detectable as early as 1.5 h after completion of embryogenesis. It is then present for the rest of the life cycle. The gonidia-specific pherophorin PhG is localized in the gonidial cellular zone 1 ("gonidial vesicle") suggesting its involvement in the protection of gonidia and developing embryos until hatching. Even if somatic cells produce the main portion of the ECM of the spheroids, ECM components produced by gonidia are also required to cooperatively assemble the total ECM. Our results provide insights into the evolution of the pherophorin protein family and convey a more detailed picture of Volvox ECM synthesis.

Molecular and cellular dynamics of early embryonic cell divisions in Volvox carteri

Cell division is fundamental to all organisms and the green alga used here exhibits both key animal and plant functions. Specifically, we analyzed the molecular and cellular dynamics of early embryonic divisions of the multicellular green alga Volvox carteri (Chlamydomonadales). Relevant proteins related to mitosis and cytokinesis were identified in silico, the corresponding genes were cloned, fused to yfp, and stably expressed in Volvox, and the tagged proteins were studied by live-cell imaging. We reveal rearrangements of the microtubule cytoskeleton during centrosome separation, spindle formation, establishment of the phycoplast, and generation of previously unknown structures. The centrosomes participate in initiation of spindle formation and determination of spindle orientation. Although the nuclear envelope does not break down during early mitosis, intermixing of cytoplasm and nucleoplasm results in loss of nuclear identity. Finally, we present a model for mitosis in Volvox. Our study reveals enormous dynamics, clarifies spatio-temporal relationships of subcellular structures, and provides insight into the evolution of cell division.

Genome sequencing of the multicellular alga Astrephomene provides insights into convergent evolution of germ-soma differentiation

Germ-soma differentiation evolved independently in many eukaryotic lineages and contributed to complex multicellular organizations. However, the molecular genetic bases of such convergent evolution remain unresolved. Two multicellular volvocine green algae, Volvox and Astrephomene, exhibit convergent evolution of germ-soma differentiation. The complete genome sequence is now available for Volvox, while genome information is scarce for Astrephomene. Here, we generated the de novo whole genome sequence of Astrephomene gubernaculifera and conducted RNA-seq analysis of isolated somatic and reproductive cells. In Volvox, tandem duplication and neofunctionalization of the ancestral transcription factor gene (RLS1/rlsD) might have led to the evolution of regA, the master regulator for Volvox germ-soma differentiation. However, our genome data demonstrated that Astrephomene has not undergone tandem duplication of the RLS1/rlsD homolog or acquisition of a regA-like gene. Our RNA-seq analysis revealed the downregulation of photosynthetic and anabolic gene expression in Astrephomene somatic cells, as in Volvox. Among genes with high expression in somatic cells of Astrephomene, we identified three genes encoding putative transcription factors, which may regulate somatic cell differentiation. Thus, the convergent evolution of germ-soma differentiation in the volvocine algae may have occurred by the acquisition of different regulatory circuits that generate a similar division of labor.

Integrated Gene Co-expression Analysis and Metabolites Profiling Highlight the Important Role of ZmHIR3 in Maize Resistance to Gibberella Stalk Rot

Gibberella stalk rot (GSR) caused by Fusarium graminearum is one of the most devastating diseases causing significant yield loss of maize, and GSR resistance is a quantitative trait controlled by multiple genes. Although a few quantitative trait loci/resistance genes have been identified, the molecular mechanisms underlying GSR resistance remain largely unexplored. To identify potential resistance genes and to better understand the molecular mechanism of GSR resistance, a joint analysis using a comparative transcriptomic and metabolomic approaches was conducted using two inbred lines with contrasting GSR resistance, K09 (resistant) and A08 (susceptible), upon infection with F. graminearum. While a substantial number of differentially expressed genes associated with various defense-related signaling pathways were identified between two lines, multiple hub genes likely associated with GSR resistance were pinpointed using Weighted Gene Correlation Network Analysis and K-means clustering. Moreover, a core set of metabolites, including anthocyanins, associated with the hub genes was determined. Among the complex co-expression networks, ZmHIR3 showed strong correlation with multiple key genes, and genetic and histological studies showed that zmhir3 mutant is more susceptible to GSR, accompanied by enhanced cell death in the stem in response to infection with F. graminearum. Taken together, our study identified differentially expressed key genes and metabolites, as well as co-expression networks associated with distinct infection stages of F. graminearum. Moreover, ZmHIR3 likely plays a positive role in disease resistance to GSR, probably through the transcriptional regulation of key genes, functional metabolites, and the control of cell death.

Enhanced nucleotide analysis enables the quantification of deoxynucleotides in plants and algae revealing connections between nucleoside and deoxynucleoside metabolism

Detecting and quantifying low-abundance (deoxy)ribonucleotides and (deoxy)ribonucleosides in plants remains difficult; this is a major roadblock for the investigation of plant nucleotide (NT) metabolism. Here, we present a method that overcomes this limitation, allowing the detection of all deoxy- and ribonucleotides as well as the corresponding nucleosides from the same plant sample. The method is characterized by high sensitivity and robustness enabling the reproducible detection and absolute quantification of these metabolites even if they are of low abundance. Employing the new method, we analyzed Arabidopsis thaliana null mutants of CYTIDINE DEAMINASE, GUANOSINE DEAMINASE, and NUCLEOSIDE HYDROLASE 1, demonstrating that the deoxyribonucleotide (dNT) metabolism is intricately interwoven with the catabolism of ribonucleosides (rNs). In addition, we discovered a function of rN catabolic enzymes in the degradation of deoxyribonucleosides in vivo. We also determined the concentrations of dNTs in several mono- and dicotyledonous plants, a bryophyte, and three algae, revealing a correlation of GC to AT dNT ratios with genomic GC contents. This suggests a link between the genome and the metabolome previously discussed but not experimentally addressed. Together, these findings demonstrate the potential of this new method to provide insight into plant NT metabolism.

Targeted migration of pherophorin-S indicates extensive extracellular matrix dynamics in Volvox carteri

Hydroxyproline-rich glycoproteins (HRGPs) constitute a major group of proteins of the extracellular matrix (ECM). The multicellular green alga Volvox carteri is a suitable model organism in which to study the evolutionary transition to multicellularity, including the basic principles and characteristics of an ECM. In Volvox, the ECM is dominated by a single HRGP family: the pherophorins. Our inventory amounts to 117 pherophorin-related genes in V. carteri. We focused on a pherophorin with an unexpected characteristic: pherophorin-S is a soluble, non-cross-linked ECM protein. Using transformants expressing a YFP-tagged pherophorin-S we observed the synthesis and secretion of pherophorin-S by somatic cells in vivo, and we then traced the protein during its conspicuous migration to the ECM around prehatching juveniles and its localized concentration there. Our results provide insights into how an ECM zone surrounding the progeny is remotely affected by distantly located parental somatic cells. In view of the properties and migration of pherophorin-S, we conclude that pherophorin-S is likely to act as an ECM plasticizer to allow for dynamic ECM remodeling.

Volvox and volvocine green algae

The transition of life from single cells to more complex multicellular forms has occurred at least two dozen times among eukaryotes and is one of the major evolutionary transitions, but the early steps that enabled multicellular life to evolve and thrive remain poorly understood. Volvocine green algae are a taxonomic group that is uniquely suited to investigating the step-wise acquisition of multicellular organization. The multicellular volvocine species Volvox carteri exhibits many hallmarks of complex multicellularity including complete germ–soma division of labor, asymmetric cell divisions, coordinated tissue-level morphogenesis, and dimorphic sexes—none of which have obvious analogs in its closest unicellular relative, the model alga Chlamydomonas reinhardtii. Here, I summarize some of the key questions and areas of study that are being addressed with Volvox carteri and how increasing genomic information and methodologies for volvocine algae are opening up the entire group as an integrated experimental system for exploring the evolution of multicellularity and more.

Babo1, formerly Vop1 and Cop1/2, is no eyespot photoreceptor but a basal body protein illuminating cell division in Volvox carteri

In photosynthetic organisms many processes are light dependent and sensing of light requires light-sensitive proteins. The supposed eyespot photoreceptor protein Babo1 (formerly Vop1) has previously been classified as an opsin due to the capacity for binding retinal. Here, we analyze Babo1 and provide evidence that it is no opsin. Due to the localization at the basal bodies, the former Vop1 and Cop1/2 proteins were renamed V.c. Babo1 and C.r. Babo1. We reveal a large family of more than 60 Babo1-related proteins from a wide range of species. The detailed subcellular localization of fluorescence-tagged Babo1 shows that it accumulates at the basal apparatus. More precisely, it is located predominantly at the basal bodies and to a lesser extent at the four strands of rootlet microtubules. We trace Babo1 during basal body separation and cell division. Dynamic structural rearrangements of Babo1 particularly occur right before the first cell division. In four-celled embryos Babo1 was exclusively found at the oldest basal bodies of the embryo and on the corresponding d-roots. The unequal distribution of Babo1 in four-celled embryos could be an integral part of a geometrical system in early embryogenesis, which establishes the anterior-posterior polarity and influences the spatial arrangement of all embryonic structures and characteristics. Due to its retinal-binding capacity, Babo1 could also be responsible for the unequal distribution of retinoids, knowing that such concentration gradients of retinoids can be essential for the correct patterning during embryogenesis of more complex organisms. Thus, our findings push the Babo1 research in another direction.

Different regulations of cell-type transcription by UV-B in multicellular green alga Volvox carteri

There is a scarcity of research reports on the effect of ultraviolet (UV)-B radiation on genome-wide transcriptional regulation in the multicellular green microalga including Volvox carteri (V. carteri). This microalga possesses only two cell types including mortal and motile somatic cells, as well as immortal and immotile reproductive cells. Therefore, the present study evaluated the effect of low-dose UV-B radiation on the cell-type-specific gene expression pattern of reproductive and somatic cells in an asexual life cycle of V. carteri using RNA sequence method. To this end, the separated reproductive and somatic cells were treated for 1 hour at an intensity of 0.056 mW/cm^-2 UV-B radiation. Then, a transcriptome analysis was conducted between the UV-B and white light treated groups in either of the cell types. Based on differential gene expression analyses, no differentially expressed genes were found in reproductive cells under the treatment as compared to the control group. This type of cell maintained its steady state. However, treating the somatic cells with UV-B radiation led to at least 126 differentially expressed genes compared to the untreated control group. In addition, the results of a direct comparison demonstrated a restricted and wide response to UV-B radiation in somatic cells as compared to reproductive cells. Based on the results, UV-B radiation could be involved in cell-type-specific regulation of biological pathways.

A halotolerant growth promoting rhizobacteria triggers induced systemic resistance in plants and defends against fungal infection

A halotolerant rhizobacteria, Klebsiella species (referred to MBE02), was identified that had a growth stimulation effect on peanut. To gain mechanistic insights into how molecular components were reprogrammed during the interaction of MBE02 and peanut roots, we performed deep RNA-sequencing. In total, 1260 genes were differentially expressed: 979 genes were up-regulated, whereas 281 were down-regulated by MBE02 treatment as compared to uninoculated controls. A large component of the differentially regulated genes were related to phytohormone signalling. This included activation of a significant proportion of genes involved in jasmonic acid, ethylene and pathogen-defense signalling, which indicated a role of MBE02 in modulating plant immunity. In vivo and in vitro pathogenesis assays demonstrated that MBE02 treatment indeed provide fitness benefits to peanut against Aspergillus infection under controlled as well as field environment. Further, MBE02 directly reduced the growth of a wide range of fungal pathogens including Aspergillus. We also identified possible molecular components involved in rhizobacteria-mediated plant protection. Our results show the potential of MBE02 as a biocontrol agent in preventing infection against several fungal phytopathogens.

CRISPR/Cas9 mutagenesis in Volvox carteri

Volvox carteri and other volvocine green algae comprise an excellent model for investigating developmental complexity and its origins. Here we describe a method for targeted mutagenesis in V. carteri using CRISPR/Cas9 components expressed from transgenes. We used V. carteri nitrate reductase gene (nitA) regulatory sequences to conditionally express Streptococcus pyogenes Cas9, and V. carteri U6 RNA gene regulatory sequences to constitutively express single-guide RNA (sgRNA) transcripts. Volvox carteri was bombarded with both Cas9 vector and one of several sgRNA vectors programmed to target different test genes (glsA, regA and invA), and transformants were selected for expression of a hygromycin-resistance marker present on the sgRNA vector. Hygromycin-resistant transformants grown with nitrate as sole nitrogen source (inducing for nitA) were tested for Cas9 and sgRNA expression, and for the ability to generate progeny with expected mutant phenotypes. Some transformants of a somatic regenerator (Reg) mutant strain receiving sgRNA plasmid with glsA protospacer sequence yielded progeny (at a rate of ~0.01%) with a gonidialess (Gls) phenotype similar to that observed for previously described glsA mutants, and sequencing of the glsA gene in independent mutants revealed short deletions within the targeted region of glsA, indicative of Cas9-directed non-homologous end joining. Similarly, bombardment of a morphologically wild-type strain with the Cas9 plasmid and sgRNA plasmids targeting regA or invA yielded regA and invA mutant transformants/progeny, respectively (at rates of 0.1-100%). The capacity to make precisely directed frameshift mutations should greatly accelerate the molecular genetic analysis of development in V. carteri, and of developmental novelty in the volvocine algae.

Two-component cyclase opsins of green algae are ATP-dependent and light-inhibited guanylyl cyclases

Background

The green algae Chlamydomonas reinhardtii and Volvox carteri are important models for studying light perception and response, expressing many different photoreceptors. More than 10 opsins were reported in C. reinhardtii, yet only two—the channelrhodopsins—were functionally characterized. Characterization of new opsins would help to understand the green algae photobiology and to develop new tools for optogenetics.

Results

Here we report the characterization of a novel opsin family from these green algae: light-inhibited guanylyl cyclases regulated through a two-component-like phosphoryl transfer, called “two-component cyclase opsins” (2c-Cyclops). We prove the existence of such opsins in C. reinhardtii and V. carteri and show that they have cytosolic N- and C-termini, implying an eight-transmembrane helix structure. We also demonstrate that cGMP production is both light-inhibited and ATP-dependent. The cyclase activity of Cr2c-Cyclop1 is kept functional by the ongoing phosphorylation and phosphoryl transfer from the histidine kinase to the response regulator in the dark, proven by mutagenesis. Absorption of a photon inhibits the cyclase activity, most likely by inhibiting the phosphoryl transfer. Overexpression of Vc2c-Cyclop1 protein in V. carteri leads to significantly increased cGMP levels, demonstrating guanylyl cyclase activity of Vc2c-Cyclop1 in vivo. Live cell imaging of YFP-tagged Vc2c-Cyclop1 in V. carteri revealed a development-dependent, layer-like structure at the immediate periphery of the nucleus and intense spots in the cell periphery.

Conclusions

Cr2c-Cyclop1 and Vc2c-Cyclop1 are light-inhibited and ATP-dependent guanylyl cyclases with an unusual eight-transmembrane helix structure of the type I opsin domain which we propose to classify as type Ib, in contrast to the 7 TM type Ia opsins. Overexpression of Vc2c-Cyclop1 protein in V. carteri led to a significant increase of cGMP, demonstrating enzyme functionality in the organism of origin. Fluorescent live cell imaging revealed that Vc2c-Cyclop1 is located in the periphery of the nucleus and in confined areas at the cell periphery.

Electronic supplementary material

The online version of this article (10.1186/s12915-018-0613-5) contains supplementary material, which is available to authorized users.