Search

Background: Scrophularia L., a genus of the Scrophulariaceae, is a group of important medicinal plants used for eliminating heat and detoxifying. East Asia has an abundance of potentially medicinal Scrophularia species, and it serves as a secondary diversity center of the genus. However, the genomic resources available for germplasm identification and pharmaceutical exploration of East Asian Scrophularia are insufficient, hindering its commercial and industrial development. Additionally, the interspecific relationships of most East Asian Scrophularia species remain unclear. Methods: In this study, we sequenced the leaves of 25 East Asian species of the genus Scrophularia, assembled and annotated the complete chloroplast genomes, and subsequently performed comparative and phylogenetic analyses on these genomes. Results and discussion: The conserved plastome length of these 25 species ranged from 151,582 bp to 153,239 bp, containing a total of 132 coding genes, including 18 duplicated genes and 114 unique genes. Through genome alignment of these 25 species, 38-53 repeated sequences and 7 shared SSRs were identified, along with regions with high nucleotide polymorphism (Pi), which could potentially serve as molecular markers for species identification. The genome structure, gene content, and arrangement showed conservation, while variations were observed in the IR boundary regions and IGS. Phylogenetic inferences based on whole plastomes or on coding sequences (CDS) only yielded congruent results. We categorized the 25 East Asian Scrophularia species into six distinct clades and further explored their interspecies relationships using morphological characteristics, such as flower color, the relative position of stamens and corolla, and plant height. This could lay a genetic basis for future resource development of Scrophularia in East Asia. [ABSTRACT FROM AUTHOR]

Alisma L. is a medicinally important genus of aquatic and wetland plants consisting of c. 10 recognized species. However, largely due to polyploidy and limited taxon and gene sampling, the phylogenomic relationships of Alisma remain challenging. In this study, we sequenced 34 accessions of Alismataceae, including eight of the ten species of Alisma, one species of Echinodorus and one species of Luronium, to perform comparative analyses of plastid genomes and phylogenetic analyses. Comparative analysis of plastid genomes revealed high sequence similarity among species within the genus. Our study analyzed structural changes and variations in the plastomes of Alisma, including IR expansion or contraction, and gene duplication or loss. Phylogenetic results suggest that Alisma is monophyletic, and constitutes four groups: (1) A. lanceolatum and A. canaliculatum; (2) the North American clade of A. subcordatum and A. triviale; (3) A. wahlenbergii and A. gramineum; and (4) A. plantago-aquatica from Eurasia and northern Africa with the eastern Asian A. orientale nested within it. Hence the results challenge the recognition of A. orientale as a distinct species and raise the possibility of treating it as a synonym of the widespread A. plantago-aquatica. The well-known Alismatis Rhizoma (Zexie) in Chinese medicine was likely derived from the morphologically variable Alisma plantago-aquatica throughout its long history of cultivation in Asia. The plastome phylogenetic results also support the tetraploid A. lanceolatum as the likely maternal parent of the hexaploid eastern Asian A. canaliculatum [ABSTRACT FROM AUTHOR]

Reconstructing evolutionary trajectories and transitions that have shaped floral diversity relies heavily on the phylogenetic framework on which traits are modelled. In this study, we focus on the angiosperm order Ranunculales, sister to all other eudicots, to unravel higher-level relationships, especially those tied to evolutionary transitions in flower symmetry within the family Papaveraceae. This family presents an astonishing array of floral diversity, with actinomorphic, disymmetric (two perpendicular symmetry axes), and zygomorphic flowers. We generated nuclear and plastid datasets using the Angiosperms353 universal probe set for target capture sequencing (of 353 single-copy nuclear ortholog genes), together with publicly available transcriptome and plastome data mined from open-access online repositories. We relied on the fossil record of the order Ranunculales to date our phylogenies and to establish a timeline of events. Our phylogenomic workflow shows that nuclear-plastid incongruence accompanies topological uncertainties in Ranunculales. A cocktail of incomplete lineage sorting, post-hybridization introgression, and extinction following rapid speciation most likely explain the observed knots in the topology. These knots coincide with major floral symmetry transitions and thus obscure the order of evolutionary events. [ABSTRACT FROM AUTHOR]

Introduction: The genus Sanicula L. is a taxonomically complicated taxa within Apiaceae, as its high variability in morphology. Although taxonomists have performed several taxonomic revisions for this genus, the interspecific relationships and species boundaries have not been satisfactorily resolved, especially for those endemic to China. This study mainly focused on S. giraldii var. ovicalycina, S. tienmuensis var. pauciflora, and S. orthacantha var. stolonifera and also described two new members of the genus. Methods: We newly sequenced sixteen plastomes from nine Sanicula species. Combined with eleven plastomes previously reported by us and one plastome downloaded, we performed a comprehensively plastid phylogenomics analysis of 21 Sanicula taxa. Results and Discussion: The comparative results showed that 21 Sanicula plastomes in their structure and features were highly conserved and further justified that two new species were indeed members of Sanicula. Nevertheless, eleven mutation hotspot regions were still identified. Phylogenetic analyses based on plastome data and the ITS sequences strongly supported that these three varieties were clearly distant from three type varieties. The results implied that these three varieties should be considered as three independent species, which were further justified by their multiple morphological characters. Therefore, revising these three varieties into three independent species was reasonable and convincing. Moreover, we also identified and described two new Sanicula species (S. hanyuanensis and S. langaoensis) from Sichuan and Shanxi, China, respectively. Based on their distinct morphological characteristics and molecular phylogenetic analysis, two new species were included in Sanicula. In summary, our study impelled the revisions of Sanicula members and improved the taxonomic system of the genus. [ABSTRACT FROM AUTHOR]

The hyperdiverse orchid genus Bulbophyllum is the second largest genus of flowering plants and exhibits a pantropical distribution with a center of diversity in tropical Asia. The only Bulbophyllum section with a center of diversity inAustralasia is sect. Adelopetalum. However, the phylogenetic placement, interspecific relationships, and spatio-temporal evolution of this section remain largely unclear. To infer broad-level relationships within Bulbophyllum, and interspecific relationshipswithinsect. Adelopetalum, a genome skimmingdatasetwasgenerated for 89 samples, which yielded 70 plastid coding regions and a nuclear ribosomal DNA cistron. For 18 additional samples, Sanger data fromtwo plastid loci (matK and ycf1) and nuclear ITSwere added using a supermatrix approach. The study provided newinsights into broad-level relationships in Bulbophyllum, including phylogenetic evidence for the non-monophyly of sections Beccariana, Brachyantha, Brachypus, Cirrhopetaloides, Cirrhopetalum, Desmosanthes, Minutissima, Oxysepala, Polymeres, and Sestochilos. Section Adelopetalum and sect. Minutissima s.s. formed a highly supported clade that was resolved as a sister group to the remainder of the genus. Divergence time estimations based on a relaxed molecular clock model placed the origin of Bulbophyllum in the Early Oligocene (ca. 33.2 Ma) and sect. Adelopetalumin the Late Oligocene (ca. 23.6 Ma). Ancestral range estimations based on a BAYAREALIKE model identified the Australian continent as the ancestral area of the sect. Adelopetalum. The section underwent crown diversification fromthemid-Miocene to the late Pleistocene, predominantly in continental Australia. At least two independent long-distance dispersal events were inferred eastward from the Australian continent to New Zealand and to New Caledonia from the early Pliocene onwards, likely mediated by predominantly westerly winds of the Southern hemisphere. Retraction and fragmentation of the eastern Australian rainforests from the early Miocene onwards are likely drivers of lineage divergence within sect. Adelopetalum facilitating allopatric speciation. [ABSTRACT FROM AUTHOR]

Equisetum is one of the oldest extant group vascular plants and is considered to be the key to understanding vascular plant evolution. Equisetum is distributed almost all over the world and has a high degree of adaptability to different environments. Despite the fossil record of horsetails (Equisetum, Equisetaceae) dating back to the Carboniferous, the phylogenetic relationship of this genus is not well, and the chloroplast evolution in Equisetum remains poorly understood. In order to fill this gap, we sequenced, assembled, and annotated the chloroplast genomes of 12 species of Equisetum, and compared them to 13 previously published vascular plants chloroplast genomes to deeply examine the plastome evolutionary dynamics of Equisetum. The chloroplast genomes have a highly conserved quadripartite structure across the genus, but these chloroplast genomes have a lower GC content than other ferns. The size of Equisetum plastomes ranges from 130,773 bp to 133,684 bp and they encode 130 genes. Contraction/expansion of IR regions and the number of simple sequences repeat regions underlie large genomic variations in size among them. Comparative analysis revealed we also identified 13 divergence hotspot regions. Additionally, the genes accD and ycf1 can be used as potential DNA barcodes for the identification and phylogeny of the genus Equisetum. Twelve photosynthesisrelated genes were specifically selected in Equisetum. Comparative genomic analyses implied divergent evolutionary patterns between Equisetum and other ferns. Phylogenomic analyses and molecular dating revealed a relatively distant phylogenetic relationship between Equisetum and other ferns, supporting the division of pteridophyte into Lycophytes, Equisetaceae and ferns. The results show that the chloroplast genome can be used to solve phylogenetic problems within or between Equisetum species, and also provide genomic resources for the study of Equisetum systematics and evolution. [ABSTRACT FROM AUTHOR]

