Gene duplication as a driver of plant morphogenetic evolution

Characteristics of duplicated gene expression and DNA methylation regulation in different tissues of allopolyploid Brassica napus

Plant polyploidization increases the complexity of epigenomes and transcriptional regulation, resulting in genome evolution and enhanced adaptability. However, few studies have been conducted on the relationship between gene expression and epigenetic modification in different plant tissues after allopolyploidization. In this study, we studied gene expression and DNA methylation modification patterns in four tissues (stems, leaves, flowers and siliques) of Brassica napusand its diploid progenitors. On this basis, the alternative splicing patterns and cis-trans regulation patterns of four tissues in B. napus and its diploid progenitors were also analyzed. It can be seen that the number of alternative splicing occurs in the B. napus is higher than that in the diploid progenitors, and the IR type increases the most during allopolyploidy. In addition, we studied the fate changes of duplicated genes after allopolyploidization in B. napus. We found that the fate of most duplicated genes is conserved, but the number of neofunctionalization and specialization is also large. The genetic fate of B. napus was classified according to five replication types (WGD, PD, DSD, TD, TRD). This study also analyzed generational transmission analysis of expression and DNA methylation patterns. Our study provides a reference for the fate differentiation of duplicated genes during allopolyploidization.

Identification and Characterization of the Glutathione S-Transferase Gene Family in Blueberry (Vaccinium corymbosum) and Their Potential Roles in Anthocyanin Intracellular Transportation

The glutathione S-transferases (GSTs, EC 2.5.1.18) constitute a versatile enzyme family with pivotal roles in plant stress responses and detoxification processes. Recent discoveries attributed the additional function of facilitating anthocyanin intracellular transportation in plants to GSTs. Our study identified 178 VcGST genes from 12 distinct subfamilies in the blueberry genome. An uneven distribution was observed among these genes across blueberry's chromosomes. Members within the same subfamily displayed homogeneity in gene structure and conserved protein motifs, whereas marked divergence was noted among subfamilies. Functional annotations revealed that VcGSTs were significantly enriched in several gene ontology and KEGG pathway categories. Promoter regions of VcGST genes predominantly contain light-responsive, MYB-binding, and stress-responsive elements. The majority of VcGST genes are subject to purifying selection, with whole-genome duplication or segmental duplication serving as key processes that drive the expansion of the VcGST gene family. Notably, during the ripening of the blueberry fruit, 100 VcGST genes were highly expressed, and the expression patterns of 24 of these genes demonstrated a strong correlation with the dynamic content of fruit anthocyanins. Further analysis identified VcGSTF8, VcGSTF20, and VcGSTF22 as prime candidates of VcGST genes involved in the anthocyanin intracellular transport. This study provides a reference for the exploration of anthocyanin intracellular transport mechanisms and paves the way for investigating the spectrum of GST functions in blueberries.

Sequence similarity networks bear out hierarchical relationships of green cytochrome P450

Land plants have diversified enzyme families. One of the most prominent is the cytochrome P450 (CYP or CYP450) family. With over 443,000 CYP proteins sequenced across the tree of life, CYPs are ubiquitous in archaea, bacteria, and eukaryotes. Here, we focused on land plants and algae to study the role of CYP diversification. CYPs, acting as monooxygenases, catalyze hydroxylation reactions crucial for specialized plant metabolic pathways, including detoxification and phytohormone production; the CYPome consists of one enormous superfamily that is divided into clans and families. Their evolutionary history speaks of high substrate promiscuity; radiation and functional diversification have yielded numerous CYP families. To understand the evolutionary relationships within the CYPs, we employed sequence similarity network analyses. We recovered distinct clusters representing different CYP families, reflecting their diversified sequences that we link to the prediction of functionalities. Hierarchical clustering and phylogenetic analysis further elucidated relationships between CYP clans, uncovering their shared deep evolutionary history. We explored the distribution and diversification of CYP subfamilies across plant and algal lineages, uncovering novel candidates and providing insights into the evolution of these enzyme families. This identified unexpected relationships between CYP families, such as the link between CYP82 and CYP74, shedding light on their roles in plant defense signaling pathways. Our approach provides a methodology that brings insights into the emergence of new functions within the CYP450 family, contributing to the evolutionary history of plants and algae. These insights can be further validated and implemented via experimental setups under various external conditions.

The Genome of Vitis zhejiang-adstricta Strengthens the Protection and Utilization of the Endangered Ancient Grape Endemic to China

Vitis zhejiang-adstricta (V. zhejiang-adstricta) is one of the most important and endangered wild grapes. It is a national key protected wild, rare and endangered ancient grape endemic to China and used as a candidate material for resistance breeding owing to its excellent significant disease resistance. Here, we present a high-quality chromosome-level assembly of V. zhejiang-adstricta (IB-VB-01), comprising 506.66 Mb assembled into 19 pseudo-chromosomes. The contig N50 length is 3.91 Mb with 31,196 annotated protein-coding genes. Comparative genome and evolutionary analyses illustrated that V. zhejiang-adstricta has a specific position in the evolution of East Asian Vitis and shared a common ancestor with Vitis vinifera during the divergence of the two species about 10.42 (between 9.34 and 11.12) Mya. The expanded gene families compared with those in plants were related to disease resistance, and constructed gene families were related to plant growth and primary metabolism. With the analysis of gene family expansion and contraction, the evolution of environmental adaptability and especially the NBS-LRR gene family of V. zhejiang-adstricta was elucidated based on the pathways of resistance genes (R genes), unique genes and structural variations. The near-complete and accurate diploid V. zhejiang-adstricta reference genome obtained herein serves as an important complement to wild grape genomes and will provide valuable genomic resources for investigating the genomic architecture of V. zhejiang-adstricta as well as for improving disease resistance breeding strategies in grape.

MpTGA, together with MpNPR, regulates sexual reproduction and independently affects oil body formation in Marchantia polymorpha

In angiosperms, basic leucine-zipper (bZIP) TGACG-motif-binding (TGA) transcription factors (TFs) regulate developmental and stress-related processes, the latter often involving NON EXPRESSOR OF PATHOGENESIS-RELATED GENES (NPR) coregulator interactions. To gain insight into their functions in an early diverging land-plant lineage, the single MpTGA and sole MpNPR genes were investigated in the liverwort Marchantia polymorpha. We generated Marchantia MpTGA and MpNPR knockout and overexpression mutants and conducted morphological, transcriptomic and expression studies. Furthermore, we investigated MpTGA interactions with wild-type and mutagenized MpNPR and expanded our analyses including TGA TFs from two streptophyte algae. Mptga mutants fail to induce the switch from vegetative to reproductive development and lack gametangiophore formation. MpTGA and MpNPR proteins interact and Mpnpr mutant analysis reveals a novel coregulatory NPR role in sexual reproduction. Additionally, MpTGA acts independently of MpNPR as a repressor of oil body (OB) formation and can thereby affect herbivory. The single MpTGA TF exerts a dual role in sexual reproduction and OB formation in Marchantia. Common activities of MpTGA/MpNPR in sexual development suggest that coregulatory interactions were established after emergence of land-plant-specific NPR genes and contributed to the diversification of TGA TF functions during land-plant evolution.

Investigation of regulatory divergence between homoeologs in the recently formed allopolyploids, Tragopogon mirus and T. miscellus (Asteraceae)

