Incongruence between gene trees and species trees and phylogenetic signal variation in plastid genes

Plastid phylogenomics and cytonuclear discordance in Rubioideae, Rubiaceae

In this study of evolutionary relationships in the subfamily Rubioideae (Rubiaceae), we take advantage of the off-target proportion of reads generated via previous target capture sequencing projects based on nuclear genomic data to build a plastome phylogeny and investigate cytonuclear discordance. The assembly of off-target reads resulted in a comprehensive plastome dataset and robust inference of phylogenetic relationships, where most intratribal and intertribal relationships are resolved with strong support. While the phylogenetic results were mostly in agreement with previous studies based on plastome data, novel relationships in the plastid perspective were also detected. For example, our analyses of plastome data provide strong support for the SCOUT clade and its sister relationship to the remaining members of the subfamily, which differs from previous results based on plastid data but agrees with recent results based on nuclear genomic data. However, several instances of highly supported cytonuclear discordance were identified across the Rubioideae phylogeny. Coalescent simulation analysis indicates that while ILS could, by itself, explain the majority of the discordant relationships, plastome introgression may be the better explanation in some cases. Our study further indicates that plastomes across the Rubioideae are, with few exceptions, highly conserved and mainly conform to the structure, gene content, and gene order present in the majority of the flowering plants.

Phylogenomic conflict analyses of the plastid and mitochondrial genomes via deep genome skimming highlight their independent evolutionary histories: A case study in the cinquefoil genus Potentilla sensu lato (Potentilleae, Rosaceae)

Phylogenomic conflicts are widespread among genomic data, with most previous studies primarily focusing on nuclear datasets instead of organellar genomes. In this study, we investigate phylogenetic conflict analyses within and between plastid and mitochondrial genomes using Potentilla as a case study. We generated three plastid datasets (coding, noncoding, and all-region) and one mitochondrial dataset (coding regions) to infer phylogenies based on concatenated and multispecies coalescent (MSC) methods. Conflict analyses were then performed using PhyParts and Quartet Sampling (QS). Both plastid and mitochondrial genomes divided the Potentilla into eight highly supported clades, two of which were newly identified in this study. While most organellar loci were uninformative for the majority of nodes (bootstrap value < 70%), PhyParts and QS detected conflicting signals within the two organellar genomes. Regression analyses revealed that conflict signals mainly occurred among shorter loci, whereas longer loci tended to be more concordant with the species tree. In addition, two significant disagreements between the two organellar genomes were detected, likely attributed to hybridization and/or incomplete lineage sorting. Our results demonstrate that mitochondrial genes can fully resolve the phylogenetic relationships among eight major clades of Potentilla and are not always linked with plastome in evolutionary history. Stochastic inferences appear to be the primary source of observed conflicts among the gene trees. We recommend that the loci with short sequence length or containing limited informative sites should be used cautiously in MSC analysis, and suggest the joint application of concatenated and MSC methods for phylogenetic inference using organellar genomes.

Phylogenomics and topological conflicts in the tribe Anthospermeae (Rubiaceae)

Genome skimming (shallow whole-genome sequencing) offers time- and cost-efficient production of large amounts of DNA data that can be used to address unsolved evolutionary questions. Here we address phylogenetic relationships and topological incongruence in the tribe Anthospermeae (Rubiaceae), using phylogenomic data from the mitochondrion, the nuclear ribosomal cistron, and the plastome. All three genomic compartments resolve relationships in the Anthospermeae; the tribe is monophyletic and consists of three major subclades. Carpacoce Sond. is sister to the remaining clade, which comprises an African subclade and a Pacific subclade. Most results, from all three genomic compartments, are statistically well supported; however, not fully consistent. Intergenomic topological incongruence is most notable in the Pacific subclade but present also in the African subclade. Hybridization and introgression followed by organelle capture may explain these conflicts but other processes, such as incomplete lineage sorting (ILS), can yield similar patterns and cannot be ruled out based on the results. Whereas the null hypothesis of congruence among all sequenced loci in the individual genomes could not be rejected for nuclear and mitochondrial data, it was rejected for plastid data. Phylogenetic analyses of three subsets of plastid loci identified using the hierarchical likelihood ratio test demonstrated statistically supported intragenomic topological incongruence. Given that plastid genes are thought to be fully linked, this result is surprising and may suggest modeling or sampling error. However, biological processes such as biparental inheritance and inter-plastome recombination have been reported and may be responsible for the observed intragenomic incongruence. Mitochondrial insertions into the plastome are rarely documented in angiosperms. Our results indicate that a mitochondrial insertion event in the plastid trnS GGA - rps4 IGS region occurred in the common ancestor of the Pacific clade of Anthospermeae. Exclusion/inclusion of this locus in phylogenetic analyses had a strong impact on topological results in the Pacific clade.

Insights into the phylogenetic relationships and species boundaries of the Myricaria squamosa complex (Tamaricaceae) based on the complete chloroplast genome

Myricaria plants are widely distributed in Eurasia and are helpful for windbreak and embankment protection. Current molecular evidence has led to controversy regarding species boundaries within the Myricaria genus and interspecific phylogenetic relationships between three specific species-M. bracteata, M. paniculata and M. squamosa-which have remained unresolved. This study treated these three unresolved taxa as a species complex, named the M. squamosa complex. The genome skimming approach was used to determine 35 complete plastome sequences and nuclear ribosomal DNA sequences for the said complex and other closely related species, followed by de novo assembly. Comparative analyses were conducted across Myricaria to identify the genome size, gene content, repeat type and number, SSR (simple sequence repeat) abundance, and codon usage bias of chloroplast genomes. Tree-based species delimitation results indicated that M. bracteata, M. paniculata and M. squamosa could not be distinguished and formed two monophyletic lineages (P1 and P2) that were clustered together. Compared to plastome-based species delimitation, the standard nuclear DNA barcode had the lowest species resolution, and the standard chloroplast DNA barcode and group-specific barcodes delimitated a maximum of four out of the five species. Plastid phylogenomics analyses indicated that the monophyletic M. squamosa complex is comprised of two evolutionarily significant units: one in the western Tarim Basin and the other in the eastern Qinghai-Tibet Plateau. This finding contradicts previous species discrimination and promotes the urgent need for taxonomic revision of the threatened genus Myricaria. Dense sampling and plastid genomes will be essential in this effort. The super-barcodes and specific barcode candidates outlined in this study will aid in further studies of evolutionary history.

The complete chloroplast genome sequence of Zygophyllum brachypterum (Zygophyllaceae) reveals its distinctive characteristics and evolutionary implication

Zygophyllum brachypterum Karelin & Kirilov belongs to Zygophyllaceae and is mainly distributed in the desert regions of Central Asia, Mongolia, and Northwest China. The species is valuable in exploring the adaptations of Zygophyllaceae plants to salt stress in ecological environments. In this study, we report the complete chloroplast (cp) genome of Z. brachypterum. The entire cp genome was 104590 bp in length, with a large single-copy region (LSC, 79170 bp), a small single-copy region (SSC, 16778 bp), and two inverted repeats (IRa/IRb) of 4321 bp each. A total of 106 genes were detected, among which seven were located in the IRs, and 65, 30, and 4 were protein-coding, tRNA, and rRNA genes, respectively. Notably, eleven genes encoding the subunits of NAD(P)H dehydrogenase complex (NDH) were absent. Phylogenetic analysis indicated that Z. brachypterum belonged to Zygophylloideae (Zygophyllaceae). Furthermore, it was closely related to Z. fabago and Z. kansuense.

More than a spiny morphology: plastome variation in the prickly pear cacti (Opuntieae)

BACKGROUND

Plastid genomes (plastomes) have long been recognized as highly conserved in their overall structure, size, gene arrangement and content among land plants. However, recent studies have shown that some lineages present unusual variations in some of these features. Members of the cactus family are one of these lineages, with distinct plastome structures reported across disparate lineages, including gene losses, inversions, boundary movements or loss of the canonical inverted repeat (IR) region. However, only a small fraction of cactus diversity has been analysed so far.

METHODS

Here, we investigated plastome features of the tribe Opuntieae, the remarkable prickly pear cacti, which represent one of the most diverse and important lineages of Cactaceae. We assembled de novo the plastome of 43 species, representing a comprehensive sampling of the tribe, including all seven genera, and analysed their evolution in a phylogenetic comparative framework. Phylogenomic analyses with different datasets (full plastome sequences and genes only) were performed, followed by congruence analyses to assess signals underlying contentious nodes.

KEY RESULTS

Plastomes varied considerably in length, from 121 to 162 kbp, with striking differences in the content and size of the IR region (contraction and expansion events), including a lack of the canonical IR in some lineages and the pseudogenization or loss of some genes. Overall, nine different types of plastomes were reported, deviating in the presence of the IR region or the genes contained in the IR. Overall, plastome sequences resolved phylogenetic relationships within major clades of Opuntieae with high bootstrap values but presented some contentious nodes depending on the dataset analysed (e.g. whole plastome vs. genes only). Congruence analyses revealed that most plastidial regions lack phylogenetic resolution, while few markers are supporting the most likely topology. Likewise, alternative topologies are driven by a handful of plastome markers, suggesting recalcitrant nodes in the phylogeny.

