Tomoaki Nishiyama

Comparative genomics of Physcomitrella patens gametophytic transcriptome and Arabidopsis thaliana: Implication for land plant evolution

The mosses and flowering plants diverged >400 million years ago. The mosses have haploid-dominant life cycles, whereas the flowering plants are diploid-dominant. The common ancestors of land plants have been inferred to be haploid-dominant, suggesting that genes used in the diploid body of flowering plants were recruited from the genes used in the haploid body of the ancestors during the evolution of land plants. To assess this evolutionary hypothesis, we constructed an EST library of the moss Physcomitrella patens, and compared the moss transcriptome to the genome of Arabidopsis thaliana. We constructed full-length enriched cDNA libraries from auxin-treated, cytokinin-treated, and untreated gametophytes of P. patens, and sequenced both ends of >40,000 clones. These data, together with the mRNA sequences in the public databases, were assembled into 15,883 putative transcripts. Sequence comparisons of A. thaliana and P. patens showed that at least 66% of the A. thaliana genes had homologues in P. patens. Comparison of the P. patens putative transcripts with all known proteins, revealed 9,907 putative transcripts with high levels of similarity to vascular plant genes, and 850 putative transcripts with high levels of similarity to other organisms. The haploid transcriptome of P. patens appears to be quite similar to the A. thaliana genome, supporting the evolutionary hypothesis. Our study also revealed that a number of genes are moss specific and were lost in the flowering plant lineage.

The GID1-Mediated Gibberellin Perception Mechanism Is Conserved in the Lycophyte Selaginella moellendorffii but Not in the Bryophyte Physcomitrella patens

In rice (Oryza sativa) and Arabidopsis thaliana, gibberellin (GA) signaling is mediated by GIBBERELLIN-INSENSITIVE DWARF1 (GID1) and DELLA proteins in collaboration with a GA-specific F-box protein. To explore when plants evolved the ability to perceive GA by the GID1/DELLA pathway, we examined these GA signaling components in the lycophyte Selaginella moellendorffii and the bryophyte Physcomitrella patens. An in silico search identified several homologs of GID1, DELLA, and GID2, a GA-specific F-box protein in rice, in both species. Sm GID1a and Sm GID1b, GID1 proteins from S. moellendorffii, showed GA binding activity in vitro and interacted with DELLA proteins from S. moellendorffii in a GA-dependent manner in yeast. Introduction of constitutively expressed Sm GID1a, Sm G1D1b, and Sm GID2a transgenes rescued the dwarf phenotype of rice gid1 and gid2 mutants. Furthermore, treatment with GA4, a major GA in S. moellendorffii, caused downregulation of Sm GID1b, Sm GA20 oxidase, and Sm GA3 oxidase and degradation of the Sm DELLA1 protein. These results demonstrate that the homologs of GID1, DELLA, and GID2 work in a similar manner in S. moellendorffii and in flowering plants. Biochemical studies revealed that Sm GID1s have different GA binding properties from GID1s in flowering plants. No evidence was found for the functional conservation of these genes in P. patens, indicating that GID1/DELLA-mediated GA signaling, if present, differs from that in vascular plants. Our results suggest that GID1/DELLA-mediated GA signaling appeared after the divergence of vascular plants from the moss lineage.

Reannotation and extended community resources for the genome of the non-seed plant Physcomitrella patens provide insights into the evolution of plant gene structures and functions.

BackgroundThe moss Physcomitrella patens as a model species provides an important reference for early-diverging lineages of plants and the release of the genome in 2008 opened the doors to genome-wide studies. The usability of a reference genome greatly depends on the quality of the annotation and the availability of centralized community resources. Therefore, in the light of accumulating evidence for missing genes, fragmentary gene structures, false annotations and a low rate of functional annotations on the original release, we decided to improve the moss genome annotation.ResultsHere, we report the complete moss genome re-annotation (designated V1.6) incorporating the increased transcript availability from a multitude of developmental stages and tissue types. We demonstrate the utility of the improved P. patens genome annotation for comparative genomics and new extensions to the cosmoss.org resource as a central repository for this plant “flagship” genome. The structural annotation of 32,275 protein-coding genes results in 8387 additional loci including 1456 loci with known protein domains or homologs in Plantae. This is the first release to include information on transcript isoforms, suggesting alternative splicing events for at least 10.8% of the loci. Furthermore, this release now also provides information on non-protein-coding loci. Functional annotations were improved regarding quality and coverage, resulting in 58% annotated loci (previously: 41%) that comprise also 7200 additional loci with GO annotations. Access and manual curation of the functional and structural genome annotation is provided via the http://www.cosmoss.org model organism database.ConclusionsComparative analysis of gene structure evolution along the green plant lineage provides novel insights, such as a comparatively high number of loci with 5’-UTR introns in the moss. Comparative analysis of functional annotations reveals expansions of moss house-keeping and metabolic genes and further possibly adaptive, lineage-specific expansions and gains including at least 13% orphan genes.

TCC: an R package for comparing tag count data with robust normalization strategies

