Phylogenomics insights into gene evolution, rapid species diversification, and morphological innovation of the apple tribe (Maleae, Rosaceae).

Maleae is one of the most widespread tribes of Rosaceae and includes several important fruit crops and ornamental plants. We used nuclear genes from 62 transcriptomes/genomes, including 26 newly generated transcriptomes, to reconstruct a well-supported phylogeny and study the evolution of fruit and leaf morphology and the possible effect of whole genome duplication (WGD). Our phylogeny recovered 11 well-supported clades and supported the monophyly of most genera (except Malus, Sorbus, and Pourthiaea) with at least two sampled species. A WGD was located to the most recent common ancestor (MRCA) of Maleae and dated to c. 54 million years ago (Ma) near the Early Eocene Climatic Optimum, supporting Gillenieae (x = 9) being a parental lineage of Maleae (x = 17) and including duplicate regulatory genes related to the origin of the fleshy pome fruit. Whole genome duplication-derived paralogs that are retained in specific lineages but lost in others are predicted to function in development, metabolism, and other processes. An upshift of diversification and innovations of fruit and leaf morphologies occurred at the MRCA of the Malinae subtribe, coinciding with the Eocene-Oligocene transition (c. 34 Ma), following a lag from the time of the WGD event. Our results provide new insights into the Maleae phylogeny, its rapid diversification, and morphological and molecular evolution.

Genome sequence and salinity adaptation of the desert shrub Nitraria sibirica (Nitrariaceae, Sapindales).

The genetic bases of halophytes for salinity tolerance are crucial for genetically breeding salt-tolerant crops. All natural Nitrariaceae species that exclusively occur in arid environments are highly tolerant to salt stress, but the underlying genomic bases to this adaptation remain unknown. Here we present a high-quality, chromosome-level genome sequence of Nitraria sibirica, with an assembled size of 456.66 Mb and 23,365 annotated genes. Phylogenomic analyses confirmed N. sibirica as the sister to all other sampled representatives from other families in Sapindales, and no lineage-specific whole-genome duplication was found except the gamma triplication event. Still, we found that the genes involved in K + retention, energy supply, and Fe absorption expanded greatly in N. sibirica. Deep transcriptome analyses showed that leaf photosynthesis and cuticular wax formation in roots were enhanced under salt treatments. Furthermore, many transcription factors involved in salt tolerance changed their expressions significantly and displayed tissue- and concentration-dependent signaling in response to salt stress. Additionally, we found vacuolar Na + compartmentalization is an ongoing process under salt treatment, while Na + exclusion tends to function at high salt concentrations. These genomic and transcriptomic changes conferred salt tolerance in N. sibirica and pave the way for future breeding of salt-tolerant crops.

Target sequence capture data shed light on the deeper evolutionary relationships of subgenus Chamaecerasus in Lonicera (Caprifoliaceae).

The genus Lonicera L. is widely distributed in the north temperate zone and is well-known for its high species richness and morphological diversity. Previous studies have suggested that many sections of Lonicera are not monophyletic and phylogenetic relationships within the genus are still poorly resolved. In this study, we sampled 37 accessions of Lonicera, covering four sections of subgenus Chamaecerasus plus six outgroup taxa, to recover the main clades of Lonicera based on sequences of nuclear loci generated by target enrichment and cpDNA from genome skimming. We found extensive cytonuclear discordance across the subgenus. Both nuclear and plastid phylogenetic analyses supported subgenus Chamaecerasus sister to subgenus Lonicera. Within subgenus Chamaecerasus, sections Isika and Nintooa were each polyphyletic. Based on the nuclear and chloroplast phylogenies, we propose to merge Lonicera korolkowii into section Coeloxylosteum and Lonicera caerulea into section Nintooa. In addition, Lonicera is estimated to have originated in the mid Oligocene (26.45 Ma). The stem age of section Nintooa was estimated to be 17.09 Ma (95% HPD: 13.30-24.45). The stem age of subgenus Lonicera was estimated to be 16.35 Ma (95% HPD: 14.12-23.66). Ancestral area reconstruction analyses indicate that subgenus Chamaecerasus originated in East Asia and Central Asia. In addition, sections Coeloxylosteum and Nintooa originated in East Asia, with subsequent dispersals into other areas. The aridification of the Asian interior likely promoted the rapid radiation of sections Coeloxylosteum and Nintooa within this region. Moreover, our biogeographic analysis fully supports the Bering and the North Atlantic Land Bridge hypotheses for the intercontinental migrations in the Northern Hemisphere. Overall, this study provides new insights into the taxonomically complex lineages of subgenus Chamaecerasus and the process of speciation.

