Sequencing and Analyzing the Transcriptomes of a Thousand Species Across the Tree of Life for Green Plants

The sexual lability hypothesis for the origin of the land plant generation cycle

The evolution of the land plant alternation of generations has been an open question for the past 150 years. Two hypotheses have dominated the discussion: the antithetic hypothesis, which posits that the diploid sporophyte generation arose de novo and gradually increased in complexity, and the homologous hypothesis, which holds that land plant ancestors had independently living sporophytes and haploid gametophytes of similar complexity. Changes in ploidy levels were unknown to early researchers. The antithetic hypothesis is contradicted by generation cycles in Lower Devonian Rhynie chert plants, whose sporophytes and gametophytes have similar morphologies and by some Silurian sporophytes whose complexity exceeds that of Rhynie chert sporophytes. The oldest unambiguous bryophyte gametophytes (thalli) are from the upper Middle Devonian, with an unconnected sporophyte nearby. Based on the 2024 discovery that conjugate algae are paraphyletic to land plants, we present a new hypothesis for the evolution of the land plant generation cycle, focusing on labile ploidy levels and types of reproduction found in conjugate algae. Our 'sexual lability' hypothesis assumes a period of unstable generation cycles (as regards ploidy), likely with predominant clonal growth, as is common in conjugate algae, resulting in sporophytes and gametophytes of similar morphology. When sexual reproduction became stabilized, the timing of gamete fusion, meiosis, and resistant wall formation, which are heterochronic in some conjugate algae, became standardized, with wall formation permanently delayed. In our scenario, independently living adult sporophytes are the land plant ancestral condition, and life-long sporophyte retention on the gametophyte is a bryophyte apomorphy.

The crane fly glycosylated triketide δ-lactone cornicinine elicits akinete differentiation of the cyanobiont in aquatic Azolla fern symbioses

The restriction of plant-symbiont dinitrogen fixation by an insect semiochemical had not been previously described. Here we report on a glycosylated triketide δ-lactone from Nephrotoma cornicina crane flies, cornicinine, that causes chlorosis in the floating-fern symbioses from the genus Azolla. Only the glycosylated trans-A form of chemically synthesized cornicinine was active: 500 nM cornicinine in the growth medium turned all cyanobacterial filaments from Nostoc azollae inside the host leaf-cavities into akinetes typically secreting CTB-bacteriocins. Cornicinine further inhibited akinete germination in Azolla sporelings, precluding re-establishment of the symbiosis during sexual reproduction. It did not impact development of the plant Arabidopsis thaliana or several free-living cyanobacteria from the genera Anabaena or Nostoc but affected the fern host without cyanobiont. Fern-host mRNA sequencing from isolated leaf cavities confirmed high NH4-assimilation and proanthocyanidin biosynthesis in this trichome-rich tissue. After cornicinine treatment, it revealed activation of Cullin-RING ubiquitin-ligase-pathways, known to mediate metabolite signaling and plant elicitation consistent with the chlorosis phenotype, and increased JA-oxidase, sulfate transport and exosome formation. The work begins to uncover molecular mechanisms of cyanobiont differentiation in a seed-free plant symbiosis important for wetland ecology or circular crop-production today, that once caused massive CO2 draw-down during the Eocene geological past.

Biosynthesis and differential spatial distribution of the 3-deoxyanthocyanidins apigenidin and luteolinidin at the interface of a plant-cyanobacteria symbiosis exposed to cold

Aquatic ferns of the genus Azolla (Azolla) form highly productive symbioses with filamentous cyanobacteria fixing N2 in their leaf cavities, Nostoc azollae. Stressed symbioses characteristically turn red due to 3-deoxyanthocyanidin (DA) accumulation, rare in angiosperms and of unknown function. To understand DA accumulation upon cold acclimation and recovery, we integrated laser-desorption-ionization mass-spectrometry-imaging (LDI-MSI), a new Azolla filiculoides genome-assembly and annotation, and dual RNA-sequencing into phenotypic analyses of the symbioses. Azolla sp. Anzali recovered even when cold-induced DA-accumulation was inhibited by abscisic acid. Cyanobacterial filaments generally disappeared upon cold acclimation and Nostoc azollae transcript profiles were unlike those of resting stages formed in cold-resistant sporocarps, yet filaments re-appeared in leaf cavities of newly formed green fronds upon cold-recovery. The high transcript accumulation upon cold acclimation of AfDFR1 encoding a flavanone 4-reductase active in vitro suggested that the enzyme of the first step in the DA-pathway may regulate accumulation of DAs in different tissues. However, LDI-MSI highlighted the necessity to describe metabolite accumulation beyond class assignments as individual DA and caffeoylquinic acid metabolites accumulated differentially. For example, luteolinidin accumulated in epithelial cells, including those lining the leaf cavity, supporting a role for the former in the symbiotic interaction during cold acclimation.

Apomixis and the paradox of sex in plants

BACKGROUND

The predominance of sex in eukaryotes, despite the high costs of meiosis and mating, remains an evolutionary enigma. Many theories have been proposed, none of them being conclusive on its own, and they are, in part, not well applicable to land plants. Sexual reproduction is obligate in embryophytes for the great majority of species.

SCOPE

This review compares the main forms of sexual and asexual reproduction in ferns and angiosperms, based on the generation cycling of sporophyte and gametophyte (leaving vegetative propagation aside). The benefits of sexual reproduction for maintenance of genomic integrity in comparison to asexuality are discussed in the light of developmental, evolutionary, genetic and phylogenetic studies.

CONCLUSIONS

