Shared xerophytic genes and their re-use in local adaptation to aridity in the desert plant Gymnocarpos przewalskii.

In order to thrive and survive, plant species need to combine stability in the long term and rapid response to environmental challenges in the short term. The former would be reflected by parallel or convergent adaptation across species, and the latter by pronounced local adaptation among populations of the same species. In the present study, we generated a high-quality genome and re-sequenced 177 individuals for Gymnocarpos przewalskii, an important desert plant species from North-West China, to detect local adaptation. We first focus on ancient adaptation to aridity at the molecular level by comparing the genomic data of 15 species that vary in their ability to withstand aridity. We found that a total of 118 genes were shared across xerophytic species but absent from non-xerophytic species. Of the 65 found in G. przewalskii, 63 were under purifying selection and two under positive selection. We then focused on local adaptation. Up to 20% of the G. przewalskii genome showed signatures of local adaptation to aridity during population divergence. Thirteen of the selected shared xerophytic genes were reused in local adaptation after population differentiation. Hence, only about 20% of the genes shared and specific to xerophytic species and associated with adaptation to aridity were later recruited for local adaptation in G. przewalskii.

Low-coverage whole genome sequencing of diverse Dioscorea bulbifera accessions for plastome resource development, polymorphic nuclear SSR identification, and phylogenetic analyses.

Dioscorea bulbifera (Dioscoreaceae), a versatile herbaceous climber native to Africa and Asia, holds significant nutritional and medicinal value. Despite extensive characterization and genetic variability analyses of African accessions, studies on the genetic variation of this species in China are limited. To address this gap, we conducted low-coverage whole genome sequencing on D. bulbifera accessions from diverse regions across mainland China and Taiwan island. Our initial investigation encompassed comprehensive comparative plastome analyses of these D. bulbifera accessions, and developing plastome resources (including plastome-derived repetitive sequences, SSRs, and divergent hotspots). We also explored polymorphic nuclear SSRs and elucidated the intraspecific phylogeny of these accessions. Comparative plastome analyses revealed that D. bulbifera plastomes exhibited a conserved quadripartite structure with minimal size variation mainly attributed to intergenic spacer regions, reinforcing prior observations of a high degree of conservation within a species. We identified 46 to 52 dispersed repeats and 151 to 163 plastome-derived SSRs, as well as highlighted eight key divergent hotspots in these D. bulbifera accessions. Furthermore, we developed 2731 high-quality candidate polymorphic nuclear SSRs for D. bulbifera. Intraspecific phylogenetic analysis revealed three distinct clades, where accessions from Southeast China formed a sister group to those from South China and Taiwan island, and collectively, these two clades formed a sister group to the remaining accessions, indicating potential regional genetic divergence. These findings not only contributed to the understanding of the genetic variation of D. bulbifera, but also offered valuable resources for future research, breeding efforts, and utilization of this economically important plant species.

Unraveling the Molecular Basis of Color Variation in Dioscorea alata Tubers: Integrated Transcriptome and Metabolomics Analysis.