BackgroundDifferential expression analysis based on “next-generation” sequencing technologies is a fundamental means of studying RNA expression. We recently developed a multi-step normalization method (called TbT) for two-group RNA-seq data with replicates and demonstrated that the statistical methods available in four R packages (edgeR, DESeq, baySeq, and NBPSeq) together with TbT can produce a well-ranked gene list in which true differentially expressed genes (DEGs) are top-ranked and non-DEGs are bottom ranked. However, the advantages of the current TbT method come at the cost of a huge computation time. Moreover, the R packages did not have normalization methods based on such a multi-step strategy.ResultsTCC (an acronym for Tag Count Comparison) is an R package that provides a series of functions for differential expression analysis of tag count data. The package incorporates multi-step normalization methods, whose strategy is to remove potential DEGs before performing the data normalization. The normalization function based on this DEG elimination strategy (DEGES) includes (i) the original TbT method based on DEGES for two-group data with or without replicates, (ii) much faster methods for two-group data with or without replicates, and (iii) methods for multi-group comparison. TCC provides a simple unified interface to perform such analyses with combinations of functions provided by edgeR, DESeq, and baySeq. Additionally, a function for generating simulation data under various conditions and alternative DEGES procedures consisting of functions in the existing packages are provided. Bioinformatics scientists can use TCC to evaluate their methods, and biologists familiar with other R packages can easily learn what is done in TCC.ConclusionDEGES in TCC is essential for accurate normalization of tag count data, especially when up- and down-regulated DEGs in one of the samples are extremely biased in their number. TCC is useful for analyzing tag count data in various scenarios ranging from unbiased to extremely biased differential expression. TCC is available at http://www.iu.a.u-tokyo.ac.jp/~kadota/TCC/ and will appear in Bioconductor (http://bioconductor.org/) from ver. 2.13.

Involvement of auxin and a homeodomain-leucine zipper I gene in rhizoid development of the moss Physcomitrella patens.

Differentiation of epidermal cells is important for plants because they are in direct contact with the environment. Rhizoids are multicellular filaments that develop from the epidermis in a wide range of plants, including pteridophytes, bryophytes, and green algae; they have similar functions to root hairs in vascular plants in that they support the plant body and are involved in water and nutrient absorption. In this study, we examined mechanisms underlying rhizoid development in the moss, Physcomitrella patens, which is the only land plant in which high-frequency gene targeting is possible. We found that rhizoid development can be split into two processes: determination and differentiation. Two types of rhizoids with distinct developmental patterns (basal and mid-stem rhizoids) were recognized. The development of basal rhizoids from epidermal cells was induced by exogenous auxin, while that of mid-stem rhizoids required an unknown factor in addition to exogenous auxin. Once an epidermal cell had acquired a rhizoid initial cell fate, expression of the homeodomain-leucine zipper I gene Pphb7 was induced. Analysis of Pphb7 disruptant lines showed that Pphb7 affects the induction of pigmentation and the increase in the number and size of chloroplasts, but not the position or number of rhizoids. This is the first report on the involvement of a homeodomain-leucine zipper I gene in epidermal cell differentiation.

Class 1 KNOX genes are not involved in shoot development in the moss Physcomitrella patens but do function in sporophyte development

SUMMARY Although the number and form of metazoan organs are determined in the embryo, plants continuously form organs via pluripotent stem cells contained within the meristem. Flowering plants have an indeterminate meristem in their diploid generation, whereas the common ancestor of land plants is inferred to have formed an indeterminate meristem in its haploid generation, as observed in the extant basal land plants, bryophytes, including mosses. It is hypothesized that the underlying gene networks for the diploid meristem were initially present in the haploid generation of the basal land plants and were eventually co‐opted for expression in the diploid generation. In flowering plants, the class 1 KNOTTED1‐LIKE HOMEOBOX (KNOX) transcription factors are essential for the function of the indeterminate apical meristem. Here, we show that the class 1 KNOX orthologs function in the diploid organ, with determinate growth in the moss Physcomitrella patens, but do not function in the haploid indeterminate meristem. We propose that the genetic networks governing the indeterminate meristem in land plants are variable, and the networks governing the diploid indeterminate meristem with the class 1 KNOX genes likely evolved de novo in the flowering plant lineage.

A polycomb repressive complex 2 gene regulates apogamy and gives evolutionary insights into early land plant evolution

Land plants have distinct developmental programs in haploid (gametophyte) and diploid (sporophyte) generations. Although usually the two programs strictly alternate at fertilization and meiosis, one program can be induced during the other program. In a process called apogamy, cells of the gametophyte other than the egg cell initiate sporophyte development. Here, we report for the moss Physcomitrella patens that apogamy resulted from deletion of the gene orthologous to the Arabidopsis thaliana CURLY LEAF (PpCLF), which encodes a component of polycomb repressive complex 2 (PRC2). In the deletion lines, a gametophytic vegetative cell frequently gave rise to a sporophyte-like body. This body grew indeterminately from an apical cell with the character of a sporophytic pluripotent stem cell but did not form a sporangium. Furthermore, with continued culture, the sporophyte-like body branched. Sporophyte branching is almost unknown among extant bryophytes. When PpCLF was expressed in the deletion lines once the sporophyte-like bodies had formed, pluripotent stem cell activity was arrested and a sporangium-like organ formed. Supported by the observed pattern of PpCLF expression, these results demonstrate that, in the gametophyte, PpCLF represses initiation of a sporophytic pluripotent stem cell and, in the sporophyte, represses that stem cell activity and induces reproductive organ development. In land plants, branching, along with indeterminate apical growth and delayed initiation of spore-bearing reproductive organs, were conspicuous innovations for the evolution of a dominant sporophyte plant body. Our study provides insights into the role of PRC2 gene regulation for sustaining evolutionary innovation in land plants.

