Phylogenomics using low-depth whole genome sequencing: A case study with the olive tribe

Phylogenomic data exploration with increased sampling provides new insights into the higher-level relationships of butterflies and moths (Lepidoptera)

A robust and stable phylogenetic framework is a fundamental goal of evolutionary biology. As the third largest insect order in the world following Coleoptera and Diptera, Lepidoptera (butterflies and moths) play a central role in almost every terrestrial ecosystem as indicators of environmental change and serve as important models for biologists exploring questions related to ecology and evolutionary biology. However, for such a charismatic insect group, the higher-level phylogenetic relationships among its superfamilies are still poorly resolved. Compared to earlier phylogenomic studies, we increased taxon sampling among Lepidoptera (37 superfamilies and 68 families containing 263 taxa) and acquired a series of large amino-acid datasets from 69,680 to 400,330 for phylogenomic reconstructions. Using these datasets, we explored the effect of different taxon sampling with significant increases in the number of included genes on tree topology by considering a series of systematic errors using maximum-likelihood (ML) and Bayesian inference (BI) methods. Moreover, we also tested the effectiveness in topology robustness among the three ML-based models. The results showed that taxon sampling is an important determinant in tree robustness of accurate lepidopteran phylogenetic estimation. Long-branch attraction (LBA) caused by site-wise heterogeneity is a significant source of bias giving rise to unstable positions of ditrysian groups in phylogenomic reconstruction. Phylogenetic inference showed the most comprehensive framework to reveal the relationships among lepidopteran superfamilies, and presented some newly relationships with strong supports (Papilionoidea was sister to Gelechioidea and Immoidea was sister to Galacticoidea, respectively), but limited by taxon sampling, the relationships within the species-rich and relatively rapid radiation Ditrysia and especially Apoditrysia remain poorly resolved, which need to increase taxon sampling for further phylogenomic reconstruction. The present study demonstrates that taxon sampling is an important determinant for an accurate lepidopteran tree of life and provides some essential insights for future lepidopteran phylogenomic studies.

Variant calling and genotyping accuracy of ddRAD-seq: Comparison with 20X WGS in layers

Whole Genome Sequencing (WGS) remains a costly or unsuitable method for routine genotyping of laying hens. Until now, breeding companies have been using or developing SNP chips. Nevertheless, alternatives methods based on sequencing have been developed. Among these, reduced representation sequencing approaches can offer sequencing quality and cost-effectiveness by reducing the genomic regions covered by sequencing. The aim of this study was to evaluate the ability of double digested Restriction site Associated DNA sequencing (ddRAD-seq) to identify and genotype SNPs in laying hens, by comparison with a presumed reliable WGS approach. Firstly, the sensitivity and precision of variant calling and the genotyping reliability of ddRADseq were determined. Next, the SNP Call Rate (CRSNP) and mean depth of sequencing per SNP (DPSNP) were compared between both methods. Finally, the effect of multiple combinations of thresholds for these parameters on genotyping reliability and amount of remaining SNPs in ddRAD-seq was studied. In raw form, the ddRAD-seq identified 349,497 SNPs evenly distributed on the genome with a CRSNP of 0.55, a DPSNP of 11X and a mean genotyping reliability rate per SNP of 80%. Considering genomic regions covered by expected enzymatic fragments (EFs), the sensitivity of the ddRAD-seq was estimated at 32.4% and its precision at 96.4%. The low CRSNP and DPSNP values were explained by the detection of SNPs outside the EFs theoretically generated by the ddRAD-seq protocol. Indeed, SNPs outside the EFs had significantly lower CRSNP (0.25) and DPSNP (1X) values than SNPs within the EFs (0.7 and 17X, resp.). The study demonstrated the relationship between CRSNP, DPSNP, genotyping reliability and the number of SNPs retained, to provide a decision-support tool for defining filtration thresholds. Severe quality control over ddRAD-seq data allowed to retain a minimum of 40% of the SNPs with a CcR of 98%. Then, ddRAD-seq was defined as a suitable method for variant calling and genotyping in layers.

The complete chloroplast genome sequence of Olea dioica Roxb, 1820 (Oleaceae)

Olea dioica Roxb, 1820 is a very important ethnomedicinal tree because of its medicinal properties and it belongs to the Oleaceae family. It is mainly distributed in evergreen and semi-evergreen forests. However, the chloroplast genome of O. dioica has not yet been reported. In this study, the chloroplast genome sequence of O. dioica was sequenced using next-generation sequencing technologies. The complete chloroplast genome of O. dioica was 155,138 bp in length (GenBank accession no. PP048999), comprising a large single-copy (LSC) region of 86,048 bp, a small single-copy (SSC) region of 17,816 bp, and two inverted repeat (IR) regions 25,637 bp each. The overall GC content was 37.8%. The complete chloroplast genome of O. dioica contains 131 complete genes, which are 88 protein-coding genes, 35 transfer RNA genes, and eight ribosomal RNA genes. A maximum-likelihood (ML) tree of O. dioica and 14 other species in the family Oleaceae suggested that O. dioica showed a close relationship with Olea brachiata.

Historical climate change and vicariance events contributed to the intercontinental disjunct distribution pattern of ash species (Fraxinus, Oleaceae)

The Northern Hemisphere temperate forests exhibit a disjunct distributional pattern in Europe, North America, and East Asia. Here, to reveal the promoter of intercontinental disjunct distribution, Fraxinus was used as a model organism to integrate abundant fossil evidence with high-resolution phylogenies in a phytogeographic analysis. We constructed a robust phylogenetic tree using genomic data, reconstructed the geographic ancestral areas, and evaluated the effect of incorporating fossil information on the reconstructed biogeographic history. The phylogenetic relationships of Fraxinus were highly resolved and divided into seven clades. Fraxinus originated in western North America during Eocene, and six intercontinental dispersal events and five intercontinental vicariance events were occured. Results suggest that climate change and vicariance contributed to the intercontinental disjunct distribution pattern of Fraxinus. Moreover, results highlight the necessity of integrating phylogenetic relationship and fossil to improve the reliability of inferred biogeographic events and our understanding of the processes underlying disjunct distributions.

A hemizygous supergene controls homomorphic and heteromorphic self-incompatibility systems in Oleaceae

Self-incompatibility (SI) has evolved independently multiple times and prevents self-fertilization in hermaphrodite angiosperms. Several groups of Oleaceae such as jasmines exhibit distylous flowers, with two compatibility groups each associated with a specific floral morph.1 Other Oleaceae species in the olive tribe have two compatibility groups without associated morphological variation.2,3,4,5 The genetic basis of both homomorphic and dimorphic SI systems in Oleaceae is unknown. By comparing genomic sequences of three olive subspecies (Olea europaea) belonging to the two compatibility groups, we first locate the genetic determinants of SI within a 700-kb hemizygous region present only in one compatibility group. We then demonstrate that the homologous hemizygous region also controls distyly in jasmine. Phylogenetic analyses support a common origin of both systems, following a segmental genomic duplication in a common ancestor. Examination of the gene content of the hemizygous region in different jasmine and olive species suggests that the mechanisms determining compatibility groups and floral phenotypes (whether homomorphic or dimorphic) in Oleaceae rely on the presence/absence of two genes involved in gibberellin and brassinosteroid regulation.