Dioscorea alata L. (Dioscoreaceae) is a widely cultivated tuber crop with variations in tuber color, offering potential value as health-promoting foods. This study focused on the comparison of D. alata tubers possessing two distinct colors, white and purple, to explore the underlying mechanisms of color variation. Flavonoids, a group of polyphenols known to influence plant color and exhibit antioxidant properties, were of particular interest. The total phenol and total flavonoid analyses revealed that purple tubers (PTs) have a significantly higher content of these metabolites than white tubers (WTs) and a higher antioxidant activity than WTs, suggesting potential health benefits of PT D. alata. The transcriptome analysis identified 108 differentially expressed genes associated with the flavonoid synthesis pathway, with 57 genes up-regulated in PTs, including CHS, CHI, DFR, FLS, F3H, F3'5'H, LAR, ANS, and ANR. The metabolomics analysis demonstrated that 424 metabolites, including 104 flavonoids and 8 tannins, accumulated differentially in PTs and WTs. Notably, five of the top ten up-regulated metabolites were flavonoids, including 6-hydroxykaempferol-7-O-glucoside, pinocembrin-7-O-(6″-O-malonyl)glucoside, 6-hydroxykaempferol-3,7,6-O-triglycoside, 6-hydroxykaempferol-7-O-triglycoside, and cyanidin-3-O-(6″-O-feruloyl)sophoroside-5-O-glucoside, with the latter being a precursor to anthocyanin synthesis. Integrating transcriptome and metabolomics data revealed that the 57 genes regulated 20 metabolites within the flavonoid synthesis pathway, potentially influencing the tubers' color variation. The high polyphenol content and antioxidant activity of PTs indicate their suitability as nutritious and health-promoting food sources. Taken together, the findings of this study provide insights into the molecular basis of tuber color variation in D. alata and underscore the potential applications of purple tubers in the food industry and human health promotion. The findings contribute to the understanding of flavonoid biosynthesis and pigment accumulation in D. alata tubers, opening avenues for future research on enhancing the nutritional quality of D. alata cultivars.

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

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

Pangenome of water caltrop reveals structural variations and asymmetric subgenome divergence after allopolyploidization.

Water caltrop (Trapa spp., Lythraceae) is a traditional but currently underutilized non-cereal crop. Here, we generated chromosome-level genome assemblies for the two diploid progenitors of allotetraploid Trapa. natans (4x, AABB), i.e., diploid T. natans (2x, AA) and Trapa incisa (2x, BB). In conjunction with four published (sub)genomes of Trapa, we used gene-based and graph-based pangenomic approaches and a pangenomic transposable element (TE) library to develop Trapa genomic resources. The pangenome displayed substantial gene-content variation with dispensable and private gene clusters occupying a large proportion (51.95%) of the total cluster sets in the six (sub)genomes. Genotyping of presence-absence variation (PAVs) identified 40 453 PAVs associated with 2570 genes specific to A- or B-lineages, of which 1428 were differentially expressed, and were enriched in organ development process, organic substance metabolic process and response to stimulus. Comparative genome analyses showed that the allotetraploid T. natans underwent asymmetric subgenome divergence, with the B-subgenome being more dominant than the A-subgenome. Multiple factors, including PAVs, asymmetrical amplification of TEs, homeologous exchanges (HEs), and homeolog expression divergence, together affected genome evolution after polyploidization. Overall, this study sheds lights on the genome architecture and evolution of Trapa, and facilitates its functional genomic studies and breeding program.

Intraspecific phylogeny and genomic resources development for an important medical plant Dioscorea nipponica, based on low-coverage whole genome sequencing data.

Dioscorea nipponica Makino, a perennial twining herb with medicinal importance, has a disjunctive distribution in the Sino-Japanese Floristic Region. It has a long history in traditional Chinese medicine, with demonstrated efficacy against various health conditions. However, the limited genomic data and knowledge of genetic variation have hindered its comprehensive exploration, utilization and conservation. In this study, we undertook low-coverage whole genome sequencing of diverse D. nipponica accessions to develop both plastome (including whole plastome sequences, plastome-derived SSRs and plastome-divergent hotspots) and nuclear genomic resources (including polymorphic nuclear SSRs and single-copy nuclear genes), as well as elucidate the intraspecific phylogeny of this species. Our research revealed 639 plastome-derived SSRs and highlighted six key mutational hotspots (namely CDS ycf1, IGS trnL-rpl32, IGS trnE-trnT, IGS rps16-trnQ, Intron 1 of clpP, and Intron trnG) within these accessions. Besides, three IGS regions (i.e., ndhD-cssA, trnL-rpl32, trnD-trnY), and the intron rps16 were identified as potential markers for distinguishing D. nipponica from its closely related species. In parallel, we successfully developed 988 high-quality candidate polymorphic nuclear SSRs and identified 17 single-copy nuclear genes for D. nipponica, all of which empower us to conduct in-depth investigations into phylogenetics and population genetics of this species. Although our phylogenetic analyses, based on plastome sequences and single-copy nuclear genes revealed cytonuclear discordance within D. nipponica, both findings challenged the current subspecies classification. In summary, this study developed a wealth of genomic resources for D. nipponica and enhanced our understanding of the intraspecific phylogeny of this species, offering valuable insights that can be instrumental in the conservation and strategic utilization of this economically significant plant.

