Green Plant Lineage Research Articles

Order in Chaos: Lesser-Conserved and Repeat Structures in Dehydrins.

Dehydrins (Dhns) are a group of intrinsically disordered land plant proteins that are closely associated with tolerance of dehydrative stress. Dhns are recognized and classified by the presence and sequence of five different conserved segments, varying in length from 8 to 15 residues, separated by highly variable disordered regions. In addition to one or more copies of the diagnostic, fifteen-residue K segment, most Dhns can be classified into one of three major groups based on the mutually exclusive presence of three other conserved segments (H, Y, or F), with all three groups typically incorporating multi-serine S segments. Many Dhns also include repeat structures. From an input library of 8675 non-redundant candidate sequences, a specialized R script identified and classified 2658 complete and 236 partial Dhn sequences in all major green plant (Viridiplantae) lineages, including a few green algal genera. An examination of the connecting segments bridging the conserved segments identified additional conserved patterns, suggesting that multi-Y, S-K, and K-S domains may act as functional units. Dhn Decoder identified 857 Dhns with repeat structures, ranging from 3 short, simple repeats to elaborate variations with up to 45 repeats or repeats of up to 85 residues comprising 1 or more of the conserved segments, suggesting that internal sequence duplication is an important mode of evolution in Dhns.

Chromosome Numbers and Reproductive Life Cycles in Green Plants: A Phylotranscriptomic Perspective.

The strong correlation between reproductive life cycle type and chromosome numbers in green plants has been a long-standing mystery in evolutionary biology. Within green plants, the derived condition of heterosporous reproduction has emerged from the ancestral condition of homospory in disparate locations on the phylogenetic tree at least 11 times, of which three lineages are extant. In all green plant lineages where heterospory has emerged, there has been a significant downsizing in chromosome numbers. This dynamic has been investigated without clear answers for many decades. In this study, we combine known ideas from existing literature with novel methods, tools, and data to generate fresh insights into an old question. Using gene family evolution models and selection analyses, we identified gene families that have undergone significant expansion, contraction, or selection in heterosporous lineages. Alongside lineage-specific genomic changes, our results revealed shared genomic changes/trends among heterosporous lineages. We found expansions in gene families related to developmental regulation, signaling pathways, and stress responses across heterosporous groups. Notably, the MATE efflux family showed consistent expansion and evidence of selection in heterosporous lineages, suggesting a potentially conserved role in heterospory evolution. These findings could provide novel avenues to investigate and probe the underlying mechanism that may underpin the association between heterospory and genomic changes. The general importance of chromosome numbers, structure, and sizes in cellular biology notwithstanding, the association between the emergence of heterosporous reproduction and chromosome number reduction/genome downsizing is not fully understood. It remains unclear why there exists an association between aspects of biology at such disparate levels as reproductive life cycles and chromosome numbers/genome size. Exploring and answering this conundrum of evolutionary biology can add to our broader understanding of life sciences and of biology at different levels. Applying the novel tools and methods emerging from ongoing progress in biotechnology and computational sciences presents an opportunity to make new inroads into this long-standing question.

Genome-Wide Analysis of the Class III Peroxidase Gene Family in Physcomitrium patens and a Search for Clues to Ancient Class III Peroxidase Functions

Plant class III peroxidases (PRXs) catalyze generation of reactive oxygen species and oxidation of various compounds including lignin precursors. PRXs function in cell wall metabolism, defense, and stress responses. However, gene redundancy and catalytic versatility have impeded detailed functional characterization of PRX genes. The genome of the model moss Physcomitrium patens harbors a relatively small number (49) of PRX genes. Conserved architecture of four exons and three ‘001’ introns, found in some algal PRX genes and in the PpPRX family, suggests that this architecture predated divergence of the green algal and land plant lineages. The PpPRX family expanded mainly through whole-genome duplications. All duplicated pairs but one were under purifying selection and generally exhibited similar expression profiles. An expanded phylogenetic tree revealed a conserved land plant-wide clade that contained PRXs implicated in stress responses in non-lignifying cells, providing a clue to ancient functions of land plant PRXs. Functional clustering was not observed, suggesting convergent evolution of specific PRX functions (e.g., lignification) in different plant lineages. With its small complement of PRXs, P. patens may be useful for functional characterization of land plant PRXs. Several PpPRXs were proposed for further study, including PpPRX34 and PpPRX39 in the ancient land plant-wide clade.