Phylogenomic Analyses of Alismatales Shed Light into Adaptations to Aquatic Environments.

Land plants first evolved from freshwater algae, and flowering plants returned to water as early as the Cretaceous and multiple times subsequently. Alismatales is the largest clade of aquatic angiosperms including all marine angiosperms, as well as terrestrial plants. We used Alismatales to explore plant adaptations to aquatic environments by analyzing a data set that included 95 samples (89 Alismatales species) covering four genomes and 91 transcriptomes (59 generated in this study). To provide a basis for investigating adaptations, we assessed phylogenetic conflict and whole-genome duplication (WGD) events in Alismatales. We recovered a relationship for the three main clades in Alismatales as (Tofieldiaceae, Araceae) + core Alismatids. We also found phylogenetic conflict among the three main clades that was best explained by incomplete lineage sorting and introgression. Overall, we identified 18 putative WGD events across Alismatales. One of them occurred at the most recent common ancestor of core Alismatids, and three occurred at seagrass lineages. We also found that lineage and life-form were both important for different evolutionary patterns for the genes related to freshwater and marine adaptation. For example, several light- or ethylene-related genes were lost in the seagrass Zosteraceae, but are present in other seagrasses and freshwater species. Stomata-related genes were lost in both submersed freshwater species and seagrasses. Nicotianamine synthase genes, which are important in iron intake, expanded in both submersed freshwater species and seagrasses. Our results advance the understanding of the adaptation to aquatic environments and WGDs using phylogenomics.

Phylogenomic analyses in Phrymaceae reveal extensive gene tree discordance in relationships among major clades.

PREMISE: Phylogenomic datasets using genomes and transcriptomes provide rich opportunities beyond resolving bifurcating phylogenetic relationships. Monkeyflower (Phrymaceae) is a model system for evolutionary ecology. However, it lacks a well-supported phylogeny as a basis for a stable taxonomy and for macroevolutionary comparisons. METHODS: We sampled 24 genomes and transcriptomes in Phrymaceae and closely related families, including eight newly sequenced transcriptomes. We reconstructed the phylogeny using IQ-TREE and ASTRAL, evaluated gene tree discordance using PhyParts, Quartet Sampling, and a cloudogram, and carried out reticulation analyses using PhyloNet and HyDe. We searched for whole genome duplication (WGD) events using chromosome numbers, synonymous distances, and gene duplication events as evidence. RESULTS: Most gene trees support the monophyly of Phrymaceae and each of its tribes. Most gene trees also support tribe Mimuleae being sister to Phrymeae + Diplaceae + Leucocarpeae, with extensive gene tree discordance among the latter three. Despite the discordance, the monophyly of Mimulus s.l. is rejected, and no individual reticulation event among the Phrymaceae tribes is well-supported. Reticulation likely occurred among Erythranthe bicolor and closely related species. No ancient WGD was detected in Phrymaceae. Instead, small-scale duplications are among potential drivers of macroevolutionary diversification of Phrymaceae. CONCLUSIONS: We show that analysis of reticulate evolution is sensitive to taxon sampling and methods used. We also demonstrate that phylogenomic datasets using genomes and transcriptomes present rich opportunities to investigate gene family evolution and genome duplication events involved in lineage diversification and adaptation.

On the hybrid origin of the C2 Salsola divaricata agg. (Amaranthaceae) from C3 and C4 parental lineages.

C2 photosynthesis is characterised using recapturing photorespiratory CO2 by RuBisCo in Kranz-like cells and is therefore physiologically intermediate between C3 and C4 photosynthesis. C2 can be interpreted as an evolutionary precursor of C4 and/or as the result of hybridisation between a C3 and C4 lineage. We compared the expression of photosynthetic traits among populations of the Salsola divaricata agg. (C2 ) from humid subtropical to arid habitats on the coasts of the Canary Islands and Morocco and subjected them to salt and drought treatments. We screened for enhanced C4 -like expression of traits related to habitat or treatment. We estimated species trees with a transcriptome dataset of Salsoleae and explored patterns of gene tree discordance. With phylogenetic networks and hybridisation analyses we tested for the hybrid origin of the Salsola divaricata agg. We observed distinct independent variation of photosynthetic traits within and among populations and no clear evidence for selection towards C4 -like trait expression in more stressful habitats or treatments. We found reticulation and gene tree incongruence in Salsoleae supporting a putative hybrid origin of the Salsola divaricata agg. C2 photosynthesis in the Salsola divaricata agg. combines traits inherited from its C3 and C4 parental lineages and seems evolutionarily stable, possibly well adapted to a wide climatic amplitude.

