Chromosome-scale genome assembly and annotation of Cotoneaster glaucophyllus

Cotoneaster glaucophyllus is a semi-evergreen plant that blossoms in late summer, producing dense, attractive, fragrant white flowers with significant ornamental and ecological value. Here, a chromosome-scale genome assembly was obtained by integrating PacBio and Illumina sequencing data with the aid of Hi-C technology. The genome assembly was 563.3 Mb in length, with contig N50 and scaffold N50 values of ~6 Mb and ~31 Mb, respectively. Most (95.59%) of the sequences were anchored onto 17 pseudochromosomes (538.4 Mb). We predicted 35,856 protein-coding genes, 1,401 miRNAs, 655 tRNAs, 425 rRNAs, and 795 snRNAs. The functions of 34,967 genes (97.52%) were predicted. The availability of this chromosome-level genome will provide valuable resources for molecular studies of this species, facilitating future research on speciation, functional genomics, and comparative genomics within the Rosaceae family.

Chromosome-level genome assembly of the diploid oat species Avena longiglumis

Diploid wild oat Avena longiglumis has nutritional and adaptive traits which are valuable for common oat (A. sativa) breeding. The combination of Illumina, Nanopore and Hi-C data allowed us to assemble a high-quality chromosome-level genome of A. longiglumis (ALO), evidenced by contig N50 of 12.68 Mb with 99% BUSCO completeness for the assembly size of 3,960.97 Mb. A total of 40,845 protein-coding genes were annotated. The assembled genome was composed of 87.04% repetitive DNA sequences. Dotplots of the genome assembly (PI657387) with two published ALO genomes were compared to indicate the conservation of gene order and equal expansion of all syntenic blocks among three genome assemblies. Two recent whole-genome duplication events were characterized in genomes of diploid Avena species. These findings provide new knowledge for the genomic features of A. longiglumis, give information about the species diversity, and will accelerate the functional genomics and breeding studies in oat and related cereal crops.

A chromosome-level genome assembly for the amphibious plant Rorippa aquatica reveals its allotetraploid origin and mechanisms of heterophylly upon submergence

The ability to respond to varying environments is crucial for sessile organisms such as plants. The amphibious plant Rorippa aquatica exhibits a striking type of phenotypic plasticity known as heterophylly, a phenomenon in which leaf form is altered in response to environmental factors. However, the underlying molecular mechanisms of heterophylly are yet to be fully understood. To uncover the genetic basis and analyze the evolutionary processes driving heterophylly in R. aquatica, we assembled the chromosome-level genome of the species. Comparative chromosome painting and chromosomal genomics revealed that allopolyploidization and subsequent post-polyploid descending dysploidy occurred during the speciation of R. aquatica. Based on the obtained genomic data, the transcriptome analyses revealed that ethylene signaling plays a central role in regulating heterophylly under submerged conditions, with blue light signaling acting as an attenuator of ethylene signal. The assembled R. aquatica reference genome provides insights into the molecular mechanisms and evolution of heterophylly.

Phylogenomic profiles of whole-genome duplications in Poaceae and landscape of differential duplicate retention and losses among major Poaceae lineages

Poaceae members shared a whole-genome duplication called rho. However, little is known about the evolutionary pattern of the rho-derived duplicates among Poaceae lineages and implications in adaptive evolution. Here we present phylogenomic/phylotranscriptomic analyses of 363 grasses covering all 12 subfamilies and report nine previously unknown whole-genome duplications. Furthermore, duplications from a single whole-genome duplication were mapped to multiple nodes on the species phylogeny; a whole-genome duplication was likely shared by woody bamboos with possible gene flow from herbaceous bamboos; and recent paralogues of a tetraploid Oryza are implicated in tolerance of seawater submergence. Moreover, rho duplicates showing differential retention among subfamilies include those with functions in environmental adaptations or morphogenesis, including ACOT for aquatic environments (Oryzoideae), CK2β for cold responses (Pooideae), SPIRAL1 for rapid cell elongation (Bambusoideae), and PAI1 for drought/cold responses (Panicoideae). This study presents a Poaceae whole-genome duplication profile with evidence for multiple evolutionary mechanisms that contribute to gene retention and losses.

Transcriptome profiling of aerial and subterranean peanut pod development

Peanut (Arachis hypogaea) showcases geocarpic behavior, transitioning from aerial flowering to subterranean seed development. We recently obtained an atavistic variant of this species, capable of producing aerial and subterranean pods on a single plant. Notably, although these pod types share similar vigor levels, they exhibit distinct differences in their physical aspects, such as pod size, color, and shell thickness. We constructed 63 RNA-sequencing datasets, comprising three biological replicates for each of 21 distinct tissues spanning six developmental stages for both pod types, providing a rich tapestry of the pod development process. This comprehensive analysis yielded an impressive 409.36 Gb of clean bases, facilitating the detection of 42,401 expressed genes. By comparing the transcriptomic data of the aerial and subterranean pods, we identified many differentially expressed genes (DEGs), highlighting their distinct developmental pathways. By providing a detailed workflow from the initial sampling to the final DEGs, this study serves as an important resource, paving the way for future research into peanut pod development and aiding transcriptome-based expression profiling and candidate gene identification.

Three de novo assembled wild cacao genomes from the Upper Amazon

Theobroma cacao, the chocolate tree, is indigenous to the Amazon basin, the greatest biodiversity hotspot on earth. Recent advancement in plant genomics highlights the importance of de novo sequencing of multiple reference genomes to capture the genome diversity present in different cacao populations. In this study, three high-quality chromosome-level genomes of wild cacao were constructed, de novo assembled with HiFi long reads sequencing, and scaffolded using a reference-free strategy. These genomes represent the three most important genetic clusters of cacao trees from the Upper Amazon region. The three wild cacao genomes were compared with two reference genomes of domesticated cacao. The five cacao genetic clusters were inferred to have diverged in the early and middle Pleistocene period, approximately 1.83-0.69 million years ago. The results shown here serve as an example of understanding how the Amazonian biodiversity was developed. The three wild cacao genomes provide valuable resources for studying genetic diversity and advancing genetic improvement of this species.

Allopolyploid origin and diversification of the Hawaiian endemic mints

Island systems provide important contexts for studying processes underlying lineage migration, species diversification, and organismal extinction. The Hawaiian endemic mints (Lamiaceae family) are the second largest plant radiation on the isolated Hawaiian Islands. We generated a chromosome-scale reference genome for one Hawaiian species, Stenogyne calaminthoides, and resequenced 45 relatives, representing 34 species, to uncover the continental origins of this group and their subsequent diversification. We further resequenced 109 individuals of two Stenogyne species, and their purported hybrids, found high on the Mauna Kea volcano on the island of Hawai'i. The three distinct Hawaiian genera, Haplostachys, Phyllostegia, and Stenogyne, are nested inside a fourth genus, Stachys. We uncovered four independent polyploidy events within Stachys, including one allopolyploidy event underlying the Hawaiian mints and their direct western North American ancestors. While the Hawaiian taxa may have principally diversified by parapatry and drift in small and fragmented populations, localized admixture may have played an important role early in lineage diversification. Our genomic analyses provide a view into how organisms may have radiated on isolated island chains, settings that provided one of the principal natural laboratories for Darwin's thinking about the evolutionary process.

