Genomic insights on the contribution of balancing selection and local adaptation to the long-term survival of a widespread living fossil tree, Cercidiphyllum japonicum

Rapidly Evolved Genes in Three Reaumuria Transcriptomes and Potential Roles of Pentatricopeptide Repeat Superfamily Proteins in Endangerment of R. trigyna

Reaumuria genus (Tamaricaceae) is widely distributed across the desert and semi-desert regions of Northern China, playing a crucial role in the restoration and protection of desert ecosystems. Previous studies mainly focused on the physiological responses to environmental stresses; however, due to the limited availability of genomic information, the underlying mechanism of morphological and ecological differences among the Reaumuria species remains poorly understood. In this study, we presented the first catalog of expressed transcripts for R. kaschgarica, a sympatric species of xerophyte R. soongorica. We further performed the pair-wise transcriptome comparison to determine the conserved and divergent genes among R. soongorica, R. kaschgarica, and the relict recretohalophyte R. trigyna. Annotation of the 600 relatively conserved genes revealed that some common genetic modules are employed by the Reaumuria species to confront with salt and drought stresses in arid environment. Among the 250 genes showing strong signs of positive selection, eight pentatricopeptide repeat (PPR) superfamily protein genes were specifically identified, including seven PPR genes in the R. soongorica vs. R. trigyna comparison and one PPR gene in the R. kaschgarica vs. R. trigyna comparison, while the cyclin D3 gene was found in the R. soongorica vs. R. trigyna comparison. These findings suggest that genetic variations in PPR genes may affect the fertility system or compromise the extent of organelle RNA editing in R. trigyna. The present study provides valuable genomic information for R. kaschgarica and preliminarily reveals the conserved genetic bases for the abiotic stress adaptation and interspecific divergent selection in the Reaumuria species. The rapidly evolved PPR and cyclin D3 genes provide new insights on the endangerment of R. trigyna and the leaf length difference among the Reaumuria species.

Integrating evolutionary genomics of forest trees to inform future tree breeding amid rapid climate change

Global climate change is leading to rapid and drastic shifts in environmental conditions, posing threats to biodiversity and nearly all life forms worldwide. Forest trees serve as foundational components of terrestrial ecosystems and play a crucial and leading role in combating and mitigating the adverse effects of extreme climate events, despite their own vulnerability to these threats. Therefore, understanding and monitoring how natural forests respond to rapid climate change is a key priority for biodiversity conservation. Recent progress in evolutionary genomics, driven primarily by cutting-edge multi-omics technologies, offers powerful new tools to address several key issues. These include precise delineation of species and evolutionary units, inference of past evolutionary histories and demographic fluctuations, identification of environmentally adaptive variants, and measurement of genetic load levels. As the urgency to deal with more extreme environmental stresses grows, understanding the genomics of evolutionary history, local adaptation, future responses to climate change, and conservation and restoration of natural forest trees will be critical for research at the nexus of global change, population genomics, and conservation biology. In this review, we explore the application of evolutionary genomics to assess the effects of global climate change using multi-omics approaches and discuss the outlook for breeding of climate-adapted trees.

Chromosome-level genome assembly of American sweetgum (Liquidambar styraciflua, Altingiaceae)

The deciduous American sweetgum (Liquidambar styraciflua, Altingiaceae) is a popular ornamental and economically valuable tree renowned for its sweet-smelling bark resin, abundant volatile substances, and spectacular fall leaf color. However, the absence of a reference genome hinders thorough investigations into the mechanisms underlying phenotypic variation, secondary metabolite synthesis and adaptation, both in this species and other Liquidambar members. In this study, we sequenced and constructed a chromosome-level assembly of the L. styraciflua genome, covering 662.48 Mb with a scaffold N50 of 39.54 Mb, by integrating PacBio, Illumina and chromosome conformation capture data. We identified 58.83% of the genome sequences as repetitive elements and 25,713 protein-coding genes, 97.28% of which were functionally annotated. The genome sequencing reads, assembly and annotation data have been deposited in publicly available repositories. This high-quality genome assembly provides valuable resources for further evolutionary and functional genomic studies in American sweetgum and other Liquidambar species.