Low-coverage whole genome sequencing of eleven species/subspecies in Dioscorea sect. Stenophora (Dioscoreaceae): comparative plastome analyses, molecular markers development and phylogenetic inference.

Dioscorea sect. Stenophora (Dioscoreaceae) comprises about 30 species that are distributed in the temperate and subtropical regions of the Northern Hemisphere. Despite being evolutionarily "primitive" and medically valuable, genomic resources and molecular studies of this section are still scarce. Here, we conducted low-coverage whole genome sequencing of 11 Stenophora species/subspecies to retrieve their plastome information (whole plastome characteristics, plastome-divergent hotspots, plastome-derived SSRs, etc.) and polymorphic nuclear SSRs, as well as performed comparative plastome and phylogenetic analyses within this section. The plastomes of Stenophora species/subspecies ranged from 153,691 bp (D. zingiberensis) to 154,149 bp (D. biformifolia) in length, and they all contained the same 114 unique genes. All these plastomes were highly conserved in gene structure, gene order and GC content, although variations at the IR/SC borders contributed to the whole length differences among them. The number of plastome-derived SSRs among Stenophora species/subspecies varied from 74 (D. futschauensis) to 93 (D. zingiberensis), with A/T found to be the most frequent one. Seven highly variable regions and 12 polymorphic nuclear SSRs were identified in this section, thereby providing important information for further taxonomical, phylogenetic and population genetic studies. Phylogenomic analyses based on whole plastome sequences and 80 common protein coding genes strongly supported D. biformifolia and D. banzhuana constituted the successive sister species to the remaining sampled species, which could be furtherly divided into three clades. Overall, this study provided a new perspective for plastome evolution of Stenophora, and proved the role of plastome phylogenomic in improving the phylogenetic resolution in this section. These results also provided an important reference for the protection and utilization of this economically important section.

Plastid Genome Assembly Using Long-read data.

Although plastid genome (plastome) structure is highly conserved across most seed plants, investigations during the past two decades have revealed several disparately related lineages that experienced substantial rearrangements. Most plastomes contain a large inverted repeat and two single-copy regions, and a few dispersed repeats; however, the plastomes of some taxa harbour long repeat sequences (>300 bp). These long repeats make it challenging to assemble complete plastomes using short-read data, leading to misassemblies and consensus sequences with spurious rearrangements. Single-molecule, long-read sequencing has the potential to overcome these challenges, yet there is no consensus on the most effective method for accurately assembling plastomes using long-read data. We generated a pipeline, plastid Genome Assembly Using Long-read data (ptGAUL), to address the problem of plastome assembly using long-read data from Oxford Nanopore Technologies (ONT) or Pacific Biosciences platforms. We demonstrated the efficacy of the ptGAUL pipeline using 16 published long-read data sets. We showed that ptGAUL quickly produces accurate and unbiased assemblies using only ~50× coverage of plastome data. Additionally, we deployed ptGAUL to assemble four new Juncus (Juncaceae) plastomes using ONT long reads. Our results revealed many long repeats and rearrangements in Juncus plastomes compared with basal lineages of Poales. The ptGAUL pipeline is available on GitHub: https://github.com/Bean061/ptgaul.

Pan-Plastome of Greater Yam (Dioscorea alata) in China: Intraspecific Genetic Variation, Comparative Genomics, and Phylogenetic Analyses.