Multiple horizontal gene transfer events have shaped plant glycosyl hydrolase diversity and function

Plant glycosyl hydrolases (GHs) play a crucial role in selectively breaking down carbohydrates and glycoconjugates during various cellular processes, such as reserve mobilization, pathogen defense, and modification/disassembly of the cell wall. In this study, we examined the distribution of GH genes in the Archaeplastida supergroup, which encompasses red algae, glaucophytes, and green plants. We identified that the GH repertoire expanded from a few tens of genes in early archaeplastidians to over 400 genes in modern angiosperms, spanning 40 GH families in land plants. Our findings reveal that major evolutionary transitions were accompanied by significant changes in the GH repertoire. Specifically, we identified at least 23 GH families acquired by green plants through multiple horizontal gene transfer events, primarily from bacteria and fungi. We found a significant shift in the subcellular localization of GH activity during green plant evolution, with a marked increase in extracellular-targeted GH proteins associated with the diversification of plant cell wall polysaccharides and defense mechanisms against pathogens. In conclusion, our study sheds light on the macroevolutionary processes that have shaped the GH repertoire in plants, highlighting the acquisition of GH families through horizontal transfer and the role of GHs in plant adaptation and defense mechanisms.

Complete organelle genomes of Korean fir, Abies koreana and phylogenomics of the gymnosperm genus Abies using nuclear and cytoplasmic DNA sequence data

Abies koreana E.H.Wilson is an endangered evergreen coniferous tree that is native to high altitudes in South Korea and susceptible to the effects of climate change. Hybridization and reticulate evolution have been reported in the genus; therefore, multigene datasets from nuclear and cytoplasmic genomes are needed to better understand its evolutionary history. Using the Illumina NovaSeq 6000 and Oxford Nanopore Technologies (ONT) PromethION platforms, we generated complete mitochondrial (1,174,803 bp) and plastid (121,341 bp) genomes from A. koreana. The mitochondrial genome is highly dynamic, transitioning from cis- to trans-splicing and breaking conserved gene clusters. In the plastome, the ONT reads revealed two structural conformations of A. koreana. The short inverted repeats (1186 bp) of the A. koreana plastome are associated with different structural types. Transcriptomic sequencing revealed 1356 sites of C-to-U RNA editing in the 41 mitochondrial genes. Using A. koreana as a reference, we additionally produced nuclear and organelle genomic sequences from eight Abies species and generated multiple datasets for maximum likelihood and network analyses. Three sections (Balsamea, Momi, and Pseudopicea) were well grouped in the nuclear phylogeny, but the phylogenomic relationships showed conflicting signals in the mitochondrial and plastid genomes, indicating a complicated evolutionary history that may have included introgressive hybridization. The obtained data illustrate that phylogenomic analyses based on sequences from differently inherited organelle genomes have resulted in conflicting trees. Organelle capture, organelle genome recombination, and incomplete lineage sorting in an ancestral heteroplasmic individual can contribute to phylogenomic discordance. We provide strong support for the relationships within Abies and new insights into the phylogenomic complexity of this genus.

Dynamics of accessible chromatin regions and subgenome dominance in octoploid strawberry

Subgenome dominance has been reported in diverse allopolyploid species, where genes from one subgenome are preferentially retained and are more highly expressed than those from other subgenome(s). However, the molecular mechanisms responsible for subgenome dominance remain poorly understood. Here, we develop genome-wide map of accessible chromatin regions (ACRs) in cultivated strawberry (2n = 8x = 56, with A, B, C, D subgenomes). Each ACR is identified as an MNase hypersensitive site (MHS). We discover that the dominant subgenome A contains a greater number of total MHSs and MHS per gene than the submissive B/C/D subgenomes. Subgenome A suffers fewer losses of MHS-related DNA sequences and fewer MHS fragmentations caused by insertions of transposable elements. We also discover that genes and MHSs related to stress response have been preferentially retained in subgenome A. We conclude that preservation of genes and their cognate ACRs, especially those related to stress responses, play a major role in the establishment of subgenome dominance in octoploid strawberry.

Chitosan combined with humic applications during sensitive growth stages to drought improves nutritional status and water relations of sweet potato

The current decline in freshwater resources presents a significant global challenge to crop production, a situation expected to intensify with ongoing climate change. This underscores the need for extensive research to enhance crop yields under drought conditions, a priority for scientists given its vital role in global food security. Our study explores the effects of using humic and chitosan treatments to alleviate drought stress during critical growth phases and their impact on crop yield and water efficiency. We employed four different irrigation strategies: full irrigation, 70% irrigation at the early vine development stage, 70% irrigation during the storage root bulking stage, and 85% irrigation across both stages, complemented by full irrigation in other periods. The plants received either humic treatments through foliar spray or soil application, or chitosan foliar applications, with tap water serving as a control. Our findings highlight that the early vine development stage is particularly vulnerable to drought, with a 42.0% decrease in yield observed under such conditions. In normal growth scenarios, foliar application of humic substances significantly improved growth parameters, resulting in a substantial increase in yield and water efficiency by 66.9% and 68.4%, respectively, compared to the control treatment under full irrigation. For sweet potatoes irrigated with 70% water at the storage root bulking stage, ground application of humic substances outperformed both foliar applications of chitosan and humic in terms of yield results. The highest tuber yield and water efficiency were attained by combining chitosan and humic ground applications, regardless of whether 70% irrigation was used at the storage root bulking stage or 85% irrigation during both the early vine development and storage root bulking stages.

Disparity of cycad leaves dispels the living fossil metaphor

The living fossil metaphor is tightly linked with the cycads. This group of gymnosperms is supposed to be characterised by long-term morphological stasis, particularly after their peak of diversity and disparity in the Jurassic. However, no formal test of this hypothesis exists. Here, we use a recent phylogenetic framework and an improved character matrix to reconstruct the Disparity Through Time for cycad leaves using a Principal Coordinate Analysis and employing Pre-Ordination Ancestral State Reconstruction to test the impact of sampling on the results. Our analysis shows that the cycad leaf morsphospace expanded up to the present, with numerous shifts in its general positioning, independently of sampling biases. Moreover, they also show that Zamiaceae expanded rapidly in the Early Cretaceous and continued to expand up to the present, while now-extinct clades experienced a slow contraction from their peak in the Triassic. We also show that rates of evolution were constantly high up to the Early Cretaceous, and then experienced a slight decrease in the Paleogene, followed by a Neogene acceleration. These results show a much more dynamic history for cycads, and suggest that the 'living fossil' metaphor is actually a hindrance to our understanding of their macroevolution.

Evolutionary insights into strategy shifts for the safe and effective accumulation of ascorbate in plants

Plants accumulate high concentrations of ascorbate, commonly in their leaves, as a redox buffer. While ascorbate levels have increased during plant evolution, the mechanisms behind this phenomenon are unclear. Moreover, has the increase in ascorbate concentration been achieved without imposing any detrimental effects on the plants? In this review, we focus on potential transitions in two regulatory mechanisms related to ascorbate biosynthesis and the availability of cellular dehydroascorbate (DHA) during plant evolution. The first transition might be that the trigger for the transcriptional induction of VTC2, which encodes the rate-limiting enzyme in ascorbate biosynthesis, has shifted from oxidative stress (in green algae) to light/photosynthesis (in land plants), probably enabling the continuous accumulation of ascorbate under illumination. This could serve as a preventive system against the unpredictable occurrence of oxidative stress. The second transition might be that DHA-degrading enzymes, which protect cells from the highly reactive DHA in green algae and mosses, have been lost in ferns or flowering plants. Instead, flowering plants may have increased glutathione concentrations to reinforce the DHA reduction capacity, possibly allowing ascorbate accumulation and avoiding the toxicity of DHA. These potential transitions may have contributed to strategies for plants' safe and effective accumulation of ascorbate.

