Phylogenetic signal detection from an ancient rapid radiation: Effects of noise reduction, long-branch attraction, and model selection in crown clade Apocynaceae

African Mountain Thistles: Three New Genera in the Carduus-Cirsium Group

The floras on the highest mountains in tropical eastern Africa are among the most unique floras in the world. Despite the exceptionally high concentration of endemic species, these floras remain understudied from an evolutionary point of view. In this study, we focus on the Carduus-Cirsium group (subtribe Carduinae) to unravel the evolutionary relationships of the species endemic to the tropical Afromontane and Afroalpine floras, aiming to improve the systematics of the group. We applied the Hyb-Seq approach using the Compositae1061 probe set on 190 samples (159 species), encompassing representatives of all genera of Carduinae. We used two recently developed pipelines that enabled the processing of raw sequence reads, identification of paralogous sequences and segregation into orthologous alignments. After the implementation of a missing data filter, we retained sequences from 986 nuclear loci and 177 plastid regions. Phylogenomic analyses were conducted using both concatenated and summary-coalescence methods. The resulting phylogenies were highly resolved and revealed three distinct evolutionary lineages consisting of the African species traditionally referred to as Carduus and Cirsium. Consequently, we propose the three new genera Afrocarduus, Afrocirsium and Nuriaea; the latter did notably not belong to the Carduus-Cirsium group. We detected some incongruences between the phylogenies based on concatenation vs. coalescence and on nuclear vs. plastid datasets, likely attributable to incomplete lineage sorting and/or hybridization.

Evolution of 101 Apocynaceae plastomes and phylogenetic implications

Apocynaceae are one of the ten species-richest angiosperm families. However, the backbone phylogeny of the family is yet less well supported, and the evolution of plastome structure has not been thoroughly studied for the whole family. Herein, a total of 101 complete plastomes including 35 newly sequenced, 24 reassembled from public raw data and the rest from the NCBI GenBank database, representing 26 of 27 tribes of Apocynaceae, were used for comparative plastome analysis. Phylogenetic analyses were conducted using a combined plastid data matrix of 77 protein-coding genes from 162 taxa, encompassing all tribes and 41 of 49 subtribes of Apocynaceae. Plastome lengths ranged from 150,897 bp in Apocynum venetum to 178,616 bp in Hoya exilis. Six types of boundaries between the inverted repeat (IR) regions and single copy (SC) regions were identified. Different sizes of IR expansion were found in three lineages, including Alyxieae, Ceropegieae and Marsdenieae, suggesting multiple expansion events of the IRs over the SC regions in Apocynaceae. The IR regions of Marsdenieae evolved in two ways: expansion towards the large single copy (LSC) region in Lygisma + Stephanotis + Ruehssia + Gymnema (Cosmopolitan clade), and expansion towards both LSC and small single copy (SSC) region in Dischidia-Hoya alliance and Marsdenia (Asia-Pacific clade). Six coding genes and five non-coding regions were identified as highly variable, including accD, ccsA-ndhD, clpP, matK, ndhF, ndhG-ndhI, trnG(GCC)-trnfM(CAU), trnH(GUG)-psbA, trnY(GUA)-trnE(UUC), ycf1, and ycf2. Maximum likelihood and Bayesian phylogenetic analyses resulted in nearly identical tree topologies and produced a well-resolved backbone comprising 15 consecutive dichotomies that subdivided Apocynaceae into 15 clades. The subfamily Periplocoideae were embedded in the Apocynoid grade and were sister to the Echiteae-Odontadenieae-Mesechiteae clade with high support values. Three tribes (Melodineae, Vinceae, and Willughbeieae), the subtribe Amphineuriinae, and four genera (Beaumontia, Ceropegia, Hoya, and Stephanotis) were not resolved as monophyletic. Our work sheds light on the backbone phylogenetic relationships in the family Apocynaceae and offers insights into the evolution of Apocynaceae plastomes using the most densely sampled plastome dataset to date.

The complete chloroplast genome of Holarrhena pubescens and its phylogenetic analysis

Holarrhena pubescens Wall. ex G. Don, 1837 is an important medicinal plant belonging to the Holarrhena genus in the Apocynaceae family. In this study, the complete chloroplast (cp) genome sequence of H. pubescens was sequenced using the Illumina NovaSeq platform. The cp genome of H. pubescens was 160,108 bp in length with 37.21% overall GC content. The cp genome of H. pubescens containing a large single-copy region (LSC, 88,685 bp), a small single-copy region (SSC, 18,671 bp), and a pair of inverted repeat regions (SSC, 26,376 bp). The cp genome encoded 129 genes, including 84 protein-coding genes, 37 tRNA genes, and eight rRNA genes. Phylogenetic analysis based on complete protein coding genes sequences revealed that H. pubescens was closest to Beaumontia murtonii.

Incomplete lineage sorting and reticulate evolution mask species relationships in Brunelliaceae, an Andean family with rapid, recent diversification

PREMISE

To date, phylogenetic relationships within the monogeneric Brunelliaceae have been based on morphological evidence, which does not provide sufficient phylogenetic resolution. Here we use target-enriched nuclear data to improve our understanding of phylogenetic relationships in the family.

METHODS

We used the Angiosperms353 toolkit for targeted recovery of exonic regions and supercontigs (exons + introns) from low copy nuclear genes from 53 of 70 species in Brunellia, and several outgroup taxa. We removed loci that indicated biased inference of relationships and applied concatenated and coalescent methods to infer Brunellia phylogeny. We identified conflicts among gene trees that may reflect hybridization or incomplete lineage sorting events and assessed their impact on phylogenetic inference. Finally, we performed ancestral-state reconstructions of morphological traits and assessed the homology of character states used to define sections and subsections in Brunellia.

RESULTS

Brunellia comprises two major clades and several subclades. Most of these clades/subclades do not correspond to previous infrageneric taxa. There is high topological incongruence among the subclades across analyses.

CONCLUSIONS

Phylogenetic reconstructions point to rapid species diversification in Brunelliaceae, reflected in very short branches between successive species splits. The removal of putatively biased loci slightly improves phylogenetic support for individual clades. Reticulate evolution due to hybridization and/or incomplete lineage sorting likely both contribute to gene-tree discordance. Morphological characters used to define taxa in current classification schemes are homoplastic in the ancestral character-state reconstructions. While target enrichment data allows us to broaden our understanding of diversification in Brunellia, the relationships among subclades remain incompletely understood.

An NGS-Based Phylogeny of Orthotricheae (Orthotrichaceae, Bryophyta) With the Proposal of the New Genus Rehubryum From Zealandia

