Uncovering ghost introgression through genomic analysis of a distinct eastern Asian hickory species.

Ghost introgression, or the transfer of genetic material from extinct or unsampled lineages to sampled species, has attracted much attention. However, conclusive evidence for ghost introgression, especially in plant species, remains scarce. Here, we newly assembled chromosome-level genomes for both Carya sinensis and Carya cathayensis, and additionally re-sequenced the whole genomes of 43 C. sinensis individuals as well as 11 individuals representing 11 diploid hickory species. These genomic datasets were used to investigate the reticulation and bifurcation patterns within the genus Carya (Juglandaceae), with a particular focus on the beaked hickory C. sinensis. By combining the D-statistic and BPP methods, we obtained compelling evidence that supports the occurrence of ghost introgression in C. sinensis from an extinct ancestral hickory lineage. This conclusion was reinforced through the phylogenetic network analysis and a genome scan method VolcanoFinder, the latter of which can detect signatures of adaptive introgression from unknown donors. Our results not only dispel certain misconceptions about the phylogenetic history of C. sinensis but also further refine our understanding of Carya's biogeography via divergence estimates. Moreover, the successful integration of the D-statistic and BPP methods demonstrates their efficacy in facilitating a more precise identification of introgression types.

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, P. 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.

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.

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.

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.

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.

Do annual and perennial populations of an insect-pollinated plant species differ in mating system?

BACKGROUND AND AIMS: Theory predicts that outcrossing should be more prevalent among perennials than annuals, a pattern confirmed by comparative evidence from diverse angiosperm families. However, intraspecific comparisons between annual and perennial populations are few because such variation is uncommon among flowering plants. Here, we test the hypothesis that perennial populations outcross more than annual populations by investigating Incarvillea sinensis, a wide-ranging insect-pollinated herb native to China. The occurrence of both allopatric and sympatric populations allows us to examine the stability of mating system differences between life histories under varying ecological conditions. METHODS: We estimated outcrossing rates and biparental inbreeding in 16 allopatric and five sympatric populations in which both life histories coexisted using 20 microsatellite loci. In each population we measured height, branch number, corolla size, tube length and herkogamy for ~30 individuals. In a sympatric population, we recorded daily flower number, pollinator visitation and the fruit and seed set of annual and perennial plants. KEY RESULTS: As predicted, outcrossing rates (t) were considerably higher in perennial (mean = 0.76) than annual (mean = 0.09) populations. This difference in mating system was also maintained at sympatric sites where plants grew intermixed. In both allopatric and sympatric populations the degree of herkogamy was consistently larger in outcrossing than selfing plants. Perennials were more branched, with more and larger flowers than in annuals. In a sympatric population, annuals had a significantly higher fruit and seed set than perennials. CONCLUSIONS: Genetically based differences in herkogamy between annuals and perennials appear to play a key role in governing outcrossing rates in populations, regardless of variation in local ecological conditions. The maintenance of mating system and life history trait differentiation between perennial and annual populations of I. sinensis probably results from correlated evolution in response to local environmental conditions.

Phylogenomics Reveals an Ancient Hybrid Origin of the Persian Walnut.

Persian walnut (Juglans regia) is cultivated worldwide for its high-quality wood and nuts, but its origin has remained mysterious because in phylogenies it occupies an unresolved position between American black walnuts and Asian butternuts. Equally unclear is the origin of the only American butternut, J. cinerea. We resequenced the whole genome of 80 individuals from 19 of the 22 species of Juglans and assembled the genome of its relatives Pterocarya stenoptera and Platycarya strobilacea. Using phylogenetic-network analysis of single-copy nuclear genes, genome-wide site pattern probabilities, and Approximate Bayesian Computation, we discovered that J. regia (and its landrace J. sigillata) arose as a hybrid between the American and the Asian lineages and that J. cinerea resulted from massive introgression from an immigrating Asian butternut into the genome of an American black walnut. Approximate Bayesian Computation modeling placed the hybrid origin in the late Pliocene, ∼3.45 My, with both parental lineages since having gone extinct in Europe.

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

