A genome-scale mining strategy for recovering novel rapidly-evolving nuclear single-copy genes for addressing shallow-scale phylogenetics in Hydrangea

Placing human gene families into their evolutionary context

Following the draft sequence of the first human genome over 20 years ago, we have achieved unprecedented insights into the rules governing its evolution, often with direct translational relevance to specific diseases. However, staggering sequence complexity has also challenged the development of a more comprehensive understanding of human genome biology. In this context, interspecific genomic studies between humans and other animals have played a critical role in our efforts to decode human gene families. In this review, we focus on how the rapid surge of genome sequencing of both model and non-model organisms now provides a broader comparative framework poised to empower novel discoveries. We begin with a general overview of how comparative approaches are essential for understanding gene family evolution in the human genome, followed by a discussion of analyses of gene expression. We show how homology can provide insights into the genes and gene families associated with immune response, cancer biology, vision, chemosensation, and metabolism, by revealing similarity in processes among distant species. We then explain methodological tools that provide critical advances and show the limitations of common approaches. We conclude with a discussion of how these investigations position us to gain fundamental insights into the evolution of gene families among living organisms in general. We hope that our review catalyzes additional excitement and research on the emerging field of comparative genomics, while aiding the placement of the human genome into its existentially evolutionary context.

Molecular Phylogeny, Character Evolution, and Biogeography of Hydrangea Section Cornidia, Hydrangeaceae

Background: Hydrangea section Cornidia consists of 26 currently accepted species and a yet undefined number of new species and erroneously synonymized taxa. This clade consists of (sub)tropical lianas occurring from northern Mexico to southern Chile and Argentina, and one species from Southeast Asia. Currently, no molecular phylogenetic hypothesis is available that includes more than a few species of this section. Hence, a resolved and well-sampled molecular phylogenetic hypothesis may help to enforce taxonomic decisions. In this study, we present a phylogenetic framework based on sequences from two low copy nuclear genes from a comprehensive taxon sampling of H. section Cornidia and a selection of outgroups. Our phylogenetic reconstructions prove the non-monophyly of the traditionally recognized subsections Monosegia and Polysegia and their corresponding series, Speciosae and Aphananthae, and Synstyleae and Chorystyleae, respectively. Three morphologically defined species were recovered with high support as monophyletic, namely, Hydrangea panamensis, Hydrangea serratifolia, and Hydrangea tarapotensis. However, statistical support for some shallow nodes did not allow to refute, with high support, the monophyly of several of the herein recognized species for which more than one individual could be analyzed. Based on the obtained phylogenetic framework, we reconstructed the evolution of selected reproductive characters. Hydrangea section Cornidia is the only genus section for which dioecism has been extensively documented. Our character reconstruction of sexual dimorphism shows that dioecism is the ancestral state in this section and that this was reversed to monoecy in Hydrangea seemannii and Hydrangea integrifolia. Character reconstruction for the enlarged marginal flowers recovered their presence as the ancestral character state in H. section Cornidia, although at least three internal lineages independently lost them; thus, losses were reconstructed to be more likely than gain. With respect to the flower color, more species exhibit white than red flowers, and white is reconstructed as the ancestral state. Cornidia also shows an unusual disjunct geographic distribution between Asia and Central Mesoamerica-South America, as it is not present in the USA and Canada. The origin of Cornidia is reconstructed to be the New World with higher probability, and the presence of one species in Asia is likely due to long-distance dispersal.

Draft genome of Puya raimondii (Bromeliaceae), the Queen of the Andes

Puya raimondii, the Queen of the Andes, is an endangered high Andean species in the Bromeliaceae family. Here, we report its first genome to promote its conservation and evolutionary study. Comparative genomics showed P. raimondii diverged from Ananas comosus about 14.8 million years ago, and the long terminal repeats were likely to contribute to the genus diversification in last 3.5 million years. The gene families related to plant reproductive development and stress responses significantly expanded in the genome. At the same time, gene families involved in disease defense, photosynthesis and carbohydrate metabolism significantly contracted, which may be an evolutionary strategy to adapt to the harsh conditions in high Andes. The demographic history analysis revealed the P. raimondii population size sharply declined in the Pleistocene and then increased in the Holocene. We also designed and tested 46 pairs of universal primers for amplifying orthologous single-copy nuclear genes in Puya species.

On Hydrangea peruviana, an endangered species from Ecuador, and Hydrangea oerstedii, very common in Costa Rica and Panama, and seven threatened Central and South American Hydrangeas, which have been confounded with these