Evolutionary Analysis and Catalytic Function of LOG Proteins in Plants.

The plant hormone cytokinin is a conserved regulator of plant development. LONELY GUY (LOG) proteins are pivotal in cytokinin biosynthesis. However, their origin, evolutionary history, and enzymatic characteristics remain largely uncharacterized. To elucidate LOG family evolution history and protein motif composition, we conducted phylogenetic and motif analyses encompassing representative species across the whole green plant lineage. Catalytic activity and structure analysis were conducted to thoroughly characterize the LOG proteins. Our phylogeny showed that LOG proteins could be divided into five groups and revealed three major duplication events giving rise to four distinct groups of vascular LOG proteins. LOG proteins share a conserved structure characterized by a canonical motif arrangement comprising motifs 1, 2, 3, 4, 5, 6, and 7. Two significant changes in LOG motif composition occurred during the transition to land plants: the emergence of motif 3 in charophyte LOG sequences and the subsequent acquisition of motif 8 at the C-terminus of LOG proteins. Enzymatic assays demonstrated that LOG proteins can be classified into two groups based on their enzyme activity. One group act as cytokinin riboside 5'-monophosphate phosphoribohydrolase and primarily convert iPRMP to iP, while the other group act as 5'-ribonucleotide phosphohydrolase, and preferentially produce iPR from the same substrates. TaLOG5-4A1, TaLOG5-4A2, TaLOG5-5B2, and TaLOG5-D1 shared conserved residues in the critical motif and were predicted to have similar protein structures, but displayed distinct enzymatic activities. Our findings provide a comprehensive overview of LOG protein phylogeny and lay a foundation for further investigations into their functional diversification.

Genome-wide patterns of selection-drift variation strongly associate with organismal traits across the green plant lineage.

There are many gaps in our knowledge of how life cycle variation and organismal body architecture associate with molecular evolution. Using the diverse range of green algal body architectures and life cycle types as a test case, we hypothesize that increases in cytomorphological complexity are likely to be associated with a decrease in the effective population size, because larger-bodied organisms typically have smaller populations, resulting in increased drift. For life cycles, we expect haploid-dominant lineages to evolve under stronger selection intensity relative to diploid-dominant life cycles owing to masking of deleterious alleles in heterozygotes. We use a genome-scale data set spanning the phylogenetic diversity of green algae and phylogenetic comparative approaches to measure the relative selection intensity across different trait categories. We show stronger signatures of drift in lineages with more complex body architectures compared with unicellular lineages, which we consider to be a consequence of smaller effective population sizes of the more complex algae. Significantly higher rates of synonymous as well as nonsynonymous substitutions relative to other algal body architectures highlight that siphonous and siphonocladous body architectures, characteristic of many green seaweeds, form an interesting test case to study the potential impacts of genome redundancy on molecular evolution. Contrary to expectations, we show that levels of selection efficacy do not show a strong association with life cycle types in green algae. Taken together, our results underline the prominent impact of body architecture on the molecular evolution of green algal genomes.

Diversification of plant SUPPRESSOR OF MAX2 1 (SMAX1)-like genes and genome-wide identification and characterization of cotton SMXL gene family

BackgroundStrigolactones (SLs) are a recently discovered class of plant hormones. SUPPRESSOR OF MAX2 1 (SMAX1)-like proteins, key component of the SL signaling pathway, have been studied extensively for their roles in regulating plant growth and development, such as plant branching. However, systematic identification and functional characterization of SMXL genes in cotton (Gossypium sp.), an important fiber and oil crop, has rarely been conducted.ResultsWe identified 210 SMXL genes from 21 plant genomes and examined their evolutionary relationships. The structural characteristics of the SMXL genes and their encoded proteins exhibited both consistency and diversity. All plant SMXL proteins possess a conserved Clp-N domain, P-loop NTPase, and EAR motif. We identified 63 SMXL genes in cotton and classified these into four evolutionary branches. Gene expression analysis revealed tissue-specific expression patterns of GhSMXL genes, with some upregulated in response to GR24 treatment. Protein co-expression network analysis showed that GhSMXL6, GhSMXL7-1, and GhSMXL7-2 mainly interact with proteins functioning in growth and development, while virus-induced gene silencing revealed that GhSMAX1-1 and GhSMAX1-2 suppress the growth and development of axillary buds.ConclusionsSMXL gene family members show evolutionary diversification through the green plant lineage. GhSMXL6/7–1/7–2 genes play critical roles in the SL signaling pathway, while GhSMXL1-1 and GhSMXL1-2 function redundantly in growth of axillary buds. Characterization of the cotton SMXL gene family provides new insights into their roles in responding to SL signals and in plant growth and development. Genes identified in this study could be used as the candidate genes for improvement of plant architecture and crop yield.