Whether species demography and diversification are driven primarily by extrinsic environmental changes such as climatic oscillations in the Quaternary or by intrinsic biological interactions like coevolution between antagonists is a matter of active debate. In fact, their relative importance can be assessed by tracking past population fluctuations over considerable time periods. We applied the pairwise sequentially Markovian coalescent approach on the genomes of 11 temperate Juglans species to estimate trajectories of changes in effective population size (Ne ) and used a Bayesian-coalescent based approach that simultaneously considers multiple genomes (G-PhoCS) to estimate divergence times between lineages. Ne curves of all study species converged 1.0 million yr ago, probably reflecting the time when the walnut genus last shared a common ancestor. This estimate was confirmed by the G-PhoCS estimates of divergence times. But all species did not react similarly to the dramatic climatic oscillations following early Pleistocene cooling, so the timing and amplitude of changes in Ne differed among species and even among conspecific lineages. The population histories of temperate walnut species were not driven by extrinsic environmental changes alone, and a key role was probably played by species-specific factors such as coevolutionary interactions with specialized pathogens.

Phylogeographic breaks within Asian butternuts indicate the existence of a phytogeographic divide in East Asia.

East Asia has been hypothesized to be subdivided into two distinct northern and southern areas, separated by a band of dry climate that was far more severe in the early Tertiary but still exists today. However, this biogeographic hypothesis has rarely been tested using a molecular phylogeographic approach. We genotyped 70 populations throughout the distributional range of Asian butternuts (Juglans section Cardiocaryon) using eight chloroplast DNA regions, one single-copy nuclear gene, and 17 nuclear microsatellite loci, supplemented with paleodistribution modeling of the major genetic clades. The genetic data consistently identified two clades, one northern, comprising Juglans mandshurica and Juglans ailantifolia, and one southern, comprising Juglans cathayensis. The two clades diverged through climate-induced vicariance of an ancestral northern range during the mid-Miocene and remained mostly separate thereafter, with geographical isolation of the Japanese Islands and refugial isolation or secondary contacts in the late Pleistocene producing further subdivision within the northern clade. But beyond all that, we also discovered a role of environmental adaptation in maintaining and/or reinforcing the north-south divergence. Asian butternuts offer a strong case for the existence of a biogeographic divide between the northern and southern parts of East Asia during the Neogene and into the Pleistocene.

Phylogeography of postglacial range expansion in Juglans mandshurica (Juglandaceae) reveals no evidence of bottleneck, loss of genetic diversity, or isolation by distance in the leading-edge populations.

The past studies of postglacial recolonization patterns in high latitude regions have revealed a significant role of dispersal capacity in shaping the genetic diversity and population structure of temperate trees. However, most of these studies have focused on species with long-distance dispersal followed by exponential population growth and were therefore unable to reveal the patterns in the case of a gradual expansion. Here we studied the impacts of postglacial range expansions on the distribution of genetic diversity in the Manchurian walnut (Juglans mandshurica), a common tree of East Asian cool-temperate deciduous forests that apparently lacks long-distance seed dispersal ability. The genetic diversity and structure of 19 natural walnut populations in Northeast China and the Korean Peninsula were examined using 17 nuclear simple sequence repeat (SSR) loci. Potential habitats under current and past climatic conditions were predicted using the ecological niche modelling (ENM) method. Bayesian clustering analysis revealed three groups, which were inferred to have diverged through multiple glacial-interglacial cycles in multiple refugia during the Quaternary Period. ENM estimated a southward range shift at the LGM, but high suitability scores still occurred in the western parts of the Changbai Mountains (Northeast China), the Korean peninsula and the exposed seafloor of the Yellow Sea. In contrast to most other cool-temperate trees co-occurring in the same region, the Manchurian walnut did not show any evidence of a population bottleneck, loss of genetic diversity or isolation by distance during the postglacial expansion. Our study clearly indicates that current northern populations originated from one glacial lineage and recolonization via a gradually advancing front due to the lack of a long-distance seed dispersal mechanism led to no latitudinal decrease in genetic diversity.

Seasonal variation in the mating system of a selfing annual with large floral displays.

