Phylogenomics and the rise of the angiosperms.

Angiosperms are the cornerstone of most terrestrial ecosystems and human livelihoods1,2. A robust understanding of angiosperm evolution is required to explain their rise to ecological dominance. So far, the angiosperm tree of life has been determined primarily by means of analyses of the plastid genome3,4. Many studies have drawn on this foundational work, such as classification and first insights into angiosperm diversification since their Mesozoic origins5-7. However, the limited and biased sampling of both taxa and genomes undermines confidence in the tree and its implications. Here, we build the tree of life for almost 8,000 (about 60%) angiosperm genera using a standardized set of 353 nuclear genes8. This 15-fold increase in genus-level sampling relative to comparable nuclear studies9 provides a critical test of earlier results and brings notable change to key groups, especially in rosids, while substantiating many previously predicted relationships. Scaling this tree to time using 200 fossils, we discovered that early angiosperm evolution was characterized by high gene tree conflict and explosive diversification, giving rise to more than 80% of extant angiosperm orders. Steady diversification ensued through the remaining Mesozoic Era until rates resurged in the Cenozoic Era, concurrent with decreasing global temperatures and tightly linked with gene tree conflict. Taken together, our extensive sampling combined with advanced phylogenomic methods shows the deep history and full complexity in the evolution of a megadiverse clade.

Pollinators contribute to the maintenance of flowering plant diversity.

Mechanisms that favour rare species are key to the maintenance of diverse communities1-3. One of the most critical tasks for conservation of flowering plant biodiversity is to understand how plant-pollinator interactions contribute to the maintenance of rare species4-7. Here we show that niche partitioning in pollinator use and asymmetric facilitation confer fitness advantage of rarer species in a biodiversity hotspot using phylogenetic structural equation modelling that integrates plant-pollinator and interspecific pollen transfer networks with floral functional traits. Co-flowering species filtered pollinators via floral traits, and rarer species showed greater pollinator specialization leading to higher pollination-mediated male and female fitness than more abundant species. When plants shared pollinator resources, asymmetric facilitation via pollen transport dynamics benefitted the rarer species at the cost of more abundant species, serving as an alternative diversity-promoting mechanism. Our results emphasize the importance of community-wide plant-pollinator interactions that affect reproduction for biodiversity maintenance.

Inferring the Evolutionary Model of Community-Structuring Traits with Convolutional Kitchen Sinks.

When communities are assembled through processes such as filtering or limiting similarity acting on phylogenetically conserved traits, the evolutionary signature of those traits may be reflected in patterns of community membership. We show how the model of trait evolution underlying community-structuring traits can be inferred from community membership data using both a variation of a traditional eco-phylogenetic metric-the mean pairwise phylogenetic distance (MPD) between taxa-and a recent machine learning tool, Convolutional Kitchen Sinks (CKS). Both methods perform well across a range of phylogenetically informative evolutionary models, but CKS outperforms MPD as tree size increases. We demonstrate CKS by inferring the evolutionary history of freeze tolerance in angiosperms. Our analysis is consistent with a late burst model, suggesting freeze tolerance evolved recently. We suggest that multiple data types that are ordered on phylogenies, such as trait values, species interactions, or community presence/absence, are good candidates for CKS modeling because the generative models produce structured differences between neighboring points that CKS is well-suited for. We introduce the R package kitchen to perform CKS for generic application of the technique.

ANAgdb: a multi-omics and taxonomy database for ANA-grade.