Bioinformatics Analysis of MSH1 Genes of Green Plants: Multiple Parallel Length Expansions, Intron Gains and Losses, Partial Gene Duplications, and Alternative Splicing.

MutS homolog 1 (MSH1) is involved in the recombining and repairing of organelle genomes and is essential for maintaining their stability. Previous studies indicated that the length of the gene varied greatly among species and detected species-specific partial gene duplications in Physcomitrella patens. However, there are critical gaps in the understanding of the gene size expansion, and the extent of the partial gene duplication of MSH1 remains unclear. Here, we screened MSH1 genes in 85 selected species with genome sequences representing the main clades of green plants (Viridiplantae). We identified the MSH1 gene in all lineages of green plants, except for nine incomplete species, for bioinformatics analysis. The gene is a singleton gene in most of the selected species with conserved amino acids and protein domains. Gene length varies greatly among the species, ranging from 3234 bp in Ostreococcus tauri to 805,861 bp in Cycas panzhihuaensis. The expansion of MSH1 repeatedly occurred in multiple clades, especially in Gymnosperms, Orchidaceae, and Chloranthus spicatus. MSH1 has exceptionally long introns in certain species due to the gene length expansion, and the longest intron even reaches 101,025 bp. And the gene length is positively correlated with the proportion of the transposable elements (TEs) in the introns. In addition, gene structure analysis indicated that the MSH1 of green plants had undergone parallel intron gains and losses in all major lineages. However, the intron number of seed plants (gymnosperm and angiosperm) is relatively stable. All the selected gymnosperms contain 22 introns except for Gnetum montanum and Welwitschia mirabilis, while all the selected angiosperm species preserve 21 introns except for the ANA grade. Notably, the coding region of MSH1 in algae presents an exceptionally high GC content (47.7% to 75.5%). Moreover, over one-third of the selected species contain species-specific partial gene duplications of MSH1, except for the conserved mosses-specific partial gene duplication. Additionally, we found conserved alternatively spliced MSH1 transcripts in five species. The study of MSH1 sheds light on the evolution of the long genes of green plants.

Occurrence and conversion of progestogens and androgens are conserved in land plants.

Progestogens and androgens have been found in many plants, but little is known about their biosynthesis and the evolution of steroidogenesis in these organisms. Here, we show that the occurrence and biosynthesis of progestogens and androgens are conserved across the viridiplantae lineage. An UHPLC-ESI-MS/MS method allowed high-throughput analysis of the occurrence and chemical conversion of progestogens and androgens in 41 species across the green plant lineage. Dehydroepiandrosterone, testosterone, and 5α-dihydrotestosterone are plants' most abundant mammalian-like steroids. Progestogens are converted into 17α-hydroxyprogesterone and 5α-pregnane-3,20-dione. Androgens are converted into testosterone and 5α-dihydrotestosterone. 17,20-Lyases, essential for converting progestogens to androgens, seem to be most effective in monocot species. Our data suggest that the occurrence of progestogens and androgens is highly conserved in plants, and their biosynthesis might favor a route using the Δ4 pathway.

Genome-Wide Identification and Expression Analysis of Calmodulin (CaM) and Calmodulin-Like (CML) Genes in the Brown Algae Saccharina japonica.

