Demographically idiosyncratic responses to climate change and rapid Pleistocene diversification of the walnut genus Juglans (Juglandaceae) revealed by whole-genome sequences

Pleistocene Glaciation Drove Shared Population Coexpansion in Eastern North American Snakes

Glacial cycles during the Pleistocene had profound impacts on local environments and climatic conditions. In North America, some regions that currently support diverse biomes were entirely covered by ice sheets, while other regions were environmentally unsuitable for the organisms that live there now. Organisms that occupy these regions in the present day must have expanded or dispersed into these regions since the last glacial maximum, leading to the possibility that species with similar geographic distributions may show temporally concordant population size changes associated with these warming trends. We examined 17 lineages from 9 eastern North American snake species and species complexes to test for a signal of temporally concordant coexpansion using a machine learning approach. We found that the majority of lineages show population size increases towards the present, with evidence for coexpansion in five out of fourteen lineages, while expansion in others was idiosyncratic. We also examined relationships between genetic distance and current environmental predictors and showed that genomic responses to environmental predictors are not consistent among species. We, therefore, conclude that Pleistocene warming resulted in population size increases in most eastern North American snake species, but variation in environmental preferences and other species-specific traits results in variance in the exact timing of expansion.

'Chimes of resilience': what makes forest trees genetically resilient?

Forest trees are foundation species of many ecosystems and are challenged by global environmental changes. We assemble genetic facts and arguments supporting or undermining resilient responses of forest trees to those changes. Genetic resilience is understood here as the capacity of a species to restore its adaptive potential following environmental changes and disturbances. Importantly, the data come primarily from European temperate tree species with large distributions and consider only marginally species with small distributions. We first examine historical trajectories of trees during repeated climatic changes. Species that survived the Pliocene-Pleistocene transition and underwent the oscillations of glacial and interglacial periods were equipped with life history traits enhancing persistence and resilience. Evidence of their resilience also comes from the maintenance of large effective population sizes across time and rapid microevolutionary responses to recent climatic events. We then review genetic mechanisms and attributes shaping resilient responses. Usually, invoked constraints to resilience, such as genetic load or generation time and overlap, have limited consequences or are offset by positive impacts. Conversely, genetic plasticity, gene flow, introgression, genetic architecture of fitness-related traits and demographic dynamics strengthen resilience by accelerating adaptive responses. Finally, we address the limitations of this review and highlight critical research gaps.

Low genetic diversity and weak population structure of Albizia odoratissima on Hainan Island

BACKGROUND

The increasing demand for wood owing to societal development has highlighted the potential of Albizia odoratissima, a valuable timber species, to address significant timber shortages in China. However, the lack of effective genetic and genomic resources has limited the development and utilization of this species.

RESULTS

In this study, we utilised 95.3 Gb of HiFi reads to assemble a draft genome of A. odoratissima, resulting in a genome size of 788 Mb, comprising 511 contigs. We conducted whole-genome resequencing on 106 individuals from 7 populations on Hainan Island to explore these resources. Our analysis identified 498,308 high-quality single nucleotide polymorphisms, which were used to assess the genetic diversity, structure, and demographic history of A. odoratissima on Hainan Island. The results indicated that the genetic diversity of A. odoratissima on Hainan Island is relatively low (observed heterozygosity = 0.189, expected heterozygosity = 0.189, genetic diversity = 1.319 × 10-4) with minimal genetic differentiation (Fst = 0.0151) among the seven populations. Furthermore, molecular variance, principal coordinate analysis, neighbour-joining tree analysis, and genetic structure analysis revealed a shallow population structure. The linkage disequilibrium (LD) decay ranged from 11.4 kb for Jianfengling (JFL) to 39.2 kb for Wuzhishan (WZS). LD decay, demographic history, and Tajima's D analyses indicated that the WZS population has experienced a bottleneck effect.

CONCLUSIONS

This study offers new insights into the genetic diversity and population structure of A. odoratissima on Hainan Island, providing a foundation for future resource utilization and genetic improvement strategies for this species.

Tracing the path from conservation to expansion evolutionary insights into NLR genes in oleaceae

The Oleaceae family, encompassing key genera such as Fraxinus (ash trees), Olea (olives), Jasminum (jasmine), Syringa (lilac), and Forsythia, plays a crucial ecological and economic role. Despite their importance, the evolutionary dynamics and immune system adaptations of their NLR (Nucleotide binding leucine-rich repeats) gene family remain largely unexplored. This study employs high-throughput comparative genomics to investigate NLR gene evolution across the Oleaceae family. The genus Fraxinus is widely distributed across both the New and Old Worlds, with 23 distinct species analyzed in this study. Our results reveal a predominant strategy of gene conservation in the evolution of the NLR gene family across these species. Geographical adaptation has played a significant role, particularly in Old World ash tree species, which exhibit dynamic patterns of gene expansion and contraction within the last 50 million years. Notably, genes acquired from an ancient whole genome duplication event (~ 35 Mya) have been retained across Fraxinus lineages. In contrast, the genus Olea (olives) has undergone extensive gene expansion driven by recent duplications and significant birth of novel NLR gene families. These differences in NLR gene evolution likely enhance Olea's ability to recognize diverse pathogens through recent expansions, while Fraxinus maintains specialized immune responses through conserved genes, with potential trade-offs in pathogen adaptation and energy efficiency. In terms of NLR distribution, all species of the Oleaceae family show an enhanced pseudogenization of TIR-NLRs and expansion in CCG10-NLR. However, the comparative RNA-seq expression analysis in olive suggests that partial NLR genes, despite their incomplete structure, have significant expression and may play important roles in plant immune responses. This study provides a comprehensive analysis of NLR gene evolution within the Oleaceae family, offering insights into the adaptive mechanisms of immune response evolution across diverse genera.

Genome Assembly and Winged Fruit Gene Regulation of Chinese Wingnut: Insights from Genomic and Transcriptomic Analyses

The genomic basis and biology of winged fruit are interesting issues in ecological and evolutionary biology. Chinese wingnut (Pterocarya stenoptera) is an important horticultural and economic tree species in China. The genomic resources of this hardwood tree could advance the genomic studies of Juglandaceae species and elucidate their evolutionary relationships. Here, we reported a high-quality reference genome of P. stenoptera (N50 = 35.15 Mb) and performed a comparative genomic analysis across Juglandaceae species. Paralogous relationships among the 16 chromosomes of P. stenoptera revealed eight main duplications representing the subgenomes. Molecular dating suggested that the most recent common ancestor of P. stenoptera and Cyclocarya paliurus diverged from Juglans species around 56.7 million years ago (MYA). The expanded and contracted gene families were associated with cutin, suberine, and wax biosynthesis, cytochrome P450, and anthocyanin biosynthesis. We identified large inversion blocks between P. stenoptera and its relatives, which were enriched with genes involved in lipid biosynthesis and metabolism, as well as starch and sucrose metabolism. Whole-genome resequencing of 28 individuals revealed clearly phylogenetic clustering into three groups corresponding to Pterocarya macroptera, Pterocarya hupehensis, and P. stenoptera. Morphological and transcriptomic analyses showed that CAD, COMT, LOX, and MADS-box play important roles during the five developmental stages of wingnuts. This study highlights the evolutionary history of the P. stenoptera genome and supports P. stenoptera as an appropriate Juglandaceae model for studying winged fruits. Our findings provide a theoretical basis for understanding the evolution, development, and diversity of winged fruits in woody plants.

Species-Specific Traits Shape Genetic Diversity During an Expansion-Contraction Cycle and Bias Demographic History Reconstruction

Species ranges are dynamic, experiencing expansions, contractions or shifts in response to habitat changes driven by extrinsic factors such as climate change or human activities. While existing research examines the genetic consequences of spatial processes, few studies integrate species-specific traits to analyse how habitat changes affect co-existing species. In this study, we address this gap by investigating how genetic diversity patterns vary among species with different traits (such as generation length, population density and dispersal) experiencing similar habitat changes. Using spatial simulations and a simpler panmictic population model, we investigate the temporal genetic diversity in refugium populations undergoing range expansion of their habitat, followed by stationary and contraction periods. By varying habitat contraction speed and species traits, we identified three distinct temporal dynamics of genetic diversity during contraction: (i) a decrease in genetic diversity, (ii) an initial increase followed by a decrease and (iii) a continuous increase throughout the contraction period. We show that genetic diversity trajectories during population decline can be predicted by comparing sampled population diversity to equilibrium values expected under expanded and contracted habitat ranges. Our study also challenges the belief that high genetic diversity in a refugium population is due to a recent and rapid habitat loss. Instead, we found contrasting effects of contraction speed on genetic diversity depending on the interaction between species-specific traits and the dynamics of habitat change. Finally, using simulated genetic data, we found that demographic histories inferred from effective population size estimates may vary across species, even when they experience similar habitat changes.