Asexual reproduction represents modifications of the sexual pathway, with various forms of facultative sexuality. For sexual land plants, meiosis provides direct DNA repair mechanisms for oxidative damage in reproductive tissues. The ploidy alternations of meiosis-syngamy cycles and prolonged multicellular stages in the haploid phase in the gametophytes provide a high efficiency of purifying selection against recessive deleterious mutations. Asexual lineages might buffer effects of such mutations via polyploidy and can purge the mutational load via facultative sexuality. The role of organelle-nuclear genome compatibility for maintenance of genome integrity is not well understood. In plants in general, the costs of mating are low because of predominant hermaphroditism. Phylogenetic patterns in the archaeplastid clade suggest that high frequencies of sexuality in land plants are concomitant with a stepwise increase of intrinsic and extrinsic stress factors. Furthermore, expansion of genome size in land plants would increase the potential mutational load. Sexual reproduction appears to be essential for keeping long-term genomic integrity, and only rare combinations of extrinsic and intrinsic factors allow for shifts to asexuality.

Facilitating taxonomy and phylogenetics: An informative and cost-effective protocol integrating long amplicon PCRs and third-generation sequencing

Phylogenetic inference has become a standard technique in integrative taxonomy and systematics, as well as in biogeography and ecology. DNA barcodes are often used for phylogenetic inference, despite being strongly limited due to their low number of informative sites. Also, because current DNA barcodes are based on a fraction of a single, fast-evolving gene, they are highly unsuitable for resolving deeper phylogenetic relationships due to saturation. In recent years, methods that analyse hundreds and thousands of loci at once have improved the resolution of the Tree of Life, but these methods require resources, experience and molecular laboratories that most taxonomists do not have. This paper introduces a PCR-based protocol that produces long amplicons of both slow- and fast-evolving unlinked mitochondrial and nuclear gene regions, which can be sequenced by the affordable and portable ONT MinION platform with low infrastructure or funding requirements. As a proof of concept, we inferred a phylogeny of a sample of 63 spider species from 20 families using our proposed protocol. The results were overall consistent with the results from approaches based on hundreds and thousands of loci, while requiring just a fraction of the cost and labour of such approaches, making our protocol accessible to taxonomists worldwide.

Plant evolution: Streptophyte multicellularity, ecology, and the acclimatisation of plants to life on land

Land plants are celebrated as one of the three great instances of complex multicellularity, but new phylogenomic and phenotypic analyses are revealing deep evolutionary roots of multicellularity among algal relatives, prompting questions about the causal basis of this major evolutionary transition.

MCPtaggR: R package for accurate genotype calling in reduced representation sequencing data by eliminating error-prone markers based on genome comparison

Reduced representation sequencing (RRS) offers cost-effective, high-throughput genotyping platforms such as genotyping-by-sequencing (GBS). RRS reads are typically mapped onto a reference genome. However, mapping reads harbouring mismatches against the reference can potentially result in mismapping and biased mapping, leading to the detection of error-prone markers that provide incorrect genotype information. We established a genotype-calling pipeline named mappable collinear polymorphic tag genotyping (MCPtagg) to achieve accurate genotyping by eliminating error-prone markers. MCPtagg was designed for the RRS-based genotyping of a population derived from a biparental cross. The MCPtagg pipeline filters out error-prone markers prior to genotype calling based on marker collinearity information obtained by comparing the genome sequences of the parents of a population to be genotyped. A performance evaluation on real GBS data from a rice F2 population confirmed its effectiveness. Furthermore, our performance test using a genome assembly that was obtained by genome sequence polishing on an available genome assembly suggests that our pipeline performs well with converted genomes, rather than necessitating de novo assembly. This demonstrates its flexibility and scalability. The R package, MCPtaggR, was developed to provide functions for the pipeline and is available at https://github.com/tomoyukif/MCPtaggR.

Angiosperm-wide analysis of fruit and ovary evolution aided by a new nuclear phylogeny supports association of the same ovary type with both dry and fleshy fruits

Fruit functions in seed protection and dispersal and belongs to many dry and fleshy types, yet their evolutionary pattern remains unclear in part due to uncertainties in the phylogenetic relationships among several orders and families. Thus we used nuclear genes of 502 angiosperm species representing 231 families to reconstruct a well supported phylogeny, with resolved relationships for orders and families with previously uncertain placements. Using this phylogeny as a framework, molecular dating supports a Triassic origin of the crown angiosperms, followed by the emergence of most orders in the Jurassic and Cretaceous and their rise to ecological dominance during the Cretaceous Terrestrial Revolution. The robust phylogeny allowed an examination of the evolutionary pattern of fruit and ovary types, revealing a trend of parallel carpel fusions during early diversifications in eudicots, monocots, and magnoliids. Moreover, taxa in the same order or family with the same ovary type can develop either dry or fleshy fruits with strong correlations between specific types of dry and fleshy fruits; such associations of ovary, dry and fleshy fruits define several ovary-fruit "modules" each found in multiple families. One of the frequent modules has an ovary containing multiple ovules, capsules and berries, and another with an ovary having one or two ovules, achenes (or other single-seeded dry fruits) and drupes. This new perspective of relationships among fruit types highlights the closeness of specific dry and fleshy fruit types, such as capsule and berry, that develop from the same ovary type and belong to the same module relative to dry and fleshy fruits of other modules (such as achenes and drupes). Further analyses of gene families containing known genes for ovary and fruit development identified phylogenetic nodes with multiple gene duplications, supporting a possible role of whole-genome duplications, in combination with climate changes and animal behaviors, in angiosperm fruit and ovary diversification.

The duplication of genomes and genetic networks and its potential for evolutionary adaptation and survival during environmental turmoil