Calmodulins (CaMs) and Calmodulin-like proteins (CMLs) are vital in plant growth, development, and stress responses. However, CaMs and CMLs have not been fully identified and characterized in brown algae, which has been evolving independently of the well-studied green plant lineage. In this study, whole-genome searches revealed one SjCaM and eight SjCMLs in Saccharina japonica, and one EsCaM and eleven EsCMLs in Ectocarpus sp. SjCaM and EsCaM encoded identical protein products and shared 88.59-89.93% amino acid identities with Arabidopsis thaliana AtCaMs, thereby indicating that brown algae CaMs retained a similar Ca2+ sensors function as in plants. The phylogenetic and gene structure analysis results showed that there was significant divergence in the gene sequences among brown algae CMLs. Furthermore, evolutionary analysis indicated that the function of brown alga CMLs was relatively conserved, which may be related to the fact that brown algae do not need to face complex environments like terrestrial plants. Regulatory elements prediction and the expression analysis revealed the probable functioning of SjCaM/CML genes in gametophyte development and the stress response in S. japonica. In addition, the SjCaM/SjCMLs interacting proteins and chemicals were preliminarily predicted, suggesting that SjCaM/SjCMLs might play putative roles in Ca2+/CaM-mediated growth and development processes and stimulus responses. Therefore, these results will facilitate our understanding of the evolution of brown algae CaMs/CMLs and the functional identification of SjCaM/SjCMLs.

Chromosome-level reference genome of tetraploid Isoetes sinensis provides insights into evolution and adaption of lycophytes.

The Lycophyta species are the extant taxa most similar to early vascular plants that were once abundant on Earth. However, their distribution has greatly diminished. So far, the absence of chromosome-level assembled lycophyte genomes has hindered our understanding of evolution and environmental adaption of lycophytes. We present the reference genome of the tetraploid aquatic quillwort, Isoetes sinensis, a lycophyte. This genome represents the first chromosome-level assembled genome of a tetraploid seed-free plant. Comparison of genomes between I. sinensis and Isoetestaiwanensis revealed conserved and different genomic features between diploid and polyploid lycophytes. Comparison of the I. sinensis genome with those of other species representing the evolutionary lineages of green plants revealed the inherited genetic tools for transcriptional regulation and most phytohormones in I. sinensis. The presence and absence of key genes related to development and stress responses provide insights into environmental adaption of lycophytes. The high-quality reference genome and genomic analysis presented in this study are crucial for future genetic and environmental studies of not only I. sinensis but also other lycophytes.

The origin and evolution of salicylic acid signaling and biosynthesis in plants

Salicylic acid (SA) plays a pivotal role in plant response to biotic and abiotic stress. Several core SA signaling regulators and key proteins in SA biosynthesis have been well characterized. However, much remains unknown about the origin, evolution, and early diversification of core elements in plant SA signaling and biosynthesis. In this study, we identified 10 core protein families in SA signaling and biosynthesis across green plant lineages. We found that the key SA signaling receptors, the nonexpresser of pathogenesis-related (NPR) proteins, originated in the most recent common ancestor (MRCA) of land plants and formed divergent groups in the ancestor of seed plants. However, key transcription factors for SA signaling, TGACG motif-binding proteins (TGAs), originated in the MRCA of streptophytes, arguing for the stepwise evolution of core SA signaling in plants. Different from the assembly of the core SA signaling pathway in the ancestor of seed plants, SA exists extensively in green plants, including chlorophytes and streptophyte algae. However, the full isochorismate synthase (ICS)-based SA synthesis pathway was first assembled in the MRCA of land plants. We further revealed that the ancient abnormal inflorescence meristem 1 (AIM1)-based β-oxidation pathway is crucial for the biosynthesis of SA in chlorophyte algae, and this biosynthesis pathway may have facilitated the adaptation of early-diverging green algae to the high-light-intensity environment on land. Taken together, our findings provide significant insights into the early evolution and diversification of plant SA signaling and biosynthesis pathways, highlighting a crucial role of SA in stress tolerance during plant terrestrialization.

Comprehensive evolutionary analysis of the TCP gene family: Further insights for its origin, expansion, and diversification.