Genome sequences and population genomics reveal climatic adaptation and genomic divergence between two closely related sweetgum species

Understanding the genetic basis of population divergence and adaptation is an important goal in population genetics and evolutionary biology. However, the relative roles of demographic history, gene flow, and/or selective regime in driving genomic divergence, climatic adaptation, and speciation in non-model tree species are not yet fully understood. To address this issue, we generated whole-genome resequencing data of Liquidambar formosana and L. acalycina, which are broadly sympatric but altitudinally segregated in the Tertiary relict forests of subtropical China. We integrated genomic and environmental data to investigate the demographic history, genomic divergence, and climatic adaptation of these two sister species. We inferred a scenario of allopatric species divergence during the late Miocene, followed by secondary contact during the Holocene. We identified multiple genomic islands of elevated divergence that mainly evolved through divergence hitchhiking and recombination rate variation, likely fostered by long-term refugial isolation and recent differential introgression in low-recombination genomic regions. We also found some candidate genes with divergent selection signatures potentially involved in climatic adaptation and reproductive isolation. Our results contribute to a better understanding of how late Tertiary/Quaternary climatic change influenced speciation, genomic divergence, climatic adaptation, and introgressive hybridization in East Asia's Tertiary relict flora. In addition, they should facilitate future evolutionary, conservation genomics, and molecular breeding studies in Liquidambar, a genus of important medicinal and ornamental values.

A Diverging Species within the Stewartia gemmata (Theaceae) Complex Revealed by RAD-Seq Data

Informed species delimitation is crucial in diverse biological fields; however, it can be problematic for species complexes. Showing a peripatric distribution pattern, Stewartia gemmata and S. acutisepala (the S. gemmata complex) provide us with an opportunity to study species boundaries among taxa undergoing nascent speciation. Here, we generated genomic data from representative individuals across the natural distribution ranges of the S. gemmata complex using restriction site-associated DNA sequencing (RAD-seq). Based on the DNA sequence of assembled loci containing 41,436 single-nucleotide polymorphisms (SNPs) and invariant sites, the phylogenetic analysis suggested strong monophyly of both the S. gemmata complex and S. acutisepala, and the latter was nested within the former. Among S. gemmata individuals, the one sampled from Mt. Tianmu (Zhejiang) showed the closest evolutionary affinity with S. acutisepala (which is endemic to southern Zhejiang). Estimated from 2996 high-quality SNPs, the genetic divergence between S. gemmata and S. acutisepala was relatively low (an Fst of 0.073 on a per-site basis). Nevertheless, we observed a proportion of genomic regions showing relatively high genetic differentiation on a windowed basis. Up to 1037 genomic bins showed an Fst value greater than 0.25, accounting for 8.31% of the total. After SNPs subject to linkage disequilibrium were pruned, the principal component analysis (PCA) showed that S. acutisepala diverged from S. gemmata along the first and the second PCs to some extent. By applying phylogenomic analysis, the present study determines that S. acutisepala is a variety of S. gemmata and is diverging from S. gemmata, providing empirical insights into the nascent speciation within a species complex.

Cryptic divergences and repeated hybridizations within the endangered "living fossil" dove tree (Davidia involucrata) revealed by whole genome resequencing