Phylogenomic data increase the possibilities of resolving the evolutionary and systematic relationships among taxa. This is especially valuable in groups with few and homoplasious morphological characters, in which systematic and taxonomical delimitations have been traditionally difficult. Such is the case of several lineages within Bryophyta, like Orthotrichaceae, the second most diverse family of mosses. Members of tribe Orthotricheae are common in temperate and cold regions, as well as in high tropical mountains. In extratropical areas, they represent one of the main components of epiphytic communities, both in dry and oceanic or hyperoceanic conditions. The epiphytic environment is considered a hostile one for plant development, mainly due to its low capacity of moisture retention. Thus, the diversification of the Orthotrichaceae in this environment could be seen as striking. Over the last two decades, great taxonomic and systematic progresses have led to a rearrangement at the generic level in this tribe, providing a new framework to link environment to patterns of diversification. Here, we use nuclear loci targeted with the GoFlag 408 enrichment probe set to generate a well-sampled phylogeny with well-supported suprageneric taxa and increasing the phylogenetic resolution within the two recognized subtribes. Specifically, we show that several genera with Ulota-like morphology jointly constitute an independent lineage. Within this lineage, the recently described Atlantichella from Macaronesia and Western Europe appears as the sister group of Ulota bellii from Zealandia. This latter species is here segregated in the new genus Rehubryum. Assessment of the ecological and biogeographical affinities of the species within the phylogenetic framework suggests that niche adaptation (including climate and substrate) may be a key evolutionary driver that shaped the high diversification of Orthotricheae.

A mitochondrial genome phylogeny of voles and lemmings (Rodentia: Arvicolinae): Evolutionary and taxonomic implications

Arvicolinae is one of the most impressive placental radiations with over 150 extant and numerous extinct species that emerged since the Miocene in the Northern Hemisphere. The phylogeny of Arvicolinae has been studied intensively for several decades using morphological and genetic methods. Here, we sequenced 30 new mitochondrial genomes to better understand the evolutionary relationships among the major tribes and genera within the subfamily. The phylogenetic and molecular dating analyses based on 11,391 bp concatenated alignment of protein-coding mitochondrial genes confirmed the monophyly of the subfamily. While Bayesian analysis provided a high resolution across the entire tree, Maximum Likelihood tree reconstruction showed weak support for the ordering of divergence and interrelationships of tribal level taxa within the most ancient radiation. Both the interrelationships among tribes Lagurini, Ellobiusini and Arvicolini, comprising the largest radiation and the position of the genus Dinaromys within it also remained unresolved. For the first time complex relationships between genus level taxa within the species-rich tribe Arvicolini received full resolution. Particularly Lemmiscus was robustly placed as sister to the snow voles Chionomys in the tribe Arvicolini in contrast with a long-held belief of its affinity with Lagurini. Molecular dating of the origin of Arvicolinae and early divergences obtained from the mitogenome data were consistent with fossil records. The mtDNA estimates for putative ancestors of the most genera within Arvicolini appeared to be much older than it was previously proposed in paleontological studies.

Sweet or salty? The origin of freshwater gastrotrichs (Gastrotricha, Chaetonotida) revealed by molecular phylogenetic analysis

Chaetonotidae is the most diverse and widely distributed family of the order Chaetonotida (Gastrotricha) and includes both marine and freshwater species. Although the family is regarded as a sister taxon to the exclusively marine Xenotrichulidae, the type of environment, marine or freshwater, where Chaetonotidae originated is still not known. Here, we reconstructed the phylogeny of the family based on molecular sequence data and mapped both morphological and ecological characters to determine the ancestral environment of the first members of the family. Our results revealed that the freshwater genus Bifidochaetus is the earliest branching lineage in the paraphyletic Chaetonotidae (encompassing Dasydytidae and Neogosseidae). Moreover, we reconstructed Lepidochaetus-Cephalionotus clade as a monophyletic sister group to the remaining chaetonotids, which supports Kisielewski's morphological based hypothesis concerning undifferentiated type of body scales as a most primary character in Chaetonotidae. We also found that reversals to marine habitats occurred independently in different Chaetonotidae lineages, thus marine species in the genera Heterolepidoderma, Halichaetonotus, Aspidiophorus and subgenera Chaetonotus (Schizochaetonotus) or Chaetonotus (Marinochaetus) should be assumed as having secondarily invaded the marine environment. Character mapping revealed a series of synapomorphies that define the clade that includes Chaetonotidae (with Dasydytidae and Neogosseidae), the most important of which may be those linked to reproduction.

Evolution of Dispersal, Habit, and Pollination in Africa Pushed Apocynaceae Diversification After the Eocene-Oligocene Climate Transition

Apocynaceae (the dogbane and milkweed family) is one of the ten largest flowering plant families, with approximately 5,350 species and diverse morphology and ecology, ranging from large trees and lianas that are emblematic of tropical rainforests, to herbs in temperate grasslands, to succulents in dry, open landscapes, and to vines in a wide variety of habitats. Despite a specialized and conservative basic floral architecture, Apocynaceae are hyperdiverse in flower size, corolla shape, and especially derived floral morphological features. These are mainly associated with the development of corolline and/or staminal coronas and a spectrum of integration of floral structures culminating with the formation of a gynostegium and pollinaria—specialized pollen dispersal units. To date, no detailed analysis has been conducted to estimate the origin and diversification of this lineage in space and time. Here, we use the most comprehensive time-calibrated phylogeny of Apocynaceae, which includes approximately 20% of the species covering all major lineages, and information on species number and distributions obtained from the most up-to-date monograph of the family to investigate the biogeographical history of the lineage and its diversification dynamics. South America, Africa, and Southeast Asia (potentially including Oceania), were recovered as the most likely ancestral area of extant Apocynaceae diversity; this tropical climatic belt in the equatorial region retained the oldest extant lineages and these three tropical regions likely represent museums of the family. Africa was confirmed as the cradle of pollinia-bearing lineages and the main source of Apocynaceae intercontinental dispersals. We detected 12 shifts toward accelerated species diversification, of which 11 were in the APSA clade (apocynoids, Periplocoideae, Secamonoideae, and Asclepiadoideae), eight of these in the pollinia-bearing lineages and six within Asclepiadoideae. Wind-dispersed comose seeds, climbing growth form, and pollinia appeared sequentially within the APSA clade and probably work synergistically in the occupation of drier and cooler habitats. Overall, we hypothesize that temporal patterns in diversification of Apocynaceae was mainly shaped by a sequence of morphological innovations that conferred higher capacity to disperse and establish in seasonal, unstable, and open habitats, which have expanded since the Eocene-Oligocene climate transition.

Taxon-specific or universal? Using target capture to study the evolutionary history of rapid radiations

Target capture has emerged as an important tool for phylogenetics and population genetics in nonmodel taxa. Whereas developing taxon‐specific capture probes requires sustained efforts, available universal kits may have a lower power to reconstruct relationships at shallow phylogenetic scales and within rapidly radiating clades. We present here a newly developed target capture set for Bromeliaceae, a large and ecologically diverse plant family with highly variable diversification rates. The set targets 1776 coding regions, including genes putatively involved in key innovations, with the aim to empower testing of a wide range of evolutionary hypotheses. We compare the relative power of this taxon‐specific set, Bromeliad1776, to the universal Angiosperms353 kit. The taxon‐specific set results in higher enrichment success across the entire family; however, the overall performance of both kits to reconstruct phylogenetic trees is relatively comparable, highlighting the vast potential of universal kits for resolving evolutionary relationships. For more detailed phylogenetic or population genetic analyses, for example the exploration of gene tree concordance, nucleotide diversity or population structure, the taxon‐specific capture set presents clear benefits. We discuss the potential lessons that this comparative study provides for future phylogenetic and population genetic investigations, in particular for the study of evolutionary radiations.

Highly Diverse Shrub Willows (Salix L.) Share Highly Similar Plastomes

