A chromosome-level genome of a Kordofan melon illuminates the origin of domesticated watermelons

Chromosome-level genome assembly and functional annotation of Citrullus colocynthis: unlocking genetic resources for drought-resilient crop development

MAIN CONCLUSION

The chromosome-level genome assembly of Citrullus colocynthis reveals its genetic potential for enhancing drought tolerance, paving the way for innovative crop improvement strategies. This study presents the first comprehensive genome assembly and annotation of Citrullus colocynthis, a drought-tolerant wild close relative of cultivated watermelon, highlighting its potential for enhancing agricultural resilience to climate change. The study achieved a chromosome-level assembly using advanced sequencing technologies, including PacBio HiFi and Hi-C, revealing a genome size of approximately 366 Mb with low heterozygosity and substantial repetitive content. Our analysis identified 23,327 gene models, that could encode stress response mechanisms for species' adaptation to arid environments. Comparative genomics with closely related species illuminated the evolutionary dynamics within the Cucurbitaceae family. In addition, resequencing of 27 accessions from the United Arab Emirates (UAE) identified genetic diversity, suggesting a foundation for future breeding programs. This genomic resource opens new avenues for the de novo domestication of C. colocynthis, offering a blueprint for developing crops with enhanced drought tolerance, disease resistance, and nutritional profiles, crucial for sustaining future food security in the face of escalating climate challenges.

Telomere-to-telomere Citrullus super-pangenome provides direction for watermelon breeding

To decipher the genetic diversity within the cucurbit genus Citrullus, we generated telomere-to-telomere (T2T) assemblies of 27 distinct genotypes, encompassing all seven Citrullus species. This T2T super-pangenome has expanded the previously published reference genome, T2T-G42, by adding 399.2 Mb and 11,225 genes. Comparative analysis has unveiled gene variants and structural variations (SVs), shedding light on watermelon evolution and domestication processes that enhanced attributes such as bitterness and sugar content while compromising disease resistance. Multidisease-resistant loci from Citrullus amarus and Citrullus mucosospermus were successfully introduced into cultivated Citrullus lanatus. The SVs identified in C. lanatus have not only been inherited from cordophanus but also from C. mucosospermus, suggesting additional ancestors beyond cordophanus in the lineage of cultivated watermelon. Our investigation substantially improves the comprehension of watermelon genome diversity, furnishing comprehensive reference genomes for all Citrullus species. This advancement aids in the exploration and genetic enhancement of watermelon using its wild relatives.

Quantitative Trait Loci Mapping and Comparative Transcriptome Analysis of Fruit Weight (FW) in Watermelon (Citrullus lanatus L.)

The watermelon (Citrullus lanatus L.) holds substantial economic value as a globally cultivated horticultural crop. However, the genetic architecture of watermelon fruit weight (FW) remains poorly understood. In this study, we used sh14-11 with small fruit and N14 with big fruit to construct 100 recombinant inbred lines (RILs). Based on whole-genome resequencing (WGR), 218,127 single nucleotide polymorphisms (SNPs) were detected to construct a high-quality genetic map. After quantitative trait loci (QTL) mapping, a candidate interval of 31-38 Mb on chromosome 2 was identified for FW. Simultaneously, the bulked segregant analysis (BSA) in the F2 population corroborated the identification of the same interval, encompassing the homologous gene linked to the known FW-related gene fas. Additionally, RNA-seq was carried out across 11 tissues from sh14-11 and N14, revealing expression profiles that identified 1695 new genes and corrected the annotation of 2941 genes. Subsequent differential expression analysis unveiled 8969 differentially expressed genes (DEGs), with 354 of these genes exhibiting significant differences across four key developmental stages. The integration of QTL mapping and differential expression analysis facilitated the identification of 14 FW-related genes, including annotated TGA and NAC transcription factors implicated in fruit development. This combined approach offers valuable insights into the genetic basis of FW, providing crucial resources for enhancing watermelon cultivation.

Molecular Mapping of Putative Genomic Regions Controlling Fruit and Seed Morphology of Watermelon