CONCLUSIONS

Our study reveals a dynamic nature of plastome evolution across closely related lineages, shedding light on peculiar features of plastomes. Variation of plastome types across Opuntieae is remarkable in size, structure and content and can be important for the recognition of species in some major clades. Unravelling connections between the causes of plastome variation and the consequences for species biology, physiology, ecology, diversification and adaptation is a promising and ambitious endeavour in cactus research. Although plastome data resolved major phylogenetic relationships, the generation of nuclear genomic data is necessary to confront these hypotheses and assess the recalcitrant nodes further.

On Pattern-Cladistic Analyses Based on Complete Plastid Genome Sequences

The fundamental Hennigian principle, grouping solely on synapomorphy, is seldom used in modern phylogenetics. In the submitted paper, we apply this principle in reanalyzing five datasets comprising 197 complete plastid genomes (plastomes). We focused on the latter because plastome-based DNA sequence data gained dramatic popularity in molecular systematics during the last decade. We show that pattern-cladistic analyses based on complete plastid genome sequences can successfully resolve affinities between plant taxa, simultaneously simplifying both the genomic and analytical frameworks of phylogenetic studies. We developed "Matrix to Newick" (M2N), a program to represent the standard molecular alignment of plastid genomes in the form of trees or relationships directly. Thus, massive plastome-based DNA sequence data can be successfully represented in a relational form rather than as a standard molecular alignment. Application of methods of median supertree construction (the Average Consensus method has been used as an example in this study) or Maximum Parsimony analysis to relational representations of plastome sequence data may help systematist to avoid the complicated assumption-based frameworks of Maximum Likelihood or Bayesian phylogenetics that are most used today in massive plastid sequence data analyses. We also found that significant amounts of pure genomic information that typically accommodate the majority of current plastid phylogenomic studies can be effectively dropped by systematists if they focus on the pattern-cladistics or relational analyses of plastome-based molecular data. The proposed pattern-cladistic approach is a powerful and straightforward heuristic alternative to modern plastome-based phylogenetics.

Structural divergence and phylogenetic relationships of Ajania (Asteraceae) from plastomes and ETS

BACKGROUND

Ajania Poljakov, an Asteraceae family member, grows mostly in Asia's arid and semi-desert areas and is a significant commercial and decorative plant. Nevertheless, the genus' classification has been disputed, and the evolutionary connections within the genus have not been thoroughly defined. Hence, we sequenced and analyzed Ajania's plastid genomes and combined them with ETS data to assess their phylogenetic relationships.

RESULTS

We obtained a total of six new Ajania plastid genomes and nine ETS sequences. The whole plastome lengths of the six species sampled ranged from 151,002 bp to 151,115 bp, showing conserved structures. Combined with publicly available data from GenBank, we constructed six datasets to reconstruct the phylogenetic relationships, detecting nucleoplasmic clashes. Our results reveal the affinities of Artemisia, Chrysanthemum and Stilpnolepis to Ajania and validate the early taxonomy reclassification. Some of the plastid genes with low phylogenetic information and gene trees with topological differences may have contributed to the ambiguous phylogenetic results of Ajania. There is extensive evolutionary rate heterogeneity in plastid genes. The psbH and ycf2 genes, which are involved in photosynthesis and ATP transport, are under selective pressure. Plastomes from Ajania species diverged, and structural aspects of plastomes may indicate some of the real evolutionary connections. We suggest the ycf1 gene as a viable plastid DNA barcode because it has significant nucleotide diversity and better reflects evolutionary connections.

CONCLUSION

Our findings validate the early Ajania taxonomy reclassification and show evolutionary rate heterogeneity, genetic variety, and phylogenetic heterogeneity of plastid genes. This research might provide new insights into the taxonomy and evolution of Ajania, as well as provide useful information for germplasm innovation and genetic enhancement in horticultural species.

IRplus: An Augmented Tool to Detect Inverted Repeats in Plastid Genomes

High-throughput sequencing methods have increased the accessibility of plastid genomes, which are crucial for clarifying phylogenetic relationships. Current large sequencing efforts require software tools for routine display of their distinctive quadripartite structure, which is denoted by four junction sites. By concentrating on these junctions and their close vicinity, IRscope has emerged as the standard tool for detection of this structure and creating simplified comparative graphical maps of plastid genomes. Here, we provide an augmented version (IRplus) that encompasses a novel set of functions such as integrated error detection, flexible color schemes, and an upgraded method to detect inverted repeats in genomic sequences. Spanning across the plant tree of life, IRplus allows the quick visualization of various sets of plastid genomes and features, next to smooth interoperability with other widely used annotation file formats and platforms. The IRplus can be accessed at https://irscope.shinyapps.io/IRplus/, and source codes are freely available at https://github.com/AmiryousefiLab/IRplus.

Genomic Analysis of Plastid-Nuclear Interactions and Differential Evolution Rates in Coevolved Genes across Juglandaceae Species

The interaction between the nuclear and chloroplast genomes in plants is crucial for preserving essential cellular functions in the face of varying rates of mutation, levels of selection, and modes of transmission. Despite this, identifying nuclear genes that coevolve with chloroplast genomes at a genome-wide level has remained a challenge. In this study, we conducted an evolutionary rate covariation analysis to identify candidate nuclear genes coevolving with chloroplast genomes in Juglandaceae. Our analysis was based on 4,894 orthologous nuclear genes and 76 genes across seven chloroplast partitions in nine Juglandaceae species. Our results indicated that 1,369 (27.97%) of the nuclear genes demonstrated signatures of coevolution, with the Ycf1/2 partition yielding the largest number of hits (765) and the ClpP1 partition yielding the fewest (13). These hits were found to be significantly enriched in biological processes related to leaf development, photoperiodism, and response to abiotic stress. Among the seven partitions, AccD, ClpP1, MatK, and RNA polymerase partitions and their respective hits exhibited a narrow range, characterized by dN/dS values below 1. In contrast, the Ribosomal, Photosynthesis, Ycf1/2 partitions and their corresponding hits, displayed a broader range of dN/dS values, with certain values exceeding 1. Our findings highlight the differences in the number of candidate nuclear genes coevolving with the seven chloroplast partitions in Juglandaceae species and the correlation between the evolution rates of these genes and their corresponding chloroplast partitions.

Target capture and genome skimming for plant diversity studies

Recent technological advances in long-read high-throughput sequencing and assembly methods have facilitated the generation of annotated chromosome-scale whole-genome sequence data for evolutionary studies; however, generating such data can still be difficult for many plant species. For example, obtaining high-molecular-weight DNA is typically impossible for samples in historical herbarium collections, which often have degraded DNA. The need to fast-freeze newly collected living samples to conserve high-quality DNA can be complicated when plants are only found in remote areas. Therefore, short-read reduced-genome representations, such as target capture and genome skimming, remain important for evolutionary studies. Here, we review the pros and cons of each technique for non-model plant taxa. We provide guidance related to logistics, budget, the genomic resources previously available for the target clade, and the nature of the study. Furthermore, we assess the available bioinformatic analyses, detailing best practices and pitfalls, and suggest pathways to combine newly generated data with legacy data. Finally, we explore the possible downstream analyses allowed by the type of data generated using each technique. We provide a practical guide to help researchers make the best-informed choice regarding reduced genome representation for evolutionary studies of non-model plants in cases where whole-genome sequencing remains impractical.

Phylogenomics and plastome evolution of Indigofera (Fabaceae)

Introduction

Indigofera L. is the third largest genus in Fabaceae and includes economically important species that are used for indigo dye-producing, medicinal, ornamental, and soil and water conservation. The genus is taxonomically difficult due to the high level of overlap in morphological characters of interspecies, fewer reliability states for classification, and extensive adaptive evolution. Previous characteristic-based taxonomy and nuclear ITS-based phylogenies have contributed to our understanding of Indigofera taxonomy and evolution. However, the lack of chloroplast genomic resources limits our comprehensive understanding of the phylogenetic relationships and evolutionary processes of Indigofera.

Methods

Here, we newly assembled 18 chloroplast genomes of Indigofera. We performed a series of analyses of genome structure, nucleotide diversity, phylogenetic analysis, species pairwise Ka/Ks ratios, and positive selection analysis by combining with allied species in Papilionoideae.

Results and discussion

The chloroplast genomes of Indigofera exhibited highly conserved structures and ranged in size from 157,918 to 160,040 bp, containing 83 protein-coding genes, 37 tRNA genes, and eight rRNA genes. Thirteen highly variable regions were identified, of which trnK-rbcL, ndhF-trnL, and ycf1 were considered as candidate DNA barcodes for species identification of Indigofera. Phylogenetic analysis using maximum likelihood (ML) and Bayesian inference (BI) methods based on complete chloroplast genome and protein-coding genes (PCGs) generated a well-resolved phylogeny of Indigofera and allied species. Indigofera monophyly was strongly supported, and four monophyletic lineages (i.e., the Pantropical, East Asian, Tethyan, and Palaeotropical clades) were resolved within the genus. The species pairwise Ka/Ks ratios showed values lower than 1, and 13 genes with significant posterior probabilities for codon sites were identified in the positive selection analysis using the branch-site model, eight of which were associated with photosynthesis. Positive selection of accD suggested that Indigofera species have experienced adaptive evolution to selection pressures imposed by their herbivores and pathogens. Our study provided insight into the structural variation of chloroplast genomes, phylogenetic relationships, and adaptive evolution in Indigofera. These results will facilitate future studies on species identification, interspecific and intraspecific delimitation, adaptive evolution, and the phylogenetic relationships of the genus Indigofera.

