Quality and quantity of data recovered from massively parallel sequencing: Examples in Asparagales and Poaceae

Young evolutionary origins of dioecy in the genus Asparagus

PREMISE

Dioecy (separate sexes) has independently evolved numerous times across the angiosperm phylogeny and is recently derived in many lineages. However, our understanding is limited regarding the evolutionary mechanisms that drive the origins of dioecy in plants. The recent and repeated evolution of dioecy across angiosperms offers an opportunity to make strong inferences about the ecological, developmental, and molecular factors influencing the evolution of dioecy, and thus sex chromosomes. The genus Asparagus (Asparagaceae) is an emerging model taxon for studying dioecy and sex chromosome evolution, yet estimates for the age and origin of dioecy in the genus are lacking.

METHODS

We use plastome sequences and fossil time calibrations in phylogenetic analyses to investigate the age and origin of dioecy in the genus Asparagus. We also review the diversity of sexual systems present across the genus to address contradicting reports in the literature.

RESULTS

We estimate that dioecy evolved once or twice approximately 2.78-3.78 million years ago in Asparagus, of which roughly 27% of the species are dioecious and the remaining are hermaphroditic with monoclinous flowers.

CONCLUSIONS

Our findings support previous work implicating a young age and the possibility of two origins of dioecy in Asparagus, which appear to be associated with rapid radiations and range expansion out of Africa. Lastly, we speculate that paleoclimatic oscillations throughout northern Africa may have helped set the stage for the origin(s) of dioecy in Asparagus approximately 2.78-3.78 million years ago.

REFMAKER: make your own reference to target nuclear loci in low coverage genome skimming libraries. Phylogenomic application in Sapotaceae

Genome skimming approach is widely used in plant systematics to infer phylogenies mostly from organelle genomes. However, organelles represent only 10% of the produced libraries, and the low coverage associated with these libraries (< 3X) prevents the capture of nuclear sequences, which are not always available in non-model organisms or limited to the ribosomal regions. We developed REFMAKER, a user-friendly pipeline, to create specific sets of nuclear loci that can next be extracted directly from the genome skimming libraries. For this, a catalogue is built from the meta-assembly of each library contigs and cleaned by selecting the nuclear regions and removing duplicates from clustering steps. Libraries are next mapped onto this catalogue and consensus sequences are generated to produce a ready-to-use phylogenetic matrix following different filtering parameters aiming at removing putative errors and paralogous sequences. REFMAKER allowed us to infer a well resolved phylogeny in Capurodendron (Sapotaceae) on 67 nuclear loci from low-coverage libraries (<1X). The resulting phylogeny is concomitant with one previously inferred on 638 nuclear genes from target enrichment libraries. While it remains preliminary because of this low sequencing depth, REFMAKER therefore opens perspectives in phylogenomics by allowing nuclear phylogeny reconstructions with genome skimming datasets.

Phylogenomics and morphological evolution of the mega-diverse genus Artemisia (Asteraceae: Anthemideae): implications for its circumscription and infrageneric taxonomy

BACKGROUND AND AIMS

Artemisia is a mega-diverse genus consisting of ~400 species. Despite its medicinal importance and ecological significance, a well-resolved phylogeny for global Artemisia, a natural generic delimitation and infrageneric taxonomy remain missing, owing to the obstructions from limited taxon sampling and insufficient information on DNA markers. Its morphological characters, such as capitulum, life form and leaf, show marked variations and are widely used in its infrageneric taxonomy. However, their evolution within Artemisia is poorly understood. Here, we aimed to reconstruct a well-resolved phylogeny for global Artemisia via a phylogenomic approach, to infer the evolutionary patterns of its key morphological characters and to update its circumscription and infrageneric taxonomy.

METHODS

We sampled 228 species (258 samples) of Artemisia and its allies from both fresh and herbarium collections, covering all the subgenera and its main geographical areas, and conducted a phylogenomic analysis based on nuclear single nucleotide polymorphisms (SNPs) obtained from genome skimming data. Based on the phylogenetic framework, we inferred the possible evolutionary patterns of six key morphological characters widely used in its previous taxonomy.

KEY RESULTS

The genus Kaschgaria was revealed to be nested in Artemisia with strong support. A well-resolved phylogeny of Artemisia consisting of eight highly supported clades was recovered, two of which were identified for the first time. Most of the previously recognized subgenera were not supported as monophyletic. Evolutionary inferences based on the six morphological characters showed that different states of these characters originated independently more than once.

CONCLUSIONS

The circumscription of Artemisia is enlarged to include the genus Kaschgaria. The morphological characters traditionally used for the infrageneric taxonomy of Artemisia do not match the new phylogenetic tree. They experienced a more complex evolutionary history than previously thought. We propose a revised infrageneric taxonomy of the newly circumscribed Artemisia, with eight recognized subgenera to accommodate the new results.

Plastid phylogenomics and molecular evolution of Thismiaceae (Dioscoreales)

PREMISE

Species in Thismiaceae can no longer photosynthesize and instead obtain carbon from soil fungi. Here we infer Thismiaceae phylogeny using plastid genome data and characterize the molecular evolution of this genome.

METHODS

We assembled five Thismiaceae plastid genomes from genome skimming data, adding to previously published data for phylogenomic inference. We investigated plastid-genome structural changes, considering locally colinear blocks (LCBs). We also characterized possible shifts in selection pressure in retained genes by considering changes in the ratio of nonsynonymous to synonymous changes (ω).

RESULTS

Thismiaceae experienced two major pulses of gene loss around the early diversification of the family, with subsequent scattered gene losses across descendent lineages. In addition to massive size reduction, Thismiaceae plastid genomes experienced occasional inversions, and there were likely two independent losses of the plastid inverted repeat (IR) region. Retained plastid genes remain under generally strong purifying selection (ω << 1), with significant and sporadic weakening or strengthening in several instances. The bifunctional trnE-UUC gene of Thismia huangii may retain a secondary role in heme biosynthesis, despite a probable loss of functionality in protein translation. Several cis-spliced group IIA introns have been retained, despite the loss of the plastid intron maturase, matK.

CONCLUSIONS

We infer that most gene losses in Thismiaceae occurred early and rapidly, following the initial loss of photosynthesis in its stem lineage. As a species-rich, fully mycoheterotrophic lineage, Thismiaceae provide a model system for uncovering the unique and divergent ways in which plastid genomes evolve in heterotrophic plants.

Phylogeny and evolution of Asparagaceae subfamily Nolinoideae: new insights from plastid phylogenomics

BACKGROUND AND AIMS

Asparagaceae subfamily Nolinoideae is an economically important plant group, but the deep relationships and evolutionary history of the lineage remain poorly understood. Based on a large data set including 37 newly sequenced samples and publicly available plastomes, this study aims to better resolve the inter-tribal relationships of Nolinoideae, and to rigorously examine the tribe-level monophyly of Convallarieae, Ophiopogoneae and Polygonateae.