Hydrangea section Cornidia, currently consisting of 19 accepted taxa, occurs from Mexico to Chile and Argentina, with one species in southeast Asia. Its representatives are root-climbing lianas which may grow up to 60 m high in the tree canopy of temperate to (sub)tropical forests. Our extensive field work throughout its distribution area, study of herbarium specimens and ongoing molecular studies have resulted in the discovery of species new to science, as well as new insights into the circumscription of many taxa. We here present amended descriptions for seven Hydrangea species of Central and South America and discuss the taxonomical situation of two Colombian Hydrangeas, including an identification key, illustrations, and distribution maps. Field work was carried out in Costa Rica, Panama, Ecuador and Peru, including exploration in areas where the genus had not been collected before. These specimens and observations were complemented with the study of specimens of 41 herbaria of North, Central and South America, as well as Europe. Detailed morphological studies of all species were carried out, based on living plants in their natural habitat, as well as on dried specimens from our own collections and all available herbarium material. Type material was studied in detail for all species concerned. Based on an extensive number of morphological characters, combined with distribution patterns, phenological differences and ecological preferences, including molecular data in most cases, Hydrangea peruviana and H. oerstedii are clearly distinct taxa, as well as the other seven species mentioned here, which had been synonymized with either of these two species. The present study results in the recognition of 26 species in section Cornidia and exemplifies the urgent need for profound taxonomic studies in plants, as in many families we do not dispose of well-circumscribed units for conservation to mitigate the already occurring unprecedented loss of biodiversity.

Targeted generation of polyploids in Hydrangea macrophylla through cross-based breeding

Background

Up to now, diploid and triploid cultivars were reported for the ornamental crop Hydrangea macrophylla. Especially, the origin of triploids and their crossing behaviors are unknown, but the underlying mechanisms are highly relevant for breeding polyploids.

Results

By screening a cultivar collection, we identified diploid, triploid, tetraploid and even aneuploid H. macrophylla varieties. The pollen viability of triploids and tetraploids was comparable to that of diploids. Systematic crosses with these cultivars resulted in viable diploid, triploid, tetraploid and aneuploid offspring. Interestingly, crosses between diploids produced diploid and 0 or 1–94% triploid offspring, depending on the cultivars used as pollen parent. This finding suggests that specific diploids form unreduced pollen, either at low or high frequencies. In contrast, crosses of triploids with diploids or tetraploids produced many viable aneuploids, whose 2C DNA contents ranged between the parental 2C values. As expected, crosses between diploid and tetraploid individuals generated triploid offspring. Putative tetraploid plants were obtained at low frequencies in crosses between diploids and in interploid crosses of triploids with either diploid or tetraploid plants. The analysis of offspring populations indicated the production of 1n = 2x gametes for tetraploid plants, whereas triploids produced obviously reduced, aneuploid gametes with chromosome numbers ranging between haploid and diploid level. While euploid offspring grew normally, aneuploid plants showed mostly an abnormal development and a huge phenotypic variation within offspring populations, most likely due to the variation in chromosome numbers. Subsequent crosses with putative diploid, triploid and aneuploid offspring plants from interploid crosses resulted in viable offspring and germination rates ranging from 21 to 100%.

Conclusions

The existence of diploids that form unreduced pollen and of tetraploids allows the targeted breeding of polyploid H. macrophylla. Different ploidy levels can be addressed by combining the appropriate crossing partners. In contrast to artificial polyploidization, cross-based polyploidization is easy, cheap and results in genetically variable offspring that allows the direct selection of more robust and stress tolerant polyploid varieties. Furthermore, the generation of polyploid H. macrophylla plants will favor interspecific breeding programs within the genus Hydrangea.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-020-00954-z.

Unraveling Reticulate Evolution in Opuntia (Cactaceae) From Southern Mexico

The process of hybridization occurs in approximately 40% of vascular plants, and this exchange of genetic material between non-conspecific individuals occurs unequally among plant lineages, being more frequent in certain groups such as Opuntia (Cactaceae). This genus is known for multiple taxonomic controversies due to widespread polyploidy and probable hybrid origin of several of its species. Southern Mexico species of this genus have been poorly studied despite their great diversity in regions such as the Tehuacán-Cuicatlán Valley which contains around 12% of recognized Mexico's native Opuntia species. In this work, we focus on testing the hybrid status of two putative hybrids from this region, Opuntia tehuacana and Opuntia pilifera, and estimate if hybridization occurs among sampled southern opuntias using two newly identified nuclear intron markers to construct phylogenetic networks with HyDe and Dsuite and perform invariant analysis under the coalescent model with HyDe and Dsuite. For the test of hybrid origin in O. tehuacana, our results could not recover hybridization as proposed in the literature, but we found introgression into O. tehuacana individuals involving O. decumbens and O. huajuapensis. Regarding O. pilifera, we identified O. decumbens as probable parental species, supported by our analysis, which sustains the previous hybridization hypothesis between Nopalea and Basilares clades. Finally, we suggest new hybridization and introgression cases among southern Mexican species involving O. tehuantepecana and O. depressa as parental species of O. velutina and O. decumbens.