Comparative analysis of dioecious Amaranthus plastomes and phylogenomic implications within Amaranthaceae s.s

BACKGROUND

The genus Amaranthus L. consists of 70-80 species distributed across temperate and tropical regions of the world. Nine species are dioecious and native to North America; two of which are agronomically important weeds of row crops. The genus has been described as taxonomically challenging and relationships among species including the dioecious ones are poorly understood. In this study, we investigated the phylogenetic relationships among the dioecious amaranths and sought to gain insights into plastid tree incongruence. A total of 19 Amaranthus species' complete plastomes were analyzed. Among these, seven dioecious Amaranthus plastomes were newly sequenced and assembled, an additional two were assembled from previously published short reads sequences and 10 other plastomes were obtained from a public repository (GenBank).

RESULTS

Comparative analysis of the dioecious Amaranthus species' plastomes revealed sizes ranged from 150,011 to 150,735 bp and consisted of 112 unique genes (78 protein-coding genes, 30 transfer RNAs and 4 ribosomal RNAs). Maximum likelihood trees, Bayesian inference trees and splits graphs support the monophyly of subgenera Acnida (7 dioecious species) and Amaranthus; however, the relationship of A. australis and A. cannabinus to the other dioecious species in Acnida could not be established, as it appears a chloroplast capture occurred from the lineage leading to the Acnida + Amaranthus clades. Our results also revealed intraplastome conflict at some tree branches that were in some cases alleviated with the use of whole chloroplast genome alignment, indicating non-coding regions contribute valuable phylogenetic signals toward shallow relationship resolution. Furthermore, we report a very low evolutionary distance between A. palmeri and A. watsonii, indicating that these two species are more genetically related than previously reported.

CONCLUSIONS

Our study provides valuable plastome resources as well as a framework for further evolutionary analyses of the entire Amaranthus genus as more species are sequenced.

Characterizing conflict and congruence of molecular evolution across organellar genome sequences for phylogenetics in land plants

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.

A comparative study across the parasitic plants of Cuscuta subgenus Grammica (Convolvulaceae) reveals a possible loss of the plastid genome in its section Subulatae

Most species in Cuscuta subgenus Grammica retain many photosynthesis-related plastid genes, generally under purifying selection. A group of holoparasitic species in section Subulatae may have lost their plastid genomes entirely. The c. 153 species of plants belonging to Cuscuta subgenus Grammica are all obligate stem parasites. However, some have completely lost the ability to conduct photosynthesis while others retain photosynthetic machinery and genes. The plastid genome that primarily encodes key photosynthesis genes functions as a bellwether for how reliant plants are on primary production. This research assembles and analyses 17 plastomes across Cuscuta subgenus Grammica with the aim of characterizing the state of the plastome in each of its sections. By comparing the structure and content of plastid genomes across the subgenus, as well as by quantifying the selection acting upon each gene, we reconstructed the patterns of plastome change within the phylogenetic context for this group. We found that species in 13 of the 15 sections that comprise Grammica retain the bulk of plastid photosynthesis genes and are thus hemiparasitic. The complete loss of photosynthesis can be traced to two clades: the entire section Subulatae and a complex of three species within section Ceratophorae. We were unable to recover any significant plastome sequences from section Subulatae, suggesting that plastomes in these species are either drastically reduced or lost entirely.

Plastid Phylogenomic Insights into the Inter-Tribal Relationships of Plantaginaceae

Plantaginaceae, consisting of 12 tribes, is a diverse, cosmopolitan family. To date, the inter-tribal relationships of this family have been unresolved, and the plastome structure and composition within Plantaginaceae have seldom been comprehensively investigated. In this study, we compared the plastomes from 41 Plantaginaceae species (including 6 newly sequenced samples and 35 publicly representative species) representing 11 tribes. To clarify the inter-tribal relationships of Plantaginaceae, we inferred phylogenic relationships based on the concatenated and coalescent analyses of 68 plastid protein-coding genes. PhyParts analysis was performed to assess the level of concordance and conflict among gene trees across the species tree. The results indicate that most plastomes of Plantaginaceae are largely conserved in terms of genome structure and gene content. In contrast to most previous studies, a robust phylogeny was recovered using plastome data, providing new insights for better understanding the inter-tribal relationships of Plantaginaceae. Both concatenated and coalescent phylogenies favored the sister relationship between Plantagineae and Digitalideae, as well as between Veroniceae and Hemiphragmeae. Sibthorpieae diverged into a separate branch which was sister to a clade comprising the four tribes mentioned above. Furthermore, the sister relationship between Russelieae and Cheloneae is strongly supported. The results of PhyParts showed gene tree congruence and conflict to varying degrees, but most plastid genes were uninformative for phylogenetic nodes, revealing the defects of previous studies using single or multiple plastid DNA sequences to infer the phylogeny of Plantaginaceae.

Phylogenomics and the flowering plant tree of life

The advances accelerated by next-generation sequencing and long-read sequencing technologies continue to provide an impetus for plant phylogenetic study. In the past decade, a large number of phylogenetic studies adopting hundreds to thousands of genes across a wealth of clades have emerged and ushered plant phylogenetics and evolution into a new era. In the meantime, a roadmap for researchers when making decisions across different approaches for their phylogenomic research design is imminent. This review focuses on the utility of genomic data (from organelle genomes, to both reduced representation sequencing and whole-genome sequencing) in phylogenetic and evolutionary investigations, describes the baseline methodology of experimental and analytical procedures, and summarizes recent progress in flowering plant phylogenomics at the ordinal, familial, tribal, and lower levels. We also discuss the challenges, such as the adverse impact on orthology inference and phylogenetic reconstruction raised from systematic errors, and underlying biological factors, such as whole-genome duplication, hybridization/introgression, and incomplete lineage sorting, together suggesting that a bifurcating tree may not be the best model for the tree of life. Finally, we discuss promising avenues for future plant phylogenomic studies.

Chloroplast phylogenomics and the taxonomy of Saxifraga section Ciliatae (Saxifragaceae)

Comprising ca. 200 species, Saxifraga sect. Ciliatae is the most species-rich section of Saxifraga s.str., whose center of diversity is in the Tibeto-Himalayan region. The infra-sectional classification of sect. Ciliatae is still in debate due to the high level of species richness, as well as remarkable variations of habitat, morphology, physiology and life cycles. Subdivisions of this section proposed in various taxonomic systems have not been adequately tested in previous phylogenetic studies, partly due to low taxonomic sampling density, but also to the use of few DNA markers. In order to achieve a more robust infra-sectional classification of sect. Ciliatae, complete chloroplast genomes of 94 taxa from this section were analyzed, of which 93 were newly sequenced, assembled and annotated. The length of the 94 plastomes of sect. Ciliatae taxa range from 143,479 to 159,938 bp, encoding 75 to 79 unique protein-coding genes (PCGs). Analyses of the 94 plastomes revealed high conservation in structural organization, gene arrangement, and gene content. Gene loss and changes of IR boundaries were detected but in extremely low frequency. The molecular phylogenetic tree from concatenated PCGs and complete chloroplast genome sequences exhibited high resolution and support values and confirms that sect. Ciliatae is monophyletic. Three well-supported clades were revealed within the section that agree relatively well with the subsectional taxonomy of Gornall (1987), but some minor modifications should be made. Firstly, the monotypic subsection Cinerascentes should be abandoned and its constituent species, S. cinerascens, assigned to subsect. Gemmiparae. Secondly, subsections Rosulares and Serpyllifoliae should be merged and become subsect. Rosulares. Section Ciliatae thus comprises: subsect. Hirculoideae Engl. & Irmsch.; subsect. Rosulares Gornall; subsect. Gemmiparae Engl. & Irmsch.; subsect. Flagellares (C. B. Clarke) Engl. & Irmsch. and subsect. Hemisphaericae (Engl. & Irmsch.) Gornall.

The complete chloroplast genome of a fast-growing tree Lophostemon confertus (Myrtaceae)

Lophostemon confertus (Myrtaceae), a fast-growing ornamental tree, is widely cultivated in tropical and subtropical regions. To determine its phylogenetic position within Myrtaceae, here we report its complete chloroplast (cp) genome, which is 160,297 bp long and contains two inverted repeats (IRs) of 26,490 bp each, separated by a small single-copy region of 18,826 bp and a large single-copy region of 88,491 bp. The cp genome contains 123 genes, including 73 unique protein-coding genes (six duplicated in the IR regions), 29 unique tRNA genes (seven duplicated in the IR regions), and four unique rRNA genes (all located in the IR regions). Phylogenetic analysis of 18 species of Myrtaceae showed that L. confertus is sister to Xanthostemon chrysanthus. The complete cp genome of L. confertus provides a valuable genetic resource for further phylogenetic studies.