The genetic regulatory basis of qualitative and quantitative phenotypes of watermelon is being investigated in different types of molecular and genetic breeding studies around the world. In this study, biparental F2 mapping populations were developed over two experimental years, and the collected datasets of fruit and seed traits exhibited highly significant correlations. Whole-genome resequencing of comparative parental lines was performed and detected single nucleotide polymorphism (SNP) loci were converted into cleaved amplified polymorphic sequence (CAPS) markers. The screened polymorphic markers were genotyped in segregating populations and two genetic linkage maps were constructed, which covered a total of 2834.28 and 2721.45 centimorgan (cM) genetic lengths, respectively. A total of 22 quantitative trait loci (QTLs) for seven phenotypic traits were mapped; among them, five stable and major-effect QTLs (PC-8-1, SL-9-1, SWi-9-1, SSi-9-1, and SW-6-1) and four minor-effect QTLs (PC-2-1 and PC-2-2; PT-2-1 and PT-2-2; SL-6-1 and SSi-6-2; and SWi-6-1 and SWi-6-2) were observed with 3.77-38.98% PVE. The adjacent QTL markers showed a good fit marker-trait association, and a significant allele-specific contribution was also noticed for genetic inheritance of traits. Further, a total of four candidate genes (Cla97C09G179150, Cla97C09G179350, Cla97C09G180040, and Cla97C09G180100) were spotted in the stable colocalized QTLs of seed size linked traits (SL-9-1 and SWi-9-1) that showed non-synonymous type mutations. The gene expression trends indicated that the seed morphology had been formed in the early developmental stage and showed the genetic regulation of seed shape formation. Hence, we think that our identified QTLs and genes would provide powerful genetic insights for marker-assisted breeding aimed at improving the quality traits of watermelon.

Evidence that the domesticated fungus Leucoagaricus gongylophorus recycles its cytoplasmic contents as nutritional rewards to feed its leafcutter ant farmers

Leafcutter ants farm a fungal cultivar (Leucoagaricus gongylophorus) that converts inedible vegetation into food that sustains colonies with up to millions of workers. Analogous to edible fruits of crops domesticated by humans, L. gongylophorus has evolved specialized nutritional rewards-swollen hyphal cells called gongylidia that package metabolites and are consumed by ant farmers. Yet, little is known about how gongylidia form, and thus how fungal physiology and ant provisioning collectively govern farming performance. We explored the process of gongylidium formation using advanced microscopy to image the cultivar at scales of nanometers, and both in vitro experiments and in silico analyses to examine the mechanisms of gongylidia formation when isolated from ant farmers. We first used transmission electron, fluorescence, and confocal microscopy imaging to see inside hyphal cells. This imaging showed that the cultivar uses a process called autophagy to recycle its own cellular material (e.g. cytosol, mitochondria) and then shuttles the resulting metabolites into a vacuole whose continual expansion displaces other organelles and causes the gongylidium cell's bulging bulb-like appearance. We next used scanning electron microscopy and light microscopy to link this intracellular rearrangement to the external branching patterns of gongylidium cells as they clump together into edible bundles called staphyla. We next confirmed that autophagy plays a critical role in gongylidium formation both: (1) in vitro as gongylidium suppression occurred when isolated fungal cultures were grown on media with autophagy inhibitors, and (2) in silico as differential transcript expression (RNA-seq) analyses showed upregulation of multiple autophagy gene isoforms in gongylidia relative to undifferentiated hyphae. While autophagy is a ubiquitous and often highly derived process across the tree of life, our study reveals a new role for autophagy as a mechanism of functional integration between ant farmers and their fungal crop, and potentially as a signifier of higher-level homeostasis between uniquely life-time committed ectosymbionts.

Rtapas: An R Package to Assess Cophylogenetic Signal between Two Evolutionary Histories

