Phylogenetic informativeness analyses to clarify past diversification processes in Cucurbitaceae

The plastid genome characters and phylogenetic status of the endemic species Trichosanthes sunhangii D. G. Zhang, Z. M. Li, Qun Liu & T. Deng 2021 (Cucurbitaceae) in the Shennongjia forestry district of China

Trichosanthes sunhangii (Cucurbitaceae) is an endemic species native to the Shennongjia forestry district of China, whose plastid genome was reported in this study. The whole genome exhibits the typical quadripartite structure with 156,906 bp in size. A total of 130 genes were identified, containing 85 protein-coding genes (CDS), 37 tRNA, and 8 rRNA genes. Phylogenetic reconstruction based on 83 shared CDS sequences reaffirmed the status of T. sunhangii within the Sect. Foliobracteola, revealing close relationships with morphologically similar species, T. kirilowii and T. rosthornii. Our findings will provide a significant foundation for future investigations into the evolution, conservation, and potential utilization of this species.

The complete chloroplast genome sequence of Melothria scabra (Cucurbitaceae)

Melothria scabra has gradually become an economically important plant worldwide. The complete chloroplast genome of M. scabra has a length of 156,744 bp, contains a large single-copy (LSC) region (86,387 bp), a small single-copy (SSC) region (18,055 bp), and two inverted repeats (IRs) with the same length of 26,151 bp. In total, 126 genes were detected, including 83 protein-encoding genes, 35 transfer RNA (tRNA) genes, and eight ribosomal RNA (rRNA) genes. For phylogenetic analysis, M. scabra has a closer genetic relationship with Cucumis sativus and Citrullus lanatus. The complete chloroplast genome sequence of M. scabra would promote the germplasm exploration, phylogenetic relationships, and molecular biology researches in Melothria.

The complete plastome sequence of Datisca cannabina L. (Datiscaceae, Cucurbitales)

In this study, we report the first complete plastome sequence of Datisca cannabina (GenBank acc. no. OP432690). The plastome had a typical quadripartite structure. Its size was 162,914 bp, consisting of 90,890 bp large single-copy (LSC), 19,296 bp small single-copy (SSC), and 26,364 bp inverted repeat (IR) regions. It contained 112 genes, including 78 protein-coding, 30 tRNA, and four rRNA genes. The infA gene was pseudogenized. Sixteen genes contain one intron and two genes (clpP and ycf3) had two introns. Our phylogenetic tree showed that D. cannabina formed a close relationship with Begoniaceae. However, further samples are required to determine the phylogenetic placement of Datiscaceae in Cucurbitales.

The complete chloroplast genome sequence of Thladiantha nudiflora Hemsl. ex F.B.Forbes & Hemsl. 1887 (Cucurbitaceae)

Thladiantha nudiflora Hemsl. ex F.B.Forbes & Hemsl. 1887 (Cucurbitaceae) has been widely known as a traditional medicine plant. In this study, we sequenced, assembled, and annotated the complete chloroplast genome of T. nudiflora. The chloroplast genome of T. nudiflora is 156,824 base pair (bp) in length, containing a large single-copy region of 86,566 bp and a small single-copy region of 18,070 bp, separated by a pair of inverted repeats of 26,094 bp. The chloroplast genome contains 132 genes, including 87 protein-coding, 37 transfer RNA, and eight ribosomal RNA genes. Phylogenetic analysis of the chloroplast genome revealed that species of the genus Thladiantha were clustered together in the phylogenetic trees. This study will not only shed light on T. nudiflora's evolutionary position but also provide valuable chloroplast genomic information for future studies into the origins and diversification of the genus Thladiantha and the Cucurbitaceae family.

Resolving the phylogeny of Thladiantha (Cucurbitaceae) with three different target capture pipelines

BACKGROUND

Despite recent advances, reliable tools to simultaneously handle different types of sequencing data (e.g., target capture, genome skimming) for phylogenomics are still scarce. Here, we evaluate the performance of the recently developed pipeline Captus in comparison with the well-known target capture pipelines HybPiper and SECAPR. As test data, we analyzed newly generated sequences for the genus Thladiantha (Cucurbitaceae) for which no well-resolved phylogeny estimate has been available so far, as well as simulated reads derived from the genome of Arabidopsis thaliana.

RESULTS

Our pipeline comparisons are based on (1) the time needed for data assembly and locus extraction, (2) locus recovery per sample, (3) the number of informative sites in nucleotide alignments, and (4) the topology of the nuclear and plastid phylogenies. Additionally, the simulated reads derived from the genome of Arabidopsis thaliana were used to evaluate the accuracy and completeness of the recovered loci. In terms of computation time, locus recovery per sample, and informative sites, Captus outperforms HybPiper and SECAPR. The resulting topologies of Captus and SECAPR are identical for coalescent trees but differ when trees are inferred from concatenated alignments. The HybPiper phylogeny is similar to Captus in both methods. The nuclear genes recover a deep split of Thladiantha in two clades, but this is not supported by the plastid data.

