Molecular phylogeny and inflorescence evolution of Prunus (Rosaceae) based on RAD-seq and genome skimming analyses

Different reference genomes determine different results: Comparing SNP calling in RAD-seq of Engelhardia roxburghiana using different reference genomes

Advances in next-generation sequencing (NGS) have significantly reduced the cost and improved the efficiency of obtaining single nucleotide polymorphism (SNP) markers, particularly through restriction site-associated DNA sequencing (RAD-seq). Meanwhile, the progression in whole genome sequencing has led to the utilization of an increasing number of reference genomes in SNP calling processes. This study utilized RAD-seq data from 242 individuals of Engelhardia roxburghiana, a tropical tree of the walnut family (Juglandaceae), with SNP calling conducted using the STACKS pipeline. We aimed to compare both reference-based approaches, namely, employing a closely related species as the reference genome versus the species itself as the reference genome, to evaluate their respective merits and limitations. Our findings indicate a substantial discrepancy in the number of obtained SNPs between using a closely related species as opposed to the species itself as reference genomes, the former yielded approximately an order of magnitude fewer SNPs compared to the latter. While the missing rate of individuals and sites of the final SNPs obtained in the two scenarios showed no significant difference. The results showed that using the reference genome of the species itself tends to be prioritized in RAD-seq studies. However, if this is unavailable, considering closely related genomes is feasible due to their wide applicability and low missing rate as alternatives. This study contributes to enrich the understanding of the impact of SNP acquisition when utilizing different reference genomes.

Assessment of phylogenetic relationships and genetic diversity of Sagittaria trifolia using phenotypic traits and SNP markers

The aquatic perennial herb Sagittaria trifolia L. commonly known as arrowhead, has been utilized in China both as a culinary vegetable and in traditional medicines. Characterizing the phylogenetic relationships and genetic diversity of arrowheads is crucial for improved management, conservation, and efficient utilization of the germplasm resources associated with this species. Herein, we presented the phenotypic traits and genome-wide DNA marker-based analyses of 111 arrowhead accessions, most of which were from China. Cluster analysis revealed that arrowhead could be categorized into two clusters based on 11 phenotypic traits, with Cluster 1 comprising two subclusters. All accessions were clustered into three sub-clusters based primarily on leaf shape and tuber weight. A set of 759,237 high-quality single-nucleotide polymorphisms was identified and used to assess the phylogenetic relationships. Population structure and phylogenetic tree analyses suggested that the accessions could be classified into two major groups, Group I was further subdivided into two subgroups, aligning with the clusters identified through morphological classification. By employing Sagittaria lichuanensis as an outgroup, the rooted tree revealed that the evolutionary relationships within the three groups followed a progression from Group I-1 to Group I-2 and finally to Group II. The landraces were clustered into one group along with the remaining wild accessions. The level of genetic diversity for Group I (π = 0.26) was slightly lower than that which was estimated for Group II (π = 0.29). The lowest pairwise differentiation levels (Fst, 0.008) were obtained from the comparison between groups I-2 and II, indicating that the two groups were the most closely related. This study provides novel insights into germplasm classification, evolutionary relationships, genomics and arrowhead breeding.

The pan-plastome of Prunus mume: insights into Prunus diversity, phylogeny, and domestication history

Backgrounds

Prunus mume in the Rosaceae and commonly referred to as mei or Chinese plum is widely used as a traditional ornamental flowering plant and fruit tree in China. Although some population and genetic analyses have been conducted for this species, no extensive comparisons of genetic variation from plastomes have yet been investigated.

Methods

We de novo assembled a total of 322 complete P. mume plastomes in this study and did a series of comparative analyses to better resolve pan-plastomic patterns of P. mume. To determine the phylogeny and domestication history of this species, we reconstructed the phylogenetic tree of Prunus genus, and resolved the population structure of P. mume. We also examined the nucleotide variation of P. mume to find potential DNA barcodes.

