The organization structure and regulatory elements of Chlamydomonas histone genes reveal features linking plant and animal genes

Characterization of activating cis-regulatory elements from the histone genes of Chlamydomonas reinhardtii

Regulation of gene expression in eukaryotes is controlled by cis-regulatory modules (CRMs). A major class of CRMs are enhancers which are composed of activating cis-regulatory elements (CREs) responsible for upregulating transcription. To date, most enhancers and activating CREs have been studied in angiosperms; in contrast, our knowledge about these key regulators of gene expression in green algae is limited. In this study, we aimed at characterizing putative activating CREs/CRMs from the histone genes of the unicellular model alga Chlamydomonas reinhardtii. To test the activity of four candidates, reporter constructs consisting of a tetramerized CRE, an established promoter, and a gene for the mCerulean3 fluorescent protein were incorporated into the nuclear genome of C. reinhardtii, and their activity was quantified by flow cytometry. Two tested candidates, Eupstr and Ehist cons, significantly upregulated gene expression and were characterized in detail. Eupstr, which originates from highly expressed genes of C. reinhardtii, is an orientation-independent CRE capable of activating both the RBCS2 and β2-tubulin promoters. Ehist cons, which is a CRM from histone genes of angiosperms, upregulates the β2-tubulin promoter in C. reinhardtii over a distance of at least 1.5 kb. The octamer motif present in Ehist cons was identified in C. reinhardtii and the related green algae Chlamydomonas incerta, Chlamydomonas schloesseri, and Edaphochlamys debaryana, demonstrating its high evolutionary conservation. The results of this investigation expand our knowledge about the regulation of gene expression in green algae. Furthermore, the characterized activating CREs/CRMs can be applied as valuable genetic tools.

The landscape of Chlamydomonas histone H3 lysine 4 methylation reveals both constant features and dynamic changes during the diurnal cycle

Chromatin modifications are epigenetic regulatory features with major roles in various cellular events, yet they remain understudied in algae. We interrogated the genome-wide distribution pattern of mono- and trimethylated histone H3 lysine 4 (H3K4) using chromatin-immunoprecipitation followed by deep-sequencing (ChIP-seq) during key phases of the Chlamydomonas cell cycle: early G1 phase, Zeitgeber Time 1 (ZT1), when cells initiate biomass accumulation, S/M phase (ZT13) when cells are replicating DNA and undergoing mitosis, and late G0 phase (ZT23) when they are quiescent. Tri-methylated H3K4 was predominantly enriched at transcription start sites of the majority of protein coding genes (85%). The likelihood of a gene being marked by H3K4me3 correlated with it being transcribed at some point during the life cycle but not necessarily by continuous active transcription, as exemplified by early zygotic genes, which may remain transcriptionally dormant for thousands of generations between sexual cycles. The exceptions to this rule were around 120 loci, some of which encode non-poly-adenylated transcripts, such as small nuclear RNAs and replication-dependent histones that had H3K4me3 peaks only when they were being transcribed. Mono-methylated H3K4 was the default state for the vast majority of histones that were bound outside of transcription start sites and terminator regions of genes. A small fraction of the genome that was depleted of any H3 lysine 4 methylation was enriched for DNA cytosine methylation and the genes within these DNA methylation islands were poorly expressed. Besides marking protein coding genes, H3K4me3 ChIP-seq data served also as a annotation tool for validation of hundreds of long non-coding RNA genes.

Chlamydomonas CHT7 is involved in repressing DNA replication and mitotic genes during synchronous growth