Introduction: Pectobacterium cacticida was identified as the causative agent of soft rot disease in cacti. Due to a high potential of spread in the face of global warming, the species poses a significant threat to horticultural and crop industry. The aim of this study was to revise the genomic, physiology and virulence characteristics of P. cacticida and update its phylogenetic position within the Pectobacterium genus. Methods: Whole genome sequences of five P. cacticida strains were obtained and subjected to comprehensive genomic and phylogenomic data analyses. We assessed the presence of virulence determinants and genes associated with host and environmental adaptation. Lipidomic analysis, as well as biochemical and phenotypic assays were performed to correlate genomic findings. Results: Phylogenomic analysis revealed that P. cacticida forms a distinct lineage within the Pectobacterium genus. Genomic evaluation uncovered 516 unique proteins, most of which were involved in cellular metabolism. They included genes of carbohydrate metabolism and transport and ABC transporters. The main differing characteristics from other Pectobacterium species were the lack of a myo-inositol degradation pathway and the presence of the malonate decarboxylase gene. All tested strains were pathogenic towards Opuntia spp., chicory, Chinese cabbage, and potato, but exhibited only mild pathogenicity towards carrot. Discussion: This study sheds light into the genomic characteristics of P. cacticida and highlights the pathogenic potential of the species. Unique genes found in P. cacticida genomes possibly enhance the species' survival and virulence. Based on phylogenomic analyses, we propose the reclassification of P. cacticida to a new genus, Alcorniella comb. nov. [ABSTRACT FROM AUTHOR]

Reconstructing evolutionary trajectories and transitions that have shaped floral diversity relies heavily on the phylogenetic framework on which traits are modelled. In this study, we focus on the angiosperm order Ranunculales, sister to all other eudicots, to unravel higher-level relationships, especially those tied to evolutionary transitions in flower symmetry within the family Papaveraceae. This family presents an astonishing array of floral diversity, with actinomorphic, disymmetric (two perpendicular symmetry axes), and zygomorphic flowers. We generated nuclear and plastid datasets using the Angiosperms353 universal probe set for target capture sequencing (of 353 single-copy nuclear ortholog genes), together with publicly available transcriptome and plastome data mined from open-access online repositories. We relied on the fossil record of the order Ranunculales to date our phylogenies and to establish a timeline of events. Our phylogenomic workflow shows that nuclear-plastid incongruence accompanies topological uncertainties in Ranunculales. A cocktail of incomplete lineage sorting, post-hybridization introgression, and extinction following rapid speciation most likely explain the observed knots in the topology. These knots coincide with major floral symmetry transitions and thus obscure the order of evolutionary events.

IntroductionPectobacterium cacticida was identified as the causative agent of soft rot disease in cacti. Due to a high potential of spread in the face of global warming, the species poses a significant threat to horticultural and crop industry. The aim of this study was to revise the genomic, physiology and virulence characteristics of P. cacticida and update its phylogenetic position within the Pectobacterium genus.MethodsWhole genome sequences of five P. cacticida strains were obtained and subjected to comprehensive genomic and phylogenomic data analyses. We assessed the presence of virulence determinants and genes associated with host and environmental adaptation. Lipidomic analysis, as well as biochemical and phenotypic assays were performed to correlate genomic findings.ResultsPhylogenomic analysis revealed that P. cacticida forms a distinct lineage within the Pectobacterium genus. Genomic evaluation uncovered 516 unique proteins, most of which were involved in cellular metabolism. They included genes of carbohydrate metabolism and transport and ABC transporters. The main differing characteristics from other Pectobacterium species were the lack of a myo-inositol degradation pathway and the presence of the malonate decarboxylase gene. All tested strains were pathogenic towards Opuntia spp., chicory, Chinese cabbage, and potato, but exhibited only mild pathogenicity towards carrot.DiscussionThis study sheds light into the genomic characteristics of P. cacticida and highlights the pathogenic potential of the species. Unique genes found in P. cacticida genomes possibly enhance the species’ survival and virulence. Based on phylogenomic analyses, we propose the reclassification of P. cacticida to a new genus, Alcorniella comb. nov.

Introduction: The recently established Linderniaceae, separated from the traditionally defined Scrophulariaceae, is a taxonomically complicated family. Although previous phylogenetic studies based on a few short DNA markers have made great contributions to the taxonomy of Linderniaceae, limited sampling and low resolution of the phylogenetic tree have failed to resolve controversies between some generic circumscriptions. The plastid genome exhibits a powerful ability to solve phylogenetic relationships ranging from shallow to deep taxonomic levels. To date, no plastid phylogenomic studies have been carried out in Linderniaceae. Methods: In this study, we newly sequenced 26 plastid genomes of Linderniaceae, including eight genera and 25 species, to explore the phylogenetic relationships and genome evolution of the family through plastid phylogenomic and comparative genomic analyses. Results: The plastid genome size of Linderniaceae ranged from 152,386 bp to 154,402 bp, exhibiting a typical quartile structure. All plastomes encoded 114 unique genes, comprising 80 protein-coding genes, 30 tRNA genes, and four rRNA genes. The inverted repeat regions were more conserved compared with the single-copy regions. A total of 1803 microsatellites and 1909 long sequence repeats were identified, and five hypervariable regions (petN-psbM, rps16-trnQ, rpl32-trnL, rpl32, and ycf1) were screened out. Most protein-coding genes were relatively conserved, with only the ycf2 gene found under positive selection in a few species. Phylogenomic analyses confirmed that Linderniaceae was a distinctive lineage and revealed that the presently circumscribed Vandellia and Torenia were non-monophyletic. Discussion: Comparative analyses showed the Linderniaceae plastomes were highly conservative in terms of structure, gene order, and gene content. Combining morphological and molecular evidence, we supported the newly established Yamazakia separating from Vandellia and the monotypic Picria as a separate genus. These findings provide further evidence to recognize the phylogenetic relationships among Linderniaceae and new insights into the evolution of the plastid genomes. [ABSTRACT FROM AUTHOR]

Introgression can produce novel genetic variation in organisms that hybridize. Sympatric species pairs in the carnivorous plant genus Sarracenia L. frequently hybridize, and all known hybrids are fertile. Despite being a desirable system for studying the evolutionary consequences of hybridization, the extent to which introgression occurs in the genus is limited to a few species in only two field sites. Previous phylogenomic analysis of Sarracenia estimated a highly resolved species tree from 199 nuclear genes, but revealed a plastid genome that is highly discordant with the species tree. Such cytonuclear discordance could be caused by chloroplast introgression (i.e. chloroplast capture) or incomplete lineage sorting (ILS). To better understand the extent to which introgression is occurring in Sarracenia, the chloroplast capture and ILS hypotheses were formally evaluated. Plastomes were assembled de-novo from sequencing reads generated from 17 individuals in addition to reads obtained from the previous study. Assemblies of 14 whole plastomes were generated and annotated, and the remaining fragmented assemblies were scaffolded to these whole-plastome assemblies. Coding sequence from 79 homologous genes were aligned and concatenated for maximum-likelihood phylogeny estimation. The plastome tree is extremely discordant with the published species tree. Plastome trees were simulated under the coalescent and tree distance from the species tree was calculated to generate a null distribution of discordance that is expected under ILS alone. A t-test rejected the null hypothesis that ILS could cause the level of discordance seen in the plastome tree, suggesting that chloroplast capture must be invoked to explain the discordance. Due to the extreme level of discordance in the plastome tree, it is likely that chloroplast capture has been common in the evolutionary history of Sarracenia. [ABSTRACT FROM AUTHOR]

The Asteraceae is the largest angiosperm family with more than 25,000 species. Individual studies have shown that MADS-box and TCP transcription factors are regulators of the development and symmetry of flowers, contributing to their iconic flower-head (capitulum) and floret. However, a systematic study of MADSbox and TCP genes across the Asteraceae is lacking. We performed a comparative analysis of genome sequences of 33 angiosperm species including our de novo assembly of diploid sexual dandelion (Taraxacum officinale) and 11 other Asteraceae to investigate the lineage-specific evolution of MADS-box and TCP genes in the Asteraceae. We compared the phylogenomic results of MADS-box and TCP genes with their expression in T. officinale floral tissues at different developmental stages to demonstrate the regulation of genes with Asteraceae-specific attributes. Here, we show that MADS-box MIKCc and TCP-CYCLOIDEA (CYC) genes have expanded in the Asteraceae. The phylogenomic analysis identified AGAMOUS-like (AG-like: SEEDSTICK [STK]-like), SEPALATA-like (SEP3-like), and TCP-PROLIFERATING CELL FACTOR (PCF)-like copies with lineage-specific genomic contexts in the Asteraceae, Cichorioideae, or dandelion. Different expression patterns of some of these gene copies suggest functional divergence. We also confirm the presence and revisit the evolutionary history of previously named "Asteraceae-Specific MADSbox genes (AS-MADS)." Specifically, we identify non-Asteraceae homologs, indicating a more ancient origin of this gene clade. Syntenic relationships support that AS-MADS is paralogous to FLOWERING LOCUS C (FLC) as demonstrated by the shared ancient duplication of FLC and SEP3. [ABSTRACT FROM AUTHOR]