Plastome phylogenomics is used in a broad range of studies where single markers do not bear enough information. Phylogenetic reconstruction in the genus Salix is difficult due to the lack of informative characters and reticulate evolution. Here, we use a genome skimming approach to reconstruct 41 complete plastomes of 32 Eurasian and North American Salix species representing different lineages, different ploidy levels, and separate geographic regions. We combined our plastomes with published data from Genbank to build a comprehensive phylogeny of 61 samples (50 species) using RAxML (Randomized Axelerated Maximum Likelihood). Additionally, haplotype networks for two observed subclades were calculated, and 72 genes were tested to be under selection. The results revealed a highly conserved structure of the observed plastomes. Within the genus, we observed a variation of 1.68%, most of which separated subg. Salix from the subgeneric Chamaetia/Vetrix clade. Our data generally confirm previous plastid phylogenies, however, within Chamaetia/Vetrix phylogenetic results represented neither taxonomical classifications nor geographical regions. Non-coding DNA regions were responsible for most of the observed variation within subclades and 5.6% of the analyzed genes showed signals of diversifying selection. A comparison of nuclear restriction site associated DNA (RAD) sequencing and plastome data on a subset of 10 species showed discrepancies in topology and resolution. We assume that a combination of (i) a very low mutation rate due to efficient mechanisms preventing mutagenesis, (ii) reticulate evolution, including ancient and ongoing hybridization, and (iii) homoplasy has shaped plastome evolution in willows.

Plastome evolution and organisation in the Hoya group (Apocynaceae)

The genus Hoya is highly diverse and many of its species are popular ornamental plants. However, the relationships between Hoya and related genera (the Hoya group) are not fully resolved. In this study, we report 20 newly sequenced plastomes of species in the Hoya group. The complete plastomes vary in length from 175,405 to 178,525 bp while the LSCs vary from 90,248 to 92,364 bp and the complete SSCs vary from 2,285 to 2,304 bp, making the SSC in the Hoya group one of the shortest known in the angiosperms. The plastome structure in the Hoya group is characterised by a massive increase in the size of the inverted repeats as compared to the outgroups. In all ingroup species, the IR/SSC boundary moved from ycf1 to ndhF while this was not observed in outgroup taxa, making it a synapomorphy for the Hoya group. We have also assembled the mitogenome of Hoya lithophytica, which, at 718,734 bp, is the longest reported in the family. The phylogenetic analysis using exons from 42 taxa in the Hoya group and three outgoups confirms that the earliest divergent genus in the Hoya group is Papuahoya, followed by Dischidia. The relationship between Dischidia and the clade which includes all Hoya and Oreosparte taxa, is not fully supported. Oreosparte is nested in Hoya making it paraphyletic unless Clemensiella is recognised as a separate genus.

Hundreds of nuclear and plastid loci yield novel insights into orchid relationships

PREMISE

The inference of evolutionary relationships in the species-rich family Orchidaceae has hitherto relied heavily on plastid DNA sequences and limited taxon sampling. Previous studies have provided a robust plastid phylogenetic framework, which was used to classify orchids and investigate the drivers of orchid diversification. However, the extent to which phylogenetic inference based on the plastid genome is congruent with the nuclear genome has been only poorly assessed.

METHODS

We inferred higher-level phylogenetic relationships of orchids based on likelihood and ASTRAL analyses of 294 low-copy nuclear genes sequenced using the Angiosperms353 universal probe set for 75 species (representing 69 genera, 16 tribes, 24 subtribes) and a concatenated analysis of 78 plastid genes for 264 species (117 genera, 18 tribes, 28 subtribes). We compared phylogenetic informativeness and support for the nuclear and plastid phylogenetic hypotheses.

RESULTS

Phylogenetic inference using nuclear data sets provides well-supported orchid relationships that are highly congruent between analyses. Comparisons of nuclear gene trees and a plastid supermatrix tree showed that the trees are mostly congruent, but revealed instances of strongly supported phylogenetic incongruence in both shallow and deep time. The phylogenetic informativeness of individual Angiosperms353 genes is in general better than that of most plastid genes.

CONCLUSIONS

Our study provides the first robust nuclear phylogenomic framework for Orchidaceae and an assessment of intragenomic nuclear discordance, plastid-nuclear tree incongruence, and phylogenetic informativeness across the family. Our results also demonstrate what has long been known but rarely thoroughly documented: nuclear and plastid phylogenetic trees can contain strongly supported discordances, and this incongruence must be reconciled prior to interpretation in evolutionary studies, such as taxonomy, biogeography, and character evolution.

Settling a family feud: a high-level phylogenomic framework for the Gentianales based on 353 nuclear genes and partial plastomes

PREMISE

Comprising five families that vastly differ in species richness-ranging from Gelsemiaceae with 13 species to the Rubiaceae with 13,775 species-members of the Gentianales are often among the most species-rich and abundant plants in tropical forests. Despite considerable phylogenetic work within particular families and genera, several alternative topologies for family-level relationships within Gentianales have been presented in previous studies.

METHODS

Here we present a phylogenomic analysis based on nuclear genes targeted by the Angiosperms353 probe set for approximately 150 species, representing all families and approximately 85% of the formally recognized tribes. We were able to retrieve partial plastomes from off-target reads for most taxa and infer phylogenetic trees for comparison with the nuclear-derived trees.

RESULTS

We recovered high support for over 80% of all nodes. The plastid and nuclear data are largely in agreement, except for some weakly to moderately supported relationships. We discuss the implications of our results for the order's classification, highlighting points of increased support for previously uncertain relationships. Rubiaceae is sister to a clade comprising (Gentianaceae + Gelsemiaceae) + (Apocynaceae + Loganiaceae).

CONCLUSIONS

The higher-level phylogenetic relationships within Gentianales are confidently resolved. In contrast to recent studies, our results support the division of Rubiaceae into two subfamilies: Cinchonoideae and Rubioideae. We do not formally recognize Coptosapelteae and Luculieae within any particular subfamily but treat them as incertae sedis. Our framework paves the way for further work on the phylogenetics, biogeography, morphological evolution, and macroecology of this important group of flowering plants.

Joining forces in Ochnaceae phylogenomics: a tale of two targeted sequencing probe kits

PREMISE

Both universal and family-specific targeted sequencing probe kits are becoming widely used for reconstruction of phylogenetic relationships in angiosperms. Within the pantropical Ochnaceae, we show that with careful data filtering, universal kits are equally as capable in resolving intergeneric relationships as custom probe kits. Furthermore, we show the strength in combining data from both kits to mitigate bias and provide a more robust result to resolve evolutionary relationships.

METHODS

We sampled 23 Ochnaceae genera and used targeted sequencing with two probe kits, the universal Angiosperms353 kit and a family-specific kit. We used maximum likelihood inference with a concatenated matrix of loci and multispecies-coalescence approaches to infer relationships in the family. We explored phylogenetic informativeness and the impact of missing data on resolution and tree support.

RESULTS

For the Angiosperms353 data set, the concatenation approach provided results more congruent with those of the Ochnaceae-specific data set. Filtering missing data was most impactful on the Angiosperms353 data set, with a relaxed threshold being the optimum scenario. The Ochnaceae-specific data set resolved consistent topologies using both inference methods, and no major improvements were obtained after data filtering. Merging of data obtained with the two kits resulted in a well-supported phylogenetic tree.