Resilience of genetic diversity in forest trees over the Quaternary

The effect of past environmental changes on the demography and genetic diversity of natural populations remains a contentious issue and has rarely been investigated across multiple, phylogenetically distant species. Here, we perform comparative population genomic analyses and demographic inferences for seven widely distributed and ecologically contrasting European forest tree species based on concerted sampling of 164 populations across their natural ranges. For all seven species, the effective population size, Ne, increased or remained stable over many glacial cycles and up to 15 million years in the most extreme cases. Surprisingly, the drastic environmental changes associated with the Pleistocene glacial cycles have had little impact on the level of genetic diversity of dominant forest tree species, despite major shifts in their geographic ranges. Based on their trajectories of Ne over time, the seven tree species can be divided into three major groups, highlighting the importance of life history and range size in determining synchronous variation in genetic diversity over time. Altogether, our results indicate that forest trees have been able to retain their evolutionary potential over very long periods of time despite strong environmental changes.

The 'queen of the Andes' (Puya raimondii) is genetically fragile and fragmented: a consequence of long generation time and semelparity?

Understanding how life history shapes genetic diversity is a fundamental issue in evolutionary biology, with important consequences for conservation. However, we still have an incomplete picture of the impact of life history on genome-wide patterns of diversity, especially in long-lived semelparous plants. Puya raimondii is a high-altitude semelparous species from the Andes that flowers at 40-100 years of age. We sequenced the whole genome and estimated the nucleotide diversity of 200 individuals sampled from nine populations. Coalescent-based approaches were then used to infer past population dynamics. Finally, these results were compared with results obtained for the iteroparous species, Puya macrura. The nine populations of P. raimondii were highly divergent, highly inbred, and carried an exceptionally high genetic load. They are genetically depauperate, although, locally in the genome, balancing selection contributed to the maintenance of genetic polymorphism. While both P. raimondii and P. macrura went through a severe bottleneck during the Pleistocene, P. raimondii did not recover from it and continuously declined, while P. macrura managed to bounce back. Our results demonstrate the importance of life history, in particular generation time and reproductive strategy, in affecting population dynamics and genomic variation, and illustrate the genetic fragility of long-lived semelparous plants.

Complete mitochondrial genome assembly of Juglans regia unveiled its molecular characteristics, genome evolution, and phylogenetic implications

BACKGROUND

The Persian walnut (Juglans regia), an economically vital species within the Juglandaceae family, has seen its mitochondrial genome sequenced and assembled in the current study using advanced Illumina and Nanopore sequencing technology.

RESULTS

The 1,007,576 bp mitogenome of J. regia consisted of three circular chromosomes with a 44.52% GC content encoding 39 PCGs, 47 tRNA, and five rRNA genes. Extensive repetitive sequences, including 320 SSRs, 512 interspersed, and 83 tandem repeats, were identified, contributing to genomic complexity. The protein-coding sequences (PCGs) favored A/T-ending codons, and the codon usage bias was primarily shaped by selective pressure. Intracellular gene transfer occurred among the mitogenome, chloroplast, and nuclear genomes. Comparative genomic analysis unveiled abundant structure and sequence variation among J. regia and related species. The results of selective pressure analysis indicated that most PCGs underwent purifying selection, whereas the atp4 and ccmB genes had experienced positive selection between many species pairs. In addition, the phylogenetic examination, grounded in mitochondrial genome data, precisely delineated the evolutionary and taxonomic relationships of J. regia and its relatives. We identified a total of 539 RNA editing sites, among which 288 were corroborated by transcriptome sequencing data. Furthermore, expression profiling under temperature stress highlighted the complex regulation pattern of 28 differently expressed PCGs, wherein NADH dehydrogenase and ATP synthase genes might be critical in the mitochondria response to cold stress.

CONCLUSIONS

Our results provided valuable molecular resources for understanding the genetic characteristics of J. regia and offered novel perspectives for population genetics and evolutionary studies in Juglans and related woody species.

Genomic divergence and mutation load in the Begonia masoniana complex from limestone karsts

Understanding genome-wide diversity, inbreeding, and the burden of accumulated deleterious mutations in small and isolated populations is essential for predicting and enhancing population persistence and resilience. However, these effects are rarely studied in limestone karst plants. Here, we re-sequenced the nuclear genomes of 62 individuals of the Begonia masoniana complex (B. liuyanii, B. longgangensis, B. masoniana and B. variegata) and investigated genomic divergence and genetic load for these four species. Our analyses revealed four distinct clusters corresponding to each species within the complex. Notably, there was only limited admixture between B. liuyanii and B. longgangensis occurring in overlapping geographic regions. All species experienced historical bottlenecks during the Pleistocene, which were likely caused by glacial climate fluctuations. We detected an asymmetric historical gene flow between group pairs within this timeframe, highlighting a distinctive pattern of interspecific divergence attributable to karst geographic isolation. We found that isolated populations of B. masoniana have limited gene flow, the smallest recent population size, the highest inbreeding coefficients, and the greatest accumulation of recessive deleterious mutations. These findings underscore the urgency to prioritize conservation efforts for these isolated population. This study is among the first to disentangle the genetic differentiation and specific demographic history of karst Begonia plants at the whole-genome level, shedding light on the potential risks associated with the accumulation of deleterious mutations over generations of inbreeding. Moreover, our findings may facilitate conservation planning by providing critical baseline genetic data and a better understanding of the historical events that have shaped current population structure of rare and endangered karst plants.

Historical spread routes of wild walnuts in Central Asia shaped by man-made and nature

Walnuts have substantial economic value and are of significant interest being a wild-cultivated species. The study has re-sequenced the entire genome of the wild walnut, aligning it with the walnut reference genome, to identify 2,021,717 single nucleotide polymorphisms (SNPs). These were used to examine the genetics of 130 wild walnut samples collected from three countries. Utilizing structural and principal component analysis, the walnut samples from Central Asia were classified into four populations: Ili ah in Xinjiang (I), Dushanbe region in Tajikistan (II), Sary-Chelek, Arslanbob in Kara-Alma regions of Kyrgyzstan (III), and Kok-Tundy region of Kyrgyzstan (IV). The 4 groups showed large differences in nucleotide diversity, population differentiation, and linkage disequilibrium decay, as well as gene flow among them. The present geographic distribution of these populations does not align with the genetic distribution pattern as the populations of Central Asian wild walnuts have experienced similar population dynamics in the past, i.e., the highest effective population size at ca. 6 Ma, two sharp population declines at 6 and 0.2 Ma, and convergence at ca. 0.2 Ma. The genetic distribution patterns are better explained by human activity, notably through archaeological findings of walnut use and the influence of the Silk Road, rather than by current geographic distributions.

Genetic diversity and selection signatures in Hainan black goats revealed by whole-genome sequencing data

Understanding the genetic characteristics of indigenous goat breeds is crucial for their conservation and breeding efforts. Hainan black goats, as a native breed of south China's tropical island province of Hainan, possess distinctive traits such as black hair, a moderate growth rate, good meat quality, and small body size. However, they exhibit exceptional resilience to rough feeding conditions, possess high-quality meat, and show remarkable resistance to stress and heat. In this study, we resequenced the whole genome of Hainan black goats to study the economic traits and genetic basis of these goats, we leveraged whole-genome sequencing data from 33 Hainan black goats to analyze single nucleotide polymorphism (SNP) density, Runs of homozygosity (ROH), Integrated Haplotype Score (iHS), effective population size (Ne), Nucleotide diversity Analysis (Pi) and selection characteristics. Our findings revealed that Hainan black goats harbor a substantial degree of genetic variation, with a total of 23 608 983 SNPs identified. Analysis of ROHs identified 53 710 segments, predominantly composed of short fragments, with inbreeding events mainly occurring in ancient ancestors, the estimates of inbreeding based on ROH in Hainan black goats typically exhibit moderate values ranging from 0.107 to 0.186. This is primarily attributed to significant declines in the effective population size over recent generations. Moreover, we identified 921 candidate genes within the intersection candidate region of ROH and iHS. Several of these genes are associated with crucial traits such as immunity (PTPRC, HYAL1, HYAL2, HYAL3, CENPE and PKN1), heat tolerance (GNG2, MAPK8, CAPN2, SLC1A1 and LEPR), meat quality (ACOX1, SSTR1, CAMK2B, PPP2CA and PGM1), cashmere production (AKT4, CHRM2, OXTR, AKT3, HMCN1 and CDK19), and stress resistance (TLR2, IFI44, ENPP1, STK3 and NFATC1). The presence of these genes may be attributed to the genetic adaptation of Hainan black goats to local climate conditions. The insights gained from this study provide valuable references and a solid foundation for the preservation, breeding, and utilization of Hainan black goats and their valuable genetic resources.