The identification and understanding of cryptic intraspecific evolutionary units (lineages) are crucial for planning effective conservation strategies aimed at preserving genetic diversity in endangered species. However, the factors driving the evolution and maintenance of these intraspecific lineages in most endangered species remain poorly understood. In this study, we conducted resequencing of 77 individuals from 22 natural populations of Davidia involucrata, a "living fossil" dove tree endemic to central and southwest China. Our analysis revealed the presence of three distinct local lineages within this endangered species, which emerged approximately 3.09 and 0.32 million years ago. These divergence events align well with the geographic and climatic oscillations that occurred across the distributional range. Additionally, we observed frequent hybridization events between the three lineages, resulting in the formation of hybrid populations in their adjacent as well as disjunct regions. These hybridizations likely arose from climate-driven population expansion and/or long-distance gene flow. Furthermore, we identified numerous environment-correlated gene variants across the total and many other genes that exhibited signals of positive evolution during the maintenance of two major local lineages. Our findings shed light on the highly dynamic evolution underlying the remarkably similar phenotype of this endangered species. Importantly, these results not only provide guidance for the development of conservation plans but also enhance our understanding of evolutionary past for this and other endangered species with similar histories.

The chromosome-scale genome of Magnolia sinica (Magnoliaceae) provides insights into the conservation of plant species with extremely small populations (PSESP)

Magnolia sinica (Magnoliaceae) is a highly threatened tree endemic to southeast Yunnan, China. In this study, we generated for the first time a high-quality chromosome-scale genome sequence from M. sinica, by combining Illumina and ONT data with Hi-C scaffolding methods. The final assembled genome size of M. sinica was 1.84 Gb, with a contig N50 of ca. 45 Mb and scaffold N50 of 92 Mb. Identified repeats constituted approximately 57% of the genome, and 43,473 protein-coding genes were predicted. Phylogenetic analysis shows that the magnolias form a sister clade with the eudicots and the order Ceratophyllales, while the monocots are sister to the other core angiosperms. In our study, a total of 21 individuals from the 5 remnant populations of M. sinica, as well as 22 specimens belonging to 8 related Magnoliaceae species, were resequenced. The results showed that M. sinica had higher genetic diversity (θw = 0.01126 and θπ = 0.01158) than other related species in the Magnoliaceae. However, population structure analysis suggested that the genetic differentiation among the 5 M. sinica populations was very low. Analyses of the demographic history of the species using different models consistently revealed that 2 bottleneck events occurred. The contemporary effective population size of M. sinica was estimated to be 10.9. The different patterns of genetic loads (inbreeding and numbers of deleterious mutations) suggested constructive strategies for the conservation of these 5 different populations of M. sinica. Overall, this high-quality genome will be a valuable genomic resource for conservation of M. sinica.

Genome sequences and population genomics provide insights into the demographic history, inbreeding, and mutation load of two 'living fossil' tree species of Dipteronia

'Living fossils', that is, ancient lineages of low taxonomic diversity, represent an exceptional evolutionary heritage, yet we know little about how demographic history and deleterious mutation load have affected their long-term survival and extinction risk. We performed whole-genome sequencing and population genomic analyses on Dipteronia sinensis and D. dyeriana, two East Asian Tertiary relict trees. We found large-scale genome reorganizations and identified species-specific genes under positive selection that are likely involved in adaptation. Our demographic analyses suggest that the wider-ranged D. sinensis repeatedly recovered from population bottlenecks over late Tertiary/Quaternary periods of adverse climate conditions, while the population size of the narrow-ranged D. dyeriana steadily decreased since the late Miocene, especially after the Last Glacial Maximum (LGM). We conclude that the efficient purging of deleterious mutations in D. sinensis facilitated its survival and repeated demographic recovery. By contrast, in D. dyeriana, increased genetic drift and reduced selection efficacy, due to recent severe population bottlenecks and a likely preponderance of vegetative propagation, resulted in fixation of strongly deleterious mutations, reduced fitness, and continuous population decline, with likely detrimental consequences for the species' future viability and adaptive potential. Overall, our findings highlight the significant impact of demographic history on levels of accumulation and purging of putatively deleterious mutations that likely determine the long-term survival and extinction risk of Tertiary relict trees.

Two chromosome-level genome assemblies of Rhodiola shed new light on genome evolution in rapid radiation and evolution of the biosynthetic pathway of salidroside