The homomorphic self-incompatibility system in Oleaceae is controlled by a hemizygous genomic region expressing a gibberellin pathway gene

In flowering plants, outcrossing is commonly ensured by self-incompatibility (SI) systems. These can be homomorphic (typically with many different allelic specificities) or can accompany flower heteromorphism (mostly with just two specificities and corresponding floral types). The SI system of the Oleaceae family is unusual, with the long-term maintenance of only two specificities but often without flower morphology differences. To elucidate the genomic architecture and molecular basis of this SI system, we obtained chromosome-scale genome assemblies of Phillyrea angustifolia individuals and related them to a genetic map. The S-locus region proved to have a segregating 543-kb indel unique to one specificity, suggesting a hemizygous region, as observed in all distylous systems so far studied at the genomic level. Only one of the predicted genes in this indel region is found in the olive tree, Olea europaea, genome, also within a segregating indel. We describe complete association between the presence/absence of this gene and the SI types determined for individuals of seven distantly related Oleaceae species. This gene is predicted to be involved in catabolism of the gibberellic acid (GA) hormone, and experimental manipulation of GA levels in developing buds modified the male and female SI responses of the two specificities in different ways. Our results provide a unique example of a homomorphic SI system, where a single conserved gibberellin-related gene in a hemizygous indel underlies the long-term maintenance of two groups of reproductive compatibility.

Twelve newly assembled jasmine chloroplast genomes: unveiling genomic diversity, phylogenetic relationships and evolutionary patterns among Oleaceae and Jasminum species

BACKGROUND

Jasmine (Jasminum), renowned for its ornamental value and captivating fragrance, has given rise to numerous species and accessions. However, limited knowledge exists regarding the evolutionary relationships among various Jasminum species.

RESULTS

In the present study, we sequenced seven distinct Jasminum species, resulting in the assembly of twelve high-quality complete chloroplast (cp) genomes. Our findings revealed that the size of the 12 cp genomes ranged from 159 to 165 kb and encoded 134-135 genes, including 86-88 protein-coding genes, 38-40 tRNA genes, and 8 rRNA genes. J. nudiflorum exhibited a larger genome size compared to other species, mainly attributed to the elevated number of forward repeats (FRs). Despite the typically conservative nature of chloroplasts, variations in the presence or absence of accD have been observed within J. sambac. The calculation of nucleotide diversity (Pi) values for 19 cp genomes indicated that potential mutation hotspots were more likely to be located in LSC regions than in other regions, particularly in genes ycf2, rbcL, atpE, ndhK, and ndhC (Pi > 0.2). Ka/Ks values revealed strong selection pressure on the genes rps2, atpA, rpoA, rpoC1, and rpl33 when comparing J. sambac with the three most closely related species (J. auriculatum, J. multiflorum, and J. dichotomum). Additionally, SNP identification, along with the results of Structure, PCA, and phylogenetic tree analyses, divided the Jasminum cp genomes into six groups. Notably, J. polyanthum showed gene flow signals from both the G5 group (J. nudiflorum) and the G3 group (J. tortuosum and J. fluminense). Phylogenetic tree analysis reflected that most species from the same genus clustered together with robust support in Oleaceae, strongly supporting the monophyletic nature of cp genomes within the genus Jasminum.

CONCLUSION

Overall, this study provides comprehensive insights into the genomic composition, variation, and phylogenetic relationships among various Jasminum species. These findings enhance our understanding of the genetic diversity and evolutionary history of Jasminum.

LocoGSE, a sequence-based genome size estimator for plants

Extensive research has focused on exploring the range of genome sizes in eukaryotes, with a particular emphasis on land plants, where significant variability has been observed. Accurate estimation of genome size is essential for various research purposes, but existing sequence-based methods have limitations, particularly for low-coverage datasets. In this study, we introduce LocoGSE, a novel genome size estimator designed specifically for low-coverage datasets generated by genome skimming approaches. LocoGSE relies on mapping the reads on single copy consensus proteins without the need for a reference genome assembly. We calibrated LocoGSE using 430 low-coverage Angiosperm genome skimming datasets and compared its performance against other estimators. Our results demonstrate that LocoGSE accurately predicts monoploid genome size even at very low depth of coverage (<1X) and on highly heterozygous samples. Additionally, LocoGSE provides stable estimates across individuals with varying ploidy levels. LocoGSE fills a gap in sequence-based plant genome size estimation by offering a user-friendly and reliable tool that does not rely on high coverage or reference assemblies. We anticipate that LocoGSE will facilitate plant genome size analysis and contribute to evolutionary and ecological studies in the field. Furthermore, at the cost of an initial calibration, LocoGSE can be used in other lineages.

Nuclear phylogenomics of angiosperms and insights into their relationships and evolution

Angiosperms (flowering plants) are by far the most diverse land plant group with over 300,000 species. The sudden appearance of diverse angiosperms in the fossil record was referred to by Darwin as the "abominable mystery," hence contributing to the heightened interest in angiosperm evolution. Angiosperms display wide ranges of morphological, physiological, and ecological characters, some of which have probably influenced their species richness. The evolutionary analyses of these characteristics help to address questions of angiosperm diversification and require well resolved phylogeny. Following the great successes of phylogenetic analyses using plastid sequences, dozens to thousands of nuclear genes from next-generation sequencing have been used in angiosperm phylogenomic analyses, providing well resolved phylogenies and new insights into the evolution of angiosperms. In this review we focus on recent nuclear phylogenomic analyses of large angiosperm clades, orders, families, and subdivisions of some families and provide a summarized Nuclear Phylogenetic Tree of Angiosperm Families. The newly established nuclear phylogenetic relationships are highlighted and compared with previous phylogenetic results. The sequenced genomes of Amborella, Nymphaea, Chloranthus, Ceratophyllum, and species of monocots, Magnoliids, and basal eudicots, have facilitated the phylogenomics of relationships among five major angiosperms clades. All but one of the 64 angiosperm orders were included in nuclear phylogenomics with well resolved relationships except the placements of several orders. Most families have been included with robust and highly supported placements, especially for relationships within several large and important orders and families. Additionally, we examine the divergence time estimation and biogeographic analyses of angiosperm on the basis of the nuclear phylogenomic frameworks and discuss the differences compared with previous analyses. Furthermore, we discuss the implications of nuclear phylogenomic analyses on ancestral reconstruction of morphological, physiological, and ecological characters of angiosperm groups, limitations of current nuclear phylogenomic studies, and the taxa that require future attention.

Reticulate evolution of the tertiary relict Osmanthus

