The liverwort Marchantia polymorpha is equipped with a wide range of molecular and genetic tools and resources that have led to its wide use to explore the evo-devo aspects of land plants. Although its diverse transcriptome data are rapidly accumulating, there is no extensive yet user-friendly tool to exploit such a compilation of data and to summarize results with the latest annotations. Here, we have developed a web-based suite of tools, MarpolBase Expression (MBEX, https://marchantia.info/mbex/), where users can visualize gene expression profiles, identify differentially expressed genes, perform co-expression and functional enrichment analyses and summarize their comprehensive output in various portable formats. Using oil body biogenesis as an example, we demonstrated that the results generated by MBEX were consistent with the published experimental evidence and also revealed a novel transcriptional network in this process. MBEX should facilitate the exploration and discovery of the genetic and functional networks behind various biological processes in M. polymorpha and promote our understanding of the evolution of land plants.

Secondary loss of photosynthesis is observed across almost all plastid-bearing branches of the eukaryotic tree of life. However, genome-based insights into the transition from a phototroph into a secondary heterotroph have so far only been revealed for parasitic species. Free-living organisms can yield unique insights into the evolutionary consequence of the loss of photosynthesis, as the parasitic lifestyle requires specific adaptations to host environments. Here, we report on the diploid genome of the free-living diatomNitzschia putrida(35 Mbp), a nonphotosynthetic osmotroph whose photosynthetic relatives contribute ca. 40% of net oceanic primary production. Comparative analyses with photosynthetic diatoms and heterotrophic algae with parasitic lifestyle revealed that a combination of gene loss, the accumulation of genes involved in organic carbon degradation, a unique secretome, and the rapid divergence of conserved gene families involved in cell wall and extracellular metabolism appear to have facilitated the lifestyle of a free-living secondary heterotroph.

The liverwort Marchantia polymorpha is a key model organism for understanding land plant evolution, development, and gene regulation. To support the growing demand for high-quality genomic resources, we present MarpolBase, a comprehensive and integrated genome database that hosts newly assembled, high-accuracy reference genomes for both the male Tak-1 and female Tak-2 accessions, designated as ver. 7.1 reference genomes. These new assemblies, generated using PacBio HiFi long-read sequencing, represent nearly telomere-to-telomere chromosome-level genomes, with improvements in assembly continuity, annotation accuracy, and structural resolution—especially for repeat-rich regions and sex chromosomes. MarpolBase offers not only access to genome sequences and gene annotations but also provides a unified platform for data exploration, comparative analysis, and community-driven gene nomenclature for M. polymorpha. It includes keyword-searchable gene pages with structural and functional annotations, expression data integration, genome browser visualization, and online analytical and utility tools. By unifying genome assembly, annotation, nomenclature, and analysis tools in a single platform, MarpolBase serves as a central resource for functional genomics and evolutionary studies in M. polymorpha, and a model for future plant genome databases. The genomic resources of MarpolBase are freely available at https://marchantia.info.

The liverwort Marchantia polymorpha is a key model organism for understanding land plant evolution, development, and gene regulation. To support the growing demand for high-quality genomic resources, we present MarpolBase, a comprehensive and integrated genome database that hosts newly assembled, high-accuracy reference genomes for both the male Tak-1 and female Tak-2 accessions, designated as ver. 7.1 reference genomes. These new assemblies, generated using PacBio HiFi long-read sequencing, represent nearly telomere-to-telomere chromosome-level genomes, with improvements in assembly continuity, annotation accuracy, and structural resolution, especially for repeat-rich regions and sex chromosomes. MarpolBase offers not only access to genome sequences and gene annotations but also provides a unified platform for data exploration, comparative analysis, and community-driven gene nomenclature. It includes keyword-searchable gene pages with structural and functional annotations, expression data integration, genome browser visualization, and online analytical and utility tools. By unifying genome assembly, annotation, nomenclature, and analysis tools in a single platform, MarpolBase serves as a central resource for functional genomics and evolutionary studies in M. polymorpha, and a model for future plant genome databases.

Since the development of transcriptome analysis systems, many expression evolution studies characterized evolutionary forces acting on gene expression, without explicit discrimination between global expression differences and tissue specific expression differences. However, different types of gene expression alteration should have different effects on an organism, the evolutionary forces that act on them might be different, and different types of genes might show different types of differential expression between species. To confirm this, we studied differentially expressed (DE) genes among closely related groups that have extensive gene expression atlases, and clarified characteristics of different types of DE genes including the identification of regulating loci for differential expression using expression quantitative loci (eQTL) analysis data. We detected differentially expressed (DE) genes between rice subspecies in five homologous tissues that were verified using japonica and indica transcriptome atlases in public databases. Using the transcriptome atlases, we classified DE genes into two types, global DE genes and changed-tissues DE genes. Global type DE genes were not expressed in any tissues in the atlas of one subspecies, however changed-tissues type DE genes were expressed in both subspecies with different tissue specificity. For the five tissues in the two japonica-indica combinations, 4.6 ± 0.8 and 5.9 ± 1.5 % of highly expressed genes were global and changed-tissues DE genes, respectively. Changed-tissues DE genes varied in number between tissues, increasing linearly with the abundance of tissue specifically expressed genes in the tissue. Molecular evolution of global DE genes was rapid, unlike that of changed-tissues DE genes. Based on gene ontology, global and changed-tissues DE genes were different, having no common GO terms. Expression differences of most global DE genes were regulated by cis-eQTLs. Expression evolution of changed-tissues DE genes was rapid in tissue specifically expressed genes and those rapidly evolved changed-tissues DE genes were regulated not by cis-eQTLs, but by complicated trans-eQTLs. Global DE genes and changed-tissues DE genes had contrasting characteristics. The two contrasting types of DE genes provide possible explanations for the previous controversial conclusions about the relationships between molecular evolution and expression evolution of genes in different species, and the relationship between expression breadth and expression conservation in evolution.

Genome annotation is critically important data that can support research. Draft genome annotations cover representative genes; however, they often do not include genes that are expressed only in limited tissues and stages, or genes with low expression levels. Neuropeptides are responsible for regulation of various physiological and biological processes. A recent study disclosed the genome draft of the two-spotted cricket Gryllus bimaculatus, which was utilized to understand the intriguing physiology and biology of crickets. Thus far, only two of the nine reported neuropeptides in G. bimaculatus were annotated in the draft genome. Even though de novo assembly using transcriptomic analyses can comprehensively identify neuropeptides, this method does not follow those annotations on the genome locus. In this study, we performed the annotations based on the reference mapping, de novo transcriptome assembly, and manual curation. Consequently, we identified 41 neuropeptides out of 43 neuropeptides, which were reported in the insects. Further, 32 of the identified neuropeptides on the genomic loci in G. bimaculatus were annotated. The present annotation methods can be applicable for the neuropeptide annotation of other insects. Furthermore, the methods will help to generate useful infrastructures for studies relevant to neuropeptides.