Phylogenomic analyses of the East Asian endemic Abelia (Caprifoliaceae) shed insights into the temporal and spatial diversification history with widespread hybridization.

BACKGROUND AND AIMS: Abelia (Caprifoliaceae) is a small genus with five species, including one artificial hybrid and several natural hybrids. The genus has a discontinuous distribution in Mainland China, Taiwan Island and the Ryukyu Islands, providing a model system to explore the mechanisms of species dispersal in the East Asian flora. However, the current phylogenetic relationships within Abelia remain uncertain. METHODS: We reconstructed the phylogenetic relationships within Abelia using nuclear loci generated by target enrichment and plastomes from genome skimming. Divergence time estimation, ancestral area reconstruction and ecological niche modelling (ENM) were used to examine the diversification history of Abelia. KEY RESULTS: We found extensive cytonuclear discordance across the genus. By integrating lines of evidence from molecular phylogenies, divergence times and morphology, we propose to merge Abelia macrotera var. zabelioides into A. uniflora. Network analyses suggested that there have been multiple widespread hybridization events among Abelia species. These hybridization events may have contributed to the speciation mechanism and resulted in the high observed morphological diversity. The diversification of Abelia began in the early Eocene, followed by A. chinensis var. ionandra colonizing Taiwan Island during the Middle Miocene. The ENM results suggested an expansion of climatically suitable areas during the Last Glacial Maximum and range contraction during the Last Interglacial. Disjunction between the Himalayan-Hengduan Mountain region and Taiwan Island is probably the consequence of topographical isolation and postglacial contraction. CONCLUSIONS: We used genomic data to reconstruct the phylogeny of Abelia and found a clear pattern of reticulate evolution in the group. In addition, our results suggest that shrinkage of postglacial range and the heterogeneity of the terrain have led to the disjunction between Mainland China and Taiwan Island. This study provides important new insights into the speciation process and taxonomy of Abelia.

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

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

Disentangling Sources of Gene Tree Discordance in Phylogenomic Data Sets: Testing Ancient Hybridizations in Amaranthaceae s.l.

Gene tree discordance in large genomic data sets can be caused by evolutionary processes such as incomplete lineage sorting and hybridization, as well as model violation, and errors in data processing, orthology inference, and gene tree estimation. Species tree methods that identify and accommodate all sources of conflict are not available, but a combination of multiple approaches can help tease apart alternative sources of conflict. Here, using a phylotranscriptomic analysis in combination with reference genomes, we test a hypothesis of ancient hybridization events within the plant family Amaranthaceae s.l. that was previously supported by morphological, ecological, and Sanger-based molecular data. The data set included seven genomes and 88 transcriptomes, 17 generated for this study. We examined gene-tree discordance using coalescent-based species trees and network inference, gene tree discordance analyses, site pattern tests of introgression, topology tests, synteny analyses, and simulations. We found that a combination of processes might have generated the high levels of gene tree discordance in the backbone of Amaranthaceae s.l. Furthermore, we found evidence that three consecutive short internal branches produce anomalous trees contributing to the discordance. Overall, our results suggest that Amaranthaceae s.l. might be a product of an ancient and rapid lineage diversification, and remains, and probably will remain, unresolved. This work highlights the potential problems of identifiability associated with the sources of gene tree discordance including, in particular, phylogenetic network methods. Our results also demonstrate the importance of thoroughly testing for multiple sources of conflict in phylogenomic analyses, especially in the context of ancient, rapid radiations. We provide several recommendations for exploring conflicting signals in such situations. [Amaranthaceae; gene tree discordance; hybridization; incomplete lineage sorting; phylogenomics; species network; species tree; transcriptomics.].

Chloroplast primers for clade-wide phylogenetic studies of Thalictrum.