BACKGROUND AND AIMS: Flowering plants display considerable variation in mating system, specifically the relative frequency of cross- and self-fertilization. The majority of estimates of outcrossing rate do not account for temporal variation, particularly during the flowering season. Here, we investigated seasonal variation in mating and fertility in Incarvillea sinensis (Bignoniaceae), an annual with showy, insect-pollinated, 'one-day' flowers capable of delayed selfing. We examined the influence of several biotic and abiotic environmental factors on day-to-day variation in fruit set, seed set and patterns of mating. METHODS: We recorded daily flower number and pollinator abundance in nine 3 × 3-m patches in a population at Mu Us Sand land, Inner Mongolia, China. From marked flowers we collected data on daily fruit and seed set and estimated outcrossing rate and biparental inbreeding using six microsatellite loci and 172 open-pollinated families throughout the flowering period. KEY RESULTS: Flower density increased significantly over most of the 50-d flowering season, but was associated with a decline in levels of pollinator service by bees, particularly on windy days. Fruit and seed set declined over time, especially during the latter third of the flowering period. Multilocus estimates of outcrossing rate were obtained using two methods (the programs MLTR and BORICE) and both indicated high selfing rates of ∼80 %. There was evidence for a significant increase in levels of selfing as the flowering season progressed and pollinator visitation declined. Biparental inbreeding also declined significantly as the flowering season progressed. CONCLUSIONS: Temporal variation in outcrossing rates may be a common feature of the mating biology of annual, insect-pollinated plants of harsh environments but our study is the first to examine seasonal mating-system dynamics in this context. Despite having large flowers and showy floral displays, I. sinensis attracted relatively few pollinators. Delayed selfing by corolla dragging largely explains the occurrence of mixed mating in I. sinensis, and this mode of self-fertilization probably functions to promote reproductive assurance when pollinator service is limited by windy environmental conditions and competition from co-occurring flowering plants.

Contrasts between the phylogeographic patterns of chloroplast and nuclear DNA highlight a role for pollen-mediated gene flow in preventing population divergence in an East Asian temperate tree.

Plant phylogeographic studies in East Asia have provided support for the biogeographic hypothesis that the complex landforms and climate of this region have provided substantial opportunities for allopatric speciation. However, most of these studies have been based on maternally inherited chloroplast DNA (cpDNA) markers and were therefore unable to reveal the role of pollen-mediated gene flow in preventing population divergence. Here, we investigate the phylogeography of the Chinese walnut Juglans cathayensis, a temperate deciduous tree widely distributed across disjunct montane sites in subtropical China. We genotyped 19 populations using seven cpDNA fragments and ten nuclear microsatellite loci and modeled the ecological niche of J. cathayensis. CpDNA analysis identified a total of nine haplotypes, and each of the 19 sampled populations was fixed for a single haplotype, displaying a prominent phylogeographic structure. The results of ecological niche modeling indicated that J. cathayensis populations survived the last glaciation in situ, although they were probably more fragmented than today. In contrast, we detected a much weaker, but nonetheless clear, genetic structure based on nuclear microsatellite data. Our study demonstrates how extensive pollen flow can erase the genetic imprint of long-term refugial isolation in maternal lineages, effectively preventing population differentiation in temperate, particularly wind-pollinated, forest trees in subtropical China.

ploidyfrost: Reference-free estimation of ploidy level from whole genome sequencing data based on de Bruijn graphs.

Polyploidy is ubiquitous and its consequences are complex and variable. A change of ploidy level generally influences genetic diversity and results in morphological, physiological and ecological differences between cells or organisms with different ploidy levels. To avoid cumbersome experiments and take advantage of the less biased information provided by the vast amounts of genome sequencing data, computational tools for ploidy estimation are urgently needed. Until now, although a few such tools have been developed, many aspects of this estimation, such as the requirement of a reference genome, the lack of informative results and objective inferences, and the influence of false positives from errors and repeats, need further improvement. We have developed ploidyfrost, a de Bruijn graph-based method, to estimate ploidy levels from whole genome sequencing data sets without a reference genome. ploidyfrost provides a visual representation of allele frequency distribution generated using the ggplot2 package as well as quantitative results using the Gaussian mixture model. In addition, it takes advantage of colouring information encoded in coloured de Bruijn graphs to analyse multiple samples simultaneously and to flexibly filter putative false positives. We evaluated the performance of ploidyfrost by analysing highly heterozygous or repetitive samples of Cyclocarya paliurus and a complex allooctoploid sample of Fragaria × ananassa. Moreover, we demonstrated that the accuracy of analysis results can be improved by constraining a threshold such as Cramér's V coefficient on variant features, which may significantly reduce the side effects of sequencing errors and annoying repeats on the graphical structure constructed.