METHODS

Maximum likelihood (ML) and Bayesian inference (BI) methods were used to infer phylogenetic relationships of Nolinoideae at the genus level and above. The diversification history of Nolinoideae was explored using molecular dating.

KEY RESULTS

Both ML and BI analyses identically recovered five clades within Nolinoideae, respectively corresponding to Dracaeneae + Rusceae, Polygonateae + Theropogon, Ophiopogoneae, Nolineae, and Convallarieae excluding Theropogon, and most deep nodes were well supported. As Theropogon was embedded in Polygonateae, the plastome phylogeny failed to resolve Convallarieae and Polygonateae as reciprocally monophyletic. Divergence time estimation showed that the origins of most Nolinoideae genera were dated to the Miocene and Pliocene. The youthfulness of Nolinoideae genera is well represented in the three herbaceous tribes (Convallarieae, Ophiopogoneae and Polygonateae) chiefly distributed in temperate areas of the Northern Hemisphere, as the median stem ages of all 14 genera currently belonging to them were estimated at <12.37 Ma.

CONCLUSIONS

This study recovered a robust backbone phylogeny, providing new insights for better understanding the evolution and classification of Nolinoideae. Compared with the deep relationships recovered by a previous study based on transcriptomic data, our data suggest that ancient hybridization or incomplete lineage sorting may have occurred in the early diversification of Nolinoideae. Our findings will provide important reference for further study of the evolutionary complexity of Nolinoideae using nuclear genomic data. The recent origin of these herbaceous genera currently belonging to Convallarieae, Ophiopogoneae and Polygonateae provides new evidence to support the hypothesis that the global expansion of temperate habitats caused by the climate cooling over the past 15 million years may have dramatically driven lineage diversification and speciation in the Northern Hemisphere temperate flora.

Mitochondrial genomic data are effective at placing mycoheterotrophic lineages in plant phylogeny

Fully mycoheterotrophic plants can be difficult to place in plant phylogeny due to elevated substitution rates associated with photosynthesis loss. This potentially limits the effectiveness of downstream analyses of mycoheterotrophy that depend on accurate phylogenetic inference. Although mitochondrial genomic data sets are rarely used in plant phylogenetics, theory predicts that they should be resilient to long-branch artefacts, thanks to their generally slow evolution, coupled with limited rate elevation in heterotrophs. We examined the utility of mitochondrial genomes for resolving contentious higher-order placements of mycoheterotrophic lineages in two test cases: monocots (focusing on Dioscoreales) and Ericaceae. We find Thismiaceae to be distantly related to Burmanniaceae in the monocot order Dioscoreales, conflicting with current classification schemes based on few gene data sets. We confirm that the unusual Afrothismia is related to Taccaceae-Thismiaceae, with a corresponding independent loss of photosynthesis. In Ericaceae we recovered the first well supported relationships among its five major lineages: mycoheterotrophic Ericaceae are not monophyletic, as pyroloids are inferred to be sister to core Ericaceae, and monotropoids to arbutoids. Genes recovered from mitochondrial genomes collectively resolved previously ambiguous mycoheterotroph higher-order relationships. We propose that mitochondrial genomic data should be considered in standardised gene panels for inferring overall plant phylogeny.

Phylogenomic resolution of order- and family-level monocot relationships using 602 single-copy nuclear genes and 1375 BUSCO genes

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

Dark-colored Oncocyclus irises in Israel analyzed by AFLP, whole chloroplast genome sequencing and species distribution modeling

The Haynei is one of seven species aggregates (clusters of species having similar flower morphology) recognized in section Oncocyclus of genus Iris. This aggregate, characterized by dark-colored flowers, is represented by six species in Israel and adjacent Jordan. There is, however, no knowledge of the genetic relationship of these species making verification of their taxonomic status impossible. We investigated genetic variation in this group using analysis of whole chloroplast genomes and amplified fragment length polymorphism (AFLP). We also used species distribution modeling (SDM) to predict species ranges under current climatic conditions. We found some population groups within the currently recognized species of section Oncocyclus to represent dramatically different genetic entities which devaluates a general trend of merging many previously recognized species of section Oncocyclus based on their flower morphology. Despite the importance of homoploid hybridization in this group’s evolution and some apparently sporadically happening inter-specific gene flow, the main evolutionary forces in Oncocyclus appear to be vicariance and spatial isolation. Our findings suggest that some of the currently recognized species in section Oncocyclus need revision. A revision must be based on genetic analyses allowing the reconstruction of ancestry and recognition of the importance of vicariance and spatial isolation in the evolution of this group. The implications of the present findings for conservation are discussed.

Comparative and phylogenetic analyses of eleven complete chloroplast genomes of Dipterocarpoideae

BACKGROUND

In South-east Asia, Dipterocarpoideae is predominant in most mature forest communities, comprising around 20% of all trees. As large quantity and high quality wood are produced in many species, Dipterocarpoideae plants are the most important and valuable source in the timber market. The d-borneol is one of the essential oil components from Dipterocarpoideae (for example, Dryobalanops aromatica or Dipterocarpus turbinatus) and it is also an important traditional Chinese medicine (TCM) formulation known as "Bingpian" in Chinese, with antibacterial, analgesic and anti-inflammatory effects and can enhance anticancer efficiency.

METHODS

In this study, we analyzed 20 chloroplast (cp) genomes characteristics of Dipterocarpoideae, including eleven newly reported genomes and nine cp genomes previously published elsewhere, then we explored the chloroplast genomic features, inverted repeats contraction and expansion, codon usage, amino acid frequency, the repeat sequences and selective pressure analyses. At last, we constructed phylogenetic relationships of Dipterocarpoideae and found the potential barcoding loci.

RESULTS

The cp genome of this subfamily has a typical quadripartite structure and maintains a high degree of consistency among species. There were slightly more tandem repeats in cp genomes of Dipterocarpus and Vatica, and the psbH gene was subjected to positive selection in the common ancestor of all the 20 species of Dipterocarpoideae compared with three outgroups. Phylogenetic tree showed that genus Shorea was not a monophyletic group, some Shorea species and genus Parashorea are placed in one clade. In addition, the rpoC2 gene can be used as a potential marker to achieve accurate and rapid species identification in subfamily Dipterocarpoideae.

CONCLUSIONS

Dipterocarpoideae had similar cp genomic features and psbM, rbcL, psbH may function in the growth of Dipterocarpoideae. Phylogenetic analysis suggested new taxon treatment is needed for this subfamily indentification. In addition, rpoC2 is potential to be a barcoding gene to TCM distinguish.

The Molecular Basis of Kale Domestication: Transcriptional Profiling of Developing Leaves Provides New Insights Into the Evolution of a Brassica oleracea Vegetative Morphotype