Phylogeography of Pterocarya hupehensis reveals the evolutionary patterns of a Cenozoic relict tree around the Sichuan Basin

Environmental factors such as mountain tectonic movements and monsoons can enhance genetic differentiation by hindering inter- and intra-specific gene flow. However, the phylogeographic breaks detected within species may differ depending on the different molecular markers used, and biological traits may be a major confounding factor. Pterocarya hupehensis is a vulnerable relict species distributed throughout the Sichuan Basin. Here, we investigated the phylogeographic patterns and evolutionary history of P. hupehensis using chloroplast DNA and restriction site-associated DNA sequencing data from 18 populations around the Sichuan Basin. The 24 chloroplast haplotypes separated into western and eastern lineages at approximately 16.7 Mya, largely coincident with a strengthening of the East Asian monsoon system during the early to middle Miocene. Both cpDNA and nuclear DNA datasets consistently identified distinct western and eastern lineages whose phylogeographic break conformed to the boundary of the Sino-Himalayan and Sino-Japanese forest sub-kingdoms. However, in contrast to the nuclear gene data, the cpDNA data revealed further divergence of the eastern lineage into northern and southern groups along the Yangtze River, a result that likely reflects differences in the extent of pollen vs seed dispersal. During the temperature decline in the penultimate (Riss) glacial period of the Pleistocene epoch, P. hupehensis experienced a genetic bottleneck event, and ecological niche modeling suggests that a subsequent population expansion occurred during the last interglacial period. Our findings not only establish a basis for conservation of this species, but also serve as a case study for the effects of geography and climate change on the evolutionary history of wind-pollinated relict plants.

High quality genomes produced from single MinION flow cells clarify polyploid and demographic histories of critically endangered Fraxinus (ash) species

With populations of threatened and endangered species declining worldwide, efforts are being made to generate high quality genomic records of these species before they are lost forever. Here, we demonstrate that data from single Oxford Nanopore Technologies (ONT) MinION flow cells can, even in the absence of highly accurate short DNA-read polishing, produce high quality de novo plant genome assemblies adequate for downstream analyses, such as synteny and ploidy evaluations, paleodemographic analyses, and phylogenomics. This study focuses on three North American ash tree species in the genus Fraxinus (Oleaceae) that were recently added to the International Union for Conservation of Nature (IUCN) Red List as critically endangered. Our results support a hexaploidy event at the base of the Oleaceae as well as a subsequent whole genome duplication shared by Syringa, Osmanthus, Olea, and Fraxinus. Finally, we demonstrate the use of ONT long-read sequencing data to reveal patterns in demographic history.

A Multiplex PCR System of Novel Microsatellite Loci for Population Genetic Application in Walnuts

Multiplex polymerase chain reaction (PCR) of microsatellite loci allows for simultaneous amplification of two or more pairs of primers in a single PCR reaction; hence, it is cost and time effective. However, very few attempts have been reported in non-model species. In this study, by combining a genome-based de novo development and cross-species application approach, a multiplex PCR system comprising 5 PCR reactions of 33 microsatellites consisting of 26 novel genomic and 7 literature-sourced loci was tested for polymorphisms, cross-species transferability, and the ability to assess genetic diversity and population structure of three walnut species (Juglans spp.). We found that the genome-based approach is more efficient than other methods. An allelic ladder was developed for each locus to enhance consistent genotyping among laboratories. The population genetic analysis results showed that all 33 loci were successfully transferred across the three species, showing high polymorphism and a strong genetic structure. Hence, the multiplex PCR system is highly applicable in walnut species. Furthermore, we propose an efficient pipeline to characterize and genotype polymorphic microsatellite loci. The novel toolbox developed here will aid future ecology and evolution studies in walnut and could serve as a model for other plant species.

Genomic insights into biased allele loss and increased gene numbers after genome duplication in autotetraploid Cyclocarya paliurus

BACKGROUND

Autopolyploidy is a valuable model for studying whole-genome duplication (WGD) without hybridization, yet little is known about the genomic structural and functional changes that occur in autopolyploids after WGD. Cyclocarya paliurus (Juglandaceae) is a natural diploid-autotetraploid species. We generated an allele-aware autotetraploid genome, a chimeric chromosome-level diploid genome, and whole-genome resequencing data for 106 autotetraploid individuals at an average depth of 60 × per individual, along with 12 diploid individuals at an average depth of 90 × per individual.

RESULTS

Autotetraploid C. paliurus had 64 chromosomes clustered into 16 homologous groups, and the majority of homologous chromosomes demonstrated similar chromosome length, gene numbers, and expression. The regions of synteny, structural variation and nonalignment to the diploid genome accounted for 81.3%, 8.8% and 9.9% of the autotetraploid genome, respectively. Our analyses identified 20,626 genes (69.18%) with four alleles and 9191 genes (30.82%) with one, two, or three alleles, suggesting post-polyploid allelic loss. Genes with allelic loss were found to occur more often in proximity to or within structural variations and exhibited a marked overlap with transposable elements. Additionally, such genes showed a reduced tendency to interact with other genes. We also found 102 genes with more than four copies in the autotetraploid genome, and their expression levels were significantly higher than their diploid counterparts. These genes were enriched in enzymes involved in stress response and plant defense, potentially contributing to the evolutionary success of autotetraploids. Our population genomic analyses suggested a single origin of autotetraploids and recent divergence (~ 0.57 Mya) from diploids, with minimal interploidy admixture.

CONCLUSIONS

Our results indicate the potential for genomic and functional reorganization, which may contribute to evolutionary success in autotetraploid C. paliurus.

Genomic Insights into Adaptation to Karst Limestone and Incipient Speciation in East Asian Platycarya spp. (Juglandaceae)

When challenged by similar environmental conditions, phylogenetically distant taxa often independently evolve similar traits (convergent evolution). Meanwhile, adaptation to extreme habitats might lead to divergence between taxa that are otherwise closely related. These processes have long existed in the conceptual sphere, yet molecular evidence, especially for woody perennials, is scarce. The karst endemic Platycarya longipes and its only congeneric species, Platycarya strobilacea, which is widely distributed in the mountains in East Asia, provide an ideal model for examining the molecular basis of both convergent evolution and speciation. Using chromosome-level genome assemblies of both species, and whole-genome resequencing data from 207 individuals spanning their entire distribution range, we demonstrate that P. longipes and P. strobilacea form two species-specific clades, which diverged around 2.09 million years ago. We find an excess of genomic regions exhibiting extreme interspecific differentiation, potentially due to long-term selection in P. longipes, likely contributing to the incipient speciation of the genus Platycarya. Interestingly, our results unveil underlying karst adaptation in both copies of the calcium influx channel gene TPC1 in P. longipes. TPC1 has previously been identified as a selective target in certain karst-endemic herbs, indicating a convergent adaptation to high calcium stress among karst-endemic species. Our study reveals the genic convergence of TPC1 among karst endemics and the driving forces underneath the incipient speciation of the two Platycarya lineages.

Pan-genome and transcriptome analyses provide insights into genomic variation and differential gene expression profiles related to disease resistance and fatty acid biosynthesis in eastern black walnut (Juglans nigra)

Walnut (Juglans) species are used as nut crops worldwide. Eastern black walnut (EBW, Juglans nigra), a diploid, horticultural important woody species is native to much of eastern North America. Although it is highly valued for its wood and nut, there are few resources for understanding EBW genetics. Here, we present a high-quality genome assembly of J. nigra based on Illumina, Pacbio, and Hi-C technologies. The genome size was 540.8 Mb, with a scaffold N50 size of 35.1 Mb, and 99.0% of the assembly was anchored to 16 chromosomes. Using this genome as a reference, the resequencing of 74 accessions revealed the effective population size of J. nigra declined during the glacial maximum. A single whole-genome duplication event was identified in the J. nigra genome. Large syntenic blocks among J. nigra, Juglans regia, and Juglans microcarpa predominated, but inversions of more than 600 kb were identified. By comparing the EBW genome with those of J. regia and J. microcarpa, we detected InDel sizes of 34.9 Mb in J. regia and 18.3 Mb in J. microcarpa, respectively. Transcriptomic analysis of differentially expressed genes identified five presumed NBS-LRR (NUCLEOTIDE BINDING SITE-LEUCINE-RICH REPEAT) genes were upregulated during the development of walnut husks and shells compared to developing embryos. We also identified candidate genes with essential roles in seed oil synthesis, including FAD (FATTY ACID DESATURASE) and OLE (OLEOSIN). Our work advances the understanding of fatty acid bioaccumulation and disease resistance in nut crops, and also provides an essential resource for conducting genomics-enabled breeding in walnut.

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.