Detecting signals of adaptive evolution in grape plastomes with a focus on the Cretaceous-Palaeogene (K/Pg) transition

BACKGROUND AND AIMS

Although plastid genes are widely used in phylogenetic studies, signals of positive selection have been scarcely investigated in the grape family. The plastomes from 91 accessions of Vitaceae were examined to understand the extent to which positive selection is present and to identify which genes are involved. Moreover, the changes through time of genes under episodic positive selection were investigated and the hypothesis of an adaptive process following the Cretaceous-Palaeogene (K/Pg) transition about 66 million years ago was tested.

METHODS

Different codon-substitution models were used to assess pervasive and episodic positive selection events on 70 candidate plastid genes. Divergence times between lineages were estimated and stochastic character mapping analysis was used to simulate variation over time of the genes found to be under episodic positive selection.

KEY RESULTS

A total of 20 plastid genes (29 %) showed positive selection. Among them, 14 genes showed pervasive signatures of positive selection and nine genes showed episodic signatures of positive selection. In particular, four of the nine genes (psbK, rpl20, rpoB, rps11) exhibited a similar pattern showing an increase in the rate of variation close to the K/Pg transition.

CONCLUSION

Multiple analyses have shown that the grape family has experienced ancient and recent positive selection events and that the targeted genes are involved in essential functions such as photosynthesis, self-replication and metabolism. Our results are consistent with the idea that the K/Pg transition has favoured an increased rate of change in some genes. Intense environmental perturbations have influenced the rapid diversification of certain lineages, and new mutations arising on some plastid genes may have been fixed by natural selection over the course of many generations.

Plastome phylogenomics provide new perspective into the phylogeny and evolution of Betulaceae (Fagales)

BACKGROUND

Betulaceae is a relatively small but morphologically diverse family, with many species having important economic and ecological values. Although plastome structure of Betulaceae has been reported sporadically, a comprehensive exploration for plastome evolution is still lacking. Besides, previous phylogenies had been constructed based on limited gene fragments, generating unrobust phylogenetic framework and hindering further studies on divergence ages, biogeography and character evolution. Here, 109 plastomes (sixteen newly assembled and 93 previously published) were subject to comparative genomic and phylogenomic analyses to reconstruct a robust phylogeny and trace the diversification history of Betulaceae.

RESULTS

All Betulaceae plastomes were highly conserved in genome size, gene order, and structure, although specific variations such as gene loss and IR boundary shifts were revealed. Ten divergent hotspots, including five coding regions (P_i > 0.02) and five noncoding regions (P_i > 0.035), were identified as candidate DNA barcodes for phylogenetic analysis and species delimitation. Phylogenomic analyses yielded high-resolution topology that supported reciprocal monophyly between Betula and Alnus within Betuloideae, and successive divergence of Corylus, Ostryopsis, and Carpinus-Ostrya within Coryloideae. Incomplete lineage sorting and hybridization may be responsible for the mutual paraphyly between Ostrya and Carpinus. Betulaceae ancestors originated from East Asia during the upper Cretaceous; dispersals and subsequent vicariance accompanied by historical environment changes contributed to its diversification and intercontinental disjunction. Ancestral state reconstruction indicated the acquisition of many taxonomic characters was actually the results of parallel or reversal evolution.

CONCLUSIONS

Our research represents the most comprehensive taxon-sampled and plastome-level phylogenetic inference for Betulaceae to date. The results clearly document global patterns of plastome structural evolution, and established a well-supported phylogeny of Betulaceae. The robust phylogenetic framework not only provides new insights into the intergeneric relationships, but also contributes to a perspective on the diversification history and evolution of the family.

Comparative chloroplast genomics provides insights into the genealogical relationships of endangered Tetraena mongolica and the chloroplast genome evolution of related Zygophyllaceae species

A comprehensive understanding of genetic background for rare species will provide an important theoretical basis for the future species management, monitoring and conservation. Tetraena mongolica is restrictedly distributed in the western Ordos plateau of China and has been listed as a national protected plant. We generated 13 chloroplast (cp) genomes of T. mongolica (size range of 106,062-106,230 bp) and conducted a series of comparative analyses of six Zygophyllaceae cp genomes. T. mongolica cp genome exhibited a quadripartite structure with drastically reduced inverted repeats (IRs, 4,315 bp) and undergone the loss of a suit of ndh genes and a copy of rRNAs. Furthermore, all the T. mongolica populations were divided into two genetic groups based on complete cp phylogenomics. In addition, notably variable genome size, gene order and structural changes had been observed among the six Zygophyllaceae cp genomes. Overall, our findings provide insights into the cp genome evolution mode and intraspecific relationships of T. mongolica, and provide a molecular basis for scientific conservation of this endangered plant.

New insights into the phylogeny of the complex thalloid liverworts (Marchantiopsida) based on chloroplast genomes

Marchantiopsida (complex thalloid liverworts) are one of the earliest lineages of embryophytes (land plants), and well-known for their air pores and chambers, pegged rhizoids, and absence of organellular RNA editing sites. Despite their importance to an understanding of early embryophyte evolution, many key nodes within this class remain poorly resolved, owing to the paucity of genetic loci previously available for phylogenetic analyses. Here, we sequenced 54 plastomes, representing 28 genera, nearly all families, and all orders of Marchantiopsida. Based on these plastomes, we present a hypothesis of deep relationships within the class, and make the first investigations of gene contents and synteny. Overall, the Marchantiopsida plastomes were well-conserved, with the exception of the genus Cyathodium that has plastomes with higher GC content, fewer single sequence repeats (SSRs), and more structural variations, implying that this genus might possess RNA editing sites. Abundant repetitive elements and six highly divergent regions were identified as suitable for future infrafamilial taxonomic studies. The phylogenetic topology of Sphaerocarpales, Neohodgsoniales and Blasiales within Marchantiopsida was essentially congruent with previous studies but generally we obtained higher support values. Based on molecular evidence and previous morphological studies, we include Lunulariales in Marchantiales and suggest the retention of narrowed delimitation of monotypic families. The phylogenetic relationships within Marchantiales were better resolved, and 13 monophyletic families were recovered. Our analyses confirmed that the loss of intron 2 of ycf3 is a synapomorphy of Marchantiidae. Finally, we propose a new genus, Asterellopsis (Aytoniaceae), and present an updated classification of Marchantiopsida. The highly supported phylogenetic backbone provided here establishes a framework for future comparative and evolutionary studies of the complex thalloid liverworts.

Untying the Gordian knot of plastid phylogenomic conflict: A case from ferns

Phylogenomic studies based on plastid genome have resolved recalcitrant relationships among various plants, yet the phylogeny of Dennstaedtiaceae at the level of family and genera remains unresolved due to conflicting plastid genes, limited molecular data and incomplete taxon sampling of previous studies. The present study generated 30 new plastid genomes of Dennstaedtiaceae (9 genera, 29 species), which were combined with 42 publicly available plastid genomes (including 24 families, 27 genera, 42 species) to explore the evolution of Dennstaedtiaceae. In order to minimize the impact of systematic errors on the resolution of phylogenetic inference, we applied six strategies to generate 30 datasets based on CDS, intergenic spacers, and whole plastome, and two tree inference methods (maximum-likelihood, ML; and multispecies coalescent, MSC) to comprehensively analyze the plastome-scale data. Besides, the phylogenetic signal among all loci was quantified for controversial nodes using ML framework, and different topologies hypotheses among all datasets were tested. The species trees based on different datasets and methods revealed obvious conflicts at the base of the polypody ferns. The topology of the "CDS-codon-align-rm3" (CDS with the removal of the third codon) matrix was selected as the primary reference or summary tree. The final phylogenetic tree supported Dennstaedtiaceae as the sister group to eupolypods, and Dennstaedtioideae was divided into four clades with full support. This robust reconstructed phylogenetic backbone establishes a framework for future studies on Dennstaedtiaceae classification, evolution and diversification. The present study suggests considering plastid phylogenomic conflict when using plastid genomes. From our results, reducing saturated genes or sites can effectively mitigate tree conflicts for distantly related taxa. Moreover, phylogenetic trees based on amino acid sequences can be used as a comparison to verify the confidence of nucleotide-based trees.

Phylotranscriptomics Illuminates the Placement of Whole Genome Duplications and Gene Retention in Ferns

Ferns are the second largest clade of vascular plants with over 10,000 species, yet the generation of genomic resources for the group has lagged behind other major clades of plants. Transcriptomic data have proven to be a powerful tool to assess phylogenetic relationships, using thousands of markers that are largely conserved across the genome, and without the need to sequence entire genomes. We assembled the largest nuclear phylogenetic dataset for ferns to date, including 2884 single-copy nuclear loci from 247 transcriptomes (242 ferns, five outgroups), and investigated phylogenetic relationships across the fern tree, the placement of whole genome duplications (WGDs), and gene retention patterns following WGDs. We generated a well-supported phylogeny of ferns and identified several regions of the fern phylogeny that demonstrate high levels of gene tree-species tree conflict, which largely correspond to areas of the phylogeny that have been difficult to resolve. Using a combination of approaches, we identified 27 WGDs across the phylogeny, including 18 large-scale events (involving more than one sampled taxon) and nine small-scale events (involving only one sampled taxon). Most inferred WGDs occur within single lineages (e.g., orders, families) rather than on the backbone of the phylogeny, although two inferred events are shared by leptosporangiate ferns (excluding Osmundales) and Polypodiales (excluding Lindsaeineae and Saccolomatineae), clades which correspond to the majority of fern diversity. We further examined how retained duplicates following WGDs compared across independent events and found that functions of retained genes were largely convergent, with processes involved in binding, responses to stimuli, and certain organelles over-represented in paralogs while processes involved in transport, organelles derived from endosymbiotic events, and signaling were under-represented. To date, our study is the most comprehensive investigation of the nuclear fern phylogeny, though several avenues for future research remain unexplored.