Chromosome-level genome assembly of Platycarya strobilacea

Platycarya strobilacea belongs to the walnut family (Juglandaceae), is commonly known as species endemic to East Asia, and is an ecologically important, wind pollinated, woody deciduous tree. To facilitate this ancient tree for the ecological value and conservation of this ancient tree, we report a new high-quality genome assembly of P. strobilacea. The genome size was 677.30 Mb, with a scaffold N50 size of 45,791,698 bp, and 98.43% of the assembly was anchored to 15 chromosomes. We annotated 32,246 protein-coding genes in the genome, of which 96.30% were functionally annotated in six databases. This new high-quality assembly of P. strobilacea provide valuable resource for the phylogenetic and evolutionary analysis of the walnut family and angiosperm.

Evolutionary and synteny analysis of HIS1, BADH2, GBSS1, and GBSS2 in rice: insights for effective introgression breeding strategies

The key genes BADH2, GBSS1, GBSS2, and HIS1 regulate the fragrance, starch synthesis, and herbicide resistance in rice. Although the molecular functions of four genes have been investigated in the Oryza sativa species, little is known regarding their evolutionary history in the Oryza genus. Here, we studied the evolution of four focal genes in 10 Oryza species using phylogenetic and syntenic approaches. The HIS1 family underwent several times of tandem duplication events in the Oryza species, resulting in copy number variation ranging from 2 to 7. At most one copy of BADH2, GBSS1, and GBSS2 orthologs were identified in each Oryza species, and gene loss events of BADH2 and GBSS2 were identified in three Oryza species. Gene transfer analysis proposed that the functional roles of GBSS1 and GBSS2 were developed in the Asian and African regions, respectively, and most allelic variations of BADH2 in japonica rice emerged after the divergence between the Asian and African rice groups. These results provide clues to determine the origin and evolution of the key genes in rice breeding as well as valuable information for molecular breeders and scientists to develop efficient strategies to simultaneously improve grain quality and yield potential in rice.

Genomic evidence for rediploidization and adaptive evolution following the whole-genome triplication

Whole-genome duplication (WGD), or polyploidy, events are widespread and significant in the evolutionary history of angiosperms. However, empirical evidence for rediploidization, the major process where polyploids give rise to diploid descendants, is still lacking at the genomic level. Here we present chromosome-scale genomes of the mangrove tree Sonneratia alba and the related inland plant Lagerstroemia speciosa. Their common ancestor has experienced a whole-genome triplication (WGT) approximately 64 million years ago coinciding with a period of dramatic global climate change. Sonneratia, adapting mangrove habitats, experienced extensive chromosome rearrangements post-WGT. We observe the WGT retentions display sequence and expression divergence, suggesting potential neo- and sub-functionalization. Strong selection acting on three-copy retentions indicates adaptive value in response to new environments. To elucidate the role of ploidy changes in genome evolution, we improve a model of the polyploidization-rediploidization process based on genomic evidence, contributing to the understanding of adaptive evolution during climate change.

Cepharanthine analogs mining and genomes of Stephania accelerate anti-coronavirus drug discovery

Cepharanthine is a secondary metabolite isolated from Stephania. It has been reported that it has anti-conronaviruses activities including severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). Here, we assemble three Stephania genomes (S. japonica, S. yunnanensis, and S. cepharantha), propose the cepharanthine biosynthetic pathway, and assess the antiviral potential of compounds involved in the pathway. Among the three genomes, S. japonica has a near telomere-to-telomere assembly with one remaining gap, and S. cepharantha and S. yunnanensis have chromosome-level assemblies. Following by biosynthetic gene mining and metabolomics analysis, we identify seven cepharanthine analogs that have broad-spectrum anti-coronavirus activities, including SARS-CoV-2, Guangxi pangolin-CoV (GX_P2V), swine acute diarrhoea syndrome coronavirus (SADS-CoV), and porcine epidemic diarrhea virus (PEDV). We also show that two other genera, Nelumbo and Thalictrum, can produce cepharanthine analogs, and thus have the potential for antiviral compound discovery. Results generated from this study could accelerate broad-spectrum anti-coronavirus drug discovery.

Repeated upslope biome shifts in Saxifraga during late-Cenozoic climate cooling

Mountains are among the most biodiverse places on Earth, and plant lineages that inhabit them have some of the highest speciation rates ever recorded. Plant diversity within the alpine zone - the elevation above which trees cannot grow-contributes significantly to overall diversity within mountain systems, but the origins of alpine plant diversity are poorly understood. Here, we quantify the processes that generate alpine plant diversity and their changing dynamics through time in Saxifraga (Saxifragaceae), an angiosperm genus that occurs predominantly in mountain systems. We present a time-calibrated molecular phylogenetic tree for the genus that is inferred from 329 low-copy nuclear loci and incorporates 73% (407) of known species. We show that upslope biome shifts into the alpine zone are considerably more prevalent than dispersal of alpine specialists between regions, and that the rate of upslope biome shifts increased markedly in the last 5 Myr, a timeframe concordant with a cooling and fluctuating climate that is likely to have increased the extent of the alpine zone. Furthermore, alpine zone specialists have lower speciation rates than generalists that occur inside and outside the alpine zone, and major speciation rate increases within Saxifraga significantly pre-date increased rates of upslope biome shifts. Specialisation to the alpine zone is not therefore associated with speciation rate increases. Taken together, this study presents a quantified and broad scale perspective of processes underpinning alpine plant diversity.

Plastid ancestors lacked a complete Entner-Doudoroff pathway, limiting plants to glycolysis and the pentose phosphate pathway

The Entner-Doudoroff (ED) pathway provides an alternative to glycolysis. It converts 6-phosphogluconate (6-PG) to glyceraldehyde-3-phosphate and pyruvate in two steps consisting of a dehydratase (EDD) and an aldolase (EDA). Here, we investigate its distribution and significance in higher plants and determine the ED pathway is restricted to prokaryotes due to the absence of EDD genes in eukaryotes. EDDs share a common origin with dihydroxy-acid dehydratases (DHADs) of the branched chain amino acid pathway (BCAA). Each dehydratase features strict substrate specificity. E. coli EDD dehydrates 6-PG to 2-keto-3-deoxy-6-phosphogluconate, while DHAD only dehydrates substrates from the BCAA pathway. Structural modeling identifies two divergent domains which account for their non-overlapping substrate affinities. Coupled enzyme assays confirm only EDD participates in the ED pathway. Plastid ancestors lacked EDD but transferred metabolically promiscuous EDA, which explains the absence of the ED pathway from the Viridiplantae and sporadic persistence of EDA genes across the plant kingdom.

Phylogenomic insights into the first multicellular streptophyte