Polyploidy is an important evolutionary process throughout eukaryotes, particularly in flowering plants. Duplicated gene pairs (homoeologs) in allopolyploids provide additional genetic resources for changes in molecular, biochemical, and physiological mechanisms that result in evolutionary novelty. Therefore, understanding how divergent genomes and their regulatory networks reconcile is vital for unraveling the role of polyploidy in plant evolution. Here, we compared the leaf transcriptomes of recently formed natural allotetraploids (Tragopogon mirus and T. miscellus) and their diploid parents (T. porrifolius X T. dubius and T. pratensis X T. dubius, respectively). Analysis of 35 400 expressed loci showed a significantly higher level of transcriptomic additivity compared to old polyploids; only 22% were non-additively expressed in the polyploids, with 5.9% exhibiting transgressive expression (lower or higher expression in the polyploids than in the diploid parents). Among approximately 7400 common orthologous regions (COREs), most loci in both allopolyploids exhibited expression patterns that were vertically inherited from their diploid parents. However, 18% and 20.3% of the loci showed novel expression bias patterns in T. mirus and T. miscellus, respectively. The expression changes of 1500 COREs were explained by cis-regulatory divergence (the condition in which the two parental subgenomes do not interact) between the diploid parents, whereas only about 423 and 461 of the gene expression changes represent trans-effects (the two parental subgenomes interact) in T. mirus and T. miscellus, respectively. The low degree of both non-additivity and trans-effects on gene expression may present the ongoing evolutionary processes of the newly formed Tragopogon polyploids (~80-90 years).

Identification, evolution, and expression of GDSL-type Esterase/Lipase (GELP) gene family in three cotton species: a bioinformatic analysis

BACKGROUND

GDSL esterase/lipases (GELPs) play important roles in plant growth, development, and response to biotic and abiotic stresses. Presently, an extensive and in-depth analysis of GELP family genes in cotton is still not clear enough, which greatly limits the further understanding of cotton GELP function and regulatory mechanism.

RESULTS

A total of 389 GELP family genes were identified in three cotton species of Gossypium hirsutum (193), G. arboreum (97), and G. raimondii (99). These GELPs could be classified into three groups and eight subgroups, with the GELPs in same group to have similar gene structures and conserved motifs. Evolutionary event analysis showed that the GELP family genes tend to be diversified at the spatial dimension and certain conservative at the time dimension, with a trend of potential continuous expansion in the future. The orthologous or paralogous GELPs among different genomes/subgenomes indicated the inheritance from genome-wide duplication during polyploidization, and the paralogous GELPs were derived from chromosomal segment duplication or tandem replication. GELP genes in the A/D subgenome underwent at least three large-scale replication events in the evolutionary process during the period of 0.6-3.2 MYA, with two large-scale evolutionary events between 0.6-1.8 MYA that were associated with tetraploidization, and the large-scale duplication between 2.6-9.1 MYA that occurred during diploidization. The cotton GELPs indicated diverse expression patterns in tissue development, ovule and fiber growth, and in response to biotic and abiotic stresses, combining the existing cis-elements in the promoter regions, suggesting the GELPs involvements of functions to be diversification and of the mechanisms to be a hormone-mediated manner.

CONCLUSIONS

Our results provide a systematic and comprehensive understanding the function and regulatory mechanism of cotton GELP family, and offer an effective reference for in-depth genetic improvement utilization of cotton GELPs.

Functional Diversification and the Plant Secondary Cell Wall

Much evidence exists suggesting the presence of genetic functional diversification in plants, though literature associated with the role of functional diversification in the evolution of the plant secondary cell wall (SCW) has sparsely been compiled and reviewed in a recent context. This review aims to elucidate, through the examination of gene phylogenies associated with its biosynthesis and maintenance, the role of functional diversification in shaping the critical, dynamic, and characteristic organelle, the secondary cell wall. It will be asserted that gene families resulting from gene duplication and subsequent functional divergence are present and are heavily involved in SCW biosynthesis and maintenance. Furthermore, diversification will be presented as a significant driver behind the evolution of the many functional characteristics of the SCW. The structure and function of the plant cell wall and its constituents will first be explored, followed by a discussion on the phenomenon of gene duplication and the resulting genetic functional divergence that can emerge. Finally, the major constituents of the SCW and their individual relationships with duplication and divergence will be reviewed to the extent of current knowledge on the subject.

Large-Scale Identification and Characterization Analysis of VQ Family Genes in Plants, Especially Gymnosperms

VQ motif-containing (VQ) proteins are a class of transcription regulatory cofactors widely present in plants, playing crucial roles in growth and development, stress response, and defense. Although there have been some reports on the member identification and functional research of VQ genes in some plants, there is still a lack of large-scale identification and clear graphical presentation of their basic characterization information to help us to better understand this family. Especially in gymnosperms, the VQ family genes and their evolutionary relationships have not yet been reported. In this study, we systematically identified 2469 VQ genes from 56 plant species, including bryophytes, gymnosperms, and angiosperms, and analyzed their molecular and evolutionary features. We found that amino acids are only highly conserved in the VQ domain, while other positions are relatively variable; most VQ genes encode relatively small proteins and do not have introns. The GC content in Poaceae plants is the highest (up to 70%); these VQ proteins can be divided into nine subgroups. In particular, we analyzed the molecular characteristics, chromosome distribution, duplication events, and expression levels of VQ genes in three gymnosperms: Ginkgo biloba, Taxus chinensis, and Pinus tabuliformis. In gymnosperms, VQ genes are classified into 11 groups, with highly similar motifs in each group; most VQ proteins have less than 300 amino acids and are predicted to be located in nucleus. Tandem duplication is an important driving force for the expansion of the VQ gene family, and the evolutionary processes of most VQ genes and duplication events are relatively independent; some candidate VQ genes are preliminarily screened, and they are likely to be involved in plant growth and stress and defense responses. These results provide detailed information and powerful references for further understanding and utilizing the VQ family genes in various plants.

Evolutionary Consequences of Functional and Regulatory Divergence of HD-Zip I Transcription Factors as a Source of Diversity in Protein Interaction Networks in Plants

The HD superfamily has been studied in detail for several decades. The plant-specific HD-Zip I subfamily attracts the most attention because of its involvement in plant development and stress responses. In this review, we provide a comprehensive insight into the evolutionary events responsible for the functional redundancy and diversification of the HD-Zip I genes in regulating various biological processes. We summarized the evolutionary history of the HD-Zip family, highlighting the important role of WGDs in its expansion and divergence of retained duplicates in the genome. To determine the relationship between the evolutionary origin and functional conservation of HD-Zip I in different species, we performed a phylogenetic analysis, compared their expression profiles in different tissues and under stress and traced the role of orthologs and paralogs in regulating developmental processes. We found that HD-Zip I from different species have similar gene structures with a highly conserved HD and Zip, bind to the same DNA sequences and are involved in similar biological processes. However, they exhibit a functional diversity, which is manifested in altered expression patterns. Some of them are involved in the regulation of species-specific leaf morphology and phenotypes. Here, we discuss the role of changes in functional domains involved in DNA binding and protein interaction of HD-Zip I and in cis-regulated regions of its target genes in promoting adaptive innovations through the formation of de novo regulatory systems. Understanding the role of the HD-Zip I subfamily in organism-environment interactions remains a challenge for evolutionary developmental biology (evo-devo).

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.

Genome-wide analysis of KIX gene family for organ size regulation in soybean (Glycine max L.)

The KIX domain, conserved among various nuclear and co-activator factors, acts as a binding site that interacts with other transcriptional activators and co-activators, playing a crucial role in gene expression regulation. In plants, the KIX domain is involved in plant hormone signaling, stress response regulation, cell cycle control, and differentiation, indicating its potential relevance to crop productivity. This study aims to identify and characterize KIX domains within the soybean (Glycine max L.) genome to predict their potential role in improving crop productivity. The conservation and evolutionary history of the KIX domains were explored in 59 plant species, confirming the presence of the KIX domains in diverse plants. Specifically, 13 KIX domains were identified within the soybean genome and classified into four main groups, namely GmKIX8/9, GmMED15, GmHAC, and GmRECQL, through sequence alignment, structural analysis, and phylogenetic tree construction. Association analysis was performed between KIX domain haplotypes and soybean seed-related agronomic traits using re-sequencing data from a core collection of 422 accessions. The results revealed correlations between SNP variations observed in GmKIX8-3 and GmMED15-4 and soybean seed phenotypic traits. Additionally, transcriptome analysis confirmed significant expression of the KIX domains during the early stages of soybean seed development. This study provides the first characterization of the structural, expression, genomic haplotype, and molecular features of the KIX domain in soybean, offering a foundation for functional analysis of the KIX domain in soybean and other plants.

Phylogenomic discovery of deleterious mutations facilitates hybrid potato breeding