New Insights Into the Relationships Within Subtribe Scorzonerinae (Cichorieae, Asteraceae) Using Hybrid Capture Phylogenomics (Hyb-Seq)

Subtribe Scorzonerinae (Cichorieae, Asteraceae) contains 12 main lineages and approximately 300 species. Relationships within the subtribe, either at inter- or intrageneric levels, were largely unresolved in phylogenetic studies to date, due to the lack of phylogenetic signal provided by traditional Sanger sequencing markers. In this study, we employed a phylogenomics approach (Hyb-Seq) that targets 1,061 nuclear-conserved ortholog loci designed for Asteraceae and obtained chloroplast coding regions as a by-product of off-target reads. Our objectives were to evaluate the potential of the Hyb-Seq approach in resolving the phylogenetic relationships across the subtribe at deep and shallow nodes, investigate the relationships of major lineages at inter- and intrageneric levels, and examine the impact of the different datasets and approaches on the robustness of phylogenetic inferences. We analyzed three nuclear datasets: exon only, excluding all potentially paralogous loci; exon only, including loci that were only potentially paralogous in 1-3 samples; exon plus intron regions (supercontigs); and the plastome CDS region. Phylogenetic relationships were reconstructed using both multispecies coalescent and concatenation (Maximum Likelihood and Bayesian analyses) approaches. Overall, our phylogenetic reconstructions recovered the same monophyletic major lineages found in previous studies and were successful in fully resolving the backbone phylogeny of the subtribe, while the internal resolution of the lineages was comparatively poor. The backbone topologies were largely congruent among all inferences, but some incongruent relationships were recovered between nuclear and plastome datasets, which are discussed and assumed to represent cases of cytonuclear discordance. Considering the newly resolved phylogenies, a new infrageneric classification of Scorzonera in its revised circumscription is proposed.

Comparative Chloroplast Genomics of Seven Endangered Cypripedium Species and Phylogenetic Relationships of Orchidaceae

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.

Plastid phylogenomic analyses of the Selaginella sanguinolenta group (Selaginellaceae) reveal conflict signatures resulting from sequence types, outlier genes, and pervasive RNA editing

Different from the generally conserved plastomes (plastid genomes) of most land plants, the Selaginellaceae plastomes exhibit dynamic structure, high GC content and high substitution rates. Previous plastome analyses identified strong conflict on several clades in Selaginella, however the factors causing the conflictions and the impact on the phylogenetic inference have not been sufficiently investigated. Here, we dissect the distribution of phylogenetic signals and conflicts in Selaginella sanguinolenta group, the plastome of which is DR (direct repeats) structure and with genome-wide RNA editing. We analyzed the data sets including 22 plastomes representing all species of the S. sanguinolenta group, covering the entire geographical distribution from the Himalayas to Siberia and the Russian Far East regions. We recovered four different topologies by applying multispecies coalescent (ASTRAL) and concatenation methods (IQ-TREE and RAxML) on four data sets of PC (protein-coding genes), NC (non-coding sequences), PCN (the concatenated PC and NC), and RC (predicted RNA editing sites "C" were corrected by "T"), respectively. Six monophyletic clades, S. nummularifolia clade, S. rossii clade, S. sajanensis clade, S. sanguinolenta I clade, S. sanguinolenta II clade, and S. sanguinolenta III clade, were consistently resolved and supported by the characteristics of GC content, RNA editing frequency, and gene content. However, the relationships among these clades varied across the four topologies. To explore the underlying causes of the uncertainty, we compared the phylogenetic signals of the four topologies. We identified that the sequence types (coding versus non-coding), outlier genes (genes with extremely high |ΔGLS| values), and C-to-U RNA editing frequency in the protein-coding genes were responsible for the unstable phylogenomic relationship. We further revealed a significant positive correlation between the |ΔGLS| values and the variation coefficient of the RNA editing number. Our results demonstrated that the coalescent method performed better than the concatenation method in overcoming the problems caused by outlier genes and extreme RNA editing events. Our study particularly focused on the importance of exploring the plastid phylogenomic conflicts and suggested conducting concatenated analyses cautiously when adopting organelle genome data.

Increased resolution in the face of conflict: phylogenomics of the Neotropical bellflowers (Campanulaceae: Lobelioideae), a rapid plant radiation

BACKGROUND AND AIMS

The centropogonid clade (Lobelioideae: Campanulaceae) is an Andean-centred rapid radiation characterized by repeated convergent evolution of morphological traits, including fruit type and pollination syndromes. While previous studies have resolved relationships of lineages with fleshy fruits into subclades, relationships among capsular species remain unresolved. This lack of resolution has impeded reclassification of non-monophyletic genera, whose current taxonomy relies heavily on traits that have undergone convergent evolution.

METHODS

Targeted sequence capture using a probe-set recently developed for the centropogonid clade was used to obtain phylogenomic data from DNA extracted from both silica-dried and herbarium leaf tissue. These data were used to infer relationships among species using concatenated and partitioned species tree methods, and to quantify gene tree discordance.

KEY RESULTS

While silica-dried leaf tissue resulted in longer assembled sequence data, the inclusion of herbarium samples improved taxonomic representation. Relationships among baccate lineages are similar to those inferred in previous studies, although they differ for lineages within and among capsular clades. We improve the phylogenetic resolution of Siphocampylus, which forms ten groups of closely related species which we informally name. Two subclades of Siphocampylus and two individual species are rogue taxa whose placement differs widely across analyses. Gene tree discordance (including cytonuclear discordance) is rampant.

CONCLUSIONS

This first phylogenomic study of the centropogonid clade considerably improves our understanding of relationships in this rapid radiation. Differences across analyses and the possibility of additional lineage discoveries still hamper a solid and stable reclassification. Rapid morphological innovation corresponds with a high degree of phylogenomic complexity, including cytonuclear discordance, nuclear gene tree conflict and well-supported differences between analyses based on different nuclear loci. Together, these results point to a potential role of hemiplasy underlying repeated convergent evolution. This hallmark of rapid radiations is probably present in many other species-rich Andean plant radiations.

Phylogenomic conflict analyses in the apple genus Malus s.l. reveal widespread hybridization and allopolyploidy driving diversification, with insights into the complex biogeographic history in the Northern Hemisphere

Phylogenomic evidence from an increasing number of studies has demonstrated that different data sets and analytical approaches often reconstruct strongly supported but conflicting relationships. In this study, 785 single-copy nuclear genes and 75 complete plastomes were used to infer the phylogenetic relationships and estimate the historical biogeography of the apple genus Malus sensu lato, an economically important lineage disjunctly distributed in the Northern Hemisphere and involved in known and suspected hybridization and allopolyploidy events. The nuclear phylogeny recovered the monophyly of Malus s.l. (including Docynia); however, the genus was supported to be biphyletic in the plastid phylogeny. An ancient chloroplast capture event in the Eocene in western North America best explains the cytonuclear discordance. Our conflict analysis demonstrated that ILS, hybridization, and allopolyploidy could explain the widespread nuclear gene tree discordance. One deep hybridization event (Malus doumeri) and one recent event (Malus coronaria) were detected in Malus s.l. Furthermore, our historical biogeographic analysis integrating living and fossil data supported a widespread East Asian-western North American origin of Malus s.l. in the Eocene, followed by several extinction and dispersal events in the Northern Hemisphere. We also propose a general workflow for assessing phylogenomic discordance and biogeographic analysis using deep genome skimming data sets.

New plastomes of eight Ipomoea species and four putative hybrids from Eastern Amazon

Ipomoea is a large pantropical genus globally distributed, which importance goes beyond the economic value as food resources or ornamental crops. This highly diverse genus has been the focus of a great number of studies, enriching the plant genomics knowledge, and challenging the plant evolution models. In the Carajás mountain range, located in Eastern Amazon, the savannah-like ferruginous ecosystem known as canga harbors highly specialized plant and animal populations, and Ipomoea is substantially representative in such restrictive habitat. Thus, to provide genetic data and insights into whole plastome phylogenetic relationships among key Ipomoea species from Eastern Amazon with little to none previously available data, we present the complete plastome sequences of twelve lineages of the genus, including the canga microendemic I. cavalcantei, the closely related I. marabaensis, and their putative hybrids. The twelve plastomes presented similar gene content as most publicly available Ipomoea plastomes, although the putative hybrids were correctly placed as closely related to the two parental species. The cavalcantei-marabaensis group was consistently grouped between phylogenetic methods. The closer relationship of the I. carnea plastome with the cavalcantei-marabaensis group, as well as the branch formed by I. quamoclit, I. asarifolia and I. maurandioides, were probably a consequence of insufficient taxonomic representativity, instead of true genetic closeness, reinforcing the importance of new plastome assemblies to resolve inconsistencies and boost statistical confidence, especially the case for South American clades of Ipomoea. The search for k-mers presenting high dispersion among the frequency distributions pointed to highly variable coding and intergenic regions, which may potentially contribute to the genetic diversity observed at species level. Our results contribute to the resolution of uncertain clades within Ipomoea and future phylogenomic studies, bringing unprecedented results to Ipomoea species with restricted distribution, such as I. cavalcantei.