CONCLUSIONS

The Angiosperms353 data set improved upon data filtering, and missing data played an important role in phylogenetic reconstruction. The Angiosperms353 data set resolved the phylogenetic backbone of Ochnaceae as equally well as the family specific data set. All analyses indicated that both Sauvagesia L. and Campylospermum Tiegh. as currently circumscribed are polyphyletic and require revised delimitation.

Enabling evolutionary studies at multiple scales in Apocynaceae through Hyb-Seq

PREMISE

Apocynaceae is the 10th largest flowering plant family and a focus for study of plant-insect interactions, especially as mediated by secondary metabolites. However, it has few genomic resources relative to its size. Target capture sequencing is a powerful approach for genome reduction that facilitates studies requiring data from the nuclear genome in non-model taxa, such as Apocynaceae.

METHODS

Transcriptomes were used to design probes for targeted sequencing of putatively single-copy nuclear genes across Apocynaceae. The sequences obtained were used to assess the success of the probe design, the intrageneric and intraspecific variation in the targeted genes, and the utility of the genes for inferring phylogeny.

RESULTS

From 853 candidate nuclear genes, 835 were consistently recovered in single copy and were variable enough for phylogenomics. The inferred gene trees were useful for coalescent-based species tree analysis, which showed all subfamilies of Apocynaceae as monophyletic, while also resolving relationships among species within the genus Apocynum. Intraspecific comparison of Elytropus chilensis individuals revealed numerous single-nucleotide polymorphisms with potential for use in population-level studies.

DISCUSSION

Community use of this Hyb-Seq probe set will facilitate and promote progress in the study of Apocynaceae across scales from population genomics to phylogenomics.

The complete plastome sequences of five Aponogeton species (Aponogetonaceae): Insights into the structural organization and mutational hotspots

Members of the aquatic plant genus Aponogeton are widely used commercially in aquariums because of their variable leaf shape and unique inflorescences. However, due to extensive similarity between species in this genus, morphological characters are generally inadequate for taxonomic classification. Currently, molecular makers available for taxonomic and phylogenetic studies of Aponogeton are limited. One approach to clarifying relationships between species in these complex groups is to use divergence hotspot regions within the genome. Here, we sequenced and analyzed the plastomes of five Aponogeton species collected from China, Zambia, and Kenya, and subsequently screened these plastomes for divergent DNA hotspots. The five plastomes are circular structures with sizes ranging from 154,167 bp to 154,860 bp. The Large and the Small Single Copies are separated by two Inverted Repeats. One hundred and thirteen unique genes were identified including 79 protein-coding, 30 tRNA, and four rRNA genes. We found that the most abundant repeats in all but one species were mononucleotide repeats (A/T) and that there were 23 potential RNA ending sites. Interestingly, a ~3 kb inversion, which includes the accD gene, was detected within the Asian species of Aponogeton. The inversion may be related to more frequent exchanges between this region and the nuclear genome. Furthermore, we detected mutational hotspot sites among the five Aponogeton species. Three of these hotspots are intergenic spacer regions (accD-psaI, rbcL-accD and trnH-GUG-psbA) that might be suitable for use as barcodes to resolve intra-generic relationships. We also identified four highly variable protein-coding genes (ccsA, rpl22, rps16 and ycf1) may be used as barcodes to resolve the higher-level phylogenies. Our study will provide valuable molecular resources for the taxonomic and phylogenomic study of the complex genus Aponogeton.

Utility of targeted sequence capture for phylogenomics in rapid, recent angiosperm radiations: Neotropical Burmeistera bellflowers as a case study

Targeted sequence capture is a promising approach for large-scale phylogenomics. However, rapid evolutionary radiations pose significant challenges for phylogenetic inference (e.g. incomplete lineages sorting (ILS), phylogenetic noise), and the ability of targeted nuclear loci to resolve species trees despite such issues remains poorly studied. We test the utility of targeted sequence capture for inferring phylogenetic relationships in rapid, recent angiosperm radiations, focusing on Burmeistera bellflowers (Campanulaceae), which diversified into ~130 species over less than 3 million years. We compared phylogenies estimated from supercontig (exons plus flanking sequences), exon-only, and flanking-only datasets with 506-546 loci (~4.7 million bases) for 46 Burmeistera species/lineages and 10 outgroup taxa. Nuclear loci resolved backbone nodes and many congruent internal relationships with high support in concatenation and coalescent-based species tree analyses, and inferences were largely robust to effects of missing taxa and base composition biases. Nevertheless, species trees were incongruent between datasets, and gene trees exhibited remarkably high levels of conflict (~4-60% congruence, ~40-99% conflict) not simply driven by poor gene tree resolution. Higher gene tree heterogeneity at shorter branches suggests an important role of ILS, as expected for rapid radiations. Phylogenetic informativeness analyses also suggest this incongruence has resulted from low resolving power at short internal branches, consistent with ILS, and homoplasy at deeper nodes, with exons exhibiting much greater risk of incorrect topologies due to homoplasy than other datasets. Our findings suggest that targeted sequence capture is feasible for resolving rapid, recent angiosperm radiations, and that results based on supercontig alignments containing nuclear exons and flanking sequences have higher phylogenetic utility and accuracy than either alone. We use our results to make practical recommendations for future target capture-based studies of Burmeistera and other rapid angiosperm radiations, including that such studies should analyze supercontigs to maximize the phylogenetic information content of loci.

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.

Molecular phylogeny of the diatom genera Amphora and Halamphora (Bacillariophyta) with a focus on morphological and ecological evolution

The taxonomic history of the diatom genus Amphora is one of a broad early morphological concept resulting in the inclusion of a diversity of taxa, followed by an extended period of revision and refinement. The introduction of molecular systematics has increased the pace of revision and has largely resolved the relationships between the major lineages, indicating homoplasy in the evolution of amphoroid symmetry. Within the two largest monophyletic lineages, the genus Halamphora and the now taxonomically refined genus Amphora, the intrageneric morphological and ecological relationships have yet to be explored within a phylogenetic framework. Critical among this is whether the range of morphological features exhibited within these diverse genera are reflective of evolutionary groupings or, as with many previously studied amphoroid features, are nonhomologous when examined phylogenetically. Presented here is a four-marker molecular phylogeny that includes 31 taxa from the genus Amphora and 77 taxa from the genus Halamphora collected from fresh, brackish, and salt waters from coastal and inland habitats of the United States and Japan. These phylogenies illustrate complex patterns in the evolution of frustule morphology and ecology within the genera and the implications of this on the taxonomy, classification, and organization of the genera are discussed.

The diversity and evolution of pollination systems in large plant clades: Apocynaceae as a case study

BACKGROUND AND AIMS

Large clades of angiosperms are often characterized by diverse interactions with pollinators, but how these pollination systems are structured phylogenetically and biogeographically is still uncertain for most families. Apocynaceae is a clade of >5300 species with a worldwide distribution. A database representing >10 % of species in the family was used to explore the diversity of pollinators and evolutionary shifts in pollination systems across major clades and regions.

METHODS

The database was compiled from published and unpublished reports. Plants were categorized into broad pollination systems and then subdivided to include bimodal systems. These were mapped against the five major divisions of the family, and against the smaller clades. Finally, pollination systems were mapped onto a phylogenetic reconstruction that included those species for which sequence data are available, and transition rates between pollination systems were calculated.