Hybrid potato breeding will transform the crop from a clonally propagated tetraploid to a seed-reproducing diploid. Historical accumulation of deleterious mutations in potato genomes has hindered the development of elite inbred lines and hybrids. Utilizing a whole-genome phylogeny of 92 Solanaceae and its sister clade species, we employ an evolutionary strategy to identify deleterious mutations. The deep phylogeny reveals the genome-wide landscape of highly constrained sites, comprising ∼2.4% of the genome. Based on a diploid potato diversity panel, we infer 367,499 deleterious variants, of which 50% occur at non-coding and 15% at synonymous sites. Counterintuitively, diploid lines with relatively high homozygous deleterious burden can be better starting material for inbred-line development, despite showing less vigorous growth. Inclusion of inferred deleterious mutations increases genomic-prediction accuracy for yield by 24.7%. Our study generates insights into the genome-wide incidence and properties of deleterious mutations and their far-reaching consequences for breeding.

Multi-Knock-a multi-targeted genome-scale CRISPR toolbox to overcome functional redundancy in plants

Plant genomes are characterized by large and complex gene families that often result in similar and partially overlapping functions. This genetic redundancy severely hampers current efforts to uncover novel phenotypes, delaying basic genetic research and breeding programmes. Here we describe the development and validation of Multi-Knock, a genome-scale clustered regularly interspaced short palindromic repeat toolbox that overcomes functional redundancy in Arabidopsis by simultaneously targeting multiple gene-family members, thus identifying genetically hidden components. We computationally designed 59,129 optimal single-guide RNAs that each target two to ten genes within a family at once. Furthermore, partitioning the library into ten sublibraries directed towards a different functional group allows flexible and targeted genetic screens. From the 5,635 single-guide RNAs targeting the plant transportome, we generated over 3,500 independent Arabidopsis lines that allowed us to identify and characterize the first known cytokinin tonoplast-localized transporters in plants. With the ability to overcome functional redundancy in plants at the genome-scale level, the developed strategy can be readily deployed by scientists and breeders for basic research and to expedite breeding efforts.

The origin and evolution of carpels and fruits from an evo-devo perspective

The flower is an evolutionary innovation in angiosperms that drives the evolution of biodiversity. The carpel is integral to a flower and develops into fruits after fertilization, while the perianth, consisting of the calyx and corolla, is decorative to facilitate pollination and protect the internal organs, including the carpels and stamens. Therefore, the nature of flower origin is carpel and stamen origin, which represents one of the greatest and fundamental unresolved issues in plant evolutionary biology. Here, we briefly summarize the main progress and key genes identified for understanding floral development, focusing on the origin and development of the carpels. Floral ABC models have played pioneering roles in elucidating flower development, but remain insufficient for resolving flower and carpel origin. The genetic basis for carpel origin and subsequent diversification leading to fruit diversity also remains elusive. Based on current research progress and technological advances, simplified floral models and integrative evolutionary-developmental (evo-devo) strategies are proposed for elucidating the genetics of carpel origin and fruit evolution. Stepwise birth of a few master regulatory genes and subsequent functional diversification might play a pivotal role in these evolutionary processes. Among the identified transcription factors, AGAMOUS (AG) and CRABS CLAW (CRC) may be the two core regulatory genes for carpel origin as they determine carpel organ identity, determinacy, and functionality. Therefore, a comparative identification of their protein-protein interactions and downstream target genes between flowering and non-flowering plants from an evo-devo perspective may be primary projects for elucidating carpel origin and development.

Genome-Wide Analysis of ATP Binding Cassette (ABC) Transporters in Peach (Prunus persica) and Identification of a Gene PpABCC1 Involved in Anthocyanin Accumulation

The ATP-binding cassette (ABC) transporter family is a large and diverse protein superfamily that plays various roles in plant growth and development. Although the ABC transporters are known to aid in the transport of a wide range of substrates across biological membranes, their role in anthocyanin transport remains elusive. In this study, we identified a total of 132 putative ABC genes in the peach genome, and they were phylogenetically classified into eight subfamilies. Variations in spatial and temporal gene expression levels resulted in differential expression patterns of PpABC family members in various tissues of peach. PpABCC1 was identified as the most likely candidate gene essential for anthocyanin accumulation in peach. Transient overexpression of PpABCC1 caused a significant increase in anthocyanin accumulation in tobacco leaves and peach fruit, whereas virus-induced gene silencing of PpABCC1 in the blood-fleshed peach resulted in a significant decrease in anthocyanin accumulation. The PpABCC1 promoter contained an MYB binding cis-element, and it could be activated by anthocyanin-activator PpMYB10.1 based on yeast one-hybrid and dual luciferase assays. Thus, it seems that PpABCC1 plays a crucial role in anthocyanin accumulation in peach. Our results provide a new insight into the vacuolar transport of anthocyanins in peach.

Predicted landscape of RETINOBLASTOMA-RELATED LxCxE-mediated interactions across the Chloroplastida

The colonization of land by a single streptophyte algae lineage some 450 million years ago has been linked to multiple key innovations such as three-dimensional growth, alternation of generations, the presence of stomata, as well as innovations inherent to the birth of major plant lineages, such as the origins of vascular tissues, roots, seeds and flowers. Multicellularity, which evolved multiple times in the Chloroplastida coupled with precise spatiotemporal control of proliferation and differentiation were instrumental for the evolution of these traits. RETINOBLASTOMA-RELATED (RBR), the plant homolog of the metazoan Retinoblastoma protein (pRB), is a highly conserved and multifunctional core cell cycle regulator that has been implicated in the evolution of multicellularity in the green lineage as well as in plant multicellularity-related processes such as proliferation, differentiation, stem cell regulation and asymmetric cell division. RBR fulfills these roles through context-specific protein-protein interactions with proteins containing the Leu-x-Cys-x-Glu (LxCxE) short-linear motif (SLiM); however, how RBR-LxCxE interactions have changed throughout major innovations in the Viridiplantae kingdom is a question that remains unexplored. Here, we employ an in silico evo-devo approach to predict and analyze potential RBR-LxCxE interactions in different representative species of key Chloroplastida lineages, providing a valuable resource for deciphering RBR-LxCxE multiple functions. Furthermore, our analyses suggest that RBR-LxCxE interactions are an important component of RBR functions and that interactions with chromatin modifiers/remodelers, DNA replication and repair machinery are highly conserved throughout the Viridiplantae, while LxCxE interactions with transcriptional regulators likely diversified throughout the water-to-land transition.

The role of CDPKs in plant development, nutrient and stress signaling

The second messenger calcium (Ca2+) is a ubiquitous intracellular signaling molecule found in eukaryotic cells. In plants, the multigene family of calcium-dependent protein kinases (CDPKs) plays an important role in regulating plant growth, development, and stress tolerance. CDPKs sense changes in intracellular Ca2+ concentration and translate them into phosphorylation events that initiate downstream signaling processes. Several functional and expression studies on different CDPKs and their encoding genes have confirmed their multifunctional role in stress. Here, we provide an overview of the signal transduction mechanisms and functional roles of CDPKs. This review includes details on the regulation of secondary metabolites, nutrient uptake, regulation of flower development, hormonal regulation, and biotic and abiotic stress responses.

Genome-Wide Analysis of the Type-B Authentic Response Regulator Gene Family in Brassica napus

The type-B authentic response regulators (type-B ARRs) are positive regulators of cytokinin signaling and involved in plant growth and stress responses. In this study, we used bioinformatics, RNA-seq, and qPCR to study the phylogenetic and expression pattern of 35 type-B ARRs in Brassica napus. The BnARRs experienced gene expansion and loss during genome polyploidization and were classified into seven groups. Whole-genome duplication (WGD) and segmental duplication were the main forces driving type-B ARR expansion in B. napus. Several BnARRs with specific expression patterns during rapeseed development were identified, including BnARR12/14/18/23/33. Moreover, we found the type-B BnARRs were involved in rapeseed development and stress responses, through participating in cytokinin and ABA signaling pathways. This study revealed the origin, evolutionary history, and expression pattern of type-B ARRs in B. napus and will be helpful to the functional characterization of BnARRs.