Cophylogeny represents a framework to understand how ecological and evolutionary process influence lineage diversification. The recently developed algorithm Random Tanglegram Partitions provides a directly interpretable statistic to quantify the strength of cophylogenetic signal and incorporates phylogenetic uncertainty into its estimation, and maps onto a tanglegram the contribution to cophylogenetic signal of individual host-symbiont associations. We introduce Rtapas, an R package to perform Random Tanglegram Partitions. Rtapas applies a given global-fit method to random partial tanglegrams of a fixed size to identify the associations, terminals, and internal nodes that maximize phylogenetic congruence. This new package extends the original implementation with a new algorithm that examines the contribution to phylogenetic incongruence of each host-symbiont association and adds ParaFit, a method designed to test for topological congruence between two phylogenies, to the list of global-fit methods than can be applied. Rtapas facilitates and speeds up cophylogenetic analysis, as it can handle large phylogenies (100+ terminals) in affordable computational time as illustrated with two real-world examples. Rtapas can particularly cater for the need for causal inference in cophylogeny in two domains: (i) Analysis of complex and intricate host-symbiont evolutionary histories and (ii) assessment of topological (in)congruence between phylogenies produced with different DNA markers and specifically identify subsets of loci for phylogenetic analysis that are most likely to reflect gene-tree evolutionary histories. [Cophylogeny; cophylogenetic signal; gene tree incongruence; phylogenetic congruence; phylogenomics.].

A chromosome-level genome assembly of radish (Raphanus sativus L.) reveals insights into genome adaptation and differential bolting regulation

High-quality radish (Raphanus sativus) genome represents a valuable resource for agronomical trait improvements and understanding genome evolution among Brassicaceae species. However, existing radish genome assembly remains fragmentary, which greatly hampered functional genomics research and genome-assisted breeding. Here, using a NAU-LB radish inbred line, we generated a reference genome of 476.32 Mb with a scaffold N50 of 56.88 Mb by incorporating Illumina, PacBio and BioNano optical mapping techniques. Utilizing Hi-C data, 448.12 Mb (94.08%) of the assembled sequences were anchored to nine radish chromosomes with 40 306 protein-coding genes annotated. In total, 249.14 Mb (52.31%) comprised the repetitive sequences, among which long terminal repeats (LTRs, 30.31%) were the most abundant class. Beyond confirming the whole-genome triplication (WGT) event in R. sativus lineage, we found several tandem arrayed genes were involved in stress response process, which may account for the distinctive phenotype of high disease resistance in R. sativus. By comparing against the existing Xin-li-mei radish genome, a total of 2 108 573 SNPs, 7740 large insertions, 7757 deletions and 84 inversions were identified. Interestingly, a 647-bp insertion in the promoter of RsVRN1 gene can be directly bound by the DOF transcription repressor RsCDF3, resulting into its low promoter activity and late-bolting phenotype of NAU-LB cultivar. Importantly, introgression of this 647-bp insertion allele, RsVRN1^In-536 , into early-bolting genotype could contribute to delayed bolting time, indicating that it is a potential genetic resource for radish late-bolting breeding. Together, this genome resource provides valuable information to facilitate comparative genomic analysis and accelerate genome-guided breeding and improvement in radish.

Comparative chloroplast genome analysis of seven extant Citrullus species insight into genetic variation, phylogenetic relationships, and selective pressure

Citrullus ecirrhosus, Citrullus rehmii, and Citrullus naudinianus are three important related wild species of watermelon in the genus Citrullus, and their morphological differences are clear, however, their chloroplast genome differences remain unknown. This study is the first to assemble, analyze, and publish the complete chloroplast genomes of C. ecirrhosus, C. rehmii, and C. naudinianus. A comparative analysis was then conducted among the complete chloroplast genomes of seven extant Citrullus species, and the results demonstrated that the average genome sizes of Citrullus is 157,005 bp, a total of 130-133 annotated genes were identified, including 8 rRNA, 37 tRNA and 85-88 protein-encoding genes. Their gene content, order, and genome structure were similar. However, noncoding regions were more divergent than coding regions, and rps16-trnQ was a hypervariable fragment. Thirty-four polymorphic SSRs, 1,271 SNPs and 234 INDELs were identified. Phylogenetic trees revealed a clear phylogenetic relationship of Citrullus species, and the developed molecular markers (SNPs and rps16-trnQ) could be used for taxonomy in Citrullus. Three genes (atpB, clpP1, and rpoC2) were identified to undergo selection and would promote the environmental adaptation of Citrullus.