TCP proteins are plant-specific transcription factors, which are involved in a broad range of physiological processes of plant growth and development. However, the origin and evolutionary history of this gene family is not fully resolved. Here, we present a genome-wide survey of TCP genes in 59 species (including 42 genomes and 17 transcriptomes) covering all main lineages of green plants, and reconstruct the evolutionary history of this gene family. Our results suggested that the origin of TCP genes predated the emergence of land plants, possibly in the common ancestor of Phragmoplastophyta. The TCP gene family gradually experienced a continuous expansion and grew from a few members in algae, moss and lycophytes to dozens, and sometimes over 50 members in angiosperms. Phylogenetic analysis indicated that at least four subclades (Class I and three subclades of Class II) have been occurred in the ancestor of spermatophyte (seed plant). Both dispersed duplication and segmental duplication or whole-genome duplication (WGD) contributed significantly to the expansion of the TCP gene family over the course of evolution. Our findings provide a comprehensive evolutionary analysis of the TCP gene family and highlight the importance of gene duplications in the evolution of this plant-specific transcription factors.

Molecular Evolution of Calcium Signaling and Transport in Plant Adaptation to Abiotic Stress.

Adaptation to unfavorable abiotic stresses is one of the key processes in the evolution of plants. Calcium (Ca2+) signaling is characterized by the spatiotemporal pattern of Ca2+ distribution and the activities of multi-domain proteins in integrating environmental stimuli and cellular responses, which are crucial early events in abiotic stress responses in plants. However, a comprehensive summary and explanation for evolutionary and functional synergies in Ca2+ signaling remains elusive in green plants. We review mechanisms of Ca2+ membrane transporters and intracellular Ca2+ sensors with evolutionary imprinting and structural clues. These may provide molecular and bioinformatics insights for the functional analysis of some non-model species in the evolutionarily important green plant lineages. We summarize the chronological order, spatial location, and characteristics of Ca2+ functional proteins. Furthermore, we highlight the integral functions of calcium-signaling components in various nodes of the Ca2+ signaling pathway through conserved or variant evolutionary processes. These ultimately bridge the Ca2+ cascade reactions into regulatory networks, particularly in the hormonal signaling pathways. In summary, this review provides new perspectives towards a better understanding of the evolution, interaction and integration of Ca2+ signaling components in green plants, which is likely to benefit future research in agriculture, evolutionary biology, ecology and the environment.

PLAZA 5.0: extending the scope and power of comparative and functional genomics in plants.

PLAZA is a platform for comparative, evolutionary, and functional plant genomics. It makes a broad set of genomes, data types and analysis tools available to researchers through a user-friendly website, an API, and bulk downloads. In this latest release of the PLAZA platform, we are integrating a record number of 134 high-quality plant genomes, split up over two instances: PLAZA Dicots 5.0 and PLAZA Monocots 5.0. This number of genomes corresponds with a massive expansion in the number of available species when compared to PLAZA 4.0, which offered access to 71 species, a 89% overall increase. The PLAZA 5.0 release contains information for 5 882 730 genes, and offers pre-computed gene families and phylogenetic trees for 5 274 684 protein-coding genes. This latest release also comes with a set of new and updated features: a new BED import functionality for the workbench, improved interactive visualizations for functional enrichments and genome-wide mapping of gene sets, and a fully redesigned and extended API. Taken together, this new version offers extended support for plant biologists working on different families within the green plant lineage and provides an efficient and versatile toolbox for plant genomics. All PLAZA releases are accessible from the portal website: https://bioinformatics.psb.ugent.be/plaza/.

Molecular Interaction and Evolution of Jasmonate Signaling With Transport and Detoxification of Heavy Metals and Metalloids in Plants.

An increase in environmental pollution resulting from toxic heavy metals and metalloids [e.g., cadmium (Cd), arsenic (As), and lead (Pb)] causes serious health risks to humans and animals. Mitigation strategies need to be developed to reduce the accumulation of the toxic elements in plant-derived foods. Natural and genetically-engineered plants with hyper-tolerant and hyper-accumulating capacity of toxic minerals are valuable for phytoremediation. However, the molecular mechanisms of detoxification and accumulation in plants have only been demonstrated in very few plant species such as Arabidopsis and rice. Here, we review the physiological and molecular aspects of jasmonic acid and the jasmonate derivatives (JAs) in response to toxic heavy metals and metalloids. Jasmonates have been identified in, limiting the accumulation and enhancing the tolerance to the toxic elements, by coordinating the ion transport system, the activity of antioxidant enzymes, and the chelating capacity in plants. We also propose the potential involvement of Ca2+ signaling in the stress-induced production of jasmonates. Comparative transcriptomics analyses using the public datasets reveal the key gene families involved in the JA-responsive routes. Furthermore, we show that JAs may function as a fundamental phytohormone that protects plants from heavy metals and metalloids as demonstrated by the evolutionary conservation and diversity of these gene families in a large number of species of the major green plant lineages. Using ATP-Binding Cassette G (ABCG) transporter subfamily of six representative green plant species, we propose that JA transporters in Subgroup 4 of ABCGs may also have roles in heavy metal detoxification. Our paper may provide guidance toward the selection and development of suitable plant and crop species that are tolerant to toxic heavy metals and metalloids.