Elimination of chemically synthesized pesticides, such as fungicides and nematicides, in agricultural products is a key to successful practice of the Vietnamese agriculture. We describe here the route for developing successful biostimulants based on members of the Bacillus subtilis species complex. A number of endospore-forming Gram-positive bacterial strains with antagonistic action against plant pathogens were isolated from Vietnamese crop plants. Based on their draft genome sequence, thirty of them were assigned to the Bacillus subtilis species complex. Most of them were assigned to the species Bacillus velezensis. Whole genome sequencing of strains BT2.4 and BP1.2A corroborated their close relatedness to B. velezensis FZB42, the model strain for Gram-positive plant growth-promoting bacteria. Genome mining revealed that at least 15 natural product biosynthesis gene clusters (BGCs) are well conserved in all B. velezensis strains. In total, 36 different BGCs were identified in the genomes of the strains representing B. velezensis, B. subtilis, Bacillus tequilensis, and Bacillus. altitudinis. In vitro and in vivo assays demonstrated the potential of the B. velezensis strains to enhance plant growth and to suppress phytopathogenic fungi and nematodes. Due to their promising potential to stimulate plant growth and to support plant health, the B. velezensis strains TL7 and S1 were selected as starting material for the development of novel biostimulants, and biocontrol agents efficient in protecting the important Vietnamese crop plants black pepper and coffee against phytopathogens. The results of the large-scale field trials performed in the Central Highlands in Vietnam corroborated that TL7 and S1 are efficient in stimulating plant growth and protecting plant health in large-scale applications. It was shown that treatment with both bioformulations resulted in prevention of the pathogenic pressure exerted by nematodes, fungi, and oomycetes, and increased harvest yield in coffee, and pepper. [ABSTRACT FROM AUTHOR]

Trichosporeae is the largest and most taxonomically difficult tribe of Gesneriaceae due to its diverse morphology. Previous studies have not clarified the phylogenetic relationships within this tribe on several DNA markers, including the generic relationships within subtribes. Recently, plastid phylogenomics have been successfully employed to resolve the phylogenetic relationships at different taxonomic levels. In this study, plastid phylogenomics were used to explore the relationships within Trichosporeae. Eleven plastomes of Hemiboea were newly reported. Comparative analyses, phylogeny and morphological character evolution within Trichosporeae were conducted on 79 species representing seven subtribes. The Hemiboea plastomes range from 152,742 bp to 153,695 bp in length. Within Trichosporeae, the sampled plastomes range from 152,196 bp to 156,614 bp and GC content from 37.2% to 37.8%. A total of 121-133 genes were annotated in each species, including 80-91 proteincoding genes, 34-37 tRNA genes, and 8 rRNA genes. The contraction and expansion of IR borders were not detected, and gene rearrangements and inversions did not occur. The 13 hypervariable regions were proposed as the potential molecular markers for species identification. A total of 24,299 SNPs and 3,378 indels were inferred, and most of the SNPs were functionally missense and silent variations. There were 1968 SSRs, 2055 tandem repeats and 2802 dispersed repeats. The RSCU and ENC values indicated that the codon usage pattern was conserved in Trichosporeae. Both the phylogenetic frameworks based on the whole plastome and 80 CDSs were basically concordant. The sister relationships between Loxocarpinae and Didymocarpinae were confirmed, and Oreocharis was a sister group of Hemiboea with high support. The morphological characters showed a complex evolutionary pattern of Trichosporeae. Our findings may contribute to future research on genetic diversity, morphological evolutionary patterns, and conservation of the tribe Trichosporeae. [ABSTRACT FROM AUTHOR]

Introgression can produce novel genetic variation in organisms that hybridize. Sympatric species pairs in the carnivorous plant genus Sarracenia L. frequently hybridize, and all known hybrids are fertile. Despite being a desirable system for studying the evolutionary consequences of hybridization, the extent to which introgression occurs in the genus is limited to a few species in only two field sites. Previous phylogenomic analysis of Sarracenia estimated a highly resolved species tree from 199 nuclear genes, but revealed a plastid genome that is highly discordant with the species tree. Such cytonuclear discordance could be caused by chloroplast introgression (i.e. chloroplast capture) or incomplete lineage sorting (ILS). To better understand the extent to which introgression is occurring in Sarracenia, the chloroplast capture and ILS hypotheses were formally evaluated. Plastomes were assembled de-novo from sequencing reads generated from 17 individuals in addition to reads obtained from the previous study. Assemblies of 14 whole plastomes were generated and annotated, and the remaining fragmented assemblies were scaffolded to these whole-plastome assemblies. Coding sequence from 79 homologous genes were aligned and concatenated for maximum-likelihood phylogeny estimation. The plastome tree is extremely discordant with the published species tree. Plastome trees were simulated under the coalescent and tree distance from the species tree was calculated to generate a null distribution of discordance that is expected under ILS alone. A t-test rejected the null hypothesis that ILS could cause the level of discordance seen in the plastome tree, suggesting that chloroplast capture must be invoked to explain the discordance. Due to the extreme level of discordance in the plastome tree, it is likely that chloroplast capture has been common in the evolutionary history of Sarracenia.

Chloroplasts and mitochondria each contain their own genomes, which have historically been and continue to be important sources of information for inferring the phylogenetic relationships among land plants. The organelles are predominantly inherited from the same parent, and therefore should exhibit phylogenetic concordance. In this study, we examine the mitochondrion and chloroplast genomes of 226 land plants to infer the degree of similarity between the organelles' evolutionary histories. Our results show largely concordant topologies are inferred between the organelles, aside from four wellsupported conflicting relationships that warrant further investigation. Despite broad patterns of topological concordance, our findings suggest that the chloroplast and mitochondrial genomes evolved with significant differences in molecular evolution. The differences result in the genes from the chloroplast and the mitochondrion preferentially clustering with other genes from their respective organelles by a program that automates selection of evolutionary model partitions for sequence alignments. Further investigation showed that changes in compositional heterogeneity are not always uniform across divergences in the land plant tree of life. These results indicate that although the chloroplast and mitochondrial genomes have coexisted for over 1 billion years, phylogenetically, they are still evolving sufficiently independently to warrant separate models of evolution. As genome sequencing becomes more accessible, research into these organelles' evolution will continue revealing insight into the ancient cellular events that shaped not only their history, but the history of plants as a whole. [ABSTRACT FROM AUTHOR]

Sapindales is an angiosperm order of high economic and ecological value comprising nine families, c. 479 genera, and c. 6570 species. However, family and subfamily relationships in Sapindales remain unclear, making reconstruction of the order's spatio-temporal and morphological evolution difficult. In this study, we used Angiosperms353 target capture data to generate the most densely sampled phylogenetic trees of Sapindales to date, with 448 samples and c. 85% of genera represented. The percentage of paralogous loci and allele divergence was characterized across the phylogeny, which was time-calibrated using 29 rigorously assessed fossil calibrations. All families were supported as monophyletic. Two core family clades subdivide the order, the first comprising Kirkiaceae, Burseraceae, and Anacardiaceae, the second comprising Simaroubaceae, Meliaceae, and Rutaceae. Kirkiaceae is sister to Burseraceae and Anacardiaceae, and, contrary to current understanding, Simaroubaceae is sister to Meliaceae and Rutaceae. Sapindaceae is placed with Nitrariaceae and Biebersteiniaceae as sister to the core Sapindales families, but the relationships between these families remain unclear, likely due to their rapid and ancient diversification. Sapindales families emerged in rapid succession, coincident with the climatic change of the Mid-Cretaceous Hothouse event. Subfamily and tribal relationships within the major families need revision, particularly in Sapindaceae, Rutaceae and Meliaceae. Much of the difficulty in reconstructing relationships at this level may be caused by the prevalence of paralogous loci, particularly in Meliaceae and Rutaceae, that are likely indicative of ancient gene duplication events such as hybridization and polyploidization playing a role in the evolutionary history of these families. This study provides key insights into factors that may affect phylogenetic reconstructions in Sapindales across multiple scales, and provides a state-of-the-art phylogenetic framework for further research. [ABSTRACT FROM AUTHOR]

The BAHD acyltransferase family is one of the largest enzyme families in flowering plants, containing dozens to hundreds of genes in individual genomes. Highly prevalent in angiosperm genomes, members of this family contribute to several pathways in primary and specialized metabolism. In this study, we performed a phylogenomic analysis of the family using 52 genomes across the plant kingdom to gain deeper insights into its functional evolution and enable function prediction. We found that BAHD expansion in land plants was associated with significant changes in various gene features. Using pre-defined BAHD clades, we identified clade expansions in different plant groups. In some groups, these expansions coincided with the prominence of metabolite classes such as anthocyanins (flowering plants) and hydroxycinnamic acid amides (monocots). Clade-wise motif-enrichment analysis revealed that some clades have novel motifs fixed on either the acceptor or the donor side, potentially reflecting historical routes of functional evolution. Co-expression analysis in rice and Arabidopsis further identified BAHDs with similar expression patterns, however, most co-expressed BAHDs belonged to different clades. Comparing BAHD paralogs, we found that gene expression diverges rapidly after duplication, suggesting that sub/neo-functionalization of duplicate genes occurs quickly via expression diversification. Analyzing co-expression patterns in Arabidopsis in conjunction with orthology-based substrate class predictions and metabolic pathway models led to the recovery of metabolic processes of most of the already-characterized BAHDs as well as definition of novel functional predictions for some uncharacterized BAHDs. Overall, this study provides new insights into the evolution of BAHD acyltransferases and sets up a foundation for their functional characterization. [ABSTRACT FROM AUTHOR]

