Transcription of damage-induced RNA in Arabidopsis was frequently initiated from DSB loci within the genic regions.

DNA double-strand breaks (DSBs) are the most severe DNA lesions and need to be removed immediately to prevent loss of genomic information. Recently, it has been revealed that DSBs induce novel transcription from the cleavage sites in various species, resulting in RNAs being referred to as damage-induced RNAs (diRNAs). While diRNA synthesis is an early event in the DNA damage response and plays an essential role in DSB repair activation, the location where diRNAs are newly generated in plants remains unclear, as does their transcriptional mechanism. Here, we performed the sequencing of polyadenylated (polyA) diRNAs that emerged around all DSB loci in Arabidopsis thaliana under the expression of the exogenous restriction enzyme Sbf I and observed 88 diRNAs transcribed via RNA polymerase II in 360 DSB loci. Most of the detected diRNAs originated within active genes and were transcribed from DSBs in a bidirectional manner. Furthermore, we found that diRNA elongation tends to terminate at the boundary of an endogenous gene located near DSB loci. Our results provide reliable evidence for understanding the importance of new transcription at DSBs and show that diRNA is a crucial factor for successful DSB repair.

Inducible Expression of the Restriction Enzyme Uncovered Genome-Wide Distribution and Dynamic Behavior of Histones H4K16ac and H2A.Z at DNA Double-Strand Breaks in Arabidopsis.

DNA double-strand breaks (DSBs) are among the most serious types of DNA damage, causing mutations and chromosomal rearrangements. In eukaryotes, DSBs are immediately repaired in coordination with chromatin remodeling for the deposition of DSB-related histone modifications and variants. To elucidate the details of DSB-dependent chromatin remodeling throughout the genome, artificial DSBs need to be reproducibly induced at various genomic loci. Recently, a comprehensive method for elucidating chromatin remodeling at multiple DSB loci via chemically induced expression of a restriction enzyme was developed in mammals. However, this DSB induction system is unsuitable for investigating chromatin remodeling during and after DSB repair, and such an approach has not been performed in plants. Here, we established a transgenic Arabidopsis plant harboring a restriction enzyme gene Sbf I driven by a heat-inducible promoter. Using this transgenic line, we performed chromatin immunoprecipitation followed by deep sequencing (ChIP-seq) of histones H4K16ac and H2A.Z and investigated the dynamics of these histone marks around the endogenous 623 Sbf I recognition sites. We also precisely quantified DSB efficiency at all cleavage sites using the DNA resequencing data obtained by the ChIP-seq procedure. From the results, Sbf I-induced DSBs were detected at 360 loci, which induced the transient deposition of H4K16ac and H2A.Z around these regions. Interestingly, we also observed the co-localization of H4K16ac and H2A.Z at some DSB loci. Overall, DSB-dependent chromatin remodeling was found to be highly conserved between plants and animals. These findings provide new insights into chromatin remodeling that occurs in response to DSBs in Arabidopsis.

Phylogeny analysis of whole protein-coding genes in metagenomic data detected an environmental gradient for the microbiota.

Environmental factors affect the growth of microorganisms and therefore alter the composition of microbiota. Correlative analysis of the relationship between metagenomic composition and the environmental gradient can help elucidate key environmental factors and establishment principles for microbial communities. However, a reasonable method to quantitatively compare whole metagenomic data and identify the primary environmental factors for the establishment of microbiota has not been reported so far. In this study, we developed a method to compare whole proteomes deduced from metagenomic shotgun sequencing data, and quantitatively display their phylogenetic relationships as metagenomic trees. We called this method Metagenomic Phylogeny by Average Sequence Similarity (MPASS). We also compared one of the metagenomic trees with dendrograms of environmental factors using a comparison tool for phylogenetic trees. The MPASS method correctly constructed metagenomic trees of simulated metagenomes and soil and water samples. The topology of the metagenomic tree of samples from the Kirishima hot springs area in Japan was highly similarity to that of the dendrograms based on previously reported environmental factors for this area. The topology of the metagenomic tree also reflected the dynamics of microbiota at the taxonomic and functional levels. Our results strongly suggest that MPASS can successfully classify metagenomic shotgun sequencing data based on the similarity of whole protein-coding sequences, and will be useful for the identification of principal environmental factors for the establishment of microbial communities. Custom Perl script for the MPASS pipeline is available at https://github.com/s0sat/MPASS.

