Multiple measures could alleviate long-branch attraction in phylogenomic reconstruction of Cupressoideae (Cupressaceae)

Is Over-parameterization a Problem for Profile Mixture Models?

Biochemical constraints on the admissible amino acids at specific sites in proteins lead to heterogeneity of the amino acid substitution process over sites in alignments. It is well known that phylogenetic models of protein sequence evolution that do not account for site heterogeneity are prone to long-branch attraction (LBA) artifacts. Profile mixture models were developed to model heterogeneity of preferred amino acids at sites via a finite distribution of site classes each with a distinct set of equilibrium amino acid frequencies. However, it is unknown whether the large number of parameters in such models associated with the many amino acid frequency vectors can adversely affect tree topology estimates because of over-parameterization. Here, we demonstrate theoretically that for long sequences, over-parameterization does not create problems for estimation with profile mixture models. Under mild conditions, tree, amino acid frequencies, and other model parameters converge to true values as sequence length increases, even when there are large numbers of components in the frequency profile distributions. Because large sample theory does not necessarily imply good behavior for shorter alignments we explore the performance of these models with short alignments simulated with tree topologies that are prone to LBA artifacts. We find that over-parameterization is not a problem for complex profile mixture models even when there are many amino acid frequency vectors. In fact, simple models with few site classes behave poorly. Interestingly, we also found that misspecification of the amino acid frequency vectors does not lead to increased LBA artifacts as long as the estimated cumulative distribution function of the amino acid frequencies at sites adequately approximates the true one. In contrast, misspecification of the amino acid exchangeability rates can severely negatively affect parameter estimation. Finally, we explore the effects of including in the profile mixture model an additional "F-class" representing the overall frequencies of amino acids in the data set. Surprisingly, the F-class does not help parameter estimation significantly and can decrease the probability of correct tree estimation, depending on the scenario, even though it tends to improve likelihood scores.

Mitochondrial genomes of four slug moths (Lepidoptera, Limacodidae): Genome description and phylogenetic implications

The family Limacodidae belongs to the superfamily Zygaenoidea, which includes 1672 species commonly referred to as slug moths. Limacodidae larvae are major pests for many economically important plant species and can cause human dermatitis. At present, the structure of the mitochondrial genome (mitogenome), phylogenetic position, and adaptive evolution of slug moths are poorly understood. Herein, the mitogenomes of Parasa lepida, Phlossa conjuncta, Thosea sinensis, and Setora sinensis were sequenced and compared with other available mitogenome sequences to better characterize the mitogenomic diversity and evolution of this moth family. The mitogenomes of P. lepida, P. conjuncta, T. sinensis, and S. sinensis were confirmed to be circular in structure with lengths of 15,575 bp, 15,553 bp, 15,535 bp, and 15,529 bp, respectively. The Limacodidae mitogenomes exhibited similar nucleotide composition, codon usage, RNA structure, and control region patterns, indicating the conservation of the mitogenome in the family Limacodidae. A sliding window, Ka/Ks, and genetic distance analyses revealed that the atp8 and nad6 genes exhibited the highest levels of variability and the most rapid evolutionary rates among the 13 protein-coding genes (PCGs) encoded in these Limacodidae mitogenomes, suggesting that they may offer value as candidate DNA markers. The phylogenetic analysis recovered the overall relationship as Tortricoidea + (Sesiidae + (Zygaenoidea + (Cossoidea/+Choreutoidea + (others)))). Within Zygaenoidea, Limacodidae was recovered as monophyletic, and the phylogenetic relationships were recovered as (Phaudidae + Zyganidae) + Limacodidae in all six phylogenetic trees. The analysis indicated that P. lepida, P. conjuncta, T. sinensis, and S. sinensis are members of the Limacodidae.

Phylogenetic studies of magnoliids: Advances and perspectives

Magnoliids are the largest flowering plant clades outside of the eudicots and monocots, which are distributed worldwide and have high economic, ornamental and ecological values. Eudicots, monocots and magnoliids are the three major clades of Mesangiospermae, and their phylogenetic relationship is one of the most interesting issues. In recent years, with the continuous accumulation of genomic information, the evolutionary status of magnoliids has become a hot spot in plant phylogenetic research. Although great efforts have been made to study the evolution of magnoliids using molecular data from several representative species such as nuclear genome, plastid genome, mitochondrial genome, and transcriptome, the results of current studies on the phylogenetic status of magnoliids are inconsistent. Here, we systematically describe the current understanding of the molecular research on magnoliid phylogeny and review the differences in the evolutionary state of magnoliids. Understanding the research approaches and limitations of magnoliid phylogeny can guide research strategies to further improve the study of the phylogenetic evolution of magnoliids.