Genome-wide identification and characterization of the bHLH gene family and analysis of their potential relevance to chlorophyll metabolism in Raphanus sativus L

BACKGROUND

Green-fleshed radish (Raphanus sativus L.) is an economically important root vegetable of the Brassicaceae family, and chlorophyll accumulates in its root tissues. It was reported that the basic helix-loop-helix (bHLH) transcription factors play vital roles in the process of chlorophyll metabolism. Nevertheless, a comprehensive study on the bHLH gene family has not been performed in Raphanus sativus L.

RESULTS

In this study, a total of 213 Raphanus sativus L. bHLH (RsbHLH) genes were screened in the radish genome, which were grouped into 22 subfamilies. 204 RsbHLH genes were unevenly distributed on nine chromosomes, and nine RsbHLH genes were located on the scaffolds. Gene structure analysis showed that 25 RsbHLH genes were intron-less. Collineation analysis revealed the syntenic orthologous bHLH gene pairs between radish and Arabidopsis thaliana/Brassica rapa/Brassica oleracea. 162 RsbHLH genes were duplicated and retained from the whole genome duplication event, indicating that the whole genome duplication contributed to the expansion of the RsbHLH gene family. RNA-seq results revealed that RsbHLH genes had a variety of expression patterns at five development stages of green-fleshed radish and white-fleshed radish. In addition, the weighted gene co-expression network analysis confirmed four RsbHLH genes closely related to chlorophyll content.

CONCLUSIONS

A total of 213 RsbHLH genes were identified, and we systematically analyzed their gene structure, evolutionary and collineation relationships, conserved motifs, gene duplication, cis-regulatory elements and expression patterns. Finally, four bHLH genes closely involved in chlorophyll content were identified, which may be associated with the photosynthesis of the green-fleshed radish. The current study would provide valuable information for further functional exploration of RsbHLH genes, and facilitate clarifying the molecular mechanism underlying photosynthesis process in green-fleshed radish.

Origin, Expansion, and Divergence of ETHYLENE-INSENSITIVE 3 (EIN3)/EIN3-LIKE Transcription Factors During Streptophytes Evolution

The transition of plants to land required several regulatory adaptive mechanisms. Little is known about these mechanisms, but they no doubt involved the evolution of transcription factor (TF) families. ETHYLENE-INSENSITIVE 3 (EIN3)/EIN3-LIKE (EIL) transcription factors (TFs) are core components of the ethylene signaling pathway that play important roles in almost every aspect of plant development and environmental responses by regulating the transcription of numerous genes. However, the evolutionary history of EIN3/EIL TFs, which are present in a wide range of streptophytes, is still not clear. Here, to explore the evolution and functions of EIN3/EIL TFs, we performed phylogenetic analysis of these TFs and investigated their gene and protein structures as well as sequence features. Our results suggest that the EIN3/EIL TF family was already was already present in the ancestor of streptophyte algae. Phylogenetic analysis divided the EIN3/EIL TFs into three groups (Group A-C). Analysis of gene and protein structure revealed that most of the structural features of these TFs had already formed in ancient lineages. Further investigation suggested that all groups have undergone several duplication events related to whole-genome duplications in plants, generating multiple, functionally diverse gene copies. Therefore, as plants colonized terrestrial habitats and evolved key traits, the EIN3/EIL TF family expanded broadly via multiple duplication events, which could have promoted their fundamental neo- and sub-functionalization to help plants adapt to terrestrial life. Our findings shed light on the functional evolution of the EIN3/EIL TF family in the streptophytes.

High-quality genome and methylomes illustrate features underlying evolutionary success of oaks

The genus Quercus, which emerged ∼55 million years ago during globally warm temperatures, diversified into ∼450 extant species. We present a high-quality de novo genome assembly of a California endemic oak, Quercus lobata, revealing features consistent with oak evolutionary success. Effective population size remained large throughout history despite declining since early Miocene. Analysis of 39,373 mapped protein-coding genes outlined copious duplications consistent with genetic and phenotypic diversity, both by retention of genes created during the ancient γ whole genome hexaploid duplication event and by tandem duplication within families, including numerous resistance genes and a very large block of duplicated DUF247 genes, which have been found to be associated with self-incompatibility in grasses. An additional surprising finding is that subcontext-specific patterns of DNA methylation associated with transposable elements reveal broadly-distributed heterochromatin in intergenic regions, similar to grasses. Collectively, these features promote genetic and phenotypic variation that would facilitate adaptability to changing environments.

Genome-Wide Identification, Expansion, and Evolution Analysis of Homeobox Gene Family Reveals TALE Genes Important for Secondary Cell Wall Biosynthesis in Moso Bamboo (Phyllostachys edulis)

The TALE gene family is a subfamily of the homeobox gene family and has been implicated in regulating plant secondary growth. However, reports about the evolutionary history and function of the TALE gene family in bamboo are limited. Here, the homeobox gene families of moso bamboo Olyra latifolia and Bonia amplexicaulis were identified and compared. Many duplication events and obvious expansions were found in the TALE family of woody bamboo. PhTALEs were found to have high syntenies with TALE genes in rice. Through gene co-expression analysis and quantitative real-time PCR analysis, the candidate PhTALEs were thought to be involved in regulating secondary cell wall development of moso bamboo during the fast-growing stage. Among these candidate PhTALEs, orthologs of OsKNAT7, OSH15, and SH5 in moso bamboo may regulate xylan synthesis by regulating the expression of IRX-like genes. These results suggested that PhTALEs may participate in the secondary cell wall deposition in internodes during the fast-growing stage of moso bamboo. The expansion of the TALE gene family may be implicated in the increased lignification of woody bamboo when divergent from herbaceous bamboos.

Genome-wide Identification and Expression Analysis of RcMYB Genes in Rhodiola crenulata

Modern research has proved that the main medicinal component of Rhodiola crenulata, which has a wide range of medicinal value, is its secondary metabolite salidroside. The MYB transcription factor family is widely involved in biosynthesis of second metabolism and other roles in the stress response in plants, so a genome-wide identification and analysis for this family in R. crenulata is worth conducting. In this research, genome-wide analysis identified 139 MYB genes based on conserved domains in the R. crenulata genome, and 137 genes were used to construct a phylogenetic tree and modified with expression files to reveal evolutionary characteristics. Physical and chemical characteristics, gene structure, and conserved motif analysis were also used to further analyze RcMYBs. Additionally, cis-acting elements related to transcription, hormone, and MYB binding were found in the promoter region of the selected RcMYBs. Four RcMYBs were cloned, sequenced, and their gene expression pattern was analyzed for further analysis of their functions. The research results lay the foundation for further research on the function of RcMYB and R. crenulata.

Genome-Wide Identification and Expression Profiling of CONSTANS-Like Genes in Pepper (Capsicum annuum): Gaining an Insight to Their Phylogenetic Evolution and Stress-Specific Roles

CONSTANS-like (COL) genes play important regulatory roles in multiple growth and development processes of plants but have rarely been studied in Capsicum annuum. This study explored the evolutionary relationship and expression patterns of COL genes from C. annuum. A total of 10 COL genes were identified in the genome of the cultivated pepper Zunla-1 and were named CaCOL01-10. These genes were unequally distributed among five chromosomes and could be divided into three groups based on differences in gene structure characteristics. During evolutionary history, duplications and retentions were divergent among different groups of COL genes. Tandem duplication caused amplification of group I genes. Genetic distance among COL genes was the largest in group III, suggesting that group III genes undergo more relaxed selection pressure compared with the other groups. Expression patterns of CaCOLs in tissues were significantly different, with CaCOL08 exhibiting the highest expression in stem and leaf. Some COL orthologous genes showed markedly different expression patterns in pepper compared with tomato, such as COL_1 orthologs, which may be involved in fruit development in pepper. In addition, CaCOLs participated in the regulation of abiotic stresses to varying degrees. Five CaCOL genes were induced by cold, and CaCOL02 and CaCOL03 were specifically upregulated by cold and downregulated by heat. This study provides a theoretical basis for the in-depth understanding of the functions of COL genes in pepper and their molecular mechanisms involved in growth and development and responses to abiotic stresses.