KEY RESULTS

Most Apocynaceae are insect pollinated with few records of bird pollination. Almost three-quarters of species are pollinated by a single higher taxon (e.g. flies or moths); 7 % have bimodal pollination systems, whilst the remaining approx. 20 % are insect generalists. The less phenotypically specialized flowers of the Rauvolfioids are pollinated by a more restricted set of pollinators than are more complex flowers within the Apocynoids + Periplocoideae + Secamonoideae + Asclepiadoideae (APSA) clade. Certain combinations of bimodal pollination systems are more common than others. Some pollination systems are missing from particular regions, whilst others are over-represented.

CONCLUSIONS

Within Apocynaceae, interactions with pollinators are highly structured both phylogenetically and biogeographically. Variation in transition rates between pollination systems suggest constraints on their evolution, whereas regional differences point to environmental effects such as filtering of certain pollinators from habitats. This is the most extensive analysis of its type so far attempted and gives important insights into the diversity and evolution of pollination systems in large clades.

Exploring data processing strategies in NGS target enrichment to disentangle radiations in the tribe Cardueae (Compositae)

Target enrichment is a cost-effective sequencing technique that holds promise for elucidating evolutionary relationships in fast-evolving lineages. However, potential biases and impact of bioinformatic sequence treatments in phylogenetic inference have not been thoroughly explored yet. Here, we investigate this issue with an ultimate goal to shed light into a highly diversified group of Compositae (Asteraceae) constituted by four main genera: Arctium, Cousinia, Saussurea, and Jurinea. Specifically, we compared sequence data extraction methods implemented in two easy-to-use workflows, PHYLUCE and HybPiper, and assessed the impact of two filtering practices intended to reduce phylogenetic noise. In addition, we compared two phylogenetic inference methods: (1) the concatenation approach, in which all loci were concatenated in a supermatrix; and (2) the coalescence approach, in which gene trees were produced independently and then used to construct a species tree under coalescence assumptions. Here we confirm the usefulness of the set of 1061 COS targets (a nuclear conserved orthology loci set developed for the Compositae) across a variety of taxonomic levels. Intergeneric relationships were completely resolved: there are two sister groups, Arctium-Cousinia and Saussurea-Jurinea, which are in agreement with a morphological hypothesis. Intrageneric relationships among species of Arctium, Cousinia, and Saussurea are also well defined. Conversely, conflicting species relationships remain for Jurinea. Methodological choices significantly affected phylogenies in terms of topology, branch length, and support. Across all analyses, the phylogeny obtained using HybPiper and the strictest scheme of removing fast-evolving sites was estimated as the optimal. Regarding methodological choices, we conclude that: (1) trees obtained under the coalescence approach are topologically more congruent between them than those inferred using the concatenation approach; (2) refining treatments only improved support values under the concatenation approach; and (3) branch support values are maximized when fast-evolving sites are removed in the concatenation approach, and when a higher number of loci is analyzed in the coalescence approach.

Plastome sequences and exploration of tree-space help to resolve the phylogeny of riceflowers (Thymelaeaceae: Pimelea)

Data sets comprising small numbers of genetic markers are not always able to resolve phylogenetic relationships. This has frequently been the case in molecular systematic studies of plants, with many analyses being based on sequence data from only two or three chloroplast genes. An example of this comes from the riceflowers Pimelea Banks & Sol. ex Gaertn. (Thymelaeaceae), a large genus of flowering plants predominantly distributed in Australia. Despite the considerable morphological variation in the genus, low sequence divergence in chloroplast markers has led to the phylogeny of Pimelea remaining largely uncertain. In this study, we resolve the backbone of the phylogeny of Pimelea in comprehensive Bayesian and maximum-likelihood analyses of plastome sequences from 41 taxa. However, some relationships received only moderate to poor support, and the Pimelea clade contained extremely short internal branches. By using topology-clustering analyses, we demonstrate that conflicting phylogenetic signals can be found across the trees estimated from individual chloroplast protein-coding genes. A relaxed-clock dating analysis reveals that Pimelea arose in the mid-Miocene, with most divergences within the genus occurring during a subsequent rapid diversification. Our new phylogenetic estimate offers better resolution and is more strongly supported than previous estimates, providing a platform for future taxonomic revisions of both Pimelea and the broader subfamily. Our study has demonstrated the substantial improvements in phylogenetic resolution that can be achieved using plastome-scale data sets in plant molecular systematics.

Employing hypothesis testing and data from multiple genomic compartments to resolve recalcitrant backbone nodes in Goodenia s.l. (Goodeniaceae)

Goodeniaceae is a primarily Australian flowering plant family with a complex taxonomy and evolutionary history. Previous phylogenetic analyses have successfully resolved the backbone topology of the largest clade in the family, Goodenia s.l., but have failed to clarify relationships within the species-rich and enigmatic Goodenia clade C, a prerequisite for taxonomic revision of the group. We used genome skimming to retrieve sequences for chloroplast, mitochondrial, and nuclear markers for 24 taxa representing Goodenia s.l., with a particular focus on Goodenia clade C. We performed extensive hypothesis tests to explore incongruence in clade C and evaluate statistical support for clades within this group, using datasets from all three genomic compartments. The mitochondrial dataset is comparable to the chloroplast dataset in providing resolution within Goodenia clade C, though backbone support values within this clade remain low. The hypothesis tests provided an additional, complementary means of evaluating support for clades. We propose that the major subclades of Goodenia clade C (C1-C3 + Verreauxia) are the result of a rapid radiation, and each represents a distinct lineage.

Plastid genomes reveal recurrent formation of allopolyploid Fragaria

PREMISE OF THE STUDY

Recurrent formation of polyploid taxa is a common observation in many plant groups. Haploid, cytoplasmic genomes like the plastid genome can be used to overcome the problem of homeologous genes and recombination in polyploid taxa. Fragaria (Rosaceae) contains several octo- and decaploid species. We use plastome sequences to infer the plastid ancestry of these taxa with special focus on the decaploid Fragaria cascadensis.

METHODS

We used genome skimming of 96 polyploid Fragaria samples on a single Illumina HiSeq 3000 lane to obtain whole plastome sequences. These sequences were used for phylogenetic reconstructions and dating analyses. Ploidy of all samples was inferred with flow cytometry, and plastid inheritance was examined in a controlled cross of F. cascadensis.

KEY RESULTS

The plastid genome phylogeny shows that only the octoploid F. chiloensis is monophyletic, all other polyploid taxa were supported to be para- or polyphyletic. The decaploid Fragaria cascadensis has biparental plastid inheritance and four different plastid donors. Diversification of the F. cascadensis clades occurred in the last 230,000 years. The southern part of its distribution range harbors considerably higher genetic diversity, suggestive of a potential refugium.

CONCLUSIONS

Fragaria cascadensis had at least four independent origins from parents with different plastomes. In contrast, para- and polyphyletic taxa of the octoploid Fragaria species are best explained by incomplete lineage sorting and/or hybridization. Biogeographic patterns in F. cascadensis are probably a result of range shift during the last glacial maximum.

Evolution of pyrrolizidine alkaloid biosynthesis in Apocynaceae: revisiting the defence de-escalation hypothesis