Dioscorea alata L. (Dioscoreaceae), commonly known as greater yam, water yam, or winged yam, is a popular tuber vegetable/food crop worldwide, with nutritional, health, and economical importance. China is an important domestication center of D. alata, and hundreds of cultivars (accessions) have been established. However, genetic variations among Chinese accessions remain ambiguous, and genomic resources currently available for the molecular breeding of this species in China are very scarce. In this study, we generated the first pan-plastome of D. alata, based on 44 Chinese accessions and 8 African accessions, and investigated the genetic variations, plastome evolution, and phylogenetic relationships within D. alata and among members of the section Enantiophyllum. The D. alata pan-plastome encoded 113 unique genes and ranged in size from 153,114 to 153,161 bp. A total of four whole-plastome haplotypes (Haps I-IV) were identified in the Chinese accessions, showing no geographical differentiation, while all eight African accessions shared the same whole-plastome haplotype (Hap I). Comparative genomic analyses revealed that all four whole plastome haplotypes harbored identical GC content, gene content, gene order, and IR/SC boundary structures, which were also highly congruent with other species of Enantiophyllum. In addition, four highly divergent regions, i.e., trnC-petN, trnL-rpl32, ndhD-ccsA, and exon 3 of clpP, were identified as potential DNA barcodes. Phylogenetic analyses clearly separated all the D. alata accessions into four distinct clades corresponding to the four haplotypes, and strongly supported that D. alata was more closely related to D. brevipetiolata and D. glabra than D. cirrhosa, D. japonica, and D. polystachya. Overall, these results not only revealed the genetic variations among Chinese D. alata accessions, but also provided the necessary groundwork for molecular-assisted breeding and industrial utilization of this species.

Genome Sequence Resource for Colletotrichum gloeosporioides, an Important Pathogenic Fungus Causing Anthracnose of Dioscorea alata.

Anthracnose disease is one of the most important diseases of Dioscorea alata and many other food yams, which is caused by Colletotrichum gloeosporioides fungus from the Glomerellaceae family of the Sordariomycetes class. In the present study, a C. gloeosporioides starin named CgDa01 was isolated from D. alata, and its genome was sequenced based on Oxford Nanopore technology (ONT) and the Illumina sequencing platform. The high-quality genome of CgDa01 was assembled with a 62.78 Mb genome size and 15,845 predicted protein-coding genes. The proteins of predicted genes were annotated using multiple public databases, including the nonredundant protein database, the InterProScan databases, and Kyoto Encyclopedia of Genes and Genomes. Among the annotated protein-coding genes, 55 were predicted as potential virulence genes by the fungal virulence factor database. The C. gloeosporioides CgDa01 genome assembly described in this study can serve as a resource for better understanding the pathogenic mechanism of C. gloeosporioides on yam hosts.

Morphological Characteristics and Comparative Chloroplast Genome Analyses between Red and White Flower Phenotypes of Pyracantha fortuneana (Maxim.) Li (Rosaceae), with Implications for Taxonomy and Phylogeny.

Pyracantha fortuneana (Maxim.) Li (Rosaceae), commonly known as Chinese firethorn, is an evergreen shrub with high nutritional, medicinal, and horticultural importance. This species typically has white flowers, but a rare red flower phenotype has been found in very few wild populations in western Hubei, China, showing great ornamental potential. In this study, the complete chloroplast genome of the red flower phenotype of P. fortuneana was reported for the first time, using high-throughput sequencing technology. The complete chloroplast genome was 160,361 bp in length and showed a typical quadripartite structure with a pair of inverted repeat (IR) regions (26,350 bp) separated by a large single-copy (LSC) region (88,316 bp) and a small single-copy (SSC) region (19,345 bp). A total of 131 functional genes were annotated in this chloroplast genome, including 86 protein-coding genes (PCGs), eight rRNA genes, and 37 tRNA genes. Comparative chloroplast genome analyses revealed that high genome similarity existed not only between red and white flower phenotypes of P. fortuneana, but also among Pyracantha species. No evidence for positive selection was found in any PCG, suggesting the evolutionary conservation of Pyracantha chloroplast genomes. Furthermore, four mutational hotspots (trnG-trnR-atpA, psbZ-trnG-trnfM-rps14, ycf3-trnS-rps4, and ndhF-rpl32) with π > 0.004 were identified as potential molecular markers for Pyracantha species. Phylogenomic analysis strongly supported that the red flower phenotype of P. fortuneana was nested within the common white flower phenotype. Based on both morphological and molecular evidence, we suggest that the red flower phenotype of P. fortuneana could be considered as a new forma. Overall, the availability of these genetic resources will not only offer valuable information for further studies on molecular taxonomy, phylogeny, and population genetics of Pyracantha species but also could be used as potential genetic resources for Chinese firethorn breeding.