Arabidopsis Mitochondrial Transcription Termination Factor mTERF2 Promotes Splicing of Group IIB Introns.

Plastid gene expression (PGE) is essential for chloroplast biogenesis and function and, hence, for plant development. However, many aspects of PGE remain obscure due to the complexity of the process. A hallmark of nuclear-organellar coordination of gene expression is the emergence of nucleus-encoded protein families, including nucleic-acid binding proteins, during the evolution of the green plant lineage. One of these is the mitochondrial transcription termination factor (mTERF) family, the members of which regulate various steps in gene expression in chloroplasts and/or mitochondria. Here, we describe the molecular function of the chloroplast-localized mTERF2 in Arabidopsis thaliana. The complete loss of mTERF2 function results in embryo lethality, whereas directed, microRNA (amiR)-mediated knockdown of MTERF2 is associated with perturbed plant development and reduced chlorophyll content. Moreover, photosynthesis is impaired in amiR-mterf2 plants, as indicated by reduced levels of photosystem subunits, although the levels of the corresponding messenger RNAs are not affected. RNA immunoprecipitation followed by RNA sequencing (RIP-Seq) experiments, combined with whole-genome RNA-Seq, RNA gel-blot, and quantitative RT-PCR analyses, revealed that mTERF2 is required for the splicing of the group IIB introns of ycf3 (intron 1) and rps12.

Transcriptional profiling reveals conserved and species-specific plant defense responses during the interaction of Physcomitrium patens with Botrytis cinerea.

Evolutionary conserved defense mechanisms present in extant bryophytes and angiosperms, as well as moss-specific defenses are part of the immune response of Physcomitrium patens. Bryophytes and tracheophytes are descendants of early land plants that evolved adaptation mechanisms to cope with different kinds of terrestrial stresses, including drought, variations in temperature and UV radiation, as well as defense mechanisms against microorganisms present in the air and soil. Although great advances have been made on pathogen perception and subsequent defense activation in angiosperms, limited information is available in bryophytes. In this study, a transcriptomic approach uncovered the molecular mechanisms underlying the defense response of the bryophyte Physcomitrium patens (previously Physcomitrella patens) against the important plant pathogen Botrytis cinerea. A total of 3.072 differentially expressed genes were significantly affected during B. cinerea infection, including genes encoding proteins with known functions in angiosperm immunity and involved in pathogen perception, signaling, transcription, hormonal signaling, metabolic pathways such as shikimate and phenylpropanoid, and proteins with diverse role in defense against biotic stress. Similarly as in other plants, B. cinerea infection leads to downregulation of genes involved in photosynthesis and cell cycle progression. These results highlight the existence of evolutionary conserved defense responses to pathogens throughout the green plant lineage, suggesting that they were probably present in the common ancestors of land plants. Moreover, several genes acquired by horizontal transfer from prokaryotes and fungi, and a high number of P. patens-specific orphan genes were differentially expressed during B. cinerea infection, suggesting that they are important players in the moss immune response.

The Transfer of the Ferredoxin Gene From the Chloroplast to the Nuclear Genome Is Ancient Within the Paraphyletic Genus Thalassiosira.