Rhodiola L. is a genus that has undergone rapid radiation in the mid-Miocene and may represent a typic case of adaptive radiation. Many species of Rhodiola have also been widely used as an important adaptogen in traditional medicines for centuries. However, a lack of high-quality chromosome-level genomes hinders in-depth study of its evolution and biosynthetic pathway of secondary metabolites. Here, we assembled two chromosome-level genomes for two Rhodiola species with different chromosome number and sexual system. The assembled genome size of R. chrysanthemifolia (2n = 14; hermaphrodite) and R. kirilowii (2n = 22; dioecious) were of 402.67 and 653.62 Mb, respectively, with approximately 57.60% and 69.22% of transposable elements (TEs). The size difference between the two genomes was mostly due to proliferation of long terminal repeat-retrotransposons (LTR-RTs) in the R. kirilowii genome. Comparative genomic analysis revealed possible gene families responsible for high-altitude adaptation of Rhodiola, including a homolog of plant cysteine oxidase 2 gene of Arabidopsis thaliana (AtPCO2), which is part of the core molecular reaction to hypoxia and contributes to the stability of Group VII ethylene response factors (ERF-VII). We found extensive chromosome fusion/fission events and structural variations between the two genomes, which might have facilitated the initial rapid radiation of Rhodiola. We also identified candidate genes in the biosynthetic pathway of salidroside. Overall, our results provide important insights into genome evolution in plant rapid radiations, and possible roles of chromosome fusion/fission and structure variation played in rapid speciation.

Phylogenomic and syntenic data demonstrate complex evolutionary processes in early radiation of the rosids

Some of the most vexing problems of deep level relationship that remain in angiosperms involve the superrosids. The superrosid clade contains a quarter of all angiosperm species, with 18 orders in three subclades (Vitales, Saxifragales and core rosids) exhibiting remarkable morphological and ecological diversity. To help resolve deep-level relationships, we constructed a high-quality chromosome-level genome assembly for Tiarella polyphylla (Saxifragaceae) thus providing broader genomic representation of Saxifragales. Whole genome microsynteny analysis of superrosids showed that Saxifragales shared more synteny clusters with core rosids than Vitales, further supporting Saxifragales as more closely related with core rosids. To resolve the ordinal phylogeny of superrosids, we screened 122 single copy nuclear genes from genomes of 36 species, representing all 18 superrosid orders. Vitales were recovered as sister to all other superrosids (Saxifragales + core rosids). Our data suggest dramatic differences in relationships compared to earlier studies within core rosids. Fabids should be restricted to the nitrogen-fixing clade, while Picramniales, the Celastrales-Malpighiales (CM) clade, Huerteales, Oxalidales, Sapindales, Malvales and Brassicales formed an "expanded" malvid clade. The Celastrales-Oxalidales-Malpighiales (COM) clade (sensu APG IV) was not monophyletic. Crossosomatales, Geraniales, Myrtales and Zygophyllales did not belong to either of our well-supported malvids or fabids. There is strong discordance between nuclear and plastid phylogenetic hypotheses for superrosid relationships; we show that this is best explained by a combination of incomplete lineage sorting and ancient reticulation.

Karyotype evolution of the Asterids insights from the first genome sequences of the family Cornaceae

Cornaceae is a core representative family in Cornales, the earliest branching lineage in the Asterids on the life tree of angiosperms. This family includes the only genus Cornus, a group of ~55 species. These species occur widely in Northern Hemisphere and have been used as resources for horticultural ornaments, medicinal and industrial manufacturing. However, no any genome sequences are available for this family. Here, we reported a chromosome-level genome for Cornus controversa. This was generated using high-fidelity plus Hi-C sequencing, and totally ~771.80 Mb assembled sequences and 39,886 protein-coding genes were obtained. We provided evidence for a whole-genome duplication event (WGD) unique to C. controversa. The evolutionary features of this genome indicated that the expanded and unique genes might have contributed to response to stress, stimulus and defense. By using chromosome-level syntenic blocks shared between eight living genomes, we found high degrees of genomic diversification from the ancestral core-eudicot genome to the present-day genomes, suggesting an important role of WGD in genomic plasticity that leads to speciation and diversification. These results provide foundational insights on the evolutionary history of Cornaceae, as well as on the Asterids diversification.

