Divergence in the Aquilegia ecalcarata complex is correlated with geography and climate oscillations: Evidence from plastid genome data

Unraveling the extensive phylogenetic discordance and evolutionary history of spurless taxa within the Aquilegia ecalcarata complex

Parallel evolution of the same, or at least very similar, phenotype(s) in different lineages is often interpreted as evidence for the action of natural selection. However, caution is required when inferring parallel evolution based on uncertain or potentially incorrect phylogenetic frameworks. Here, by conducting extensive phylogenomic and population genetic analyses, we aim to clarify the evolutionary history of spurless taxa within the Aquilegia ecalcarata complex. We observed substantial discordance in the phylogenetic patterns across the entire genome, primarily attributed to ancient introgression and incomplete lineage sorting. Additionally, we identified several spurless lineages whose phylogenetic positions were distorted by admixture events. Using a backbone tree and demographic modeling, we determined that these spurless taxa independently originated twice within this group. Intriguingly, our investigation revealed that the spurless taxa experienced population expansion during global cooling, while their spurred sister groups underwent population contraction. The parallel losses of petal spurs, therefore, may be linked to adaptations for low-temperature conditions. These findings emphasize the importance of comprehensive population-level analyses in phylogenetic inference and provide valuable insights into the dynamics of trait loss and its implications for the adaptive strategies.

Standing genetic variation and introgression shape the cryptic radiation of Aquilegia in the mountains of Southwest China

Cryptic diversity in evolutionary radiation offers an excellent system for investigating the intricacies of evolutionary progress. Understanding the evolution of cryptic diversity is imperative for unraveling the hidden complexities of biodiversity. However, empirical evidence elucidating the mechanisms behind cryptic radiation remains limited, particularly in plants. Here, we focus on a monophyletic group of Aquilegia species mainly distributed in the mountains of Southwest China, one of the world's biodiversity hotspots. Using whole-genome resequencing of 158 individuals from 23 natural populations, we identify three to four paraphyletic lineages within each morphological species. Our findings reveal that 39 out of 43 detected instances of introgression occurred post-lineage formation. Identifying shared genomic regions indicates that the divergence of fixed singletons in lineages from morphological species A. kansuensis and A. rockii predates lineage formation, supporting a scenario where incomplete lineage sorting of standing variation contributes to morphological parallelism. Furthermore, strong positive correlations among genomic differentiation, divergence, and introgression suggest that standing variations and introgression from non-sister lineages contribute to the rapid genetic divergence. Our study illuminates the important roles of standing variations and introgression in plant cryptic radiation, advancing our understanding of the complex mechanisms behind the evolution of biodiversity in recent radiation events.

Multiple Ecological Niche Modeling Reveals Niche Conservatism and Divergence in East Asian Yew (Taxus)

Understanding ecological niche evolution patterns is crucial for elucidating biogeographic history and guiding biodiversity conservation. Taxus is a Tertiary relict gymnosperm with 11 lineages mainly distributed across East Asia, spanning from tropical to subarctic regions. However, the spatiotemporal dynamics of its ecological niche evolution and the roles of ecological and geographical factors in lineage diversification, remain unclear. Using occurrence records, environmental data, and reconstructed phylogenies, we employed ensemble ecological niche models (eENMs), environmental principle components analysis (PCA-env), and phyloclimatic modeling to analyze niche similarity and evolution among 11 Taxus lineages. Based on reconstructed Bayesian trees and geographical distribution characteristics, we classified the eleven lineages into four clades: Northern (T. cuspidata), Central (T. chinensis, T. qinlingensis, and the Emei type), Western (T. wallichiana, T. florinii, and T. contorta), and Southern (T. calcicola, T. phytonii, T. mairei, and the Huangshan type). Orogenic activities and climate changes in the Tibetan Plateau since the Late Miocene likely facilitated the local adaptation of ancestral populations in Central China, the Hengduan Mountains, and the Yunnan-Guizhou Plateau, driving their expansion and diversification towards the west and south. Key environmental variables, including extreme temperature, temperature and precipitation variability, light, and altitude, were identified as major drivers of current niche divergence. Both niche conservatism and divergence were observed, with early conservatism followed by recent divergence. The Southern clade exhibits high heat and moisture tolerance, suggesting an adaptive shift, while the Central and Western clades retain ancestral drought and cold tolerance, displaying significant phylogenetic niche conservatism (PNC). We recommend prioritizing the conservation of T. qinlingensis, which exhibits the highest PNC level, particularly in the Qinling, Daba, and Taihang Mountains, which are highly degraded and vulnerable to future climate fluctuations.