When interspecific gene flow is common, species relationships are more accurately represented by a phylogenetic network than by a bifurcating tree. This study aimed to uncover the role of introgression in the evolution of Osmanthus, the only genus of the subtribe Oleinae (Oleaceae) with its distribution center in East Asia. We built species trees, detected introgression, and constructed networks using multiple kinds of sequencing data (whole genome resequencing, transcriptome sequencing, and Sanger sequencing of nrDNA) combined with concatenation and coalescence approaches. Then, based on well-understood species relationships, historical biogeographic analyses and diversification rate estimates were employed to reveal the history of Osmanthus. Osmanthus originated in mid-Miocene Europe and dispersed to the eastern Tibetan Plateau in the late Miocene. Thereafter, it continued to spread eastwards. Phylogenetic conflict is common within the 'Core Osmanthus' clade and is seen at both early and late stages of diversification, leading to hypotheses of net-like species relationships. Incomplete lineage sorting proved ineffective in explaining phylogenetic conflicts and thus supported introgression as the main cause of conflicts. This study elucidates the diversification history of a relict genus in the subtropical regions of eastern Asia and reveals that introgression had profound effects on its evolutionary history.

One hundred single-copy nuclear sequence markers for olive variety identification: a case of fingerprinting database construction in China

Olive is an ancient oil-producing tree, widely cultivated in Mediterranean countries, and now spread to other areas of the world, including China. Recently, several molecular databases were constructed in different countries and platforms for olive identification using simple sequence repeats (SSRs) or single-nucleotide polymorphisms (SNPs). However, comparing their results across laboratories was difficult. Herein, hundreds of polymorphic single-copy nuclear sequence markers were developed from the olive genome. Using the advantage of multiplex PCR amplification and high-throughput sequencing, a fingerprint database was constructed for the majority of olives cultivated in China. We used 100 high-quality sequence loci and estimated the genetic diversity and structure among all these varieties. We found that compared with that based on SSRs, the constructed fingerprint database based on these 100 sequences or a few of them, could provide a reliable olive variety identification platform in China, with high discrimination among different varieties using the principle of BLAST algorithm. An example of such identification platform based on this study was displayed on the web for the olive database in China (http://olivedb.cn/jianding). After resolving redundant genotypes, we identified 126 olive varieties with distinct genotypes in China. These varieties could be divided into two clusters, and it was revealed that the grouping of the varieties has a certain relationship with their origin. Herein, it is concluded that these single-copy orthologous nuclear sequences could be used to construct a universal fingerprint database of olives across different laboratories and platforms inexpensively. Based on such a database, variety identification can be performed easily by any laboratory, which would further facilitate olive breeding and variety exchange globally.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01434-9.

Genealogical tracing of Olea europaea species and pedigree relationships of var. europaea using chloroplast and nuclear markers

BACKGROUND

Olive is one of the most cultivated species in the Mediterranean Basin and beyond. Despite being extensively studied for its commercial relevance, the origin of cultivated olive and the history of its domestication remain open questions. Here, we present a genealogical and kinship relationships analysis by mean of chloroplast and nuclear markers of different genera, subgenus, species, subspecies, ecotypes, cultivated, ancient and wild types, which constitutes one of the most inclusive research to date on the diversity within Olea europaea species. A complete survey of the variability across the nuclear and plastid genomes of different genotypes was studied through single nucleotide polymorphisms, indels (insertions and deletions), and length variation.

RESULTS

Fifty-six different chlorotypes were identified among the Oleaceae family including Olea europaea, other species and genera. The chloroplast genome evolution, within Olea europaea subspecies, probably started from subsp. cuspidata, which likely represents the ancestor of all the other subspecies and therefore of wild types and cultivars. Our study allows us to hypothesize that, inside the subspecies europaea containing cultivars and the wild types, the ancestral selection from var. sylvestris occurred both in the eastern side of the Mediterranean and in the central-western part of Basin. Moreover, it was elucidated the origin of several cultivars, which depends on the introduction of eastern cultivars, belonging to the lineage E1, followed by crossing and replacement of the autochthonous olive germplasm of central-western Mediterranean Basin. In fact, our study highlighted that two main 'founders' gave the origin to more than 60% of analyzed olive cultivars. Other secondary founders, which strongly contributed to give origin to the actual olive cultivar diversity, were already detected.

CONCLUSIONS

The application of comparative genomics not only paves the way for a better understanding of the phylogenetic relationships within the Olea europaea species but also provides original insights into other elusive evolutionary processes, such as chloroplast inheritance and parentage inside olive cultivars, opening new scenarios for further research such as the association studies and breeding programs.

Museomics and phylogenomics of lovebirds (Psittaciformes, Psittaculidae, Agapornis) using low-coverage whole-genome sequencing

Natural history collections contain specimens that provide important insights into studies of ecology and evolution. With the advancement of high-throughput sequencing, historical DNA (hDNA) from museum specimens has become a valuable source of genomic data to study the evolutionary history of organisms. Low-coverage whole genome sequencing (WGS) has been increasingly applied to museum specimens for analyzing organelle genomes, but is still uncommon for genotyping the nuclear DNA fraction. In this study, we applied low-coverage WGS to phylogenomic analyses of parrots in the genus Agapornis by including both modern samples and historical specimens of ∼100-year-old. Agapornis are small-sized African and Malagasy parrots with diverse characters. Earlier phylogenetic studies failed to resolve the positions of some key lineages, prohibiting a robust interpretation of the biogeography and evolution of these African parrots. Here, we demonstrated the use of low-coverage WGS for generating both mitochondrial and nuclear genomic data, and evaluated data quality differences between modern and historical samples. Our resolved Agapornis phylogeny indicates the ancestor of Agapornis likely colonized Madagascar from Australasia by trans-oceanic dispersal events before dispersing to the African continent. Genome-wide SNPs also allowed us to identify the parental origins of hybrid Agapornis individuals. This study demonstrates the potential of applying low-coverage WGS to phylogenomics and population genomics analyses and illustrates how including historical museum specimens can address outstanding questions regarding the evolutionary history of contemporary lineages.

Leaf anatomy explains the strength of C4 activity within the grass species Alloteropsis semialata

C₄ photosynthesis results from anatomical and biochemical characteristics that together concentrate CO₂ around ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco), increasing productivity in warm conditions. This complex trait evolved through the gradual accumulation of components, and particular species possess only some of these, resulting in weak C₄ activity. The consequences of adding C₄ components have been modelled and investigated through comparative approaches, but the intraspecific dynamics responsible for strengthening the C₄ pathway remain largely unexplored. Here, we evaluate the link between anatomical variation and C₄ activity, focusing on populations of the photosynthetically diverse grass Alloteropsis semialata that fix various proportions of carbon via the C₄ cycle. The carbon isotope ratios in these populations range from values typical of C₃ to those typical of C₄ plants. This variation is statistically explained by a combination of leaf anatomical traits linked to the preponderance of bundle sheath tissue. We hypothesize that increased investment in bundle sheath boosts the strength of the intercellular C₄ pump and shifts the balance of carbon acquisition towards the C₄ cycle. Carbon isotope ratios indicating a stronger C₄ pathway are associated with warmer, drier environments, suggesting that incremental anatomical alterations can lead to the emergence of C₄ physiology during local adaptation within metapopulations.

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.