Molecular evolution of the Pi-d2 gene conferring resistance to rice blast in Oryza.

The exploitation of plant disease resistance (R) genes in breeding programs is an effective strategy for coping with pathogens. An understanding of R gene variation is the basis for this strategy. Rice blast disease, caused by the Magnaporthe oryzae fungus, is a destructive disease of rice. The rice blast resistance gene Pi-d2 represents a new class of plant R gene because of its novel extracellular domain. We investigated the nucleotide polymorphism, phylogenetic topology and evolution patterns of the Pi-d2 gene among 67 cultivated and wild rice relatives. The Pi-d2 gene originated early in the basal Poales and has remained as a single gene without expansion. The striking finding is that susceptible Pi-d2 alleles might be derived from a single nucleotide substitution of the resistant alleles after the split of Oryza subspecies. Functional pleiotropy and linkage effects are proposed for the evolution and retention of the disease-susceptible alleles in rice populations. One set of DNA primers was developed from the polymorphic position to detect the functional nucleotide polymorphism for disease resistance of the Pi-d2 gene based on conventional Polymerase Chain Reaction. The nucleotide diversity level varied between different domains of the Pi-d2 gene, which might be related to distinct functions of each domain in the disease defense response. Directional (or purifying) selection appears dominant in the molecular evolution of the Pi-d2 gene and has shaped its conserved variation pattern.

Chromosome-Level Genome Assembly of a Fragrant Japonica Rice Cultivar 'Changxianggeng 1813' Provides Insights into Genomic Variations between Fragrant and Non-Fragrant Japonica Rice.

East Asia has an abundant resource of fragrant japonica rice that is gaining increasing interest among both consumers and producers. However, genomic resources and in particular complete genome sequences currently available for the breeding of fragrant japonica rice are still scarce. Here, integrating Nanopore long-read sequencing, Illumina short-read sequencing, and Hi-C methods, we presented a high-quality chromosome-level genome assembly (~378.78 Mb) for a new fragrant japonica cultivar 'Changxianggeng 1813', with 31,671 predicated protein-coding genes. Based on the annotated genome sequence, we demonstrated that it was the badh2-E2 type of deletion (a 7-bp deletion in the second exon) that caused fragrance in 'Changxianggeng 1813'. Comparative genomic analyses revealed that multiple gene families involved in the abiotic stress response were expanded in the 'Changxianggeng 1813' genome, which further supported the previous finding that no generalized loss of abiotic stress tolerance associated with the fragrance phenotype. Although the 'Changxianggeng 1813' genome showed high genomic synteny with the genome of the non-fragrant japonica rice cultivar Nipponbare, a total of 289,970 single nucleotide polymorphisms (SNPs), 96,093 small insertion-deletion polymorphisms (InDels), and 8690 large structure variants (SVs, >1000 bp) were identified between them. Together, these genomic resources will be valuable for elucidating the mechanisms underlying economically important traits and have wide-ranging implications for genomics-assisted breeding in fragrant japonica rice.

A Comparative Study of Genetic Responses to Short- and Long-Term Habitat Fragmentation in a Distylous Herb Hedyotis chyrsotricha (Rubiaceae).