Evolution of a novel cell type in Dictyostelia required gene duplication of a cudA-like transcription factor

The evolution of novel cell types has been proposed to result from duplication of gene regulatory networks, but proven examples are rare. In addition to stalk cells and spores that make up the fruiting bodies of three major groups of Dictyostelia, those in group 4 additionally evolved basal disc and cup cells that respectively anchor the stalk to the substratum and the spore mass to the stalk. We noted a putative group-4-specific duplication of a cudA-like transcription factor (TF) in a comparative analysis of group-representative genomes. Using increased taxon sampling, we here confirmed that this TF, cdl1, duplicated into cdl1a and cdl1b in the common ancestor to group 4. cdl1a, but not cdl1b, showed signatures of positive selection, indicative of functional innovation. Deletion of cdl1a in Dictyostelium discoideum resulted in fruiting bodies with sagging spore heads that lacked the supporting cup cells and expression of cup-specific genes. Deletion of cdl1b resulted in thinner fruiting body stalks, while a cdl1b-cdl1a- double knockout showed more severe stalk defects, suggesting an ancestral role of cdl1 in stalk formation. This was confirmed in a closely related non-group 4 species, Polysphondylium violaceum, where cdl1 knockout caused defective stalk formation. These data indicate that the group-specific duplication of cdl1 and subsequent diversification of cdl1a played a pivotal role in the evolution of a novel somatic cell type in group 4 Dictyostelia.

Genome-wide identification and expression analysis of the VQ gene family in Cucurbita pepo L

VQ protein is a plant specific protein, which plays an important role in plant growth and development and biological and abiotic stress response. This study aimed to systematically analyze for the first time the VQ of Cucurbita pepo and understand their expression patterns in response to different stimuli. Herein, 44 VQ genes were identified, which were divided into eight groups (I-VIII) based on phylogenetic analysis. Two genes (CpVQ1 and CpVQ2) could not be located on the chromosome, whereas the remaining CpVQ genes were randomly distributed on the chromosomes, except for chromosomes 15 and 18. Noteworthy, the main event driving the expansion of the VQ gene family was chromosome fragment duplication. Based on qRT-PCR analysis, VQ genes are expressed in different tissues, and VQ genes are differentially regulated under a variety of abiotic stresses and powdery mildew stress, indicating that they play an important role in plant stress response and other aspects. This report presents the first systematic analysis of VQ genes from C. pepo and provides a solid foundation for further research of the specific functions of VQ proteins.

Genomic Analysis of the Glutathione S-Transferase Family in Pear (Pyrus communis) and Functional Identification of PcGST57 in Anthocyanin Accumulation

Anthocyanin accumulation in vacuoles results in red coloration in pear peels. Glutathione S-transferase (GST) proteins have emerged as important regulators of anthocyanin accumulation. Here, a total of 57 PcGST genes were identified in the European pear ‘Bartlett’ (Pyrus communis) through comprehensive genomic analysis. Phylogenetic analysis showed that PcGST genes were divided into 10 subfamilies. The gene structure, chromosomal localization, collinearity relationship, cis-elements in the promoter region, and conserved motifs of PcGST genes were analyzed. Further research indicated that glutamic acid (Glu) can significantly improve anthocyanin accumulation in pear peels. RNA sequencing (RNA-seq) analysis showed that Glu induced the expression of most PcGST genes, among which PcGST57 was most significantly induced. Further phylogenetic analysis indicated that PcGST57 was closely related to GST genes identified in other species, which were involved in anthocyanin accumulation. Transcript analysis indicated that PcGST57 was expressed in various tissues, other than flesh, and associated with peel coloration at different developmental stages. Silencing of PcGST57 by virus-induced gene silencing (VIGS) inhibited the expression of PcGST57 and reduced the anthocyanin content in pear fruit. In contrast, overexpression of PcGST57 improved anthocyanin accumulation. Collectively, our results demonstrated that PcGST57 was involved in anthocyanin accumulation in pear and provided candidate genes for red pear breeding.

Genomic and Meiotic Changes Accompanying Polyploidization

Hybridization and polyploidy have been considered as significant evolutionary forces in adaptation and speciation, especially among plants. Interspecific gene flow generates novel genetic variants adaptable to different environments, but it is also a gene introgression mechanism in crops to increase their agronomical yield. An estimate of 9% of interspecific hybridization has been reported although the frequency varies among taxa. Homoploid hybrid speciation is rare compared to allopolyploidy. Chromosome doubling after hybridization is the result of cellular defects produced mainly during meiosis. Unreduced gametes, which are formed at an average frequency of 2.52% across species, are the result of altered spindle organization or orientation, disturbed kinetochore functioning, abnormal cytokinesis, or loss of any meiotic division. Meiotic changes and their genetic basis, leading to the cytological diploidization of allopolyploids, are just beginning to be understood especially in wheat. However, the nature and mode of action of homoeologous recombination suppressor genes are poorly understood in other allopolyploids. The merger of two independent genomes causes a deep modification of their architecture, gene expression, and molecular interactions leading to the phenotype. We provide an overview of genomic changes and transcriptomic modifications that particularly occur at the early stages of allopolyploid formation.

Identification and Comprehensive Genome-Wide Analysis of Glutathione S-Transferase Gene Family in Sweet Cherry (Prunus avium) and Their Expression Profiling Reveals a Likely Role in Anthocyanin Accumulation

Glutathione S-transferases (GSTs) in plants are multipurpose enzymes that are involved in growth and development and anthocyanins transportation. However, members of the GST gene family were not identified in sweet cherry (Prunus avium). To identify the GST genes in sweet cherry, a genome-wide analysis was conducted. In this study, we identified 67 GST genes in P. avium genome and nomenclature according to chromosomal distribution. Phylogenetic tree analysis revealed that PavGST genes were classified into seven chief subfamily: TCHQD, Theta, Phi, Zeta, Lambda, DHAR, and Tau. The majority of the PavGST genes had a relatively well-maintained exon-intron and motif arrangement within the same group, according to gene structure and motif analyses. Gene structure (introns-exons) and conserved motif analysis revealed that the majority of the PavGST genes showed a relatively well-maintained motif and exons-introns configuration within the same group. The chromosomal localization, GO enrichment annotation, subcellular localization, syntenic relationship, Ka/Ks analysis, and molecular characteristics were accomplished using various bioinformatics tools. Mode of gene duplication showed that dispersed duplication might play a key role in the expansion of PavGST gene family. Promoter regions of PavGST genes contain numerous cis-regulatory components, which are involved in multiple stress responses, such as abiotic stress and phytohormones responsive factors. Furthermore, the expression profile of sweet cherry PavGSTs showed significant results under LED treatment. Our findings provide the groundwork for future research into induced LED anthocyanin and antioxidants deposition in sweet cherries.

Comprehensive Analysis of the SUV Gene Family in Allopolyploid Brassica napus and Its Diploid Ancestors

SUV (the Suppressor of variegation [Su(var)] homologs and related) gene family is a subgroup of the SET gene family. According to the SRA domain and WIYLD domain distributions, it can be divided into two categories, namely SUVH (the Suppressor of variegation [Su(var)] homologs) and SUVR (the Suppressor of variegation [Su(var)] related). In this study, 139 SUV genes were identified in allopolyploid Brassica napus and its diploid ancestors, and their evolutionary relationships, protein properties, gene structures, motif distributions, transposable elements, cis-acting elements and gene expression patterns were analyzed. Our results showed that the SUV gene family of B. napus was amplified during allopolyploidization, in which the segmental duplication and TRD played critical roles. After the separation of Brassica and Arabidopsis lineages, orthologous gene analysis showed that many SUV genes were lost during the evolutionary process in B. rapa, B. oleracea and B. napus. The analysis of the gene and protein structures and expression patterns of 30 orthologous gene pairs which may have evolutionary relationships showed that most of them were conserved in gene structures and protein motifs, but only four gene pairs had the same expression patterns.