Phylogeny and evolution of Cupressaceae: Updates on intergeneric relationships and new insights on ancient intergeneric hybridization

After the merger of the former Taxodiaceae and Cupressaceae s.s., currently the conifer family Cupressaceae (sensu lato) comprises seven subfamilies and 32 genera, most of which are important components of temperate and mountainous forests. With the exception of a recently published genus-level phylogeny of gymnosperms inferred from sequence analysis of 790 orthologs, previous phylogenetic studies of Cupressaceae were based mainly on morphological characters or a few molecular markers, and did not completely resolve the intergeneric relationships. In this study, we reconstructed a robust and well-resolved phylogeny of Cupressaceae represented by all 32 genera, using 1944 genes (Orthogroups) generated from transcriptome sequencing. Reticulate evolution analyses detected a possible ancient hybridization that occurred between ancestors of two subclades of Cupressoideae, including Microbiota-Platycladus-Tetraclinis (MPT) and Juniperus-Cupressus-Hesperocyparis-Callitropsis-Xanthocyparis (JCHCX), although both concatenation and coalescent trees are highly supported. Moreover, divergence time estimation and ancestral area reconstruction indicate that Cupressaceae very likely originated in Asia in the Triassic, and geographic isolation caused by continental separation drove the vicariant evolution of the two subfamilies Cupressoideae and Callitroideae in the northern and southern hemispheres, respectively. Evolutionary analyses of some morphological characters suggest that helically arranged linear-acicular leaves and imbricate bract-scale complexes represent ancestral states, and the shift from linear-acicular leaves to scale-like leaves was associated with the shift from helical to decussate arrangement. Our study sheds new light on phylogeny and evolutionary history of Cupressaceae, and strongly suggests that both dichotomous phylogenetic and reticulate evolution analyses be conducted in phylogenomic studies.

hlh-12, a gene that is necessary and sufficient to promote migration of gonadal regulatory cells in Caenorhabditis elegans, evolved within the Caenorhabditis clade

Specialized cells of the somatic gonad primordium of nematodes play important roles in the final form and function of the mature gonad. Caenorhabditis elegans hermaphrodites are somatic females that have a two-armed, U-shaped gonad that connects to the vulva at the midbody. The outgrowth of each gonad arm from the somatic gonad primordium is led by two female distal tip cells (fDTCs), while the anchor cell (AC) remains stationary and central to coordinate uterine and vulval development. The bHLH protein HLH-2 and its dimerization partners LIN-32 and HLH-12 had previously been shown to be required for fDTC specification. Here, we show that ectopic expression of both HLH-12 and LIN-32 in cells with AC potential transiently transforms them into fDTC-like cells. Furthermore, hlh-12 was known to be required for the fDTCs to sustain gonad arm outgrowth. Here, we show that ectopic expression of HLH-12 in the normally stationary AC causes displacement from its normal position and that displacement likely results from activation of the leader program of fDTCs because it requires genes necessary for gonad arm outgrowth. Thus, HLH-12 is both necessary and sufficient to promote gonadal regulatory cell migration. As differences in female gonadal morphology of different nematode species reflect differences in the fate or migratory properties of the fDTCs or of the AC, we hypothesized that evolutionary changes in the expression of hlh-12 may underlie the evolution of such morphological diversity. However, we were unable to identify an hlh-12 ortholog outside of Caenorhabditis. Instead, by performing a comprehensive phylogenetic analysis of all Class II bHLH proteins in multiple nematode species, we found that hlh-12 evolved within the Caenorhabditis clade, possibly by duplicative transposition of hlh-10. Our analysis suggests that control of gene regulatory hierarchies for gonadogenesis can be remarkably plastic during evolution without adverse phenotypic consequence.

Molecular evolution and SSRs analysis based on the chloroplast genome of Callitropsis funebris