Morphotypes of Brassica oleracea are the result of a dynamic interaction between genes that regulate the transition between vegetative and reproductive stages and those that regulate leaf morphology and plant architecture. In kales, ornate leaves, extended vegetative phase, and nutritional quality are some of the characters potentially selected by humans during domestication. We used a combination of developmental studies and transcriptomics to understand the vegetative domestication syndrome of kale. To identify candidate genes that are responsible for the evolution of domestic kale, we searched for transcriptome-wide differences among three vegetative B. oleracea morphotypes. RNA-seq experiments were used to understand the global pattern of expressed genes during a mixture of stages at one time in kale, cabbage, and the rapid cycling kale line TO1000. We identified gene expression patterns that differ among morphotypes and estimate the contribution of morphotype-specific gene expression that sets kale apart (3958 differentially expressed genes). Differentially expressed genes that regulate the vegetative to reproductive transition were abundant in all morphotypes. Genes involved in leaf morphology, plant architecture, defense, and nutrition were differentially expressed in kale. This allowed us to identify a set of candidate genes we suggest may be important in the kale domestication syndrome. Understanding candidate genes responsible for kale domestication is of importance to ultimately improve Cole crop production.

Complete chloroplast genomes shed light on phylogenetic relationships, divergence time, and biogeography of Allioideae (Amaryllidaceae)

Allioideae includes economically important bulb crops such as garlic, onion, leeks, and some ornamental plants in Amaryllidaceae. Here, we reported the complete chloroplast genome (cpDNA) sequences of 17 species of Allioideae, five of Amaryllidoideae, and one of Agapanthoideae. These cpDNA sequences represent 80 protein-coding, 30 tRNA, and four rRNA genes, and range from 151,808 to 159,998 bp in length. Loss and pseudogenization of multiple genes (i.e., rps2, infA, and rpl22) appear to have occurred multiple times during the evolution of Alloideae. Additionally, eight mutation hotspots, including rps15-ycf1, rps16-trnQ-UUG, petG-trnW-CCA, psbA upstream, rpl32-trnL-UAG, ycf1, rpl22, matK, and ndhF, were identified in the studied Allium species. Additionally, we present the first phylogenomic analysis among the four tribes of Allioideae based on 74 cpDNA coding regions of 21 species of Allioideae, five species of Amaryllidoideae, one species of Agapanthoideae, and five species representing selected members of Asparagales. Our molecular phylogenomic results strongly support the monophyly of Allioideae, which is sister to Amaryllioideae. Within Allioideae, Tulbaghieae was sister to Gilliesieae-Leucocoryneae whereas Allieae was sister to the clade of Tulbaghieae- Gilliesieae-Leucocoryneae. Molecular dating analyses revealed the crown age of Allioideae in the Eocene (40.1 mya) followed by differentiation of Allieae in the early Miocene (21.3 mya). The split of Gilliesieae from Leucocoryneae was estimated at 16.5 mya. Biogeographic reconstruction suggests an African origin for Allioideae and subsequent spread to Eurasia during the middle Eocene. Cool and arid conditions during the late Eocene led to isolation between African and Eurasian species. African Allioideae may have diverged to South American taxa in the late Oligocene. Rather than vicariance, long-distance dispersal is the most likely explanation for intercontinental distribution of African and South American Allioideae species.

Phylogenomics With Hyb-Seq Unravels Korean Hosta Evolution

The genus Hosta (Agavoideae and Asparagaceae) is one of the most popular landscaping and ornamental plants native to temperate East Asia. Their popularity has led to extensive hybridization to develop various cultivars. However, their long history of hybridization, cultivation, and selection has brought about taxonomic confusion in the Hosta species delimitation along with their indistinguishable morphology. Here, we conducted the first broad phylogenetic analyses of Hosta species based on the most comprehensive genomic data set to date. To do so, we captured 246 nuclear gene sequences and plastomes from 55 accessions of Korean Hosta species using the Hyb-Seq method. As a result, this study provides the following novel and significant findings: (1) phylogenetic analyses of the captured sequences retrieved six species of Hosta in South Korea compared to five to eleven species based on the previous studies, (2) their phylogenetic relationships suggested that the large genome size was ancestral and the diversification of Korean Hosta species was accompanied by decreases in genome sizes, (3) comparison between nuclear genes and plastome revealed several introgressive hybridization events between Hosta species, and (4) divergence times estimated here showed that Hosta diverged 35.59 million years ago, while Korean Hosta species rapidly diversified during the late Miocene. Last, we explored whether these genomic data could be used to infer the origin of cultivars. In summary, this study provides the most comprehensive genomic resources to be used in phylogenetic, population, and conservation studies of Hosta, as well as for unraveling the origin of many cultivars.

The Mitochondrial Genome of Eleusine indica and Characterization of Gene Content within Poaceae

Plant mitochondrial (mt) genome assembly provides baseline data on size, structure, and gene content, but resolving the sequence of these large and complex organelle genomes remains challenging due to fragmentation, frequent recombination, and transfers of DNA from neighboring plastids. The mt genome for Eleusine indica (Poaceae: goosegrass) is comprehensibly analyzed here, providing key reference data for an economically significant invasive species that is also the maternal parent of the allotetraploid crop Finger millet (Eleusine coracana). The assembled E. indica genome contains 33 protein coding genes, 6 rRNA subunits, 24 tRNA, 8 large repetitive regions 15 kb of transposable elements across a total of 520,691 bp. Evidence of RNA editing and loss of rpl2, rpl5, rps14, rps11, sdh4, and sdh3 genes is evaluated in the context of an updated survey of mt genomic gene content across the grasses through an analysis of publicly available data. Hypothesized patterns of Poaceae mt gene loss are examined in a phylogenetic context to clarify timing, showing that rpl2 was transferred to the nucleus from the mitochondrion prior to the origin of the PACMAD clade.

Evolution of Lomandroideae: Multiple origins of polyploidy and biome occupancy in Australia

Asparagaceae: Lomandroideae are a species-rich and economically important subfamily in the monocot order Asparagales, with a center of diversity in Australia. Lomandroideae are ecologically diverse, occupying mesic and arid biomes in Australia and possessing an array of key traits, including sexual dimorphism, storage organs and polyploidy that are potentially adaptive for survival in seasonally arid and fire-dependent habitats. The Lomandroideae phylogeny was reconstructed using maximum likelihood and Bayesian inference criteria, based on plastome data from genome-skimming to infer relationships. A fossil-calibrated chronogram provided a temporal framework for understanding trait transitions. Ancestral state reconstructions and phylogenetic comparative trait correlation analyses provided insights into the evolutionary and ecological drivers associated with Lomandroideae diversification. Lomandroideae diverged from the other Asparagaceae ca. 56.61 million years ago (95% highest posterior density values 70.31-45.34 million years) and the major lineages diversified since the Oligocene. The most recent common ancestor of the clade likely occupied the mesic biome, was hermaphroditic and geophytic. Biome occupancy transitions were correlated with polyploidy and the presence of storage roots. Polyploidy potentially serves as an "enabler" trait, generating novel phenotypes, which may confer tolerance to climatic ranges and soil conditions putatively required for expansion into and occupation of new arid biomes. Storage roots, as a key factor driving biome transitions, may have been associated with fire rather than with aridification events in the Australian flora. This study contributes significantly to our understanding of biome evolution by identifying polyploidy and storage organs as key factors associated with transitions in biome occupancy in this lineage.

