Opening a door to the spatiotemporal history of plants from the tropical Indochina Peninsula to subtropical China

Plastome evolution of Engelhardia facilitates phylogeny of Juglandaceae

BACKGROUND

Engelhardia (Juglandaceae) is a genus of significant ecological and economic importance, prevalent in the tropics and subtropics of East Asia. Although previous efforts based on multiple molecular markers providing profound insights into species delimitation and phylogeography of Engelhardia, the maternal genome evolution and phylogeny of Engelhardia in Juglandaceae still need to be comprehensively evaluated. In this study, we sequenced plastomes from 14 samples of eight Engelhardia species and the outgroup Rhoiptelea chiliantha, and incorporated published data from 36 Juglandaceae and six outgroup species to test phylogenetic resolution. Moreover, comparative analyses of the plastomes were conducted to investigate the plastomes evolution of Engelhardia and the whole Juglandaceae family.

RESULTS

The 13 Engelhardia plastomes were highly similar in genome size, gene content, and order. They exhibited a typical quadripartite structure, with lengths from 161,069 bp to 162,336 bp. Three mutation hotspot regions (TrnK-rps16, ndhF-rpl32, and ycf1) could be used as effective molecular markers for further phylogenetic analyses and species identification. Insertion and deletion (InDels) may be an important driving factor for the evolution of plastomes in Juglandoideae and Engelhardioideae. A total of ten codons were identified as the optimal codons in Juglandaceae. The mutation pressure mostly contributed to shaping codon usage. Seventy-eight protein-coding genes in Juglandaceae experienced relaxed purifying selection, only rpl22 and psaI genes showed positive selection (Ka/Ks > 1). Phylogenetic results fully supported Engelhardia as a monophyletic group including two sects and the division of Juglandaceae into three subfamilies. The Engelhardia originated in the Late Cretaceous and diversified in the Late Eocene, and Juglandaceae originated in the Early Cretaceous and differentiated in Middle Cretaceous. The phylogeny and divergence times didn't support rapid radiation occurred in the evolution history of Engelhardia.

CONCLUSION

Our study fully supported the taxonomic treatment of at the section for Engelhardia species and three subfamilies for Juglandaceae and confirmed the power of phylogenetic resolution using plastome sequences. Moreover, our results also laid the foundation for further studying the course, tempo and mode of plastome evolution of Engelhardia and the whole Juglandaceae family.

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.

Distribution and Population Genetic Structure of the Hau Giang Medaka, Oryzias haugiangensis, Along the East Coast of the Indochinese Peninsula

The east coast of the Indochinese Peninsula is a well-known transition zone from subtropical to tropical systems, yet only a small number of studies have been conducted on the biogeography and phylogeography of aquatic organisms in this region. The Hau Giang medaka, Oryzias haugiangensis, was originally described from the Mekong Delta in southern Vietnam, and later reported also from southeastern Thailand, west of the Mekong Delta region. However, the species' full geographic range and population genetic structures remain unknown. Field surveys showed a widespread distribution of this species along the east coast of the Indochinese Peninsula, as far as northern Vietnam. A mitochondrial gene phylogeny and population genetic structure analysis using genome-wide single nucleotide polymorphisms revealed that the populations of O. haugiangensis are highly structuralized along the east coast of Vietnam, with the southernmost Mekong Delta population clearly separated from three populations north of central Vietnam. Further field collections are necessary to determine the boundary between the southern and northern populations, and the presence or absence of a hybrid zone.

Adaptive divergence and genetic vulnerability of relict species under climate change: a case study of Pterocarya macroptera

BACKGROUND AND AIMS

Understanding adaptive genetic variation and whether it can keep pace with predicted future climate change is critical in assessing the genetic vulnerability of species and developing conservation management strategies. The lack of information on adaptive genetic variation in relict species carrying abundant genetic resources hinders the assessment of genetic vulnerability. Using a landscape genomics approach, this study aimed to determine how adaptive genetic variation shapes population divergence and to predict the adaptive potential of Pterocarya macroptera (a vulnerable relict species in China) under future climate scenarios.

METHODS

We applied restriction site-associated DNA sequencing (RAD-seq) to obtain 8244 single-nucleotide polymorphisms (SNPs) from 160 individuals across 28 populations. We examined the pattern of genetic diversity and divergence, and then identified outliers by genetic differentiation (FST) and genotype-environment association (GEA) methods. We further dissected the effect of geographical/environmental gradients on genetic variation. Finally, we predicted genetic vulnerability and adaptive risk under future climate scenarios.

KEY RESULTS

We identified three genetic lineages within P. macroptera: the Qinling-Daba-Tianmu Mountains (QDT), Western Sichuan (WS) and Northwest Yunnan (NWY) lineages, which showed significant signals of isolation by distance (IBD) and isolation by environment (IBE). IBD and IBE explained 3.7-5.7 and 8.6-12.8 % of the genetic structure, respectively. The identified GEA SNP-related genes were involved in chemical defence and gene regulation and may exhibit higher genetic variation to adapt to the environment. Gradient forest analysis revealed that the genetic variation was mainly shaped by temperature-related variables, indicating its adaptation to local thermal environments. A limited adaptive potential was suggested by the high levels of genetic vulnerability in marginal populations.