Chloroplast genome sequences have been used to understand evolutionary events and to infer efficiently phylogenetic relationships. Callitropsis funebris (Cupressaceae) is an endemic species in China. Its phylogenetic position is controversial due to morphological characters similar to those of Cupressus, Callitropsis, and Chamaecyparis. This study used next-generation sequencing technology to sequence the complete chloroplast genome of Ca. funebris and then constructed the phylogenetic relationship between Ca. funebris and its related species based on a variety of data sets and methods. Simple sequence repeats (SSRs) and adaptive evolution analysis were also conducted. Our results showed that the monophyletic branch consisting of Ca. funebris and Cupressus tonkinensis is a sister to Cupressus, while Callitropsis is not monophyletic; Ca. nootkatensis and Ca. vietnamensis are nested in turn at the base of the monophyletic group Hesperocyparis. The statistical results of SSRs supported the closest relationship between Ca. funebris and Cupressus. By performing adaptive evolution analysis under the phylogenetic background of Cupressales, the Branch model detected three genes and the Site model detected 10 genes under positive selection; and the Branch-Site model uncovered that rpoA has experienced positive selection in the Ca. funebries branch. Molecular analysis from the chloroplast genome highly supported that Ca. funebris is at the base of Cupressus. Of note, SSR features were found to be able to shed some light on phylogenetic relationships. In short, this chloroplast genomic study has provided new insights into the phylogeny of Ca. funebris and revealed multiple chloroplast genes possibly undergoing adaptive evolution.

Testing the utility of alternative metrics of branch support to address the ancient evolutionary radiation of tunas, stromateoids, and allies (Teleostei: Pelagiaria)

The use of high-throughput sequencing technologies to produce genome-scale datasets was expected to settle some long-standing controversies across the Tree of Life, particularly in areas where short branches occur at deep timescales. Instead, these datasets have often yielded many well-supported but conflicting topologies, and highly variable gene-tree distributions. A variety of branch-support metrics beyond the nonparametric bootstrap are now available to assess how robust a phylogenetic hypothesis may be, as well as new methods to quantify gene-tree discordance. We applied multiple branch support metrics to an ancient group of marine fishes (Teleostei: Pelagiaria) whose interfamilial relationships have proven difficult to resolve due to a rapid accumulation of lineages very early in its history. We analyzed hundreds of loci including published UCE data and newly generated exonic data along with their flanking regions to represent all 16 extant families for more than 150 out of 284 valid species in the group. Branch support was lower for interfamilial relationships (except the SH-like aLRT and aBayes methods) regardless of the type of marker used. Several nodes that were highly supported with bootstrap had very low site and gene-tree concordance, revealing underlying conflict. Despite this conflict, we were able to identify four consistent interfamilial clades, each comprised of two or three families. Combining exons with their flanking regions also produced increased branch lengths in the deep branches of the pelagiarian tree. Our results demonstrate the limitations of employing current metrics of branch support and species-tree estimation when assessing the confidence of ancient evolutionary radiations and emphasize the necessity to embrace alternative measurements to explore phylogenetic uncertainty and discordance in phylogenomic datasets.

Plastome Phylogenomic and Biogeographical Study on Thuja (Cupressaceae)

Investigating the biogeographical disjunction of East Asian and North American flora is key to understanding the formation and dynamics of biodiversity in the Northern Hemisphere. The small Cupressaceae genus Thuja, comprising five species, exhibits a typical disjunct distribution in East Asia and North America. Owing to obscure relationships, the biogeographical history of the genus remains controversial. Here, complete plastomes were employed to investigate the plastome evolution, phylogenetic relationships, and biogeographic history of Thuja. All plastomes of Thuja share the same gene content arranged in the same order. The loss of an IR was evident in all Thuja plastomes, and the B-arrangement as previously recognized was detected. Phylogenomic analyses resolved two sister pairs, T. standishii-T. koraiensis and T. occidentalis-T. sutchuenensis, with T. plicata sister to T. occidentalis-T. sutchuenensis. Molecular dating and biogeographic results suggest the diversification of Thuja occurred in the Middle Miocene, and the ancestral area of extant species was located in northern East Asia. Incorporating the fossil record, we inferred that Thuja likely originated from the high-latitude areas of North America in the Paleocene with a second diversification center in northern East Asia. The current geographical distribution of Thuja was likely shaped by dispersal events attributed to the Bering Land Bridge in the Miocene and subsequent vicariance events accompanying climate cooling. The potential effect of extinction may have profound influence on the biogeographical history of Thuja.

The Origin and Evolution of Plastid Genome Downsizing in Southern Hemispheric Cypresses (Cupressaceae)