Streptophytes are best known as the clade containing the teeming diversity of embryophytes (land plants).1,2,3,4 Next to embryophytes are however a range of freshwater and terrestrial algae that bear important information on the emergence of key traits of land plants. Among these, the Klebsormidiophyceae stand out. Thriving in diverse environments-from mundane (ubiquitous occurrence on tree barks and rocks) to extreme (from the Atacama Desert to the Antarctic)-Klebsormidiophyceae can exhibit filamentous body plans and display remarkable resilience as colonizers of terrestrial habitats.5,6 Currently, the lack of a robust phylogenetic framework for the Klebsormidiophyceae hampers our understanding of the evolutionary history of these key traits. Here, we conducted a phylogenomic analysis utilizing advanced models that can counteract systematic biases. We sequenced 24 new transcriptomes of Klebsormidiophyceae and combined them with 14 previously published genomic and transcriptomic datasets. Using an analysis built on 845 loci and sophisticated mixture models, we establish a phylogenomic framework, dividing the six distinct genera of Klebsormidiophyceae in a novel three-order system, with a deep divergence more than 830 million years ago. Our reconstructions of ancestral states suggest (1) an evolutionary history of multiple transitions between terrestrial-aquatic habitats, with stem Klebsormidiales having conquered land earlier than embryophytes, and (2) that the body plan of the last common ancestor of Klebsormidiophyceae was multicellular, with a high probability that it was filamentous whereas the sarcinoids and unicells in Klebsormidiophyceae are likely derived states. We provide evidence that the first multicellular streptophytes likely lived about a billion years ago.

Chromosome level genome assembly of endangered medicinal plant Anisodus tanguticus

Anisodus tanguticus is a medicinal herb that belongs to the Anisodus genus of the Solanaceae family. This endangered herb is mainly distributed in Qinghai-Tibet Plateau. In this study, we combined the Illumina short-read, Nanopore long-read and high-throughput chromosome conformation capture (Hi-C) sequencing technologies to de novo assemble the A. tanguticus genome. A high-quality chromosomal-level genome assembly was obtained with a genome size of 1.26 Gb and a contig N50 of 25.07 Mb. Of the draft genome sequences, 97.47% were anchored to 24 pseudochromosomes with a scaffold N50 of 51.28 Mb. In addition, 842.14 Mb of transposable elements occupying 66.70% of the genome assembly were identified and 44,252 protein-coding genes were predicted. The genome assembly of A. tanguticus will provide genetic repertoire to understand the adaptation strategy of Anisodus species in the plateau, which will further promote the conservation of endangered A. tanguticus resources.

Cryogenian Origins of Multicellularity in Archaeplastida

Earth was impacted by global glaciations during the Cryogenian (720 to 635 million years ago; Ma), events invoked to explain both the origins of multicellularity in Archaeplastida and radiation of the first land plants. However, the temporal relationship between these environmental and biological events is poorly established, due to a paucity of molecular and fossil data, precluding resolution of the phylogeny and timescale of archaeplastid evolution. We infer a time-calibrated phylogeny of early archaeplastid evolution based on a revised molecular dataset and reappraisal of the fossil record. Phylogenetic topology testing resolves deep archaeplastid relationships, identifying two clades of Viridiplantae and placing Bryopsidales as sister to the Chlorophyceae. Our molecular clock analysis infers an origin of Archaeplastida in the late-Paleoproterozoic to early-Mesoproterozoic (1712 to 1387 Ma). Ancestral state reconstruction of cytomorphological traits on this time-calibrated tree reveals many of the independent origins of multicellularity span the Cryogenian, consistent with the Cryogenian multicellularity hypothesis. Multicellular rhodophytes emerged 902 to 655 Ma while crown-Anydrophyta (Zygnematophyceae and Embryophyta) originated 796 to 671 Ma, broadly compatible with the Cryogenian plant terrestrialization hypothesis. Our analyses resolve the timetree of Archaeplastida with age estimates for ancestral multicellular archaeplastids coinciding with the Cryogenian, compatible with hypotheses that propose a role of Snowball Earth in plant evolution.

Induced responses contribute to rapid adaptation of Spirodela polyrhiza to herbivory by Lymnaea stagnalis

Herbivory-induced responses in plants are typical examples of phenotypic plasticity, and their evolution is thought to be driven by herbivory. However, direct evidence of the role of induced responses in plant adaptive evolution to herbivores is scarce. Here, we experimentally evolve populations of an aquatic plant (Spirodela polyrhiza, giant duckweed) and its native herbivore (Lymnaea stagnalis, freshwater snail), testing whether herbivory drives rapid adaptive evolution in plant populations using a combination of bioassays, pool-sequencing, metabolite analyses, and amplicon metagenomics. We show that snail herbivory drove rapid phenotypic changes, increased herbivory resistance, and altered genotype frequencies in the plant populations. Additional bioassays suggest that evolutionary changes of induced responses contributed to the rapid increase of plant resistance to herbivory. This study provides direct evidence that herbivory-induced responses in plants can be subjected to selection and have an adaptive role by increasing resistance to herbivores.

A chromosome-level genome assembly for Onobrychis viciifolia reveals gene copy number gain underlying enhanced proanthocyanidin biosynthesis

Sainfoin (Onobrychis viciifolia), which belongs to subfamily Papilionoideae of Leguminosae, is a vital perennial forage known as "holy hay" due to its high contents of crude proteins and proanthocyanidins (PAs, also called condensed tannins) that have various pharmacological properties in animal feed, such as alleviating rumen tympanic disease in ruminants. In this study, we select an autotetraploid common sainfoin (2n = 4x = 28) and report its high-quality chromosome-level genome assembly with 28 pseudochromosomes and four haplotypes (~1950.14 Mb, contig N50 = 10.91 Mb). The copy numbers of genes involved in PA biosynthesis in sainfoin are significantly greater than those in four selected Fabales species, namely, autotetraploid Medicago sativa and three other diploid species, Lotus japonicus, Medicago truncatula, and Glycine max. Furthermore, gene expansion is confirmed to be the key contributor to the increased expression of these genes and subsequent PA enhancement in sainfoin. Transcriptomic analyses reveal that the expression of genes involved in the PA biosynthesis pathway is significantly increased in the lines with high PA content compared to the lines with medium and low PA content. The sainfoin genome assembly will improve our understanding of leguminous genome evolution and biosynthesis of secondary metabolites in sainfoin.

RAF-like protein kinases mediate a deeply conserved, rapid auxin response

The plant-signaling molecule auxin triggers fast and slow cellular responses across land plants and algae. The nuclear auxin pathway mediates gene expression and controls growth and development in land plants, but this pathway is absent from algal sister groups. Several components of rapid responses have been identified in Arabidopsis, but it is unknown if these are part of a conserved mechanism. We recently identified a fast, proteome-wide phosphorylation response to auxin. Here, we show that this response occurs across 5 land plant and algal species and converges on a core group of shared targets. We found conserved rapid physiological responses to auxin in the same species and identified rapidly accelerated fibrosarcoma (RAF)-like protein kinases as central mediators of auxin-triggered phosphorylation across species. Genetic analysis connects this kinase to both auxin-triggered protein phosphorylation and rapid cellular response, thus identifying an ancient mechanism for fast auxin responses in the green lineage.