Two Is Company, but Four Is a Party-Challenges of Tetraploidization for Cell Wall Dynamics and Efficient Tip-Growth in Pollen

Some cells grow by an intricately coordinated process called tip-growth, which allows the formation of long tubular structures by a remarkable increase in cell surface-to-volume ratio and cell expansion across vast distances. On a broad evolutionary scale, tip-growth has been extraordinarily successful, as indicated by its recurrent 're-discovery' throughout evolutionary time in all major land plant taxa which allowed for the functional diversification of tip-growing cell types across gametophytic and sporophytic life-phases. All major land plant lineages have experienced (recurrent) polyploidization events and subsequent re-diploidization that may have positively contributed to plant adaptive evolutionary processes. How individual cells respond to genome-doubling on a shorter evolutionary scale has not been addressed as elaborately. Nevertheless, it is clear that when polyploids first form, they face numerous important challenges that must be overcome for lineages to persist. Evidence in the literature suggests that tip-growth is one of those processes. Here, I discuss the literature to present hypotheses about how polyploidization events may challenge efficient tip-growth and strategies which may overcome them: I first review the complex and multi-layered processes by which tip-growing cells maintain their cell wall integrity and steady growth. I will then discuss how they may be affected by the cellular changes that accompany genome-doubling. Finally, I will depict possible mechanisms polyploid plants may evolve to compensate for the effects caused by genome-doubling to regain diploid-like growth, particularly focusing on cell wall dynamics and the subcellular machinery they are controlled by.

Gene duplication drove the loss of awn in sorghum

Loss of the awn in some cereals, including sorghum, is a key transition during cereal domestication or improvement that has facilitated grain harvest and storage. The genetic basis of awn loss in sorghum during domestication or improvement remains unknown. Here, we identified the awn1 gene encoding a transcription factor with the ALOG domain that is responsible for awn loss during sorghum domestication or improvement. awn1 arose from a gene duplication on chromosome 10 that translocated to chromosome 3, recruiting a new promoter from the neighboring intergenic region filled with "noncoding DNA" and recreating the first exon and intron. awn1 acquired high expression after duplication and represses the elongation of awns in domesticated sorghum. Comparative mapping revealed high collinearity at the awn1 paralog locus on chromosome 10 across cereals, and awn growth and development were successfully reactivated on the rice spikelet by inactivating the rice awn1 ortholog. RNA-seq and DAP-seq revealed that as a transcriptional repressor, AWN1 bound directly to a motif in the regulatory regions of three MADS genes related to flower development and two genes, DL and LKS2, involved in awn development. AWN1 downregulates the expression of these genes, thereby repressing awn elongation. The preexistence of regulatory elements in the neighboring intergenic region of awn1 before domestication implicates that noncoding DNA may serve as a treasure trove for evolution during sorghum adaptation to a changing world. Taken together, our results suggest that gene duplication can rapidly drive the evolution of gene regulatory networks in plants.

Ectopic expression of OsJAZ6, which interacts with OsJAZ1, alters JA signaling and spikelet development in rice

Jasmonates (JAs) are key phytohormones that regulate plant responses and development. JASMONATE-ZIM DOMAIN (JAZ) proteins safeguard JA signaling by repressing JA-responsive gene expression in the absence of JA. However, the interaction and cooperative roles of JAZ repressors remain unclear during plant development. Here, we found that OsJAZ6 interacts with OsJAZ1 depending on a single amino acid in the so-called ZIM domain of OsJAZ6 in rice JA signaling transduction and JA-regulated rice spikelet development. In vivo protein distribution analysis revealed that the OsJAZ6 content is efficiently regulated during spikelet development, and biochemical and genetic evidence showed that OsJAZ6 is more sensitive to JA-mediated degradation than OsJAZ1. Through over- and mis-expression experiments, we further showed that the protein stability and levels of OsJAZ6 orchestrate the output of JA signaling during rice spikelet development. A possible mechanism, which outlines how OsJAZ repressors interact and function synergistically in specifying JA signaling output through degradation titration, is also discussed.

Physalis floridana CRABS CLAW mediates neofunctionalization of GLOBOSA genes in carpel development

Floral B-function MADS-box genes, such as GLOBOSA (GLO), function in corolla and stamen organ identity specification. The functions of these genes outside these floral whorls are rarely reported. DOLL1 is a GLO gene controlling corolla and androecium organ identity. In this study we found that, in Physalis floridana double-layered-lantern 1 (doll1) mutant pollinated with wild-type pollen, fruit set was extremely low, indicating that doll1 females are dysfunctional. Stigma and style structure, stigma receptivity, pollen tube guidance, and embryo sac development were also impaired in doll1. P. floridana CRABS CLAW (PFCRC), predominantly expressed in carpels, was repressed in doll1 native carpels. Loss-of-function of PFCRC altered carpel meristem determinacy, carpel closure, and ovule number, and the resultant 'pistil' consisted of multiple spirally-arranged dorsiventral carpels occasionally with 1-2 naked ovules on the margin and trichomes at each mutated carpel tip, implying an alteration of carpel organ identity. Regulatory and genetic interactions between B-class MADS-box genes and PFCRC were revealed in a context-dependent manner in floral development. Our work reveals a new role for the B-function genes in carpel and ovule development via regulating PFCRC, providing a new understanding of genetic regulatory networks between MADS-domain and CRC transcription factors in mediating carpel organ specification, functionality, and origin.

Lipid Transfer Proteins (LTPs)-Structure, Diversity and Roles beyond Antimicrobial Activity

Lipid transfer proteins (LTPs) are among the most promising plant-exclusive antimicrobial peptides (AMPs). They figure among the most challenging AMPs from the point of view of their structural diversity, functions and biotechnological applications. This review presents a current picture of the LTP research, addressing not only their structural, evolutionary and further predicted functional aspects. Traditionally, LTPs have been identified by their direct isolation by biochemical techniques, whereas omics data and bioinformatics deserve special attention for their potential to bring new insights. In this context, new possible functions have been identified revealing that LTPs are actually multipurpose, with many additional predicted roles. Despite some challenges due to the toxicity and allergenicity of LTPs, a systematic review and search in patent databases, indicate promising perspectives for the biotechnological use of LTPs in human health and also plant defense.

Genome-wide identification of ascorbate-glutathione cycle gene families in soybean (Glycine max) reveals gene duplication events and specificity of gene members linked to development and stress conditions

Ascorbate-glutathione (AsA-GSH) cycle plays an important role in tuning beneficial ROS accumulation for intracellular signals and imparts plant tolerance to oxidative stress by detoxifying excess of ROS. Here, we present genome-wide identification of AsA-GSH cycle genes (APX, MDHAR, DHAR, and GR) in several leguminous species and expression analyses in G. max during stress, germination and tissue development. Our data revealed 24 genes in Glycine genus against the maximum of 15 in other leguminous species, which was due to 9 pars of duplicated genes mostly originated from sub/neofunctionalization. Cytosolic APX and MDHAR genes were highly expressed in different tissues and physiological conditions. Germination induced genes encoding AsA-GSH proteins from different cell compartments, whereas vegetative phase (leaves) stimulated predominantly genes related to chloroplast/mitochondria proteins. Moreover, cytosolic APX-1, 2, MDHAR-1a, 1b and GR genes were the primary genes linked to senescence and biotic stresses, while stAPX-a, b and GR (from organelles) were the most abiotic stress related genes. Biotic and abiotic stress tolerant genotypes generally showed increased MDHAR, DHAR and/or GR mRNA levels compared to susceptible genotypes. Overall, these data clarified evolutionary events in leguminous plants and point to the functional specificity of duplicated genes of the AsA-GSH cycle in G. max.

Ancient gene duplications, rather than polyploidization, facilitate diversification of petal pigmentation patterns in Clarkia gracilis (Onagraceae)