BACKGROUND: The ANA-grade, encompassing early-diverging angiosperm lineages, Amborellales, Nymphaeales, and Austrobaileyales, represents a fundamental phase in the evolutionary history of flowering plants. Since the completion of key assembly of the Amborella genome, the continuous influx of omics data from the lineage underscores the need for a specialized database. RESULTS: Here, we introduce the ANA-grade Genome Database (ANAgdb, https://anagenome.cn/ ), which integrates multi-omics data including 11 genomes, 167 transcriptomes, and 10 miRNAomes, as well as extensive taxonomic details specific to the ANA-grade. Designed with an array of user-friendly tools, ANAgdb not only facilitates the effective storage, querying, and analysis of data but also enables the integration and dissemination of crucial genomic and taxonomic information. CONCLUSION: By integrating the comprehensive resources and tools, ANAgdb aims to significantly advance research in phylogenomics and taxonomic studies, providing a robust platform for researchers to explore the genetic and morphological diversities of these ancient plant lineages.

Cretaceous and Paleocene fossils reveal an extinct higher clade within Cornales, the dogwood order.

PREMISE: Characterization and phylogenetic integration of fossil angiosperms with uncertain affinities is relatively limited, which may obscure the diversity of extinct higher taxa in the flowering plant tree of life. The order Cornales contains a diversity of extinct taxa with uncertain familial affinities that make it an ideal group for studying turnover in angiosperms. Here, we describe a new extinct genus of Cornales unassignable to an extant family and conduct a series of phylogenetic analyses to reconstruct relationships of fossils across the order. METHODS: Two permineralized endocarps were collected from the Cedar District Formation (Campanian, 82-80 Ma) of Sucia Island, State of Washington, United States. Fossils were sectioned with the cellulose acetate peel technique and incorporated into a morphological dataset. To assess the utility of this dataset to accurately place taxa in their respective clades, we used a series of phylogenetic pseudofossilization analyses. We then conducted a total-evidence analysis and a scaffold-based approach to determine relationships of fossils. RESULTS: Based on their unique combination of characters, the fossils represent a new genus, Fenestracarpa washingtonensis gen. nov. et sp. nov. Pseudofossilization analyses indicate that our morphological dataset can be used to accurately recover taxa at the major clade to family level, generally with moderate to high support. The total-evidence and scaffold-based analyses recovered Fenestracarpa and other fossil genera in an entirely extinct clade within Cornales. CONCLUSIONS: Our findings increase the reported diversity of extinct Cornales and indicate that the order's initial radiation likely included the divergence of an extinct higher clade that endured the end-Cretaceous Mass extinction but perished during the Cenozoic.

Evolutionary history of the angiosperm flora of China.

High species diversity may result from recent rapid speciation in a 'cradle' and/or the gradual accumulation and preservation of species over time in a 'museum'. China harbours nearly 10% of angiosperm species worldwide and has long been considered as both a museum, owing to the presence of many species with hypothesized ancient origins, and a cradle, as many lineages have originated as recent topographic changes and climatic shifts-such as the formation of the Qinghai-Tibetan Plateau and the development of the monsoon-provided new habitats that promoted remarkable radiation. However, no detailed phylogenetic study has addressed when and how the major components of the Chinese angiosperm flora assembled to form the present-day vegetation. Here we investigate the spatio-temporal divergence patterns of the Chinese flora using a dated phylogeny of 92% of the angiosperm genera for the region, a nearly complete species-level tree comprising 26,978 species and detailed spatial distribution data. We found that 66% of the angiosperm genera in China did not originate until early in the Miocene epoch (23 million years ago (Mya)). The flora of eastern China bears a signature of older divergence (mean divergence times of 22.04-25.39 Mya), phylogenetic overdispersion (spatial co-occurrence of distant relatives) and higher phylogenetic diversity. In western China, the flora shows more recent divergence (mean divergence times of 15.29-18.86 Mya), pronounced phylogenetic clustering (co-occurrence of close relatives) and lower phylogenetic diversity. Analyses of species-level phylogenetic diversity using simulated branch lengths yielded results similar to genus-level patterns. Our analyses indicate that eastern China represents a floristic museum, and western China an evolutionary cradle, for herbaceous genera; eastern China has served as both a museum and a cradle for woody genera. These results identify areas of high species richness and phylogenetic diversity, and provide a foundation on which to build conservation efforts in China.

A Collaborative Classroom Investigation of the Evolution of SABATH Methyltransferase Substrate Preference Shifts over 120 My of Flowering Plant History.

Next-generation sequencing has resulted in an explosion of available data, much of which remains unstudied in terms of biochemical function; yet, experimental characterization of these sequences has the potential to provide unprecedented insight into the evolution of enzyme activity. One way to make inroads into the experimental study of the voluminous data available is to engage students by integrating teaching and research in a college classroom such that eventually hundreds or thousands of enzymes may be characterized. In this study, we capitalize on this potential to focus on SABATH methyltransferase enzymes that have been shown to methylate the important plant hormone, salicylic acid (SA), to form methyl salicylate. We analyze data from 76 enzymes of flowering plant species in 23 orders and 41 families to investigate how widely conserved substrate preference is for SA methyltransferase orthologs. We find a high degree of conservation of substrate preference for SA over the structurally similar metabolite, benzoic acid, with recent switches that appear to be associated with gene duplication and at least three cases of functional compensation by paralogous enzymes. The presence of Met in active site position 150 is a useful predictor of SA methylation preference in SABATH methyltransferases but enzymes with other residues in the homologous position show the same substrate preference. Although our dense and systematic sampling of SABATH enzymes across angiosperms has revealed novel insights, this is merely the "tip of the iceberg" since thousands of sequences remain uncharacterized in this enzyme family alone.

Network-based microsynteny analysis identifies major differences and genomic outliers in mammalian and angiosperm genomes.

A comprehensive analysis of relative gene order, or microsynteny, can provide valuable information for understanding the evolutionary history of genes and genomes, and ultimately traits and species, across broad phylogenetic groups and divergence times. We have used our network-based phylogenomic synteny analysis pipeline to first analyze the overall patterns and major differences between 87 mammalian and 107 angiosperm genomes. These two important groups have both evolved and radiated over the last ∼170 MYR. Secondly, we identified the genomic outliers or "rebel genes" within each clade. We theorize that rebel genes potentially have influenced trait and lineage evolution. Microsynteny networks use genes as nodes and syntenic relationships between genes as edges. Networks were decomposed into clusters using the Infomap algorithm, followed by phylogenomic copy-number profiling of each cluster. The differences in syntenic properties of all annotated gene families, including BUSCO genes, between the two clades are striking: most genes are single copy and syntenic across mammalian genomes, whereas most genes are multicopy and/or have lineage-specific distributions for angiosperms. We propose microsynteny scores as an alternative and complementary metric to BUSCO for assessing genome assemblies. We further found that the rebel genes are different between the two groups: lineage-specific gene transpositions are unusual in mammals, whereas single-copy highly syntenic genes are rare for flowering plants. We illustrate several examples of mammalian transpositions, such as brain-development genes in primates, and syntenic conservation across angiosperms, such as single-copy genes related to photosynthesis. Future experimental work can test if these are indeed rebels with a cause.

Comprehensive analysis of polygalacturonase genes offers new insights into their origin and functional evolution in land plants.

Polygalacturonase (PG) is a hydrolase that participates in pectin degradation, pod shattering and fruit softening. Here, we identified 2786 PG genes across 54 plants, which could be divided into three groups. Evolutionary analysis suggested that PG family originated from the charophyte green algae, and Subgroups A2-A4 evolved from the Subgroup A1 after the tracheophyte-angiosperm split. Whole-genome duplication was the major force leading to PG gene expansion. Interestingly, the PG genes continuously expanded in eudicots, whereas it contracted in monocots after the eudicot-monocot split. PG genes in Group A are expressed at high levels in floral organs, whereas genes in Groups B and C are expressed at high levels in various tissues. Moreover, three BnaPG15 members were found for their potential possibility in pod shattering in Brassica napus. Our results provide new insight into the evolutionary history of PG family, and their potentially functional role in plants.

New insights on the function of plant acyl carrier proteins from comparative and evolutionary analysis.

Acyl carrier proteins (ACPs) play a central role in both plastidial and mitochondrial Type II fatty acid synthesis in plant cells. However, a large proportion of plant ACPs remain functionally uncharacterized, and their evolutionary history remains elusive. In present study, 97 putative ACPs were identified from ten angiosperm species examined. Based on phylogenetic analysis, ACP genes were grouped into plastidial (cpACP: ACP1/2/3/4/5) and mitochondrial (mtACP: mtACP1/mtACP2/mtACP3) ACPs. Protein sequence (motifs and length), tertiary structure, and gene structure (exon number, average intron length, and intron phase) were highly conserved in different ACP subclades. The differentiation of ACPs into distinct types occurred 85-98 and 45-57 million years ago. A limited proportion of ACP genes experience tandem or segmental duplication, corresponding to two rounds of whole genome duplication. Ka/Ks ratios revealed that duplicated ACP genes underwent a purifying selection. Regarding expression patterns, most ACPs were expressed constitutively and tissue-specifically. Notably, the average expression levels of ACP1, mtACP3, and mtACP1 were positively correlated with those of ACP3, ACP4, and mtACP2, respectively. Analysis of cis-elements showed that seven motifs (CACTFTPPCA1, DOFCOREZM, GT1CONSENSUS, CAATBOX1, ARR1AT, POLLEN1LELAT52, and GATABOX) related to tissue-specific, ABA, and light-mediated gene regulation were ubiquitous in all ACPs investigated, which shed new light on the regulation patterns of these central enzymatic partners of the FAS system. This study presents a thorough overview of angiosperm ACP gene families and provides informative clues for the functional characterization of plant ACPs in the future.

On Deleterious Mutations in Perennials: Inbreeding Depression, Mutation Load, and Life-History Evolution.

AbstractIn angiosperms, perennials typically present much higher levels of inbreeding depression than annuals. One hypothesis to explain this pattern stems from the observation that inbreeding depression is expressed across multiple life stages in angiosperms. It posits that increased inbreeding depression in more long-lived species could be explained by differences in the way mutations affect fitness, through the life stages at which they are expressed. In this study, we investigate this hypothesis. We combine a physiological growth model and multilocus population genetics approaches to describe a full genotype-to-phenotype-to-fitness map. We study the behavior of mutations affecting growth or survival and explore their consequences in terms of inbreeding depression and mutation load. Although our results agree with empirical data only within a narrow range of conditions, we argue that they may point us toward the type of traits capable of generating high inbreeding depression in long-lived species-that is, traits under sufficiently strong selection, on which selection decreases sharply as life expectancy increases. Then we study the role deleterious mutations maintained at mutation-selection balance may play in the joint evolution of growth and survival strategies.

Interactions between bacteria and fungi in macrophyte leaf litter decomposition.

Microbes play an important role in decomposition of macrophytes in shallow lakes, and the process can be greatly affected by bacteria-fungi interactions in response to material composition and environmental conditions. In this study, microbes involved in the decomposition of leaf litter from three macrophyte species, Zizania latifolia, Hydrilla verticillata and Nymphoides peltata, were analysed at temperatures of 5, 15 and 25 °C. Results indicate that the decomposition rate was affected by temperature. Bacterial alpha diversity increased significantly along the time, while both temperature and plant species had a significant impact on the bacterial community, and plant type was shown to be the most important driving factor for the fungal community. The cosmopolitan bacterial taxa affiliated with Gammaproteobacteria, Bacteroidetes, Deltaproteobacteria, Firmicutes and Spirochaetes were key species in the investigated ecological networks, demonstrating significant co-occurrence or co-exclusion relationships with Basidiomycota and Ascomycota, according to different macrophyte species. This study indicates that bacteria involved in the decomposition of macrophyte leaf litter are more sensitive to temperature variance, and that fungi have a higher specificity to the composition of plant materials. The nutrient content of Hydrilla verticillata promoted a positive bacteria-fungi interaction, thereby accelerating the decomposition and re-circulation of leaf litter.

Evidence for Middle Miocene origin and morphological evolutionary stasis in a Barbacenia Inselberg clade (Velloziaceae).

Atlantic Forest Inselbergs (AFI) and Campos Rupestres (CR) are mountains and highlands of eastern South America, relatively poorly studied and highly threatened, which display extraordinary levels of plant endemism and richness. In spite of their geographical and environmental differences, the origin of the flora of CR and AFI are likely linked to each other, because several plant clades are distributed across both ecosystems. In addition to these studies, little has been investigated about the historical biogeographical connections between AFI and CR and most evolutionary studies are restricted to CR. Barbacenia (Velloziaceae) is widely spread and nearly endemic to the AFI and CR outcrops and thus represent an ideal system to study the biogeographical connections between CR and AFI. Besides, given the remarkable diversity of Barbacenia in CR compared to AFI, it appears that different factors were important drivers in the diversification of Barbacenia lineages, likely leading to different patterns of morphological diversification. Here, we integrate phylogenetic, biogeographic and morphological approaches to: (i) address whether AFI species of Barbacenia are monophyletic and thus a single colonization of AFI can be inferred; (ii) understand the timing and geographical origin of CR and AFI clades; (iii) compare morphological diversity between Barbacenia from AFI and CR under the hypothesis that these two systems have experienced similar levels of morphological diversification during their evolutionary history. To this end, we presented a phylogeny inferred using plastid (atpB-rbcL, trnH-psbA and trnL-trnF) and nuclear (ITS) markers and a complete sampling of AFI Barbacenia, estimated divergence times, reconstructed the ancestral areas of Barbacenia clades and compared their morphological diversity based on a dataset of 16 characters. Our results provided evidence for a diversification of Barbacenia from the Middle Miocene to Pleistocene, as suggested in previous studies. We suggest that stepping-stone dispersal across mountaintops in interplay with paleovegetation dynamics during the global Miocene cooling and Pleistocene climatic oscillations may played an important role in the range expansion of modern AFI Barbacenia lineages. Finally, our results also showed a significant differences in morphological diversity between AFI and CR clades, suggesting a long-term morphological stasis in AFI species.

Phylogeny, divergence times, and diversification in Calophyllaceae: Linking key characters and habitat changes to the evolution of Neotropical Calophylleae.

The clusioid clade comprises five monophyletic families: Bonnetiaceae, Calophyllaceae, Clusiaceae s.s., Hypericaceae, and Podostemaceae. Even though the circumscription of these families is well established, phylogenetic relationships within some families remain unresolved. This study aims to infer phylogenetic relationships within the Neotropical Calophylleae based on a broad sampling of taxa and a multilocus approach. We then use our phylogenetic framework as basis to investigate the evolution and biogeography of Calophylleae and diversification shifts in Calophyllaceae. To reconstruct the phylogeny of the Neotropical Calophylleae, we used five plastid (matK, ndhF, rbcL, psbA-trnH, and trnK), two mitochondrial (matR and rps3), and two nuclear (EMB2765 and ITS) markers, including previously published and newly generated sequences. We sampled 74 species, increasing sampling of Neotropical taxa by 500%. Our phylogenetic hypothesis for Calophyllaceae provides additional support for the monophyly of all genera and allowed us to identify four main clades: Calophyllum, Kayea, Mammea, and the Neotropical clade. The Neotropical clade includes three main lineages, a small clade composed of Clusiella and Marila, and a large HaCaKi clade (i.e., Haplocarpa, Caraipa, and Kilmeyera) that is sister to Mahurea exstipulata. The evolution of three morphological traits (i.e., fleshy fruits, anther glands, and winged seeds) were shown to be associated with changes in evolutionary dynamics in Calophyllaceae, while a biome shift was detected in Kielmeyera, affecting net diversification within this genus. Major geological and climatic events such as the Andean uplift and a gradual decrease in temperatures seem to have influenced diversification rates within the Neotropical Calophylleae.

Phylogenomics of Gesneriaceae using targeted capture of nuclear genes.

Gesneriaceae (ca. 3400 species) is a pantropical plant family with a wide range of growth form and floral morphology that are associated with repeated adaptations to different environments and pollinators. Although Gesneriaceae systematics has been largely improved by the use of Sanger sequencing data, our understanding of the evolutionary history of the group is still far from complete due to the limited number of informative characters provided by this type of data. To overcome this limitation, we developed here a Gesneriaceae-specific gene capture kit targeting 830 single-copy loci (776,754 bp in total), including 279 genes from the Universal Angiosperms-353 kit. With an average of 557,600 reads and 87.8% gene recovery, our target capture was successful across the family Gesneriaceae and also in other families of Lamiales. From our bait set, we selected the most informative 418 loci to resolve phylogenetic relationships across the entire Gesneriaceae family using maximum likelihood and coalescent-based methods. Upon testing the phylogenetic performance of our baits on 78 taxa representing 20 out of 24 subtribes within the family, we showed that our data provided high support for the phylogenetic relationships among the major lineages, and were able to provide high resolution within more recent radiations. Overall, the molecular resources we developed here open new perspectives for the study of Gesneriaceae phylogeny at different taxonomical levels and the identification of the factors underlying the diversification of this plant group.

Species and geographic specificity between endophytic fungi and host supported by parasitic Cynomorium songaricum and its host Nitraria tangutorum distributed in desert.

This study was aimed to investigate whether host plant species and lifestyles, and environmental conditions in the desert affect endophytic fungi composition. Endophytic fungal communities from parasitic plant Cynomorium songaricum and its host Nitraria tangutorum were investigated from three sites including Tonggu Naoer, Xilin Gaole, and Guazhou in Tengger and Badain Jaran Deserts in China using the next-generation sequencing of a ribosomal RNA gene region. Similarity and difference in endophytic fungal composition from different geographic locations were evaluated through multivariate statistical analysis. It showed that plant genetics was a deciding factor affecting endophytic fungal composition even when C. songaricum and N. tangutorum grow together tightly. Not only that, the fungal composition was also greatly affected by the local environment and rainfall. However, the distribution and richness of fungal species indicated that the geographical distance exerted little influence on characterizing the fungal composition. Overall, the findings suggested that plant species, parasitic or non-parasitic lifestyles of the plant, and local environment strongly affected the number and diversity of the endophytic fungal species, which may provide valuable insights into the microbe ecology, symbiosis specificity, and the tripartite relationship among parasitic plant, host, and endophytic fungi, especially under desert environment.

Simultaneous diversification of Polypodiales and angiosperms in the Mesozoic.

Comprising about 82% of the extant fern species diversity, Polypodiales are generally believed to have diversified in the Late Cretaceous. We estimated the divergence times of Polypodiales using both penalized likelihood and Bayesian methods, based on a dataset consisting of 208 plastomes representing all 28 families and 14 fossil constraints reflecting current interpretations of fossil record. Our plastome phylogeny recovered the same six major lineages as a recent nuclear phylogeny, but the position of Dennstaedtiineae was different. The present phylogeny showed high resolution of relationships among the families of Polypodiales, especially among those forming the Aspleniineae. The divergence time estimates supported the most recent common ancestor of Polypodiales and its closest relative dating back to the Triassic, establishment of the major lineages in the Jurassic, and a likely accelerated radiation during the late Jurassic and the Early Cretaceous. The estimated divergence patterns of Polypodiales and angiosperms converge to a scenario in which their main lineages were established simultaneously shortly before the onset of the Cretaceous Terrestrial Revolution, and further suggest a pre-Cretaceous hidden history for both lineages. The pattern of simultaneous diversifications shown here elucidate an important gap in our understanding of the Terrestrial Revolution that shaped today's ecosystems.

Re"CYC"ling molecular regulators in the evolution and development of flower symmetry.

Flower forms are both highly diverse and multifaceted. As well as varying in colour, size, organ number, and much more, flowers show different types of symmetry. Floral symmetry can be grouped into three main categories: asymmetry, bilateral symmetry and radial symmetry, characterised by zero, one, and multiple planes of symmetry, respectively. This review will first explore floral symmetry from a classical morphological view, then from a modern molecular perspective. The recent molecular work on symmetry in monocots and eudicots will be discussed, followed by an in-depth discussion into the evolution of CYC genes, particularly in the capitulum of the sunflower family (Asteraceae). Whilst recent studies on non-model species are helping to bring new light to this field, more species coverage is required to understand how traits such as bilateral symmetry have evolved so many times, and whether the same molecular regulators were recruited for this function.

Implications of region-specific gene expression for development of the partially fused petunia corolla.

Angiosperm petal fusion (sympetaly) has evolved multiple times independently and is associated with increased specificity between plants and their pollinators. To uncover developmental genetic changes that might have led to the evolution of sympetaly in the asterid core eudicot genus Petunia (Solanaceae), we carried out global and fine-scale gene expression analyses in different regions of the corolla. We found that, despite several similarities with the choripetalous model species Arabidopsis thaliana in the proximal-distal transcriptome, the Petunia axillaris fused and proximal corolla tube expresses several genes that in A. thaliana are associated with the distal petal region. This difference aligns with variation in petal shape and fusion across ontogeny of the two species. Moreover, differential gene expression between the unfused lobes and fused tube of P. axillaris petals revealed three strong candidate genes for sympetaly based on functional annotation in organ boundary specification. Partial silencing of one of these, the HANABA TARANU (HAN)-like gene PhGATA19, resulted in reduced fusion of Petunia hybrida petals, with silencing of both PhGATA19 and its close paralog causing premature plant senescence. Finally, detailed expression analyses for the previously characterized organ boundary gene candidate NO APICAL MERISTEM (NAM) supports the hypothesis that it establishes boundaries between most P. axillaris floral organs, with the exception of boundaries between petals.

Evolutionary significance of amino acid permease transporters in 17 plants from Chlorophyta to Angiospermae.

BACKGROUND: Nitrogen is an indispensable nutrient for plant growth. It is used and transported in the form of amino acids in living organisms. Transporting amino acids to various parts of plants requires relevant transport proteins, such as amino acid permeases (AAPs), which were our focus in this study. RESULTS: We found that 5 AAP genes were present in Chlorophyte species and more AAP genes were predicted in Bryophyta and Lycophytes. Two main groups were defined and group I comprised 5 clades. Our phylogenetic analysis indicated that the origin of clades 2, 3, and 4 is Gymnospermae and that these clades are closely related. The members of clade 1 included Chlorophyta to Gymnospermae. Group II, as a new branch consisting of non-seed plants, is first proposed in our research. Our results also indicated that the AAP family was already present in Chlorophyta and then expanded accompanying the development of vasculature. Concurrently, the AAP family experienced multiple duplication events that promoted the generation of new functions and differentiation of sub-functions. CONCLUSIONS: Our findings suggest that the AAP gene originated in Chlorophyta, and some non-seed AAP genes clustered in one group. A second group, which contained plants of all evolutionary stages, indicated the evolution of AAPs. These new findings can be used to guide future research.