CONCLUSIONS

Environmental gradient mainly shaped the population differentiation of P. macroptera. Marginal populations may be at high risk of extinction, and thus proactive management measures, such as assisted gene flow, are required to ensure the survival of these populations.

Genome structure-based Juglandaceae phylogenies contradict alignment-based phylogenies and substitution rates vary with DNA repair genes

In lineages of allopolyploid origin, sets of homoeologous chromosomes may coexist that differ in gene content and syntenic structure. Presence or absence of genes and microsynteny along chromosomal blocks can serve to differentiate subgenomes and to infer phylogenies. We here apply genome-structural data to infer relationships in an ancient allopolyploid lineage, the walnut family (Juglandaceae), by using seven chromosome-level genomes, two of them newly assembled. Microsynteny and gene-content analyses yield identical topologies that place Platycarya with Engelhardia as did a 1980s morphological-cladistic study. DNA-alignment-based topologies here and in numerous earlier studies instead group Platycarya with Carya and Juglans, perhaps misled by past hybridization. All available data support a hybrid origin of Juglandaceae from extinct or unsampled progenitors nested within, or sister to, Myricaceae. Rhoiptelea chiliantha, sister to all other Juglandaceae, contains proportionally more DNA repair genes and appears to evolve at a rate 2.6- to 3.5-times slower than the remaining species.

Monsoon intensification in East Asia triggered the evolution of its flora

INTRODUCTION

East Asia (EA), which falls within the region of the Asian monsoon that is composed of the East Asia monsoon (EAM) and the Indian monsoon (IM), is known for its high species diversity and endemism. This has been attributed to extreme physiographical heterogeneity in conjunction with climate and sea-level changes during the Pleistocene, this hypothesis has been widely proven by phylogeographic studies. Recently, dated phylogenies have indicated that the origins (stem age) of the flora occurred after the Oligocene-Miocene boundary and are related to the establishment of the EAM.

METHODS

Hence, this study further examined whether the strengthening of the monsoons triggered floral evolution via a meta-analysis of the tempo-spatial pattern of evolutionary radiation dates (crown ages) of 101 endemic seed plant genera.

RESULTS

Taxonomic diversification began during the late Eocene, whereas the accumulated number of diversifications did not significantly accelerate until the late Miocene. The distribution of the weighted mean and the average divergence times in the EAM, IM, or transitional regions all fall within the mid-late Miocene. Fossils of the Tertiary relict genera are mostly and widely distributed outside EA and only half of the earliest fossils in the EA region are not older than Miocene, while their divergence times are mostly after the late Miocene. The pattern of divergence time of monotypic and polytypic taxa suggest the climatic changes after the late Pliocene exert more influence on monotypic taxa.

DISCUSSION

The two key stages of floral evolution coincide with the intensifications of the EAM and IM, especially the summer monsoon which brings a humid climate. An integrated review of previous studies concerning flora, genus, and species levels further supports our suggestion that monsoon intensification in EA triggered the evolution of its flora.

Latitudinal Diversity Gradient in the Changing World: Retrospectives and Perspectives

The latitudinal diversity gradient (LDG) is one of the most extensive and important biodiversity patterns on the Earth. Various studies have established that species diversity increases with higher taxa numbers from the polar to the tropics. Studies of multicellular biotas have supported the LDG patterns from land (e.g., plants, animals, forests, wetlands, grasslands, fungi, and so forth) to oceans (e.g., marine organisms from freshwater invertebrates, continental shelve, open ocean, even to the deep sea invertebrates). So far, there are several hypotheses proposed to explore the diversity patterns and mechanisms of LDG, however, there has been no consensus on the underlying causes of LDG over the past few decades. Thus, we reviewed the progress of LDG studies in recent years. Although several explanations for the LDG have been proposed, these hypotheses are only based on species richness, evolution and the ecosystems. In this review, we summarize the effects of evolution and ecology on the LDG patterns to synthesize the formation mechanisms of the general biodiversity distribution patterns. These intertwined factors from ecology and evolution in the LDG are generally due to the wider distribution of tropical areas, which hinders efforts to distinguish their relative contributions. However, the mechanisms of LDG always engaged controversies, especially in such a context that the human activity and climate change has affected the biodiversity. With the development of molecular biology, more genetic/genomic data are available to facilitate the estimation of global biodiversity patterns with regard to climate, latitude, and other factors. Given that human activity and climate change have inevitably impacted on biodiversity loss, biodiversity conservation should focus on the change in LDG pattern. Using large-scale genetic/genomic data to disentangle the diversity mechanisms and patterns of LDG, will provide insights into biodiversity conservation and management measures. Future perspectives of LDG with integrative genetic/genomic, species, evolution, and ecosystem diversity patterns, as well as the mechanisms that apply to biodiversity conservation, are discussed. It is imperative to explore integrated approaches for recognizing the causes of LDG in the context of rapid loss of diversity in a changing world.