Plastome downsizing is rare in photosynthetic seed plants. However, a large-scale study of five cupressophyte families (conifers II) indicated that the plastomes of some Cupressaceous genera are notably reduced and compact. Here, we enriched taxon sampling in Cupressaceae by adding plastomes of ten previously unreported genera to determine the origin, evolution, and consequences of plastome reduction in this family. We discovered that plastome downsizing is specific to Callitroideae (a Southern Hemispheric subfamily). Their plastomes are the smallest, encode the fewest plastid genes, and contain the fewest GC-end codons among Cupressaceae. We show that repeated tRNA losses and shrinkage of intergenic spacers together contributed to the plastome downsizing in Callitroideae. Moreover, our absolute nucleotide substitution rate analyses suggest relaxed functional constraints in translation-related plastid genes (clpP, infA, rpl, and rps), but not in photosynthesis- or transcription-related ones, of Callitris (the most diverse genus in Callitroideae). We hypothesize that the small and low-GC plastomes of Callitroideae emerged ca. 112-75 million years ago as an adaptation to increased competition with angiosperms on the Gondwana supercontinent. Our findings highlight Callitroideae as another case of plastome downsizing in photosynthetic seed plant lineages.

Chloroplast genome of an extremely endangered conifer Thuja sutchuenensis Franch.: gene organization, comparative and phylogenetic analysis

Thuja sutchuenensis is a critically endangered tertiary relict species of Cupressaceae from southwestern China. We sequenced the complete chloroplast (cp) genome of T. sutchuenensis, showing the genome content of 129,776 bp, 118 unique genes including 82 unique protein-coding genes, 32 tRNA genes, and 4 rRNA genes. The genome structures, gene order, and GC content are similar to other typical gymnosperm cp genomes. Thirty-eight simple sequence repeats were identified in the T. sutchuenensis cp genome. We also found an apparent inversion between trnT and psbK between genera Thuja and Thujopsis. In addition, positive selection signals were detected in seven genes with high Ka/Ks ratios. The reconstructed phylogeny based on locally collinear blocks of cp genomes among 21 gymnosperms species is similar to previous inferences. We also inferred a Late-Miocene divergence between T. sutchuenensis and T. standishii, according to the dating of ~ 11.05 Mya by cp genomes. These results will be helpful for future studies of Cupressaceae phylogeny as well as studies in population genetics, systematics, and cp genetic engineering.

An empirical assessment of a single family-wide hybrid capture locus set at multiple evolutionary timescales in Asteraceae

PREMISE

Hybrid capture with high-throughput sequencing (Hyb-Seq) is a powerful tool for evolutionary studies. The applicability of an Asteraceae family-specific Hyb-Seq probe set and the outcomes of different phylogenetic analyses are investigated here.

METHODS

Hyb-Seq data from 112 Asteraceae samples were organized into groups at different taxonomic levels (tribe, genus, and species). For each group, data sets of non-paralogous loci were built and proportions of parsimony informative characters estimated. The impacts of analyzing alternative data sets, removing long branches, and type of analysis on tree resolution and inferred topologies were investigated in tribe Cichorieae.

RESULTS

Alignments of the Asteraceae family-wide Hyb-Seq locus set were parsimony informative at all taxonomic levels. Levels of resolution and topologies inferred at shallower nodes differed depending on the locus data set and the type of analysis, and were affected by the presence of long branches.

DISCUSSION

The approach used to build a Hyb-Seq locus data set influenced resolution and topologies inferred in phylogenetic analyses. Removal of long branches improved the reliability of topological inferences in maximum likelihood analyses. The Astereaceae Hyb-Seq probe set is applicable at multiple taxonomic depths, which demonstrates that probe sets do not necessarily need to be lineage-specific.

A Critical Appraisal of the Placement of Xiphosura (Chelicerata) with Account of Known Sources of Phylogenetic Error