The Treasure Vault Can be Opened: Large-Scale Genome Skimming Works Well Using Herbarium and Silica Gel Dried Material

Genome skimming has the potential for generating large data sets for DNA barcoding and wider biodiversity genomic studies, particularly via the assembly and annotation of full chloroplast (cpDNA) and nuclear ribosomal DNA (nrDNA) sequences. We compare the success of genome skims of 2051 herbarium specimens from Norway/Polar regions with 4604 freshly collected, silica gel dried specimens mainly from the European Alps and the Carpathians. Overall, we were able to assemble the full chloroplast genome for 67% of the samples and the full nrDNA cluster for 86%. Average insert length, cover and full cpDNA and rDNA assembly were considerably higher for silica gel dried than herbarium-preserved material. However, complete plastid genomes were still assembled for 54% of herbarium samples compared to 70% of silica dried samples. Moreover, there was comparable recovery of coding genes from both tissue sources (121 for silica gel dried and 118 for herbarium material) and only minor differences in assembly success of standard barcodes between silica dried (89% ITS2, 96% matK and rbcL) and herbarium material (87% ITS2, 98% matK and rbcL). The success rate was > 90% for all three markers in 1034 of 1036 genera in 160 families, and only Boraginaceae worked poorly, with 7 genera failing. Our study shows that large-scale genome skims are feasible and work well across most of the land plant families and genera we tested, independently of material type. It is therefore an efficient method for increasing the availability of plant biodiversity genomic data to support a multitude of downstream applications.

Complete Chloroplast Genomes of Chlorophytum comosum and Chlorophytum gallabatense: Genome Structures, Comparative and Phylogenetic Analysis

The genus Chlorophytum includes many economically important species well-known for medicinal, ornamental, and horticultural values. However, to date, few molecular genomic resources have been reported for this genus. Therefore, there is limited knowledge of phylogenetic studies, and the available chloroplast (cp) genome of Chlorophytum (C. rhizopendulum) does not provide enough information on this genus. In this study, we present genomic resources for C. comosum and C. gallabatense, which had lengths of 154,248 and 154,154 base pairs (bp), respectively. They had a pair of inverted repeats (IRa and IRb) of 26,114 and 26,254 bp each in size, separating the large single-copy (LSC) region of 84,004 and 83,686 bp from the small single-copy (SSC) region of 18,016 and 17,960 bp in C. comosum and C. gallabatense, respectively. There were 112 distinct genes in each cp genome, which were comprised of 78 protein-coding genes, 30 tRNA genes, and four rRNA genes. The comparative analysis with five other selected species displayed a generally high level of sequence resemblance in structural organization, gene content, and arrangement. Additionally, the phylogenetic analysis confirmed the previous phylogeny and produced a phylogenetic tree with similar topology. It showed that the Chlorophytum species (C. comosum, C. gallabatense and C. rhizopendulum) were clustered together in the same clade with a closer relationship than other plants to the Anthericum ramosum. This research, therefore, presents valuable records for further molecular evolutionary and phylogenetic studies which help to fill the gap in genomic resources and resolve the taxonomic complexes of the genus.

Museomics Clarifies the Classification of Aloidendron (Asphodelaceae), the Iconic African Tree Aloes

Arborescent succulent plants are regarded as keystone and indicator species in desert ecosystems due to their large stature and long lifespans. Tree aloes, the genus Aloidendron, are icons of the southern African deserts yet have proved elusive subjects due to the difficulty of obtaining material of known provenance for comparative study. Consequently, evolutionary relationships among representatives of the unusual arborescent life form have remained unclear until now. We used a museomics approach to overcome this challenge. Chloroplast genomes of six Aloidendron species and 12 other members of Asphodelaceae were sequenced from modern living collections and herbarium specimens, including the type specimens of all but two Aloidendron species, the earliest of which was collected 130 years ago. Maximum-likelihood trees estimated from full chloroplast genomes and the nuclear internal transcribed spacer (ITS) region show that Aloidendron sabaeum, from the Arabian Peninsula, is nested within Aloe while the Madagascar endemic Aloestrela suzannae is most closely related to the Somalian Aloidendron eminens. We observed phylogenetic conflicts between the plastid and nuclear topologies, which may be indicative of recurrent hybridisation or incomplete lineage sorting events in Aloe and in Aloidendron. Comparing species ecology in the context provided by our phylogeny suggests that habitat preference to either xeric deserts or humid forests/thickets evolved repeatedly in Aloidendron. Our findings demonstrate the value of botanical collections for the study and classification of taxonomically challenging succulent plants.

Cytogenetic and genomic organization analyses of chloroplast DNA invasions in the nuclear genome of Asparagus officinalis L. provides signatures of evolutionary complexity and informativity in sex chromosome evolution

BACKGROUND

The transfer of chloroplast DNA into nuclear genome is a common process in plants. These transfers form nuclear integrants of plastid DNAs (NUPTs), which are thought to be driving forces in genome evolution, including sex chromosome evolution. In this study, NUPTs in the genome of a dioecious plant Asparagus officinalis L. were systematically analyzed, in order to investigate the characteristics of NUPTs in the nuclear genome and the relationship between NUPTs and sex chromosome evolution in this species.

RESULTS

A total of 3155 NUPT insertions were detected, and they represented approximated 0.06% of the nuclear genome. About 45% of the NUPTs were organized in clusters. These clusters were derived from various evolutionary events. The Y chromosome contained the highest number and largest proportion of NUPTs, suggesting more accumulation of NUPTs on sex chromosomes. NUPTs were distributed widely in all of the chromosomes, and some regions preferred these insertions. The highest density of NUPTs was found in a 47 kb region in the Y chromosome; more than 75% of this region was occupied by NUPTs. Further cytogenetic and sequence alignment analysis revealed that this region was likely the centromeric region of the sex chromosomes. On the other hand, the male-specific region of the Y chromosome (MSY) and the adjacent regions did not have NUPT insertions.

CONCLUSIONS

These results indicated that NUPTs were involved in shaping the genome of A. officinalis through complicated process. NUPTs may play important roles in the centromere shaping of the sex chromosomes of A. officinalis, but were not implicated in MSY formation.

Phylogeny of Chinese Allium Species in Section Daghestanica and Adaptive Evolution of Allium (Amaryllidaceae, Allioideae) Species Revealed by the Chloroplast Complete Genome