Cryptic promoter activation occurs by at least two different mechanisms in the Arabidopsis genome.

In gene-trap screening of plant genomes, promoterless reporter constructs are often expressed without trapping of annotated gene promoters. The molecular basis of this phenomenon, which has been interpreted as the trapping of cryptic promoters, is poorly understood. Here, we found that cryptic promoter activation occurs by at least two different mechanisms using Arabidopsis gene-trap lines in which a firefly luciferase (LUC) open reading frame (ORF) without an apparent promoter sequence was expressed from intergenic regions: one mechanism is 'cryptic promoter capturing', in which the LUC ORF captured pre-existing promoter-like chromatin marked by H3K4me3 and H2A.Z, and the other is 'promoter de novo origination', in which the promoter chromatin was newly formed near the 5' end of the inserted LUC ORF. The latter finding raises a question as to how the inserted LUC ORF sequence is involved in this phenomenon. To examine this, we performed a model experiment with chimeric LUC genes in transgenic plants. Using Arabidopsis psaH1 promoter-LUC constructs, we found that the functional core promoter region, where transcription start sites (TSSs) occur, cannot simply be determined by the upstream nor core promoter sequences; rather, its positioning proximal to the inserted LUC ORF sequence was more critical. This result suggests that the insertion of the coding sequence alters the local distribution of TSSs in the plant genome. The possible impact of the two types of cryptic promoter activation mechanisms on plant genome evolution and endosymbiotic gene transfer is discussed.

De novo activated transcription of inserted foreign coding sequences is inheritable in the plant genome.

The manner in which inserted foreign coding sequences become transcriptionally activated and fixed in the plant genome is poorly understood. To examine such processes of gene evolution, we performed an artificial evolutionary experiment in Arabidopsis thaliana. As a model of gene-birth events, we introduced a promoterless coding sequence of the firefly luciferase (LUC) gene and established 386 T2-generation transgenic lines. Among them, we determined the individual LUC insertion loci in 76 lines and found that one-third of them were transcribed de novo even in the intergenic or inherently unexpressed regions. In the transcribed lines, transcription-related chromatin marks were detected across the newly activated transcribed regions. These results agreed with our previous findings in A. thaliana cultured cells under a similar experimental scheme. A comparison of the results of the T2-plant and cultured cell experiments revealed that the de novo-activated transcription concomitant with local chromatin remodelling was inheritable. During one-generation inheritance, it seems likely that the transcription activities of the LUC inserts trapped by the endogenous genes/transcripts became stronger, while those of de novo transcription in the intergenic/untranscribed regions became weaker. These findings may offer a clue for the elucidation of the mechanism by which inserted foreign coding sequences become transcriptionally activated and fixed in the plant genome.

Kozak Sequence Acts as a Negative Regulator for De Novo Transcription Initiation of Newborn Coding Sequences in the Plant Genome.

The manner in which newborn coding sequences and their transcriptional competency emerge during the process of gene evolution remains unclear. Here, we experimentally simulated eukaryotic gene origination processes by mimicking horizontal gene transfer events in the plant genome. We mapped the precise position of the transcription start sites (TSSs) of hundreds of newly introduced promoterless firefly luciferase (LUC) coding sequences in the genome of Arabidopsis thaliana cultured cells. The systematic characterization of the LUC-TSSs revealed that 80% of them occurred under the influence of endogenous promoters, while the remainder underwent de novo activation in the intergenic regions, starting from pyrimidine-purine dinucleotides. These de novo TSSs obeyed unexpected rules; they predominantly occurred ∼100 bp upstream of the LUC inserts and did not overlap with Kozak-containing putative open reading frames (ORFs). These features were the output of the immediate responses to the sequence insertions, rather than a bias in the screening of the LUC gene function. Regarding the wild-type genic TSSs, they appeared to have evolved to lack any ORFs in their vicinities. Therefore, the repulsion by the de novo TSSs of Kozak-containing ORFs described above might be the first selection gate for the occurrence and evolution of TSSs in the plant genome. Based on these results, we characterized the de novo type of TSS identified in the plant genome and discuss its significance in genome evolution.