Horseshoe crabs (Xiphosura) are traditionally regarded as sister group to the clade of terrestrial chelicerates (Arachnida). This hypothesis has been challenged by recent phylogenomic analyses, but the non-monophyly of Arachnida has consistently been disregarded as artifactual. We re-evaluated the placement of Xiphosura among chelicerates using the most complete phylogenetic data set to date, expanding outgroup sampling, and including data from whole genome sequencing projects. In spite of uncertainty in the placement of some arachnid clades, all analyses show Xiphosura consistently nested within Arachnida as the sister group to Ricinulei (hooded tick spiders). It is apparent that the radiation of arachnids is an old one and occurred over a brief period of time, resulting in several consecutive short internodes, and thus is a potential case for the confounding effects of incomplete lineage sorting (ILS). We simulated coalescent gene trees to explore the effects of increasing levels of ILS on the placement of horseshoe crabs. In addition, common sources of systematic error were evaluated, as well as the effects of fast-evolving partitions and the dynamics of problematic long branch orders. Our results indicated that the placement of horseshoe crabs cannot be explained by missing data, compositional biases, saturation, or ILS. Interrogation of the phylogenetic signal showed that the majority of loci favor the derived placement of Xiphosura over a monophyletic Arachnida. Our analyses support the inference that horseshoe crabs represent a group of aquatic arachnids, comparable to aquatic mites, breaking a long-standing paradigm in chelicerate evolution and altering previous interpretations of the ancestral transition to the terrestrial habitat. Future studies testing chelicerate relationships should approach the task with a sampling strategy where the monophyly of Arachnida is not held as the premise.

Statistical Inconsistency of Maximum Parsimony for k-Tuple-Site Data

One of the main aims of phylogenetics is to reconstruct the "Tree of Life." In this respect, different methods and criteria are used to analyze DNA sequences of different species and to compare them in order to derive the evolutionary relationships of these species. Maximum parsimony is one such criterion for tree reconstruction, and it is the one which we will use in this paper. However, it is well known that tree reconstruction methods can lead to wrong relationship estimates. One typical problem of maximum parsimony is long branch attraction, which can lead to statistical inconsistency. In this work, we will consider a blockwise approach to alignment analysis, namely the so-called k-tuple analyses. For four taxa, it has already been shown that k-tuple-based analyses are statistically inconsistent if and only if the standard character-based (site-based) analyses are statistically inconsistent. So, in the four-taxon case, going from individual sites to k-tuples does not lead to any improvement. However, real biological analyses often consider more than only four taxa. Therefore, we analyze the case of five taxa for 2- and 3-tuple-site data and consider alphabets with two and four elements. We show that the equivalence of single-site data and k-tuple-site data then no longer holds. Even so, we can show that maximum parsimony is statistically inconsistent for k-tuple-site data and five taxa.

A transcriptome-based resolution for a key taxonomic controversy in Cupressaceae

BACKGROUND AND AIMS

Rapid evolutionary divergence and reticulate evolution may result in phylogenetic relationships that are difficult to resolve using small nucleotide sequence data sets. Next-generation sequencing methods can generate larger data sets that are better suited to solving these puzzles. One major and long-standing controversy in conifers concerns generic relationships within the subfamily Cupressoideae (105 species, approx. 1/6 of all conifers) of Cupressaceae, in particular the relationship between Juniperus, Cupressus and the Hesperocyparis-Callitropsis-Xanthocyparis (HCX) clade. Here we attempt to resolve this question using transcriptome-derived data.

METHODS

Transcriptome sequences of 20 species from Cupressoideae were collected. Using MarkerMiner, single-copy nuclear (SCN) genes were extracted. These were applied to estimate phylogenies based on concatenated data, species trees and a phylogenetic network. We further examined the effect of alternative backbone topologies on downstream analyses, including biogeographic inference and dating analysis.

RESULTS

Based on the 73 SCN genes (>200 000 bp total alignment length) we considered, all tree-building methods lent strong support for the relationship (HCX, (Juniperus, Cupressus)); however, strongly supported conflicts among individual gene trees were also detected. Molecular dating suggests that these three lineages shared a most recent common ancestor approx. 60 million years ago (Mya), and that Juniperus and Cupressus diverged about 56 Mya. Ancestral area reconstructions (AARs) suggest an Asian origin for the entire clade, with subsequent dispersal to North America, Europe and Africa.

CONCLUSIONS

Our analysis of SCN genes resolves a controversial phylogenetic relationship in the Cupressoideae, a major clade of conifers, and suggests that rapid evolutionary divergence and incomplete lineage sorting probably acted together as the source for conflicting phylogenetic inferences between gene trees and between our robust results and recently published studies. Our updated backbone topology has not substantially altered molecular dating estimates relative to previous studies; however, application of the latest AAR approaches has yielded a clearer picture of the biogeographic history of Cupressoideae.

A Malvaceae mystery: A mallow maelstrom of genome multiplications and maybe misleading methods?