Ferredoxins are iron–sulfur proteins essential for a wide range of organisms because they are an electron transfer mediator involved in multiple metabolic pathways. In phytoplankton, these proteins are active in the mature chloroplasts, but the petF gene, encoding for ferredoxin, has been found either to be in the chloroplast genome or transferred to the nuclear genome as observed in the green algae and higher plant lineage. We experimentally determined the location of the petF gene in 12 strains of Thalassiosira covering three species using DNA sequencing and qPCR assays. The results showed that petF gene is located in the nuclear genome of all confirmed Thalassiosira oceanica strains (CCMP0999, 1001, 1005, and 1006) tested. In contrast, all Thalassiosira pseudonana (CCMP1012, 1013, 1014, and 1335) and Thalassiosira weissflogii (CCMP1010, 1049, and 1052) strains studied retained the gene in the chloroplast genome, as generally observed for Bacillariophyceae. Our evolutionary analyses further extend the dataset on the localization of the petF gene in the Thalassiosirales. The realization that the petF gene is nuclear-encoded in the Skeletonema genus allowed us to trace the petF gene transfer back to a single event that occurred within the paraphyletic genus Thalassiosira. Phylogenetic analyses revealed the need to reassess the taxonomic assignment of the Thalassiosira strain CCMP1616, since the genes used in our study did not cluster within the T. oceanica lineage. Our results suggest that this strains’ diversification occurred prior to the ferredoxin gene transfer event. The functional transfer of petF genes provides insight into the evolutionary processes leading to chloroplast genome reduction and suggests ecological adaptation as a driving force for such chloroplast to nuclear gene transfer.

Independent recurrent evolution of MICRORNA genes converging onto similar non-canonical organisation across green plant lineages is driven by local and segmental duplication events in species, family and lineages

The relationship between evolutionary history, organisation and transcriptional regulation of genes are intrinsically linked. These have been well studied in canonically organised protein-coding genes but not of MIRNAs. In the present study, we investigated the non-canonical arrangement of MIRNAs across taxonomic boundaries from algae to angiosperms employing a combination of genome organization, phylogeny and synteny. We retrieved the complete dataset of MIRNA from twenty-five species to identify and classify based on organisational patterns. The median size of cluster was between 2-5 kb and between 1-20 % of all MIRNAs are organized in head-to-head (with bidirectional promoter), head-to-tail (tandem), and overlapping manner. Although majority of the clusters are composed of MIRNA homologs, 25% of all clusters comprises of non-homologous genes with a potential of generating functional and regulatory complexity. A comparison of phylogeny and organizational patterns revealed that multiple independent events, some of which are species-specific, and some ancient, in different lineages, are responsible for non-canonical organization. Detailed investigation of MIR395 family across the plants revealed a complex origin of non-canonical arrangement through ancient and recent, segmental and local duplications; analysis of MIR399 family revealed major expansion occurred prior to monocot-dicot split, with few lineage-specific events. Evolution of “convergent” organization pattern of non-canonical arrangement originating from independent loci through recurrent event highlights our poor understanding of evolutionary process of MIRNA genes. The present investigation thus paves way for comparative functional genomics to understand the role of non-canonical organization on transcriptional regulation and regulatory diversity in MIRNA gene families.

Phylogenomics reveals convergent evolution of red-violet coloration in land plants and the origins of the anthocyanin biosynthetic pathway

The flavonoids, one of the largest classes of plant secondary metabolites, are found in lineages that span the land plant phylogeny and play important roles in stress responses and as pigments. Perhaps the most well-studied flavonoids are the anthocyanins that have human health benefits and help plants attract pollinators, regulate hormone production, and confer resistance to abiotic and biotic stresses. The canonical biochemical pathway responsible for the production of these pigments is well-characterized for flowering plants yet its conservation across deep divergences in land plants remains debated and poorly understood. Many early land plants such as mosses, liverworts, and ferns produce flavonoid pigments, but their biosynthetic origins and homologies to the anthocyanin pathway remain uncertain. We conducted phylogenetic analyses using full genome sequences representing nearly all major green plant lineages to reconstruct the evolutionary history of the anthocyanin biosynthetic pathway then test the hypothesis that genes in this pathway are present in early land plants. We found that the entire pathway was not intact until the most recent common ancestor of seed plants and that orthologs of many downstream enzymes are absent from seedless plants including mosses, liverworts, and ferns. Our results also highlight the utility of phylogenetic inference, as compared to pairwise sequence similarity, in orthology assessment within large gene families that have complex duplication-loss histories. We suggest that the production of red-violet flavonoid pigments widespread in seedless plants, including the 3-deoxyanthocyanins, requires the activity of novel, as-yet discovered enzymes, and represents convergent evolution of red-violet coloration across land plants.