Identification and assessment of variable single-copy orthologous (SCO) nuclear loci for low-level phylogenomics: a case study in the genus Rosa (Rosaceae)

BACKGROUND

With an ever-growing number of published genomes, many low levels of the Tree of Life now contain several species with enough molecular data to perform shallow-scale phylogenomic studies. Moving away from using just a few universal phylogenetic markers, we can now target thousands of other loci to decipher taxa relationships. Making the best possible selection of informative sequences regarding the taxa studied has emerged as a new issue. Here, we developed a general procedure to mine genomic data, looking for orthologous single-copy loci capable of deciphering phylogenetic relationships below the generic rank. To develop our strategy, we chose the genus Rosa, a rapid-evolving lineage of the Rosaceae family in which several species genomes have recently been sequenced. We also compared our loci to conventional plastid markers, commonly used for phylogenetic inference in this genus.

RESULTS

We generated 1856 sequence tags in putative single-copy orthologous nuclear loci. Associated in silico primer pairs can potentially amplify fragments able to resolve a wide range of speciation events within the genus Rosa. Analysis of parsimony-informative site content showed the value of non-coding genomic regions to obtain variable sequences despite the fact that they may be more difficult to target in less related species. Dozens of nuclear loci outperform the conventional plastid phylogenetic markers in terms of phylogenetic informativeness, for both recent and ancient evolutionary divergences. However, conflicting phylogenetic signals were found between nuclear gene tree topologies and the species-tree topology, shedding light on the many patterns of hybridization and/or incomplete lineage sorting that occur in the genus Rosa.

CONCLUSIONS

With recently published genome sequence data, we developed a set of single-copy orthologous nuclear loci to resolve species-level phylogenomics in the genus Rosa. This genome-wide scale dataset contains hundreds of highly variable loci which phylogenetic interest was assessed in terms of phylogenetic informativeness and topological conflict. Our target identification procedure can easily be reproduced to identify new highly informative loci for other taxonomic groups and ranks.

Plant phylogenomics based on genome-partitioning strategies: Progress and prospects

The rapid expansion of next-generation sequencing (NGS) has generated a powerful array of approaches to address fundamental questions in biology. Several genome-partitioning strategies to sequence selected subsets of the genome have emerged in the fields of phylogenomics and evolutionary genomics. In this review, we summarize the applications, advantages and limitations of four NGS-based genome-partitioning approaches in plant phylogenomics: genome skimming, transcriptome sequencing (RNA-seq), restriction site associated DNA sequencing (RAD-Seq), and targeted capture (Hyb-seq). Of these four genome-partitioning approaches, targeted capture (especially Hyb-seq) shows the greatest promise for plant phylogenetics over the next few years. This review will aid researchers in their selection of appropriate genome-partitioning approaches to address questions of evolutionary scale, where we anticipate continued development and expansion of whole-genome sequencing strategies in the fields of plant phylogenomics and evolutionary biology research.

Topological support and data quality can only be assessed through multiple tests in reviewing Blattodea phylogeny

Assessing support for molecular phylogenies is difficult because the data is heterogeneous in quality and overwhelming in quantity. Traditionally, node support values (bootstrap frequency, Bayesian posterior probability) are used to assess confidence in tree topologies. Other analyses to assess the quality of phylogenetic data (e.g. Lento plots, saturation plots, trait consistency) and the resulting phylogenetic trees (e.g. internode certainty, parameter permutation tests, topological tests) exist but are rarely applied. Here we argue that a single qualitative analysis is insufficient to assess support of a phylogenetic hypothesis and relate data quality to tree quality. We use six molecular markers to infer the phylogeny of Blattodea and apply various tests to assess relationship support, locus quality, and the relationship between the two. We use internode-certainty calculations in conjunction with bootstrap scores, alignment permutations, and an approximately unbiased (AU) test to assess if the molecular data unambiguously support the phylogenetic relationships found. Our results show higher support for the position of Lamproblattidae, high support for the termite phylogeny, and low support for the position of Anaplectidae, Corydioidea and phylogeny of Blaberoidea. We use Lento plots in conjunction with mutation-saturation plots, calculations of locus homoplasy to assess locus quality, identify long branch attraction, and decide if the tree's relationships are the result of data biases. We conclude that multiple tests and metrics need to be taken into account to assess tree support and data robustness.