Previous research suggests that Gossypium has undergone a 5- to 6-fold multiplication following its divergence from Theobroma. However, the number of events, or where they occurred in the Malvaceae phylogeny remains unknown. We analyzed transcriptomic and genomic data from representatives of eight of the nine Malvaceae subfamilies. Phylogenetic analysis of nuclear data placed Dombeya (Dombeyoideae) as sister to the rest of Malvadendrina clade, but the plastid DNA tree strongly supported Durio (Helicteroideae) in this position. Intraspecific Ks plots indicated that all sampled taxa, except Theobroma (Byttnerioideae), Corchorus (Grewioideae), and Dombeya (Dombeyoideae), have experienced whole genome multiplications (WGMs). Quartet analysis suggested WGMs were shared by Malvoideae-Bombacoideae and Sterculioideae-Tilioideae, but did not resolve whether these are shared with each other or Helicteroideae (Durio). Gene tree reconciliation and Bayesian concordance analysis suggested a complex history. Alternative hypotheses are suggested, each involving two independent autotetraploid and one allopolyploid event. They differ in that one entails an allopolyploid origin for the Durio lineage, whereas the other invokes an allopolyploid origin for Malvoideae-Bombacoideae. We highlight the need for more genomic information in the Malvaceae and improved methods to resolve complex evolutionary histories that may include allopolyploidy, incomplete lineage sorting, and variable rates of gene and genome evolution.

Phylogenomic evidence for ancient recombination between plastid genomes of the Cupressus-Juniperus-Xanthocyparis complex (Cupressaceae)

BACKGROUND

Phylogenetic relationships among Eastern Hemisphere cypresses, Western Hemisphere cypresses, junipers, and their closest relatives are controversial, and generic delimitations have been in flux for the past decade. To address relationships and attempt to produce a more robust classification, we sequenced 11 new plastid genomes (plastomes) from the five variously described genera in this complex (Callitropsis, Cupressus, Hesperocyparis, Juniperus, and Xanthocyparis) and compared them with additional plastomes from diverse members of Cupressaceae.

RESULTS

Phylogenetic analysis of protein-coding genes recovered a topology in which Juniperus is sister to Cupressus, whereas a tree based on whole plastomes indicated that the Callitropsis-Hesperocyparis-Xanthocyparis (CaHX) clade is sister to Cupressus. A sliding window analysis of site-specific phylogenetic support identified a ~ 15 kb region, spanning the genes ycf1 and ycf2, which harbored an anomalous signal relative to the rest of the genome. After excluding these genes, trees based on the remainder of the genes and genome consistently recovered a topology grouping the CaHX clade and Cupressus with strong bootstrap support. In contrast, trees based on the ycf1 and ycf2 region strongly supported a sister relationship between Cupressus and Juniperus.

CONCLUSIONS

These results demonstrate that standard phylogenomic analyses can result in strongly supported but conflicting trees. We suggest that the conflicting plastomic signals result from an ancient introgression event involving ycf1 and ycf2 that occurred in an ancestor of this species complex. The introgression event was facilitated by plastomic recombination in an ancestral heteroplasmic individual carrying distinct plastid haplotypes, offering further evidence that recombination occurs between plastomes. Finally, we provide strong support for previous proposals to recognize five genera in this species complex: Callitropsis, Cupressus, Hesperocyparis, Juniperus, and Xanthocyparis.

Phylogenomic analysis on the exceptionally diverse fish clade Gobioidei (Actinopterygii: Gobiiformes) and data-filtering based on molecular clocklikeness

The use of genome-scale data to infer phylogenetic relationships has gained in popularity in recent years due to the progress made in target-gene capture and sequencing techniques. Data filtering, the approach of excluding data inconsistent with the model from analyses, presumably could alleviate problems caused by systematic errors in phylogenetic inference. Different data filtering criteria, such as those based on evolutionary rate and molecular clocklikeness as well as others have been proposed for selecting useful phylogenetic markers, yet few studies have tested these criteria using phylogenomic data. We developed a novel set of single-copy nuclear coding markers to capture thousands of target genes in gobioid fishes, a species-rich lineages of vertebrates, and tested the effects of data-filtering methods based on substitution rate and molecular clocklikeness while attempting to control for the compounding effects of missing data and variation in locus length. We found that molecular clocklikeness was a better predictor than overall substitution rate for phylogenetic usefulness of molecular markers in our study. In addition, when the 100 best ranked loci for our predictors were concatenated and analyzed using maximum likelihood, or combined in a coalescent-based species-tree analysis, the resulting trees showed a well-resolved topology of Gobioidei that mostly agrees with previous studies. However, trees generated from the 100 least clocklike frequently recovered conflicting, and in some cases clearly erroneous topologies with strong support, thus indicating strong systematic biases in those datasets. Collectively these results suggest that data filtering has the potential improve the performance of phylogenetic inference when using both a concatenation approach as well as methods that rely on input from individual gene trees (i.e. coalescent species-tree approaches), which may be preferred in scenarios where incomplete lineage sorting is likely to be an issue.