Hidden functional complexity in the flora of an early land ecosystem

The first land ecosystems were composed of organisms considered simple in nature, yet the morphological diversity of their flora was extraordinary. The biological significance of this diversity remains a mystery largely due to the absence of feasible study approaches. To study the functional biology of Early Devonian flora, we have reconstructed extinct plants from fossilised remains in silico. We explored the morphological diversity of sporangia in relation to their mechanical properties using finite element method. Our approach highlights the impact of sporangia morphology on spore dispersal and adaptation. We discovered previously unidentified innovations among early land plants, discussing how different species might have opted for different spore dispersal strategies. We present examples of convergent evolution for turgor pressure resistance, achieved by homogenisation of stress in spherical sporangia and by torquing force in Tortilicaulis-like specimens. In addition, we show a potential mechanism for stress-assisted sporangium rupture. Our study reveals the deceptive complexity of this seemingly simple group of organisms. We leveraged the quantitative nature of our approach and constructed a fitness landscape to understand the different ecological niches present in the Early Devonian Welsh Borderland flora. By connecting morphology to functional biology, these findings facilitate a deeper understanding of the diversity of early land plants and their place within their ecosystem.

A chromosome-scale assembly of the early-flowering Prunus campanulata and comparative genomics of cherries

Prunus campanulata is an important flowering cherry germplasm of high ornamental value. Given its early-flowering phenotypes, P. campanulata could be used for molecular breeding of ornamental species and fruit crops belonging to the subgenus Cerasus. Here, we report a chromosome-scale assembly of P. campanulata with a genome size of 282.6 Mb and a contig N50 length of 12.04 Mb. The genome contained 24,861 protein-coding genes, of which 24,749 genes (99.5%) were functionally annotated, and 148.20 Mb (52.4%) of the assembled sequences are repetitive sequences. A combination of genomic and population genomic analyses revealed a number of genes under positive selection or accelerated molecular evolution in P. campanulata. Our study provides a reliable genome resource, and lays a solid foundation for genetic improvement of flowering cherry germplasm.

RHO of plant signaling was established early in streptophyte evolution

The algal ancestors of land plants underwent a transition from a unicellular to a multicellular body plan.1 This transition likely took place early in streptophyte evolution, sometime after the divergence of the Chlorokybophyceae/Mesostigmatophyceae lineage, but before the divergence of the Klebsormidiophyceae lineage.2 How this transition was brought about is unknown; however, it was likely facilitated by the evolution of novel mechanisms to spatially regulate morphogenesis. In land plants, RHO of plant (ROP) signaling plays a conserved role in regulating polarized cell growth and cell division orientation to orchestrate morphogenesis.3,4,5,6,7,8 ROP constitutes a plant-specific subfamily of the RHO GTPases, which are more widely conserved throughout eukaryotes.9,10 Although the RHO family originated in early eukaryotes,11,12 how and when the ROP subfamily originated had remained elusive. Here, we demonstrate that ROP signaling was established early in the streptophyte lineage, sometime after the divergence of the Chlorokybophyceae/Mesostigmatophyceae lineage, but before the divergence of the Klebsormidiophyceae lineage. This period corresponds to when the unicellular-to-multicellular transition likely took place in the streptophytes. In addition to being critical for the complex morphogenesis of extant land plants, we speculate that ROP signaling contributed to morphological evolution in early streptophytes.

Haplotype-resolved chromosomal-level genome assembly of Buzhaye (Microcos paniculata)

Microcos paniculata is a shrub used traditionally as folk medicine and to make herbal teas. Previous research into this species has mainly focused on its chemical composition and medicinal value. However, the lack of a reference genome limits the study of the molecular mechanisms of active compounds in this species. Here, we assembled a haplotype-resolved chromosome-level genome of M. paniculata based on PacBio HiFi and Hi-C data. The assembly contains two haploid genomes with sizes 399.43 Mb and 393.10 Mb, with contig N50 lengths of 43.44 Mb and 30.17 Mb, respectively. About 99.93% of the assembled sequences could be anchored to 18 pseudo-chromosomes. Additionally, a total of 482 Mb repeat sequences were identified, accounting for 60.76% of the genome. A total of 49,439 protein-coding genes were identified, of which 48,979 (99%) were functionally annotated. This haplotype-resolved chromosome-level assembly and annotation of M. paniculata will serve as a valuable resource for investigating the biosynthesis and genetic basis of active compounds in this species, as well as advancing evolutionary phylogenomic studies in Malvales.

A chromosome-scale assembly of the quinoa genome provides insights into the structure and dynamics of its subgenomes

Quinoa (Chenopodium quinoa Willd.) is an allotetraploid seed crop with the potential to help address global food security concerns. Genomes have been assembled for four accessions of quinoa; however, all assemblies are fragmented and do not reflect known chromosome biology. Here, we use in vitro and in vivo Hi-C data to produce a chromosome-scale assembly of the Chilean accession PI 614886 (QQ74). The final assembly spans 1.326 Gb, of which 90.5% is assembled into 18 chromosome-scale scaffolds. The genome is annotated with 54,499 protein-coding genes, 96.9% of which are located on the 18 largest scaffolds. We also report an updated genome assembly for the B-genome diploid C. suecicum and use it, together with the A-genome diploid C. pallidicaule, to identify genomic rearrangements within the quinoa genome, including a large pericentromeric inversion representing 71.7% of chromosome Cq3B. Repetitive sequences comprise 65.2%, 48.6%, and 57.9% of the quinoa, C. pallidicaule, and C. suecicum genomes, respectively. Evidence suggests that the B subgenome is more dynamic and has expanded more than the A subgenome. These genomic resources will enable more accurate assessments of genome evolution within the Amaranthaceae and will facilitate future efforts to identify variation in genes underlying important agronomic traits in quinoa.

Cell geometry regulates tissue fracture

In vascular plants, the epidermal surfaces of leaves and flower petals often display cells with wavy geometries forming intricate jigsaw puzzle patterns. The prevalence and diversity of these complex epidermal patterns, originating from simple polyhedral progenitor cells, suggest adaptive significance. However, despite multiple efforts to explain the evolutionary drivers behind these geometrical features, compelling validation remains elusive. Employing a multidisciplinary approach that integrates microscopic and macroscopic fracture experiments with computational fracture mechanics, we demonstrate that wavy epidermal cells toughen the plants' protective skin. Through a multi-scale framework, we demonstrate that this energy-efficient patterning mechanism is universally applicable for toughening biological and synthetic materials. Our findings reveal a tunable structural-mechanical strategy employed in the microscale design of plants to protect them from deleterious surface fissures while facilitating and strategically directing beneficial ones. These findings hold implications for targeted plant breeding aimed at enhancing resilience in fluctuating environmental conditions. From an engineering perspective, our work highlights the sophisticated design principles the plant kingdom offers to inspire metamaterials.

Evolution of a plant growth-regulatory protein interaction specificity