The genetic effects of habitat fragmentation are complex and are influenced by both species traits and landscape features. For plants with strong seed or pollen dispersal capabilities, the question of whether the genetic erosion of an isolated population becomes stronger or is counterbalanced by sufficient gene flow across landscapes as the timescales of fragmentation increase has been less studied. In this study, we compared the population structure and genetic diversity of a distylous herb, Hedyotis chyrsotricha (Rubiaceae), in two contrasting island systems of southeast China. Based on RAD-Seq data, our results showed that populations from the artificially created Thousand-Island Lake (TIL) harbored significantly higher levels of genetic diversity than those from the Holocene-dated Zhoushan Archipelago (ZA) (π = 0.247 vs. 0.208, HO = 0.307 vs. 0.256, HE = 0.228 vs. 0.190), while genetic differences between island and mainland populations were significant in neither the TIL region nor the ZA region. A certain level of population substructure was found in TIL populations, and the level of gene flow among TIL populations was also lower than in ZA populations (m = 0.019 vs. 0.027). Overall, our comparative study revealed that genetic erosion has not become much stronger for the island populations of either the TIL or ZA regions. Our results emphasized that the matrix of water in the island system may facilitate the seed (fruit) dispersal of H. chrysotricha, thus maintaining population connectivity and providing ongoing resilience to the effects of habitat fragmentation over thousands of years.

Genome-wide analyses of introgression between two sympatric Asian oak species.

Introgression can be an important source of new alleles for adaption under rapidly changing environments, perhaps even more important than standing variation. Though introgression has been extensively studied in many plants and animals, key questions on the underlying mechanisms of introgression still remain unanswered. In particular, we are yet to determine the genomic distribution of introgressed regions along the genome; whether the extent and patterns of introgression are influenced by ecological factors; and when and how introgression contributes to adaptation. Here, we generated high-quality genomic resources for two sympatric widespread Asian oak species, Quercus acutissima and Q. variabilis, sampled in multiple forests to study introgression between them. We show that introgressed regions are broadly distributed across the genome. Introgression was affected by genetic divergence between pairs of populations and by the similarity of the environments in which they live-populations occupying similar ecological sites tended to share the same introgressed regions. Introgressed genomic footprints of adaptation were preferentially located in regions with suppressed recombination rate. Introgression probably confers adaptation in these oak populations by introducing allelic variation in cis-regulatory elements, in particular through transposable element insertions, thereby altering the regulation of genes related to stress. Our results provide new avenues of research for uncovering mechanisms of adaptation due to hybridization in sympatric species.

Fitness benefits play a vital role in the retention of the Pi-ta susceptible alleles.

In plants, large numbers of R genes, which segregate as loci with alternative alleles conferring different levels of disease resistance to pathogens, have been maintained over a long period of evolution. The reason why hosts harbor susceptible alleles in view of their null contribution to resistance is unclear. In rice, a single copy gene, Pi-ta, segregates for 2 expressed clades of alleles, 1 resistant and the other susceptible. We simulated loss-of-function of the Pi-ta susceptible allele using the CRISPR/Cas9 system to detect subsequent fitness changes and obtained insights into fitness effects related to the retention of the Pi-ta susceptible allele. Our creation of an artificial knockout of the Pi-ta susceptible allele suffered fitness-related trait declines of up to 49% in terms of filled grain yield upon the loss of Pi-ta function. The Pi-ta susceptible alleles might serve as an off-switch to downstream immune signaling, thus contributing to the fine-tuning of plant defense responses. The results demonstrated that the susceptible Pi-ta alleles should have evolved pleiotropic functions, facilitating their retention in populations. As Pi-ta is a single copy gene with no paralogs in the genome, its function cannot be compensated by an alternative gene; whereas most other R genes form gene clusters by tandem duplications, and the function could be compensated by paralogs with high sequence similarity. This attempt to evaluate the fitness effects of the R gene in crops indicates that not all disease resistance genes incur fitness costs, which also provides a plausible explanation for how host genomes can tolerate the possible genetic load associated with a vast repertoire of R genes.

Genome sequencing and transcriptome analyses provide insights into the origin and domestication of water caltrop (Trapa spp., Lythraceae).