Genomic Analysis of Plastid-Nuclear Interactions and Differential Evolution Rates in Coevolved Genes across Juglandaceae Species.

The interaction between the nuclear and chloroplast genomes in plants is crucial for preserving essential cellular functions in the face of varying rates of mutation, levels of selection, and modes of transmission. Despite this, identifying nuclear genes that coevolve with chloroplast genomes at a genome-wide level has remained a challenge. In this study, we conducted an evolutionary rate covariation analysis to identify candidate nuclear genes coevolving with chloroplast genomes in Juglandaceae. Our analysis was based on 4,894 orthologous nuclear genes and 76 genes across seven chloroplast partitions in nine Juglandaceae species. Our results indicated that 1,369 (27.97%) of the nuclear genes demonstrated signatures of coevolution, with the Ycf1/2 partition yielding the largest number of hits (765) and the ClpP1 partition yielding the fewest (13). These hits were found to be significantly enriched in biological processes related to leaf development, photoperiodism, and response to abiotic stress. Among the seven partitions, AccD, ClpP1, MatK, and RNA polymerase partitions and their respective hits exhibited a narrow range, characterized by dN/dS values below 1. In contrast, the Ribosomal, Photosynthesis, Ycf1/2 partitions and their corresponding hits, displayed a broader range of dN/dS values, with certain values exceeding 1. Our findings highlight the differences in the number of candidate nuclear genes coevolving with the seven chloroplast partitions in Juglandaceae species and the correlation between the evolution rates of these genes and their corresponding chloroplast partitions.

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.

Influence of asymmetrical mating patterns and male reproductive success on the maintenance of sexual polymorphism in Acer pictum subsp. mono (Aceraceae).

Populations of Acer species often contain more than three sex phenotypes with complex sexual polymorphism including duodichogamy, protandry and protogyny. We identified the mechanisms that maintain sexual polymorphism in Acer pictum subsp. mono, a temperate tree from northern China, by investigating maternal mating patterns and male reproductive success. We used paternity analyses to estimate rates of outcrossing and disassortative mating, as well as male outcrossed siring success, in a population of A. pictum subsp. mono with uneven sex phenotype ratios (duodichogamous 69.1%, protandrous 19.6%, protogynous 11.3%). We used a pollen-transfer model to investigate whether the unequal ratios of sex phenotypes could be explained by the observed patterns of mating. Most progeny resulted from outcrossing, particularly disassortative among the sex phenotypes. Although the duodichogamous phenotype showed a significant amount of intraphenotypic mating, the frequency did not exceed that of disassortative mating. We detected no significant differences in male outcrossed siring success among the sex phenotypes. The pollen-transfer model demonstrated that sex phenotype ratios could be maintained by the observed mating pattern in the population. Our results indicate that disassortative mating among the sex phenotypes can maintain sexual polymorphism in A. pictum subsp. mono and that ratios biased towards duodichogamy can result from frequent intraphenotypic mating in this phenotype.

The rediscovery of Caryapoilanei (Juglandaceae) after 63 years reveals a new record from China.

Despite having been first published in 1941, Caryapoilanei (A.Chev.) J.-F.Leroy is only known from three collections in Vietnam, Laos and Thailand. It has not been recollected since then and was long suspected to have become extinct through repeated deforestation events. Here, we report the rediscovery, and meanwhile the first new record in China, of this extremely rare gigantic hickory species at Yunnan province 63 years after its last collection in 1958. Three small patchy subpopulations were found with a total of about 50 adult trees having diameter at breast height (DBH) larger than 60 cm, together with some seedlings and saplings, but the fruit set was low. Based on new and fresh material, we present a revised morphological description of C.poilanei, and an updated distribution map for the species. In addition, we also provide a key for the hickories in China. Lastly, we suggest C.poilanei should be listed as a Critically Endangered (CR) species according to the latest IUCN Red List Categories and Criteria.

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.