Plants produce specialized metabolites for their defence. However, specialist herbivores adapt to these compounds and use them for their own benefit. Plants attacked predominantly by specialists may be under selection to reduce or eliminate production of co-opted chemicals: the defence de-escalation hypothesis. We studied the evolution of pyrrolizidine alkaloids (PAs) in Apocynaceae, larval host plants for PA-adapted butterflies (Danainae, milkweed and clearwing butterflies), to test if the evolutionary pattern is consistent with de-escalation. We used the first PA biosynthesis specific enzyme (homospermidine synthase, HSS) as tool for reconstructing PA evolution. We found hss orthologues in diverse Apocynaceae species, not all of them known to produce PAs. The phylogenetic analysis showed a monophyletic origin of the putative hss sequences early in the evolution of one Apocynaceae lineage (the APSA clade). We found an hss pseudogene in Asclepias syriaca, a species known to produce cardiac glycosides but no PAs, and four losses of an HSS amino acid motif. APSA clade species are significantly more likely to be Danainae larval host plants than expected if all Apocynaceae species were equally likely to be exploited. Our findings are consistent with PA de-escalation as an adaptive response to specialist attack.

Evolution on the backbone: Apocynaceae phylogenomics and new perspectives on growth forms, flowers, and fruits

PREMISE OF THE STUDY

We provide the largest phylogenetic analyses to date of Apocynaceae in terms of taxa and molecular data as a framework for analyzing the evolution of vegetative and reproductive traits.

METHODS

We produced maximum-likelihood phylogenies of Apocynaceae using 21 plastid loci sampled from 1045 species (nearly 25% of the family) and complete plastomes from 73 species. We reconstructed ancestral states and used model comparisons in a likelihood framework to analyze character evolution across Apocynaceae.

KEY RESULTS

We obtained a well-supported phylogeny of Apocynaceae, resolving poorly understood tribal and subtribal relationships (e.g., among Amsonieae and Hunterieae, within Asclepiadeae), rejecting monophyly of Melodineae and Odontadenieae, and placing previously unsampled and enigmatic taxa (e.g., Pycnobotrya). We provide new insights into the evolution of Apocynaceae, including frequent shifts between herbaceousness and woodiness, reversibility of twining, integrated evolution of the corolla and gynostegium, and ancestral baccate fruits.

CONCLUSIONS

Increased sampling and selection of best-fitting models of evolution provide more resolved and robust estimates of phylogeny and character evolution than obtained in previous studies. Evolutionary inferences are sensitive to choice of phylogenetic frameworks and models.

The unexpected depths of genome-skimming data: A case study examining Goodeniaceae floral symmetry genes

PREMISE OF THE STUDY

The use of genome skimming allows systematists to quickly generate large data sets, particularly of sequences in high abundance (e.g., plastomes); however, researchers may be overlooking data in low abundance that could be used for phylogenetic or evo-devo studies. Here, we present a bioinformatics approach that explores the low-abundance portion of genome-skimming next-generation sequencing libraries in the fan-flowered Goodeniaceae.

METHODS

Twenty-four previously constructed Goodeniaceae genome-skimming Illumina libraries were examined for their utility in mining low-copy nuclear genes involved in floral symmetry, specifically the CYCLOIDEA (CYC)-like genes. De novo assemblies were generated using multiple assemblers, and BLAST searches were performed for CYC1, CYC2, and CYC3 genes.

RESULTS

Overall Trinity, SOAPdenovo-Trans, and SOAPdenovo implementing lower k-mer values uncovered the most data, although no assembler consistently outperformed the others. Using SOAPdenovo-Trans across all 24 data sets, we recovered four CYC-like gene groups (CYC1, CYC2, CYC3A, and CYC3B) from a majority of the species. Alignments of the fragments included the entire coding sequence as well as upstream and downstream regions.

DISCUSSION

Genome-skimming data sets can provide a significant source of low-copy nuclear gene sequence data that may be used for multiple downstream applications.

The East Asian origin of the giant lobelias

PREMISE OF THE STUDY

Rapid radiations are difficult to reconstruct when organismal diversification and biogeographic movement outpace the evolution of genes typically used in phylogenetic analyses. The 125 kb of unique sequence from complete plastid genomes (= plastomes) largely solves the molecular sampling problem, and taxon sampling that triangulates the base of each major subclade largely solves the long-branch attraction problem. This combination of molecular and phylogenetic sampling is used to reconstruct the cosmopolitan radiation of lobeliads, with special focus on the origin of the giant lobelias.

METHODS

An alignment of 18 previously generated and 61 new plastomes was analyzed to produce the phylogenetic estimate upon which the biogeographic reconstruction was based.

KEY RESULTS

Originating in southern Africa, the Lobeliaceae underwent a spectacular cosmopolitan radiation about 20 million years ago. One lineage colonized Madagascar and eastern Asia, which was the source area for the evolution of the giant lobelias. A second lineage colonized the Mediterranean and North America, in quick succession. South America and Australia were also colonized from South Africa, most likely as independent events, but detailed biogeographic reconstruction is limited by inferred extinction events. The south Pacific segregate genera Apetahia and Sclerotheca are inferred to have Hawaiian ancestry. The East African radiation independently reached Ethiopia, West Africa, and Brazil.

CONCLUSIONS

With adequate molecular and taxon sampling, many details of rapid radiations can be accurately inferred. However, not all lineages survived, and analyses of extant species cannot recover details that have been lost due to extinction.

Draft Genome Sequence of Mentha longifolia and Development of Resources for Mint Cultivar Improvement

The genus Mentha encompasses mint species cultivated for their essential oils, which are formulated into a vast array of consumer products. Desirable oil characteristics and resistance to the fungal disease Verticillium wilt are top priorities for the mint industry. However, cultivated mints have complex polyploid genomes and are sterile. Breeding efforts, therefore, require the development of genomic resources for fertile mint species. Here, we present draft de novo genome and plastome assemblies for a wilt-resistant South African accession of Mentha longifolia (L.) Huds., a diploid species ancestral to cultivated peppermint and spearmint. The 353 Mb genome contains 35 597 predicted protein-coding genes, including 292 disease resistance gene homologs, and nine genes determining essential oil characteristics. A genetic linkage map ordered 1397 genome scaffolds on 12 pseudochromosomes. More than two million simple sequence repeats were identified, which will facilitate molecular marker development. The M. longifolia genome is a valuable resource for both metabolic engineering and molecular breeding. This is exemplified by employing the genome sequence to clone and functionally characterize the promoters in a peppermint cultivar, and demonstrating the utility of a glandular trichome-specific promoter to increase expression of a biosynthetic gene, thereby modulating essential oil composition.

The complete plastome sequence of Carissa macrocarpa (Eckl.) A. DC. (Apocynaceae)

In this study, we determined the complete plastome sequence of Carissa macrocarpa (Eckl.) A. DC. (Apocynaceae) (NCBI acc. no. KX364402). The gene order and structure of the C. macrocarpa plastome are similar to those of a typical angiosperm. The complete plastome is 155,297 bp in length, and consists of a large single-copy region of 85,586 bp and a small single-copy region of 18,131 bp, which are separated by two inverted repeats of 25,792 bp. The plastome contains 113 genes, of which 79 are protein-coding genes, 30 are tRNA genes and 4 are rRNA genes. Sixteen genes contained one intron and two genes have two introns. The average A–T content of the plastome is 62.0%. A total of 31 simple sequence repeat loci were identified within the genome. Phylogenetic analysis revealed that C. macrocarpa is a member of the paraphyletic subfamily Rauvolfioideae of Apocynaceae. The sister group relationship of C. macrocarpa to the Apocynoideae–Asclepiadoideae clade is supported by 100% bootstrap values.