Towards the plastome evolution and phylogeny of Cycas L. (Cycadaceae): molecular-morphology discordance and gene tree space analysis

BACKGROUND

Plastid genomes (plastomes) present great potential in resolving multiscale phylogenetic relationship but few studies have focused on the influence of genetic characteristics of plastid genes, such as genetic variation and phylogenetic discordance, in resolving the phylogeny within a lineage. Here we examine plastome characteristics of Cycas L., the most diverse genus among extant cycads, and investigate the deep phylogenetic relationships within Cycas by sampling 47 plastomes representing all major clades from six sections.

RESULTS

All Cycas plastomes shared consistent gene content and structure with only one gene loss detected in Philippine species C. wadei. Three novel plastome regions (psbA-matK, trnN-ndhF, chlL-trnN) were identified as containing the highest nucleotide variability. Molecular evolutionary analysis showed most of the plastid protein-coding genes have been under purifying selection except ndhB. Phylogenomic analyses that alternatively included concatenated and coalescent methods, both identified four clades but with conflicting topologies at shallow nodes. Specifically, we found three species-rich Cycas sections, namely Stangerioides, Indosinenses and Cycas, were not or only weakly supported as monophyly based on plastomic phylogeny. Tree space analyses based on different tree-inference methods both revealed three gene clusters, of which the cluster with moderate genetic properties showed the best congruence with the favored phylogeny.

CONCLUSIONS

Our exploration in plastomic data for Cycas supports the idea that plastid protein-coding genes may exhibit discordance in phylogenetic signals. The incongruence between molecular phylogeny and morphological classification reported here may largely be attributed to the uniparental attribute of plastid, which cannot offer sufficient information to resolve the phylogeny. Contrasting to a previous consensus that genes with longer sequences and a higher proportion of variances are superior for phylogeny reconstruction, our result implies that the most effective phylogenetic signals could come from loci that own moderate variation, GC content, sequence length, and underwent modest selection.

Highly Resolved Papilionoid Legume Phylogeny Based on Plastid Phylogenomics

Comprising 501 genera and around 14,000 species, Papilionoideae is not only the largest subfamily of Fabaceae (Leguminosae; legumes), but also one of the most extraordinarily diverse clades among angiosperms. Papilionoids are a major source of food and forage, are ecologically successful in all major biomes, and display dramatic variation in both floral architecture and plastid genome (plastome) structure. Plastid DNA-based phylogenetic analyses have greatly improved our understanding of relationships among the major groups of Papilionoideae, yet the backbone of the subfamily phylogeny remains unresolved. In this study, we sequenced and assembled 39 new plastomes that are covering key genera representing the morphological diversity in the subfamily. From 244 total taxa, we produced eight datasets for maximum likelihood (ML) analyses based on entire plastomes and/or concatenated sequences of 77 protein-coding sequences (CDS) and two datasets for multispecies coalescent (MSC) analyses based on individual gene trees. We additionally produced a combined nucleotide dataset comprising CDS plus matK gene sequences only, in which most papilionoid genera were sampled. A ML tree based on the entire plastome maximally supported all of the deep and most recent divergences of papilionoids (223 out of 236 nodes). The Swartzieae, ADA (Angylocalyceae, Dipterygeae, and Amburaneae), Cladrastis, Andira, and Exostyleae clades formed a grade to the remainder of the Papilionoideae, concordant with nine ML and two MSC trees. Phylogenetic relationships among the remaining five papilionoid lineages (Vataireoid, Dermatophyllum, Genistoid s.l., Dalbergioid s.l., and Baphieae + Non-Protein Amino Acid Accumulating or NPAAA clade) remained uncertain, because of insufficient support and/or conflicting relationships among trees. Our study fully resolved most of the deep nodes of Papilionoideae, however, some relationships require further exploration. More genome-scale data and rigorous analyses are needed to disentangle phylogenetic relationships among the five remaining lineages.

Incompatibility Phylogenetic Signals between Double-Digest Restriction Site-Associated DNA Sequencing and Plastid Genomes in Chinese Curcuma (Zingiberaceae)—A Recent Qinghai–Tibetan Plateau Diversification Genera

Curcuma is of high economic value, credited to its medicinal, edible, and ornamental properties, which possess all signatures of adaptability, and rapid radiation, especially species of Curcuma (Chinese Curcuma, a recent Qinghai–Tibetan Plateau diversification genera) scattered in China. However, little is known about the incongruent phylogenetic signals within this genera from different inheritance patterns that will militate against the further development of this genera. In this research, we applied complete chloroplast genome data together with double-digest restriction site-associated DNA sequencing data (ddRAD-seq) strategy to investigate phylogenetic signals of Chinese Curcuma species, clustering using two RAD analysis pipelines (STACKS and pyRAD). Phylogenetic trees were obtained from each locus based on the maximum likelihood (ML) and multispecies coalescent (BEAST) methods. For visual comparison, multi-method and different datasets were used to infer the phylogeny. We discovered inconsistent relationships for the Chinese Curcuma with varying degrees of support using different methods and datasets.

Phylogenetic Relationships Within the Hyper-Diverse Genus Eugenia (Myrtaceae: Myrteae) Based on Target Enrichment Sequencing

Eugenia is one of the most taxonomically challenging lineages of flowering plants, in which morphological delimitation has changed over the last few years resulting from recent phylogenetic study based on molecular data. Efforts, until now, have been limited to Sanger sequencing of mostly plastid markers. These phylogenetic studies indicate 11 clades formalized as infrageneric groups. However, relationships among these clades are poorly supported at key nodes and inconsistent between studies, particularly along the backbone and within Eugenia sect. Umbellatae encompasses ca. 700 species. To resolve and better understand systematic discordance, 54 Eugenia taxa were subjected to phylogenomic Hyb-Seq using 353 low-copy nuclear genes. Twenty species trees based on coding and non-coding loci of nuclear and plastid datasets were recovered using coalescent and concatenated approaches. Concordant and conflicting topologies were assessed by comparing tree landscapes, topology tests, and gene and site concordance factors. The topologies are similar except between nuclear and plastid datasets. The coalescent trees better accommodate disparity in the intron dataset, which contains more parsimony informative sites, while concatenated trees recover more conservative topologies, as they have narrower distribution in the tree landscape. This suggests that highly supported phylogenetic relationships determined in previous studies do not necessarily indicate overwhelming concordant signal. Congruence must be interpreted carefully especially in concatenated datasets. Despite this, the congruence between the multi-species coalescent (MSC) approach and concatenated tree topologies found here is notable. Our analysis does not support Eugenia subg. Pseudeugenia or sect. Pilothecium, as currently circumscribed, suggesting necessary taxonomic reassessment. Five clades are further discussed within Eugenia sect. Umbellatae progress toward its division into workable clades. While targeted sequencing provides a massive quantity of data that improves phylogenetic resolution in Eugenia, uncertainty still remains in Eugenia sect. Umbellatae. The general pattern of higher site coefficient factor (CF) than gene CF in the backbone of Eugenia suggests stochastic error from limited signal. Tree landscapes in combination with concordance factor scores, as implemented here, provide a comprehensive approach that incorporates several phylogenetic hypotheses. We believe the protocols employed here will be of use for future investigations on the evolutionary history of Myrtaceae.

Organelle Phylogenomics and Extensive Conflicting Phylogenetic Signals in the Monocot Order Poales

The Poales is one of the largest orders of flowering plants with significant economic and ecological values. Reconstructing the phylogeny of the Poales is important for understanding its evolutionary history that forms the basis for biological studies. However, due to sparse taxon sampling and limited molecular data, previous studies have resulted in a variety of contradictory topologies. In particular, there are three nodes surrounded by incongruence: the phylogenetic ambiguity near the root of the Poales tree, the sister family of Poaceae, and the delimitation of the xyrid clade. We conducted a comprehensive sampling and reconstructed the phylogenetic tree using plastid and mitochondrial genomic data from 91 to 66 taxa, respectively, representing all the 16 families of Poales. Our analyses support the finding of Bromeliaceae and Typhaceae as the earliest diverging groups within the Poales while having phylogenetic relationships with the polytomy. The clade of Ecdeiocoleaceae and Joinvilleaceae is recovered as the sister group of Poaceae. The three families, Mayacaceae, Eriocaulaceae, and Xyridaceae, of the xyrid assembly diverged successively along the backbone of the Poales phylogeny, and thus this assembly is paraphyletic. Surprisingly, we find substantial phylogenetic conflicts within the plastid genomes of the Poales, as well as among the plastid, mitochondrial, and nuclear data. These conflicts suggest that the Poales could have a complicated evolutionary history, such as rapid radiation and polyploidy, particularly allopolyploidy through hybridization. In sum, our study presents a new perspicacity into the complex phylogenetic relationships and the underlying phylogenetic conflicts within the Poales.