The importance of whole-genome duplication (WGD) for evolution is controversial. Whereas some view WGD mainly as detrimental and an evolutionary dead end, there is growing evidence that polyploidization can help overcome environmental change, stressful conditions, or periods of extinction. However, despite much research, the mechanistic underpinnings of why and how polyploids might be able to outcompete or outlive nonpolyploids at times of environmental upheaval remain elusive, especially for autopolyploids, in which heterosis effects are limited. On the longer term, WGD might increase both mutational and environmental robustness due to redundancy and increased genetic variation, but on the short-or even immediate-term, selective advantages of WGDs are harder to explain. Here, by duplicating artificially generated Gene Regulatory Networks (GRNs), we show that duplicated GRNs-and thus duplicated genomes-show higher signal output variation than nonduplicated GRNs. This increased variation leads to niche expansion and can provide polyploid populations with substantial advantages to survive environmental turmoil. In contrast, under stable environments, GRNs might be maladaptive to changes, a phenomenon that is exacerbated in duplicated GRNs. We believe that these results provide insights into how genome duplication and (auto)polyploidy might help organisms to adapt quickly to novel conditions and to survive ecological uproar or even cataclysmic events.

Channelrhodopsins: From Phototaxis to Optogenetics

Channelrhodopsins stand out among other retinal proteins because of their capacity to generate passive ionic currents following photoactivation. Owing to that, channelrhodopsins are widely used in neuroscience and cardiology as instruments for optogenetic manipulation of the activity of excitable cells. Photocurrents generated by channelrhodopsins were first discovered in the cells of green algae in the 1970s. In this review we describe this discovery and discuss the current state of research in the field.

A plastome phylogeny of Rumex (Polygonaceae) illuminates the divergent evolutionary histories of docks and sorrels

The genus Rumex L. (Polygonaceae) provides a unique system for investigating the evolutionary development of sex determination and molecular rate evolution. Historically, Rumex has been divided, both taxonomically and colloquially into two groups: 'docks' and 'sorrels'. A well-resolved phylogeny can help evaluate a genetic basis for this division. Here we present a plastome phylogeny for 34 species of Rumex, inferred using maximum likelihood criteria. The historical 'docks' (Rumex subgenus Rumex) were resolved as monophyletic. The historical 'sorrels' (Rumex subgenera Acetosa and Acetosella) were resolved together, though not monophyletic due to the inclusion of R. bucephalophorus (Rumex subgenus Platypodium). Emex is supported as its own subgenus within Rumex, instead of resolved as sister taxa. We found remarkably low nucleotide diversity among the docks, consistent with recent diversification in that group, especially as compared to the sorrels. Fossil calibration of the phylogeny suggested that the common ancestor for Rumex (including Emex) has origins in the lower Miocene (22.13 MYA). The sorrels appear to have subsequently diversified at a relatively constant rate. The origin of the docks, however, was placed in the upper Miocene, but with most speciation occurring in the Plio-Pleistocene.

MALDI Imaging Assisted Discovery of a Di-O-glycosyltransferase from Platycodon grandiflorum Root

A matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI MSI) assisted genome mining strategy was developed for the discovery of glycosyltransferase (GT) from the root of Platycodon grandiflorum. A di-O-glycosyltransferase PgGT1 was discovered and characterized that is capable of catalyzing platycoside E (PE) synthesis through the attachment of two β-1,6-linked glucosyl residues sequentially to the glucosyl residue at the C3 position of platycodin D (PD). Although UDP-glucose is the preferred sugar donor for PgGT1, it could also utilize UDP-xylose and UDP-N-acetylglucosamine as weak donors. Residues S273, E274, and H350 played important roles in stabilizing the glucose donor and positioning the glucose in the optimal orientation for the glycosylation reaction. This study clarified two key steps involved in the biosynthetic pathway of PE and could greatly contribute to improving its industrial biotransformation.

Polyploidy: its consequences and enabling role in plant diversification and evolution

BACKGROUND

Most, if not all, green plant (Virdiplantae) species including angiosperms and ferns are polyploids themselves or have ancient polyploid or whole genome duplication signatures in their genomes. Polyploids are not only restricted to our major crop species such as wheat, maize, potato and the brassicas, but also occur frequently in wild species and natural habitats. Polyploidy has thus been viewed as a major driver in evolution, and its influence on genome and chromosome evolution has been at the centre of many investigations. Mechanistic models of the newly structured genomes are being developed that incorporate aspects of sequence evolution or turnover (low-copy genes and regulatory sequences, as well as repetitive DNAs), modification of gene functions, the re-establishment of control of genes with multiple copies, and often meiotic chromosome pairing, recombination and restoration of fertility.

SCOPE

World-wide interest in how green plants have evolved under different conditions - whether in small, isolated populations, or globally - suggests that gaining further insight into the contribution of polyploidy to plant speciation and adaptation to environmental changes is greatly needed. Forward-looking research and modelling, based on cytogenetics, expression studies, and genomics or genome sequencing analyses, discussed in this Special Issue of the Annals of Botany, consider how new polyploids behave and the pathways available for genome evolution. They address fundamental questions about the advantages and disadvantages of polyploidy, the consequences for evolution and speciation, and applied questions regarding the spread of polyploids in the environment and challenges in breeding and exploitation of wild relatives through introgression or resynthesis of polyploids.

CONCLUSION

Chromosome number, genome size, repetitive DNA sequences, genes and regulatory sequences and their expression evolve following polyploidy - generating diversity and possible novel traits and enabling species diversification. There is the potential for ever more polyploids in natural, managed and disturbed environments under changing climates and new stresses.

Inference of Ancient Polyploidy Using Transcriptome Data