In the green alga Chlamydomonas reinhardtii, regulation of the cell cycle in response to external cues is critical for survival in a changing environment. The loss of the nuclear COMPROMISED HYDROLYSIS OF TRIACYLGLYCEROLS 7 (CHT7) protein affects the expression of many genes especially in response to nitrogen availability. Cells lacking CHT7 exhibit abnormal cell morphology following nitrogen deprivation and fail to resume normal cell division after N resupply. To investigate the function of CHT7 in the regulation of cell cycle-related pathways, cells were synchronized, and RNA-seq analysis was performed during various stages of the cell cycle. In the cht7 mutant following nitrogen deprivation, the cells were not dividing, but a subset of cell cycle genes involved in DNA replication and mitosis were found to be derepressed, suggesting that the CHT7 protein plays a role in cell cycle regulation that is opposite to that of the mitotic cyclin-dependent kinases. Furthermore, genes for cell wall synthesis and remodeling were found to be abnormally induced in nondividing cht7 cells; this misregulation may deplete cellular resources and thus contribute to cell death following nitrogen deprivation. Lastly, 43 minimally characterized kinases were found to be highly misregulated in cht7. Further analysis suggested that some of these CHT7-regulated kinases may be related to the MAP3K and Aurora-like kinases, while others are unique. Together, these results suggest a role of CHT7 in transcriptional regulation of the cell cycle and reveal several pathways and genes whose expression appears to be subject to a CHT7-mediated regulatory network.

A three-dimensional RNA motif mediates directional trafficking of Potato spindle tuber viroid from epidermal to palisade mesophyll cells in Nicotiana benthamiana

Potato spindle tuber viroid (PSTVd) is a circular non-coding RNA of 359 nucleotides that replicates and spreads systemically in host plants, thus all functions required to establish an infection are mediated by sequence and structural elements in the genome. The PSTVd secondary structure contains 26 Watson-Crick base-paired stems and 27 loops. Most of the loops are believed to form three-dimensional (3D) structural motifs through non-Watson-Crick base pairing, base stacking, and other local interactions. Homology-based prediction using the JAR3D online program revealed that loop 27 (nucleotides 177-182) most likely forms a 3D structure similar to the loop of a conserved hairpin located in the 3' untranslated region of histone mRNAs in animal cells. This stem-loop, which is involved in 3'-end maturation, is not found in polyadenylated plant histone mRNAs. Mutagenesis showed that PSTVd genomes containing base substitutions in loop 27 predicted by JAR3D to disrupt the 3D structure were unable to replicate in Nicotiana benthamiana leaves following mechanical rub inoculation, with one exception: a U178G/U179G double mutant was replication-competent and able to spread within the upper epidermis of inoculated leaves, but was confined to this cell layer. Remarkably, direct delivery of the U178G/U179G mutant into the vascular system by needle puncture inoculation allowed it to spread systemically and enter mesophyll cells and epidermal cells of upper leaves. These findings highlight the importance of RNA 3D structure for PSTVd replication and intercellular trafficking and indicate that loop 27 is required for epidermal exit, but not epidermal entry or transit between other cell types. Thus, requirements for RNA trafficking between epidermal and underlying palisade mesophyll cells are unique and directional. Our findings further suggest that 3D structure and RNA-protein interactions constrain RNA sequence evolution, and validate JAR3D as a tool to predict RNA 3D structure.

Multiomics resolution of molecular events during a day in the life of Chlamydomonas

The unicellular green alga Chlamydomonas reinhardtii displays metabolic flexibility in response to a changing environment. We analyzed expression patterns of its three genomes in cells grown under light-dark cycles. Nearly 85% of transcribed genes show differential expression, with different sets of transcripts being up-regulated over the course of the day to coordinate cellular growth before undergoing cell division. Parallel measurements of select metabolites and pigments, physiological parameters, and a subset of proteins allow us to infer metabolic events and to evaluate the impact of the transcriptome on the proteome. Among the findings are the observations that Chlamydomonas exhibits lower respiratory activity at night compared with the day; multiple fermentation pathways, some oxygen-sensitive, are expressed at night in aerated cultures; we propose that the ferredoxin, FDX9, is potentially the electron donor to hydrogenases. The light stress-responsive genes PSBS, LHCSR1, and LHCSR3 show an acute response to lights-on at dawn under abrupt dark-to-light transitions, while LHCSR3 genes also exhibit a later, second burst in expression in the middle of the day dependent on light intensity. Each response to light (acute and sustained) can be selectively activated under specific conditions. Our expression dataset, complemented with coexpression networks and metabolite profiling, should constitute an excellent resource for the algal and plant communities.