Comparative chloroplast genomics provides insights into the genealogical relationships of endangered Tetraena mongolica and the chloroplast genome evolution of related Zygophyllaceae species

A comprehensive understanding of genetic background for rare species will provide an important theoretical basis for the future species management, monitoring and conservation. Tetraena mongolica is restrictedly distributed in the western Ordos plateau of China and has been listed as a national protected plant. We generated 13 chloroplast (cp) genomes of T. mongolica (size range of 106,062-106,230 bp) and conducted a series of comparative analyses of six Zygophyllaceae cp genomes. T. mongolica cp genome exhibited a quadripartite structure with drastically reduced inverted repeats (IRs, 4,315 bp) and undergone the loss of a suit of ndh genes and a copy of rRNAs. Furthermore, all the T. mongolica populations were divided into two genetic groups based on complete cp phylogenomics. In addition, notably variable genome size, gene order and structural changes had been observed among the six Zygophyllaceae cp genomes. Overall, our findings provide insights into the cp genome evolution mode and intraspecific relationships of T. mongolica, and provide a molecular basis for scientific conservation of this endangered plant.

On the hidden temporal dynamics of plant adaptation

Adaptation to a wide range of environments is a major driver of plant diversity. It is now possible to catalog millions of potential adaptive genomic differences segregating between environments within a plant species in a single experiment. Understanding which of these changes contributes to adaptive phenotypic divergence between plant populations is a major goal of evolutionary biologists and crop breeders. In this review, we briefly highlight the approaches frequently used to understand the genetic basis of adaptive phenotypes in plants, and we discuss some of the limitations of these methods. We propose that direct observation of the process of adaptation using multigenerational studies and whole genome sequencing is a crucial missing component of recent studies of plant adaptation because it complements several shortcomings of sampling-based techniques.

Humoral response to SARS-CoV-2 and seasonal coronaviruses in COVID-19 patients

We used enzyme-linked immunoassay methods to measure the prevalence and the levels of antibody responses to the nucleocapsid (N) protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and four seasonal human coronaviruses (HCoV-OC43, HCoV-HKU1, HCoV 229E, and HCoV-NL63) in a cohort of 115 convalescent plasma donors infected with SARS-CoV-2 (1-61 days after symptom onset) compared to antibody levels in 114 individuals with no evidence of a recent infection with SARS-CoV-2. In the humoral response to the four seasonal coronaviruses, only HCoV-HKU1- and HCoV-229E-assays showed slightly elevated antibody levels in the COVID group compared to the control group. While in the COVID-group the levels of SARS-CoV-2 antibodies correlated significantly with disease severity, no association was found in the levels of antibodies against the seasonal coronaviruses. The most striking result in both groups was that the levels of antibodies against all tested coronaviruses, including the new SARS-CoV-2 showed a highly significant correlation with each other. There seems to be an individual predisposition to a weaker or stronger humoral immune response against all known seasonal human coronaviruses including the new SARS-CoV-2, which could lead to a definition of low and high responders against human coronaviruses with potential impact on the assessment of postinfection antibody levels and protection.

The genome of the Paleogene relic tree Bretschneidera sinensis: insights into trade-offs in gene family evolution, demographic history, and adaptive SNPs