High-quality genome assembly and genetic mapping reveal a gene regulating flesh color in watermelon (Citrullus lanatus)

The unique color and type characteristics of watermelon fruits are regulated by many molecular mechanisms. However, it still needs to be combined with more abundant genetic data to fine-tune the positioning. We assembled genomes of two Korean inbred watermelon lines (cv. 242-1 and 159-1) with unique color and fruit-type characteristics and identified 23,921 and 24,451 protein-coding genes in the two genomes, respectively. To obtain more precise results for further study, we resequenced one individual of each parental line and an F₂ population composed of 87 individuals. This identified 1,539 single-nucleotide polymorphisms (SNPs) and 80 InDel markers that provided a high-density genetic linkage map with a total length of 3,036.9 cM. Quantitative trait locus mapping identified 15 QTLs for watermelon fruit quality-related traits, including β-carotene and lycopene content in fruit flesh, fruit shape index, skin thickness, flesh color, and rind color. By investigating the mapping intervals, we identified 33 candidate genes containing variants in the coding sequence. Among them, Cla97C01G008760 was annotated as a phytoene synthase with a single-nucleotide variant (A → G) in the first exon at 9,539,129 bp of chromosome 1 that resulted in the conversion of a lysine to glutamic acid, indicating that this gene might regulate flesh color changes at the protein level. These findings not only prove the importance of a phytoene synthase gene in pigmentation but also explain an important reason for the color change of watermelon flesh.

CuGenDBv2: an updated database for cucurbit genomics