Modeled distribution shifts of North American birds over four decades based on suitable climate alone do not predict observed shifts

As climate change alters the global environment, it is critical to understand the relationship between shifting climate suitability and species distributions. Key questions include whether observed changes in population abundance are aligned with the velocity and direction of shifts predicted by climate suitability models and if the responses are consistent among species with similar ecological traits. We examined the direction and velocity of the observed abundance-based distribution centroids compared with the model-predicted bioclimatic distribution centroids of 250 bird species across the United States from 1969 to 2011. We hypothesized that there is a significant positive correlation in both direction and velocity between the observed and the modeled shifts. We then tested five additional hypotheses that predicted differential shifting velocity based on ecological adaptability and climate change exposure. Contrary to our hypotheses, we found large differences between the observed and modeled shifts among all studied bird species and within specific ecological guilds. However, temperate migrants and habitat generalist species tended to have higher velocity of observed shifts than other species. Neotropical migratory and wetland birds also had significantly different observed velocities than their counterparts, which may be due to their climate change exposure. The velocity based on modeled bioclimatic suitability did not exhibit significant differences among most guilds. Boreal forest birds were the only guild with significantly faster modeled-shifts than the other groups, suggesting an elevated conservation risk for high latitude and altitude species. The highly idiosyncratic species responses to climate and the mismatch between shifts in modeled and observed distribution centroids highlight the challenge of predicting species distribution change based solely on climate suitability and the importance of non-climatic factors traits in shaping species distributions.

Characteristics and potential functional effects of long insertions in Asian butternuts

Background

Structural variants (SVs) play important roles in adaptation evolution and species diversification. Especially, in plants, many phenotypes of response to the environment were found to be associated with SVs. Despite the prevalence and significance of SVs, long insertions remain poorly detected and studied in all but model species.

Results