Correction to: Space-time analysis of gravitropism in etiolated Arabidopsis hypocotyls using bioluminescence imaging of the IAA19 promoter fusion with a destabilized luciferase reporter.

The article Space-time analysis of gravitropism in etiolated Arabidopsis hypocotyls using bioluminescence imaging of the IAA19 promoter fusion with a destabilized luciferase reporter, written by Kotaro T. Yamamoto, Masaaki K. Watahiki, Jun Matsuzaki, Soichirou Satoh and Hisayo Shimizu, was originally published electronically on the publisher's internet portal (currently SpringerLink) on 10 April 2017 without open access.

Characterization of spliced leader trans-splicing in a photosynthetic rhizarian amoeba, Paulinella micropora, and its possible role in functional gene transfer.

Paulinella micropora is a rhizarian thecate amoeba, belonging to a photosynthetic Paulinella species group that has a unique organelle termed chromatophore, whose cyanobacterial origin is distinct from that of plant and algal chloroplasts. Because acquisition of the chromatophore was quite a recent event compared with that of the chloroplast ancestor, the Paulinella species are thought to be model organisms for studying the early process of primary endosymbiosis. To obtain insight into how endosymbiotically transferred genes acquire expression competence in the host nucleus, here we analyzed the 5' end sequences of the mRNAs of P. micropora MYN1 strain with the aid of a cap-trapper cDNA library. As a result, we found that mRNAs of 27 genes, including endosymbiotically transferred genes, possessed the common 5' end sequence of 28-33 bases that were posttranscriptionally added by spliced leader (SL) trans-splicing. We also found two subtypes of SL RNA genes encoded by the P. micropora MYN1 genome. Differing from the other SL trans-splicing organisms that usually possess poly(A)-less SL RNAs, this amoeba has polyadenylated SL RNAs. In this study, we characterize the SL trans-splicing of this unique organism and discuss the putative merits of SL trans-splicing in functional gene transfer and genome evolution.

Comparative analyses of whole-genome protein sequences from multiple organisms.

Phylogenies based on entire genomes are a powerful tool for reconstructing the Tree of Life. Several methods have been proposed, most of which employ an alignment-free strategy. Average sequence similarity methods are different than most other whole-genome methods, because they are based on local alignments. However, previous average similarity methods fail to reconstruct a correct phylogeny when compared against other whole-genome trees. In this study, we developed a novel average sequence similarity method. Our method correctly reconstructs the phylogenetic tree of in silico evolved E. coli proteomes. We applied the method to reconstruct a whole-proteome phylogeny of 1,087 species from all three domains of life, Bacteria, Archaea, and Eucarya. Our tree was automatically reconstructed without any human decisions, such as the selection of organisms. The tree exhibits a concentric circle-like structure, indicating that all the organisms have similar total branch lengths from their common ancestor. Branching patterns of the members of each phylum of Bacteria and Archaea are largely consistent with previous reports. The topologies are largely consistent with those reconstructed by other methods. These results strongly suggest that this approach has sufficient taxonomic resolution and reliability to infer phylogeny, from phylum to strain, of a wide range of organisms.

Space-time analysis of gravitropism in etiolated Arabidopsis hypocotyls using bioluminescence imaging of the IAA19 promoter fusion with a destabilized luciferase reporter.

Imaging analysis was carried out during the gravitropic response of etiolated Arabidopsis hypocotyls, using an IAA19 promoter fusion of destabilized luciferase as a probe. From the bright-field images we obtained the local deflection angle to the vertical, A, local curvature, C, and the partial derivative of C with respect to time, [Formula: see text]. These were determined every 19.9 µm along the curvilinear length of the hypocotyl, at ~10 min intervals over a period of ~6 h after turning hypocotyls through 90° to the horizontal. Similarly from the luminescence images we measured the luminescence intensity of the convex and concave flanks of the hypocotyl as well as along the median of the hypocotyl, to determine differential expression of auxin-inducible IAA19. Comparison of these parameters as a function of time and curvilinear length shows that the gravitropic response is composed of three successive elements: the first and second curving responses and a decurving response (autostraightening). The maximum of the first curving response occurs when A is 76° along the entire length of the hypocotyl, suggesting that A is the sole determinant in this response; in contrast, the decurving response is a function of both A and C, as predicted by the newly-proposed graviproprioception model (Bastien et al., Proc Natl Acad Sci USA 110:755-760, 2013). Further, differential expression of IAA19, with higher expression in the convex flank, is observed at A = 44°, and follows the Sachs' sine law. This also suggests that IAA19 is not involved in the first curving response. In summary, the gravitropic response of Arabidopsis hypocotyls consists of multiple elements that are each determined by separate principles.