Humans have domesticated diverse species from across the plant kingdom; however, our current understanding of plant domestication is largely founded on major cereal crops. Here, we examine the evolutionary processes and genetic basis underlying the domestication of water caltrop (Trapa spp., Lythraceae), a traditional, yet presently underutilized non-cereal crop that sustained early Chinese agriculturalists. We generated a chromosome-level genome assembly of tetraploid T. natans, and then divided the allotetraploid genome into two subgenomes. Based on resequencing data from 57 accessions, representing cultivated diploid T. natans, wild T. natans (2x and 4x) and diploid T. incisa, we showed that water caltrop was likely first domesticated in the Yangtze River Valley as early as 6300 yr BP, and experienced a second improvement c. 800 years ago. We also provided strong support for an allotetraploid origin of T. natans within the past 230 000-310 000 years. By integrating selective sweep and transcriptome profiling analyses, we identified a number of genes potentially selected and/or differentially expressed during domestication, some of which likely contributed not only to larger fruit sizes but also to a more vigorous root system, facilitating nutrient uptake, environmental stress response and underwater photosynthesis. Our results shed light on the evolutionary and domestication history of water caltrop, one of the earliest domesticated crops in China. This study has implications for genomic-assisted breeding of this presently underutilized aquatic plant, and improves our general understanding of plant domestication.

Genome-wide Identification and Evolutionary Analysis of NBS-LRR Genes From Secale cereale.

Secale cereale is an important crop in the Triticeae tribe of the Poaceae family, and it has unique agronomic characteristics and genome properties. It possesses resistance to many diseases and serves as an important resource for the breeding of other Triticeae crops. We performed a genome-wide study on S. cereale to identify the largest group of plant disease resistance genes (R genes), the nucleotide-binding site-leucine-rich repeat receptor (NBS-LRR) genes. In its genome, 582 NBS-LRR genes were identified, including one from the RNL subclass and 581 from the CNL subclass. The NBS-LRR gene number in the S. cereale genome is greater than that in barley and the diploid wheat genomes. S. cereale chromosome 4 contains the largest number of NBS-LRR genes among the seven chromosomes, which is different from the pattern in barley and the genomes B and D of wheat but similar to that in the genome A of wheat. Further synteny analysis suggests that more NBS-LRR genes on chromosome 4 have been inherited from a common ancestor by S. cereale and the wheat genome A than the wheat genomes B and D. Phylogenetic analysis revealed that at least 740 NBS-LRR lineages are present in the common ancestor of S. cereale, Hordeum vulgare and Triticum urartu. However, most of them have only been inherited by one or two species, with only 65 of them preserved in all three species. The S. cereale genome inherited 382 of these ancestral NBS-LRR lineages, but 120 of them have been lost in both H. vulgare and T. urartu. This study provides the full NBS-LRR profile of the S. cereale genome, which is a resource for S. cereale breeding and indicates that S. cereale can be an important material for the molecular breeding of other Triticeae crops.

Comparative Chloroplast Genomes of Four Lycoris Species (Amaryllidaceae) Provides New Insight into Interspecific Relationship and Phylogeny.

The genus Lycoris (Amaryllidaceae) consists of about 20 species, which is endemic to East Asia. Although the Lycoris species is of great horticultural and medical importance, challenges in accurate species identification persist due to frequent natural hybridization and large-scale intraspecific variation. In this study, we sequenced chloroplast genomes of four Lycoris species and retrieved seven published chloroplast (cp) genome sequences in this genus for comparative genomic and phylogenetic analyses. The cp genomes of these four newly sequenced species were found to be 158,405-158,498 bp with the same GC content of 37.8%. The structure of the genomes exhibited the typical quadripartite structure with conserved gene order and content. A total of 113 genes (20 duplicated) were identified, including 79 protein-coding genes (PCGs), 30 tRNAs, and 4 rRNAs. Phylogenetic analysis showed that the 11 species were clustered into three main groups, and L. sprengeri locate at the base of Lycoriss. The L. radiata was suggested to be the female donor of the L. incarnata, L. shaanxiensis, and L. squamigera. The L. straminea and L. houdyshelii may be derived from L. anhuiensis, L. chinensis, or L. longituba. These results could not only offer a genome-scale platform for identification and utilization of Lycoris but also provide a phylogenomic framework for future studies in this genus.

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.