Specific protein-protein interactions (PPIs) enable biological regulation. However, the evolution of PPI specificity is little understood. Here we trace the evolution of the land-plant growth-regulatory DELLA-SLY1/GID2 PPI, revealing progressive increase in specificity of affinity of SLY1/GID2 for a particular DELLA form. While early-diverging SLY1s display relatively broad-range DELLA affinity, later-diverging SLY1s tend towards increasingly stringent affinity for a specific DELLA A' form generated by the growth-promoting phytohormone gibberellin (GA). Our novel mutational strategy reveals amino acid substitutions contributing to the evolution of Arabidopsis thaliana SLY1 A' specificity, also showing that routes permitting reversion to broader affinity became increasingly constrained over evolutionary time. We suggest that progressive affinity narrowing may be an important evolutionary driver of PPI specificity and that increase in SLY1/GID2-DELLA specificity enabled the enhanced flexibility of plant physiological environmental adaptation conferred by the GA-DELLA growth-regulatory mechanism.

Landscape genomics reveals adaptive genetic differentiation driven by multiple environmental variables in naked barley on the Qinghai-Tibetan Plateau

Understanding the local adaptation of crops has long been a concern of evolutionary biologists and molecular ecologists. Identifying the adaptive genetic variability in the genome is crucial not only to provide insights into the genetic mechanism of local adaptation but also to explore the adaptation potential of crops. This study aimed to identify the climatic drivers of naked barley landraces and putative adaptive loci driving local adaptation on the Qinghai-Tibetan Plateau (QTP). To this end, a total of 157 diverse naked barley accessions were genotyped using the genotyping-by-sequencing approach, which yielded 3123 high-quality SNPs for population structure analysis and partial redundancy analysis, and 37,636 SNPs for outlier analysis. The population structure analysis indicated that naked barley landraces could be divided into four groups. We found that the genomic diversity of naked barley landraces could be partly traced back to the geographical and environmental diversity of the landscape. In total, 136 signatures associated with temperature, precipitation, and ultraviolet radiation were identified, of which 13 had pleiotropic effects. We mapped 447 genes, including a known gene HvSs1. Some genes involved in cold stress and regulation of flowering time were detected near eight signatures. Taken together, these results highlight the existence of putative adaptive loci in naked barley on QTP and thus improve our current understanding of the genetic basis of local adaptation.

Phylogenetic, population structure, and population demographic analyses reveal that Vicia sepium in Japan is native and not introduced

Vicia sepium (bush vetch) is a perennial legume widely distributed throughout the Eurasian continent. However, its distribution in Japan is limited to Mt. Ibuki and small parts of central and southern Hokkaido. Therefore, each Japanese V. sepium lineage has been considered to have been introduced separately from Europe. Here, we examined whether the species was introduced or not on the basis of cpDNA sequences and genome-wide SNPs from Japanese and overseas samples. Both the cpDNA haplotype network and the nuclear DNA phylogenetic tree showed that Japanese V. sepium is monophyletic. Furthermore, although the nuclear DNA phylogenetic tree also showed that each lineage is clearly monophyletic, genetic admixture of the genetic cluster dominated in the Hokkaido lineage was also detected in the Mt. Ibuki lineage. Population divergence analysis showed that the two lineages diverged during the last glacial period. The Mt. Ibuki lineage showed a sudden population decline 300-400 years ago, indicating that some anthropogenic activity might be involved, while the Hokkaido lineage showed a gradual population decline from 5000 years ago. Consequently, these two lineages show low current genetic diversity compared with overseas lineages. These results show that the Japanese V. sepium is not introduced but is native.

LHP1-mediated epigenetic buffering of subgenome diversity and defense responses confers genome plasticity and adaptability in allopolyploid wheat

Polyploidization is a major driver of genome diversification and environmental adaptation. However, the merger of different genomes may result in genomic conflicts, raising a major question regarding how genetic diversity is interpreted and regulated to enable environmental plasticity. By analyzing the genome-wide binding of 191 trans-factors in allopolyploid wheat, we identified like heterochromatin protein 1 (LHP1) as a master regulator of subgenome-diversified genes. Transcriptomic and epigenomic analyses of LHP1 mutants reveal its role in buffering the expression of subgenome-diversified defense genes by controlling H3K27me3 homeostasis. Stripe rust infection releases latent subgenomic variations by eliminating H3K27me3-related repression. The simultaneous inactivation of LHP1 homoeologs by CRISPR-Cas9 confers robust stripe rust resistance in wheat seedlings. The conditional repression of subgenome-diversified defenses ensures developmental plasticity to external changes, while also promoting neutral-to-non-neutral selection transitions and adaptive evolution. These findings establish an LHP1-mediated buffering system at the intersection of genotypes, environments, and phenotypes in polyploid wheat. Manipulating the epigenetic buffering capacity offers a tool to harness cryptic subgenomic variations for crop improvement.

A toxin-antidote system contributes to interspecific reproductive isolation in rice

Breakdown of reproductive isolation facilitates flow of useful trait genes into crop plants from their wild relatives. Hybrid sterility, a major form of reproductive isolation exists between cultivated rice (Oryza sativa) and wild rice (O. meridionalis, Mer). Here, we report the cloning of qHMS1, a quantitative trait locus controlling hybrid male sterility between these two species. Like qHMS7, another locus we cloned previously, qHMS1 encodes a toxin-antidote system, but differs in the encoded proteins, their evolutionary origin, and action time point during pollen development. In plants heterozygous at qHMS1, ~ 50% of pollens carrying qHMS1-D (an allele from cultivated rice) are selectively killed. In plants heterozygous at both qHMS1 and qHMS7, ~ 75% pollens without co-presence of qHMS1-Mer and qHMS7-D are selectively killed, indicating that the antidotes function in a toxin-dependent manner. Our results indicate that different toxin-antidote systems provide stacked reproductive isolation for maintaining species identity and shed light on breakdown of hybrid male sterility.

Taming the perils of photosynthesis by eukaryotes: constraints on endosymbiotic evolution in aquatic ecosystems

An ancestral eukaryote acquired photosynthesis by genetically integrating a cyanobacterial endosymbiont as the chloroplast. The chloroplast was then further integrated into many other eukaryotic lineages through secondary endosymbiotic events of unicellular eukaryotic algae. While photosynthesis enables autotrophy, it also generates reactive oxygen species that can cause oxidative stress. To mitigate the stress, photosynthetic eukaryotes employ various mechanisms, including regulating chloroplast light absorption and repairing or removing damaged chloroplasts by sensing light and photosynthetic status. Recent studies have shown that, besides algae and plants with innate chloroplasts, several lineages of numerous unicellular eukaryotes engage in acquired phototrophy by hosting algal endosymbionts or by transiently utilizing chloroplasts sequestrated from algal prey in aquatic ecosystems. In addition, it has become evident that unicellular organisms engaged in acquired phototrophy, as well as those that feed on algae, have also developed mechanisms to cope with photosynthetic oxidative stress. These mechanisms are limited but similar to those employed by algae and plants. Thus, there appear to be constraints on the evolution of those mechanisms, which likely began by incorporating photosynthetic cells before the establishment of chloroplasts by extending preexisting mechanisms to cope with oxidative stress originating from mitochondrial respiration and acquiring new mechanisms.

Discovery of isoflavone phytoalexins in wheat reveals an alternative route to isoflavonoid biosynthesis