Evolution of Green Plants Accompanied Changes in Light-Harvesting Systems.

Photosynthetic organisms have various pigments enabling them to adapt to various light environments. Green plants are divided into two groups: streptophytes and chlorophytes. Streptophytes include some freshwater green algae and land plants, while chlorophytes comprise the other freshwater green algae and seawater green algae. The environmental conditions driving the divergence of green plants into these two groups and the changes in photosynthetic properties accompanying their evolution remain unknown. Here, we separated the core antennae of PSI and the peripheral antennae [light-harvesting complexes (LHCs)] in green plants by green-native gel electrophoresis and determined their pigment compositions. Freshwater green algae and land plants have high Chl a/b ratios, with most Chl b existing in LHCs. In contrast, seawater green algae have low Chl a/b ratios. In addition, Chl b exists not only in LHCs but also in PSI core antennae in these organisms, a situation beneficial for survival in deep seawater, where blue-green light is the dominant light source. Finally, low-energy Chl (red Chl) of PSI was detected in freshwater green algae and land plants, but not in seawater green algae. We thus conclude that the different level of Chl b accumulation in core antennae and differences in PSI red Chl between freshwater and seawater green algae are evolutionary adaptations of these algae to their habitats, especially to high- or low-light environments.

Construction of a phylogenetic tree of photosynthetic prokaryotes based on average similarities of whole genome sequences.

Phylogenetic trees have been constructed for a wide range of organisms using gene sequence information, especially through the identification of orthologous genes that have been vertically inherited. The number of available complete genome sequences is rapidly increasing, and many tools for construction of genome trees based on whole genome sequences have been proposed. However, development of a reasonable method of using complete genome sequences for construction of phylogenetic trees has not been established. We have developed a method for construction of phylogenetic trees based on the average sequence similarities of whole genome sequences. We used this method to examine the phylogeny of 115 photosynthetic prokaryotes, i.e., cyanobacteria, Chlorobi, proteobacteria, Chloroflexi, Firmicutes and nonphotosynthetic organisms including Archaea. Although the bootstrap values for the branching order of phyla were low, probably due to lateral gene transfer and saturated mutation, the obtained tree was largely consistent with the previously reported phylogenetic trees, indicating that this method is a robust alternative to traditional phylogenetic methods.

Pentatricopeptide repeat proteins involved in plant organellar RNA editing.

C-to-U RNA editing has been widely observed in organellar RNAs in terrestrial plants. Recent research has revealed the significance of a large, plant-specific family of pentatricopeptide repeat (PPR) proteins for RNA editing and other RNA processing events in plant mitochondria and chloroplasts. PPR protein is a sequence-specific RNA-binding protein that identifies specific C residues for editing. Discovery of the RNA recognition code for PPR motifs, including verification and prediction of the individual RNA editing site and its corresponding PPR protein, expanded our understanding of the molecular function of PPR proteins in plant organellar RNA editing. Using this knowledge and the co-expression database, we have identified two new PPR proteins that mediate chloroplast RNA editing. Further, computational target assignment using the PPR RNA recognition codes suggests a distinct, unknown mode-of-action, by which PPR proteins serve a function beyond site recognition in RNA editing.

Niche adaptation and genome expansion in the chlorophyll d-producing cyanobacterium Acaryochloris marina.