The genus Allium (Amaryllidaceae, Allioideae) is one of the largest monocotyledonous genera and it includes many economically important crops that are cultivated for consumption or medicinal uses. Recent advances in molecular phylogenetics have revolutionized our understanding of Allium taxonomy and evolution. However, the phylogenetic relationships in some Allium sections (such as the Allium section Daghestanica) and the genetic bases of adaptative evolution, remain poorly understood. Here, we newly assembled six chloroplast genomes from Chinese endemic species in Allium section Daghestanica and by combining these genomes with another 35 allied species, we performed a series of analyses including genome structure, GC content, species pairwise Ka/Ks ratios, and the SSR component, nucleotide diversity and codon usage. Positively selected genes (PSGs) were detected in the Allium lineage using the branch-site model. Comparison analysis of Bayesian and ML phylogeny on CCG (complete chloroplast genome), SCG (single copy genes) and CDS (coding DNA sequences) produced a well-resolved phylogeny of Allioideae plastid lineages, which illustrated several novel relationships with the section Daghestanica. In addition, six species in section Daghestanica showed highly conserved structures. The GC content and the GC3s content in Allioideae species exhibited lower values than studied non-Allioideae species, along with elevated pairwise Ka/Ks ratios. The rps2 gene was lost in all examined Allioideae species, and 10 genes with significant posterior probabilities for codon sites were identified in the positive selection analysis, seven of them are associated with photosynthesis. Our study uncovered a new species relationship in section Daghestanica and suggested that the selective pressure has played an important role in Allium adaptation and evolution, these results will facilitate our further understanding of evolution and adaptation of species in the genus Allium.

Extremely low levels of chloroplast genome sequence variability in Astelia pumila (Asteliaceae, Asparagales)

Astelia pumila (G.Forst.) Gaudich. (Asteliaceae, Asparagales) is a major element of West Patagonian cushion peat bog vegetation. With the aim to identify appropriate chloroplast markers for the use in a phylogeographic study, the complete chloroplast genomes of five A. pumila accessions from almost the entire geographical range of the species were assembled and screened for variable positions. The chloroplast genome sequence was obtained via a mapping approach, using Eustrephus latifolius (Asparagaceae) as a reference. The chloroplast genome of A. pumila varies in length from 158,215 bp to 158,221 bp, containing a large single copy region of 85,981–85,983 bp, a small single copy region of 18,182–18,186 bp and two inverted repeats of 27,026 bp. Genome annotation predicted a total of 113 genes, including 30 tRNA and four rRNA genes. Sequence comparisons revealed a very low degree of intraspecific genetic variability, as only 37 variable sites (18 indels, 18 single nucleotide polymorphisms, one 3-bp mutation)—most of them autapomorphies—were found among the five assembled chloroplast genomes. A Maximum Likelihood analysis, based on whole chloroplast genome sequences of several Asparagales accessions representing six of the currently recognized 14 families (sensu APG IV), confirmed the phylogenetic position of A. pumila. The chloroplast genome of A. pumila is the first to be reported for a member of the astelioid clade (14 genera with c. 215 species), a basally branching group within Asparagales.

Monocot plastid phylogenomics, timeline, net rates of species diversification, the power of multi-gene analyses, and a functional model for the origin of monocots

PREMISE OF THE STUDY

We present the first plastome phylogeny encompassing all 77 monocot families, estimate branch support, and infer monocot-wide divergence times and rates of species diversification.

METHODS

We conducted maximum likelihood analyses of phylogeny and BAMM studies of diversification rates based on 77 plastid genes across 545 monocots and 22 outgroups. We quantified how branch support and ascertainment vary with gene number, branch length, and branch depth.

KEY RESULTS

Phylogenomic analyses shift the placement of 16 families in relation to earlier studies based on four plastid genes, add seven families, date the divergence between monocots and eudicots+Ceratophyllum at 136 Mya, successfully place all mycoheterotrophic taxa examined, and support recognizing Taccaceae and Thismiaceae as separate families and Arecales and Dasypogonales as separate orders. Only 45% of interfamilial divergences occurred after the Cretaceous. Net species diversification underwent four large-scale accelerations in PACMAD-BOP Poaceae, Asparagales sister to Doryanthaceae, Orchidoideae-Epidendroideae, and Araceae sister to Lemnoideae, each associated with specific ecological/morphological shifts. Branch ascertainment and support across monocots increase with gene number and branch length, and decrease with relative branch depth. Analysis of entire plastomes in Zingiberales quantifies the importance of non-coding regions in identifying and supporting short, deep branches.

CONCLUSIONS

We provide the first resolved, well-supported monocot phylogeny and timeline spanning all families, and quantify the significant contribution of plastome-scale data to resolving short, deep branches. We outline a new functional model for the evolution of monocots and their diagnostic morphological traits from submersed aquatic ancestors, supported by convergent evolution of many of these traits in aquatic Hydatellaceae (Nymphaeales).

Comparative analysis of the complete plastomes of garlic Allium sativum and bulb onion Allium cepa

Sequencing and comparative characterization of plant plastid genomes, or plastomes, is an important tool for modern phylogenetic and taxonomic studies, as well as for understanding the plastome evolution. The genusAlliumL. (family Amaryllidaceae) incorporates more than 900 species, includes economically signifi­cant vegetable crops such as garlicA. sativum, onionA. cepa, leekA. porrum, etc. In this work, the plastome of garlicA. sativumhas been completely sequenced. TheA. sativumplastome is 153172 bp in size. It consists of a large unique (LSC, 82035 bp) and small unique (SSC, 18015 bp) copies, separated by inverted repeats (IRa and IRb) of 26561 bp each. In the garlic plastome, 134 genes have been annotated: 82 protein-coding genes, 38 tRNA genes, 8 rRNA genes, and 6 pseudogenes. Comparative analysis ofA. sativumandA. cepaplastomes reveals differences in the sizes of structural elements and spacers at the inverted repeat bound­aries. The total numbers of genes inA. sativumandA. cepaare the same, but the gene composition is dif­ferent: therpl22gene is functional inA. sativum, being a pseudogene inA. cepa; conversely, therps16gene is a pseudogene inA. sativumand a protein-coding gene inA. cepa. In theA. sativumandA. cepaplastomes, 32 SSR sequences have been identified. More than half of them are dinucleotides, and the remaining are tetra-, penta-, and hexanucleotides at the same time, trinucleotides were absent. The compared plastomes differ in the numbers of certain SSRs, and some are present in only one of the species.

Phylogenomics and barcoding of Panax: toward the identification of ginseng species

BACKGROUND