Chloroplasts and mitochondria each contain their own genomes, which have historically been and continue to be important sources of information for inferring the phylogenetic relationships among land plants. The organelles are predominantly inherited from the same parent, and therefore should exhibit phylogenetic concordance. In this study, we examine the mitochondrion and chloroplast genomes of 226 land plants to infer the degree of similarity between the organelles’ evolutionary histories. Our results show largely concordant topologies are inferred between the organelles, aside from four well-supported conflicting relationships that warrant further investigation. Despite broad patterns of topological concordance, our findings suggest that the chloroplast and mitochondrial genomes evolved with significant differences in molecular evolution. The differences result in the genes from the chloroplast and the mitochondrion preferentially clustering with other genes from their respective organelles by a program that automates selection of evolutionary model partitions for sequence alignments. Further investigation showed that changes in compositional heterogeneity are not always uniform across divergences in the land plant tree of life. These results indicate that although the chloroplast and mitochondrial genomes have coexisted for over 1 billion years, phylogenetically, they are still evolving sufficiently independently to warrant separate models of evolution. As genome sequencing becomes more accessible, research into these organelles’ evolution will continue revealing insight into the ancient cellular events that shaped not only their history, but the history of plants as a whole.

Sapindales is an angiosperm order of high economic and ecological value comprising nine families, c. 479 genera, and c. 6570 species. However, family and subfamily relationships in Sapindales remain unclear, making reconstruction of the order’s spatio-temporal and morphological evolution difficult. In this study, we used Angiosperms353 target capture data to generate the most densely sampled phylogenetic trees of Sapindales to date, with 448 samples and c. 85% of genera represented. The percentage of paralogous loci and allele divergence was characterized across the phylogeny, which was time-calibrated using 29 rigorously assessed fossil calibrations. All families were supported as monophyletic. Two core family clades subdivide the order, the first comprising Kirkiaceae, Burseraceae, and Anacardiaceae, the second comprising Simaroubaceae, Meliaceae, and Rutaceae. Kirkiaceae is sister to Burseraceae and Anacardiaceae, and, contrary to current understanding, Simaroubaceae is sister to Meliaceae and Rutaceae. Sapindaceae is placed with Nitrariaceae and Biebersteiniaceae as sister to the core Sapindales families, but the relationships between these families remain unclear, likely due to their rapid and ancient diversification. Sapindales families emerged in rapid succession, coincident with the climatic change of the Mid-Cretaceous Hothouse event. Subfamily and tribal relationships within the major families need revision, particularly in Sapindaceae, Rutaceae and Meliaceae. Much of the difficulty in reconstructing relationships at this level may be caused by the prevalence of paralogous loci, particularly in Meliaceae and Rutaceae, that are likely indicative of ancient gene duplication events such as hybridization and polyploidization playing a role in the evolutionary history of these families. This study provides key insights into factors that may affect phylogenetic reconstructions in Sapindales across multiple scales, and provides a state-of-the-art phylogenetic framework for further research.

The BAHD acyltransferase family is one of the largest enzyme families in flowering plants, containing dozens to hundreds of genes in individual genomes. Highly prevalent in angiosperm genomes, members of this family contribute to several pathways in primary and specialized metabolism. In this study, we performed a phylogenomic analysis of the family using 52 genomes across the plant kingdom to gain deeper insights into its functional evolution and enable function prediction. We found that BAHD expansion in land plants was associated with significant changes in various gene features. Using pre-defined BAHD clades, we identified clade expansions in different plant groups. In some groups, these expansions coincided with the prominence of metabolite classes such as anthocyanins (flowering plants) and hydroxycinnamic acid amides (monocots). Clade-wise motif-enrichment analysis revealed that some clades have novel motifs fixed on either the acceptor or the donor side, potentially reflecting historical routes of functional evolution. Co-expression analysis in rice and Arabidopsis further identified BAHDs with similar expression patterns, however, most co-expressed BAHDs belonged to different clades. Comparing BAHD paralogs, we found that gene expression diverges rapidly after duplication, suggesting that sub/neo-functionalization of duplicate genes occurs quickly via expression diversification. Analyzing co-expression patterns in Arabidopsis in conjunction with orthology-based substrate class predictions and metabolic pathway models led to the recovery of metabolic processes of most of the already-characterized BAHDs as well as definition of novel functional predictions for some uncharacterized BAHDs. Overall, this study provides new insights into the evolution of BAHD acyltransferases and sets up a foundation for their functional characterization.

Convolvulaceae, the morning glories or bindweeds, is a large family containing species of economic value, including crops, traditional medicines, ornamentals, and vegetables. However, not only are the phylogenetic relationships within this group still debated at the intertribal and intergeneric levels, but also plastid genome (plastome) complexity within Convolvulaceae is not well surveyed. We gathered 78 plastomes representing 17 genera across nine of the 12 Convolvulaceae tribes. Our plastid phylogenomic trees confirm the monophyly of Convolvulaceae, place the genus Jacquemontia within the subfamily Dicranostyloideae, and suggest that the tribe Merremieae is paraphyletic. In contrast, positions of the two genera Cuscuta and Erycibe are uncertain as the bootstrap support of the branches leading to them is moderate to weak. We show that nucleotide substitution rates are extremely variable among Convolvulaceae taxa and likely responsible for the topological uncertainty. Numerous plastomic rearrangements are detected in Convolvulaceae, including inversions, duplications, contraction and expansion of inverted repeats (IRs), and losses of genes and introns. Moreover, integrated foreign DNA of mitochondrial origin was found in the Jacquemontia plastome, adding a rare example of gene transfer from mitochondria to plastids in angiosperms. In the IR of Dichondra, we discovered an extra copy of rpl16 containing a direct repeat of ca. 200 bp long. This repeat was experimentally demonstrated to trigger effective homologous recombination, resulting in the coexistence of introncontaining and -lacking rpl16 duplicates. Therefore, we propose a hypothetical model to interpret intron loss accompanied by invasion of direct repeats at appropriate positions. Our model complements the intron loss model driven by retroprocessing when genes have lost introns but contain abundant RNA editing sites adjacent to former splicing sites. [ABSTRACT FROM AUTHOR]

The genus Sabia is a woody climber belonging to the family Sabiaceae, order Proteales. Several species of this genus have been utilized as medicines for treating diseases, such as rheumatic arthritis, traumatism, hepatitis, etc. However, the lack of molecular data has prevented the accurate identification and refinement of taxonomic relationships in this genus. In this study, chloroplast genomes of 11 samples of the genus Sabia were assembled and analyzed. These chloroplast genomes showed a typical quadripartite structure and ranged in length from 160,956 to 162,209 bp. The structure of the genomes was found to be relatively conserved, with 130 genes annotated, including 85 coding genes, 37 tRNA genes, and eight rRNA genes. A total of 78– 98 simple sequence repeats and 52–61 interspersed repeats were detected. Sequence alignment revealed 11 highly variable loci in chloroplast genomes. Among these loci, ndhF-ndhD achieved a remarkably higher resolution than the other regions. In addition, phylogenetic analysis indicated that Sect. Pachydiscus and Sect. Sabia of Sabia did not form two separate monophyletic groups. The divergence time calculated based on the Reltime method indicated that the evolutionary branches of Sabia and Meliosma started to form approximately 85.95 million years ago (Mya), and the species within Sabia began to diverge approximately 7.65 Mya. In conclusion, our study provides a basis for comprehensively exploring the phylogenetic relationships of Sabia. It also provides a methodological basis and data support for establishing a standardized and scientific identification system for this genus. [ABSTRACT FROM AUTHOR]

We assess relationships among 192 species in all 12 monocot orders and 72 of 77 families, using 602 conserved single-copy (CSC) genes and 1375 benchmarking single-copy ortholog (BUSCO) genes extracted from genomic and transcriptomic datasets. Phylogenomic inferences based on these data, using both coalescentbased and supermatrix analyses, are largely congruent with the most comprehensive plastome-based analysis, and nuclear-gene phylogenomic analyses with less comprehensive taxon sampling. The strongest discordance between the plastome and nuclear gene analyses is the monophyly of a clade comprising Asparagales and Liliales in our nuclear gene analyses, versus the placement of Asparagales and Liliales as successive sister clades to the commelinids in the plastome tree. Within orders, around six of 72 families shifted supported inferred relationships in the plastome-based tree. In Poales, the nuclear data place a clade comprising Ecdeiocoleaceae+Joinvilleaceae as sister to the grasses (Poaceae); Typhaceae, (rather thanBromeliaceae)are resolved as sister to all other Poales. In Commelinales, nuclear data place Philydraceae sister to all other families rather than to a clade comprising Haemodoraceae+Pontederiaceae as seen in the plastome tree. In Liliales, nuclear data place Liliaceae sister to Smilacaceae, and Melanthiaceae are placed sister to all other Liliales except Campynemataceae. Finally, in Alismatales, nuclear data strongly place Tofieldiaceae, rather than Araceae, as sister to all the other families, providing an alternative resolution of what has been themost problematic node to resolve using plastid data, outside of those involving achlorophyllousmycoheterotrophs. As seen in numerous prior studies, the placement of orders Acorales and Alismatales as successive sister lineages to allother extantmonocots. Only 21.2% of BUSCOgenes were demonstrably single-copy, yet phylogenomic inferences based on BUSCO andCSCgenesdid notdiffer, and overall functional annotationsof thetwosetswere very similar. Our analyses also reveal significant gene tree-species tree discordance despite high support values, as expected given incomplete lineage sorting (ILS) related to rapid diversification. Our study advances understanding of monocot relationships and the robustness of phylogenetic inferences based on large numbers of nuclear single-copy genes that can be obtained from transcriptomes and genomes. [ABSTRACT FROM AUTHOR]