Acaryochloris marina is a unique cyanobacterium that is able to produce chlorophyll d as its primary photosynthetic pigment and thus efficiently use far-red light for photosynthesis. Acaryochloris species have been isolated from marine environments in association with other oxygenic phototrophs, which may have driven the niche-filling introduction of chlorophyll d. To investigate these unique adaptations, we have sequenced the complete genome of A. marina. The DNA content of A. marina is composed of 8.3 million base pairs, which is among the largest bacterial genomes sequenced thus far. This large array of genomic data is distributed into nine single-copy plasmids that code for >25% of the putative ORFs. Heavy duplication of genes related to DNA repair and recombination (primarily recA) and transposable elements could account for genetic mobility and genome expansion. We discuss points of interest for the biosynthesis of the unusual pigments chlorophyll d and alpha-carotene and genes responsible for previously studied phycobilin aggregates. Our analysis also reveals that A. marina carries a unique complement of genes for these phycobiliproteins in relation to those coding for antenna proteins related to those in Prochlorococcus species. The global replacement of major photosynthetic pigments appears to have incurred only minimal specializations in reaction center proteins to accommodate these alternate pigments. These features clearly show that the genus Acaryochloris is a fitting candidate for understanding genome expansion, gene acquisition, ecological adaptation, and photosystem modification in the cyanobacteria.

Identification of chlorophyllide a oxygenase in the Prochlorococcus genome by a comparative genomic approach.

Chl b is a major photosynthetic pigment of peripheral antenna complexes in chlorophytes and prochlorophytes. Chl b is synthesized by chlorophyllide a oxygenase (CAO), an enzyme that has been identified from higher plants, moss, green algae and two groups of prochlorophytes, Prochlorothrix and Prochloron. Based on these results, we previously proposed the hypothesis that all of the Chl b synthesis genes have a common origin. However, the CAO gene is not found in whole genome sequences of Prochlorococcus although a gene which is distantly related to CAO was reported. If Prochlorococcus employs a different enzyme, a Chl synthesis gene should have evolved several times on the different phylogenetic lineages of Prochlorococcus and other Chl b-containing organisms. To examine these hypotheses, we identified a Prochlorococcus Chl b synthesis gene by using a combination of bioinformatics and molecular genetics techniques. We first identified Prochlorococcus-specific genes by comparing the whole genome sequences of Prochlorococcus marinus MED4, MIT9313 and SS120 with Synechococcus sp. WH8102. Synechococcus is closely related to Prochlorococcus phylogenetically, but it does not contain a Chl b synthesis gene. By examining the sequences of Prochlorococcus-specific genes, we found a candidate for the Chl b synthesis gene and introduced it into Synechocystis sp. PCC6803. The transformant cells accumulated Chl b, indicating that the gene product catalyzes Chl b synthesis. In this study, we discuss the evolution of CAO based upon the molecular phylogenetic studies we performed.

Pigment shuffling in antenna systems achieved by expressing prokaryotic chlorophyllide a oxygenase in Arabidopsis.

The organization of pigment molecules in photosystems is strictly determined. The peripheral antennae have both chlorophyll a and b, but the core antennae consist of only chlorophyll a in green plants. Furthermore, according to the recent model obtained from the crystal structure of light-harvesting chlorophyll a/b-protein complexes II (LHCII), individual chlorophyll-binding sites are occupied by either chlorophyll a or chlorophyll b. In this study, we succeeded in altering these pigment organizations by introducing a prokaryotic chlorophyll b synthesis gene (chlorophyllide a oxygenase (CAO)) into Arabidopsis. In these transgenic plants (Prochlirothrix hollandica CAO plants), approximately 40% of chlorophyll a of the core antenna complexes was replaced by chlorophyll b in both photosystems. Chlorophyll a/b ratios of LHCII also decreased from 1.3 to 0.8 in PhCAO plants. Surprisingly, these transgenic plants were capable of photosynthetic growth similar to wild type under low light conditions. These results indicate that chlorophyll organizations are not solely determined by the binding affinities, but they are also controlled by CAO. These data also suggest that strict organizations of chlorophyll molecules are not essential for photosynthesis under low light conditions.

Identification of a vinyl reductase gene for chlorophyll synthesis in Arabidopsis thaliana and implications for the evolution of Prochlorococcus species.