It has been suggested that gene duplication and polyploidization create opportunities for the evolution of novel characters. However, the connections between the effects of polyploidization and morphological novelties have rarely been examined. In this study, we investigated whether petal pigmentation patterning in an allotetraploid Clarkia gracilis has evolved as a result of polyploidization. C. gracilis is thought to be derived through a recent polyploidization event with two diploid species, C. amoena huntiana and an extinct species that is closely related to C. lassenensis. We reconstructed phylogenetic relationships of the R2R3-MYBs (the regulators of petal pigmentation) from two subspecies of C. gracilis and the two purported progenitors, C. a. huntiana and C. lassenensis. The gene tree reveals that these R2R3-MYB genes have arisen through duplications that occurred before the divergence of the two progenitor species, i.e., before polyploidization. After polyploidization and subsequent gene loss, only one of the two orthologous copies inherited from the progenitors was retained in the polyploid, turning it to diploid inheritance. We examined evolutionary changes in these R2R3-MYBs and in their expression, which reveals that the changes affecting patterning (including expression domain contraction, loss-of-function mutation, cis-regulatory mutation) occurred after polyploidization within the C. gracilis lineages. Our results thus suggest that polyploidization itself is not necessary in producing novel petal color patterns. By contrast, duplications of R2R3-MYB genes in the common ancestor of the two progenitors have apparently facilitated diversification of petal pigmentation patterns.

The evolution of the phenylpropanoid pathway entailed pronounced radiations and divergences of enzyme families

Land plants constantly respond to fluctuations in their environment. Part of their response is the production of a diverse repertoire of specialized metabolites. One of the foremost sources for metabolites relevant to environmental responses is the phenylpropanoid pathway, which was long thought to be a land-plant-specific adaptation shaped by selective forces in the terrestrial habitat. Recent data have, however, revealed that streptophyte algae, the algal relatives of land plants, have candidates for the genetic toolkit for phenylpropanoid biosynthesis and produce phenylpropanoid-derived metabolites. Using phylogenetic and sequence analyses, we here show that the enzyme families that orchestrate pivotal steps in phenylpropanoid biosynthesis have independently undergone pronounced radiations and divergence in multiple lineages of major groups of land plants; sister to many of these radiated gene families are streptophyte algal candidates for these enzymes. These radiations suggest a high evolutionary versatility in the enzyme families involved in the phenylpropanoid-derived metabolism across embryophytes. We suggest that this versatility likely translates into functional divergence, and may explain the key to one of the defining traits of embryophytes: a rich specialized metabolism.

Gene duplications and phylogenomic conflict underlie major pulses of phenotypic evolution in gymnosperms

Inferring the intrinsic and extrinsic drivers of species diversification and phenotypic disparity across the tree of life is a major challenge in evolutionary biology. In green plants, polyploidy (or whole-genome duplication, WGD) is known to play a major role in microevolution and speciation, but the extent to which WGD has shaped macroevolutionary patterns of diversification and phenotypic innovation across plant phylogeny remains an open question. Here, we examine the relationship of various facets of genomic evolution-including gene and genome duplication, genome size, and chromosome number-with macroevolutionary patterns of phenotypic innovation, species diversification, and climatic occupancy in gymnosperms. We show that genomic changes, such as WGD and genome-size shifts, underlie the origins of most major extant gymnosperm clades, and notably, our results support an ancestral WGD in the gymnosperm lineage. Spikes of gene duplication typically coincide with major spikes of phenotypic innovation, while increased rates of phenotypic evolution are typically found at nodes with high gene-tree conflict, representing historic population-level dynamics during speciation. Most shifts in gymnosperm diversification since the rise of angiosperms are decoupled from putative WGDs and instead are associated with increased rates of climatic occupancy evolution, particularly in cooler and/or more arid climatic conditions, suggesting that ecological opportunity, especially in the later Cenozoic, and environmental heterogeneity have driven a resurgence of gymnosperm diversification. Our study provides critical insight on the processes underlying diversification and phenotypic evolution in gymnosperms, with important broader implications for the major drivers of both micro- and macroevolution in plants.

A plant-like mechanism coupling m6A reading to polyadenylation safeguards transcriptome integrity and developmental gene partitioning in Toxoplasma

Correct 3’end processing of mRNAs is one of the regulatory cornerstones of gene expression. In a parasite that must adapt to the regulatory requirements of its multi-host life style, there is a need to adopt additional means to partition the distinct transcriptional signatures of the closely and tandemly arranged stage-specific genes. In this study, we report our findings in T. gondii of an m6A-dependent 3’end polyadenylation serving as a transcriptional barrier at these loci. We identify the core polyadenylation complex within T. gondii and establish CPSF4 as a reader for m6A-modified mRNAs, via a YTH domain within its C-terminus, a feature which is shared with plants. We bring evidence of the specificity of this interaction both biochemically, and by determining the crystal structure at high resolution of the T. gondii CPSF4-YTH in complex with an m6A-modified RNA. We show that the loss of m6A, both at the level of its deposition or its recognition is associated with an increase in aberrantly elongated chimeric mRNAs emanating from impaired transcriptional termination, a phenotype previously noticed in the plant model Arabidopsis thaliana. Nanopore direct RNA sequencing shows the occurrence of transcriptional read-through breaching into downstream repressed stage-specific genes, in the absence of either CPSF4 or the m6A RNA methylase components in both T. gondii and A. thaliana. Taken together, our results shed light on an essential regulatory mechanism coupling the pathways of m6A metabolism directly to the cleavage and polyadenylation processes, one that interestingly seem to serve, in both T. gondii and A. thaliana, as a guardian against aberrant transcriptional read-throughs.

Signatures of Transcription Factor Evolution and the Secondary Gain of Red Algae Complexity

Red algae (Rhodophyta) belong to the superphylum Archaeplastida, and are a species-rich group exhibiting diverse morphologies. Theory has it that the unicellular red algal ancestor went through a phase of genome contraction caused by adaptation to extreme environments. More recently, the classes Porphyridiophyceae, Bangiophyceae, and Florideophyceae experienced genome expansions, coinciding with an increase in morphological complexity. Transcription-associated proteins (TAPs) regulate transcription, show lineage-specific patterns, and are related to organismal complexity. To better understand red algal TAP complexity and evolution, we investigated the TAP family complement of uni- and multi-cellular red algae. We found that the TAP family complement correlates with gain of morphological complexity in the multicellular Bangiophyceae and Florideophyceae, and that abundance of the C2H2 zinc finger transcription factor family may be associated with the acquisition of morphological complexity. An expansion of heat shock transcription factors (HSF) occurred within the unicellular Cyanidiales, potentially as an adaption to extreme environmental conditions.

A New Pipeline for Removing Paralogs in Target Enrichment Data

Target enrichment (such as Hyb-Seq) is a well-established high throughput sequencing method that has been increasingly used for phylogenomic studies. Unfortunately, current widely used pipelines for analysis of target enrichment data do not have a vigorous procedure to remove paralogs in target enrichment data. In this study, we develop a pipeline we call Putative Paralogs Detection (PPD) to better address putative paralogs from enrichment data. The new pipeline is an add-on to the existing HybPiper pipeline, and the entire pipeline applies criteria in both sequence similarity and heterozygous sites at each locus in the identification of paralogs. Users may adjust the thresholds of sequence identity and heterozygous sites to identify and remove paralogs according to the level of phylogenetic divergence of their group of interest. The new pipeline also removes highly polymorphic sites attributed to errors in sequence assembly and gappy regions in the alignment. We demonstrated the value of the new pipeline using empirical data generated from Hyb-Seq and the Angiosperms353 kit for two woody genera Castanea (Fagaceae, Fagales) and Hamamelis (Hamamelidaceae, Saxifragales). Comparisons of data sets showed that the PPD identified many more putative paralogs than the popular method HybPiper. Comparisons of tree topologies and divergence times showed evident differences between data from HybPiper and data from our new PPD pipeline. We further evaluated the accuracy and error rates of PPD by BLAST mapping of putative paralogous and orthologous sequences to a reference genome sequence of Castanea mollissima. Compared to HybPiper alone, PPD identified substantially more paralogous gene sequences that mapped to multiple regions of the reference genome (31 genes for PPD compared with 4 genes for HybPiper alone). In conjunction with HybPiper, paralogous genes identified by both pipelines can be removed resulting in the construction of more robust orthologous gene data sets for phylogenomic and divergence time analyses. Our study demonstrates the value of Hyb-Seq with data derived from the Angiosperms353 probe set for elucidating species relationships within a genus, and argues for the importance of additional steps to filter paralogous genes and poorly aligned regions (e.g., as occur through assembly errors), such as our new PPD pipeline described in this study. [Angiosperms353; Castanea; divergence time; Hamamelis; Hyb-Seq, paralogs, phylogenomics.]