The economic value of ginseng in the global medicinal plant trade is estimated to be in excess of US$2.1 billion. At the same time, the evolutionary placement of ginseng (Panax ginseng) and the complex evolutionary history of the genus is poorly understood despite several molecular phylogenetic studies. In this study, we use a full plastome phylogenomic framework to resolve relationships in Panax and to identify molecular markers for species discrimination.

RESULTS

We used high-throughput sequencing of MBD2-Fc fractionated Panax DNA to supplement publicly available plastid genomes to create a phylogeny based on fully assembled and annotated plastid genomes from 60 accessions of 8 species. The plastome phylogeny based on a 163 kbp matrix resolves the sister relationship of Panax ginseng with P. quinquefolius. The closely related species P. vietnamensis is supported as sister of P. japonicus. The plastome matrix also shows that the markers trnC-rps16, trnS-trnG, and trnE-trnM could be used for unambiguous molecular identification of all the represented species in the genus.

CONCLUSIONS

MBD2 depletion reduces the cost of plastome sequencing, which makes it a cost-effective alternative to Sanger sequencing based DNA barcoding for molecular identification. The plastome phylogeny provides a robust framework that can be used to study the evolution of morphological characters and biosynthesis pathways of ginsengosides for phylogenetic bioprospecting. Molecular identification of ginseng species is essential for authenticating ginseng in international trade and it provides an incentive for manufacturers to create authentic products with verified ingredients.

Practical considerations for plant phylogenomics

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

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.

Timing of rapid diversification and convergent origins of active pollination within Agavoideae (Asparagaceae)

PREMISE OF THE STUDY

Yucca species are ideal candidates for the study of coevolution due to the obligate mutualism they form with yucca moth pollinators (genera Tegeticula and Parategeticula). Yuccas are not the only species to exhibit a mutualism with yucca moths; the genus Hesperoyucca is pollinated by the California yucca moth (Tegeticula maculata). Relationships among yuccas, Hesperoyucca, and other members of subfamily Agavoideae are necessary to understand the evolution of this unique pollination syndrome. Here, we investigate evolutionary relationships of yuccas and closely related genera looking at the timing and origin of yucca moth pollination.

METHODS

In this study, we sequenced the chloroplast genomes of 20 species in the subfamily Agavoideae (Asparagaceae) and three confamilial outgroup taxa to resolve intergeneric phylogenetic relationships of Agavoideae. We estimated divergence times using protein-coding genes from 67 chloroplast genomes sampled across monocots to determine the timing of the yucca moth pollination origin.

KEY RESULTS

We confidently resolved intergeneric relationships in Agavoideae, demonstrating the origin of the yucca-yucca moth mutualism on two distinct lineages that diverged 27 million years ago. Comparisons of Yucca and Hesperoyucca divergence time to those of yucca moths (Tegeticula and Parategeticula, Prodoxidae) indicate overlapping ages for the origin of pollinating behavior in the moths and pollination by yucca moths in the two plant lineages.

CONCLUSION

Whereas pollinating yucca moths have been shown to have a single origin within the Prodoxidae, there were independent acquisitions of active pollination on lineages leading to Yucca and Hesperoyucca within the Agavoideae.

Verdant: automated annotation, alignment and phylogenetic analysis of whole chloroplast genomes

Motivation

Chloroplast genomes are now produced in the hundreds for angiosperm phylogenetics projects, but current methods for annotation, alignment and tree estimation still require some manual intervention reducing throughput and increasing analysis time for large chloroplast systematics projects.

Results

Verdant is a web-based software suite and database built to take advantage a novel annotation program, annoBTD. Using annoBTD, Verdant provides accurate annotation of chloroplast genomes without manual intervention. Subsequent alignment and tree estimation can incorporate newly annotated and publically available plastomes and can accommodate a large number of taxa. Verdant sharply reduces the time required for analysis of assembled chloroplast genomes and removes the need for pipelines and software on personal hardware.

Availability and Implementation

Verdant is available at: http://verdant.iplantcollaborative.org/plastidDB/. It is implemented in PHP, Perl, MySQL, Javascript, HTML and CSS with all major browsers supported.

Supplementary information

Supplementary data are available at Bioinformatics online.

Strongly asymmetric hybridization barriers shape the origin of a new polyploid species and its hybrid ancestor

PREMISE OF THE STUDY

Hybridization between diploids and tetraploids can lead to new allopolyploid species, often via a triploid intermediate. Viable triploids are often produced asymmetrically, with greater success observed for "maternal-excess" crosses where the mother has a higher ploidy than the father. Here we investigated the evolutionary origins of Mimulus peregrinus, an allohexaploid recently derived from the triploid M. ×robertsii, to determine whether reproductive asymmetry has shaped the formation of this new species.

METHODS

We used reciprocal crosses between the diploid (M. guttatus) and tetraploid (M. luteus) progenitors to determine the viability of triploid M. ×robertsii hybrids resulting from paternal- vs. maternal-excess crosses. To investigate whether experimental results predict patterns seen in the field, we performed parentage analyses comparing natural populations of M. peregrinus to its diploid, tetraploid, and triploid progenitors. Organellar sequences obtained from pre-existing genomic data, supplemented with additional genotyping was used to establish the maternal ancestry of multiple M. peregrinus and M. ×robertsii populations.

KEY RESULTS

We found strong evidence for asymmetric origins of M. peregrinus, but opposite to the common pattern, with paternal-excess crosses significantly more successful than maternal-excess crosses. These results successfully predicted hybrid formation in nature: 111 of 114 M. ×robertsii individuals, and 27 of 27 M. peregrinus, had an M. guttatus maternal haplotype.

CONCLUSION

This study, which includes the first Mimulus chloroplast genome assembly, demonstrates the utility of parentage analysis through genome skimming. We highlight the benefits of complementing genomic analyses with experimental approaches to understand asymmetry in allopolyploid speciation.

Multilocus phylogeny and phylogenomics of Eriochrysis P. Beauv. (Poaceae-Andropogoneae): Taxonomic implications and evidence of interspecific hybridization

Species delimitation is a vital issue concerning evolutionary biology and conservation of biodiversity. However, it is a challenging task for several reasons, including the low interspecies variability of markers currently used in phylogenetic reconstructions and the occurrence of reticulate evolution and polyploidy in many lineages of flowering plants. The first phylogeny of the grass genus Eriochrysis is presented here, focusing on the New World species, in order to examine its relationships to other genera of the subtribe Saccharinae/tribe Andropogoneae and to define the circumscriptions of its taxonomically complicated species. Molecular cloning and sequencing of five regions of four low-copy nuclear genes (apo1, d8, ep2-ex7 and ep2-ex8, kn1) were performed, as well as complete plastome sequencing. Trees were reconstructed using maximum parsimony, maximum likelihood, and Bayesian inference analyses. The present phylogenetic analyses indicate that Eriochrysis is monophyletic and the Old World E. pallida is sister to the New World species. Subtribe Saccharinae is polyphyletic, as is the genus Eulalia. Based on nuclear and plastome sequences plus morphology, we define the circumscriptions of the New World species of Eriochrysis: E. laxa is distinct from E. warmingiana, and E. villosa is distinct from E. cayennensis. Natural hybrids occur between E. laxa and E. villosa. The hybrids are probably tetraploids, based on the number of paralogues in the nuclear gene trees. This is the first record of a polyploid taxon in the genus Eriochrysis. Some incongruities between nuclear genes and plastome analyses were detected and are potentially caused by incomplete lineage sorting and/or ancient hybridization. The set of low-copy nuclear genes used in this study seems to be sufficient to resolve phylogenetic relationships and define the circumscriptions of other species complexes in the grass family and relatives, even in the presence of polyploidy and reticulate evolution. Complete plastome sequencing is also a promising tool for phylogenetic inference.