Clematis is one of the largest genera of Ranunculaceae with many phylogenetic problems left to be resolved. Clematis species have considerable genome size of more than 7 Gbp, and there was no whole-genome reference sequence published in this genus. This raises difficulties in acquiring nuclear genome data for its phylogenetic analysis. Previous studies based on Sanger sequencing data, plastid genome data, and nrDNA sequences did not well resolve the phylogeny of Clematis. In this study, we used genome skimming and transcriptome data to assemble the plastid genome sequences, nuclear single nucleotide polymorphisms (SNPs) datasets, and single-copy nuclear orthologous genes (SCOGs) to reconstruct the phylogenetic backbone of Clematis, and test effectiveness of these genome partitioning methods. We also further analyzed the discordance among nuclear gene trees and between plastid and nuclear phylogenies. The results showed that the SCOGs datasets, assembled from transcriptome method, well resolved the phylogenetic backbone of Clematis. The nuclear SNPs datasets from genome skimming method can also produce similar results with the SCOGs data. In contrast to the plastid phylogeny, the phylogeny resolved by nuclear genome data is more robust and better corresponds to morphological characters. Our results suggested that rapid species radiation may have generated high level of incomplete lineage sorting, which was the major cause of nuclear gene discordance. Our simulation also showed that there may have been frequent interspecific hybridization events, which led to some of the cyto-nuclear discordances in Clematis. This study not only provides the first robust phylogenetic backbone of Clematis based on nuclear genome data, but also provides suggestions of genome partitioning strategies for the phylogenomic study of other plant taxa. [ABSTRACT FROM AUTHOR]

The core particle represents the catalytic portions of the 26S proteasomal complex. The genes encoding a- and b-subunits play a crucial role in protecting plants against various environmental stresses by controlling the quality of newly produced proteins. The 20S proteasome gene family has already been reported in model plants such as Arabidopsis and rice; however, they have not been studied in oilseed crops such as rapeseed (Brassica napus L.). In the present study, we identified 20S proteasome genes for a- (PA) and bsubunits (PB) in B. napus through systematically performed gene structure analysis, chromosomal location, conserved motif, phylogenetic relationship, and expression patterns. A total of 82 genes, comprising 35 BnPA and 47 BnPB of the 20S proteasome, were revealed in the B. napus genome. These genes were distributed on all 20 chromosomes of B. napus and most of these genes were duplicated on homoeologous chromosomes. The BnPA (a1-7) and BnPB (b1-7) genes were phylogenetically placed into seven clades. The pattern of expression of all the BnPA and BnPB genes was also studied using RNA-seq datasets under biotic and abiotic stress conditions. Out of 82 BnPA/PB genes, three exhibited high expression under abiotic stresses, whereas two genes were overexpressed in response to biotic stresses at both the seedling and flowering stages. Moreover, an additional eighteen genes were expressed under normal conditions. Overall, the current findings developed our understanding of the organization of the 20S proteasome genes in B. napus, and provided specific BnPA/PB genes for further functional research in response to abiotic and biotic stresses. [ABSTRACT FROM AUTHOR]

Sexually deceptive plants achieve pollination by enticing specific male insects as pollinators using a combination of olfactory, visual, and morphological mimicry. The sexually deceptive orchid genus Chiloglottis is comprised of some 30 species with predominantly dull green-red flowers except for the dark insectiform calli/callus structure from the labellum lamina. This unique structure mimics the female of the pollinator and potentially enhances the visibility of the mimic. However, the chemical and genetic basis for the color of these structures remains poorly understood across the genus. The goal of this study was to investigate the flower color biochemistry and patterns of gene expression across the anthocyanin and flavonol glycoside biosynthetic pathway within the calli structures across the three distinct clades of Chiloglottis (Formicifera, Reflexa, and Valida) using chemical and transcriptome analysis. Our phylogenomic analysis confirmed the close sister relationship between the Reflexa/Formicifera clades and reaffirms the basal position of the Valida clade. Additionally, the biochemical basis of the dark calli/callus structures is conserved across the genus. Nonetheless, the proportion of methoxylated anthocyanin and flavonol glycoside derivatives and the mean gene expression levels appear to differentiate the Reflexa and Formicifera clades from the Valida clade. In future studies, it will be of interest to tease apart the role of phylogeny, environment, pollinators, and other factors as potential drivers of the observed biochemistry and gene expression differences. It will also be important to characterize the function of candidate genes such as DFR, LDOX, and FLS in this fascinating case of flower color mimicry. [ABSTRACT FROM AUTHOR]

We assess relationships among 192 species in all 12 monocot orders and 72 of 77 families, using 602 conserved single-copy (CSC) genes and 1375 benchmarking single-copy ortholog (BUSCO) genes extracted from genomic and transcriptomic datasets. Phylogenomic inferences based on these data, using both coalescent-based and supermatrix analyses, are largely congruent with the most comprehensive plastome-based analysis, and nuclear-gene phylogenomic analyses with less comprehensive taxon sampling. The strongest discordance between the plastome and nuclear gene analyses is the monophyly of a clade comprising Asparagales and Liliales in our nuclear gene analyses, versus the placement of Asparagales and Liliales as successive sister clades to the commelinids in the plastome tree. Within orders, around six of 72 families shifted positions relative to the recent plastome analysis, but four of these involve poorly supported inferred relationships in the plastome-based tree. In Poales, the nuclear data place a clade comprising Ecdeiocoleaceae+Joinvilleaceae as sister to the grasses (Poaceae); Typhaceae, (rather than Bromeliaceae) are resolved as sister to all other Poales. In Commelinales, nuclear data place Philydraceae sister to all other families rather than to a clade comprising Haemodoraceae+Pontederiaceae as seen in the plastome tree. In Liliales, nuclear data place Liliaceae sister to Smilacaceae, and Melanthiaceae are placed sister to all other Liliales except Campynemataceae. Finally, in Alismatales, nuclear data strongly place Tofieldiaceae, rather than Araceae, as sister to all the other families, providing an alternative resolution of what has been the most problematic node to resolve using plastid data, outside of those involving achlorophyllous mycoheterotrophs. As seen in numerous prior studies, the placement of orders Acorales and Alismatales as successive sister lineages to all other extant monocots. Only 21.2% of BUSCO genes were demonstrably single-copy, yet phylogenomic inferences based on BUSCO and CSC genes did not differ, and overall functional annotations of the two sets were very similar. Our analyses also reveal significant gene tree-species tree discordance despite high support values, as expected given incomplete lineage sorting (ILS) related to rapid diversification. Our study advances understanding of monocot relationships and the robustness of phylogenetic inferences based on large numbers of nuclear single-copy genes that can be obtained from transcriptomes and genomes.

The core particle represents the catalytic portions of the 26S proteasomal complex. The genes encoding α- and β-subunits play a crucial role in protecting plants against various environmental stresses by controlling the quality of newly produced proteins. The 20S proteasome gene family has already been reported in model plants such as Arabidopsis and rice; however, they have not been studied in oilseed crops such as rapeseed (Brassica napus L.). In the present study, we identified 20S proteasome genes for α- (PA) and β-subunits (PB) in B. napus through systematically performed gene structure analysis, chromosomal location, conserved motif, phylogenetic relationship, and expression patterns. A total of 82 genes, comprising 35 BnPA and 47 BnPB of the 20S proteasome, were revealed in the B. napus genome. These genes were distributed on all 20 chromosomes of B. napus and most of these genes were duplicated on homoeologous chromosomes. The BnPA (α1-7) and BnPB (β1-7) genes were phylogenetically placed into seven clades. The pattern of expression of all the BnPA and BnPB genes was also studied using RNA-seq datasets under biotic and abiotic stress conditions. Out of 82 BnPA/PB genes, three exhibited high expression under abiotic stresses, whereas two genes were overexpressed in response to biotic stresses at both the seedling and flowering stages. Moreover, an additional eighteen genes were expressed under normal conditions. Overall, the current findings developed our understanding of the organization of the 20S proteasome genes in B. napus, and provided specific BnPA/PB genes for further functional research in response to abiotic and biotic stresses.