Among relic species, genomic information may provide the key to inferring their long-term survival. Therefore, in this study, we investigated the genome of the Paleogene relic tree species, Bretschneidera sinensis, which is a rare endemic species within southeastern Asia. Specifically, we assembled a high-quality genome for B. sinensis using PacBio high-fidelity and high-throughput chromosome conformation capture reads and annotated it with long and short RNA sequencing reads. Using the genome, we then detected a trade-off between active and passive disease defences among the gene families. Gene families involved in salicylic acid and MAPK signalling pathways expanded as active defence mechanisms against disease, but families involved in terpene synthase activity as passive defences contracted. When inferring the long evolutionary history of B. sinensis, we detected population declines corresponding to historical climate change around the Eocene–Oligocene transition and to climatic fluctuations in the Quaternary. Additionally, based on this genome, we identified 388 single nucleotide polymorphisms (SNPs) that were likely under selection, and showed diverse functions in growth and stress responses. Among them, we further found 41 climate-associated SNPs. The genome of B. sinensis and the SNP dataset will be important resources for understanding extinction/diversification processes using comparative genomics in different lineages.

Landscape Genomics in Tree Conservation Under a Changing Environment

Understanding the genetic basis of how species respond to changing environments is essential to the conservation of species. However, the molecular mechanisms of adaptation remain largely unknown for long-lived tree species which always have large population sizes, long generation time, and extensive gene flow. Recent advances in landscape genomics can reveal the signals of adaptive selection linking genetic variations and landscape characteristics and therefore have created novel insights into tree conservation strategies. In this review article, we first summarized the methods of landscape genomics used in tree conservation and elucidated the advantages and disadvantages of these methods. We then highlighted the newly developed method "Risk of Non-adaptedness," which can predict the genetic offset or genomic vulnerability of species via allele frequency change under multiple scenarios of climate change. Finally, we provided prospects concerning how our introduced approaches of landscape genomics can assist policymaking and improve the existing conservation strategies for tree species under the ongoing global changes.

Species Delimitation of Asteropyrum (Ranunculaceae) Based on Morphological, Molecular, and Ecological Variation

Objectively evaluating different lines of evidence within a formalized framework is the most efficient and theoretically grounded approach for defining robust species hypotheses. Asteropyrum Drumm. et Hutch. is a small genus of perennial herb containing two species, A. cavaleriei and A. peltatum. The distinction of these two species mainly lies in the shape and size of leaf blades. However, these characters have been considered labile and could not differentiate the two species reliably. In this study, we investigated the variation of the leaf blades of 28 populations across the whole range of Asteropyrum using the landmark-based geometric morphometrics (GMM), sought genetic gaps within this genus using DNA barcoding, phylogenetic reconstruction and population genetic methods, and compared the predicted ecological niches of the two species. The results showed that the leaf form (shape and size) was overlapped between the two species; barcode gap was not detected within the genus Asteropyrum; and little ecological and geographical differentiation was found between the two taxa. Two genetic clusters detected by population genetic analysis did not match the two morphospecies. The results suggest that there are no distinct boundaries between the two species of Asteropyrum in terms of morphology, genetics and ecology and this present classification should be abandoned. We anticipate that range-wide population genomic studies would properly delineate the species boundaries and help to understand the evolution and speciation within Asteropyrum.

The global significance of biodiversity science in China: an overview

Biodiversity science in China has seen rapid growth over recent decades, ranging from baseline biodiversity studies to understanding the processes behind evolution across dynamic regions such as the Qinghai-Tibetan Plateau. We review research, including species catalogues; biodiversity monitoring; the origins, distributions, maintenance and threats to biodiversity; biodiversity-related ecosystem function and services; and species and ecosystems' responses to global change. Next, we identify priority topics and offer suggestions and priorities for future biodiversity research in China. These priorities include (i) the ecology and biogeography of the Qinghai-Tibetan Plateau and surrounding mountains, and that of subtropical and tropical forests across China; (ii) marine and inland aquatic biodiversity; and (iii) effective conservation and management to identify and maintain synergies between biodiversity and socio-economic development to fulfil China's vision for becoming an ecological civilization. In addition, we propose three future strategies: (i) translate advanced biodiversity science into practice for biodiversity conservation; (ii) strengthen capacity building and application of advanced technologies, including high-throughput sequencing, genomics and remote sensing; and (iii) strengthen and expand international collaborations. Based on the recent rapid progress of biodiversity research, China is well positioned to become a global leader in biodiversity research in the near future.