Dissecting plant genomes with the PLAZA comparative genomics platform

With the arrival of low-cost, next-generation sequencing, a multitude of new plant genomes are being publicly released, providing unseen opportunities and challenges for comparative genomics studies. Here, we present PLAZA 2.5, a user-friendly online research environment to explore genomic information from different plants. This new release features updates to previous genome annotations and a substantial number of newly available plant genomes as well as various new interactive tools and visualizations. Currently, PLAZA hosts 25 organisms covering a broad taxonomic range, including 13 eudicots, five monocots, one lycopod, one moss, and five algae. The available data consist of structural and functional gene annotations, homologous gene families, multiple sequence alignments, phylogenetic trees, and colinear regions within and between species. A new Integrative Orthology Viewer, combining information from different orthology prediction methodologies, was developed to efficiently investigate complex orthology relationships. Cross-species expression analysis revealed that the integration of complementary data types extended the scope of complex orthology relationships, especially between more distantly related species. Finally, based on phylogenetic profiling, we propose a set of core gene families within the green plant lineage that will be instrumental to assess the gene space of draft or newly sequenced plant genomes during the assembly or annotation phase.

Targeted gene silencing by RNA interference in Chlamydomonas

Small RNA-guided gene silencing is an evolutionarily conserved process that operates by a variety of molecular mechanisms and plays an essential role in developmental pathways and defense responses against genomic parasites in eukaryotes. Double-stranded RNA (dsRNA) triggered posttranscriptional gene silencing, termed RNA interference (RNAi), is also becoming a powerful tool for reverse genetics studies. Stable RNAi, induced by the expression of long dsRNAs or duplex small RNAs from genome-integrated transgenes, has been achieved in multiple organisms, including the green alga Chlamydomonas reinhardtii. However, the level of gene repression is often quite variable, depending on the type of construct, transgene copy number, site of integration, and target gene. Moreover, unintended transcripts partly complementary to a trigger dsRNA can also be silenced, making difficult the interpretation of observed phenotypes. To obviate some of these problems we have developed a tandem inverted repeat system that consistently induces cosilencing of a gene with a selectable RNAi-induced phenotype (encoding tryptophan synthase beta-subunit) and any other (nonessential) gene of interest. In addition, to circumvent off-target effects, for each tested gene, RNAi lines are generated with at least two transgenes, homologous to distinct and nonoverlapping sequences of the target transcript. A common phenotype among these independent RNAi strains is expected to result from suppression of expression of the gene of interest. We demonstrate this approach for the characterization of a gene of unknown function in Chlamydomonas, encoding a predicted exoribonuclease with weak similarity to 3'hExo/ERI-1.

PLAZA: a comparative genomics resource to study gene and genome evolution in plants