Isoflavones are a group of phenolic compounds mostly restricted to plants of the legume family, where they mediate important interactions with plant-associated microbes, including in defense from pathogens and in nodulation. Their well-studied health promoting attributes have made them a prime target for metabolic engineering, both for bioproduction of isoflavones as high-value molecules, and in biofortification of food crops. A key gene in their biosynthesis, isoflavone synthase, was identified in legumes over two decades ago, but little is known about formation of isoflavones outside of this family. Here we identify a specialized wheat-specific isoflavone synthase, TaCYP71F53, which catalyzes a different reaction from the leguminous isoflavone synthases, thus revealing an alternative path to isoflavonoid biosynthesis and providing a non-transgenic route for engineering isoflavone production in wheat. TaCYP71F53 forms part of a biosynthetic gene cluster that produces a naringenin-derived O-methylated isoflavone, 5-hydroxy-2',4',7-trimethoxyisoflavone, triticein. Pathogen-induced production and in vitro antimicrobial activity of triticein suggest a defense-related role for this molecule in wheat. Genomic and metabolic analyses of wheat ancestral grasses further show that the triticein gene cluster was introduced into domesticated emmer wheat through natural hybridization ~9000 years ago, and encodes a pathogen-responsive metabolic pathway that is conserved in modern bread wheat varieties.

A high-quality chromosome-level wild rice genome of Oryza coarctata

Oryza coarctata (2n = 4X = 48, KKLL) is an allotetraploid, undomesticated relative of rice and the only species in the genus Oryza with tolerance to high salinity and submergence. Therefore, it contains important stress and tolerance genes/factors for rice. The initial draft genome published was limited by data and technical restrictions, leading to an incomplete and highly fragmented assembly. This study reports a new, highly contiguous chromosome-level genome assembly and annotation of O. coarctata. PacBio high-quality HiFi reads generated 460 contigs with a total length of 573.4 Mb and an N50 of 23.1 Mb, which were assembled into scaffolds with Hi-C data, anchoring 96.99% of the assembly onto 24 chromosomes. The genome assembly comprises 45,571 genes, and repetitive content contributes 25.5% of the genome. This study provides the novel identification of the KK and LL genome types of the genus Oryza, leading to valuable insights into rice genome evolution. The chromosome-level genome assembly of O. coarctata is a valuable resource for rice research and molecular breeding.

The pan-genome and local adaptation of Arabidopsis thaliana

Arabidopsis thaliana serves as a model species for investigating various aspects of plant biology. However, the contribution of genomic structural variations (SVs) and their associate genes to the local adaptation of this widely distribute species remains unclear. Here, we de novo assemble chromosome-level genomes of 32 A. thaliana ecotypes and determine that variable genes expand the gene pool in different ecotypes and thus assist local adaptation. We develop a graph-based pan-genome and identify 61,332 SVs that overlap with 18,883 genes, some of which are highly involved in ecological adaptation of this species. For instance, we observe a specific 332 bp insertion in the promoter region of the HPCA1 gene in the Tibet-0 ecotype that enhances gene expression, thereby promotes adaptation to alpine environments. These findings augment our understanding of the molecular mechanisms underlying the local adaptation of A. thaliana across diverse habitats.

Asymmetric Evolution of Protein Domains in the Leucine-Rich Repeat Receptor-Like Kinase Family of Plant Signaling Proteins

The coding sequences of developmental genes are expected to be deeply conserved, with cis-regulatory change driving the modulation of gene function. In contrast, proteins with roles in defense are expected to evolve rapidly, in molecular arms races with pathogens. However, some gene families include both developmental and defense genes. In these families, does the tempo and mode of evolution differ between genes with divergent functions, despite shared ancestry and structure? The leucine-rich repeat receptor-like kinase (LRR-RLKs) protein family includes members with roles in plant development and defense, thus providing an ideal system for answering this question. LRR-RLKs are receptors that traverse plasma membranes. LRR domains bind extracellular ligands; RLK domains initiate intracellular signaling cascades in response to ligand binding. In LRR-RLKs with roles in defense, LRR domains evolve faster than RLK domains. To determine whether this asymmetry extends to LRR-RLKs that function primarily in development, we assessed evolutionary rates and tested for selection acting on 11 subfamilies of LRR-RLKs, using deeply sampled protein trees. To assess functional evolution, we performed heterologous complementation assays in Arabidopsis thaliana (Arabidopsis). We found that the LRR domains of all tested LRR-RLK proteins evolved faster than their cognate RLK domains. All tested subfamilies of LRR-RLKs had strikingly similar patterns of molecular evolution, despite divergent functions. Heterologous transformation experiments revealed that multiple mechanisms likely contribute to the evolution of LRR-RLK function, including escape from adaptive conflict. Our results indicate specific and distinct evolutionary pressures acting on LRR versus RLK domains, despite diverse organismal roles for LRR-RLK proteins.

Interpreting shoot flammability experiments: from fire behaviour to plant evolution

This article is a Commentary on Cui et al. (2023), 240: 105–113.

Evolution of phenotypic disparity in the plant kingdom

The plant kingdom exhibits diverse bodyplans, from single-celled algae to complex multicellular land plants, but it is unclear how this phenotypic disparity was achieved. Here we show that the living divisions comprise discrete clusters within morphospace, separated largely by reproductive innovations, the extinction of evolutionary intermediates and lineage-specific evolution. Phenotypic complexity correlates not with disparity but with ploidy history, reflecting the role of genome duplication in plant macroevolution. Overall, the plant kingdom exhibits a pattern of episodically increasing disparity throughout its evolutionary history that mirrors the evolutionary floras and reflects ecological expansion facilitated by reproductive innovations. This pattern also parallels that seen in the animal and fungal kingdoms, suggesting a general pattern for the evolution of multicellular bodyplans.

Identification and evolution of nsLTPs in the root nodule nitrogen fixation clade and molecular response of Frankia to AgLTP24

Non-specific lipid transfer proteins (nsLTPs) are antimicrobial peptides, involved in several plant biological processes including root nodule nitrogen fixation (RNF). Nodulating plants belonging to the RNF clade establish symbiosis with the nitrogen-fixing bacteria rhizobia (legumes symbiosis model) and Frankia (actinorhizal symbiosis model) leading to root nodule formation. nsLTPs are involved in processes active in early step of symbiosis and functional nodule in both models. In legumes, nsLTPs have been shown to regulate symbiont entry, promote root cortex infection, membrane biosynthesis, and improve symbiosis efficiency. More recently, a nsLTP, AgLTP24 has been described in the context of actinorhizal symbiosis between Alnus glutinosa and Frankia alni ACN14a. AgLTP24 is secreted at an early step of symbiosis on the deformed root hairs and targets the symbiont in the nitrogen-fixing vesicles in functional nodules. nsLTPs are involved in RNF, but their functions and evolutionary history are still largely unknown. Numerous putative nsLTPs were found up-regulated in functional nodules compared to non-infected roots in different lineages within the RNF clade. Here, results highlight that nodulating plants that are co-evolving with their nitrogen-fixing symbionts appear to have independently specialized nsLTPs for this interaction, suggesting a possible convergence of function, which opens perspectives to investigate nsLTPs functions in RNF.

A de novo evolved gene contributes to rice grain shape difference between indica and japonica