Comparative plastomes sheds light on phylogeny of Weigela

Weigela Thunb. is a genus in the family Caprifoliaceae. All species in this genus have high ornamental and medicinal value. However, the genetic divergence between species and the phylogeny within Weigela is still unclear. Therefore, we sequenced and analyzed four plastomes from four different Weigela species to reveal the genetic divergence among species of this genus, and the phylogeny within Weigela. The four plastomes from Weigela ranged from 156,909 bp to 157,739 bp in size, and presented a typical circular quadripartite structure. Each complete plastome contained a pair of inverted repeat regions (23,592~24,957 bp), a larger single-copy (LSC) region (89,922~90,229 bp), and a small single-copy (SSC) region (17,668~20,429 bp). We identified three types of repeats, corresponding to 268 forward repeats, 128 palindromic repeats, and 867 tandem repeats, for a total of 1,263 long repeats. A total of 352 SSRs were identified from the four plastomes, and most of them were concentrated in the LSC region and the noncoding regions. Mononucleotide repeat units were the most frequently detected types of repeats, of which A/T repeat units were the most abundant. Three mutational hotspots (trnH-psbA, trnR-ndhF, and trnN-ndhF) were identified as candidate barcodes for Weigela species. Weigela belongs to Diervilloideae located at an early diverging position in the Caprifoliaceae. Within Weigela, W. japonica and W. floribunda were sister with W. subsessilis and W. florida. This study revealed the plastome structure and variation of four well-known Weigela species, and found three candidate barcodes for further study of four well-known Weigela species. In addition, the phylogenetic location of Weigela within the Caprifoliaceae was identified.

Phylogeny, Genetic Diversity and Population Structure of Fritillaria cirrhosa and Its Relatives Based on Chloroplast Genome Data

Fritillaria cirrhosa and its relatives have been utilized in traditional Chinese medicine for many years and are under priority protection in China. Despite their medicinal and protective value, research on their phylogeny, genetic diversity, and divergence remains limited. Here, we investigate the chloroplast genome variation architecture of 46 samples of F. cirrhosa and its relatives collected from various regions, encompassing the majority of wild populations across diverse geographical areas. The results indicate abundant variations in 46 accessions including 1659 single-nucleotide polymorphisms and 440 indels. Six variable markers (psbJ, ndhD, ycf1, ndhG, trnT-trnL, and rpl32-trnL) were identified. Phylogenetic and network analysis, population structure analysis, and principal component analysis showed that the 46 accessions formed five clades with significant divergence, which were related to their geographical distribution. The regions spanning from the southern Hengduan Mountains to the Qinghai-Tibet Plateau exhibited the highest levels of genetic diversity. F. cirrhosa and its relatives may have suffered a genetic bottleneck and have a relatively low genetic diversity level. Moreover, geographical barriers and discrete patches may have accelerated population divergence. The study offers novel perspectives on the phylogeny, genetic diversity, and population structure of F. cirrhosa and its relatives, information that can inform conservation and utilization strategies in the future.

Metal nanoparticle-induced effects on green toads (Amphibia, Anura) under climate change: conservation implications