PREMISE: Chloroplast primers were developed for phylogenetic and comparative studies in Thalictrum (Ranunculaceae). METHODS AND RESULTS: We assembled and annotated the complete plastome sequence of T. thalictroides by combining multiple whole genome sequencing libraries. Using transcriptome-sequencing libraries, we also assembled a partial plastome of the related species T. hernandezii. From the newly assembled plastomes and one previously sequenced plastome, we designed and validated 28 primer pairs to target variable portions of the chloroplast genome in Thalictrum. Furthermore, we tested the validated primers in 62 species of Thalictrum. The total alignment length of the 28 regions was 15,268 bp with 2443 variable sites and 92% character occupancy. CONCLUSIONS: The newly developed chloroplast primer pairs improve the phylogenetic resolution (bootstrap support and tree certainty) in Thalictum and will be a useful resource for future phylogenetic and evolutionary studies for species in the genus and in close relatives in Thalictroideae.

Performance of gene expression analyses using de novo assembled transcripts in polyploid species.

MOTIVATION: Quality of gene expression analyses using de novo assembled transcripts in species that experienced recent polyploidization remains unexplored. RESULTS: Differential gene expression (DGE) analyses using putative genes inferred by Trinity, Corset and Grouper performed slightly differently across five plant species that experienced various polyploidy histories. In species that lack recent polyploidy events that occurred in the past several millions of years, DGE analyses using de novo assembled transcriptomes identified 54-82% of the differentially expressed genes recovered by mapping reads to the reference genes. However, in species that experienced more recent polyploidy events, the percentage decreased to 21-65%. Gene co-expression network analyses using de novo assemblies versus mapping to the reference genes recovered the same module that significantly correlated with treatment in one species that lacks recent polyploidization. AVAILABILITY AND IMPLEMENTATION: Commands and scripts used in this study are available at https://bitbucket.org/lychen83/chen_et_al_2018_benchmark_dge/; Analysis files are available at Dryad doi: 10.5061/dryad.4p6n481. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Three new species of Lachemilla (Rosaceae) from South America.

Three new species of Lachemilla (Rosaceae), two from Colombia and one from Peru, are described and illustrated. Lachemillarothmaleriana is characterized by its stout stems, sericeous-villous indumentum, and wide ascending sheaths with trilobate lateral lobes. Lachemillaargentea presents a unique combination of tripartite basal leaves with an adaxial silvery villous indumentum, and decumbent branches with verticillate lobed sheaths. Finally, Lachemillacyanea has distinctly basal reniform leaves with a blue-green color and hirsute pubescence. Phylogenetic analyses of the nuclear ribosomal cistron and multiple regions of the plastid genome revealed the allopolyploid origin of the three new taxa.

ResumenTres nuevas especies de Lachemilla (Rosaceae), dos de Colombia y una de Perú, se describen e ilustran. Lachemillarothmaleriana se caracteriza por tener tallos robustos, indumento sericeo-villoso, y verticilos anchos ascendentes con lóbulos laterales trilobados. Lachemillaargentea presenta una combinación única de hojas basales tripartitas con indumento villoso plateado en el haz y ramas decumbentes con verticilos lobulados. Finalmente, Lachemillacyanea presenta hojas basales distintivas de forma reniforme, color azul­verdoso y pubescencia hirsuta. Análisis filogenéticos del cistrón ribosomal nuclear y múltiples regiones del genoma cloroplástico revelaron el origen alopoliploide de los tres nuevos taxones.

Extensive allopolyploidy in the neotropical genus Lachemilla (Rosaceae) revealed by PCR-based target enrichment of the nuclear ribosomal DNA cistron and plastid phylogenomics.

PREMISE OF THE STUDY: Polyploidy has been long recognized as an important force in plant evolution. Previous studies had suggested widespread occurrence of polyploidy and the allopolyploid origin of several species in the diverse neotropical genus Lachemilla (Rosaceae). Nonetheless, this evidence has relied mostly on patterns of cytonuclear discordance, and direct evidence from nuclear allelic markers is still needed. METHODS: Here we used PCR target enrichment in combination with high throughput sequencing to obtain multiple copies of the nuclear ribosomal (nr) DNA cistron and 45 regions of the plastid genome (cpDNA) from 219 accessions representing 48 species of Lachemilla and to explore the allopolyploid origin of species in this group. KEY RESULTS: We were able to identify multiple nrDNA ribotypes and establish clear evidence of allopolyploidy in 33 species of Lachemilla, showing that this condition is common and widespread in the genus. Additionally, we found evidence for three autopolyploid species. We also established multiple, independent origins of several allopolyploid species. Finally, based solely on the cpDNA phylogeny, we identified that the monotypic genus Farinopsis is the sister group of Lachemilla and allied genera within subtribe Fragariinae. CONCLUSIONS: Our study demonstrates the utility of the nuclear ribosomal DNA cistron to detect allopolyploidy when concerted evolution of this region is not complete. Additionally, with a robust chloroplast phylogeny in place, the direction of hybridization events can be established, and multiple, independent origins of allopolyploid species can be identified.