Polyploidizations, or whole-genome duplications (WGDs), in plants have increased biological complexity, facilitated evolutionary innovation, and likely enabled adaptation under harsh conditions. Besides genomic data, transcriptome data have been widely employed to detect WGDs, due to their efficient accessibility to the gene space of a species. Age distributions based on synonymous substitutions (so-called K_S age distributions) for paralogs assembled from transcriptome data have identified numerous WGDs in plants, paving the way for further studies on the importance of WGDs for the evolution of seed and flowering plants. However, it is still unclear how transcriptome-based age distributions compare to those based on genomic data. In this chapter, we implemented three different de novo transcriptome assembly pipelines with two popular assemblers, i.e., Trinity and SOAPdenovo-Trans. We selected six plant species with published genomes and transcriptomes to evaluate how assembled transcripts from different pipelines perform when using K_S distributions to detect previously documented WGDs in the six species. Further, using genes predicted in each genome as references, we evaluated the effects of missing genes, gene family clustering, and de novo assembled transcripts on the transcriptome-based K_S distributions. Our results show that, although the transcriptome-based K_S distributions differ from the genome-based ones with respect to their shapes and scales, they are still reasonably reliable for unveiling WGDs, except in species where most duplicates originated from a recent WGD. We also discuss how to overcome some possible pitfalls when using transcriptome data to identify WGDs.

Molecular characterizations of genes in chloroplast genomes of the genus Arachis L. (Fabaceae) based on the codon usage divergence

Studies on the molecular characteristics of chloroplast genome are generally important for clarifying the evolutionary processes of plant species. The base composition, the effective number of codons, the relative synonymous codon usage, the codon bias index, and their correlation coefficients of a total of 41 genes in 21 chloroplast genomes of the genus Arachis were investigated to further perform the correspondence and clustering analyses, revealing significantly higher variations in genomes of wild species than those of the cultivated taxa. The codon usage patterns of all 41 genes in the genus Arachis were AT-rich, suggesting that the natural selection was the main factor affecting the evolutionary history of these genomes. Five genes (i.e., ndhC, petD, atpF, rpl14, and rps11) and five genes (i.e., atpE, psbD, psaB, ycf2, and rps12) showed higher and lower base usage divergences, respectively. This study provided novel insights into our understanding of the molecular evolution of chloroplast genomes in the genus Arachis.

Diversity, phylogeny, and adaptation of bryophytes: insights from genomic and transcriptomic data

Bryophytes including mosses, liverworts, and hornworts are among the earliest land plants, and occupy a crucial phylogenetic position to aid in the understanding of plant terrestrialization. Despite their small size and simple structure, bryophytes are the second largest group of extant land plants. They live ubiquitously in various habitats and are highly diversified, with adaptive strategies to modern ecosystems on Earth. More and more genomes and transcriptomes have been assembled to address fundamental questions in plant biology. Here, we review recent advances in bryophytes associated with diversity, phylogeny, and ecological adaptation. Phylogenomic studies have provided increasing supports for the monophyly of bryophytes, with hornworts sister to the Setaphyta clade including liverworts and mosses. Further comparative genomic analyses revealed that multiple whole-genome duplications might have contributed to the species richness and morphological diversity in mosses. We highlight that the biological changes through gene gain or neofunctionalization that primarily evolved in bryophytes have facilitated the adaptation to early land environments; among the strategies to adapt to modern ecosystems in bryophytes, desiccation tolerance is the most remarkable. More genomic information for bryophytes would shed light on key mechanisms for the ecological success of these 'dwarfs' in the plant kingdom.

Using interdisciplinary, phylogeny-guided approaches to understand the evolution of plant metabolism

To cope with relentless environmental pressures, plants produce an arsenal of structurally diverse chemicals, often called specialized metabolites. These lineage-specific compounds are derived from the simple building blocks made by ubiquitous core metabolic pathways. Although the structures of many specialized metabolites are known, the underlying metabolic pathways and the evolutionary events that have shaped the plant chemical diversity landscape are only beginning to be understood. However, with the advent of multi-omics data sets and the relative ease of studying pathways in previously intractable non-model species, plant specialized metabolic pathways are now being systematically identified. These large datasets also provide a foundation for comparative, phylogeny-guided studies of plant metabolism. Comparisons of metabolic traits and features like chemical abundances, enzyme activities, or gene sequences from phylogenetically diverse plants provide insights into how metabolic pathways evolved. This review highlights the power of studying evolution through the lens of comparative biochemistry, particularly how placing metabolism into a phylogenetic context can help a researcher identify the metabolic innovations enabling the evolution of structurally diverse plant metabolites.

Novel Approaches for Species Concepts and Delimitation in Polyploids and Hybrids

Hybridization and polyploidization are important processes for plant evolution. However, classification of hybrid or polyploid species has been notoriously difficult because of the complexity of processes and different evolutionary scenarios that do not fit with classical species concepts. Polyploid complexes are formed via combinations of allopolyploidy, autopolyploidy and homoploid hybridization with persisting sexual reproduction, resulting in many discrete lineages that have been classified as species. Polyploid complexes with facultative apomixis result in complicated net-work like clusters, or rarely in agamospecies. Various case studies illustrate the problems that apply to traditional species concepts to hybrids and polyploids. Conceptual progress can be made if lineage formation is accepted as an inevitable consequence of meiotic sex, which is established already in the first eukaryotes as a DNA restoration tool. The turnaround of the viewpoint that sex forms species as lineages helps to overcome traditional thinking of species as "units". Lineage formation and self-sustainability is the prerequisite for speciation and can also be applied to hybrids and polyploids. Species delimitation is aided by the improved recognition of lineages via various novel -omics methods, by understanding meiosis functions, and by recognizing functional phenotypes by considering morphological-physiological-ecological adaptations.

Spatial and Temporal Transcriptomic Heredity and Asymmetry in an Artificially Constructed Allotetraploid Wheat (AADD)