Chlorophyll metabolism has been extensively studied with various organisms, and almost all of the chlorophyll biosynthetic genes have been identified in higher plants. However, only the gene for 3,8-divinyl protochlorophyllide a 8-vinyl reductase (DVR), which is indispensable for monovinyl chlorophyll synthesis, has not been identified yet. In this study, we isolated an Arabidopsis thaliana mutant that accumulated divinyl chlorophyll instead of monovinyl chlorophyll by ethyl methanesulfonate mutagenesis. Map-based cloning of this mutant resulted in the identification of a gene (AT5G18660) that shows sequence similarity with isoflavone reductase genes. The mutant phenotype was complemented by the transformation with the wild-type gene. A recombinant protein encoded by AT5G18660 was expressed in Escherichia coli and found to catalyze the conversion of divinyl chlorophyllide to monovinyl chlorophyllide, thereby demonstrating that the gene encodes a functional DVR. DVR is encoded by a single copy gene in the A. thaliana genome. With the identification of DVR, finally all genes required for chlorophyll biosynthesis have been identified in higher plants. Analysis of the complete genome of A. thaliana showed that it has 15 enzymes encoded by 27 genes for chlorophyll biosynthesis from glutamyl-tRNA(glu) to chlorophyll b. Furthermore, identification of the DVR gene helped understanding the evolution of Prochlorococcus marinus, a marine cyanobacterium that is dominant in the open ocean and is uncommon in using divinyl chlorophylls. A DVR homolog was not found in the genome of P. marinus but found in the Synechococcus sp WH8102 genome, which is consistent with the distribution of divinyl chlorophyll in marine cyanobacteria of the genera Prochlorococcus and Synechococcus.

Unique constitution of photosystem I with a novel subunit in the cyanobacterium Gloeobacter violaceus PCC 7421.

Constitution of the photosystem I complex isolated from the cyanobacterium Gloeobacter violaceus PCC 7421 was investigated by tricine-urea-SDS-PAGE, followed by peptide mass fingerprinting or N-terminal sequencing. Eight subunits (PsaA, PsaB, PsaC, PsaD, PsaE, PsaF, PsaL and PsaM) were identified as predicted from the genome sequence. A novel subunit (PsaZ) was discovered, but PsaI, PsaJ, PsaK and PsaX were absent. PsaB has a C-terminal extension with 155 amino acids in addition to the conserved region and this domain is similar to the peptidoglycan-binding domain. These results suggest that PS I complexes of G. violaceus have unique structural properties.

The Arabidopsis-accelerated cell death gene ACD1 is involved in oxygenation of pheophorbide a: inhibition of the pheophorbide a oxygenase activity does not lead to the "stay-green" phenotype in Arabidopsis.

Oxygenation of pheophorbide a is a key step in chlorophyll breakdown. Several biochemical studies have implicated that this step was catalyzed by an iron-containing and ferredoxin-dependent monooxygenase, pheophorbide a oxygenase (PaO). It has been proposed that inhibition of its activity arrests the chlorophyll breakdown and leads to the "stay-green" phenotype. We searched the Arabidopsis genome for a possible PaO-encoding gene and hypothesized that it has homology to known iron-containing Rieske-type monooxygenase sequences. We identified three such open reading frames, Tic55, ACD1 and ACD1-like. We produced transgenic Arabidopsis plants which expressed antisense RNA as a method to inhibit the expression of these genes. The appearance of these antisense plants were indistinguishable from that of the wild type under illumination. However, after they were kept under darkness for 5 d and again illuminated, the leaves of the antisense ACD1 plants (AsACD1) were bleached. Leaves of AsACD1 accumulated 387 nmol (g FW)(-1) pheophorbide a which corresponded to 60% of chlorophyll a degraded. The rate of decrease in chlorophyll a was not influenced in senesced AsACD1 leaves. These results demonstrated that ACD1 is involved in PaO activity, and its inhibition led to photooxidative destruction of the cell instead of the "stay-green" phenotype.

Chlorophyll b inhibits the formation of photosystem I trimer in Synechocystis sp. PCC6803.

Chlorophyllide a oxygenase (CAO) catalyzes two-step oxygenation reactions and converts chlorophyllide a to chlorophyllide b. When CAO was introduced into the Synechocystis sp. PCC6803 genome, chlorophyll b was synthesized and incorporated into P700-chlorophyll a-protein complexes. Curve analysis of photosystem I particles showed that Ca687 was decreased with a concomitant increase in Cb652 suggesting that chlorophyll b was incorporated into Ca687-binding sites. When the level of chlorophyll b was high, Ca704, which is known as red chlorophyll and photosystem I trimers were decreased. Formation of photosystem I trimers is discussed in relation to red chlorophyll and chlorophyll b accumulation.