The first complete plastid genomes of Melastomataceae are highly structurally conserved

Background

In the past three decades, several studies have predominantly relied on a small sample of the plastome to infer deep phylogenetic relationships in the species-rich Melastomataceae. Here, we report the first full plastid sequences of this family, compare general features of the sampled plastomes to other sequenced Myrtales, and survey the plastomes for highly informative regions for phylogenetics.

Methods

Genome skimming was performed for 16 species spread across the Melastomataceae. Plastomes were assembled, annotated and compared to eight sequenced plastids in the Myrtales. Phylogenetic inference was performed using Maximum Likelihood on six different data sets, where putative biases were taken into account. Summary statistics were generated for all introns and intergenic spacers with suitable size for polymerase chain reaction (PCR) amplification and used to rank the markers by phylogenetic information.

Results

The majority of the plastomes sampled are conserved in gene content and order, as well as in sequence length and GC content within plastid regions and sequence classes. Departures include the putative presence of rps16 and rpl2 pseudogenes in some plastomes. Phylogenetic analyses of the majority of the schemes analyzed resulted in the same topology with high values of bootstrap support. Although there is still uncertainty in some relationships, in the highest supported topologies only two nodes received bootstrap values lower than 95%.

Discussion

Melastomataceae plastomes are no exception for the general patterns observed in the genomic structure of land plant chloroplasts, being highly conserved and structurally similar to most other Myrtales. Despite the fact that the full plastome phylogeny shares most of the clades with the previously widely used and reduced data set, some changes are still observed and bootstrap support is higher. The plastome data set presented here is a step towards phylogenomic analyses in the Melastomataceae and will be a useful resource for future studies.

Chloroplast genomes: diversity, evolution, and applications in genetic engineering

Chloroplasts play a crucial role in sustaining life on earth. The availability of over 800 sequenced chloroplast genomes from a variety of land plants has enhanced our understanding of chloroplast biology, intracellular gene transfer, conservation, diversity, and the genetic basis by which chloroplast transgenes can be engineered to enhance plant agronomic traits or to produce high-value agricultural or biomedical products. In this review, we discuss the impact of chloroplast genome sequences on understanding the origins of economically important cultivated species and changes that have taken place during domestication. We also discuss the potential biotechnological applications of chloroplast genomes.

Development and evolution of extreme synorganization in angiosperm flowers and diversity: a comparison of Apocynaceae and Orchidaceae

BACKGROUND AND AIMS

Apocynaceae and Orchidaceae are two angiosperm families with extreme flower synorganization. They are unrelated, the former in eudicots, the latter in monocots, but they converge in the formation of pollinia and pollinaria, which do not occur in any other angiosperm family, and for which extreme synorganization of floral organs is a precondition. In each family extensive studies on flower development and evolution have been performed; however, newer comparative studies focusing on flower synorganization and involving both families together are lacking.

SCOPE

For this study an extensive search through the morphological literature has been conducted. Based on this and my own studies on flowers in various Apocynaceae and Orchidaceae and complex flowers in other angiosperms with scanning electron microscopy and with microtome section series, a review on convergent floral traits in flower development and architecture in the two families is presented.

KEY FINDINGS

There is a tendency of protracted development of synorganized parts in Apocynaceae and Orchidaceae (development of synorganization of two or more organs begins earlier the more accentuated it is at anthesis). Synorganization (or complexity) also paves the way for novel structures. One of the most conspicuous such novel structures in Apocynaceae is the corona, which is not the product of synorganization of existing organs; however, it is probably enhanced by synorganization of other, existing, floral parts. In contrast to synorganized parts, the corona appears developmentally late.

CONCLUSIONS

Synorganization of floral organs may lead to a large number of convergences in clades that are only very distantly related. The convergences that have been highlighted in this comparative study should be developmentally investigated directly in parallel in future studies.

Phylogenomics of a rapid radiation: is chromosomal evolution linked to increased diversification in north american spiny lizards (Genus Sceloporus)?

Background

Resolving the short phylogenetic branches that result from rapid evolutionary diversification often requires large numbers of loci. We collected targeted sequence capture data from 585 nuclear loci (541 ultraconserved elements and 44 protein-coding genes) to estimate the phylogenetic relationships among iguanian lizards in the North American genus Sceloporus. We tested for diversification rate shifts to determine if rapid radiation in the genus is correlated with chromosomal evolution.

Results

The phylogenomic trees that we obtained for Sceloporus using concatenation and coalescent-based species tree inference provide strong support for the monophyly and interrelationships among nearly all major groups. The diversification analysis supported one rate shift on the Sceloporus phylogeny approximately 20–25 million years ago that is associated with the doubling of the speciation rate from 0.06 species/million years (Ma) to 0.15 species/Ma. The posterior probability for this rate shift occurring on the branch leading to the Sceloporus species groups exhibiting increased chromosomal diversity is high (posterior probability = 0.997).

Conclusions

Despite high levels of gene tree discordance, we were able to estimate a phylogenomic tree for Sceloporus that solves some of the taxonomic problems caused by previous analyses of fewer loci. The taxonomic changes that we propose using this new phylogenomic tree help clarify the number and composition of the major species groups in the genus. Our study provides new evidence for a putative link between chromosomal evolution and the rapid divergence and radiation of Sceloporus across North America.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0628-x) contains supplementary material, which is available to authorized users.

Utilizing next-generation sequencing to resolve the backbone of the Core Goodeniaceae and inform future taxonomic and floral form studies

Though considerable progress has been made in inferring phylogenetic relationships of many plant lineages, deep unresolved nodes remain a common problem that can impact downstream efforts, including taxonomic decision-making and character reconstruction. The Core Goodeniaceae is a group affected by this issue: data from the plastid regions trnL-trnF and matK have been insufficient to generate adequate support at key nodes along the backbone of the phylogeny. We performed genome skimming for 24 taxa representing major clades within Core Goodeniaceae. The plastome coding regions (CDS) and nuclear ribosomal repeats (NRR) were assembled and complemented with additional accessions sequenced for nuclear G3PDH and plastid trnL-trnF and matk. The CDS, NRR, and G3PDH alignments were analyzed independently and topology tests were used to detect the alignments' ability to reject alternative topologies. The CDS, NRR, and G3PDH alignments independently supported a Brunonia (Scaevola s.l. (Coopernookia (Goodenia s.l.))) backbone topology, but within Goodenia s.l., the strongly-supported plastome topology (Goodenia A (Goodenia B (Velleia+Goodenia C))) contrasts with the poorly supported nuclear topology ((Goodenia A+Goodenia B) (Velleia+Goodenia C)). A fully resolved and maximally supported topology for Core Goodeniaceae was recovered from the plastome CDS, and there is excellent support for most of the major clades and relationships among them in all alignments. The composition of these seven major clades renders many of the current taxonomic divisions non-monophyletic, prompting us to suggest that Goodenia may be split into several segregate genera.