Polyploidy, or whole-genome duplication (WGD), often induces dramatic changes in gene expression due to "transcriptome shock. " However, questions remain about how allopolyploidy (the merging of multiple nuclear genomes in the same nucleus) affects gene expression within and across multiple tissues and developmental stages during the initial foundation of allopolyploid plants. Here, we systematically investigated the immediate effect of allopolyploidy on gene expression variation in an artificial allopolyploidy system consisting of a constructed allotetraploid wheat (AADD genome, accession AT2) and its diploid progenitors Triticum urartu and Aegilops tauschii. We performed comprehensive RNA sequencing of 81 samples from different genotypes, tissues, and developmental stages. First, we found that intrinsic interspecific differences between the diploid parents played a major role in establishing the expression architecture of the allopolyploid. Nonetheless, allopolyploidy per se also induced dramatic and asymmetric patterns of differential gene expression between the subgenomes, and genes from the D subgenome exhibited a more drastic response. Second, analysis of homoeolog expression bias (HEB) revealed that the D subgenome exhibited significant expression bias and that de novo-generated HEB was attributed mainly to asymmetrical differential gene expression. Homoeolog-specific expression (HSE) analyses showed that the cis-only regulatory pattern was predominant in AT2, reflecting significant divergence between the parents. Co-expression network analysis revealed that homoeolog expression connectivity (HEC) was significantly correlated with sequence divergence in cis elements between subgenomes. Interestingly, allopolyploidy-induced reconstruction of network modules was also associated with different HSE patterns. Finally, a transcriptome atlas of spike development demonstrated that the phenotypic similarity of AT2 to T. urartu may be attributed to the combination of relatively stable expression of A-subgenome genes and drastic downregulation of their D-subgenome homoeologs. These findings provide a broad, multidimensional characterization of allopolyploidy-induced transcriptomic responses and suggest that allopolyploidy can have immediate and complex regulatory effects on the expression of nuclear genes.

Insights into the molecular evolution of fertilization mechanism in land plants

Comparative genetics and genomics among green plants, including algae, provide deep insights into the evolution of land plant sexual reproduction. Land plants have evolved successive changes during their conquest of the land and innovations in sexual reproduction have played a major role in their terrestrialization. Recent years have seen many revealing dissections of the molecular mechanisms of sexual reproduction and much new genomics data from the land plant lineage, including early diverging land plants, as well as algae. This new knowledge is being integrated to further understand how sexual reproduction in land plants evolved, identifying highly conserved factors and pathways, but also molecular changes that underpinned the emergence of new modes of sexual reproduction. Here, we review recent advances in the knowledge of land plant sexual reproduction from an evolutionary perspective and also revisit the evolution of angiosperm double fertilization.

Taxus yunnanensis genome offers insights into gymnosperm phylogeny and taxol production

Taxol, a natural product derived from Taxus, is one of the most effective natural anticancer drugs and the biosynthetic pathway of Taxol is the basis of heterologous bio-production. Here, we report a high-quality genome assembly and annotation of Taxus yunnanensis based on 10.7 Gb sequences assembled into 12 chromosomes with contig N50 and scaffold N50 of 2.89 Mb and 966.80 Mb, respectively. Phylogenomic analyses show that T. yunnanensis is most closely related to Sequoiadendron giganteum among the sampled taxa, with an estimated divergence time of 133.4-213.0 MYA. As with most gymnosperms, and unlike most angiosperms, there is no evidence of a recent whole-genome duplication in T. yunnanensis. Repetitive sequences, especially long terminal repeat retrotransposons, are prevalent in the T. yunnanensis genome, contributing to its large genome size. We further integrated genomic and transcriptomic data to unveil clusters of genes involved in Taxol synthesis, located on the chromosome 12, while gene families encoding hydroxylase in the Taxol pathway exhibited significant expansion. Our study contributes to the further elucidation of gymnosperm relationships and the Taxol biosynthetic pathway.

The evolutionary emergence of land plants

There can be no doubt that early land plant evolution transformed the planet but, until recently, how and when this was achieved was unclear. Coincidence in the first appearance of land plant fossils and formative shifts in atmospheric oxygen and CO2 are an artefact of the paucity of earlier terrestrial rocks. Disentangling the timing of land plant bodyplan assembly and its impact on global biogeochemical cycles has been precluded by uncertainty concerning the relationships of bryophytes to one another and to the tracheophytes, as well as the timescale over which these events unfolded. New genome and transcriptome sequencing projects, combined with the application of sophisticated phylogenomic modelling methods, have yielded increasing support for the Setaphyta clade of liverworts and mosses, within monophyletic bryophytes. We consider the evolution of anatomy, genes, genomes and of development within this phylogenetic context, concluding that many vascular plant (tracheophytes) novelties were already present in a comparatively complex last common ancestor of living land plants (embryophytes). Molecular clock analyses indicate that embryophytes emerged in a mid-Cambrian to early Ordovician interval, compatible with hypotheses on their role as geoengineers, precipitating early Palaeozoic glaciations.

Early "Rootprints" of Plant Terrestrialization: Selaginella Root Development Sheds Light on Root Evolution in Vascular Plants

Roots provide multiple key functions for plants, including anchorage and capturing of water and nutrients. Evolutionarily, roots represent a crucial innovation that enabled plants to migrate from aquatic to terrestrial environment and to grow in height. Based on fossil evidence, roots evolved at least twice independently, once in the lycophyte clade and once in the euphyllophyte (ferns and seed plants) clade. In lycophytes, roots originated in a stepwise manner. Despite their pivotal position in root evolution, it remains unclear how root development is controlled in lycophytes. Getting more insight into lycophyte root development might shed light on how genetic players controlling the root meristem and root developmental processes have evolved. Unfortunately, genetic studies in lycophytes are lagging behind, lacking advanced biotechnological tools, partially caused by the limited economic value of this clade. The technology of RNA sequencing (RNA-seq) at least enabled transcriptome studies, which could enhance the understanding or discovery of genes involved in the root development of this sister group of euphyllophytes. Here, we provide an overview of the current knowledge on root evolution followed by a survey of root developmental events and how these are genetically and hormonally controlled, starting from insights obtained in the model seed plant Arabidopsis and where possible making a comparison with lycophyte root development. Second, we suggest possible key genetic regulators in root development of lycophytes mainly based on their expression profiles in Selaginella moellendorffii and phylogenetics. Finally, we point out challenges and possible future directions for research on root evolution.