The toxicity of aluminum oxide (Al2O3), copper oxide (CuO), iron oxide (Fe3O4), nickel oxide (NiO), zinc oxide (ZnO), and titanium dioxide (TiO2) nanoparticles (NPs) on amphibians and their interaction with high temperatures, remain unknown. In this study, we investigated the survival, developmental, behavioral, and histological reactions of Bufotes viridis embryos and larvae exposed to different NPs for a duration of 10 days, using lethal concentrations (LC25%, LC50%, and LC75% mg/L) under both ambient (AT: 18 °C) and high (HT: 21 °C) temperatures. Based on LC, NiONPs > ZnONPs > CuONPs > Al2O3NPs > TiO2NPs > Fe3O4NPs showed the highest mortality at AT. A similar pattern was observed at HT, although mortality occurred at lower concentrations and Fe3O4NPs were more toxic than TiO2NPs. The results indicated that increasing concentrations of NPs significantly reduced hatching rates, except for TiO2NPs. Survival rates decreased, abnormality rates increased, and developmental processes slowed down, particularly for NiONPs and ZnONPs, under HT conditions. However, exposure to low concentrations of Fe3O4NPs for up to 7 days, CuONPs for up to 72 h, and NiO, ZnONPs, and TiO2NPs for up to 96 h did not have a negative impact on survival compared with the control group under AT. In behavioral tests with larvae, NPs generally induced hypoactivity at AT and hyperactivity at HT. Histological findings revealed liver and internal gill tissue lesions, and an increase in the number of melanomacrophage centers at HT. These results suggest that global warming may exacerbate the toxicity of metal oxide NPs to amphibians, emphasizing the need for further research and conservation efforts in this context.

Genomic divergence and demographic history of Quercus aliena populations

BACKGROUND

Quercus aliena is a major montane tree species of subtropical and temperate forests in China, with important ecological and economic value. In order to reveal the species' population dynamics, genetic diversity, genetic structure, and association with mountain habitats during the evolutionary process, we re-sequenced the genomes of 72 Q. aliena individuals.

RESULTS

The whole chloroplast and nuclear genomes were used for this study. Phylogenetic analysis using the chloroplast genome dataset supported four clades of Q. aliena, while the nuclear dataset supported three major clades. Sex-biased dispersal had a critical role in causing discordance between the chloroplast and nuclear genomes. Population structure analysis showed two groups in Q. aliena. The effective population size sharply declined 1 Mya, coinciding with the Poyang Glaciation in Eastern China. Using genotype-climate association analyses, we found a positive correlation between allele frequency variation in SNPs and temperature, suggesting the species has the capacity to adapt to changing temperatures.

CONCLUSION

Overall, this study illustrates the genetic divergence, genomic variation, and evolutionary processes behind the demographic history of Q. aliena.

Multiple paternally inherited chloroplast capture events associated with Taxus speciation in the Hengduan Mountains

Mountainous regions provide a multitude of habitats and opportunities for complex speciation scenarios. Hybridization leading to chloroplast capture, which can be revealed by incongruent phylogenetic trees, is one possible outcome. Four allopatric Taxus lineages (three species and an undescribed lineage) from the Hengduan Mountains, southwest China, exhibit conflicting phylogenetic relationships between nuclear and chloroplast phylogenies. Here, we use multi-omic data at the population level to investigate their historical speciation processes. Population genomic analysis based on ddRAD-seq data revealed limited contemporary inter-specific gene flow involving only populations located close to another species. In a historical context, chloroplast and nuclear data (transcriptome) consistently showed conflicting phylogenetic relationships for T. florinii and the Emei type lineage. ILS and chloroplast recombination were excluded as possible causes, and transcriptome and ddRAD-seq data revealed an absence of the mosaic nuclear genomes that characterize hybrid origin scenarios. Therefore, T. florinii appears to have originated when a lineage of T. florinii captured the T. chinensis plastid type, whereas plastid introgression in the opposite direction generated the Emei Type. All four species have distinct ecological niche based on community investigations and ecological niche analyses. We propose that the origins of both species represent very rare examples of chloroplast capture events despite the paternal cpDNA inheritance of gymnosperms. Specifically, allopatrically and/or ecologically diverged parental species experienced a rare secondary contact, subsequent hybridization and reciprocal chloroplast capture, generating two new lineages, each of which acquired a unique ecological niche. These events might have been triggered by orogenic activities of the Hengduan Mountains and an intensification of the Asian monsoon in the late Miocene, and may represent a scenario more common in these mountains than presently known.

Pueraria montana Population Structure and Genetic Diversity Based on Chloroplast Genome Data