Genome skimming provides new insight into the relationships in Ludwigia section Macrocarpon, a polyploid complex

PREMISE OF THE STUDY

Interpreting relationships within groups containing polyploids, which are frequent in angiosperms, can be greatly assisted by genomic techniques. In this study, we used a genome-skimming approach to investigate the evolutionary relationships and origins of polyploids in the monophyletic group, Ludwigia section Macrocarpon (Onagraceae), which includes diploid, tetraploid, and hexaploid taxa.

METHODS

We sampled all known taxa and ploidy levels in the section and conducted shotgun sequencing. We assembled plastomes, mitochondrial sequences, and completed nuclear ribosomal regions, reconstructed phylogenies, and conducted comparative genomic analyses for plastomes to gain insights into the relationships among studied taxa.

KEY RESULTS

Within the section, results showed that the South American diploid taxa L. bonariensis and L. lagunae were closely related. We reported the first chromosome count (2n = 4× = 32) for L. neograndiflora, which is closely related to the two South American diploid taxa, although its exact origin remains unclear. The samples of the widespread, polyploid taxon L. octovalvis do not form a monophyletic group. Both tetraploid and hexaploid L. octovalvis lineages have originated more than once. At least one tetraploid in the L. octovalvis lineage may have been involved in the origins of hexaploids. One or more extinct/unsampled intermediate tetraploids in the L. octovalvis lineages had also likely been involved in the origins of hexaploids.

CONCLUSIONS

Genome skimming provided important insights into the complex evolutionary relationships within sect. Macrocarpon, but additional sampling and data from single-copy nuclear regions are necessary to further elucidate the origins of the polyploids in this section.

Comparison of taxon-specific versus general locus sets for targeted sequence capture in plant phylogenomics

PREMISE OF THE STUDY

Targeted sequence capture can be used to efficiently gather sequence data for large numbers of loci, such as single-copy nuclear loci. Most published studies in plants have used taxon-specific locus sets developed individually for a clade using multiple genomic and transcriptomic resources. General locus sets can also be developed from loci that have been identified as single-copy and have orthologs in large clades of plants.

METHODS

We identify and compare a taxon-specific locus set and three general locus sets (conserved ortholog set [COSII], shared single-copy nuclear [APVO SSC] genes, and pentatricopeptide repeat [PPR] genes) for targeted sequence capture in Buddleja (Scrophulariaceae) and outgroups. We evaluate their performance in terms of assembly success, sequence variability, and resolution and support of inferred phylogenetic trees.

RESULTS

The taxon-specific locus set had the most target loci. Assembly success was high for all locus sets in Buddleja samples. For outgroups, general locus sets had greater assembly success. Taxon-specific and PPR loci had the highest average variability. The taxon-specific data set produced the best-supported tree, but all data sets showed improved resolution over previous non-sequence capture data sets.

DISCUSSION

General locus sets can be a useful source of sequence capture targets, especially if multiple genomic resources are not available for a taxon.

More on the Best Evolutionary Rate for Phylogenetic Analysis

The accumulation of genome-scale molecular data sets for nonmodel taxa brings us ever closer to resolving the tree of life of all living organisms. However, despite the depth of data available, a number of studies that each used thousands of genes have reported conflicting results. The focus of phylogenomic projects must thus shift to more careful experimental design. Even though we still have a limited understanding of what are the best predictors of the phylogenetic informativeness of a gene, there is wide agreement that one key factor is its evolutionary rate; but there is no consensus as to whether the rates derived as optimal in various analytical, empirical, and simulation approaches have any general applicability. We here use simulations to infer optimal rates in a set of realistic phylogenetic scenarios with varying tree sizes, numbers of terminals, and tree shapes. Furthermore, we study the relationship between the optimal rate and rate variation among sites and among lineages. Finally, we examine how well the predictions made by a range of experimental design methods correlate with the observed performance in our simulations.We find that the optimal level of divergence is surprisingly robust to differences in taxon sampling and even to among-site and among-lineage rate variation as often encountered in empirical data sets. This finding encourages the use of methods that rely on a single optimal rate to predict a gene's utility. Focusing on correct recovery either of the most basal node in the phylogeny or of the entire topology, the optimal rate is about 0.45 substitutions from root to tip in average Yule trees and about 0.2 in difficult trees with short basal and long-apical branches, but all rates leading to divergence levels between about 0.1 and 0.5 perform reasonably well.Testing the performance of six methods that can be used to predict a gene's utility against our simulation results, we find that the probability of resolution, signal-noise analysis, and Fisher information are good predictors of phylogenetic informativeness, but they require specification of at least part of a model tree. Likelihood quartet mapping also shows very good performance but only requires sequence alignments and is thus applicable without making assumptions about the phylogeny. Despite them being the most commonly used methods for experimental design, geometric quartet mapping and the integration of phylogenetic informativeness curves perform rather poorly in our comparison. Instead of derived predictors of phylogenetic informativeness, we suggest that the number of sites in a gene that evolve at near-optimal rates (as inferred here) could be used directly to prioritize genes for phylogenetic inference. In combination with measures of model fit, especially with respect to compositional biases and among-site and among-lineage rate variation, such an approach has the potential to greatly improve marker choice and should be tested on empirical data.