Plastid phylogenomics of tribe Perseeae (Lauraceae) yields insights into the evolution of East Asian subtropical evergreen broad-leaved forests

BACKGROUND

The East Asian subtropical evergreen broad-leaved forests (EBLFs) harbor remarkable biodiversity. However, their historical assembly remains unclear. To gain new insights into the assembly of this biome, we generated a molecular phylogeny of one of its essential plant groups, the tribe Perseeae (Lauraceae).

RESULTS

Our plastid tree topologies were robust to analyses based on different plastid regions and different strategies for data partitioning, nucleotide substitution saturation, and gap handling. We found that tribe Perseeae comprised six major clades and began to colonize the subtropical EBLFs of East Asia in the early Miocene. The diversification rates of tribe Perseeae accelerated twice in the late Miocene.

CONCLUSIONS

Our findings suggest that the intensified precipitation in East Asia in the early Miocene may have facilitated range expansions of the subtropical EBLFs and establishment of tribe Perseeae within this biome. By the late Miocene, species assembly and diversification within the EBLFs had become rapid.

Plastome structure, evolution, and phylogeny of Selaginella

As one of the earliest land plant lineages, Selaginella is important for studying land plant evolution. It is the largest genus of lycophytes containing 700-800 species. Some unique characters of Selaginella plastomes have been reported, but based only on 20 species. There have been no plastome phylogenies of Selaginella based on a relatively large sampling, and no efforts have been made to resolve the phylogeny of the enigmatic Sinensis group whose relationships have been unclear based on small datasets. Here we investigated the structures of 59 plastomes representing 51 species covering all six subgenera and 18 sections of Selaginella except two sections and including the intriguing Sinensis group for the first time. Our major results include: (1) the plastome size of Selaginella ranges tremendously from 78,492 bp to 187,632 bp; (2) there are numerous gene losses in Selaginella comparing with other lycophytes, Isoëtaceae and Lycopodiaceae; (3) the gene contents and plastome structures in Selaginella vary lineage-specifically and all infrageneric taxa are well supported in the plastome phylogeny; (4) the ndh gene family tends to lose or pseudogenize in those species with DR structure and without other short or medium repeats; (5) the short and medium repeat regions in SC mediate many conformations causing diverse and complex plastome structures, and six new conformations are discovered; (6) forty-eight species sampled have high GC content (>50%) but three species in the Sinensis group have ∼30% GC content in plastomes, similar to most vascular plants; (7) the Sinensis group is monophyletic, includes at least two subgroups, and has the smallest plastomes in land plants except some parasitic plants, and their plastomes do not contain any tRNAs; (8) the younger lineages in Selaginella tend to have higher GC content, whereas the older lineages tend to have lower GC content; and (9) because of incomplete genomic data and abnormal structures or some unknown reasons, even the concatenated plastomes could not well resolve the phylogenetic relationships in Selaginella with confidence, highlighting the difficulty in resolving the phylogeny and evolution of this particularly important land plant lineage.

Perspectives in Myrtaceae evolution from plastomes and nuclear phylogenies

Myrtaceae is a large and species-rich family of woody eudicots, with prevalent distribution in the Southern Hemisphere. Classification and taxonomy of species belonging to this family is quite challenging, sometimes with difficulty in species identification and producing phylogenies with low support for species relationships. Most of the current knowledge comes from few molecular markers, such as plastid genes and intergenic regions, which can be difficult to handle and produce conflicting results. Based on plastid protein-coding sequences and nuclear markers, we present a topology for the phylogenetic relationships among Myrtaceae tribes. Our phylogenetic estimate offers a contrasting topology over previous analysis with fewer markers. Plastome phylogeny groups the tribes Syzygieae and Eucalypteae and individual chloroplast genes produce divergent topologies, especially among species within Myrteae tribe, but also in regard to the grouping of Syzygieae and Eucalypteae. Results are consistent and reproducible with both nuclear and organellar datasets. It confronts previous data about the deep nodes of Myrtaceae phylogeny.

Reassessment of Annamocarya sinesis (Carya sinensis) Taxonomy through Concatenation and Coalescence Phylogenetic Analysis

Due to its peculiar morphological characteristics, there is dispute as to whether the genus of Annamocarya sinensis, a species of Juglandaceae, is Annamocarya or Carya. Most morphologists believe it should be distinguished from the Carya genus while genomicists suggest that A. sinensis belongs to the Carya genus. To explore the taxonomic status of A. sinensis using chloroplast genes, we collected chloroplast genomes of 16 plant species and assembled chloroplast genomes of 10 unpublished Carya species. We analyzed all 26 species' chloroplast genomes through two analytical approaches (concatenation and coalescence), using the entire and unique chloroplast coding sequence (CDS) and entire and protein sequences. Our results indicate that the analysis of the CDS and protein sequences or unique CDS and unique protein sequence of chloroplast genomes shows that A. sinensis indeed belongs to the Carya genus. In addition, our analysis shows that, compared to single chloroplast genes, the phylogeny trees constructed using numerous genes showed higher consistency. Moreover, the phylogenetic analysis calculated with the coalescence method and unique gene sequences was more robust than that done with the concatenation method, particularly for analyzing phylogenetically controversial species. Through the analysis, our results concluded that A. sinensis should be called C. sinensis.

Analysis of Paralogs in Target Enrichment Data Pinpoints Multiple Ancient Polyploidy Events in Alchemilla s.l. (Rosaceae)

Target enrichment is becoming increasingly popular for phylogenomic studies. Although baits for enrichment are typically designed to target single-copy genes, paralogs are often recovered with increased sequencing depth, sometimes from a significant proportion of loci, especially in groups experiencing whole-genome duplication (WGD) events. Common approaches for processing paralogs in target enrichment data sets include random selection, manual pruning, and mainly, the removal of entire genes that show any evidence of paralogy. These approaches are prone to errors in orthology inference or removing large numbers of genes. By removing entire genes, valuable information that could be used to detect and place WGD events is discarded. Here, we used an automated approach for orthology inference in a target enrichment data set of 68 species of Alchemilla s.l. (Rosaceae), a widely distributed clade of plants primarily from temperate climate regions. Previous molecular phylogenetic studies and chromosome numbers both suggested ancient WGDs in the group. However, both the phylogenetic location and putative parental lineages of these WGD events remain unknown. By taking paralogs into consideration and inferring orthologs from target enrichment data, we identified four nodes in the backbone of Alchemilla s.l. with an elevated proportion of gene duplication. Furthermore, using a gene-tree reconciliation approach, we established the autopolyploid origin of the entire Alchemilla s.l. and the nested allopolyploid origin of four major clades within the group. Here, we showed the utility of automated tree-based orthology inference methods, previously designed for genomic or transcriptomic data sets, to study complex scenarios of polyploidy and reticulate evolution from target enrichment data sets.[Alchemilla; allopolyploidy; autopolyploidy; gene tree discordance; orthology inference; paralogs; Rosaceae; target enrichment; whole genome duplication.].

Plastid phylogenomics provides novel insights into the infrafamilial relationship of Polypodiaceae

The polygrammoids (Polypodiaceae) are the most species-rich and diversified epiphytic fern lineages, and hold an important role to understand the deep diverging events and rapid adaptation to changing environments in the plant tree of life. Despite progress in the phylogeny of this group of ferns in previous multilocus phylogenetic studies, uncertainty remains especially in backbone relationships among closely related clades, and the phylogenetic placement of recalcitrant species or lineages. Here, we investigated the deep phylogenetic relationships within Polypodiaceae by sampling all major lineages and using 81 plastid genomes (plastomes), of which 70 plastomes were newly sequenced with high-throughput sequencing technology. Based on parsimony, maximum-likelihood, Bayesian and multispecies coalescent analyses of genome skimming data, we achieved a better resolution of the backbone phylogeny of Polypodiaceae. Using simulated data matrices, we detected that potential phylogenetic artefacts, such as long-branch attraction and insufficient taxonomic sampling, may have a confounding impact on the incongruence of phylogenetic inferences. Furthermore, our phylogenetic analyses offer greater resolution than previous multilocus studies, providing a robust framework for future phylogenetic implications on the subfamilial taxonomy of Polypodiaceae. Our phylogenomic study not only demonstrates the advantage of a character-rich plastome dataset for resolving the recalcitrant lineages that have undergone rapid radiation, but also sheds new light on integrative explorations understanding the evolutionary history of large fern groups in the genomic era.

Phylogeny of Crataegus (Rosaceae) based on 257 nuclear loci and chloroplast genomes: evaluating the impact of hybridization

Background