CONCLUSIONS

Captus is the best choice among the three pipelines in terms of computation time and locus recovery. Even though there is no significant topological difference between the Thladiantha species trees produced by the three pipelines, Captus yields a higher number of gene trees in agreement with the topology of the species tree (i.e., fewer genes in conflict with the species tree topology).

The complete plastome sequence of Momordica cochinchinensis (Cucurbitaceae)

Momordica cochinchinensis (Lour.) Spreng. is an important medicinal plant that is used to treat various diseases in South and Southeast Asia. In this study, the complete plastome of M. cochinchinensis was sequenced and found to exhibit a total length of 158,955 bp, with a large single copy (LSC) region of 87,924 bp and a small single copy (SSC) region of 18,479 bp, as well as with two inverted repeats (IRs) that were both 26,726 bp in length. In total, 129 genes were detected, comprising 86 protein-encoding genes, 8 ribosomal RNA (rRNA) genes, and 35 transfer RNA (tRNA) genes. Furthermore, the inferred phylogenetic tree confirmed that M. cochinchinensis belongs to the genus Momordica in the Cucurbitaceae family. The research results will be used for authenticating M. cochinchinensis plant materials and for analyzing the genetic diversity and phylogenetic relationships in Momordica.

A Practical Guide to Design and Assess a Phylogenomic Study

Over the last decade, molecular systematics has undergone a change of paradigm as high-throughput sequencing now makes it possible to reconstruct evolutionary relationships using genome-scale datasets. The advent of "big data" molecular phylogenetics provided a battery of new tools for biologists but simultaneously brought new methodological challenges. The increase in analytical complexity comes at the price of highly specific training in computational biology and molecular phylogenetics, resulting very often in a polarized accumulation of knowledge (technical on one side and biological on the other). Interpreting the robustness of genome-scale phylogenetic studies is not straightforward, particularly as new methodological developments have consistently shown that the general belief of "more genes, more robustness" often does not apply, and because there is a range of systematic errors that plague phylogenomic investigations. This is particularly problematic because phylogenomic studies are highly heterogeneous in their methodology, and best practices are often not clearly defined. The main aim of this article is to present what I consider as the ten most important points to take into consideration when planning a well-thought-out phylogenomic study and while evaluating the quality of published papers. The goal is to provide a practical step-by-step guide that can be easily followed by nonexperts and phylogenomic novices in order to assess the technical robustness of phylogenomic studies or improve the experimental design of a project.

A critical review on cytogenetics of Cucurbitaceae with updates on Indian taxa

The cytogenetic relationships in the species of Cucurbitaceae are becoming immensely important to answer questions pertaining to genome evolution. Here, a simplified and updated data resource on cytogenetics of Cucurbitaceae is presented on the basis of foundational parameters (basic, zygotic and gametic chromosome numbers, ploidy, genome size, karyotype) and molecular cytogenetics. We have revised and collated our own findings on seven agriculturally important Indian cucurbit species in a comparative account with the globally published reports. Chromosome count (of around 19% species) shows nearly three-fold differences while genome size (of nearly 5% species) shows 5.84-fold differences across the species. There is no significant correlation between chromosome numbers and nuclear genome sizes. The possible trend of evolution is discussed here based on molecular cytogenetics data, especially the types and distribution of nucleolus organizer regions (NORs). The review supersedes the scopes of general chromosome databases and invites scopes for continuous updates. The offline resource serves as an exclusive toolkit for research and breeding communities across the globe and also opens scope for future establishment of web-database on Cucurbitaceae cytogenetics.

Plastid phylogenomics of the Gynoxoid group (Senecioneae, Asteraceae) highlights the importance of motif-based sequence alignment amid low genetic distances

PREMISE

The genus Gynoxys and relatives form a species-rich lineage of Andean shrubs and trees with low genetic distances within the sunflower subtribe Tussilaginineae. Previous molecular phylogenetic investigations of the Tussilaginineae have included few, if any, representatives of this Gynoxoid group or reconstructed ambiguous patterns of relationships for it.

METHODS

We sequenced complete plastid genomes of 21 species of the Gynoxoid group and related Tussilaginineae and conducted detailed comparisons of the phylogenetic relationships supported by the gene, intron, and intergenic spacer partitions of these genomes. We also evaluated the impact of manual, motif-based adjustments of automatic DNA sequence alignments on phylogenetic tree inference.