Asymmetric expansions of FT and TFL1 lineages characterize differential evolution of the EuPEBP family in the major angiosperm lineages

Background

In flowering plants, precise timing of the floral transition is crucial to maximize chances of reproductive success, and as such, this process has been intensively studied. FLOWERING LOCUS T (FT) and TERMINAL FLOWER1 (TFL1) have been identified as closely related eukaryotic phosphatidylethanolamine-binding proteins (‘EuPEBPs’) that integrate multiple environmental stimuli, and act antagonistically to determine the optimal timing of the floral transition. Extensive research has demonstrated that FT acts similar to hormonal signals, being transported in the phloem from its primary site of expression in leaves to its primary site of action in the shoot meristem; TFL1 also appears to act as a mobile signal. Recent work implicates FT, TFL1, and the other members of the EuPEBP family, in the control of other important processes, suggesting that the EuPEBP family may be key general regulators of developmental transitions in flowering plants. In eudicots, there are a small number of EuPEBP proteins, but in monocots, and particularly grasses, there has been a large, but uncharacterized expansion of EuPEBP copy number, with unknown consequences for the EuPEBP function.

Results

To systematically characterize the evolution of EuPEBP proteins in flowering plants, and in land plants more generally, we performed a high-resolution phylogenetic analysis of 701 PEBP sequences from 208 species. We refine previous models of EuPEBP evolution in early land plants, demonstrating the algal origin of the family, and pin-pointing the origin of the FT/TFL1 clade at the base of monilophytes. We demonstrate how a core set of genes (MFT1, MFT2, FT, and TCB) at the base of flowering plants has undergone differential evolution in the major angiosperm lineages. This includes the radical expansion of the FT family in monocots into 5 core lineages, further re-duplicated in the grass family to 12 conserved clades.

Conclusions

We show that many grass FT proteins are strongly divergent from other FTs and are likely neo-functional regulators of development. Our analysis shows that monocots and eudicots have strongly divergent patterns of EuPEBP evolution.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01128-8.

The thionin family of antimicrobial peptides

Thionins are antimicrobial peptides found only in plants. They are first produced as preproproteins and then processed to yield the usually 5 kDa, basic thionin peptide with three or four disulfide bridges. So far, thionins had only been found in some plant families of angiosperms. The One Thousand Plant Transcriptomes Initiative (1KP project) has sequenced the transcriptomes of more than 1000 plant species. We have used these data to search for new thionin sequences which gave 225 hits. After removing doublets these resulted in 133 new thionins. No sequences were found in algae and mosses. The phylogenetically earliest hits were from Selaginella species and from conifers. Many hits were from angiosperm plant families which were previously not known to contain thionins. A large gene family for thionins was found in Papaver. We isolated a genomic clone from Papaver somniferum which confirmed the general genomic structure with two small introns within the acidic domain. We also expressed the thionin encoded by the genomic clone and found that it had antimicrobial activity in vitro, especially against fungi. Previously, we had grouped thionins into four classes. The new data reported here led us to revise this classification. We now recognize only class 1 thionins with eight cysteine residues and class 2 thionins with six cysteine residues. The different variants that we found (and also previously known variants) can all be traced back to one of these two classes. Some of the variants had an uneven number of cysteine residues and it is not clear at the moment what that means for their threedimensional structure.

Disparate Dynamics of Gene Body and cis-Regulatory Element Evolution Illustrated for the Senescence-Associated Cysteine Protease Gene SAG12 of Plants

Gene regulation networks precisely orchestrate the expression of genes that are closely associated with defined physiological and developmental processes such as leaf senescence in plants. The Arabidopsis thaliana senescence-associated gene 12 (AtSAG12) encodes a cysteine protease that is (i) involved in the degradation of chloroplast proteins and (ii) almost exclusively expressed during senescence. Transcription factors, such as WRKY53 and WRKY45, bind to W-boxes in the promoter region of AtSAG12 and play key roles in its activation. Other transcription factors, such as bZIPs, might have accessory functions in their gene regulation, as several A-boxes have been identified and appear to be highly overrepresented in the promoter region compared to the whole genome distribution but are not localized within the regulatory regions driving senescence-associated expression. To address whether these two regulatory elements exhibiting these different properties are conserved in other closely related species, we constructed phylogenetic trees of the coding sequences of orthologs of AtSAG12 and screened their respective 2000 bp promoter regions for the presence of conserved cis-regulatory elements, such as bZIP and WRKY binding sites. Interestingly, the functional relevant upstream located W-boxes were absent in plant species as closely related as Arabidopsis lyrata, whereas an A-box cluster appeared to be conserved in the Arabidopsis species but disappeared in Brassica napus. Several orthologs were present in other species, possibly because of local or whole genome duplication events, but with distinct cis-regulatory sites in different locations. However, at least one gene copy in each family analyzed carried one W-box and one A-box in its promoter. These gene differences in SAG12 orthologs are discussed in the framework of cis- and trans-regulatory factors, of promoter and gene evolution, of genetic variation, and of the enhancement of the adaptability of plants to changing environmental conditions.

Hundreds of nuclear and plastid loci yield novel insights into orchid relationships

PREMISE

The inference of evolutionary relationships in the species-rich family Orchidaceae has hitherto relied heavily on plastid DNA sequences and limited taxon sampling. Previous studies have provided a robust plastid phylogenetic framework, which was used to classify orchids and investigate the drivers of orchid diversification. However, the extent to which phylogenetic inference based on the plastid genome is congruent with the nuclear genome has been only poorly assessed.