Algal ancestor of land plants was preadapted for symbiosis

Significance Colonization of land by plants was a critical event for the emergence of extant ecosystems. The innovations that allowed the algal ancestor of land plants to succeed in such a transition remain unknown. Beneficial interaction with symbiotic fungi has been proposed as one of these innovations. Here we show that the genes required for this interaction appeared in a stepwise manner: Some evolved before the colonization of land by plants and others first appeared in land plants. We thus propose that the algal ancestor of land plants was preadapted for interaction with beneficial fungi and employed these gene networks to colonize land successfully. Colonization of land by plants was a major transition on Earth, but the developmental and genetic innovations required for this transition remain unknown. Physiological studies and the fossil record strongly suggest that the ability of the first land plants to form symbiotic associations with beneficial fungi was one of these critical innovations. In angiosperms, genes required for the perception and transduction of diffusible fungal signals for root colonization and for nutrient exchange have been characterized. However, the origin of these genes and their potential correlation with land colonization remain elusive. A comprehensive phylogenetic analysis of 259 transcriptomes and 10 green algal and basal land plant genomes, coupled with the characterization of the evolutionary path leading to the appearance of a key regulator, a calcium- and calmodulin-dependent protein kinase, showed that the symbiotic signaling pathway predated the first land plants. In contrast, downstream genes required for root colonization and their specific expression pattern probably appeared subsequent to the colonization of land. We conclude that the most recent common ancestor of extant land plants and green algae was preadapted for symbiotic associations. Subsequent improvement of this precursor stage in early land plants through rounds of gene duplication led to the acquisition of additional pathways and the ability to form a fully functional arbuscular mycorrhizal symbiosis.

Evolution of MADS-Box Gene Induction by FLO/LFY Genes

Abstract. Some MADS-box genes function as floral homeotic genes. The Arabidopsis LFY gene is a positive regulator of floral homeotic genes, and homologs of the FLO/LFY gene family in other angiosperms and gymnosperms are likely to have a similar function. To investigate the origin of the floral homeotic gene regulatory cascade involving the FLO/LFY gene, FLO/LFY homologs were cloned from a leptosporangiate fern (Ceratopteris richardii), two eusporangiate ferns (Angiopteris lygodiifolia and Botrychium multifidum var. robustum), three fern allies (Psilotum nudum, Equisetum arvense, and Isoetes asiatica), and a moss (Physcomitrella patens). The FLO/LFY gene phylogenetic tree indicates that both duplication and loss of FLO/LFY homologs occurred during the course of vascular plant evolution. The expression patterns of the Ceratopteris LFY genes (CrLFY1 and 2) were assessed. CrLFY1 expression was prominent in tissues including shoot tips and circinate reproductive leaves, but very weak in other tissues examined. Expression of CrLFY2 was also prominent in tissues, including shoot tips and circinate reproductive leaves. These patterns of expression are dissimilar to that of any Ceratopteris MADS-box gene previously reported, suggesting that the induction of MADS-box genes by FLO/LFY is not established at the stage of ferns.

A normalization strategy for comparing tag count data

BackgroundHigh-throughput sequencing, such as ribonucleic acid sequencing (RNA-seq) and chromatin immunoprecipitation sequencing (ChIP-seq) analyses, enables various features of organisms to be compared through tag counts. Recent studies have demonstrated that the normalization step for RNA-seq data is critical for a more accurate subsequent analysis of differential gene expression. Development of a more robust normalization method is desirable for identifying the true difference in tag count data.ResultsWe describe a strategy for normalizing tag count data, focusing on RNA-seq. The key concept is to remove data assigned as potential differentially expressed genes (DEGs) before calculating the normalization factor. Several R packages for identifying DEGs are currently available, and each package uses its own normalization method and gene ranking algorithm. We compared a total of eight package combinations: four R packages (edgeR, DESeq, baySeq, and NBPSeq) with their default normalization settings and with our normalization strategy. Many synthetic datasets under various scenarios were evaluated on the basis of the area under the curve (AUC) as a measure for both sensitivity and specificity. We found that packages using our strategy in the data normalization step overall performed well. This result was also observed for a real experimental dataset.ConclusionOur results showed that the elimination of potential DEGs is essential for more accurate normalization of RNA-seq data. The concept of this normalization strategy can widely be applied to other types of tag count data and to microarray data.