RESULTS

Our results indicate that the inclusion of all plastid genome partitions is needed to infer well-supported phylogenetic trees of the Gynoxoid group. Whole plastome-based tree inference suggests that the genera Gynoxys and Nordenstamia are polyphyletic and form the core clade of the Gynoxoid group. This clade is sister to a clade of Aequatorium and Paragynoxys and also includes some but not all representatives of Paracalia.

CONCLUSIONS

The concatenation and combined analysis of all plastid genome partitions and the construction of manually-curated, motif-based DNA sequence alignments are found to be instrumental in the recovery of well-supported relationships of the Gynoxoid group. We demonstrate that the correct assessment of homology in genome-level plastid sequence data sets is crucial for subsequent phylogeny reconstruction and that the manual post-processing of multiple sequence alignments improves the reliability of such reconstructions amid low genetic distances between taxa.

Phylogenomic Subsampling and the Search for Phylogenetically Reliable Loci

Phylogenomic subsampling is a procedure by which small sets of loci are selected from large genome-scale data sets and used for phylogenetic inference. This step is often motivated by either computational limitations associated with the use of complex inference methods or as a means of testing the robustness of phylogenetic results by discarding loci that are deemed potentially misleading. Although many alternative methods of phylogenomic subsampling have been proposed, little effort has gone into comparing their behavior across different data sets. Here, I calculate multiple gene properties for a range of phylogenomic data sets spanning animal, fungal, and plant clades, uncovering a remarkable predictability in their patterns of covariance. I also show how these patterns provide a means for ordering loci by both their rate of evolution and their relative phylogenetic usefulness. This method of retrieving phylogenetically useful loci is found to be among the top performing when compared with alternative subsampling protocols. Relatively common approaches such as minimizing potential sources of systematic bias or increasing the clock-likeness of the data are found to fare worse than selecting loci at random. Likewise, the general utility of rate-based subsampling is found to be limited: loci evolving at both low and high rates are among the least effective, and even those evolving at optimal rates can still widely differ in usefulness. This study shows that many common subsampling approaches introduce unintended effects in off-target gene properties and proposes an alternative multivariate method that simultaneously optimizes phylogenetic signal while controlling for known sources of bias.

Phylogenomic Subsampling and the Search for Phylogenetically Reliable Loci

Phylogenomic subsampling is a procedure by which small sets of loci are selected from large genome-scale data sets and used for phylogenetic inference. This step is often motivated by either computational limitations associated with the use of complex inference methods or as a means of testing the robustness of phylogenetic results by discarding loci that are deemed potentially misleading. Although many alternative methods of phylogenomic subsampling have been proposed, little effort has gone into comparing their behavior across different data sets. Here, I calculate multiple gene properties for a range of phylogenomic data sets spanning animal, fungal, and plant clades, uncovering a remarkable predictability in their patterns of covariance. I also show how these patterns provide a means for ordering loci by both their rate of evolution and their relative phylogenetic usefulness. This method of retrieving phylogenetically useful loci is found to be among the top performing when compared with alternative subsampling protocols. Relatively common approaches such as minimizing potential sources of systematic bias or increasing the clock-likeness of the data are found to fare worse than selecting loci at random. Likewise, the general utility of rate-based subsampling is found to be limited: loci evolving at both low and high rates are among the least effective, and even those evolving at optimal rates can still widely differ in usefulness. This study shows that many common subsampling approaches introduce unintended effects in off-target gene properties and proposes an alternative multivariate method that simultaneously optimizes phylogenetic signal while controlling for known sources of bias.

Plastid Phylogenomic Data Offers Novel Insights Into the Taxonomic Status of the Trichosanthes kirilowii Complex (Cucurbitaceae) in South Korea

Trichosanthes is a genus in Cucurbitaceae comprising 90-100 species. Trichosanthes species are valuable as herbaceous medicinal ingredients. The fruits, seeds, and roots of species such as T. kirilowii and T. rosthornii are used in Korean traditional herbal medicines. T. rosthornii is only found in China, whereas in South Korea two varieties, T. kirilowii var. kirilowii and T. kirilowii var. japonica, are distributed. T. kirilowii var. kirilowii and T. kirilowii var. japonica have different fruit and leaf shapes but are recognized as belonging to the same species. Furthermore, although its members have herbal medicine applications, genomic information of the genus is still limited. The broad goals of this study were (i) to evaluate the taxonomy of Trichosanthes using plastid phylogenomic data and (ii) provide molecular markers specific for T. kirilowii var. kirilowii and T. kirilowii var. japonica, as these have differences in their pharmacological effectiveness and thus should not be confused and adulterated. Comparison of five Trichosanthes plastid genomes revealed locally divergent regions, mainly within intergenic spacer regions (trnT-UGU-trnL-UAA: marker name Tri, rrn4.5-rrn5: TRr, trnE-UUC-trnT-GGU: TRtt). Using these three markers as DNA-barcodes for important herbal medicine species in Trichosanthes, the identity of Trichosanthes material in commercial medicinal products in South Korea could be successfully determined. Phylogenetic analysis of the five Trichosanthes species revealed that the species are clustered within tribe Sicyoeae. T. kirilowii var. kirilowii and T. rosthornii formed a clade with T. kirilowii var. japonica as their sister group. As T. kirilowii in its current circumscription is paraphyletic and as the two varieties can be readily distinguished morphologically (e.g., in leaf shape), T. kirilowii var. japonica should be treated (again) as an independent species, T. japonica.

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