Phylogeny and Historical Biogeography of the East Asian Clematis Group, Sect. Tubulosae, Inferred from Phylogenomic Data

The evolutionary history of Clematis section Tubulosae, an East Asian endemic lineage, has not been comprehensively studied. In this study, we reconstruct the phylogeny of this section with a complete sampling using a phylogenomic approach. The genome skimming method was applied to obtain the complete plastome sequence, the nuclear ribosomal DNA (nrDNA), and the nuclear SNPs data for phylogenetic reconstruction. Using a Bayesian molecular clock approach and ancestral range reconstruction, we reconstruct biogeographical history and discuss the biotic and abiotic factors that may have shaped the distribution patterns of the section. Both nuclear datasets better resolved the phylogeny of the sect. Tubulosae than the plastome sequence. Sect. Tubulosae was resolved as a monophyletic group sister to a clade mainly containing species from the sect. Clematis and sect. Aspidanthera. Within sect. Tubulosae, two major clades were resolved by both nuclear datasets. Two continental taxa, C. heracleifolia and C. tubulosa var. ichangensis, formed one clade. One continental taxon, C. tubulosa, and all the other species from Taiwan island, the Korean peninsula, and the Japanese archipelago formed the other clade. Molecular dating results showed that sect. Tubulosae diverged from its sister clade in the Pliocene, and all the current species diversified during the Pleistocene. Our biogeographical reconstruction suggested that sect. Tubulosae evolved and began species diversification, most likely in mainland China, then dispersed to the Korean peninsula, and then expanded its range through the Japanese archipelago to Taiwan island. Island species diversity may arise through allopatric speciation by vicariance events following the range fragmentation triggered by the climatic oscillation and sea level change during the Pleistocene epoch. Our results highlight the importance of climatic oscillation during the Pleistocene to the spatial-temporal diversification patterns of the sect. Tubulosae.

The complete chloroplast genome sequence of Syringa oblata Lindl. var. alba Hort. ex Rehd. 1763 (Oleaceae)

Syringa oblata var. alba is a shrub or a small tree from China with high ornamental, medicinal, and edible value. Here, we present its first complete chloroplast genome. The entire circular genome is 155,648 bp in length, with large single-copy (LSC) length of 86,247, small single-copy (SSC) length of 17,937, inverted repeat (IR) length of 25,732, and GC content of 37.9%. One hundred and thirty-two genes, including 88 protein-coding, 36 tRNA, and eight rRNA genes were predicted. A phylogenetic tree of 25 plant species was constructed based on the maximum-likelihood method, indicating that S. oblata var. alba, S. vulgaris, and S. oblata form a sister group. This study will provide valuable basic information for phylogeny, species identification, and varieties breeding of this species.

A pipeline for assembling low copy nuclear markers from plant genome skimming data for phylogenetic use

Background

Genome skimming is a popular method in plant phylogenomics that do not include a biased enrichment step, relying on random shallow sequencing of total genomic DNA. From these data the plastome is usually readily assembled and constitutes the bulk of phylogenetic information generated in these studies. Despite a few attempts to use genome skims to recover low copy nuclear loci for direct phylogenetic use, such endeavor remains neglected. Causes might include the trade-off between libraries with few reads and species with large genomes (i.e., missing data caused by low coverage), but also might relate to the lack of pipelines for data assembling.

Methods

A pipeline and its companion R package designed to automate the recovery of low copy nuclear markers from genome skimming libraries are presented. Additionally, a series of analyses aiming to evaluate the impact of key assembling parameters, reference selection and missing data are presented.

Results

A substantial amount of putative low copy nuclear loci was assembled and proved useful to base phylogenetic inference across the libraries tested (4 to 11 times more data than previously assembled plastomes from the same libraries).

Discussion