Development of single-copy nuclear intron markers for species-level phylogenetics: Case study with Paullinieae (Sapindaceae)

PREMISE OF THE STUDY

We developed a bioinformatic pipeline that leverages a publicly available genome and published transcriptomes to design primers in conserved coding sequences flanking targeted introns of single-copy nuclear loci. Paullinieae (Sapindaceae) is used to demonstrate the pipeline.

METHODS AND RESULTS

Transcriptome reads phylogenetically closer to the lineage of interest are aligned to the closest genome. Single-nucleotide polymorphisms are called, generating a "pseudoreference" closer to the lineage of interest. Several filters are applied to meet the criteria of single-copy nuclear loci with introns of a desired size. Primers are designed in conserved coding sequences flanking introns. Using this pipeline, we developed nine single-copy nuclear intron markers for Paullinieae.

CONCLUSIONS

This pipeline is highly flexible and can be used for any group with available genomic and transcriptomic resources. This pipeline led to the development of nine variable markers for phylogenetic study without generating sequence data de novo.

Multilocus coalescent species delimitation to evaluate traditionally defined morphotypes in Hydrangea sect. Asperae (Hydrangeaceae)

The number of species recognized in section Asperae of the flowering plant genus Hydrangea differs widely between subsequent revisions. This variation is largely centered around the H. aspera species complex, with numbers of recognized species varying from one to nearly a dozen. Despite indications of molecular variation in this complex, no sequence-based species delimitation methods have been employed to evaluate the primarily morphology-based species boundaries. In the present study, a multi-locus coalescent-based approach to species delimitation is employed in order to identify separate evolutionary lines within H. sect. Asperae, using four chloroplast and four nuclear molecular markers. Eight lineages were recovered within the focal group, of which five correspond with named morphotypes. The other three lineages illustrate types of conflict between molecular species delimitation and traditional morphology-based taxonomy. One molecular lineage comprises two named morphotypes, which possibly diverged recently enough to not have developed sufficient molecular divergence. A second conflict is found in H. strigosa. This morphotype is recovered as a separate lineage when occurring in geographic isolation, but when occurring in sympatry with two other morphotypes (H. aspera and H. robusta), the coalescent species delimitation lumps these taxa into a single putative species.

Evolution of Epiphytism and Fruit Traits Act Unevenly on the Diversification of the Species-Rich Genus Peperomia (Piperaceae)

The species-rich genus Peperomia (Black Pepper relatives) is the only genus among early diverging angiosperms where epiphytism evolved. The majority of fruits of Peperomia release sticky secretions or exhibit hook-shaped appendages indicative of epizoochorous dispersal, which is in contrast to other flowering plants, where epiphytes are generally characterized by fruit morphological adaptations for anemochory or endozoochory. We investigate fruit characters using Cryo-SEM. Comparative phylogenetic analyses are applied for the first time to include life form and fruit character information to study diversification in Peperomia. Likelihood ratio tests uncover correlated character evolution. We demonstrate that diversification within Peperomia is not homogenous across its phylogeny, and that net diversification rates increase by twofold within the most species-rich subgenus. In contrast to former land plant studies that provide general evidence for increased diversification in epiphytic lineages, we demonstrate that the evolution of epiphytism within Peperomia predates the diversification shift. An epiphytic-dependent diversification is only observed for the background phylogeny. An elevated frequency of life form transitions between epiphytes and terrestrials and thus evolutionary flexibility of life forms is uncovered to coincide with the diversification shift. The evolution of fruits showing dispersal related structures is key to diversification in the foreground region of the phylogeny and postdates the evolution of epiphytism. We conclude that the success of Peperomia, measured in species numbers, is likely the result of enhanced vertical and horizontal dispersal ability and life form flexibility but not the evolution of epiphytism itself.