PREMISE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Whole plastomes are not enough: phylogenomic and morphometric exploration at multiple demographic levels of the bee orchid clade Ophrys sect. Sphegodes

Plastid sequences have long dominated phylogeny reconstruction at all time depths, predicated on a usually untested assumption that they accurately represent the evolutionary histories of phenotypically circumscribed species. We combined detailed in situ morphometrics (124 plants) and whole-plastome sequencing through genome skimming (71 plants) in order to better understand species-level diversity and speciation in arguably the most challenging monophyletic group within the taxonomically controversial, pseudo-copulatory bee orchid genus Ophrys. Using trees and ordinations, we interpreted the data at four nested demographic levels-macrospecies, mesospecies, microspecies, and local population-seeking the optimal level for bona fide species. Neither morphological nor molecular discontinuities are evident at any level below macrospecies, the observed overlap among taxa suggesting that both mesospecies and microspecies reflect arbitrary division of a continuum of variation. Plastomes represent geographic location more strongly than taxonomic assignment and correlate poorly with morphology, suggesting widespread plastid capture and possibly post-glacial expansion from multiple southern refugia. As they are rarely directly involved in the speciation process, plastomes depend on extinction of intermediate lineages to provide phylogenetic signal and so cannot adequately document evolutionary radiations. The popular 'ethological' evolutionary model recognizes as numerous 'ecological species' (microspecies) lineages perceived as actively diverging as a result of density-dependent selection on very few features that immediately dictate extreme pollinator specificity. However, it is assumed rather than demonstrated that the many microspecies are genuinely diverging. We conversely envisage a complex four-dimensional reticulate network of lineages, generated locally and transiently through a wide spectrum of mechanisms, but each unlikely to maintain an independent evolutionary trajectory long enough to genuinely speciate by escaping ongoing gene flow. The frequent but localized microevolution that characterizes the Ophrys sphegodes complex is often convergent and rarely leads to macroevolution. Choosing between the contrasting 'discontinuity' and 'ethology' models will require next-generation sequencing of nuclear genomes plus ordination of corresponding morphometric matrices, seeking the crucial distinction between retained ancestral polymorphism-consistent with lineage divergence-and polymorphisms reflecting gene flow through 'hybridization'-more consistent with lineage convergence.

Transcriptome-wide analysis of the AP2/ERF transcription factor gene family involved in the regulation of gypenoside biosynthesis in Gynostemma pentaphyllum

Gynostemma pentaphyllum is a traditional Chinese medicinal herb, serving as natural source of gypenosides (triterpene saponins). The APETALA2/ethylene response factor (AP2/ERF) transcription factors, playing essential regulation roles in plant biotic and abiotic stress responses and secondary metabolism biosynthesis. However, the regulation roles of AP2/ERF transcription factors in gypenosides biosynthesis in G. pentaphyllum remains little understood. In the present study, 125 AP2/ERF genes were identified from G. pentaphyllum transcriptome datasets. Phylogenetic, conserved motifs and expression pattern were employed to comprehensively analyze the 125 GpAP2/ERF genes. Based on the sequence similarity and phylogeny tree, the 125 GpAP2/ERF genes can be classified into 10 groups. Moreover, the distribution of conserved motifs among GpAP2/ERF proteins in phylogenetic trees was consistent with previous studies, thus supporting the classification. Expression profiling indicated that the 125 GpAP2/ERF genes exhibited distinct tissue-specific expression patterns. As confirmed by qRT-PCR, the four candidate GpAP2/ERF genes and gypenoside biosynthetic genes were highly expressed in leaves and/or flowers, and show similar expression patterns in response to MeJA. Base on the expression patterns and phylogenetic relationships, two GpAP2/ERF genes were considered as potential regulatory genes for gypenoside biosynthesis. Our study enhances understanding roles of GpAP2/ERF genes in regulation of gypenosides biosynthesis.