Critical aspects of assembling low copy nuclear markers from genome skims include the minimum coverage and depth of a sequence to be used. More stringent values of these parameters reduces the amount of assembled data and increases the relative amount of missing data, which can compromise phylogenetic inference, in turn relaxing the same parameters might increase sequence error. These issues are discussed in the text, and parameter tuning through multiple comparisons tracking their effects on support and congruence is highly recommended when using this pipeline. The skimmingLoci pipeline (https://github.com/mreginato/skimmingLoci) might stimulate the use of genome skims to recover nuclear loci for direct phylogenetic use, increasing the power of genome skimming data to resolve phylogenetic relationships, while reducing the amount of sequenced DNA that is commonly wasted.

Genome-partitioning strategy, plastid and nuclear phylogenomic discordance, and its evolutionary implications of Clematis (Ranunculaceae)

Clematis is one of the largest genera of Ranunculaceae with many phylogenetic problems left to be resolved. Clematis species have considerable genome size of more than 7 Gbp, and there was no whole-genome reference sequence published in this genus. This raises difficulties in acquiring nuclear genome data for its phylogenetic analysis. Previous studies based on Sanger sequencing data, plastid genome data, and nrDNA sequences did not well resolve the phylogeny of Clematis. In this study, we used genome skimming and transcriptome data to assemble the plastid genome sequences, nuclear single nucleotide polymorphisms (SNPs) datasets, and single-copy nuclear orthologous genes (SCOGs) to reconstruct the phylogenetic backbone of Clematis, and test effectiveness of these genome partitioning methods. We also further analyzed the discordance among nuclear gene trees and between plastid and nuclear phylogenies. The results showed that the SCOGs datasets, assembled from transcriptome method, well resolved the phylogenetic backbone of Clematis. The nuclear SNPs datasets from genome skimming method can also produce similar results with the SCOGs data. In contrast to the plastid phylogeny, the phylogeny resolved by nuclear genome data is more robust and better corresponds to morphological characters. Our results suggested that rapid species radiation may have generated high level of incomplete lineage sorting, which was the major cause of nuclear gene discordance. Our simulation also showed that there may have been frequent interspecific hybridization events, which led to some of the cyto-nuclear discordances in Clematis. This study not only provides the first robust phylogenetic backbone of Clematis based on nuclear genome data, but also provides suggestions of genome partitioning strategies for the phylogenomic study of other plant taxa.

Hybridisation and chloroplast capture between distinct Themeda triandra lineages in Australia

Ecotypes are distinct populations within a species that are adapted to specific environmental conditions. Understanding how these ecotypes become established, and how they interact when reunited, is fundamental to elucidating how ecological adaptations are maintained. This study focuses on Themeda triandra, a dominant grassland species across Asia, Africa and Australia. It is the most widespread plant in Australia, where it has distinct ecotypes that are usually restricted to either wetter and cooler coastal regions or the drier and hotter interior. We generate a reference genome for T. triandra and use whole genome sequencing for over 80 Themeda accessions to reconstruct the evolutionary history of T. triandra and related taxa. Organelle phylogenies confirm that Australia was colonized by T. triandra twice, with the division between ecotypes predating their arrival in Australia. The nuclear genome provides evidence of differences in the dominant ploidal level and gene-flow among the ecotypes. In northern Queensland there appears to be a hybrid zone between ecotypes with admixed nuclear genomes and shared chloroplast haplotypes. Conversely, in the cracking claypans of Western Australia, there is cytonuclear discordance with individuals possessing the coastal chloroplast and interior clade nuclear genome. This chloroplast capture is potentially a result of adaptive introgression, with selection detected in the rpoC2 gene which is associated with water use efficiency. The reason that T. triandra is the most widespread plant in Australia appears to be a result of distinct ecotypic genetic variation and genome duplication, with the importance of each depending on the geographic scale considered.

Simultaneously collecting coding and non-coding phylogenomic data using homemade full-length cDNA probes, tested by resolving the high-level relationships of Colubridae

Resolving intractable phylogenetic relationships often requires simultaneously analyzing a large number of coding and non-coding orthologous loci. To gather both coding and non-coding data, traditional sequence capture methods require custom-designed commercial probes. Here, we present a cost-effective sequence capture method based on homemade probes, to capture thousands of coding and non-coding orthologous loci simultaneously, suitable for all organisms. This approach, called “FLc-Capture,” synthesizes biotinylated full-length cDNAs from mRNA as capture probes, eliminates the need for costly commercial probe design and synthesis. To demonstrate the utility of FLc-Capture, we prepared full-length cDNA probes from mRNA extracted from a common colubrid snake. We performed capture experiments with these homemade cDNA probes and successfully obtained thousands of coding and non-coding genomic loci from 24 Colubridae species and 12 distantly related snake species of other families. The average capture specificity of FLc-Capture across all tested snake species is 35%, similar to the previously published EecSeq method. We constructed two phylogenomic data sets, one including 1,075 coding loci (∼817,000 bp) and the other including 1,948 non-coding loci (∼1,114,000 bp), to study the phylogeny of Colubridae. Both data sets yielded highly similar and well-resolved trees, with 85% of nodes having &gt;95% bootstrap support. Our experimental tests show that FLc-Capture is a flexible, fast, and cost-effective sequence capture approach for simultaneously gathering coding and non-coding phylogenomic data sets to study intractable phylogenetic questions. We hope that this method will serve as a new data collection tool for evolutionary biologists working in the era of phylogenomics.

Phylogenomic approaches untangle early divergences and complex diversifications of the olive plant family

BACKGROUND

Deep-branching phylogenetic relationships are often difficult to resolve because phylogenetic signals are obscured by the long history and complexity of evolutionary processes, such as ancient introgression/hybridization, polyploidization, and incomplete lineage sorting (ILS). Phylogenomics has been effective in providing information for resolving both deep- and shallow-scale relationships across all branches of the tree of life. The olive family (Oleaceae) is composed of 25 genera classified into five tribes with tribe Oleeae consisting of four subtribes. Previous phylogenetic analyses showed that ILS and/or hybridization led to phylogenetic incongruence in the family. It was essential to distinguish phylogenetic signal conflicts, and explore mechanisms for the uncertainties concerning relationships of the olive family, especially at the deep-branching nodes.

RESULTS

We used the whole plastid genome and nuclear single nucleotide polymorphism (SNP) data to infer the phylogenetic relationships and to assess the variation and rates among the main clades of the olive family. We also used 2608 and 1865 orthologous nuclear genes to infer the deep-branching relationships among tribes of Oleaceae and subtribes of tribe Oleeae, respectively. Concatenated and coalescence trees based on the plastid genome, nuclear SNPs and multiple nuclear genes suggest events of ILS and/or ancient introgression during the diversification of Oleaceae. Additionally, there was extreme heterogeneity in the substitution rates across the tribes. Furthermore, our results supported that introgression/hybridization, rather than ILS, is the main factor for phylogenetic discordance among the five tribes of Oleaceae. The tribe Oleeae is supported to have originated via ancient hybridization and polyploidy, and its most likely parentages are the ancestral lineage of Jasmineae or its sister group, which is a "ghost lineage," and Forsythieae. However, ILS and ancient introgression are mainly responsible for the phylogenetic discordance among the four subtribes of tribe Oleeae.

CONCLUSIONS

This study showcases that using multiple sequence datasets (plastid genomes, nuclear SNPs and thousands of nuclear genes) and diverse phylogenomic methods such as data partition, heterogeneous models, quantifying introgression via branch lengths (QuIBL) analysis, and species network analysis can facilitate untangling long and complex evolutionary processes of ancient introgression, paleopolyploidization, and ILS.

Plastid genomes reveal evolutionary shifts in elevational range and flowering time of Osmanthus (Oleaceae)

Species of Osmanthus are economically important ornamental trees, yet information regarding their plastid genomes (plastomes) have rarely been reported, thus hindering taxonomic and evolutionary studies of this small but enigmatic genus. Here, we performed comparative genomics and evolutionary analyses on plastomes of 16 of the 28 currently accepted species, with 11 plastomes newly sequenced. Phylogenetic studies identified four main lineages within the genus that are here designated the: “Caucasian Osmanthus” (corresponding to O. decorus), “Siphosmanthus” (corresponding to O. sect. Siphosmanthus), “O. serrulatus + O. yunnanensis,” and “Core Osmanthus: (corresponding to O. sect. Osmanthus + O. sect. Linocieroides). Molecular clock analysis suggested that Osmanthus split from its sister clade c. 15.83 Ma. The estimated crown ages of the lineages were the following: genus Osmanthus at 12.66 Ma; “Siphosmanthus” clade at 5.85 Ma; “O. serrulatus + O. yunnanensis” at 4.89 Ma; and “Core Osmanthus: clade at 6.2 Ma. Ancestral state reconstructions and trait mapping showed that ancestors of Osmanthus were spring flowering and originated at lower elevations. Phylogenetic principal component analysis clearly distinguished spring‐flowering species from autumn‐flowering species, suggesting that flowering time differentiation is related to the difference in ecological niches. Nucleotide substitution rates of 80 common genes showed slow evolutionary pace and low nucleotide variations, all genes being subjected to purifying selection.

The phylogeny of Salix revealed by whole genome re-sequencing suggests different sex-determination systems in major groups of the genus

BACKGROUND AND AIMS

The largest genus of Salicaceae sensu lato, Salix, has been shown to consist of two main clades: clade Salix, in which species have XY sex-determination systems (SDSs) on chromosome 7, and clade Vetrix including species with ZW SDSs on chromosome 15. Here, we test the utility of whole genome re-sequencing (WGR) for phylogenomic reconstructions of willows to infer changes between different SDSs.

METHODS

We used more than 1 TB of WGR data from 70 Salix taxa to ascertain single nucleotide polymorphisms on the autosomes, the sex-linked regions (SLRs) and the chloroplast genomes, for phylogenetic and species tree analyses. To avoid bias, we chose reference genomes from both groups, Salix dunnii from clade Salix and S. purpurea from clade Vetrix.

KEY RESULTS

Two main largely congruent groups were recovered: the paraphyletic Salix grade and the Vetrix clade. The autosome dataset trees resolved four subclades (C1-C4) in Vetrix. C1 and C2 comprise species from the Hengduan Mountains and adjacent areas and from Eurasia, respectively. Section Longifoliae (C3) grouped within the Vetrix clade but fell into the Salix clade in trees based on the chloroplast dataset analysis. Salix triandra from Eurasia (C4) was revealed as sister to the remaining species of clade Vetrix. In Salix, the polyploid group C5 is paraphyletic to clade Vetrix and subclade C6 is consistent with Argus's subgenus Protitea. Chloroplast datasets separated both Vetrix and Salix as monophyletic, and yielded C5 embedded within Salix. Using only diploid species, both the SLR and autosomal datasets yielded trees with Vetrix and Salix as well-supported clades.

CONCLUSION

WGR data are useful for phylogenomic analyses of willows. The different SDSs may contribute to the isolation of the two major groups, but the reproductive barrier between them needs to be studied.

Hybridization boosts dispersal of two contrasted ecotypes in a grass species

Genetic exchanges between closely related groups of organisms with different adaptations have well-documented beneficial and detrimental consequences. In plants, pollen-mediated exchanges affect the sorting of alleles across physical landscapes and influence rates of hybridization. How these dynamics affect the emergence and spread of novel phenotypes remains only partially understood. Here, we use phylogenomics and population genomics to retrace the origin and spread of two geographically overlapping ecotypes of the African grass Alloteropsis angusta. In addition to an ecotype inhabiting wetlands, we report the existence of a previously undescribed ecotype inhabiting Miombo woodlands and grasslands. The two ecotypes are consistently associated with different nuclear groups, which represent an advanced stage of divergence with secondary low-level gene flow. However, the seed-transported chloroplast genomes are consistently shared by distinct ecotypes inhabiting the same region. These patterns suggest that the nuclear genome of one ecotype can enter the seeds of the other via occasional pollen movements with sorting of nuclear groups in subsequent generations. The contrasting ecotypes of A. angusta can thus use each other as a gateway to new locations across a large part of Africa, showing that hybridization can facilitate the geographical dispersal of distinct ecotypes of the same grass species.

A bioinformatic platform to integrate target capture and whole genome sequences of various read depths for phylogenomics

The increasing availability of short‐read whole genome sequencing (WGS) provides unprecedented opportunities to study ecological and evolutionary processes. Although loci of interest can be extracted from WGS data and combined with target sequence data, this requires suitable bioinformatic workflows. Here, we test different assembly and locus extraction strategies and implement them into secapr, a pipeline that processes short‐read data into multilocus alignments for phylogenetics and molecular ecology analyses. We integrate the processing of data from low‐coverage WGS (<30×) and target sequence capture into a flexible framework, while optimizing de novo contig assembly and loci extraction. Specifically, we test different assembly strategies by contrasting their ability to recover loci from targeted butterfly protein‐coding genes, using four data sets: a WGS data set across different average coverages (10×, 5× and 2×) and a data set for which these loci were enriched prior to sequencing via target sequence capture. Using the resulting de novo contigs, we account for potential errors within contigs and infer phylogenetic trees to evaluate the ability of each assembly strategy to recover species relationships. We demonstrate that choosing multiple sizes of kmer simultaneously for assembly results in the highest yield of extracted loci from de novo assembled contigs, while data sets derived from sequencing read depths as low as 5× recovers the expected species relationships in phylogenetic trees. By making the tested assembly approaches available in the secapr pipeline, we hope to inspire future studies to incorporate complementary data and make an informed choice on the optimal assembly strategy.

Widespread coexistence of self-compatible and self-incompatible phenotypes in a diallelic self-incompatibility system in Ligustrum vulgare (Oleaceae)

The breakdown of self-incompatibility (SI) in angiosperms is one of the most commonly observed evolutionary transitions. While multiple examples of SI breakdown have been documented in natural populations, there is strikingly little evidence of stable within-population polymorphism with both inbreeding (self-compatible) and outcrossing (self-incompatible) individuals. This absence of breeding system polymorphism corroborates theoretical expectations that predict that in/outbreeding polymorphism is possible only under very restricted conditions. However, theory also predicts that a diallelic sporophytic SI system should facilitate the maintenance of such polymorphism. We tested this prediction by studying the breeding system of Ligustrum vulgare L., an insect-pollinated hermaphroditic species of the Oleaceae family. Using stigma tests with controlled pollination and paternity assignment of open-pollinated progenies, we confirmed the existence of two self-incompatibility groups in this species. We also demonstrated the occurrence of self-compatible individuals in different populations of Western Europe arising from a mutation affecting the functioning of the pollen component of SI. Our results show that the observed low frequency of self-compatible individuals in natural populations is compatible with theoretical predictions only if inbreeding depression is very high.

How to Choose a Good Marker to Analyze the Olive Germplasm (Olea europaea L.) and Derived Products

The olive tree (Olea europaea L.) is one of the most cultivated crops in the Mediterranean basin. Its economic importance is mainly due to the intense production of table olives and oil. Cultivated varieties are characterized by high morphological and genetic variability and present a large number of synonyms and homonyms. This necessitates the introduction of a rapid and accurate system for varietal identification. In the past, the recognition of olive cultivars was based solely on analysis of the morphological traits, however, these are highly influenced by environmental conditions. Therefore, over the years, several methods based on DNA analysis were developed, allowing a more accurate and reliable varietal identification. This review aims to investigate the evolving history of olive tree characterization approaches, starting from the earlier morphological methods to the latest technologies based on molecular markers, focusing on the main applications of each approach. Furthermore, we discuss the impact of the advent of next generation sequencing and the recent sequencing of the olive genome on the strategies used for the development of new molecular markers.

Genome skims analysis of betel palms (Areca spp., Arecaceae) and development of a profiling method to assess their plastome diversity

The betel nut (Areca catechu L., Arecaceae) is a monoecious cultivated palm tree that is widespread in tropical regions. It is mainly cultivated for producing areca nuts, from which seeds are extracted and chewed by local populations principally in the Indo-Pacific region. Seeds contain alkaloids which are central nervous system stimulants and are highly addictive. Wild relatives of the betel nut are distributed in South Asia and Australasia, with ca. 40-50 Areca species currently recognized. The geographic origin(s) of the betel nut and its subsequent diffusion and diversification remains poorly documented. Here, a genome skimming approach was applied to screen nucleotidic variation in the most abundant genomic regions. Low coverage sequencing data allowed us to assemble full plastomes, mitochondrial regions (either full mitogenomes or the full set of mitochondrial genes) and the nuclear ribosomal DNA cluster for nine representatives of the Areca genus collected in the field and herbarium collections (including a 182-years old specimen collected during the Dumont d'Urville's expedition). These three genomic compartments provided similar phylogenetic signals, and revealed very low genomic diversity in our sample of cultivated betel nut. We finally developed a genotyping method targeting 34 plastid DNA microsatellites. This plastome profiling approach is useful for tracing the spread of matrilineages, and in combination with nuclear genomic data, can resolve the history of the betel nut. Our method also proves to be efficient for analyzing herbarium specimens, even those collected >100 years ago.

Characteristic and phylogenetic analyses of chloroplast genome for Syringa komarowii C.K.Schneid. (Oleaceae) from Huoditang, China, an important horticultural plant

The genus Syringa (Oleaceae) comprises some of the most important cultivated horticultural trees worldwide. Syringa komarowii, is one of these important horticultural plants. We sequenced the complete chloroplast (cp) genomes of S. komarowii (previously identified as var. reflexa) using Illumina Hiseq X Ten platform. The cp genome exhibits a typical quadripartite structure with 158,020 bp in length, which consists of two copies of inverted repeat (IR) regions (25,676 bp) separated by a large single copy (LSC) region (87,628 bp) and a small single copy (SSC) region (19,040 bp). The cp genome of S. komarowii encodes a total of 132 genes, including 88 protein-coding genes, 36 tRNA genes and 8 rRNA genes. A maximum likelihood (ML) phylogenetic analysis resolved that S. komarowii in a clade with S. wolfii and S. yunnanensis. The ML tree also showed Syringa appeared more closely related to Ligustrum than to the other genera in the Oleaceae.

Resolving the Phylogeny of the Olive Family (Oleaceae): Confronting Information from Organellar and Nuclear Genomes

The olive family, Oleaceae, is a group of woody plants comprising 28 genera and ca. 700 species, distributed on all continents (except Antarctica) in both temperate and tropical environments. It includes several genera of major economic and ecological importance such as olives, ash trees, jasmines, forsythias, osmanthuses, privets and lilacs. The natural history of the group is not completely understood yet, but its diversification seems to be associated with polyploidisation events and the evolution of various reproductive and dispersal strategies. In addition, some taxonomical issues still need to be resolved, particularly in the paleopolyploid tribe Oleeae. Reconstructing a robust phylogenetic hypothesis is thus an important step toward a better comprehension of Oleaceae's diversity. Here, we reconstructed phylogenies of the olive family using 80 plastid coding sequences, 37 mitochondrial genes, the complete nuclear ribosomal cluster and a small multigene family encoding phytochromes (phyB and phyE) of 61 representative species. Tribes and subtribes were strongly supported by all phylogenetic reconstructions, while a few Oleeae genera are still polyphyletic (Chionanthus, Olea, Osmanthus, Nestegis) or paraphyletic (Schrebera, Syringa). Some phylogenetic relationships among tribes remain poorly resolved with conflicts between topologies reconstructed from different genomic regions. The use of nuclear data remains an important challenge especially in a group with ploidy changes (both paleo- and neo-polyploids). This work provides new genomic datasets that will assist the study of the biogeography and taxonomy of the whole Oleaceae.

Continued Adaptation of C4 Photosynthesis After an Initial Burst of Changes in the Andropogoneae Grasses

C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} photosynthesis is a complex trait that sustains fast growth and high productivity in tropical and subtropical conditions and evolved repeatedly in flowering plants. One of the major C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} lineages is Andropogoneae, a group of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\sim $\end{document}1200 grass species that includes some of the world’s most important crops and species dominating tropical and some temperate grasslands. Previous efforts to understand C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} evolution in the group have compared a few model C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} plants to distantly related C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{3}$\end{document} species so that changes directly responsible for the transition to C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} could not be distinguished from those that preceded or followed it. In this study, we analyze the genomes of 66 grass species, capturing the earliest diversification within Andropogoneae as well as their C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{3}$\end{document} relatives. Phylogenomics combined with molecular dating and analyses of protein evolution show that many changes linked to the evolution of C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} photosynthesis in Andropogoneae happened in the Early Miocene, between 21 and 18 Ma, after the split from its C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{3}$\end{document} sister lineage, and before the diversification of the group. This initial burst of changes was followed by an extended period of modifications to leaf anatomy and biochemistry during the diversification of Andropogoneae, so that a single C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} origin gave birth to a diversity of C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} phenotypes during 18 million years of speciation events and migration across geographic and ecological spaces. Our comprehensive approach and broad sampling of the diversity in the group reveals that one key transition can lead to a plethora of phenotypes following sustained adaptation of the ancestral state. [Adaptive evolution; complex traits; herbarium genomics; Jansenelleae; leaf anatomy; Poaceae; phylogenomics.]