METHODS

We inferred higher-level phylogenetic relationships of orchids based on likelihood and ASTRAL analyses of 294 low-copy nuclear genes sequenced using the Angiosperms353 universal probe set for 75 species (representing 69 genera, 16 tribes, 24 subtribes) and a concatenated analysis of 78 plastid genes for 264 species (117 genera, 18 tribes, 28 subtribes). We compared phylogenetic informativeness and support for the nuclear and plastid phylogenetic hypotheses.

RESULTS

Phylogenetic inference using nuclear data sets provides well-supported orchid relationships that are highly congruent between analyses. Comparisons of nuclear gene trees and a plastid supermatrix tree showed that the trees are mostly congruent, but revealed instances of strongly supported phylogenetic incongruence in both shallow and deep time. The phylogenetic informativeness of individual Angiosperms353 genes is in general better than that of most plastid genes.

CONCLUSIONS

Our study provides the first robust nuclear phylogenomic framework for Orchidaceae and an assessment of intragenomic nuclear discordance, plastid-nuclear tree incongruence, and phylogenetic informativeness across the family. Our results also demonstrate what has long been known but rarely thoroughly documented: nuclear and plastid phylogenetic trees can contain strongly supported discordances, and this incongruence must be reconciled prior to interpretation in evolutionary studies, such as taxonomy, biogeography, and character evolution.

Recent Advances and Future Perspectives in Cotton Research

Cotton is not only the world's most important natural fiber crop, but it is also an ideal system in which to study genome evolution, polyploidization, and cell elongation. With the assembly of five different cotton genomes, a cotton-specific whole-genome duplication with an allopolyploidization process that combined the A- and D-genomes became evident. All existing A-genomes seemed to originate from the A₀-genome as a common ancestor, and several transposable element bursts contributed to A-genome size expansion and speciation. The ethylene production pathway is shown to regulate fiber elongation. A tip-biased diffuse growth mode and several regulatory mechanisms, including plant hormones, transcription factors, and epigenetic modifications, are involved in fiber development. Finally, we describe the involvement of the gossypol biosynthetic pathway in the manipulation of herbivorous insects, the role of GoPGF in gland formation, and host-induced gene silencing for pest and disease control. These new genes, modules, and pathways will accelerate the genetic improvement of cotton.

The Enigma of Interspecific Plasmodesmata: Insight From Parasitic Plants

Parasitic plants live in intimate physical connection with other plants serving as their hosts. These host plants provide the inorganic and organic compounds that the parasites need for their propagation. The uptake of the macromolecular compounds happens through symplasmic connections in the form of plasmodesmata. In contrast to regular plasmodesmata, which connect genetically identical cells of an individual plant, the plasmodesmata that connect the cells of host and parasite join separate individuals belonging to different species and are therefore termed "interspecific". The existence of such interspecific plasmodesmata was deduced either indirectly using molecular approaches or observed directly by ultrastructural analyses. Most of this evidence concerns shoot parasitic Cuscuta species and root parasitic Orobanchaceae, which can both infect a large range of phylogenetically distant hosts. The existence of an interspecific chimeric symplast is both striking and unique and, with exceptions being observed in closely related grafted plants, exist only in these parasitic relationships. Considering the recent technical advances and upcoming tools for analyzing parasitic plants, interspecific plasmodesmata in parasite/host connections are a promising system for studying secondary plasmodesmata. For open questions like how their formation is induced, how their positioning is controlled and if they are initiated by one or both bordering cells simultaneously, the parasite/host interface with two adjacent distinguishable genetic systems provides valuable advantages. We summarize here what is known about interspecific plasmodesmata between parasitic plants and their hosts and discuss the potential of the intriguing parasite/host system for deepening our insight into plasmodesmatal structure, function, and development.

Polyploidy: an evolutionary and ecological force in stressful times

Polyploidy has been hypothesized to be both an evolutionary dead-end and a source for evolutionary innovation and species diversification. Although polyploid organisms, especially plants, abound, the apparent nonrandom long-term establishment of genome duplications suggests a link with environmental conditions. Whole-genome duplications seem to correlate with periods of extinction or global change, while polyploids often thrive in harsh or disturbed environments. Evidence is also accumulating that biotic interactions, for instance, with pathogens or mutualists, affect polyploids differently than nonpolyploids. Here, we review recent findings and insights on the effect of both abiotic and biotic stress on polyploids versus nonpolyploids and propose that stress response in general is an important and even determining factor in the establishment and success of polyploidy.

Ancestral sequence reconstruction - An underused approach to understand the evolution of gene function in plants?

Whilst substantial research effort has been placed on understanding the interactions of plant proteins with their molecular partners, relatively few studies in plants - by contrast to work in other organisms - address how these interactions evolve. It is thought that ancestral proteins were more promiscuous than modern proteins and that specificity often evolved following gene duplication and subsequent functional refining. However, ancestral protein resurrection studies have found that some modern proteins have evolved de novo from ancestors lacking those functions. Intriguingly, the new interactions evolved as a consequence of just a few mutations and, as such, acquisition of new functions appears to be neither difficult nor rare, however, only a few of them are incorporated into biological processes before they are lost to subsequent mutations. Here, we detail the approach of ancestral sequence reconstruction (ASR), providing a primer to reconstruct the sequence of an ancestral gene. We will present case studies from a range of different eukaryotes before discussing the few instances where ancestral reconstructions have been used in plants. As ASR is used to dig into the remote evolutionary past, we will also present some alternative genetic approaches to investigate molecular evolution on shorter timescales. We argue that the study of plant secondary metabolism is particularly well suited for ancestral reconstruction studies. Indeed, its ancient evolutionary roots and highly diverse landscape provide an ideal context in which to address the focal issue around the emergence of evolutionary novelties and how this affects the chemical diversification of plant metabolism.