The number of sequenced genomes of representatives within the green lineage is rapidly increasing. Consequently, comparative sequence analysis has significantly altered our view on the complexity of genome organization, gene function, and regulatory pathways. To explore all this genome information, a centralized infrastructure is required where all data generated by different sequencing initiatives is integrated and combined with advanced methods for data mining. Here, we describe PLAZA, an online platform for plant comparative genomics (http://bioinformatics.psb.ugent.be/plaza/). This resource integrates structural and functional annotation of published plant genomes together with a large set of interactive tools to study gene function and gene and genome evolution. Precomputed data sets cover homologous gene families, multiple sequence alignments, phylogenetic trees, intraspecies whole-genome dot plots, and genomic colinearity between species. Through the integration of high confidence Gene Ontology annotations and tree-based orthology between related species, thousands of genes lacking any functional description are functionally annotated. Advanced query systems, as well as multiple interactive visualization tools, are available through a user-friendly and intuitive Web interface. In addition, detailed documentation and tutorials introduce the different tools, while the workbench provides an efficient means to analyze user-defined gene sets through PLAZA's interface. In conclusion, PLAZA provides a comprehensible and up-to-date research environment to aid researchers in the exploration of genome information within the green plant lineage.

Metabolism and regulation of canonical histone mRNAs: life without a poly(A) tail

The canonical histone proteins are encoded by replication-dependent genes and must rapidly reach high levels of expression during S phase. In metazoans the genes that encode these proteins produce mRNAs that, instead of being polyadenylated, contain a unique 3' end structure. By contrast, the synthesis of the variant, replication-independent histones, which are encoded by polyadenylated mRNAs, persists outside of S phase. Accurate positioning of both histone types in chromatin is essential for proper transcriptional regulation, the demarcation of heterochromatic boundaries and the epigenetic inheritance of gene expression patterns. Recent results suggest that the coordinated synthesis of replication-dependent and variant histone mRNAs is achieved by signals that affect formation of the 3' end of the replication-dependent histone mRNAs.

U7 snRNAs: a computational survey

U7 small nuclear RNA (snRNA) sequences have been described only for a handful of animal species in the past. Here we describe a computational search for functional U7 snRNA genes throughout vertebrates including the upstream sequence elements characteristic for snRNAs transcribed by polymerase II. Based on the results of this search, we discuss the high variability of U7 snRNAs in both sequence and structure, and report on an attempt to find U7 snRNA sequences in basal deuterostomes and non-drosophilids insect genomes based on a combination of sequence, structure, and promoter features. Due to the extremely short sequence and the high variability in both sequence and structure, no unambiguous candidates were found. These results cast doubt on putative U7 homologs in even more distant organisms that are reported in the most recent release of the Rfam database.

Early evolution of histone mRNA 3' end processing

The replication-dependent histone mRNAs in metazoa are not polyadenylated, in contrast to the bulk of mRNA. Instead, they contain an RNA stem-loop (SL) structure close to the 3' end of the mature RNA, and this 3' end is generated by cleavage using a machinery involving the U7 snRNP and protein factors such as the stem-loop binding protein (SLBP). This machinery of 3' end processing is related to that of polyadenylation as protein components are shared between the systems. It is commonly believed that histone 3' end processing is restricted to metazoa and green algae. In contrast, polyadenylation is ubiquitous in Eukarya. However, using computational approaches, we have now identified components of histone 3' end processing in a number of protozoa. Thus, the histone mRNA stem-loop structure as well as the SLBP protein are present in many different protozoa, including Dictyostelium, alveolates, Trypanosoma, and Trichomonas. These results show that the histone 3' end processing machinery is more ancient than previously anticipated and can be traced to the root of the eukaryotic phylogenetic tree. We also identified histone mRNAs from both metazoa and protozoa that are polyadenylated but also contain the signals characteristic of histone 3' end processing. These results provide further evidence that some histone genes are regulated at the level of 3' end processing to produce either polyadenylated RNAs or RNAs with the 3' end characteristic of replication-dependent histone mRNAs.

Histone H3 and H3.3 variants in the protozoan pathogens Plasmodium falciparum and Toxoplasma gondii

Histones constitute the fundamental component of chromatin and participate in the regulation of gene expression by virtue of covalent modifications to their N-terminal domains. The discovery that histone-modifying enzymes are targeted by the antiprotozoal agent apicidin has prompted further investigation of gene expression regulation in protozoan parasites; consequently, several chromatin remodeling homologues with unusual features have been isolated. To facilitate investigation of these chromatin remodeling homologues using parasite-specific substrates, we sought to clone and characterize histone H3 from two medically significant pathogens in the phylum Apicomplexa: Plasmodium falciparum (malaria) and Toxoplasma gondii (opportunistic pathogen of immunocompromised individuals). Like most eukaryotic organisms, these parasites each contain at least two histone H3 variants, termed H3 and H3.3. Sequence analysis reveals the Apicomplexan H3 proteins harbor novel and rare features. Expression and purification of recombinant H3 variants will provide species-specific substrate for the analysis of the histone-modifying machinery of these parasites.

DC3, the 21-kDa subunit of the outer dynein arm-docking complex (ODA-DC), is a novel EF-hand protein important for assembly of both the outer arm and the ODA-DC

The outer dynein arm-docking complex (ODA-DC) is a microtubule-associated structure that targets the outer dynein arm to its binding site on the flagellar axoneme (Takada et al. 2002. Mol. Biol. Cell 13, 1015-1029). The ODA-DC of Chlamydomonas contains three proteins, referred to as DC1, DC2, and DC3. We here report the isolation and sequencing of genomic and full-length cDNA clones encoding DC3. The sequence predicts a 21,341 Da protein with four EF-hands that is a member of the CTER (calmodulin, troponin C, essential and regulatory myosin light chains) group and is most closely related to a predicted protein from Plasmodium. The DC3 gene, termed ODA14, is intronless. Chlamydomonas mutants that lack DC3 exhibit slow, jerky swimming because of loss of some but not all outer dynein arms. Some outer doublet microtubules without arms had a "partial" docking complex, indicating that DC1 and DC2 can assemble in the absence of DC3. In contrast, DC3 cannot assemble in the absence of DC1 or DC2. Transformation of a DC3-deletion strain with the wild-type DC3 gene rescued both the motility phenotype and the structural defect, whereas a mutated DC3 gene was incompetent to rescue. The results indicate that DC3 is important for both outer arm and ODA-DC assembly.

Core histones of the amitochondriate protist, Giardia lamblia

Genes coding for the core histones H2a, H2b, H3, and H4 of Giardia lamblia were sequenced. A conserved organism- and gene-specific element, GRGCGCAGATTTVGG, was found upstream of the coding region in all core histone genes. The derived amino acid sequences of all four histones were similar to their homologs in other eukaryotes, although they were among the most divergent members of this protein family. Comparative protein structure modeling combined with energy evaluation of the resulting models indicated that the G. lamblia core histones individually and together can assume the same three-dimensional structures that were established by X-ray crystallography for Xenopus laevis histones and the nucleosome core particle. Since G. lamblia represents one of the earliest-diverging eukaryotes in many different molecular trees, the structure of its histones is potentially of relevance to understanding histone evolution. The G. lamblia proteins do not represent an intermediate stage between archaeal and eukaryotic histones.

Formation of the 3' end of histone mRNA

All metazoan messenger RNAs, with the exception of the replication-dependent histone mRNAs, terminate at the 3' end with a poly(A) tail. Replication-dependent histone mRNAs end instead in a conserved 26-nucleotide sequence that contains a 16-nucleotide stem-loop. Formation of the 3' end of histone mRNA occurs by endonucleolytic cleavage of pre-mRNA releasing the mature mRNA from the chromatin template. Cleavage requires several trans-acting factors, including a protein, the stem-loop binding protein (SLBP), which binds the 26-nucleotide sequence; and a small nuclear RNP, U7 snRNP. There are probably additional factors also required for cleavage. One of the functions of the SLBP is to stabilize binding of the U7 snRNP to the histone pre-mRNA. In the nucleus, both U7 snRNP and SLBP are present in coiled bodies, structures that are associated with histone genes and may play a direct role in histone pre-mRNA processing in vivo. One of the major regulatory events in the cell cycle is regulation of histone pre-mRNA processing, which is at least partially mediated by cell-cycle regulation of the levels of the SLBP protein.

Dynamics of histone acetylation in Chlamydomonas reinhardtii

The dynamic character of core histone post-translational acetylation in the unicellular green alga Chlamydomonas reinhardtii was studied by tritiated acetate incorporation. Histone H3 is the major target of acetylation, steady state, and in pulse and pulse-chase analyses. Acetylation turnover rates were measured by tracer labeling under steady-state conditions. Half-lives of 1.5-3 min were found for penta- to mono-acetylation of H3, dynamically acetylated to the 30% level. Twenty percent of H3 was multi-acetylated, on average with 3. 2 acetyl-lysines, all with rapid turnover. Deacetylase inhibitor trichostatin A (TSA) caused doubling of average acetylation levels, primarily as penta-acetylated H3, but half of H3 was not acetylated at all. The level of histone H4 acetylation was only half that of H3 and a major fraction of mono- and di-acetylated forms appeared static. The dynamic fraction had an average half-life of 3.5 min with higher turnover rates for more highly acetylated H4 forms. TSA, inhibiting less effectively deacetylases active on H4, strongly increased multi-acetylated H4 levels and doubled average acetylation. As for H3, half of histone H4 remained unacetylated. Acetylation of histone H2B was low and of H2A was barely measurable. Despite turnover with half-lives of approximately 2 min, no increase beyond di-acetylation was seen upon TSA treatment.

Translational control of cellular and viral mRNAs

We are becoming increasingly aware of the role that translational control plays in regulating gene expression in plants. There are now many examples in which specific mechanisms have evolved at the translational level that directly impact the amount of protein produced from an mRNA. All regions of an mRNA, i.e., the 5' leader, the coding region, and the 3'-untranslated region, have the potential to influence translation. The 5'-terminal cap structure and the poly(A) tail at the 3' terminus serve as additional elements controlling translation. Many viral mRNAs have evolved alternatives to the cap and poly(A) tail that are functionally equivalent. Nevertheless, for both cellular and viral mRNAs, a co-dependent interaction between the terminal controlling elements appears to be the universal basis for efficient translation.

Common features of analogous replacement histone H3 genes in animals and plants

Phylogenetic analysis of histone H3 protein sequences demonstrates the independent origin of the replacement histone H3 genes in animals and in plants. Multiple introns in the replacement histone H3 genes of animals in a pattern distinct from that in plant replacement H3 genes supports this conclusion. It is suggested that replacement H3 genes arose at the same time that, independently, multicellular forms of animals and of plants evolved. Judged by the degree of invariant and functionally constrained amino acid positions, histones H3 and H4, which form together the tetramer kernel of the nucleosome, have co-evolved with equal rates of sequence divergence. Residues 31 and 87 in histone H3 are the only residues that consistently changed across each gene duplication event that created functional replacement histone H3 variant forms. Once changed, these residues have remained invariant across divergent speciation. This suggests that they are required to allow replacement histone H3 to participate in the assembly of nucleosomes in non-S-phase cells. The abundant occurrence of polypyrimidine sequences in the introns of all replacement H3 genes, and the replacement of an intron by a polypyrimidine motif upstream of the alfalfa replacement H3 gene, suggests a function. It is speculated that they may contribute to the characteristic cell-cycle-independent pattern of replacement histone H3 genes by binding nucleosome-excluding proteins.

The histone 3'-terminal stem-loop is necessary for translation in Chinese hamster ovary cells

The metazoan cell cycle-regulated histone mRNAs are the only known cellular mRNAs that do not terminate in a poly(A) tall. Instead, mammalian histone mRNAs terminate in a highly conserved stem-loop structure which is required for 3'-end processing and regulates mRNA stability. The poly(A) tail not only regulates translational efficiency and mRNA stability but is required for the function of the cap in translation (m(7)GpppN). We show that the histone terminal stem-loop is functionally similar to a poly(A) tail in that it enhances translational efficiency and is co-dependent on a cap in order to establish an efficient level of translation. The histone stem-loop is sufficient and necessary to increase the translation of reporter mRNA in transfected Chinese hamster ovary cells but must be positioned at the 3'-terminus in order to function optimally. Mutations within the conserved stem or loop regions reduced its ability to facilitate translation. All histone mRNAs in higher plants are polyadenylated. The histone stem-loop did not function to influence translational efficiency or mRNA stability in plant protoplasts. These data demonstrate that the histone stem/loop directs efficient translation and that it is functionally analogous to a poly(A) tail.