Deciphering the organelle genomes and transcriptomes of a common ornamental plant Ligustrum quihoui reveals multiple fragments of transposable elements in the mitogenome

Ligustrum quihoui (L. quihoui) is an important hedge material for landscaping and also possesses medicinal value. To generate genomic resources for better understanding the evolutionary history of this important plant, the organelle genomes of L. quihoui are de novo assembled and functionally annotated. Compared with other Oleaceae species, the 163,069 bp chloroplast genome of L. quihoui exhibits a typical quadripartite structure with highly conserved gene content and gene order, while the 848,451 bp mitochondrial genome of L. quihoui exhibits highly divergent genome size and gene content. Codon usage analyses show that genes related with photosynthesis and mitochondrial respiratory chain show inconsistent codon biases. A total of 48,760 bp transposable elements (TEs) fragments and 41,887 bp chloroplast-like sequences are found in the L. quihoui mitochondrial genome. A striking discrepancy of RNA editing between the two organelle genomes is found in L. quihoui, in which 146 mitochondrial editing sites coexist with only 43 such sites in chloroplast. Based on DNA and RNA-Seq data, we propose that GTG may act as the start codon of mitochondrial rpl16 in Oleaceae species. Phylogenetic analysis based on chloroplast genome shows that L. quihoui and L. japonicum form a sister clade within the genus Ligustrum.