Plastome characteristics of Cannabaceae

Cannabaceae is an economically important family that includes ten genera and ca. 117 accepted species. To explore the structure and size variation of their plastomes, we sequenced ten plastomes representing all ten genera of Cannabaceae. Each plastome possessed the typical angiosperm quadripartite structure and contained a total of 128 genes. The Inverted Repeat (IR) regions in five plastomes had experienced small expansions (330-983 bp) into the Large Single-Copy (LSC) region. The plastome of Chaetachme aristata has experienced a 942-bp IR contraction and lost rpl22 and rps19 in its IRs. The substitution rates of rps19 and rpl22 decreased after they shifted from the LSC to IR. A 270-bp inversion was detected in the Parasponia rugosa plastome, which might have been mediated by 18-bp inverted repeats. Repeat sequences, simple sequence repeats, and nucleotide substitution rates varied among these plastomes. Molecular markers with more than 13% variable sites and 5% parsimony-informative sites were identified, which may be useful for further phylogenetic analysis and species identification. Our results show strong support for a sister relationship between Gironniera and Lozanell (BS = 100). Celtis, Cannabis-Humulus, Chaetachme-Pteroceltis, and Trema-Parasponia formed a strongly supported clade, and their relationships were well resolved with strong support (BS = 100). The availability of these ten plastomes provides valuable genetic information for accurately identifying species, clarifying taxonomy and reconstructing the intergeneric phylogeny of Cannabaceae.

Practical considerations for plant phylogenomics

The past decade has seen a major breakthrough in our ability to easily and inexpensively sequence genome-scale data from diverse lineages. The development of high-throughput sequencing and long-read technologies has ushered in the era of phylogenomics, where hundreds to thousands of nuclear genes and whole organellar genomes are routinely used to reconstruct evolutionary relationships. As a result, understanding which options are best suited for a particular set of questions can be difficult, especially for those just starting in the field. Here, we review the most recent advances in plant phylogenomic methods and make recommendations for project-dependent best practices and considerations. We focus on the costs and benefits of different approaches in regard to the information they provide researchers and the questions they can address. We also highlight unique challenges and opportunities in plant systems, such as polyploidy, reticulate evolution, and the use of herbarium materials, identifying optimal methodologies for each. Finally, we draw attention to lingering challenges in the field of plant phylogenomics, such as reusability of data sets, and look at some up-and-coming technologies that may help propel the field even further.

Plastid Genome Evolution in the Early-Diverging Legume Subfamily Cercidoideae (Fabaceae)

The subfamily Cercidoideae is an early-branching legume lineage, which consists of 13 genera distributed in the tropical and warm temperate Northern Hemisphere. A previous study detected two plastid genomic variations in this subfamily, but the limited taxon sampling left the overall plastid genome (plastome) diversification across the subfamily unaddressed, and phylogenetic relationships within this clade remained unresolved. Here, we assembled eight plastomes from seven Cercidoideae genera and conducted phylogenomic-comparative analyses in a broad evolutionary framework across legumes. The plastomes of Cercidoideae all exhibited a typical quadripartite structure with a conserved gene content typical of most angiosperm plastomes. Plastome size ranged from 151,705 to 165,416 bp, mainly due to the expansion and contraction of inverted repeat (IR) regions. The order of genes varied due to the occurrence of several inversions. In Tylosema species, a plastome with a 29-bp IR-mediated inversion was found to coexist with a canonical-type plastome, and the abundance of the two arrangements of isomeric molecules differed between individuals. Complete plastome data were much more efficient at resolving intergeneric relationships of Cercidoideae than the previously used selection of only a few plastid or nuclear loci. In sum, our study revealed novel insights into the structural diversification of plastomes in an early-branching legume lineage, and, thus, into the evolutionary trajectories of legume plastomes in general.