Despite having a generally conserved structure, chloroplast genome data have been helpful for plant population genetics and evolution research. To mine Pueraria montana chloroplast genome variation architecture and phylogeny, we investigated the chloroplast variation architecture of 104 P. montana accessions from across China. P. montana's chloroplast genome showed high diversity levels, with 1674 variations, including 1118 single nucleotide polymorphisms and 556 indels. The intergenic spacers, psbZ-trnS and ccsA-ndhD, are the two mutation hotspot regions in the P. montana chloroplast genome. Phylogenetic analysis based on the chloroplast genome dataset supported four P. montana clades. P. montana variations were conserved among and within clades, which showed high gene flow levels. Most P. montana clades were estimated to have diverged at 3.82-5.17 million years ago. Moreover, the East Asian summer monsoon and South Asian summer monsoon may have accelerated population divergence. Our results show that chloroplast genome sequences were highly variable and can be used as molecular markers to assess genetic variation and relationships in P. montana.

Evolutionary and phylogenetic analyses of 11 Cerasus species based on the complete chloroplast genome

The subgenus Cerasus, one of the most important groups in the genus Prunus sensu lato, comprises over 100 species; however, the taxonomic classification and phylogenetic relationships of Cerasus remain controversial. Therefore, it is necessary to reconstruct the phylogenetic tree for known Cerasus species. Here, we report the chloroplast (cp) genome sequences of 11 Cerasus species to provide insight into evolution of the plastome. The cp genomes of the 11 Cerasus species (157,571-158,830 bp) displayed a typical quadripartite circular structure. The plastomes contain 115 unique genes, including 80 protein-coding genes, four ribosomal RNAs, and 31 transfer RNAs. Twenty genes were found to be duplicated in inverted repeats as well as at the boundary. The conserved non-coding sequences showed significant divergence compared with the coding regions. We found 12 genes and 14 intergenic regions with higher nucleotide diversity and more polymorphic sites, including matK, rps16, rbcL, rps16-trnQ, petN-psbM, and trnL-trnF. During cp plastome evolution, the codon profile has been strongly biased toward the use of A/T at the third base, and leucine and isoleucine codons appear the most frequently. We identified strong purifying selection on the rpoA, cemA, atpA, and petB genes; whereas ccsA, rps19, matK, rpoC2, ycf2 and ndhI showed a signature of possible positive selection during the course of Cerasus evolution. In addition, we further analyzed the phylogenetic relationships of these species with 57 other congenic related species.Through reconstructing the Cerasus phylogeny tree, we found that true cherry is similar to the flora of China forming a distinct group, from which P. mahaleb was separated as an independent subclade. Microcerasus was genetically closer to Amygdalus, Armeniaca, and Prunus (sensu stricto) than to members of true cherry, whereas P. japonica and P. tomentosa were most closely related to P. triloba and P. pedunculata. However, P. tianshanica formed a clade with P. cerasus, P. fruticosa, P. cerasus × P. canescens 'Gisela 6', and P. avium as a true cherry group. These results provide new insights into the plastome evolution of Cerasus, along with potential molecular markers and candidate DNA barcodes for further phylogenetic and phylogeographic analyses of Cerasus species.

Maternal Donor and Genetic Variation of Lagerstroemia indica Cultivars

Lagerstroemia indica L. is a well-known ornamental plant with large pyramidal racemes, long flower duration, and diverse colors and cultivars. It has been cultivated for nearly 1600 years and is essential for investigating the germplasm and assessing genetic variation to support international cultivar identification and breeding programs. In this study, 20 common Lagerstroemia indica cultivars from different varietal groups and flower morphologies, as well as multiple wild relative species, were analyzed to investigate the maternal donor of Lagerstroemia indica cultivars and to discover the genetic variation and relationships among cultivars based on plastome and nuclear ribosomal DNA (nrDNA) sequences. A total of 47 single nucleotide polymorphisms (SNPs) and 24 insertion/deletions (indels) were identified in the 20 L. indica cultivars' plastome and 25 SNPs were identified in the nrDNA. Phylogenetic analysis based on the plastome sequences showed that all the cultivars formed a clade with the species of L. indica, indicating that L. indica was the maternal donor of the cultivars. Population structure and PCA analyses supported two clades of cultivars, which exhibited significant genetic differences according to the plastome dataset. The results of the nrDNA supported that all 20 cultivars were divided into three clades and most of the cultivars had at least two genetic backgrounds and higher gene flow. Our results suggest that the plastome and nrDNA sequences can be used as molecular markers for assessing the genetic variation and relationships of L. indica cultivars.

Characterization of the Dicranostigma leptopodum chloroplast genome and comparative analysis within subfamily Papaveroideae

BACKGROUND