Results

The assembled plastomes exhibited a typical quadripartite structure and ranged from 157,871 bp to 158,213 bp in total size with a GC content ranging from 36.73 to 36.75%. A total of 112 unique genes were identified. Single nucleotide variants (SNVs) were the most common variants found among the plastomes, followed by nucleotide insertions/deletions (InDels), and block substitutions with the intergenic spacer (IGS) regions containing the greatest number of variants. From the pan-plastome data six well-supported genetic clusters were resolved using multiple different population structure analyses. The different cultivars were unevenly distributed among multiple clades. We also reconstructed a phylogeny for multiple species of Prunus to better understand genus level diversity and history from which a complex introgressive relationship between mei and other apricots/plums was resolved.

Conclusion

This study constructed the pan-plastome of P. mume, which indicated the domestication of P. mume involved multiple genetic origins and possible matrilineal introgression from other species. The phylogenetic analysis in Prunus and the population structure of P. mume provide an important maternal history for Prunus and the groundwork for future studies on intergenomic sequence transfers, cytonuclear incompatibility, and conservation genetics.

Central Asia revealed as a key area in evolution of Eremurus (Asphodelaceae)

Eremurus was described at the beginning of the 19th century. However, due to limited sampling and the small number of gene markers to date, its phylogeny and evolution are largely unknown. In this study, we analyzed plastomes from 27 species belonging to 2 subgenera and 3 sections of Eremurus, which are found in Central Asia (its center of diversity) and China. We also analyzed nuclear DNA ITS of 33 species, encompassing all subgenera and sections of the genus in Central Asia, southwest Asia and China. Our findings revealed that the genus was monophyletic, although both subgenera Eremurus and Henningia were found to be paraphyletic. Both plastome and nrDNA-based phylogenetic trees had three clades that did not reflect the current taxonomy of the genus. Our biogeographical and time-calibrated trees suggest that Eremurus originated in the ancient Tethyan area in the second half of the Eocene. Diversification of Eremurus occurred from the early Oligocene to the late Miocene. Paratethys Sea retreat and several orogenetic events, such as the progressive uplift of the Qinghai-Tibet Plateau and surrounding mountain belts (Altai, Pamir, Tian Shan), caused serious topographic and climate (aridification) changes in Central Asia that may have triggered a split of clades and speciation. In this transformed Central Asia, speciation proceeded rapidly driven mainly by vicariance caused by numerous mountain chains and specialization to a variety of climatic, topographic and soil conditions that exist in this region.

Molecular characterization and evolutionary relationships of DOFs in four cherry species and functional analysis in sweet cherry

The DOF (DNA binding with one finger) has multiple functions in plants. However, it has received little attention in the research field of cherries. In this study, the evolutionary relationship and molecular characterization of DOF in four cherry species were analyzed, revealing its expression pattern in sweet cherry. There are 23 members in Prunus avium cv. 'Tieton', 88 in Prunus cerasus, 53 in Cerasus × yedoensis, and 27 in Cerasus serrulata. Most of these genes are intron-less or non-intron, with a conserved C2-C2 domain. Due to heterozygosity and chromosomal ploidy, whole-genome duplication (WGD) events occur to varying degrees, and DOF genes are contracted during evolution. Furthermore, these genes are affected by purifying selection pressure. Under low-temperature treatment, the expression of PavDOF2 and PavDOF18 were significantly up-regulated, while that of PavDOF16 is significantly down-regulated. The expression of PavDOF9, PavDOF12, PavDOF14, PavDOF16, PavDOF17, PavDOF18, and PavDOF19 exhibits an increasing trend during flower development and varies during sweet cherry fruit development. PavDOF1, PavDOF8, PavDOF9, and PavDOF15 are localized in the nucleus but is not transcriptionally active. The findings systemically demonstrate the molecular characteristics of DOF in different cherry varieties, providing a basis for further research on the functions of these genes.