Utility of characters evolving at diverse rates of evolution to resolve quartet trees with unequal branch lengths: analytical predictions of long-branch effects

BACKGROUND

The detection and avoidance of "long-branch effects" in phylogenetic inference represents a longstanding challenge for molecular phylogenetic investigations. A consequence of parallelism and convergence, long-branch effects arise in phylogenetic inference when there is unequal molecular divergence among lineages, and they can positively mislead inference based on parsimony especially, but also inference based on maximum likelihood and Bayesian approaches. Long-branch effects have been exhaustively examined by simulation studies that have compared the performance of different inference methods in specific model trees and branch length spaces.

RESULTS

In this paper, by generalizing the phylogenetic signal and noise analysis to quartets with uneven subtending branches, we quantify the utility of molecular characters for resolution of quartet phylogenies via parsimony. Our quantification incorporates contributions toward the correct tree from either signal or homoplasy (i.e. "the right result for either the right reason or the wrong reason"). We also characterize a highly conservative lower bound of utility that incorporates contributions to the correct tree only when they correspond to true, unobscured parsimony-informative sites (i.e. "the right result for the right reason"). We apply the generalized signal and noise analysis to classic quartet phylogenies in which long-branch effects can arise due to unequal rates of evolution or an asymmetrical topology. Application of the analysis leads to identification of branch length conditions in which inference will be inconsistent and reveals insights regarding how to improve sampling of molecular loci and taxa in order to correctly resolve phylogenies in which long-branch effects are hypothesized to exist.

CONCLUSIONS

The generalized signal and noise analysis provides analytical prediction of utility of characters evolving at diverse rates of evolution to resolve quartet phylogenies with unequal branch lengths. The analysis can be applied to identifying characters evolving at appropriate rates to resolve phylogenies in which long-branch effects are hypothesized to occur.

Employing Phylogenomics to Resolve the Relationships among Cnidarians, Ctenophores, Sponges, Placozoans, and Bilaterians

Despite an explosion in the amount of sequence data, phylogenomics has failed to settle controversy regarding some critical nodes on the animal tree of life. Understanding relationships among Bilateria, Ctenophora, Cnidaria, Placozoa, and Porifera is essential for studying how complex traits such as neurons, muscles, and gastrulation have evolved. Recent studies have cast doubt on the historical viewpoint that sponges are sister to all other animal lineages with recent studies recovering ctenophores as sister. However, the ctenophore-sister hypothesis has been criticized as unrealistic and caused by systematic error. We review past phylogenomic studies and potential causes of systematic error in an effort to identify areas that can be improved in future studies. Increased sampling of taxa, less missing data, and a priori removal of sequences and taxa that may cause systematic error in phylogenomic inference will likely be the most fruitful areas of focus when assembling future datasets. Ultimately, we foresee metazoan relationships being resolved with higher support in the near future, and we caution against dismissing novel hypotheses merely because they conflict with historical viewpoints of animal evolution.

Evaluation of a Phylogenetic Marker Based on Genomic Segment B of Infectious Bursal Disease Virus: Facilitating a Feasible Incorporation of this Segment to the Molecular Epidemiology Studies for this Viral Agent

Background

Infectious bursal disease (IBD) is a highly contagious and acute viral disease, which has caused high mortality rates in birds and considerable economic losses in different parts of the world for more than two decades and it still represents a considerable threat to poultry. The current study was designed to rigorously measure the reliability of a phylogenetic marker included into segment B. This marker can facilitate molecular epidemiology studies, incorporating this segment of the viral genome, to better explain the links between emergence, spreading and maintenance of the very virulent IBD virus (vvIBDV) strains worldwide.

Methodology/Principal Findings

Sequences of the segment B gene from IBDV strains isolated from diverse geographic locations were obtained from the GenBank Database; Cuban sequences were obtained in the current work. A phylogenetic marker named B-marker was assessed by different phylogenetic principles such as saturation of substitution, phylogenetic noise and high consistency. This last parameter is based on the ability of B-marker to reconstruct the same topology as the complete segment B of the viral genome. From the results obtained from B-marker, demographic history for both main lineages of IBDV regarding segment B was performed by Bayesian skyline plot analysis. Phylogenetic analysis for both segments of IBDV genome was also performed, revealing the presence of a natural reassortant strain with segment A from vvIBDV strains and segment B from non-vvIBDV strains within Cuban IBDV population.

Conclusions/Significance

This study contributes to a better understanding of the emergence of vvIBDV strains, describing molecular epidemiology of IBDV using the state-of-the-art methodology concerning phylogenetic reconstruction. This study also revealed the presence of a novel natural reassorted strain as possible manifest of change in the genetic structure and stability of the vvIBDV strains. Therefore, it highlights the need to obtain information about both genome segments of IBDV for molecular epidemiology studies.

Error, signal, and the placement of Ctenophora sister to all other animals

Elucidating relationships among early animal lineages has been difficult, and recent phylogenomic analyses place Ctenophora sister to all other extant animals, contrary to the traditional view of Porifera as the earliest-branching animal lineage. To date, phylogenetic support for either ctenophores or sponges as sister to other animals has been limited and inconsistent among studies. Lack of agreement among phylogenomic analyses using different data and methods obscures how complex traits, such as epithelia, neurons, and muscles evolved. A consensus view of animal evolution will not be accepted until datasets and methods converge on a single hypothesis of early metazoan relationships and putative sources of systematic error (e.g., long-branch attraction, compositional bias, poor model choice) are assessed. Here, we investigate possible causes of systematic error by expanding taxon sampling with eight novel transcriptomes, strictly enforcing orthology inference criteria, and progressively examining potential causes of systematic error while using both maximum-likelihood with robust data partitioning and Bayesian inference with a site-heterogeneous model. We identified ribosomal protein genes as possessing a conflicting signal compared with other genes, which caused some past studies to infer ctenophores and cnidarians as sister. Importantly, biases resulting from elevated compositional heterogeneity or elevated substitution rates are ruled out. Placement of ctenophores as sister to all other animals, and sponge monophyly, are strongly supported under multiple analyses, herein.

Tropical Refuges with Exceptionally High Phylogenetic Diversity Reveal Contrasting Phylogenetic Structures

Loss of phylogenetic diversity (PD) has gained increasing attention in conservation biology. However, PD is not equally distributed in a phylogeny and can be better assessed when species relatedness (phylogenetic structure: PS) is also considered. Here, we investigate PD and PS in two refuges of biodiversity in northeastern Brazil: the Bahia Costal Forest (BCF) in the Atlantic Forest domain and Chapada Diamantina (CD) in the Caatinga domain. We used geographic data of 205 species at two spatial scales and a chronogram of Apocynaceae based on matK sequences to estimate PD and PS. Our results show an exceptionally high PD in both refuges, overdispersed in BCF and clustered in CD, although this difference is less evident or absent for recent relationships, especially at a smaller spatial scale. Overall, PS suggests long-term competitive exclusion under climatic stability, currently balanced by habitat filtering, in BCF, and biome conservatism and limited dispersal leading to in situ diversification and high density of microendemics in CD. The phylogenetically clustered flora in CD, also threatened by climate changes, are naturally more vulnerable than BCF. Therefore, while in situ conservation may ensure protection of biodiversity in BCF, emergency ex situ conservation is strongly recommended in CD.