The Cucurbitaceae (cucurbit) family consists of about 1,000 species in 95 genera, including many economically important and popular fruit and vegetable crops. During the past several years, reference genomes have been generated for >20 cucurbit species, and variome and transcriptome profiling data have been rapidly accumulated for cucurbits. To efficiently mine, analyze and disseminate these large-scale datasets, we have developed an updated version of Cucurbit Genomics Database. The updated database, CuGenDBv2 (http://cucurbitgenomics.org/v2), currently hosts 34 reference genomes from 27 cucurbit species/subspecies belonging to 10 different genera. Protein-coding genes from these genomes have been comprehensively annotated by comparing their protein sequences to various public protein and domain databases. A novel 'Genotype' module has been implemented to facilitate mining and analysis of the functionally annotated variome data including SNPs and small indels from large-scale genome sequencing projects. An updated 'Expression' module has been developed to provide a comprehensive gene expression atlas for cucurbits. Furthermore, synteny blocks between any two and within each of the 34 genomes, representing a total of 595 pair-wise genome comparisons, have been identified and can be explored and visualized in the database.

Diploid chromosome-level reference genome and population genomic analyses provide insights into Gypenoside biosynthesis and demographic evolution of Gynostemma pentaphyllum (Cucurbitaceae)

Gynostemma pentaphyllum (Thunb.) Makino is a perennial creeping herbaceous plant in the family Cucurbitaceae, which has great medicinal value and commercial potential, but urgent conservation efforts are needed due to the gradual decreases and fragmented distribution of its wild populations. Here, we report the high-quality diploid chromosome-level genome of G. pentaphyllum obtained using a combination of next-generation sequencing short reads, Nanopore long reads, and Hi-C sequencing technologies. The genome is anchored to 11 pseudo-chromosomes with a total size of 608.95 Mb and 26 588 predicted genes. Comparative genomic analyses indicate that G. pentaphyllum is estimated to have diverged from Momordica charantia 60.7 million years ago, with no recent whole-genome duplication event. Genomic population analyses based on genotyping-by-sequencing and ecological niche analyses indicated low genetic diversity but a strong population structure within the species, which could classify 32 G. pentaphyllum populations into three geographical groups shaped jointly by geographic and climate factors. Furthermore, comparative transcriptome analyses showed that the genes encoding enzyme involved in gypenoside biosynthesis had higher expression levels in the leaves and tendrils. Overall, the findings obtained in this study provide an effective molecular basis for further studies of demographic genetics, ecological adaption, and systematic evolution in Cucurbitaceae species, as well as contributing to molecular breeding, and the biosynthesis and biotransformation of gypenoside.

An allelic variant in the ACS7 gene promotes primary root growth in watermelon

Gene mining in a C. lanatus × C. amarus population revealed one gene, ACS7, linked to primary root elongation in watermelon. Watermelon is a xerophytic crop characterized by a long primary root and robust lateral roots. Therefore, watermelon serves as an excellent model for studying root elongation and development. However, the genetic mechanism underlying the primary root elongation in watermelon remains unknown. Herein, through bulk segregant analysis we identified a genetic locus, qPRL.Chr03, controlling primary root length (PRL) using two different watermelon species (Citrullus lanatus and Citrullus amarus) that differ in their root architecture. Fine mapping revealed that xaa-Pro dipeptidase and 1-aminocyclopropane-1-carboxylate synthase 7 (ACS7) are candidate regulators of the primary root growth. Allelic variation in the delimited region among 193 watermelon accessions indicated that the long-root alleles might only exist in C. amarus. Interestingly, the discrepancy in PRL among the C. amarus accessions was clearly associated with a nonsynonymous single nucleotide polymorphism variant within the ACS7 gene. The ACS7 expression and ethylene levels in the primary root tips suggested that ethylene is a negative regulator of root elongation in watermelon, as supported by the application of 1-aminocyclopropane-1-carboxylate (ACC, the ethylene precursor) or 2-aminoethoxyvinyl glycine (AVG, an ACS inhibitor). To the best of our knowledge, these findings provide the first description of the genetic basis of root elongation in watermelon. The detected markers of the ACS7 gene will facilitate marker-assisted selection for the PRL trait to improve water and nutrient use efficacy in watermelon and beyond.

Phylogenomic assessment prompts recognition of the Serianthes clade and confirms the monophyly of Serianthes and its relationship with Falcataria and Wallaceodendron in the wider ingoid clade (Leguminosae, Caesalpinioideae)

The Indo-Pacific legume genus Serianthes was recently placed in the Archidendron clade (sensu Koenen et al. 2020), a subclade of the mimosoid clade in subfamily Caesalpinioideae, which also includes Acacia, Archidendron, Archidendropsis, Falcataria, Pararchidendron, Paraserianthes and Wallaceodendron. Serianthes comprises ca. 18 species, five subspecies and two varieties that are characterised by bipinnately compound leaves with alternate sessile leaflets, branched axillary corymbiform panicles and woody indehiscent pods. Generic relationships, as well as species relationships within genera in the Archidendron clade, remain uncertain. While the sister relationship between Serianthes and the genus Falcataria is strongly supported by molecular data, the distinction between Serianthes and the monotypic genus Wallaceodendron has been questioned, based on their similar flower and fruit morphologies. We combined three gene-enriched hybrid capture DNA sequence datasets (generated from the 964 mimobaits v1 probe set, the expanded 997 mimobaits v2 probe set and the GoFlag angiosperm 408 probe set) and used their overlapping markers (77 loci of the target exonic and flanking regions) to test the monophyly of Serianthes and to investigate generic relationships within the Archidendron clade using 55 ingoid plus two outgroup taxa. We show that Serianthes is monophyletic, confirm the Serianthes + Falcataria sister relationship to Wallaceodendron and recognise this combined clade as the Serianthes clade within the Archidendron clade. We also evaluated the use of overlapping loci across datasets in combination with concordance analyses to test generic relationships and further investigate previously unresolved relationships across the wider ingoid clade. Concordance analysis revealed limited gene tree conflicts near the tips of the Archidendron clade, but increased discordance at the base of the clade, which could be attributed to rapid lineage divergence (radiation) and/or incomplete lineage sorting.

A telomere-to-telomere gap-free reference genome of watermelon and its mutation library provide important resources for gene discovery and breeding

Watermelon, Citrullus lanatus, is the world's third largest fruit crop. Reference genomes with gaps and a narrow genetic base hinder functional genomics and genetic improvement of watermelon. Here, we report the assembly of a telomere-to-telomere gap-free genome of the elite watermelon inbred line G42 by incorporating high-coverage and accurate long-read sequencing data with multiple assembly strategies. All 11 chromosomes have been assembled into single-contig pseudomolecules without gaps, representing the highest completeness and assembly quality to date. The G42 reference genome is 369 321 829 bp in length and contains 24 205 predicted protein-coding genes, with all 22 telomeres and 11 centromeres characterized. Furthermore, we established a pollen-EMS mutagenesis protocol and obtained over 200 000 M1 seeds from G42 . In a sampling pool, 48 monogenic phenotypic mutations, selected from 223 M1 and 78 M2 mutants with morphological changes, were confirmed. The average mutation density was 1 SNP/1.69 Mb and 1 indel/4.55 Mb per M1 plant and 1 SNP/1.08 Mb and 1 indel/6.25 Mb per M2 plant. Taking advantage of the gap-free G42 genome, 8039 mutations from 32 plants sampled from M1 and M2 families were identified with 100% accuracy, whereas only 25% of the randomly selected mutations identified using the 97103v2 reference genome could be confirmed. Using this library and the gap-free genome, two genes responsible for elongated fruit shape and male sterility (ClMS1) were identified, both caused by a single base change from G to A. The validated gap-free genome and its EMS mutation library provide invaluable resources for functional genomics and genetic improvement of watermelon.

Genome sequencing of up to 6,000-yr-old Citrullus seeds reveals use of a bitter-fleshed species prior to watermelon domestication

Iconographic evidence from Egypt suggests that watermelon pulp was consumed there as a dessert by 4,360 BP. Earlier archaeobotanical evidence comes from seeds from Neolithic settlements in Libya, but whether these were watermelons with sweet pulp or other forms is unknown. We generated genome sequences from 6,000- and 3,300-year-old seeds from Libya and Sudan, and from worldwide herbarium collections made between 1824 and 2019, and analyzed these data together with resequenced genomes from important germplasm collections for a total of 131 accessions. Phylogenomic and population-genomic analyses reveal that (1) much of the nuclear genome of both ancient seeds is traceable to West African seed-use “egusi-type” watermelon (Citrullus mucosospermus) rather than domesticated pulp-use watermelon (Citrullus lanatus ssp. vulgaris); (2) the 6,000-year-old watermelon likely had bitter pulp and greenish-white flesh as today found in C. mucosospermus, given alleles in the bitterness regulators ClBT and in the red color marker LYCB; and (3) both ancient genomes showed admixture from C. mucosospermus, C. lanatus ssp. cordophanus, C. lanatus ssp. vulgaris, and even South African Citrullus amarus, and evident introgression between the Libyan seed (UMB-6) and populations of C. lanatus. An unexpected new insight is that Citrullus appears to have initially been collected or cultivated for its seeds, not its flesh, consistent with seed damage patterns induced by human teeth in the oldest Libyan material.

Grafting Watermelon Onto Pumpkin Increases Chilling Tolerance by Up Regulating Arginine Decarboxylase to Increase Putrescine Biosynthesis

Low temperature is a major environmental factor that severely impairs plant growth and productivity. Watermelon (Citrullus lanatus) is a chilling-sensitive crop. Grafting of watermelon onto pumpkin rootstock is an effective technique to increase the chilling tolerance of watermelon when exposure to short-time chilling stress. However, the mechanism by which pumpkin rootstock increases chilling tolerance remains poorly understood. Under 10°C/5°C (day/night) chilling stress treatment, pumpkin-grafted watermelon seedlings showed higher chilling tolerance than self-grafted watermelon plants with significantly reduced lipid peroxidation and chilling injury (CI) index. Physiological analysis revealed that pumpkin rootstock grafting led to the notable accumulation of putrescine in watermelon seedlings under chilling conditions. Pre-treat foliar with 1 mM D-arginine (inhibitor of arginine decarboxylase, ADC) increased the electrolyte leakage (EL) of pumpkin-grafted watermelon leaves under chilling stress. This result can be ascribed to the decrease in transcript levels of ADC, ornithine decarboxylase, spermidine synthase, and polyamine oxidase genes involved in the synthesis and metabolism of polyamines. Transcriptome analysis showed that pumpkin rootstock improved chilling tolerance in watermelon seedlings by regulating differential gene expression under chilling stress. Pumpkin-grafted seedling reduced the number and expression level of differential genes in watermelon scion under chilling stress. It specifically increased the up-regulated expression of ADC (Cla97C11G210580), a key gene in the polyamine metabolism pathway, and ultimately promoted the accumulation of putrescine. In conclusion, pumpkin rootstock grafting increased the chilling tolerance of watermelon through transcription adjustments, up regulating the expression level of ADC, and promoting the synthesis of putrescine, which ultimately improved the chilling tolerance of pumpkin-grafted watermelon plants.

The Collaboration Between Art History and Genetics - An Unlikely Marriage of Disciplines

Our fruits, vegetables, and cereal crops stem from a wild ancestor and have undergone major changes through millennia of domestication and selection. There are various ways to reveal plant diversity over time, and one of these is through the combination of art history and genetics (also known as #ArtGenetics). Here, we discuss this approach from the art historian's point of view and flag the advantages and caveats of such an approach. We also advocate for the development of an integrated, global art database to facilitate such analyses.

Genetic Resources and Vulnerabilities of Major Cucurbit Crops

The Cucurbitaceae family provides numerous important crops including watermelons (Citrullus lanatus), melons (Cucumis melo), cucumbers (Cucumis sativus), and pumpkins and squashes (Cucurbita spp.). Centers of domestication in Africa, Asia, and the Americas were followed by distribution throughout the world and the evolution of secondary centers of diversity. Each of these crops is challenged by multiple fungal, oomycete, bacterial, and viral diseases and insects that vector disease and cause feeding damage. Cultivated varieties are constrained by market demands, the necessity for climatic adaptations, domestication bottlenecks, and in most cases, limited capacity for interspecific hybridization, creating narrow genetic bases for crop improvement. This analysis of crop vulnerabilities examines the four major cucurbit crops, their uses, challenges, and genetic resources. ex situ germplasm banks, the primary strategy to preserve genetic diversity, have been extensively utilized by cucurbit breeders, especially for resistances to biotic and abiotic stresses. Recent genomic efforts have documented genetic diversity, population structure, and genetic relationships among accessions within collections. Collection size and accessibility are impacted by historical collections, current ability to collect, and ability to store and maintain collections. The biology of cucurbits, with insect-pollinated, outcrossing plants, and large, spreading vines, pose additional challenges for regeneration and maintenance. Our ability to address ongoing and future cucurbit crop vulnerabilities will require a combination of investment, agricultural, and conservation policies, and technological advances to facilitate collection, preservation, and access to critical Cucurbitaceae diversity.

Hundreds of nuclear and plastid loci yield novel insights into orchid relationships

PREMISE

The inference of evolutionary relationships in the species-rich family Orchidaceae has hitherto relied heavily on plastid DNA sequences and limited taxon sampling. Previous studies have provided a robust plastid phylogenetic framework, which was used to classify orchids and investigate the drivers of orchid diversification. However, the extent to which phylogenetic inference based on the plastid genome is congruent with the nuclear genome has been only poorly assessed.

METHODS

We inferred higher-level phylogenetic relationships of orchids based on likelihood and ASTRAL analyses of 294 low-copy nuclear genes sequenced using the Angiosperms353 universal probe set for 75 species (representing 69 genera, 16 tribes, 24 subtribes) and a concatenated analysis of 78 plastid genes for 264 species (117 genera, 18 tribes, 28 subtribes). We compared phylogenetic informativeness and support for the nuclear and plastid phylogenetic hypotheses.

RESULTS

Phylogenetic inference using nuclear data sets provides well-supported orchid relationships that are highly congruent between analyses. Comparisons of nuclear gene trees and a plastid supermatrix tree showed that the trees are mostly congruent, but revealed instances of strongly supported phylogenetic incongruence in both shallow and deep time. The phylogenetic informativeness of individual Angiosperms353 genes is in general better than that of most plastid genes.

CONCLUSIONS

Our study provides the first robust nuclear phylogenomic framework for Orchidaceae and an assessment of intragenomic nuclear discordance, plastid-nuclear tree incongruence, and phylogenetic informativeness across the family. Our results also demonstrate what has long been known but rarely thoroughly documented: nuclear and plastid phylogenetic trees can contain strongly supported discordances, and this incongruence must be reconciled prior to interpretation in evolutionary studies, such as taxonomy, biogeography, and character evolution.

Molecular clocks and archaeogenomics of a Late Period Egyptian date palm leaf reveal introgression from wild relatives and add timestamps on the domestication

The date palm, Phoenix dactylifera, has been a cornerstone of Middle Eastern and North African agriculture for millennia. It was first domesticated in the Persian Gulf, and its evolution appears to have been influenced by gene flow from two wild relatives, P. theophrasti, currently restricted to Crete and Turkey, and P. sylvestris, widespread from Bangladesh to the West Himalayas. Genomes of ancient date palm seeds show that gene flow from P. theophrasti to P. dactylifera may have occurred by ∼2,200 years ago, but traces of P. sylvestris could not be detected. We here integrate archeogenomics of a ∼2,100-year-old P. dactylifera leaf from Saqqara (Egypt), molecular-clock dating, and coalescence approaches with population genomic tests, to probe the hybridization between the date palm and its two closest relatives and provide minimum and maximum timestamps for its reticulated evolution. The Saqqara date palm shares a close genetic affinity with North African date palm populations, and we find clear genomic admixture from both P. theophrasti, and P. sylvestris, indicating that both had contributed to the date palm genome by 2,100 years ago. Molecular-clocks placed the divergence of P. theophrasti from P. dactylifera/P. sylvestris and that of P. dactylifera from P. sylvestris in the Upper Miocene, but strongly supported, conflicting topologies point to older gene flow between P. theophrasti and P. dactylifera, and P. sylvestris and P. dactylifera. Our work highlights the ancient hybrid origin of the date palms, and prompts the investigation of the functional significance of genetic material introgressed from both close relatives, which in turn could prove useful for modern date palm breeding.

Plastome Evolution in the Hyperdiverse Genus Euphorbia (Euphorbiaceae) Using Phylogenomic and Comparative Analyses: Large-Scale Expansion and Contraction of the Inverted Repeat Region

With c. 2,000 species, Euphorbia is one of the largest angiosperm genera, yet a lack of chloroplast genome (plastome) resources impedes a better understanding of its evolution. In this study, we assembled and annotated 28 plastomes from Euphorbiaceae, of which 15 were newly sequenced. Phylogenomic and comparative analyses of 22 plastome sequences from all four recognized subgenera within Euphorbia revealed that plastome length in Euphorbia is labile, presenting a range of variation c. 42 kb. Large-scale expansions of the inverted repeat (IR) region were identified, and at the extreme opposite, the near-complete loss of the IR region (with only 355 bp left) was detected for the first time in Euphorbiaceae. Other structural variations, including gene inversion and duplication, and gene loss/pseudogenization, were also observed. We screened the most promising molecular markers from both intergenic and coding regions for phylogeny-based utilities, and estimated maximum likelihood and Bayesian phylogenies from four datasets including whole plastome sequences. The monophyly of Euphorbia is supported, and its four subgenera are recovered in a successive sister relationship. Our study constitutes the first comprehensive investigation on the plastome structural variation in Euphorbia and it provides resources for phylogenetic research in the genus, facilitating further studies on its taxonomy, evolution, and conservation.