We used whole-genome resequencing of paired reads from 80 Asian butternuts to detect long insertions and further analyse their characteristics and potential functional effects. By combining of mapping-based and de novo assembly-based methods, we obtained a multiple related species pangenome representing higher taxonomic groups. We obtained 89,312 distinct contigs totaling 147,773,999 base pair (bp) of new sequences, of which 347 were putative long insertions placed in the reference genome. Most of the putative long insertions appeared in multiple species; in contrast, only 62 putative long insertions appeared in one species, which may be involved in the response to the environment. 65 putative long insertions fell into 61 distinct protein-coding genes involved in plant development, and 105 putative long insertions fell into upstream of 106 distinct protein-coding genes involved in cellular respiration. 3,367 genes were annotated in 2,606 contigs. We propose PLAINS (https://github.com/CMB-BNU/PLAINS.git), a streamlined, comprehensive pipeline for the prediction and analysis of long insertions using whole-genome resequencing.

Conclusions

Our study lays down an important foundation for further whole-genome long insertion studies, allowing the investigation of their effects by experiments.

Population-genomic analyses reveal bottlenecks and asymmetric introgression from Persian into iron walnut during domestication

Background

Persian walnut, Juglans regia, occurs naturally from Greece to western China, while its closest relative, the iron walnut, Juglans sigillata, is endemic in southwest China; both species are cultivated for their nuts and wood. Here, we infer their demographic histories and the time and direction of possible hybridization and introgression between them.

Results

We use whole-genome resequencing data, different population-genetic approaches (PSMC and GONE), and isolation-with-migration models (IMa3) on individuals from Europe, Iran, Kazakhstan, Pakistan, and China. IMa3 analyses indicate that the two species diverged from each other by 0.85 million years ago, with unidirectional gene flow from eastern J. regia and its ancestor into J. sigillata, including the shell-thickness gene. Within J. regia, a western group, located from Europe to Iran, and an eastern group with individuals from northern China, experienced dramatically declining population sizes about 80 generations ago (roughly 2400 to 4000 years), followed by an expansion at about 40 generations, while J. sigillata had a constant population size from about 100 to 20 generations ago, followed by a rapid decline.

Conclusions

Both J. regia and J. sigillata appear to have suffered sudden population declines during their domestication, suggesting that the bottleneck scenario of plant domestication may well apply in at least some perennial crop species. Introgression from introduced J. regia appears to have played a role in the domestication of J. sigillata.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02720-z.

Genetic diversity and population dynamic of Ziziphus jujuba var. spinosa (Bunge) Hu ex H. F. Chow in Central China

Phylogeographic research concerning Central China has been rarely conducted. Population genetic and phylogeography of Ziziphus jujuba var. spinosa (also called sour jujube) were investigated to improve our understanding of plant phylogeographic patterns in Central China. Single-copy nuclear gene markers and complete chloroplast genome data were applied to 328 individuals collected from 21 natural populations of sour jujube in China. Nucleotide variation of sour jujube was relatively high (π = 0.00720, θ w = 0.00925), which resulted from the mating system and complex population dynamics. Analysis of molecular variation analysis revealed that most of the total variation was attributed to variation within populations, and a high level of genetic differentiation among populations was detected (F st = 0.197). Relatively low long-distance dispersal capability and vitality of pollen contributed to high genetic differentiation among populations. Differences in the environmental conditions and long distance among populations further restricted gene flow. Structure clustering analysis uncovered intraspecific divergence between central and marginal populations. Migrate analysis found a high level of gene flow between these two intraspecific groups. Bayesian skyline plot detected population expansion of these two intraspecific groups. Network and phylogeny analysis of chloroplast haplotypes also found intraspecific divergence, and the divergence time was estimated to occur at about 55.86 Ma. Haplotype native to the Loess Plateau was more ancient, and multiple glacial refugia of sour jujube were found to locate at the Loess Plateau, areas adjacent to the Qinling Mountains and Tianmu Mountains. Species distribution model analysis found a typical contraction-expansion model corresponding to the Quaternary climatic oscillations. In the future, the distribution of sour jujube may shift to high-latitude areas. This study provides new insights for phylogeographic research of temperate plant species distributed in Central China and sets a solid foundation for the application of the scientific management strategy of Z. jujuba var. spinosa.

Genetic Diversity of Juglans mandshurica Populations in Northeast China Based on SSR Markers

Juglans mandshurica is a native tree species in Northeast China. Due to habitat destruction and human disturbance, its population size has sharply decreased. Currently, information on molecular markers of J. mandshurica is limited and cannot meet the needs of germplasm resource evaluation and molecular marker-assisted breeding of J. mandshurica. Based on transcriptomic data from three tissues (leaves, bark, and fruit pericarp), we developed expressed sequence tag-simple sequence repeats (EST-SSRs) for J. mandshurica, and 15 polymorphic EST-SSR primers were initially selected. The average number of alleles (Na), expected heterozygosity (He), and the polymorphic information content (PIC) at different loci were 18.27, 0.670, and 0.797, respectively. Population genetic diversity analysis revealed that the average Na, He, and Shannon information indices (I) for 15 J. mandshurica populations were 6.993, 0.670, and 1.455, respectively. Among them, population Hunchun exhibited the highest genetic diversity (Na = 7.933, He = 0.723, and I = 1.617), while population Heihe exhibited the lowest genetic diversity (Na = 4.200, He = 0.605, and I = 1.158). STRUCTURE analysis, neighbor-joining method cluster analysis, and principal coordinate analysis showed that the 343 individuals of J. mandshurica from 15 populations were clustered into three categories. Category 1 (green) had 147 individuals from eight populations in Qingyuan, Caohekou, Jian, Ningan, Yongji, Baishishan, Helong, and Maoershan; category 2 (blue) had 81 individuals from three populations in Hulin, Boli, and Sanchazi; and category 3 (red) had 115 individuals from four populations in Heihe, Hunchun, Fangzheng, and Liangshui. Analysis of molecular variance (AMOVA) showed that genetic variations among and within individuals accounted for 16.22% and 21.10% of the total genetic variation, respectively, indicating that genetic variations within populations were greater than genetic variations among populations. The average genetic differentiation coefficient (Fst) and gene flow (Nm) between different populations were 0.109 and 4.063, respectively, implying moderate levels of genetic differentiation and gene flow. Based on the genetic diversity characteristics of different populations, we proposed various genetic conservation strategies for J. mandshurica.

The Manchurian Walnut Genome: Insights into Juglone and Lipid Biosynthesis

BACKGROUND

Manchurian walnut (Juglans mandshurica Maxim.) is a tree with multiple industrial uses and medicinal properties in the Juglandaceae family (walnuts and hickories). J. mandshurica produces juglone, which is a toxic allelopathic agent and has potential utilization value. Furthermore, the seed of J. mandshurica is rich in various unsaturated fatty acids and has high nutritive value.

FINDINGS

Here, we present a high-quality chromosome-scale reference genome assembly and annotation for J. mandshurica (n = 16) with a contig N50 of 21.4 Mb by combining PacBio high-fidelity reads with high-throughput chromosome conformation capture data. The assembled genome has an estimated sequence size of 548.7 Mb and consists of 657 contigs, 623 scaffolds, and 40,453 protein-coding genes. In total, 60.99% of the assembled genome consists of repetitive sequences. Sixteen super-scaffolds corresponding to the 16 chromosomes were assembled, with a scaffold N50 length of 33.7 Mb and a BUSCO complete gene percentage of 98.3%. J. mandshurica displays a close sequence relationship with Juglans cathayensis, with a divergence time of 13.8 million years ago. Combining the high-quality genome, transcriptome, and metabolomics data, we constructed a gene-to-metabolite network and identified 566 core and conserved differentially expressed genes, which may be involved in juglone biosynthesis. Five CYP450 genes were found that may contribute to juglone accumulation. NAC, bZip, NF-YA, and NF-YC are positively correlated with the juglone content. Some candidate regulators (e.g., FUS3, ABI3, LEC2, and WRI1 transcription factors) involved in the regulation of lipid biosynthesis were also identified.

CONCLUSIONS

Our genomic data provide new insights into the evolution of the walnut genome and create a new platform for accelerating molecular breeding and improving the comprehensive utilization of these economically important tree species.

Genomic variation reveals demographic history and biological adaptation of the ancient relictual, lotus (Nelumbo Adans)

Lotus (Nelumbo Adans.), a relict plant, is the testimony of long-term sustained ecological success, but the underlying genetic changes related to its survival strategy remains unclear. Here, we assembled the high-quality lotus genome, investigated genome variation of lotus mutation accumulation (MA) lines and reconstructed the demographic history of wild Asian lotus, respectively. We identified and validated 43 base substitutions fixed in MA lines, implying a spontaneous mutation rate of 1.4 × 10-9 base/generation in lotus shoot stem cells. The past history of lotus revealed that the ancestors of lotus in eastern and southern Asia could be traced back ~20 million years ago (Mya) and experienced twice significant bottlenecks and population splits. We further identified the selected genes among three lotus groups in different habitats, suggesting that 453 genes between tropical and temperate group and 410 genes between two subgroups from Northeastern China and the Yangtze River - Yellow River Basin might play important roles in natural selection in lotus's adaptation and resilience. Our findings not only improve an understanding of the lotus evolutionary history and the genetic basis of its survival advantages, but also provide valuable data for addressing various questions in evolution and protection for the relict plants.

Dead-End Hybridization in Walnut Trees Revealed by Large-Scale Genomic Sequence Data

Although hybridization plays a large role in speciation, some unknown fraction of hybrid individuals never reproduces, instead remaining as genetic dead-ends. We investigated a morphologically distinct and culturally important Chinese walnut, Juglans hopeiensis, suspected to have arisen from hybridization of Persian walnut (J. regia) with Asian butternuts (J. cathayensis, J. mandshurica, and hybrids between J. cathayensis and J. mandshurica). Based on 151 whole-genome sequences of the relevant taxa, we discovered that all J. hopeiensis individuals are first-generation hybrids, with the time for the onset of gene flow estimated as 370,000 years, implying both strong postzygotic barriers and the presence of J. regia in China by that time. Six inversion regions enriched for genes associated with pollen germination and pollen tube growth may be involved in the postzygotic barriers that prevent sexual reproduction in the hybrids. Despite its long-recurrent origination and distinct traits, J. hopeiensis does not appear on the way to speciation.

Chloroplast and Nuclear Genetic Diversity Explain the Limited Distribution of Endangered and Endemic Thuja sutchuenensis in China

Narrow-ranged species face challenges from natural disasters and human activities, and to address why species distributes only in a limited region is of great significance. Here we investigated the genetic diversity, gene flow, and genetic differentiation in six wild and three cultivated populations of Thuja sutchuenensis, a species that survive only in the Daba mountain chain, using chloroplast simple sequence repeats (cpSSR) and nuclear restriction site-associated DNA sequencing (nRAD-seq). Wild T. sutchuenensis populations were from a common ancestral population at 203 ka, indicating they reached the Daba mountain chain before the start of population contraction at the Last Interglacial (LIG, ∼120-140 ka). T. sutchuenensis populations showed relatively high chloroplast but low nuclear genetic diversity. The genetic differentiation of nRAD-seq in any pairwise comparisons were low, while the cpSSR genetic differentiation values varied with pairwise comparisons of populations. High gene flow and low genetic differentiation resulted in a weak isolation-by-distance effect. The genetic diversity and differentiation of T. sutchuenensis explained its survival in the Daba mountain chain, while its narrow ecological niche from the relatively isolated and unique environment in the "refugia" limited its distribution.

Whole genome based insights into the phylogeny and evolution of the Juglandaceae

BACKGROUND

The walnut family (Juglandaceae) contains commercially important woody trees commonly called walnut, wingnut, pecan and hickory. Phylogenetic relationships and diversification within the Juglandaceae are classic and hot scientific topics that have been elucidated by recent fossil, morphological, molecular, and (paleo) environmental data. Further resolution of relationships among and within genera is still needed and can be achieved by analysis of the variation of chloroplast, mtDNA, and nuclear genomes.

RESULTS

We reconstructed the backbone phylogenetic relationships of Juglandaceae using organelle and nuclear genome data from 27 species. The divergence time of Juglandaceae was estimated to be 78.7 Mya. The major lineages diversified in warm and dry habitats during the mid-Paleocene and early Eocene. The plastid, mitochondrial, and nuclear phylogenetic analyses all revealed three subfamilies, i.e., Juglandoideae, Engelhardioideae, Rhoipteleoideae. Five genera of Juglandoideae were strongly supported. Juglandaceae were estimated to have originated during the late Cretaceous, while Juglandoideae were estimated to have originated during the Paleocene, with evidence for rapid diversification events during several glacial and geological periods. The phylogenetic analyses of organelle sequences and nuclear genome yielded highly supported incongruence positions for J. cinerea, J. hopeiensis, and Platycarya strobilacea. Winged fruit were the ancestral condition in the Juglandoideae, but adaptation to novel dispersal and regeneration regimes after the Cretaceous-Paleogene boundary led to the independent evolution of zoochory among several genera of the Juglandaceae.

CONCLUSIONS

A fully resolved, strongly supported, time-calibrated phylogenetic tree of Juglandaceae can provide an important framework for studying classification, diversification, biogeography, and comparative genomics of plant lineages. Our addition of new, annotated whole chloroplast genomic sequences and identification of their variability informs the study of their evolution in walnuts (Juglandaceae).

Whole genome based insights into the phylogeny and evolution of the Juglandaceae

BACKGROUND

The walnut family (Juglandaceae) contains commercially important woody trees commonly called walnut, wingnut, pecan and hickory. Phylogenetic relationships and diversification within the Juglandaceae are classic and hot scientific topics that have been elucidated by recent fossil, morphological, molecular, and (paleo) environmental data. Further resolution of relationships among and within genera is still needed and can be achieved by analysis of the variation of chloroplast, mtDNA, and nuclear genomes.

RESULTS

We reconstructed the backbone phylogenetic relationships of Juglandaceae using organelle and nuclear genome data from 27 species. The divergence time of Juglandaceae was estimated to be 78.7 Mya. The major lineages diversified in warm and dry habitats during the mid-Paleocene and early Eocene. The plastid, mitochondrial, and nuclear phylogenetic analyses all revealed three subfamilies, i.e., Juglandoideae, Engelhardioideae, Rhoipteleoideae. Five genera of Juglandoideae were strongly supported. Juglandaceae were estimated to have originated during the late Cretaceous, while Juglandoideae were estimated to have originated during the Paleocene, with evidence for rapid diversification events during several glacial and geological periods. The phylogenetic analyses of organelle sequences and nuclear genome yielded highly supported incongruence positions for J. cinerea, J. hopeiensis, and Platycarya strobilacea. Winged fruit were the ancestral condition in the Juglandoideae, but adaptation to novel dispersal and regeneration regimes after the Cretaceous-Paleogene boundary led to the independent evolution of zoochory among several genera of the Juglandaceae.

CONCLUSIONS

A fully resolved, strongly supported, time-calibrated phylogenetic tree of Juglandaceae can provide an important framework for studying classification, diversification, biogeography, and comparative genomics of plant lineages. Our addition of new, annotated whole chloroplast genomic sequences and identification of their variability informs the study of their evolution in walnuts (Juglandaceae).

Different rates of pollen and seed gene flow cause branch-length and geographic cytonuclear discordance within Asian butternuts

Topological cytonuclear discordance is commonly observed in plant phylogenetic and phylogeographic studies, yet few studies have attempted to detect two other forms of cytonuclear discordance (branch length and geographical) and to uncover the causes of the discordance. We used the whole nuclear and chloroplast genome data from 80 individual Asian butternuts to reveal the pattern and processes of cytonuclear discordance. Our findings indicate that the chloroplast genome had substantially deeper divergence (branch-length discordance) and a steeper cline in the contact zone (geographic discordance) compared with the nuclear genome. After various hypothesis have been tested, the results suggest that incomplete lineage sorting, positive selection and cytonuclear incompatibility are probably insufficient to explain this pattern. However, isolation-by-distance analysis and gene flow estimation point to a much higher level of gene flow by pollen compared with by seeds, which may have slowed down lineage divergence and mediated wider contact for nuclear genome compared with the chloroplast genome. Altogether, this study highlights a critical role of sex-biased dispersal in causing discordance between the nuclear and plastid genome of Asian butternuts. Given its ubiquity among plants, asymmetric gene flow should be given a high priority in future studies of cytonuclear discordance.

Population dynamics of caribou shaped by glacial cycles before the last glacial maximum

Pleistocene glacial cycles influenced the diversification of high‐latitude wildlife species through recurrent periods of range contraction, isolation, divergence, and expansion from refugia and subsequent admixture of refugial populations. We investigate population size changes and the introgressive history of caribou (Rangifer tarandus) in western Canada using 33 whole genome sequences coupled with larger‐scale mitochondrial data. We found that a major population expansion of caribou occurred starting around 110,000 years ago (kya), the start of the last glacial period. Additionally, we found effective population sizes of some caribou reaching ~700,000 to 1,000,000 individuals, one of the highest recorded historical effective population sizes for any mammal species thus far. Mitochondrial analyses dated introgression events prior to the LGM dating to 20–30 kya and even more ancient at 60 kya, coinciding with colder periods with extensive ice coverage, further demonstrating the importance of glacial cycles and events prior to the LGM in shaping demographic history. Reconstructing the origins and differential introgressive history has implications for predictions on species responses under climate change. Our results have implications for other whole genome analyses using pairwise sequentially Markovian coalescent (PSMC) analyses, as well as highlighting the need to investigate pre‐LGM demographic patterns to fully reconstruct the origin of species diversity, especially for high‐latitude species.

Chromosome-level genome assembly and population genetic analysis of a critically endangered rhododendron provide insights into its conservation

Rhododendrons are woody plants, famous throughout the world as having high horticultural value. However, many wild species are currently threatened with extinction. Here, we report for the first time a high-quality, chromosome-level genome of Rhododendron griersonianum, which has contributed to approximately 10% of all horticultural rhododendron varieties but which in its wild form has been evaluated as critically endangered. The final genome assembly, which has a contig N50 size of approximately 34 M and a total length of 677 M, is the highest-quality genome sequenced within the genus to date, in part due to its low heterozygosity (0.18%). Identified repeats constitute approximately 57% of the genome, and 38 280 protein-coding genes were predicted with high support. We further resequenced 31 individuals of R. griersonianum as well as 30 individuals of its widespread relative R. delavayi, and performed additional conservation genomic analysis. The results showed that R. griersonianum had lower genetic diversity (θ = 2.58e-3; π = 1.94e-3) when compared not only to R. delavayi (θ = 11.61e-3, π = 12.97e-3), but also to most other woody plants. Furthermore, three severe genetic bottlenecks were detected using both the Stairway plot and fastsimcoal2 analysis, which are thought to have occurred in the late Middle Pleistocene and the Last Glacial Maximum (LGM) period. After these bottlenecks, R. griersonianum recovered and maintained a constant effective population size (>25 000) until now. Intriguingly, R. griersonianum has accumulated significantly more deleterious mutations in the homozygous state than R. delavayi, and several deleterious mutations (e.g., in genes involved in the response to heat stress) are likely to have harmed the adaptation of this plant to its surroundings. This high-quality, chromosome-level genome and the population genomic analysis of the critically endangered R. griersonianum will provide an invaluable resource as well as insights for future study in this species to facilitate conservation and in the genus Rhododendron in general.

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

Quaternary climate oscillations and geographical heterogeneity play important roles in determining species and genetic diversity distribution patterns, but how these factors affect the migration and differentiation of East Asian plants species at the population level remains poorly understood. The Aquilegia ecalcarata complex, a group that originated in the Late Tertiary and is widely distributed throughout East Asia, displays high genetic variation that is suitable for studying elaborate phylogeographic patterns and demographic history related to the impact of Quaternary climate and geography. We used plastid genome data from 322 individuals in 60 populations of the A. ecalcarata complex to thoroughly explore the impact of Quaternary climate oscillations and geography on the phylogeographic patterns and demographic history of the A. ecalcarata complex through a series of phylogenetic, divergence time estimation, and demographic history analyses. The dry, cold climate and frequent climate oscillations that occurred during the early Pleistocene and the Mid-Pleistocene transition led to the differentiation of the A. ecalcarata complex, which was isolated in various areas. Geographically, the A. ecalcarata complex can be divided into Eastern and Western Clades and five subclades, which conform to the divergence of the East Asian flora. Our results clearly show the impact of Quaternary climate and geography on evolutionary history at the population level. These findings promote the understanding of the relationship between plant genetic differentiation and climate and geographical factors of East Asia at the population level.

Repetitive genomic regions and the inference of demographic history

Inference of demographic histories using whole-genome datasets has provided insights into diversification, adaptation, hybridization, and plant-pathogen interactions, and stimulated debate on the impact of anthropogenic interventions and past climate on species demography. However, the impact of repetitive genomic regions on these inferences has mostly been ignored by masking of repeats. We use the Populus trichocarpa genome (Pop_tri_v3) to show that masking of repeat regions leads to lower estimates of effective population size (Ne) in the distant past in contrast to an increase in Ne estimates in recent times. However, in human datasets, masking of repeats resulted in lower estimates of Ne at all time points. We demonstrate that repeats affect demographic inferences using diverse methods like PSMC, MSMC, SMC++, and the Stairway plot. Our genomic analysis revealed that the biases in Ne estimates were dependent on the repeat class type and its abundance in each atomic interval. Notably, we observed a weak, yet consistently significant negative correlation between the repeat abundance of an atomic interval and the Ne estimates for that interval, which potentially reflects the recombination rate variation within the genome. The rationale for the masking of repeats has been that variants identified within these regions are erroneous. We find that polymorphisms in some repeat classes occur in callable regions and reflect reliable coalescence histories (e.g., LTR Gypsy, LTR Copia). The current demography inference methods do not handle repeats explicitly, and hence the effect of individual repeat classes needs careful consideration in comparative analysis. Deciphering the repeat demographic histories might provide a clear understanding of the processes involved in repeat accumulation.

Recent demographic histories of temperate deciduous trees inferred from microsatellite markers

Background

Accurate inference of demographic histories for temperate tree species can aid our understanding of current climate change as a driver of evolution. Microsatellites are more suitable for inferring recent historical events due to their high mutation rates. However, most programs analyzing microsatellite data assume a strict stepwise mutation model (SMM), which could cause false detection of population shrinkage when microsatellite mutation does not follow SMM.

Results

This study aims to reconstruct the recent demographic histories of five cool-temperate tree species in Eastern Asia, Quercus mongolica, Q. liaotungensis, Juglans cathayensis, J. mandshurica and J. ailantifolia, by using 19 microsatellite markers with two methods considering generalized stepwise mutation model (GSM) (MIGRAINE and VarEff). Both programs revealed that all the five species experienced expansions after the Last Glacial Maximum (LGM). Within butternuts, J. cathayensis experienced a more serious bottleneck than the other species, and within oaks, Q. mongolica showed a moderate increase in population size and remained stable after the expansion. In addition, the point estimates of the multistep mutation proportion in the GSM model (p_GSM) for all five species were between 0.50 and 0.65, indicating that when inferring population demographic history of the cool-temperate forest species using microsatellite markers, it is better to assume a GSM rather than a SMM.

Conclusions

This study provides the first direct evidence that five cool-temperate tree species in East Asia have experienced expansions after the LGM with microsatellite data. Considering the mutation model of microsatellite has a vital influence on demographic inference, combining multiple programs such as MIGRAINE and VarEff can effectively reduce errors caused by inappropriate model selection and prior setting.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01805-w.

Improved de novo chromosome-level genome assembly of the vulnerable walnut tree Juglans mandshurica reveals gene family evolution and possible genome basis of resistance to lesion nematode

Manchurian walnut (Juglans mandshurica Maxim.) is a synonym of J. cathayensis, a diploid, vulnerable, temperate deciduous tree valued for its wood and nut. It is also valued as a rootstock for Juglans regia because of its reported tolerance of lesion nematode. Reference genomes are available for several Juglans species, our goal was to produce a de novo, chromosome-level assembly of the J. mandshurica genome. Here, we reported an improved assembly of J. mandshurica with a contig N50 size of 6.49 Mb and a scaffold N50 size of 36.1 Mb. The total genome size was 548 Mb encoding 29,032 protein coding genes which were annotated. The collinearity analysis showed that J. mandshurica and J. regia originated from a common ancestor, with both species undergoing two WGD events. A genomic comparison showed that J. mandshurica was missing 1657 genes found in J. regia, and J. mandshurica includes 2827 genes not found in of the J. regia genome. The J. mandshurica contained 1440 unique paralogues that were highly enriched for flavonoid biosynthesis, phenylpropanoid biosynthesis, and plant-pathogen interaction. Four gene families related to disease resistance notable contraction (rapidly evolving; LEA, WAK, PPR, and PR) in J. mandshurica compared to eight species. JmaPR10 and JmaPR8 contained three orthologous gene pairs with J. regia that were highly expressed in root bark. JmaPR10 is a strong candidate gene for lesion nematodes resistance in J. mandshurica. The J. mandshurica genome should be a useful resource for study of the evolution, breeding, and genetic variation in walnuts (Juglans).

Inferring historical survivals of climate relicts: the effects of climate changes, geography, and population-specific factors on herbaceous hydrangeas

Climate relicts hold considerable importance because they have resulted from numerous historical changes. However, there are major interspecific variations among the ways by which they survived climate changes. Therefore, investigating the factors and timing that affected population demographics can expand our understanding of how climate relicts responded to historical environmental changes. Here, we examined herbaceous hydrangeas of genus Deinanthe in East Asia, which show limited distributions and a remarkable disjunction between Japan and central China. Chloroplast genome and restriction site-associated DNA sequencing revealed that speciation event occurred in the late Miocene (ca. 7-9 Mya) in response to global climate change. Two lineages apparently remained not branched until the middle Quaternary, and afterwards started to diverge to regional population groups. The narrow endemic species in central China showed lower genetic diversity (He = 0.082), as its population size rapidly decreased during the Holocene due to isolation in montane refugia. Insular populations in the three Japanese islands (He = 0.137-0.160) showed a genetic structure that was inconsistent with sea barriers, indicating that it was shaped in the glacial period when its range retreated to coastal refugia on the exposed sea floor. Demographic modelling by stairway-plot analysis reconstructed variable responses of Japanese populations: some experienced glacial bottlenecks in refugial isolation, while post-glacial range expansion seemingly exerted founder effects on other populations. Overall, this study demonstrated the involvement of not just one, but multiple factors, such as the interplay between climate changes, geography, and other population-specific factors, that determine the demographics of climate relicts.

Whole-Genome Diversification Analysis of the Hornbeam Species Reveals Speciation and Adaptation Among Closely Related Species

Speciation is the key evolutionary process for generating biological diversity and has a central place in evolutionary and ecological research. How species diverge and adapt to different habitats is one of the most exciting areas in speciation studies. Here, we sequenced 55 individuals from three closely related species in the genus Carpinus: Carpinus tibetana, Carpinus monbeigiana, and Carpinus mollicoma to understand the strength and direction of gene flow and selection during the speciation process. We found low genetic diversity in C. tibetana, which reflects its extremely small effective population size. The speciation analysis between C. monbeigiana and C. mollicoma revealed that both species diverged ∼1.2 Mya with bidirectional gene flow. A total of 291 highly diverged genes, 223 copy number variants genes, and 269 positive selected genes were recovered from the two species. Genes associated with the diverged and positively selected regions were mainly involved in thermoregulation, plant development, and response to stress, which included adaptations to their habitats. We also found a great population decline and a low genetic divergence of C. tibetana, which suggests that this species is extremely vulnerable. We believe that the current diversification and adaption study and the important genomic resource sequenced herein will facilitate the speciation studies and serve as an important methodological reference for future research.

Climatic Refugia and Geographical Isolation Contribute to the Speciation and Genetic Divergence in Himalayan-Hengduan Tree Peonies (Paeonia delavayi and Paeonia ludlowii)

Himalaya and Hengduan Mountains (HHM) is a biodiversity hotspot, and very rich in endemic species. Previous phylogeographical studies proposed different hypotheses (vicariance and climate-driven speciation) in explaining diversification and the observed pattern of extant biodiversity, but it is likely that taxa are forming in this area in species-specific ways. Here, we reexplored the phylogenetic relationship and tested the corresponding hypotheses within Paeonia subsect. Delavayanae composed of one widespread species (Paeonia delavayi) and the other geographically confined species (Paeonia ludlowii). We gathered genetic variation data at three chloroplast DNA fragments and one nuclear gene from 335 individuals of 34 populations sampled from HHM. We performed a combination of population genetic summary statistics, isolation-with-migration divergence models, isolation by environment, and demographic history analyses. We found evidence for the current taxonomic treatment that P. ludlowii and P. delavayi are two different species with significant genetic differentiation. The significant isolation by environment was revealed within all sampled populations but genetic distances only explained by geographical distances within P. delavayi populations. The results of population divergence models and demographic history analyses indicated a progenitor–derivative relationship and the Late Quaternary divergence without gene flow between them. The coalescence of all sampled cpDNA haplotypes could date to the Late Miocene, and P. delavayi populations probably underwent a severe bottleneck in population size during the last glacial period. Genetic variation in Paeonia subsect. Delavayanae is associated with geographical and environmental distances. These findings point to the importance of geological and climatic changes as causes of the speciation event and lineage diversification within Paeonia subsect. Delavayanae.

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

'Living fossils' are testimonies of long-term sustained ecological success, but how demographic history and natural selection contributed to their survival, resilience, and persistence in the face of Quaternary climate fluctuations remains unclear. To better understand the interplay between demographic history and selection in shaping genomic diversity and evolution of such organisms, we assembled the whole genome of Cercidiphyllum japonicum, a widespread East Asian Tertiary relict tree, and resequenced 99 individuals of C. japonicum and its sister species, Cercidiphyllum magnificum (Central Japan). We dated this speciation event to the mid-Miocene, and the intraspecific lineage divergence of C. japonicum (China vs Japan) to the Early Pliocene. Throughout climatic upheavals of the late Tertiary/Quaternary, population bottlenecks greatly reduced the genetic diversity of C. japonicum. However, this polymorphism loss was likely counteracted by, first, long-term balancing selection at multiple chromosomal and heterozygous gene regions, potentially reflecting overdominance, and, second, selective sweeps at stress response and growth-related genes likely involved in local adaptation. Our findings contribute to a better understanding of how living fossils have survived climatic upheaval and maintained an extensive geographic range; that is, both types of selection could be major factors contributing to the species' survival, resilience, and persistence.

Computational Identification and Comparative Analysis of Conserved miRNAs and Their Putative Target Genes in the Juglans regia and J. microcarpa Genomes

MicroRNAs (miRNAs) are important factors for the post-transcriptional regulation of protein-coding genes in plants and animals. They are discovered either by sequencing small RNAs or computationally. We employed a sequence-homology-based computational approach to identify conserved miRNAs and their target genes in Persian (English) walnut, Juglans regia, and its North American wild relative, J. microcarpa. A total of 119 miRNA precursors (pre-miRNAs) were detected in the J. regia genome and 121 in the J. microcarpa genome and miRNA target genes were predicted and their functional annotations were performed in both genomes. In the J. regia genome, 325 different genes were targets; 87.08% were regulated by transcript cleavage and 12.92% by translation repression. In the J. microcarpa genome, 316 different genes were targets; 88.92% were regulated by transcript cleavage and 11.08% were regulated by translation repression. Totals of 1.3% and 2.0% of all resistance gene analogues (RGA) and 2.7% and 2.6% of all transcription factors (TFs) were regulated by miRNAs in the J. regia and J. microcarpa genomes, respectively. Juglans genomes evolved by a whole genome duplication (WGD) and consist of eight pairs of fractionated homoeologous chromosomes. Within each pair, the chromosome that has more genes with greater average transcription also harbors more pre-miRNAs and more target genes than its homoeologue. While only minor differences were detected in pre-miRNAs between the J. regia and J. microcarpa genomes, about one-third of the pre-miRNA loci were not conserved between homoeologous chromosome within each genome. Pre-miRNA and their corresponding target genes showed a tendency to be collocated within a subgenome.

The genome sequence of Mesua ferrea and comparative demographic histories of forest trees

Mesua ferrea (Family: Calophyllaceae) is a tropical forest plant used for timber, biofuel, and traditional medicine. Colloquially, it is known as Nagkesar (Cobra saffron) and is the state flower of Tripura (India). In this study, we perform the whole-genome assembly of Mesua ferrea using ~180X coverage paired-end Illumina data. Our de novo assembly is 614 Mega-base pair (Mbp), has an N50 of 392 Kilo-base pairs (Kbp), and an assembly quality comparable to other published Malpighiales genomes. Further, we collate the genomic datasets of 14 additional forest tree species to compare the temporal dynamics of Effective Population Size (N_e) and find evidence of a substantial bottleneck in all tropical forest plants during Mid-Pleistocene glaciations. The availability of this high-quality draft genome assembly will prove to be a useful resource for functional and comparative genomic studies.

Major QTL with pleiotropic effects controlling time of leaf budburst and flowering-related traits in walnut (Juglans regia L.)

Breeding studies in walnut (Juglans regia L.) are usually time consuming due to the long juvenile period and therefore, this study aimed to determine markers associated with time of leaf budburst and flowering-related traits by performing a genome-wide association study (GWAS). We investigated genotypic variation and its association with time of leaf budburst and flowering-related traits in 188 walnut accessions. Phenotypic data was obtained from 13 different traits during 3 consecutive years. We used DArT-seq for genotyping with a total of 33,519 (14,761 SNP and 18,758 DArT) markers for genome-wide associations to identify marker underlying these traits. Significant correlations were determined among the 13 different traits. Linkage disequilibrium decayed very quickly in walnut in comparison with other plants. Sixteen quantitative trait loci (QTL) with major effects (R² between 0.08 and 0.23) were found to be associated with a minimum of two phenotypic traits each. Of these QTL, QTL05 had the maximum number of associated traits (seven). Our study is GWAS for time of leaf budburst and flowering-related traits in Juglans regia L. and has a strong potential to efficiently implement the identified QTL in walnut breeding programs.

Comparative genomics of six Juglans species reveals disease-associated gene family contractions

Juglans (walnuts), the most speciose genus in the walnut family (Juglandaceae), represents most of the family's commercially valuable fruit and wood-producing trees. It includes several species used as rootstock for their resistance to various abiotic and biotic stressors. We present the full structural and functional genome annotations of six Juglans species and one outgroup within Juglandaceae (Juglans regia, J. cathayensis, J. hindsii, J. microcarpa, J. nigra, J. sigillata and Pterocarya stenoptera) produced using BRAKER2 semi-unsupervised gene prediction pipeline and additional tools. For each annotation, gene predictors were trained using 19 tissue-specific J. regia transcriptomes aligned to the genomes. Additional functional evidence and filters were applied to multi-exonic and mono-exonic putative genes to yield between 27 000 and 44 000 high-confidence gene models per species. Comparison of gene models to the BUSCO embryophyta dataset suggested that, on average, genome annotation completeness was 85.6%. We utilized these high-quality annotations to assess gene family evolution within Juglans, and among Juglans and selected Eurosid species. We found notable contractions in several gene families in J. hindsii, including disease resistance-related wall-associated kinase (WAK), Catharanthus roseus receptor-like kinase (CrRLK1L) and others involved in abiotic stress response. Finally, we confirmed an ancient whole-genome duplication that took place in a common ancestor of Juglandaceae using site substitution comparative analysis.

Morpho-Physiological and Genomic Evaluation of Juglans Species Reveals Regional Maladaptation to Cold Stress

Climate change may have unpredictable effects on the cold hardiness of woody species planted outside of their range of origin. Extreme undulations in temperatures may exacerbate susceptibility to cold stress, thereby interfering with productivity and ecosystem functioning. Juglans L. and their naturally occurring interspecific F1 hybrids, are distributed natively across many temperate regions, and J. regia has been extensively introduced. Cold hardiness, an environmental and genetic factor yet to be evaluated in many native and introduced Juglans species, may be a limiting factor under future climate change and following species introductions. We evaluated cold hardiness of native North American and Eastern Asian Juglans along with J. regia genotypes using field data from the Midwestern United States (Indiana), controlled freezing tests, and genome sequencing with close assessment of Juglans cold hardy genes. Many Juglans species previously screened for cold-hardiness were genotypes derived from the Midwest, California, and Europe. In 2014, despite general climate adaptation, Midwestern winter temperatures of -30°C killed J. regia originating from California; however, naturalized Midwestern J. regia survived and displayed low damage. Hybridization of J. regia with black walnut (J. nigra) and butternut (J. cinerea) produced F1s displaying greater cold tolerance than pure J. regia. Cold hardiness and growth are variable in Midwestern J. regia compared to native Juglans, East Asian Juglans, and F1 hybrids. Phylogeny analyses revealed that J. cinerea sorted with East Asian species using the nuclear genome but with North American species using the organellar genome. Investigation of selected cold hardy genes revealed that J. regia was distinct from other species and exhibited less genetic diversity than native Juglans species Average whole genome heterozygosity and Tajima's D for cold hardy genes was low within J. regia samples and significantly higher for hybrid as well as J. nigra. We confirmed that molecular and morpho-physiological data were highly correlated and thus can be used effectively to characterize cold hardiness in Juglans species. We conclude that the genetic diversity within local J. regia populations is low and additional germplasm is needed for development of more regionally adapted J. regia varieties.

Evidence for widespread selection in shaping the genomic landscape during speciation of Populus

Increasing our understanding of how evolutionary processes drive the genomic landscape of variation is fundamental to a better understanding of the genomic consequences of speciation. However, genome-wide patterns of within- and between- species variation have not been fully investigated in most forest tree species despite their global ecological and economic importance. Here, we use whole-genome resequencing data from four Populus species spanning the speciation continuum to reconstruct their demographic histories and investigate patterns of diversity and divergence within and between species. Using Populus trichocarpa as an outgroup species, we further infer the genealogical relationships and estimate the extent of ancient introgression among the three aspen species (Populus tremula, Populus davidiana and Populus tremuloides) throughout the genome. Our results show substantial variation in these patterns along the genomes with this variation being strongly predicted by local recombination rates and the density of functional elements. This implies that the interaction between recurrent selection and intrinsic genomic features has dramatically sculpted the genomic landscape over long periods of time. In addition, our findings provide evidence that, apart from background selection, recent positive selection and long-term balancing selection have also been crucial components in shaping patterns of genome-wide variation during the speciation process.

Development of 32 novel microsatellite loci in Juglans sigillata using genomic data

PREMISE

A novel set of microsatellite markers was developed for Juglans sigillata (Juglandaceae), an endemic walnut species in southwestern China, to facilitate cultivar identification and future investigations into the genetic structure and domestication history of this species and its close relatives.

METHODS AND RESULTS

We developed 32 microsatellite loci for J. sigillata using genomic data and used them to examine 60 individuals from three natural populations. A high level of polymorphism was detected by these primers, with up to eight alleles observed per locus, and an average of four alleles across populations. The levels of observed and expected heterozygosity ranged from 0.000-1.000 and 0.000-0.785, respectively. All but two of the loci were also successfully amplified in three closely related Eurasian Juglans species (J. regia, J. cathayensis, and J. mandshurica).

CONCLUSIONS

The microsatellite loci identified here provide a powerful resource for examining the genetic structure and domestication history of Juglans, as well as identification of its cultivars.