Sexually deceptive plants achieve pollination by enticing specific male insects as pollinators using a combination of olfactory, visual, and morphological mimicry. The sexually deceptive orchid genus Chiloglottis is comprised of some 30 species with predominantly dull green-red flowers except for the dark insectiform calli/callus structure from the labellum lamina. This unique structure mimics the female of the pollinator and potentially enhances the visibility of the mimic. However, the chemical and genetic basis for the color of these structures remains poorly understood across the genus. The goal of this study was to investigate the flower color biochemistry and patterns of gene expression across the anthocyanin and flavonol glycoside biosynthetic pathway within the calli structures across the three distinct clades of Chiloglottis (Formicifera, Reflexa, and Valida) using chemical and transcriptome analysis. Our phylogenomic analysis confirmed the close sister relationship between the Reflexa/Formicifera clades and reaffirms the basal position of the Valida clade. Additionally, the biochemical basis of the dark calli/callus structures is conserved across the genus. Nonetheless, the proportion of methoxylated anthocyanin and flavonol glycoside derivatives and the mean gene expression levels appear to differentiate the Reflexa and Formicifera clades from the Valida clade. In future studies, it will be of interest to tease apart the role of phylogeny, environment, pollinators, and other factors as potential drivers of the observed biochemistry and gene expression differences. It will also be important to characterize the function of candidate genes such as DFR, LDOX, and FLS in this fascinating case of flower color mimicry.

Subfamily Rubioideae is the largest of the main lineages in the coffee family (Rubiaceae), with over 8,000 species and 29 tribes. Phylogenetic relationships among tribes and other major clades within this group of plants are still only partly resolved despite considerable efforts. While previous studies have mainly utilized data from the organellar genomes and nuclear ribosomal DNA, we here use a large number of low-copy nuclear genes obtained via a target capture approach to infer phylogenetic relationships within Rubioideae. We included 101 Rubioideae species representing all but two (the monogeneric tribes Foonchewieae and Aitchinsonieae) of the currently recognized tribes, and all but one non-monogeneric tribe were represented by more than one genus. Using data from the 353 genes targeted with the universal Angiosperms353 probe set we investigated the impact of data type, analytical approach, and potential paralogs on phylogenetic reconstruction. We inferred a robust phylogenetic hypothesis of Rubioideae with the vast majority (or all) nodes being highly supported across all analyses and datasets and few incongruences between the inferred topologies. The results were similar to those of previous studies but novel relationships were also identified. We found that supercontigs [coding sequence (CDS) C non-coding sequence] clearly outperformed CDS data in levels of support and gene tree congruence. The full datasets (353 genes) outperformed the datasets with potentially paralogous genes removed (186 genes) in levels of support but increased gene tree incongruence slightly. The pattern of gene tree conflict at short internal branches were often consistent with high levels of incomplete lineage sorting (ILS) due to rapid speciation in the group. While concatenationand coalescence-based trees mainly agreed, the observed phylogenetic discordance between the two approaches may be best explained by their differences in accounting for ILS. The use of target capture data. greatly improved our confidence and understanding of the Rubioideae phylogeny, highlighted by the increased support for previously uncertain relationships and the increased possibility to explore sources of underlying phylogenetic discordance. [ABSTRACT FROM AUTHOR]

Pseudomonas sp. strain 1008 was isolated from the rhizosphere of field grown wheat plants at the tillering stage in an agricultural plot near Pergamino city, Argentina. Based on its in vitro phosphate solubilizing capacity and the production of IAA, strain 1008 was formulated as an inoculant for bacterization of wheat seeds and subjected to multiple field assays within the period 2010-2017. Pseudomonas sp. strain 1008 showed a robust positive impact on the grain yield (C8% on average) across a number of campaigns, soil properties, seed genotypes, and with no significant influence of the simultaneous seed treatment with a fungicide, strongly supporting the use of this biostimulant bacterium as an agricultural input for promoting the yield of wheat. Full genome sequencing revealed that strain 1008 has the capacity to access a number of sources of inorganic and organic phosphorus, to compete for iron scavenging, to produce auxin, 2,3-butanediol and acetoin, and to metabolize GABA. Additionally, the genome of strain 1008 harbors several loci related to rhizosphere competitiveness, but it is devoid of biosynthetic gene clusters for production of typical secondary metabolites of biocontrol representatives of the Pseudomonas genus. Finally, the phylogenomic, phenotypic, and chemotaxonomic comparative analysis of strain 1008 with related taxa strongly suggests that this wheat rhizospheric biostimulant isolate is a representative of a novel species within the genus Pseudomonas, for which the name Pseudomonas pergaminensis sp. nov. (type strain 1008T = DSM 113453T = ATCC TSD-287T) is proposed. [ABSTRACT FROM AUTHOR]

Subfamily Nandinoideae Heintze (Berberidaceae), comprising four genera and ca. 19 species, is disjunctively distributed in eastern North America vs. Eurasia (eastern Asia, Central Asia, Middle East, and southeastern Europe), and represents an ideal taxon to explore plastid phylogenomics and plastome evolution in Berberidaceae. Many species of this subfamily have been listed as national or international rare and endangered plants. In this study, we sequenced and assembled 20 complete plastomes, representing three genera and 13 species of Nandinoideae. Together with six plastomes from GenBank, a total of 26 plastomes, representing all four genera and 16 species of Nandinoideae, were used for comparative genomic and phylogenomic analyses. These plastomes showed significant differences in overall size (156,626-161,406 bp), which is mainly due to the expansion in inverted repeat (IR) regions and/or insertion/deletion (indel) events in intergenic spacer (IGS) regions. A 75-bp deletion in the ndhF gene occurred in Leontice and Gymnospermium when compared with Nandina and Caulophyllum. We found a severe truncation at the 5' end of ycf 1 in three G. altaicum plastomes, and a premature termination of ropC1 in G. microrrhynchum. Our phylogenomic results support the topology of {Nandina, [Caulophyllum, (Leontice, Gymnospermium)]}. Within the core genus Gymnospermium, we identified G. microrrhynchum from northeastern Asia (Clade A) as the earliest diverging species, followed by G. kiangnanense from eastern China (Clade B), while the rest species clustered into the two sister clades (C and D). Clade C included three species from West Tianshan (G. albertii, G. darwasicum, G. vitellinum). Clade D consisted of G. altaicum from northern Central Asia, plus one species from the Caucasus Mountains (G. smirnovii) and three from southeastern Europe (G. odessanum, G. peloponnesiacum, G. scipetarum). Overall, we identified 21 highly variable plastome regions, including two coding genes (rpl22, ycf1) and 19 intergenic spacer (IGS) regions, all with nucleotide diversity (Pi) values > 0.02. These molecular markers should serve as powerful tools (including DNA barcodes) for future phylogenetic, phylogeographic and conservation genetic studies. [ABSTRACT FROM AUTHOR]

Subfamily Rubioideae is the largest of the main lineages in the coffee family (Rubiaceae), with over 8,000 species and 29 tribes. Phylogenetic relationships among tribes and other major clades within this group of plants are still only partly resolved despite considerable efforts. While previous studies have mainly utilized data from the organellar genomes and nuclear ribosomal DNA, we here use a large number of low-copy nuclear genes obtained via a target capture approach to infer phylogenetic relationships within Rubioideae. We included 101 Rubioideae species representing all but two (the monogeneric tribes Foonchewieae and Aitchinsonieae) of the currently recognized tribes, and all but one non-monogeneric tribe were represented by more than one genus. Using data from the 353 genes targeted with the universal Angiosperms353 probe set we investigated the impact of data type, analytical approach, and potential paralogs on phylogenetic reconstruction. We inferred a robust phylogenetic hypothesis of Rubioideae with the vast majority (or all) nodes being highly supported across all analyses and datasets and few incongruences between the inferred topologies. The results were similar to those of previous studies but novel relationships were also identified. We found that supercontigs [coding sequence (CDS) + non-coding sequence] clearly outperformed CDS data in levels of support and gene tree congruence. The full datasets (353 genes) outperformed the datasets with potentially paralogous genes removed (186 genes) in levels of support but increased gene tree incongruence slightly. The pattern of gene tree conflict at short internal branches were often consistent with high levels of incomplete lineage sorting (ILS) due to rapid speciation in the group. While concatenation- and coalescence-based trees mainly agreed, the observed phylogenetic discordance between the two approaches may be best explained by their differences in accounting for ILS. The use of target capture data greatly improved our confidence and understanding of the Rubioideae phylogeny, highlighted by the increased support for previously uncertain relationships and the increased possibility to explore sources of underlying phylogenetic discordance.