Hawthorn species (Crataegus L.; Rosaceae tribe Maleae) form a well-defined clade comprising five subgeneric groups readily distinguished using either molecular or morphological data. While multiple subsidiary groups (taxonomic sections, series) are recognized within some subgenera, the number of and relationships among species in these groups are subject to disagreement. Gametophytic apomixis and polyploidy are prevalent in the genus, and disagreement concerns whether and how apomictic genotypes should be recognized taxonomically. Recent studies suggest that many polyploids arise from hybridization between members of different infrageneric groups.

Methods

We used target capture and high throughput sequencing to obtain nucleotide sequences for 257 nuclear loci and nearly complete chloroplast genomes from a sample of hawthorns representing all five currently recognized subgenera. Our sample is structured to include two examples of intersubgeneric hybrids and their putative diploid and tetraploid parents. We queried the alignment of nuclear loci directly for evidence of hybridization, and compared individual gene trees with each other, and with both the maximum likelihood plastome tree and the nuclear concatenated and multilocus coalescent-based trees. Tree comparisons provided a promising, if challenging (because of the number of comparisons involved) method for visualizing variation in tree topology. We found it useful to deploy comparisons based not only on tree-tree distances but also on a metric of tree-tree concordance that uses extrinsic information about the relatedness of the terminals in comparing tree topologies.

Results

We obtained well-supported phylogenies from plastome sequences and from a minimum of 244 low copy-number nuclear loci. These are consistent with a previous morphology-based subgeneric classification of the genus. Despite the high heterogeneity of individual gene trees, we corroborate earlier evidence for the importance of hybridization in the evolution of Crataegus. Hybridization between subgenus Americanae and subgenus Sanguineae was documented for the origin of Sanguineae tetraploids, but not for a tetraploid Americanae species. This is also the first application of target capture probes designed with apple genome sequence. We successfully assembled 95% of 257 loci in Crataegus, indicating their potential utility across the genera of the apple tribe.

Plastid phylogenomics of the Gynoxoid group (Senecioneae, Asteraceae) highlights the importance of motif-based sequence alignment amid low genetic distances

PREMISE

The genus Gynoxys and relatives form a species-rich lineage of Andean shrubs and trees with low genetic distances within the sunflower subtribe Tussilaginineae. Previous molecular phylogenetic investigations of the Tussilaginineae have included few, if any, representatives of this Gynoxoid group or reconstructed ambiguous patterns of relationships for it.

METHODS

We sequenced complete plastid genomes of 21 species of the Gynoxoid group and related Tussilaginineae and conducted detailed comparisons of the phylogenetic relationships supported by the gene, intron, and intergenic spacer partitions of these genomes. We also evaluated the impact of manual, motif-based adjustments of automatic DNA sequence alignments on phylogenetic tree inference.

RESULTS

Our results indicate that the inclusion of all plastid genome partitions is needed to infer well-supported phylogenetic trees of the Gynoxoid group. Whole plastome-based tree inference suggests that the genera Gynoxys and Nordenstamia are polyphyletic and form the core clade of the Gynoxoid group. This clade is sister to a clade of Aequatorium and Paragynoxys and also includes some but not all representatives of Paracalia.

CONCLUSIONS

The concatenation and combined analysis of all plastid genome partitions and the construction of manually-curated, motif-based DNA sequence alignments are found to be instrumental in the recovery of well-supported relationships of the Gynoxoid group. We demonstrate that the correct assessment of homology in genome-level plastid sequence data sets is crucial for subsequent phylogeny reconstruction and that the manual post-processing of multiple sequence alignments improves the reliability of such reconstructions amid low genetic distances between taxa.

Calophyllaceae plastomes, their structure and insights in relationships within the clusioids

A complete chloroplast genome is not yet available for numerous species of plants. Among the groups that lack plastome information is the clusioid clade (Malpighiales), which includes five families: Bonnetiaceae, Calophyllaceae, Clusiaceae, Hypericaceae, and Podostemaceae. With around 2200 species, it has few published plastomes and most of them are from Podostemaceae. Here we assembled and compared six plastomes from members of the clusioids: five from Calophyllaceae (newly sequenced) and one from Clusiaceae. Putative regions for evolutionary studies were identified and the newly assembled chloroplasts were analyzed with other available chloroplasts for the group, focusing on Calophyllaceae. Our results mostly agree with recent studies which found a general conserved structure, except for the two Podostemaceae species that have a large inversion (trnK-UUU–rbcL) and lack one intron from ycf3. Within Calophyllaceae we observed a longer LSC and reduced IRs in Mahurea exstipulata, resulting in some genic rearrangement, and a short inversion (psbJ–psbE) in Kielmeyera coriacea. Phylogenetic analyses recovered the clusioids and the five families as monophyletic and revealed that conflicts in relationships reported in the literature for the group agree with nodes concentrating uninformative or conflicting gene trees. Our study brings new insights about clusioid plastome architecture and its evolution.

Disentangling Complex Inheritance Patterns of Plant Organellar Genomes: An Example From Carrot

Plant mitochondria and plastids display an array of inheritance patterns and varying levels of heteroplasmy, where individuals harbor more than 1 version of a mitochondrial or plastid genome. Organelle inheritance in plants has the potential to be quite complex and can vary with plant growth, development, and reproduction. Few studies have sought to investigate these complicated patterns of within-individual variation and inheritance using experimental crosses in plants. We carried out crosses in carrot, Daucus carota L. (Apiaceae), which has previously been shown to exhibit organellar heteroplasmy. We used mitochondrial and plastid markers to begin to disentangle the patterns of organellar inheritance and the fate of heteroplasmic variation, with special focus on cases where the mother displayed heteroplasmy. We also investigated heteroplasmy across the plant, assaying leaf samples at different development stages and ages. Mitochondrial and plastid paternal leakage was rare and offspring received remarkably similar heteroplasmic mixtures to their heteroplasmic mothers, indicating that heteroplasmy is maintained over the course of maternal inheritance. When offspring did differ from their mother, they were likely to exhibit a loss of the genetic variation that was present in their mother. Finally, we found that mitochondrial variation did not vary significantly over plant development, indicating that substantial vegetative sorting did not occur. Our study is one of the first to quantitatively investigate inheritance patterns and heteroplasmy in plants using controlled crosses, and we look forward to future studies making use of whole genome information to study the complex evolutionary dynamics of plant organellar genomes.

Divergence in the Aquilegia ecalcarata complex is correlated with geography and climate oscillations: Evidence from plastid genome data

Quaternary climate oscillations and geographical heterogeneity play important roles in determining species and genetic diversity distribution patterns, but how these factors affect the migration and differentiation of East Asian plants species at the population level remains poorly understood. The Aquilegia ecalcarata complex, a group that originated in the Late Tertiary and is widely distributed throughout East Asia, displays high genetic variation that is suitable for studying elaborate phylogeographic patterns and demographic history related to the impact of Quaternary climate and geography. We used plastid genome data from 322 individuals in 60 populations of the A. ecalcarata complex to thoroughly explore the impact of Quaternary climate oscillations and geography on the phylogeographic patterns and demographic history of the A. ecalcarata complex through a series of phylogenetic, divergence time estimation, and demographic history analyses. The dry, cold climate and frequent climate oscillations that occurred during the early Pleistocene and the Mid-Pleistocene transition led to the differentiation of the A. ecalcarata complex, which was isolated in various areas. Geographically, the A. ecalcarata complex can be divided into Eastern and Western Clades and five subclades, which conform to the divergence of the East Asian flora. Our results clearly show the impact of Quaternary climate and geography on evolutionary history at the population level. These findings promote the understanding of the relationship between plant genetic differentiation and climate and geographical factors of East Asia at the population level.

Combining molecular and geographical data to infer the phylogeny of Lamiales and its dispersal patterns in and out of the tropics

Lamiales is one of the most intractable orders of flowering plants, with several changes in family composition, and circumscription throughout history. The order is worldwide distributed, occurring in tropical forests and frozen habitats. In this study, a comprehensive phylogeny of Lamiales was reconstructed using DNA sequences. The tree was used to infer dispersal patterns, focusing on the tropics and extratropics. Molecular and species geographic data available from public repositories were combined to address both objectives. A total of 6,910 species, and 842 genera of Lamiales were sampled using the Python tool PyPHLAWD. The tree was inferred using RAxML, and recovered a monophyletic Lamiales. All 26 families were recovered as monophyletic with high support. The families Bignoniaceae, and Plantaginaceae are remarkable examples. The first emerged as monophyletic and included tribe Jacarandeae, while the later emerged as monophyletic in its sensu lato and included both the tribes Angelonieae, and Gratioleae. Distribution points for all species were retrieved from GBIF. After filtering, 1,136,425 records were retained. Species were coded as present in extratropical or tropical environments. The in and out of the tropics dispersal patterns were inferred using a maximum likelihood approach that identifies hidden rate changes. The model recovered higher rates of transition from extratropics to tropics, estimating two rates of state transitions. When ancestral states are considered, more discrete transitions from extratropics to tropics were observed. The extratropical state was also inferred for the crown node of Lamiales and old nested nodes, revealing a rare pattern of transitions to the tropics throughout the upper Cretaceous and Tertiary. A significant phylogenetic signal was recovered for the in and out of the tropics dispersal patterns, showing that state transitions are not frequent enough to erase the effect of tree structure on the data.