Extremely low nucleotide diversity among thirty-six new chloroplast genome sequences from Aldama (Heliantheae, Asteraceae) and comparative chloroplast genomics analyses with closely related genera

Aldama (Heliantheae, Asteraceae) is a diverse genus in the sunflower family. To date, nearly 200 Asteraceae chloroplast genomes have been sequenced, but the plastomes of Aldama remain undescribed. Plastomes in Asteraceae usually show little sequence divergence, consequently, our hypothesis is that species of Aldama will be overall conserved. In this study, we newly sequenced 36 plastomes of Aldama and of five species belonging to other Heliantheae genera selected as outgroups (i.e., Dimerostemma asperatum, Helianthus tuberosus, Iostephane heterophylla, Pappobolus lanatus var. lanatus, and Tithonia diversifolia). We analyzed the structure and gene content of the assembled plastomes and performed comparative analyses within Aldama and with other closely related genera. As expected, Aldama plastomes are very conserved, with the overall gene content and orientation being similar in all studied species. The length of the plastome is also consistent and the junction between regions usually contain the same genes and have similar lengths. A large ∼20 kb and a small ∼3 kb inversion were detected in the Large Single Copy (LSC) regions of all assembled plastomes, similarly to other Asteraceae species. The nucleotide diversity is very low, with only 1,509 variable sites in 127,466 bp (i.e., 1.18% of the sites in the alignment of 36 Aldama plastomes, with one of the IRs removed, is variable). Only one gene, rbcL, shows signatures of positive selection. The plastomes of the selected outgroups feature a similar gene content and structure compared to Aldama and also present the two inversions in the LSC region. Deletions of different lengths were observed in the gene ycf2. Multiple SSRs were identified for the sequenced Aldama and outgroups. The phylogenetic analysis shows that Aldama is not monophyletic due to the position of the Mexican species A. dentata. All Brazilian species form a strongly supported clade. Our results bring new understandings into the evolution and diversity of plastomes at the species level.

Homoeolog expression bias and expression level dominance (ELD) in four tissues of natural allotetraploid Brassica napus

Background

Allopolyploidy is widespread in angiosperms, and they can coordinate two or more different genomes through genetic and epigenetic modifications to exhibit stronger vigor and adaptability. To explore the changes in homologous gene expression patterns in the natural allotetraploid Brassica napus (A_nA_nC_nC_n) relative to its two diploid progenitors, B. rapa (A_rA_r) and B. oleracea (C_oC_o), after approximately 7500 years of domestication, the global gene pair expression patterns in four major tissues (stems, leaves, flowers and siliques) of these three species were analyzed using an RNA sequencing approach.

Results

The results showed that the ‘transcriptomic shock’ phenomenon was alleviated in natural B. napus after approximately 7500 years of natural domestication, and most differentially expressed genes (DEGs) in B. napus were downregulated relative to those in its two diploid progenitors. The KEGG analysis indicated that three pathways related to photosynthesis were enriched in both comparison groups (A_nA_nC_nC_n vs A_rA_r and A_nA_nC_nC_n vs C_oC_o), and these pathways were all downregulated in four tissues of B. napus. In addition, homoeolog expression bias and expression level dominance (ELD) in B. napus were thoroughly studied through analysis of expression levels of 27,609 B. rapa-B. oleracea orthologous gene pairs. The overwhelming majority of gene pairs (an average of 86.7%) in B. napus maintained their expression pattern in two diploid progenitors, and approximately 78.1% of the gene pairs showed expression bias with a preference toward the A subgenome. Overall, an average of 48, 29.7 and 22.3% homologous gene pairs exhibited additive expression, ELD and transgressive expression in B. napus, respectively. The ELD bias varies from tissue to tissue; specifically, more gene pairs in stems and siliques showed ELD-A, whereas the opposite was observed in leaves and flowers. More transgressive upregulation, rather than downregulation, was observed in gene pairs of B. napus.

Conclusions

In general, these results may provide a comprehensive understanding of the changes in homologous gene expression patterns in natural B. napus after approximately 7500 years of evolution and domestication and may enhance our understanding of allopolyploidy.

The Peptidoglycan Biosynthesis Gene murC in Frankia: Actinorhizal vs. Plant Type

Nitrogen-fixing Actinobacteria of the genus Frankia can be subdivided into four phylogenetically distinct clades; members of clusters one to three engage in nitrogen-fixing root nodule symbioses with actinorhizal plants. Mur enzymes are responsible for the biosynthesis of the peptidoglycan layer of bacteria. The four Mur ligases,MurC, MurD, MurE, and MurF, catalyse the addition of a short polypeptide to UDP-N-acetylmuramic acid. Frankia strains of cluster-2 and cluster-3 contain two copies of murC, while the strains of cluster-1 and cluster-4 contain only one. Phylogenetically, the protein encoded by the murC gene shared only by cluster-2 and cluster-3, termed MurC1, groups with MurC proteins of other Actinobacteria. The protein encoded by the murC gene found in all Frankia strains, MurC2, shows a higher similarity to the MurC proteins of plants than of Actinobacteria. MurC2 could have been either acquired via horizontal gene transfer or via gene duplication and convergent evolution, while murC1 was subsequently lost in the cluster-1 and cluster-4 strains. In the nodules induced by the cluster-2 strains, the expression levels of murC2 were significantly higher than those of murC1. Thus, there is clear sequence divergence between both types of Frankia MurC, and Frankia murC1 is in the process of being replaced by murC2, indicating selection in favour of murC2. Nevertheless, protein modelling showed no major structural differences between the MurCs from any phylogenetic group examined.