Subfamily Nandinoideae Heintze (Berberidaceae), comprising four genera and ca. 19 species, is disjunctively distributed in eastern North America vs. Eurasia (eastern Asia, Central Asia, Middle East, and southeastern Europe), and represents an ideal taxon to explore plastid phylogenomics and plastome evolution in Berberidaceae. Many species of this subfamily have been listed as national or international rare and endangered plants. In this study, we sequenced and assembled 20 complete plastomes, representing three genera and 13 species of Nandinoideae. Together with six plastomes from GenBank, a total of 26 plastomes, representing all four genera and 16 species of Nandinoideae, were used for comparative genomic and phylogenomic analyses. These plastomes showed significant differences in overall size (156,626–161,406 bp), which is mainly due to the expansion in inverted repeat (IR) regions and/or insertion/deletion (indel) events in intergenic spacer (IGS) regions. A 75-bp deletion in the ndh F gene occurred in Leontice and Gymnospermium when compared with Nandina and Caulophyllum. We found a severe truncation at the 5' end of ycf 1 in three G. altaicum plastomes, and a premature termination of ropC 1 in G. microrrhynchum. Our phylogenomic results support the topology of { Nandina , [ Caulophyllum , (Leontice , Gymnospermium)]}. Within the core genus Gymnospermium , we identified G. microrrhynchum from northeastern Asia (Clade A) as the earliest diverging species, followed by G. kiangnanense from eastern China (Clade B), while the rest species clustered into the two sister clades (C and D). Clade C included three species from West Tianshan (G. albertii , G. darwasicum , G. vitellinum). Clade D consisted of G. altaicum from northern Central Asia, plus one species from the Caucasus Mountains (G. smirnovii) and three from southeastern Europe (G. odessanum , G. peloponnesiacum , G. scipetarum). Overall, we identified 21 highly variable plastome regions, including two coding genes (rpl 22, ycf 1) and 19 intergenic spacer (IGS) regions, all with nucleotide diversity (Pi) values > 0.02. These molecular markers should serve as powerful tools (including DNA barcodes) for future phylogenetic, phylogeographic and conservation genetic studies. [ABSTRACT FROM AUTHOR]

The topographic gradients of the Tropical Andes may have triggered species divergence by different mechanisms. Topography separates species' geographical ranges and offers climatic heterogeneity, which could potentially foster local adaptation to specific climatic conditions and result in narrowly distributed endemic species. Such a pattern is found in the Andean centered palm genus Aiphanes. To test the extent to which geographic barriers and climatic heterogeneity can explain distribution patterns in Aiphanes , we sampled 34 out of 36 currently recognized species in that genus and sequenced them by Sanger sequencing and/or sequence target capture sequencing. We generated Bayesian, likelihood, and species-tree phylogenies, with which we explored climatic trait evolution from current climatic occupation. We also estimated species distribution models to test the relative roles of geographical and climatic divergence in their evolution. We found that Aiphanes originated in the Miocene in Andean environments and possibly in mid-elevation habitats. Diversification is related to the occupation of the adjacent high and low elevation habitats tracking high annual precipitation and low precipitation seasonality (moist habitats). Different species in different clades repeatedly occupy all the different temperatures offered by the elevation gradient from 0 to 3,000 m in different geographically isolated areas. A pattern of conserved adaptation to moist environments is consistent among the clades. Our results stress the evolutionary roles of niche truncation of wide thermal tolerance by physical range fragmentation, coupled with water-related niche conservatism, to colonize the topographic gradient. [ABSTRACT FROM AUTHOR]

The species in the genus Cypripedium (Orchidaceae) are considered endangered, mainly distributed in the temperate regions of the Northern Hemisphere, with high ornamental and economic value. Despite previous extensive studies based on both morphology and molecular data, species and sections relationships within Cypripedium remain controversial. Here, we employed two newly generated Cypripedium chloroplast genomes with five other published genomes to elucidate their genomic characteristics. The two genomes were 162,773–207,142 bp in length and contained 128–130 genes, including 82–84 protein-coding genes, 38 tRNA genes, and 8 rRNA genes. We identified 2,192 simple sequence repeats, 786 large repeat sequences, and 7,929 variable loci. The increase of repeat sequences (simple sequence repeats and large repeat sequences) causes a significant amplification in the chloroplast genome size of Cypripedium. The expansion of the IR region led to the pseudogenization or loss of genes in the SSC region. In addition, we identified 12 highly polymorphic loci (Pi > 0.09) suitable for inferring the phylogeny of Cypripedium species. Based on data sets of whole chloroplast genomes (IRa excluded) and protein-coding sequences, a well-supported phylogenetic tree was reconstructed, strongly supporting the five subfamilies of Orchidaceae and the genus Cypripedium as monophyletic taxa. Our findings also supported that C. palangshanense belonged to sect. Palangshanensia rather than sect. Retinervia. This study also enriched the genomic resources of Cypripedium , which may help to promote the conservation efforts of these endangered species. [ABSTRACT FROM AUTHOR]

Phytolacca is the largest genus of Phytolaccaceae. Owing to interspecific hybridization, infraspecific variation, and apparent weak genetic control of many qualitative characters, which have obscured boundaries between species, the classification and phylogenetic relationships of this genus are unclear. Native Phytolacca is disjunctly distributed in America, eastern Asia, and Africa, and the biogeographic history of the genus remained unresolved. In this study, we used the whole chloroplast genome and three markers (nrDNA, rbcL , and matK) to reconstruct phylogenetic relationships within Phytolacca , analyze divergence times, and infer biogeographic histories. The phylogenetic results indicate that Phytolacca is monophyletic, which is inconsistent with the infrageneric classification based on morphology. According to the divergence time estimation, Phytolacca began to diversify at approximately 20.30 Ma during the early Miocene. Central America, including Mexico, Costa Rica, and Colombia, is the center of species diversity. Biogeographical analysis indicated five main dispersal events and Phytolacca originated from Central and South America. Birds may be the primary agents of dispersal because of the fleshy fruiting of Phytolacca. This study extended sampling and added more genetic characteristics to infer the evolutionary history of Phytolacca , providing new insights for resolving the classification and elucidating the dispersal events of Phytolacca. [ABSTRACT FROM AUTHOR]

Arnebiae Radix is a traditional medicine with pleiotropic properties that has been used for several 100 years. There are five species of Arnebia in China, and the two species Arnebia euchroma and Arnebia guttata are the source plants of Arnebiae Radix according to the Chinese Pharmacopoeia. Molecular markers that permit species identification and facilitate studies of the genetic diversity and divergence of the wild populations of these two source plants have not yet been developed. Here, we sequenced the chloroplast genomes of 56 samples of five Arnebia species using genome skimming methods. The Arnebia chloroplast genomes exhibited quadripartite structures with lengths from 149,539 and 152,040 bp. Three variable markers (rps16-trnQ , ndhF-rpl32 , and ycf1b) were identified, and these markers exhibited more variable sites than universal chloroplast markers. The phylogenetic relationships among the five Arnebia species were completely resolved using the whole chloroplast genome sequences. Arnebia arose during the Oligocene and diversified in the middle Miocene; this coincided with two geological events during the late Oligocene and early Miocene: warming and the progressive uplift of Tianshan and the Himalayas. Our analyses revealed that A. euchroma and A. guttata have high levels of genetic diversity and comprise two and three subclades, respectively. The two clades of A. euchroma exhibited significant genetic differences and diverged at 10.18 Ma in the middle Miocene. Three clades of A. guttata diverged in the Pleistocene. The results provided new insight into evolutionary history of Arnebia species and promoted the conservation and exploitation of A. euchroma and A. guttata. [ABSTRACT FROM AUTHOR]

The Orchidaceae is rivaled only by the Asteraceae as the largest plant family, with the estimated number of species exceeding 25,000 and encompassing more than 700 genera. To gain insights into the mechanisms driving species diversity across both global and local scales, well-supported phylogenies targeting different taxonomic groups and/or geographical regions will be crucial. High-throughput sequencing technologies have revolutionized the field of molecular phylogenetics by simplifying the process of obtaining genome-scale sequence data. Consequently, there has been an explosive growth of such data in public repositories. Here we took advantage of this unprecedented access to transcriptome data from predominantly non-phylogenetic studies to assess if it can be repurposed to gain rapid and accurate phylogenetic insights across the orchids. Exhaustive searches revealed transcriptomic data for more than 100 orchid species spanning 5 subfamilies, 13 tribes, 21 subtribes, and 50 genera that were amendable for exploratory phylotranscriptomic analysis. Next, we performed reassembly of the transcriptomes before strategic selection of the final samples based on a gene completeness evaluation. Drawing on these data, we report phylogenetic analyses at both deep and shallow evolutionary scales via maximum likelihood and shortcut coalescent species tree methods. In this perspective, we discuss some key outcomes of this study and conclude by highlighting other complementary, albeit rarely explored, insights beyond phylogenetic analysis that repurposed multi-tissue transcriptome can offer. [ABSTRACT FROM AUTHOR]

The genus Campylotropis Bunge (Desmodieae, Papilionoideae) comprises about 37 species distributed in temperate and tropical Asia. Despite the great potential in soil conservation, horticulture, and medicine usage, little is known about the evolutionary history and phylogenetic relationships of Campylotropis due to insufficient genetic resources. Here, we sequenced and assembled 21 complete chloroplast genomes of Campylotropis species. In combination with the previously published chloroplast genomes of C. macrocarpa and closely related species, we conducted comparative genomics and phylogenomic analysis on these data. Comparative analysis of the genome size, structure, expansion and contraction of inverted repeat (IR) boundaries, number of genes, GC content, and pattern of simple sequence repeats (SSRs) revealed high similarities among the Campylotropis chloroplast genomes. The activities of long sequence repeats contributed to the variation in genome size and gene content in Campylotropis chloroplast genomes. The Campylotropis chloroplast genomes showed moderate sequence variation, and 13 highly variable regions were identified for species identification and further phylogenetic studies. We also reported one more case of matK pseudogene in the legume family. The phylogenetic analysis confirmed the monophyly of Campylotropis and the sister relationship between Lespedeza and Kummerowia , the latter two genera were then sister to Campylotropis. The intrageneric relationships of Campylotropis based on genomic scale data were firstly reported in this study. The two positively selected genes (atpF and rps19) and eight fast-evolving genes identified in this study may help us to understand the adaptation of Campylotropis species. Overall, this study enhances our understanding of the chloroplast genome evolution and phylogenetic relationships of Campylotropis. [ABSTRACT FROM AUTHOR]