Plastid genomes reveal support for deep phylogenetic relationships and extensive rate variation among palms and other commelinid monocots

Despite progress based on multilocus, phylogenetic studies of the palms (order Arecales, family Arecaceae), uncertainty remains in resolution/support among major clades and for the placement of the palms among the commelinid monocots. Palms and related commelinids represent a classic case of substitution rate heterogeneity that has not been investigated in the genomic era. To address questions of relationships, support and rate variation among palms and commelinid relatives, 39 plastomes representing the palms and related family Dasypogonaceae were generated via genome skimming and integrated within a monocot-wide matrix for phylogenetic and molecular evolutionary analyses. Support was strong for 'deep' relationships among the commelinid orders, among the five palm subfamilies, and among tribes of the subfamily Coryphoideae. Additionally, there was extreme heterogeneity in the plastid substitution rates across the commelinid orders indicated by model based analyses, with c. 22 rate shifts, and significant departure from a global clock. To date, this study represents the most comprehensively sampled matrix of plastomes assembled for monocot angiosperms, providing genome-scale support for phylogenetic relationships of monocot angiosperms, and lays the phylogenetic groundwork for comparative analyses of the drivers and correlates of such drastic differences in substitution rates across a diverse and significant clade.

Complete Chloroplast Genome Sequence of Aquilaria sinensis (Lour.) Gilg and Evolution Analysis within the Malvales Order

Aquilaria sinensis (Lour.) Gilg is an important medicinal woody plant producing agarwood, which is widely used in traditional Chinese medicine. High-throughput sequencing of chloroplast (cp) genomes enhanced the understanding about evolutionary relationships within plant families. In this study, we determined the complete cp genome sequences for A. sinensis. The size of the A. sinensis cp genome was 159,565 bp. This genome included a large single-copy region of 87,482 bp, a small single-copy region of 19,857 bp, and a pair of inverted repeats (IRa and IRb) of 26,113 bp each. The GC content of the genome was 37.11%. The A. sinensis cp genome encoded 113 functional genes, including 82 protein-coding genes, 27 tRNA genes, and 4 rRNA genes. Seven genes were duplicated in the protein-coding genes, whereas 11 genes were duplicated in the RNA genes. A total of 45 polymorphic simple-sequence repeat loci and 60 pairs of large repeats were identified. Most simple-sequence repeats were located in the noncoding sections of the large single-copy/small single-copy region and exhibited high A/T content. Moreover, 33 pairs of large repeat sequences were located in the protein-coding genes, whereas 27 pairs were located in the intergenic regions. Aquilaria sinensis cp genome bias ended with A/T on the basis of codon usage. The distribution of codon usage in A. sinensis cp genome was most similar to that in the Gonystylus bancanus cp genome. Comparative results of 82 protein-coding genes from 29 species of cp genomes demonstrated that A. sinensis was a sister species to G. bancanus within the Malvales order. Aquilaria sinensis cp genome presented the highest sequence similarity of >90% with the G. bancanus cp genome by using CGView Comparison Tool. This finding strongly supports the placement of A. sinensis as a sister to G. bancanus within the Malvales order. The complete A. sinensis cp genome information will be highly beneficial for further studies on this traditional medicinal plant. Moreover, the results will enhance our understanding about the evolution of cp genomes of the Malvales order, particularly with regard to the role of A. sinensis in plant systematics and evolution.

Seven New Complete Plastome Sequences Reveal Rampant Independent Loss of the ndh Gene Family across Orchids and Associated Instability of the Inverted Repeat/Small Single-Copy Region Boundaries

Earlier research has revealed that the ndh loci have been pseudogenized, truncated, or deleted from most orchid plastomes sequenced to date, including in all available plastomes of the two most species-rich subfamilies, Orchidoideae and Epidendroideae. This study sought to resolve deeper-level phylogenetic relationships among major orchid groups and to refine the history of gene loss in the ndh loci across orchids. The complete plastomes of seven orchids, Oncidium sphacelatum (Epidendroideae), Masdevallia coccinea (Epidendroideae), Sobralia callosa (Epidendroideae), Sobralia aff. bouchei (Epidendroideae), Elleanthus sodiroi (Epidendroideae), Paphiopedilum armeniacum (Cypripedioideae), and Phragmipedium longifolium (Cypripedioideae) were sequenced and analyzed in conjunction with all other available orchid and monocot plastomes. Most ndh loci were found to be pseudogenized or lost in Oncidium, Paphiopedilum and Phragmipedium, but surprisingly, all ndh loci were found to retain full, intact reading frames in Sobralia, Elleanthus and Masdevallia. Character mapping suggests that the ndh genes were present in the common ancestor of orchids but have experienced independent, significant losses at least eight times across four subfamilies. In addition, ndhF gene loss was correlated with shifts in the position of the junction of the inverted repeat (IR) and small single-copy (SSC) regions. The Orchidaceae have unprecedented levels of homoplasy in ndh gene presence/absence, which may be correlated in part with the unusual life history of orchids. These results also suggest that ndhF plays a role in IR/SSC junction stability.

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.

SMRT sequencing only de novo assembly of the sugar beet (Beta vulgaris) chloroplast genome

BACKGROUND

Third generation sequencing methods, like SMRT (Single Molecule, Real-Time) sequencing developed by Pacific Biosciences, offer much longer read length in comparison to Next Generation Sequencing (NGS) methods. Hence, they are well suited for de novo- or re-sequencing projects. Sequences generated for these purposes will not only contain reads originating from the nuclear genome, but also a significant amount of reads originating from the organelles of the target organism. These reads are usually discarded but they can also be used for an assembly of organellar replicons. The long read length supports resolution of repetitive regions and repeats within the organelles genome which might be problematic when just using short read data. Additionally, SMRT sequencing is less influenced by GC rich areas and by long stretches of the same base.

RESULTS