Dicranostigma leptopodum (Maxim.) Fedde is a perennial herb with bright yellow flowers, well known as "Hongmao Cao" for its medicinal properties, and is an excellent early spring flower used in urban greening. However, its molecular genomic information remains largely unknown. Here, we sequenced and analyzed the chloroplast genome of D. leptopodum to discover its genome structure, organization, and phylogenomic position within the subfamily Papaveroideae.

RESULTS

The chloroplast genome size of D. leptopodum was 162,942 bp, and D. leptopodum exhibited a characteristic circular quadripartite structure, with a large single-copy (LSC) region (87,565 bp), a small single-copy (SSC) region (18,759 bp) and a pair of inverted repeat (IR) regions (28,309 bp). The D. leptopodum chloroplast genome encoded 113 genes, including 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. The dynamics of the genome structures, genes, IR contraction and expansion, long repeats, and single sequence repeats exhibited similarities, with slight differences observed among the eight Papaveroideae species. In addition, seven interspace regions and three coding genes displayed highly variable divergence, signifying their potential to serve as molecular markers for phylogenetic and species identification studies. Molecular evolution analyses indicated that most of the genes were undergoing purifying selection. Phylogenetic analyses revealed that D. leptopodum formed a clade with the tribe Chelidonieae.

CONCLUSIONS

Our study provides detailed information on the D. leptopodum chloroplast genome, expanding the available genomic resources that may be used for future evolution and genetic diversity studies.

Variations in genetic diversity in cultivated Pistacia chinensis

Identification of the evolution history and genetic diversity of a species is important in the utilization of novel genetic variation in this species, as well as for its conservation. Pistacia chinensis is an important biodiesel tree crop in China, due to the high oil content of its fruit. The aim of this study was to uncover the genetic structure of P. chinensis and to investigate the influence of intraspecific gene flow on the process of domestication and the diversification of varieties. We investigated the genetic structure of P. chinensis, as well as evolution and introgression in the subpopulations, through analysis of the plastid and nuclear genomes of 39 P. chinensis individuals from across China. High levels of variation were detected in the P. chinensis plastome, and 460 intraspecific polymorphic sites, 104 indels and three small inversions were identified. Phylogenetic analysis and population structure using the plastome dataset supported five clades of P. chinensis. Population structure analysis based on the nuclear SNPs showed two groups, clearly clustered together, and more than a third of the total individuals were classified as hybrids. Discordance between the plastid and nuclear genomes suggested that hybridization events may have occurred between highly divergent samples in the P. chinensis subclades. Most of the species in the P. chinensis subclade diverged between the late Miocene and the mid-Pliocene. The processes of domestication and cultivation have decreased the genetic diversity of P. chinensis. The extensive variability and structuring of the P. chinensis plastid together with the nuclear genomic variation detected in this study suggests that much unexploited genetic diversity is available for improvement in this recently domesticated species.

Whole-genome survey and phylogenetic analysis of Gadus macrocephalus

Gadus macrocephalus (Pacific cod) is an economically important species on the northern coast of the Pacific. Although numerous studies on G. macrocephalus exist, there are few reports on its genomic data. Here, we used whole-genome sequencing data to elucidate the genomic characteristics and phylogenetic relationship of G. macrocephalus. From the 19-mer frequency distribution, the genome size was estimated to be 658.22 Mb. The heterozygosity, repetitive sequence content and GC content were approximately 0.62%, 27.50% and 44.73%, respectively. The draft genome sequences were initially assembled, yielding a total of 500,760 scaffolds (N50 = 3565 bp). A total of 789,860 microsatellite motifs were identified from the genomic data, and dinucleotide repeat was the most dominant simple sequence repeat motif. As a byproduct of whole-genome sequencing, the mitochondrial genome was assembled to investigate the evolutionary relationships between G. macrocephalus and its relatives. On the basis of 13 protein-coding gene sequences of the mitochondrial genome of Gadidae species, the maximum likelihood phylogenetic tree showed that complicated relationships and divergence times among Gadidae species. Demographic history analysis revealed changes in the G. macrocephalus population during the Pleistocene by using the pairwise sequentially Markovian coalescent model. These findings supplement the genomic data of G. macrocephalus, and make a valuable contribution to the whole-genome studies on G. macrocephalus.