The Cellulose synthase superfamily synthesizes cellulose and different hemicellulosic polysaccharides in plant cell walls. While much has been discovered about the evolution and function of these genes, their genomic architecture and relationship with gene (sub-)functionalization and evolution remains unclear. By using 242 genomes covering plant evolution from green algae to eudicots, we performed a large-scale analysis of synteny, phylogenetic, and functional data of the CesA superfamily. Results revealed considerable gene copy number variation across species and gene families, and also two patterns – singletons vs. tandem arrays – in chromosomic gene arrangement. Synteny analysis revealed exceptional conservation of gene architecture across species, but also lineage-specific patterns across gene (sub-)families. Synteny patterns correlated with gene sub-functionalization into primary and secondary CesAs and distinct CslD functional isoforms. Furthermore, a genomic context shift of a group of cotton secondary CesAs was associated with peculiar properties of cotton fiber synthesis. Finally, phylogenetics suggested that primary CesA sequences appeared before the secondary CesAs , while phylogenomic analyses unveiled the genomic trace of the CslD duplication that initiated the CslF family. Our results describe in detail the genomic architecture of the CesA superfamily in plants, highlighting its crucial relevance for gene diversification and sub-functionalization, and for understanding their evolution. [ABSTRACT FROM AUTHOR]

The species in the genus Cypripedium (Orchidaceae) are considered endangered, mainly distributed in the temperate regions of the Northern Hemisphere, with high ornamental and economic value. Despite previous extensive studies based on both morphology and molecular data, species and sections relationships within Cypripedium remain controversial. Here, we employed two newly generated Cypripedium chloroplast genomes with five other published genomes to elucidate their genomic characteristics. The two genomes were 162,773–207,142 bp in length and contained 128–130 genes, including 82–84 protein-coding genes, 38 tRNA genes, and 8 rRNA genes. We identified 2,192 simple sequence repeats, 786 large repeat sequences, and 7,929 variable loci. The increase of repeat sequences (simple sequence repeats and large repeat sequences) causes a significant amplification in the chloroplast genome size of Cypripedium. The expansion of the IR region led to the pseudogenization or loss of genes in the SSC region. In addition, we identified 12 highly polymorphic loci (Pi > 0.09) suitable for inferring the phylogeny of Cypripedium species. Based on data sets of whole chloroplast genomes (IRa excluded) and protein-coding sequences, a well-supported phylogenetic tree was reconstructed, strongly supporting the five subfamilies of Orchidaceae and the genus Cypripedium as monophyletic taxa. Our findings also supported that C. palangshanense belonged to sect. Palangshanensia rather than sect. Retinervia. This study also enriched the genomic resources of Cypripedium, which may help to promote the conservation efforts of these endangered species.

Phytolacca is the largest genus of Phytolaccaceae. Owing to interspecific hybridization, infraspecific variation, and apparent weak genetic control of many qualitative characters, which have obscured boundaries between species, the classification and phylogenetic relationships of this genus are unclear. Native Phytolacca is disjunctly distributed in America, eastern Asia, and Africa, and the biogeographic history of the genus remained unresolved. In this study, we used the whole chloroplast genome and three markers (nrDNA, rbcL, and matK) to reconstruct phylogenetic relationships within Phytolacca, analyze divergence times, and infer biogeographic histories. The phylogenetic results indicate that Phytolacca is monophyletic, which is inconsistent with the infrageneric classification based on morphology. According to the divergence time estimation, Phytolacca began to diversify at approximately 20.30 Ma during the early Miocene. Central America, including Mexico, Costa Rica, and Colombia, is the center of species diversity. Biogeographical analysis indicated five main dispersal events and Phytolacca originated from Central and South America. Birds may be the primary agents of dispersal because of the fleshy fruiting of Phytolacca. This study extended sampling and added more genetic characteristics to infer the evolutionary history of Phytolacca, providing new insights for resolving the classification and elucidating the dispersal events of Phytolacca.

The genus Campylotropis Bunge (Desmodieae, Papilionoideae) comprises about 37 species distributed in temperate and tropical Asia. Despite the great potential in soil conservation, horticulture, and medicine usage, little is known about the evolutionary history and phylogenetic relationships of Campylotropis due to insufficient genetic resources. Here, we sequenced and assembled 21 complete chloroplast genomes of Campylotropis species. In combination with the previously published chloroplast genomes of C. macrocarpa and closely related species, we conducted comparative genomics and phylogenomic analysis on these data. Comparative analysis of the genome size, structure, expansion and contraction of inverted repeat (IR) boundaries, number of genes, GC content, and pattern of simple sequence repeats (SSRs) revealed high similarities among the Campylotropis chloroplast genomes. The activities of long sequence repeats contributed to the variation in genome size and gene content in Campylotropis chloroplast genomes. The Campylotropis chloroplast genomes showed moderate sequence variation, and 13 highly variable regions were identified for species identification and further phylogenetic studies. We also reported one more case of matK pseudogene in the legume family. The phylogenetic analysis confirmed the monophyly of Campylotropis and the sister relationship between Lespedeza and Kummerowia, the latter two genera were then sister to Campylotropis. The intrageneric relationships of Campylotropis based on genomic scale data were firstly reported in this study. The two positively selected genes (atpF and rps19) and eight fast-evolving genes identified in this study may help us to understand the adaptation of Campylotropis species. Overall, this study enhances our understanding of the chloroplast genome evolution and phylogenetic relationships of Campylotropis.

The Cellulose synthase superfamily synthesizes cellulose and different hemicellulosic polysaccharides in plant cell walls. While much has been discovered about the evolution and function of these genes, their genomic architecture and relationship with gene (sub-)functionalization and evolution remains unclear. By using 242 genomes covering plant evolution from green algae to eudicots, we performed a large-scale analysis of synteny, phylogenetic, and functional data of the CesA superfamily. Results revealed considerable gene copy number variation across species and gene families, and also two patterns – singletons vs. tandem arrays – in chromosomic gene arrangement. Synteny analysis revealed exceptional conservation of gene architecture across species, but also lineage-specific patterns across gene (sub-)families. Synteny patterns correlated with gene sub-functionalization into primary and secondary CesAs and distinct CslD functional isoforms. Furthermore, a genomic context shift of a group of cotton secondary CesAs was associated with peculiar properties of cotton fiber synthesis. Finally, phylogenetics suggested that primary CesA sequences appeared before the secondary CesAs, while phylogenomic analyses unveiled the genomic trace of the CslD duplication that initiated the CslF family. Our results describe in detail the genomic architecture of the CesA superfamily in plants, highlighting its crucial relevance for gene diversification and sub-functionalization, and for understanding their evolution.

Multidrug and toxic compound extrusion (MATE) transporters comprise a multigene family that mediates multiple functions in plants through the efflux of diverse substrates including organic molecules, specialized metabolites, hormones, and xenobiotics. MATE classification based on genome-wide studies remains ambiguous, likely due to a lack of large-scale phylogenomic studies and/or reference sequence datasets. To resolve this, we established a phylogeny of the plant MATE gene family using a comprehensive kingdom-wide phylogenomic analysis of 74 diverse plant species. We identified more than 4,000 MATEs, which were classified into 14 subgroups based on a systematic bioinformatics pipeline using USEARCH, blast+ and synteny network tools. Our classification was performed using a four-step process, whereby MATEs sharing ≥ 60% protein sequence identity with a ≤ 1E-05 threshold at different sequence lengths (either full-length, ≥ 60% length, or ≥ 150 amino acids) or retaining in the similar synteny blocks were assigned to the same subgroup. In this way, we assigned subgroups to 95.8% of the identified MATEs, which we substantiated using synteny network clustering analysis. The subgroups were clustered under four major phylogenetic groups and named according to their clockwise appearance within each group. We then generated a reference sequence dataset, the usefulness of which was demonstrated in the classification of MATEs in additional species not included in the original analysis. Approximately 74% of the plant MATEs exhibited synteny relationships with angiosperm-wide or lineage-, order/family-, and species-specific conservation. Most subgroups evolved independently, and their distinct evolutionary trends were likely associated with the development of functional novelties or the maintenance of conserved functions. Together with the systematic classification and synteny network profiling analyses, we identified all the major evolutionary events experienced by the MATE gene family in plants. We believe that our findings and the reference dataset provide a valuable resource to guide future functional studies aiming to explore the key roles of MATEs in different aspects of plant physiology. Our classification framework can also be readily extendable to other (super) families. [ABSTRACT FROM AUTHOR]