The role of de novo evolved genes from non-coding sequences in regulating morphological differentiation between species/subspecies remains largely unknown. Here, we show that a rice de novo gene GSE9 contributes to grain shape difference between indica/xian and japonica/geng varieties. GSE9 evolves from a previous non-coding region of wild rice Oryza rufipogon through the acquisition of start codon. This gene is inherited by most japonica varieties, while the original sequence (absence of start codon, gse9) is present in majority of indica varieties. Knockout of GSE9 in japonica varieties leads to slender grains, whereas introgression to indica background results in round grains. Population evolutionary analyses reveal that gse9 and GSE9 are derived from wild rice Or-I and Or-III groups, respectively. Our findings uncover that the de novo GSE9 gene contributes to the genetic and morphological divergence between indica and japonica subspecies, and provide a target for precise manipulation of rice grain shape.

Seasonal pigment fluctuation in diploid and polyploid Arabidopsis revealed by machine learning-based phenotyping method PlantServation

Long-term field monitoring of leaf pigment content is informative for understanding plant responses to environments distinct from regulated chambers but is impractical by conventional destructive measurements. We developed PlantServation, a method incorporating robust image-acquisition hardware and deep learning-based software that extracts leaf color by detecting plant individuals automatically. As a case study, we applied PlantServation to examine environmental and genotypic effects on the pigment anthocyanin content estimated from leaf color. We processed >4 million images of small individuals of four Arabidopsis species in the field, where the plant shape, color, and background vary over months. Past radiation, coldness, and precipitation significantly affected the anthocyanin content. The synthetic allopolyploid A. kamchatica recapitulated the fluctuations of natural polyploids by integrating diploid responses. The data support a long-standing hypothesis stating that allopolyploids can inherit and combine the traits of progenitors. PlantServation facilitates the study of plant responses to complex environments termed "in natura".

Characterization and organelle genome sequencing of Pyropia species from Myanmar

Pyropia is a genus comprising red algae of the Bangiaceae family that is commonly found in intertidal zones worldwide. However, understanding of Pyropia species that are prone to tropical regions remains limited despite recent breakthroughs in genomic research. Within the realm of Pyropia species thriving in tropical regions, P. vietnamensis stands out as a widely recognized species. In this study, we aimed to investigate Pyropia species in the southwest coast of Myanmar using physiological and molecular approaches, culture-based analyses, chloroplast rbcL and nuclear SSU gene sequencing, and whole chloroplast and mitochondrial genome sequencing. Physiological analysis showed that the Myanmar samples were more heat-tolerant than their Japanese counterparts, including those of subtropical origin. Additionally, molecular characterization revealed that the Myanmar samples were closely related to P. vietnamensis from India. This study is the first to sequence the chloroplast and mitochondrial genomes of Pyropia species from tropical regions. A unique deletion event was observed within a ribosomal RNA gene cluster in the chloroplast genome of the studied Pyropia species, which is a deviation from the usual characteristics of most Pyropia species. This study improves current understanding of the physiological and molecular characteristics of this comparatively understudied Pyropia species that grows in tropical regions.

Comparative and phylogenetic analysis of Asparagus meioclados Levl. and Asparagus munitus Wang et S. C. Chen plastomes and utility of plastomes mutational hotspots

Tiandong is a vital traditional Chinese herbal medicine. It is derived from the tuber root of the Asparagus cochinchinensis according to the Pharmacopoeia of the people's republic of China (2020 Edition). On account of the similar morphology, Asparagus meioclados and Asparagus munitus were used as Tian-Dong in southwest China. Chloroplast (cp) genomes are highly active genetic components of plants and play an extremely important role in improving the efficiency of the identification of plant species. To differentiate the medicinal plants belonging to the genus Asparagus, we sequenced and analyzed the complete plastomes (plastid genomes) of A. meioclados and A. munitus and obtained two plastomes whose length changed to 156,515 bp and 156,381 bp, respectively. A total of 111 unique genes have been detected in plastome, which included 78 protein-coding genes, 29 tRNA genes and 4 rRNA genes. In plastomes of A. meioclados and A. munitus, 14,685 and 14,987 codons were detected, among which 9942 and 10,207 had the relative synonymous codon usage (RSCU) values higher than 1, respectively. A. meioclados and A. munitus have 26 SSRs patterns, among which A. meioclados was 25 and A. munitus 21. The average Ka/Ks value was 0.36, and positive selection was detected in genes of the photosynthetic system (ndhF and rbcL) in Asparagus species. To perform the comparative analysis of plastomes, the two newly sequenced plastomes of the A. meioclados and A. munitus species were compared with that of A. cochinchinensis, and 12 hotspots, including 5 coding regions and 7 inter-genomic regions, were identified. Based on the whole plastome of Asparagus, 2 divergent hotspots (accD and rpl32-trnL-UAG) and 1 international barcode fragment (rbcL) were screened, which may be used as particular molecular markers for the identification of Asparagus species. In addition, we determined the phylogenetic relationship between A. meioclados and A. munitus in the genus Asparagus. This study enriches our knowledge of the molecular evolutionary relationships of the Asparagus genus and provides treasured data records for species identification, molecular breeding, and evolutionary analysis of this genus.

Phylotranscriptomics unveil a Paleoproterozoic-Mesoproterozoic origin and deep relationships of the Viridiplantae

The Viridiplantae comprise two main clades, the Chlorophyta (including a diverse array of marine and freshwater green algae) and the Streptophyta (consisting of the freshwater charophytes and the land plants). Lineages sister to core Chlorophyta, informally refer to as prasinophytes, form a grade of mainly planktonic green algae. Recently, one of these lineages, Prasinodermophyta, which is previously grouped with prasinophytes, has been identified as the sister lineage to both Chlorophyta and Streptophyta. Resolving the deep relationships among green plants is crucial for understanding the historical impact of green algal diversity on marine ecology and geochemistry, but has been proven difficult given the ancient timing of the diversification events. Through extensive taxon and gene sampling, we conduct large-scale phylogenomic analyses to resolve deep relationships and reveal the Prasinodermophyta as the lineage sister to Chlorophyta, raising questions about the necessity of classifying the Prasinodermophyta as a distinct phylum. We unveil that incomplete lineage sorting is the main cause of discordance regarding the placement of Prasinodermophyta. Molecular dating analyses suggest that crown-group green plants and crown-group Prasinodermophyta date back to the Paleoproterozoic-Mesoproterozoic. Our study establishes a plausible link between oxygen levels in the Paleoproterozoic-Mesoproterozoic and the origin of Viridiplantae.

Chromosome level genome assembly of colored calla lily (Zantedeschia elliottiana)

The colored calla lily is an ornamental floral plant native to southern Africa, belonging to the Zantedeschia genus of the Araceae family. We generated a high-quality chromosome-level genome of the colored calla lily, with a size of 1,154 Mb and a contig N50 of 42 Mb. We anchored 98.5% of the contigs (1,137 Mb) into 16 pseudo-chromosomes, and identified 60.18% of the sequences (694 Mb) as repetitive sequences. Functional annotations were assigned to 95.1% of the predicted protein-coding genes (36,165). Additionally, we annotated 469 miRNAs, 1,652 tRNAs, 10,033 rRNAs, and 1,677 snRNAs. Furthermore, Gypsy-type LTR retrotransposons insertions in the genome are the primary factor causing significant genome size variation in Araceae species. This high-quality genome assembly provides valuable resources for understanding genome size differences within the Araceae family and advancing genomic research on colored calla lily.