HAG1 and SWI3A/B control of male germ line development in P. patens suggests conservation of epigenetic reproductive control across land plants

KEY MESSAGE

Bryophytes as models to study the male germ line: loss-of-function mutants of epigenetic regulators HAG1 and SWI3a/b demonstrate conserved function in sexual reproduction. With the water-to-land transition, land plants evolved a peculiar haplodiplontic life cycle in which both the haploid gametophyte and the diploid sporophyte are multicellular. The switch between these phases was coined alternation of generations. Several key regulators that control the bauplan of either generation are already known. Analyses of such regulators in flowering plants are difficult due to the highly reduced gametophytic generation, and the fact that loss of function of such genes often is embryo lethal in homozygous plants. Here we set out to determine gene function and conservation via studies in bryophytes. Bryophytes are sister to vascular plants and hence allow evolutionary inferences. Moreover, embryo lethal mutants can be grown and vegetatively propagated due to the dominance of the bryophyte gametophytic generation. We determined candidates by selecting single copy orthologs that are involved in transcriptional control, and of which flowering plant mutants show defects during sexual reproduction, with a focus on the under-studied male germ line. We selected two orthologs, SWI3a/b and HAG1, and analyzed loss-of-function mutants in the moss P. patens. In both mutants, due to lack of fertile spermatozoids, fertilization and hence the switch to the diploid generation do not occur. Pphag1 additionally shows arrested male and impaired female gametangia development. We analyzed HAG1 in the dioecious liverwort M. polymorpha and found that in Mphag1 the development of gametangiophores is impaired. Taken together, we find that involvement of both regulators in sexual reproduction is conserved since the earliest divergence of land plants.

Novel Structural Variation and Evolutionary Characteristics of Chloroplast tRNA in Gossypium Plants

Cotton is one of the most important fiber and oil crops in the world. Chloroplast genomes harbor their own genetic materials and are considered to be highly conserved. Transfer RNAs (tRNAs) act as "bridges" in protein synthesis by carrying amino acids. Currently, the variation and evolutionary characteristics of tRNAs in the cotton chloroplast genome are poorly understood. Here, we analyzed the structural variation and evolution of chloroplast tRNA (cp tRNA) based on eight diploid and two allotetraploid cotton species. We also investigated the nucleotide evolution of chloroplast genomes in cotton species. We found that cp tRNAs in cotton encoded 36 or 37 tRNAs, and 28 or 29 anti-codon types with lengths ranging from 60 to 93 nucleotides. Cotton chloroplast tRNA sequences possessed specific conservation and, in particular, the Ψ-loop contained the conserved U-U-C-X3-U. The cp tRNAs of Gossypium L. contained introns, and cp tRNA^Ile contained the anti-codon (C-A-U), which was generally the anti-codon of tRNA^Met. The transition and transversion analyses showed that cp tRNAs in cotton species were iso-acceptor specific and had undergone unequal rates of evolution. The intergenic region was more variable than coding regions, and non-synonymous mutations have been fixed in cotton cp genomes. On the other hand, phylogeny analyses indicated that cp tRNAs of cotton were derived from several inferred ancestors with greater gene duplications. This study provides new insights into the structural variation and evolution of chloroplast tRNAs in cotton plants. Our findings could contribute to understanding the detailed characteristics and evolutionary variation of the tRNA family.

Molecular mechanisms involved in functional macroevolution of plant transcription factors

Transcription factors (TFs) are key components of the transcriptional regulation machinery. In plants, they accompanied the evolution from unicellular aquatic algae to complex flowering plants that dominate the land environment. The adaptations of the body plan and physiological responses required changes in the biological functions of TFs. Some ancestral gene regulatory networks are highly conserved, while others evolved more recently and only exist in particular lineages. The recent emergence of novel model organisms provided the opportunity for comparative studies, producing new insights to infer these evolutionary trajectories. In this review, we comprehensively revisit the recent literature on TFs of nonseed plants and algae, focusing on the molecular mechanisms driving their functional evolution. We discuss the particular contribution of changes in DNA-binding specificity, protein-protein interactions and cis-regulatory elements to gene regulatory networks. Current advances have shown that these evolutionary processes were shaped by changes in TF expression pattern, not through great innovation in TF protein sequences. We propose that the role of TFs associated with environmental and developmental regulation was unevenly conserved during land plant evolution.

Evolutionary Contribution of Duplicated Genes to Genome Evolution in the Ginseng Species Complex

Genes duplicated by whole genome duplication (WGD) and small-scale duplication (SSD) have played important roles in adaptive evolution of all flowering plants. However, it still remains underinvestigated how the distinct models of duplication events and their contending evolutionary patterns have shaped the genome and epigenomes of extant plant species. In this study, we investigated the contribution of the WGD- and SSD-derived duplicate genes to the genome evolution of one diploid and three closely related allotetraploid Panax species based on genome, methylome, and proteome data sets. Our genome-wide comparative analyses revealed that although the ginseng species complex was recently diverged, they have evolved distinct overall patterns of nucleotide variation, cytosine methylation, and protein-level expression. In particular, genetic and epigenetic asymmetries observed in the recent WGD-derived genes are largely consistent across the ginseng species complex. In addition, our results revealed that gene duplicates generated by ancient WGD and SSD mechanisms exhibited distinct evolutionary patterns. We found the ancient WGD-derived genes (i.e., ancient collinear gene) are genetically more conserved and hypomethylated at the cytosine sites. In contrast, some of the SSD-derived genes (i.e., dispersal duplicated gene) showed hypermethylation and high variance in nucleotide variation pattern. Functional enrichment analyses of the duplicated genes indicated that adaptation-related traits (i.e., photosynthesis) created during the distant ancient WGDs are further strengthened by both the more recent WGD and SSD. Together, our findings suggest that different types of duplicated genes may have played distinct but relaying evolutionary roles in the polyploidization and speciation processes in the ginseng species complex.

The eukaryotic MEP-pathway genes are evolutionarily conserved and originated from Chlaymidia and cyanobacteria

Background

Isoprenoids are the most ancient and essential class of metabolites produced in all organisms, either via mevalonate (MVA)-and/or methylerythritol phosphate (MEP)-pathways. The MEP-pathway is present in all plastid-bearing organisms and most eubacteria. However, no comprehensive study reveals the origination and evolutionary characteristics of MEP-pathway genes in eukaryotes.

Results

Here, detailed bioinformatics analyses of the MEP-pathway provide an in-depth understanding the evolutionary history of this indispensable biochemical route, and offer a basis for the co-existence of the cytosolic MVA- and plastidial MEP-pathway in plants given the established exchange of the end products between the two isoprenoid-biosynthesis pathways. Here, phylogenetic analyses establish the contributions of both cyanobacteria and Chlamydiae sequences to the plant’s MEP-pathway genes. Moreover, Phylogenetic and inter-species syntenic block analyses demonstrate that six of the seven MEP-pathway genes have predominantly remained as single-copy in land plants in spite of multiple whole-genome duplication events (WGDs). Substitution rate and domain studies display the evolutionary conservation of these genes, reinforced by their high expression levels. Distinct phenotypic variation among plants with reduced expression levels of individual MEP-pathway genes confirm the indispensable function of each nuclear-encoded plastid-targeted MEP-pathway enzyme in plant growth and development.

Conclusion

Collectively, these findings reveal the polyphyletic origin and restrict conservation of MEP-pathway genes, and reinforce the potential function of the individual enzymes beyond production of the isoprenoids intermediates.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07448-x.