We describe a workflow for a de novo assembly of the sugar beet (Beta vulgaris ssp. vulgaris) chloroplast genome sequence only based on data originating from a SMRT sequencing dataset targeted on its nuclear genome. We show that the data obtained from such an experiment are sufficient to create a high quality assembly with a higher reliability than assemblies derived from e.g. Illumina reads only. The chloroplast genome is especially challenging for de novo assembling as it contains two large inverted repeat (IR) regions. We also describe some limitations that still apply even though long reads are used for the assembly.

CONCLUSIONS

SMRT sequencing reads extracted from a dataset created for nuclear genome (re)sequencing can be used to obtain a high quality de novo assembly of the chloroplast of the sequenced organism. Even with a relatively small overall coverage for the nuclear genome it is possible to collect more than enough reads to generate a high quality assembly that outperforms short read based assemblies. However, even with long reads it is not always possible to clarify the order of elements of a chloroplast genome sequence reliantly which we could demonstrate with Fosmid End Sequences (FES) generated with Sanger technology. Nevertheless, this limitation also applies to short read sequencing data but is reached in this case at a much earlier stage during finishing.

Phylogeny and photosynthesis of the grass tribe Paniceae

PREMISE OF THE STUDY

The grass tribe Paniceae includes important food, forage, and bioenergy crops such as switchgrass, napiergrass, various millet species, and economically important weeds. Paniceae are also valuable for answering scientific and evolutionary questions about C4 photosynthetic evolution, drought tolerance, and spikelet variation. However, the phylogeny of the tribe remains incompletely resolved.

METHODS

Forty-five taxa were selected from across the tribe Paniceae and outgroups for genome survey sequencing (GSS). These data were used to build a phylogenetic tree of the Paniceae based on 102 markers (78 chloroplast, 22 mitochondrial, 2 nrDNA). Ancestral state reconstruction analyses were also performed within the Paniceae using both the traditional and two subtype classification systems to test hypotheses of C4 subtype evolution.

KEY RESULTS

The phylogenetic tree resolves many areas of the Paniceae with high support and provides insight into the origin and number of C4 evolution events within the tribe. The recovered phylogeny and ancestral state reconstructions support between four and seven independent origins of C4 photosynthesis within the tribe and indicate which species are potentially the closest C3 sister taxa of each of these events.

CONCLUSIONS

Although the sequence of evolutionary events that produced multiple C4 subtypes within the Paniceae remains undetermined, the results presented here are consistent with only a subset of currently proposed models. The species used in this study constitute a panel of C3 and C4 grasses that are suitable for further studies on C4 photosynthesis, bioenergy, food and forage crops, and various developmental features of the Paniceae.

Biodiversity comparison among phylogenetic diversity metrics and between three North American prairies

Protection of Earth's ecosystems requires identification of geographical areas of greatest biodiversity. Assessment of biodiversity begins with knowledge of the evolutionary histories of species in a geographic area. Multiple phylogenetic diversity (PD) metrics have been developed to describe biodiversity beyond species counts, but sufficient empirical studies, particularly at fine phylogenetic scales, have not been conducted to provide conservation planners with evidence for incorporating PD metrics into selection of priority regions. We review notable studies that are contributing to a growing database of empirical results, we report on the effect of using high-throughput sequencing to estimate the phylogenies used to calculate PD metrics, and we discuss difficulties in selecting appropriate diversity indices. We focused on two of the most speciose angiosperm families in prairies-Asteraceae and Fabaceae-and compared 12 PD metrics and four traditional measures of biodiversity between three North American prairie sites. The varying results from the literature and from the current data reveal the wide range of applications of PD metrics and the necessity for many more empirical studies. The accumulation of results from further investigations will eventually lead to a scientific understanding upon which conservation planners can make informed decisions about where to apply limited preservation funds.

Resolving relationships within the palm subfamily Arecoideae (Arecaceae) using plastid sequences derived from next-generation sequencing

PREMISE OF THE STUDY

Several studies have incorporated molecular and morphological data to study the phylogeny of the palms (Arecaceae), but some relationships within the family remain ambiguous-particularly those within Arecoideae, the most diverse subfamily including coconut and oil palm. Here, two next-generation, targeted plastid-enrichment methods were compared and used to elucidate Arecoideae phylogeny.

METHODS

Next-generation sequencing techniques were used to generate a plastid genome data set. Long range PCR and hybrid gene capture were used to enrich for chloroplast targets. Ten taxa were enriched using both methods for comparison. Chloroplast sequence data were generated for 31 representatives of the 14 Arecoideae tribes and five outgroup taxa. The phylogeny was reconstructed using maximum likelihood, maximum parsimony, and Bayesian analyses.

KEY RESULTS

Long range PCR and hybrid gene capture both enriched the plastid genome and provided similar sequencing coverage. Subfamily Arecoideae was resolved as monophyletic with tribe Chamaedoreeae as the earliest-diverging lineage, implying that the development of flowers in triads defines a synapomorphy for the Arecoideae clade excluding Chamaedoreeae. Three major clades within this group were recovered: Roystoneeae/Reinhardtieae/Cocoseae (RRC), Areceae/Euterpeae/Geonomateae/Leopoldinieae/Manicarieae/Pelagodoxeae (core arecoids), and Podococceae/Oranieae/Sclerospermeae (POS). An Areceae + Euterpeae clade was resolved within the core arecoids. The POS clade was sister to a RRC + core arecoids clade, implying a shared ancestral area in South America for these three clades.

CONCLUSIONS

The plastome phylogeny recovered here provides robust resolution of previously ambiguous studies and new insights into palm evolution.

Are differences in genomic data sets due to true biological variants or errors in genome assembly: an example from two chloroplast genomes

DNA sequencing has been revolutionized by the development of high-throughput sequencing technologies. Plummeting costs and the massive throughput capacities of second and third generation sequencing platforms have transformed many fields of biological research. Concurrently, new data processing pipelines made rapid de novo genome assemblies possible. However, high quality data are critically important for all investigations in the genomic era. We used chloroplast genomes of one Oryza species (O. australiensis) to compare differences in sequence quality: one genome (GU592209) was obtained through Illumina sequencing and reference-guided assembly and the other genome (KJ830774) was obtained via target enrichment libraries and shotgun sequencing. Based on the whole genome alignment, GU592209 was more similar to the reference genome (O. sativa: AY522330) with 99.2% sequence identity (SI value) compared with the 98.8% SI values in the KJ830774 genome; whereas the opposite result was obtained when the SI values in coding and noncoding regions of GU592209 and KJ830774 were compared. Additionally, the junctions of two single copies and repeat copies in the chloroplast genome exhibited differences. Phylogenetic analyses were conducted using these sequences, and the different data sets yielded dissimilar topologies: phylogenetic replacements of the two individuals were remarkably different based on whole genome sequencing or SNP data and insertions and deletions (indels) data. Thus, we concluded that the genomic composition of GU592209 was heterogeneous in coding and non-coding regions. These findings should impel biologists to carefully consider the quality of sequencing and assembly when working with next-generation data.