Fluidigm2PURC: Automated processing and haplotype inference for double-barcoded PCR amplicons.

PREMISE OF THE STUDY: Targeted enrichment strategies for phylogenomic inference are a time- and cost-efficient way to collect DNA sequence data for large numbers of individuals at multiple, independent loci. Automated and reproducible processing of these data is a crucial step for researchers conducting phylogenetic studies. METHODS AND RESULTS: We present Fluidigm2PURC, an open source Python utility for processing paired-end Illumina data from double-barcoded PCR amplicons. In combination with the program PURC (Pipeline for Untangling Reticulate Complexes), our scripts process raw FASTQ files for analysis with PURC and use its output to infer haplotypes for diploids, polyploids, and samples with unknown ploidy. We demonstrate the use of the pipeline with an example data set from the genus Thalictrum (Ranunculaceae). CONCLUSIONS: Fluidigm2PURC is freely available for Unix-like operating systems on GitHub (https://github.com/pblischak/fluidigm2purc) and for all operating systems through Docker (https://hub.docker.com/r/pblischak/fluidigm2purc).

Phylogenomic analyses reveal a deep history of hybridization and polyploidy in the Neotropical genus Lachemilla (Rosaceae).

Hybridization, incomplete lineage sorting, and phylogenetic error produce similar incongruence patterns, representing a great challenge for phylogenetic reconstruction. Here, we use sequence capture data and multiple species tree and species network approaches to resolve the backbone phylogeny of the Neotropical genus Lachemilla, while distinguishing among sources of incongruence. We used 396 nuclear loci and nearly complete plastome sequences from 27 species to clarify the relationships among the major groups of Lachemilla, and explored multiple sources of conflict between gene trees and species trees inferred with a plurality of approaches. All phylogenetic methods recovered the four major groups previously proposed for Lachemilla, but species tree methods recovered different topologies for relationships between these four clades. Species network analyses revealed that one major clade, Orbiculate, is likely of ancient hybrid origin, representing one of the main sources of incongruence among the species trees. Additionally, we found evidence for a potential whole genome duplication event shared by Lachemilla and allied genera. Lachemilla shows clear evidence of ancient and recent hybridization throughout the evolutionary history of the group. Also, we show the necessity to use phylogenetic network approaches that can simultaneously accommodate incomplete lineage sorting and gene flow when studying groups that show patterns of reticulation.

Primers for Castilleja and their utility across Orobanchaceae: I. Chloroplast primers.

PREMISE OF THE STUDY: Chloroplast primers were developed from genomic data for the taxonomically challenging genus Castilleja. We further tested the broader utility of these primers across Orobanchaceae, identifying a core set of chloroplast primers amplifying across the clade. METHODS AND RESULTS: Using a combination of three low-coverage Castilleja genomes and sequence data from 12 Castilleja plastomes, 76 primer combinations were specifically designed and tested for Castilleja. The primers targeted the most variable portions of the plastome and were validated for their applicability across the clade. Of these, 38 primer combinations were subsequently evaluated in silico and then validated across other major clades in Orobanchaceae. CONCLUSIONS: These results demonstrate the utility of these primers, not only across Castilleja, but for other clades in Orobanchaceae-particularly hemiparasitic lineages-and will contribute to future phylogenetic studies of this important clade of parasitic plants.

Lachemilla mexiquense (Rosaceae), a new species from Mexico.

A new species of Lachemilla (Rosaceae), Lachemilla mexiquense D.F. Morales-B., from Mexico is described and illustrated. This species is similar to Lachemilla aphanoides by its tripartite leaves and glomerulate inflorescence with entirely glabrous flowers, but it differs by its stonoliferous habit, persistent basal leaves and basal stipules, and smaller flowers with a campanulate-elongate hypanthium and single carpel. A key to the species of Lachemilla in Mexico is provided.

ResumenUna nueva especie de Lachemilla (Rosaceae), Lachemilla mexiquense D.F. Morales-B., de México se describe e ilustra. Esta especie es similar a Lachemilla aphanoides por sus hojas tripartitas e inflorescencias glomeruladas con flores completamente glabras, pero difiere por su hábito estolonífero, hojas basales y estipulas basales persistentes y flores de menor tamaño con hipantio campanulado-elongado con un solo carpelo. Se provee una clave para las especies de Lachemilla en México.