Resolving the Mortierellaceae phylogeny through synthesis of multi-gene phylogenetics and phylogenomics

Early efforts to classify Mortierellaceae were based on macro- and micromorphology, but sequencing and phylogenetic studies with ribosomal DNA (rDNA) markers have demonstrated conflicting taxonomic groupings and polyphyletic genera. Although some taxonomic confusion in the family has been clarified, rDNA data alone is unable to resolve higher level phylogenetic relationships within Mortierellaceae. In this study, we applied two parallel approaches to resolve the Mortierellaceae phylogeny: low coverage genome (LCG) sequencing and high-throughput, multiplexed targeted amplicon sequencing to generate sequence data for multi-gene phylogenetics. We then combined our datasets to provide a well-supported genome-based phylogeny having broad sampling depth from the amplicon dataset. Resolving the Mortierellaceae phylogeny into monophyletic groups led to the definition of 14 genera, 7 of which are newly proposed. Low-coverage genome sequencing proved to be a relatively cost-effective means of generating a well-resolved phylogeny. The multi-gene phylogenetics approach enabled much greater sampling depth and breadth than the LCG approach, but was unable to resolve higher-level organization of groups. We present this work to resolve some of the taxonomic confusion and provide a genus-level framework to empower future studies on Mortierellaceae diversity, biology, and evolution.