Historical biogeography and evolutionary diversification of Lilium (Liliaceae): New insights from plastome phylogenomics

Here, we infer the historical biogeography and evolutionary diversification of the genus Lilium. For this purpose, we used the complete plastomes of 64 currently accepted species in the genus Lilium (14 plastomes were newly sequenced) to recover the phylogenetic backbone of the genus and a time-calibrated phylogenetic framework to estimate biogeographical history scenarios and evolutionary diversification rates of Lilium. Our results suggest that ancient climatic changes and geological tectonic activities jointly shaped the distribution range and drove evolutionary radiation of Lilium, including the Middle Miocene Climate Optimum (MMCO), the late Miocene global cooling, as well as the successive uplift of the Qinghai-Tibet Plateau (QTP) and the strengthening of the monsoon climate in East Asia during the late Miocene and the Pliocene. This case study suggests that the unique geological and climatic events in the Neogene of East Asia, in particular the uplift of QTP and the enhancement of monsoonal climate, may have played an essential role in formation of uneven distribution of plant diversity in the Northern Hemisphere.

The complete chloroplast genome sequence of Prunus salicina cultivar 'Zuili' (Rosaceae)

'Zuili' is a distinguished plum (Prunus salicina Lindley 1830) originating from Tongxiang City, Zhejiang Province, China, and a nationally recognized geographical indication product. In this study, we reported the complete chloroplast genome of P. salicina cultivar 'Zuili'. The genome has a circular structure of 157,935 bp containing a large single-copy (LSC) region of 86,133 bp, a small single-copy (SSC) region of 19,028 bp, and two inverted repeats (IRs) of 26,387 by each. It harbors 130 genes (111 unique genes), including 85 protein-coding genes (78 are unique), eight ribosomal RNA genes (four are unique), and 37 transfer genes (29 are unique). The phylogenetic analysis based on whole chloroplast genomes showed 'Zuili' was clustered with Prunus salicina cultivar 'Wuyuecui' (MW406461.1) and 'No. 2 Guofeng' (MW406472.1). This study provides valuable information that can contribute to the identification and further evolutionary analysis of Prunus salicina cultivar 'Zuili'.

Evolution of Cherries (Prunus Subgenus Cerasus) Based on Chloroplast Genomes

Cherries (Prunus Subgenus Cerasus) have economic value and ecological significance, yet their phylogeny, geographic origin, timing, and dispersal patterns remain challenging to understand. To fill this gap, we conducted a comprehensive analysis of the complete chloroplast genomes of 54 subg. Cerasus individuals, along with 36 additional genomes from the NCBI database, resulting in a total of 90 genomes for comparative analysis. The chloroplast genomes of subg. Cerasus exhibited varying sizes and consisted of 129 genes, including protein-coding, transfer RNA, and ribosomIal RNA genes. Genomic variation was investigated through InDels and SNPs, showcasing distribution patterns and impact levels. A comparative analysis of chloroplast genome boundaries highlighted variations in inverted repeat (IR) regions among Cerasus and other Prunus species. Phylogeny based on whole-chloroplast genome sequences supported the division of Prunus into three subgenera, I subg. Padus, II subg. Prunus and III subg. Cerasus. The subg. Cerasus was subdivided into seven lineages (IIIa to IIIg), which matched roughly to taxonomic sections. The subg. Padus first diverged 51.42 Mya, followed by the separation of subg. Cerasus from subg. Prunus 39.27 Mya. The subg. Cerasus started diversification at 15.01 Mya, coinciding with geological and climatic changes, including the uplift of the Qinghai-Tibet Plateau and global cooling. The Himalayans were the refuge of cherries, from which a few species reached Europe through westward migration and another species reached North America through northeastward migration. The mainstage of cherry evolution was on the Qing-Tibet Plateau and later East China and Japan as well. These findings strengthen our understanding of the evolution of cherry and provide valuable insights into the conservation and sustainable utilization of cherry's genetic resources.