Phylogenomics indicates the "living fossil" Isoetes diversified in the Cenozoic

The fossil record provides an invaluable insight into the temporal origins of extant lineages of organisms. However, establishing the relationships between fossils and extant lineages can be difficult in groups with low rates of morphological change over time. Molecular dating can potentially circumvent this issue by allowing distant fossils to act as calibration points, but rate variation across large evolutionary scales can bias such analyses. In this study, we apply multiple dating methods to genome-wide datasets to infer the origin of extant species of Isoetes, a group of mostly aquatic and semi-aquatic isoetalean lycopsids, which closely resemble fossil forms dating back to the Triassic. Rate variation observed in chloroplast genomes hampers accurate dating, but genome-wide nuclear markers place the origin of extant diversity within this group in the mid-Paleogene, 45-60 million years ago. Our genomic analyses coupled with a careful evaluation of the fossil record indicate that despite resembling forms from the Triassic, extant Isoetes species do not represent the remnants of an ancient and widespread group, but instead have spread around the globe in the relatively recent past.

Museomics for reconstructing historical floristic exchanges: Divergence of stone oaks across Wallacea

Natural history collections and tropical tree diversity are both treasure troves of biological and evolutionary information, but their accessibility for scientific study is impeded by a number of properties. DNA in historical specimens is generally highly fragmented, complicating the recovery of high-grade genetic material. Furthermore, our understanding of hyperdiverse, wide-spread tree assemblages is obstructed by extensive species ranges, fragmented knowledge of tropical tree diversity and phenology, and a widespread lack of species-level diagnostic characters, prohibiting the collecting of readily identifiable specimens which can be used to build, revise or strengthen taxonomic frameworks. This, in turn, delays the application of downstream conservation action. A sizable component of botanical collections are sterile-thus eluding identification and are slowing down progress in systematic treatments of tropical biodiversity. With rapid advances in genomics and bioinformatic approaches to biodiversity research, museomics is emerging as a new field breathing life into natural collections that have been built up over centuries. Using MIGseq (multiplexed ISSR genotyping by sequencing), we generated 10,000s of short loci, for both freshly collected materials and museum specimens (aged >100 years) of Lithocarpus-a widespread tropical tree genus endemic to the Asian tropics. Loci recovery from historical and recently collected samples was not affected by sample age and preservation history of the study material, underscoring the reliability and flexibility of the MIGseq approach. Phylogenomic inference and biogeographic reconstruction across insular Asia, highlights repeated migration and diversification patterns between continental regions and islands. Results indicate that co-occurring insular species at the extremity of the distribution range are not monophyletic, raising the possibility of multiple independent dispersals along the outer edge of Wallacea. This suggests that dispersal of large seeded tree genera throughout Malesia and across Wallacea may have been less affected by large geographic distances and the presence of marine barriers than generally assumed. We demonstrate the utility of MIGseq in museomic studies using non-model taxa, presenting the first range-wide genomic assessment of Lithocarpus and tropical Fagaceae as a proof-of-concept. Our study shows the potential for developing innovative genomic approaches to improve the capture of novel evolutionary signals using valuable natural history collections of hyperdiverse taxa.

Ancient Plant Genomics in Archaeology, Herbaria, and the Environment

The ancient DNA revolution of the past 35 years has driven an explosion in the breadth, nuance, and diversity of questions that are approachable using ancient biomolecules, and plant research has been a constant, indispensable facet of these developments. Using archaeological, paleontological, and herbarium plant tissues, researchers have probed plant domestication and dispersal, plant evolution and ecology, paleoenvironmental composition and dynamics, and other topics across related disciplines. Here, we review the development of the ancient DNA discipline and the role of plant research in its progress and refinement. We summarize our understanding of long-term plant DNA preservation and the characteristics of degraded DNA. In addition, we discuss challenges in ancient DNA recovery and analysis and the laboratory and bioinformatic strategies used to mitigate them. Finally, we review recent applications of ancient plant genomic research.

Paternity tests support a diallelic self-incompatibility system in a wild olive (Olea europaea subsp. laperrinei, Oleaceae)

Self-incompatibility (SI) is the main mechanism that favors outcrossing in plants. By limiting compatible matings, SI interferes in fruit production and breeding of new cultivars. In the Oleeae tribe (Oleaceae), an unusual diallelic SI system (DSI) has been proposed for three distantly related species including the olive (Olea europaea), but empirical evidence has remained controversial for this latter. The olive domestication is a complex process with multiple origins. As a consequence, the mixing of S-alleles from two distinct taxa, the possible artificial selection of self-compatible mutants and the large phenological variation of blooming may constitute obstacles for deciphering SI in olive. Here, we investigate cross-genotype compatibilities in the Saharan wild olive (O. e. subsp. laperrinei). As this taxon was geographically isolated for thousands of years, SI should not be affected by human selection. A population of 37 mature individuals maintained in a collection was investigated. Several embryos per mother were genotyped with microsatellites in order to identify compatible fathers that contributed to fertilization. While the pollination was limited by distance inside the collection, our results strongly support the DSI hypothesis, and all individuals were assigned to two incompatibility groups (G1 and G2). No self-fertilization was observed in our conditions. In contrast, crosses between full or half siblings were frequent (ca. 45%), which is likely due to a nonrandom assortment of related trees in the collection. Finally, implications of our results for orchard management and the conservation of olive genetic resources are discussed.

Late Miocene origin and recent population collapse of the Malagasy savanna olive tree (Noronhia lowryi)

Debates regarding the origin of tropical savannas have attempted to disentangle the role of human, biotic and abiotic factors. Understanding the origins of savanna remains essential to identifying processes that gave rise to habitat mosaics, particularly those found in the Central Plateau of Madagascar. Documenting the evolutionary history and demography of native trees occurring in open habitats may reveal footprints left by past and recent environmental changes. We conducted a population genetic analysis of an endangered Malagasy shrub (Noronhia lowryi, Oleaceae) of the Central Plateau. Seventy-seven individuals were sampled from three sites and genotyped at 14 nuclear and 24 chloroplast microsatellites. We found a highly contrasting nuclear and plastid genetic structure, suggesting that pollen-mediated gene flow allows panmixia, while seed-based dispersal may rarely exceed tens of metres. From a phylogeny based on full plastomes, we dated the surprisingly old crown age of maternal lineages back to ~6.2 Mya, perhaps co-occurring with the global expansion of savanna. In contrast, recent demographic history inferred from nuclear data shows a bottleneck signature ~350 generations ago, probably reflecting an environmental shift during the Late Pleistocene or the Holocene. Ancient in situ adaptation and recent demographic collapse of an endangered woody plant highlight the unique value and vulnerability of the Malagasy savannas.