The genetic basis of sex determination in grapes

Multi-omics analysis reveals promiscuous O-glycosyltransferases involved in the diversity of flavonoid glycosides in Periploca forrestii (Apocynaceae)

Flavonoid glycosides are widespread in plants, and are of great interest owing to their diverse biological activities and effectiveness in preventing chronic diseases. Periploca forrestii, a renowned medicinal plant of the Apocynaceae family, contains diverse flavonoid glycosides and is clinically used to treat rheumatoid arthritis and traumatic injuries. However, the mechanisms underlying the biosynthesis of these flavonoid glycosides have not yet been elucidated. In this study, we used widely targeted metabolomics and full-length transcriptome sequencing to identify flavonoid diversity and biosynthetic genes in P. forrestii. A total of 120 flavonoid glycosides, including 21 C-, 96 O-, and 3 C/O-glycosides, were identified and annotated. Based on 24,123 full-length coding sequences, 99 uridine diphosphate sugar-utilizing glycosyltransferases (UGTs) were identified and classified into 14 groups. Biochemical assays revealed that four UGTs exhibited O-glycosyltransferase activity toward apigenin and luteolin. Among them, PfUGT74B4 and PfUGT92A8 were highly promiscuous and exhibited multisite O-glycosylation or consecutive glycosylation activities toward various flavonoid aglycones. These four glycosyltransferases may significantly contribute to the diversity of flavonoid glycosides in P. forrestii. Our findings provide a valuable genetic resource for further studies on P. forrestii and insights into the metabolic engineering of bioactive flavonoid glycosides.

Population comparative genomics discovers gene gain and loss during grapevine domestication

Plant domestication are evolutionary experiments conducted by early farmers since thousands years ago, during which the crop wild progenitors are artificially selected for desired agronomic traits along with dramatic genomic variation in the course of moderate to severe bottlenecks. However, previous investigations are mainly focused on small-effect variants, while changes in gene contents are rarely investigated due to the lack of population-level assemblies for both the crop and its wild relatives. Here, we applied comparative genomic analyses to discover gene gain and loss during grapevine domestication using long-read assemblies of representative population samples for both domesticated grapevines (V. vinifera ssp. vinifera) and their wild progenitors (V. vinifera ssp. sylvestris). Only ∼7% of gene families were shared by 16 Vitis genomes while ∼8% of gene families were specific to each accession, suggesting dramatic variations of gene contents in grapevine genomes. Compared to wild progenitors, the domesticated accessions exhibited an increased presence of genes associated with asexual reproduction, while the wild progenitors showcased a higher abundance of genes related to pollination, revealing the transition from sexual reproduction to clonal propagation during domestication processes. Moreover, the domesticated accessions harbored fewer disease-resistance genes than wild progenitors. The SVs occurred frequently in aroma and disease-resistance related genes between domesticated grapevines and wild progenitors, indicating the rapid diversification of these genes during domestication. Our study provides insights and resources for biological studies and breeding programs in grapevine.

The wild side of grape genomics

With broad genetic diversity and as a source of key agronomic traits, wild grape species (Vitis spp.) are crucial to enhance viticulture's climatic resilience and sustainability. This review discusses how recent breakthroughs in the genome assembly and analysis of wild grape species have led to discoveries on grape evolution, from wild species' adaptation to environmental stress to grape domestication. We detail how diploid chromosome-scale genomes from wild Vitis spp. have enabled the identification of candidate disease-resistance and flower sex determination genes and the creation of the first Vitis graph-based pangenome. Finally, we explore how wild grape genomics can impact grape research and viticulture, including aspects such as data sharing, the development of functional genomics tools, and the acceleration of genetic improvement.

Diploid genome assembly of the Malbec grapevine cultivar enables haplotype-aware analysis of transcriptomic differences underlying clonal phenotypic variation

To preserve their varietal attributes, established grapevine cultivars (Vitis vinifera L. ssp. vinifera) must be clonally propagated, due to their highly heterozygous genomes. Malbec is a France-originated cultivar appreciated for producing high-quality wines and is the offspring of cultivars Prunelard and Magdeleine Noire des Charentes. Here, we have built a diploid genome assembly of Malbec, after trio binning of PacBio long reads into the two haploid complements inherited from either parent. After haplotype-aware deduplication and corrections, complete assemblies for the two haplophases were obtained with a very low haplotype switch-error rate (<0.025). The haplophase alignment identified > 25% of polymorphic regions. Gene annotation including RNA-seq transcriptome assembly and ab initio prediction evidence resulted in similar gene model numbers for both haplophases. The annotated diploid assembly was exploited in the transcriptomic comparison of four clonal accessions of Malbec that exhibited variation in berry composition traits. Analysis of the ripening pericarp transcriptome using either haplophases as a reference yielded similar results, although some differences were observed. Particularly, among the differentially expressed genes identified only with the Magdeleine-inherited haplotype as reference, we observed an over-representation of hypothetically hemizygous genes. The higher berry anthocyanin content of clonal accession 595 was associated with increased abscisic acid responses, possibly leading to the observed overexpression of phenylpropanoid metabolism genes and deregulation of genes associated with abiotic stress response. Overall, the results highlight the importance of producing diploid assemblies to fully represent the genomic diversity of highly heterozygous woody crop cultivars and unveil the molecular bases of clonal phenotypic variation.

Why should we study plant sex chromosomes?

Understanding plant sex chromosomes involves studying interactions between developmental and physiological genetics, genome evolution, and evolutionary ecology. We focus on areas of overlap between these. Ideas about how species with separate sexes (dioecious species, in plant terminology) can evolve are even more relevant to plants than to most animal taxa because dioecy has evolved many times from ancestral functionally hermaphroditic populations, often recently. One aim of studying plant sex chromosomes is to discover how separate males and females evolved from ancestors with no such genetic sex-determining polymorphism, and the diversity in the genetic control of maleness vs femaleness. Different systems share some interesting features, and their differences help to understand why completely sex-linked regions may evolve. In some dioecious plants, the sex-determining genome regions are physically small. In others, regions without crossing over have evolved sometimes extensive regions with properties very similar to those of the familiar animal sex chromosomes. The differences also affect the evolutionary changes possible when the environment (or pollination environment, for angiosperms) changes, as dioecy is an ecologically risky strategy for sessile organisms. Dioecious plants have repeatedly reverted to cosexuality, and hermaphroditic strains of fruit crops such as papaya and grapes are desired by plant breeders. Sex-linked regions are predicted to become enriched in genes with sex differences in expression, especially when higher expression benefits one sex function but harms the other. Such trade-offs may be important for understanding other plant developmental and physiological processes and have direct applications in plant breeding.

From TgO/GABA-AT, GABA, and T-263 Mutant to Conception of Toxoplasma

Toxoplasma gondii causes morbidity, mortality, and disseminates widely via cat sexual stages. Here, we find T. gondii ornithine aminotransferase (OAT) is conserved across phyla. We solve TgO/GABA-AT structures with bound inactivators at 1.55 Å and identify an inactivator selective for TgO/GABA-AT over human OAT and GABA-AT. However, abrogating TgO/GABA-AT genetically does not diminish replication, virulence, cyst-formation, or eliminate cat's oocyst shedding. Increased sporozoite/merozoite TgO/GABA-AT expression led to our study of a mutagenized clone with oocyst formation blocked, arresting after forming male and female gametes, with "Rosetta stone"-like mutations in genes expressed in merozoites. Mutations are similar to those in organisms from plants to mammals, causing defects in conception and zygote formation, affecting merozoite capacitation, pH/ionicity/sodium-GABA concentrations, drawing attention to cyclic AMP/PKA, and genes enhancing energy or substrate formation in TgO/GABA-AT-related-pathways. These candidates potentially influence merozoite's capacity to make gametes that fuse to become zygotes, thereby contaminating environments and causing disease.

Sex-limited experimental evolution drives transcriptomic divergence in a hermaphrodite

The evolution of gonochorism from hermaphroditism is linked with the formation of sex chromosomes, as well as the evolution of sex-biased and sex-specific gene expression to allow both sexes to reach their fitness optimum. There is evidence that sexual selection drives the evolution of male-biased gene expression in particular. However, previous research in this area in animals comes from either theoretical models or comparative studies of already old sex chromosomes. We therefore investigated changes in gene expression under 3 different selection regimes for the simultaneous hermaphrodite Macrostomum lignano subjected to sex-limited experimental evolution (i.e. selection for fitness via eggs, sperm, or a control regime allowing both). After 21 and 22 generations of selection for male-specific or female-specific fitness, we characterized changes in whole-organism gene expression. We found that female-selected lines had changed the most in their gene expression. Although annotation for this species is limited, gene ontology term and Kyoto Encyclopedia of Genes and Genomes pathway analyses suggest that metabolic changes (e.g. biosynthesis of amino acids and carbon metabolism) are an important adaptive component. As predicted, we found that the expression of genes previously identified as testis-biased candidates tended to be downregulated in the female-selected lines. We did not find any significant expression differences for previously identified candidates of other sex-specific organs, but this may simply reflect that few transcripts have been characterized in this way. In conclusion, our experiment suggests that changes in testis-biased gene expression are important in the early evolution of sex chromosomes and gonochorism.

A telomere-to-telomere reference genome of ficus (Ficus hispida) provides new insights into sex determination

A high-quality reference genome is indispensable for resolving biologically essential traits. Ficus hispida is a dioecious plant. A complete Ficus reference genome will be crucial for understanding their sex evolution and important biological characteristics, such as aerial roots, mutualistic symbiosis with ficus-wasps, and fruiting from old stems. Here, we generated a telomere-to-telomere (T2T) genome for F. hispida using PacBio HiFi and Oxford Nanopore Ultra-long sequencing technologies. The genome contiguity and completeness has shown improvement compared with the previously released genome, with the annotation of six centromeres and 28 telomeres. We have refined our previously reported 2-Mb male-specific region into a 7.2-Mb genomic region containing 51 newly predicted genes and candidate sex-determination genes AG2 and AG3. Many of these genes showed extremely low expression, likely attributed to hypermethylation in the gene body and promoter regions. Gene regulatory networks (GRNs) revealed that AG2 and AG3 are related to the regulation of stamen development in male flowers, while the AG1 gene is responsible for regulating female flowers' defense responses and secondary metabolite processes. Comparative analysis of GRNs showed that the NAC, WRKY, and MYB transcription factor families dominate the female GRN, whereas the MADS and MYB transcription factor families are prevalent in the male GRN.

A super-pangenome of the North American wild grape species

BACKGROUND

Capturing the genetic diversity of wild relatives is crucial for improving crops because wild species are valuable sources of agronomic traits that are essential to enhance the sustainability and adaptability of domesticated cultivars. Genetic diversity across a genus can be captured in super-pangenomes, which provide a framework for interpreting genomic variations.

RESULTS

Here we report the sequencing, assembly, and annotation of nine wild North American grape genomes, which are phased and scaffolded at chromosome scale. We generate a reference-unbiased super-pangenome using pairwise whole-genome alignment methods, revealing the extent of the genomic diversity among wild grape species from sequence to gene level. The pangenome graph captures genomic variation between haplotypes within a species and across the different species, and it accurately assesses the similarity of hybrids to their parents. The species selected to build the pangenome are a great representation of the genus, as illustrated by capturing known allelic variants in the sex-determining region and for Pierce's disease resistance loci. Using pangenome-wide association analysis, we demonstrate the utility of the super-pangenome by effectively mapping short reads from genus-wide samples and identifying loci associated with salt tolerance in natural populations of grapes.

CONCLUSIONS

This study highlights how a reference-unbiased super-pangenome can reveal the genetic basis of adaptive traits from wild relatives and accelerate crop breeding research.

Boosting grapevine breeding for climate-smart viticulture: from genetic resources to predictive genomics

The multifaceted nature of climate change is increasing the urgency to select resilient grapevine varieties, or generate new, fitter cultivars, to withstand a multitude of new challenging conditions. The attainment of this goal is hindered by the limiting pace of traditional breeding approaches, which require decades to result in new selections. On the other hand, marker-assisted breeding has proved useful when it comes to traits governed by one or few genes with great effects on the phenotype, but its efficacy is still restricted for complex traits controlled by many loci. On these premises, innovative strategies are emerging which could help guide selection, taking advantage of the genetic diversity within the Vitis genus in its entirety. Multiple germplasm collections are also available as a source of genetic material for the introgression of alleles of interest via adapted and pioneering transformation protocols, which present themselves as promising tools for future applications on a notably recalcitrant species such as grapevine. Genome editing intersects both these strategies, not only by being an alternative to obtain focused changes in a relatively rapid way, but also by supporting a fine-tuning of new genotypes developed with other methods. A review on the state of the art concerning the available genetic resources and the possibilities of use of innovative techniques in aid of selection is presented here to support the production of climate-smart grapevine genotypes.

Genealogical Analyses of 3 Cultivated and 1 Wild Specimen of Vitis vinifera from Greece

Grapevine (Vitis vinifera) has been an important crop with considerable cultural and economic significance for over 2,500 years, and Greece has been an important entry point into Europe for lineages that were domesticated in Western Asia and the Caucasus. However, whole-genome-based investigation of the demographic history of Greek cultivars relative to other European lineages has only started recently. To understand how Greek cultivars relate to Eurasian domesticated and wild populations, we sequenced 3 iconic domesticated strains ('Xinomavro,' 'Agiorgitiko,' 'Mavrotragano') along with 1 wild accession (the vinetree of Pausanias-a historically important wild specimen) and analyzed their genomic diversity together with a large sample of publicly available domesticated and wild strains. We also reconstructed genealogies by leveraging the powerful tsinfer methodology which has not previously been used in this system. We show that cultivated strains from Greece differ genetically from other strains in Europe. Interestingly, all the 3 cultivated Greek strains clustered with cultivated and wild accessions from Transcaucasia, South Asia, and the Levant and are amongst the very few cultivated European strains belonging to this cluster. Furthermore, our results indicate that 'Xinomavro' shares close genealogical proximity with European elite cultivars such as 'Chardonnay,' 'Riesling,' and 'Gamay' but not 'Pinot.' Therefore, the proximity of 'Xinomavro' to Gouais/Heunisch Weiss is confirmed and the utility of ancestral recombination graph reconstruction approaches to study genealogical relationships in crops is highlighted.

Genetic insights into the dissolution of dioecy in diploid persimmon Diospyros oleifera Cheng

BACKGROUND

Dioecy, a sexual system of single-sexual (gynoecious/androecious) individuals, is rare in flowering plants. This rarity may be a result of the frequent transition from dioecy into systems with co-sexual individuals.

RESULTS

In this study, co-sexual expression (monoecy and hermaphroditic development), previously thought to be polyploid-specific in Diospyros species, was identified in the diploid D. oleifeara historically. We characterized potential genetic mechanisms that underlie the dissolution of dioecy to monoecy and andro(gyno)monoecy, based on multiscale genome-wide investigations of 150 accessions of Diospyros oleifera. We found all co-sexual plants, including monoecious and andro(gyno)monoecious individuals, possessed the male determinant gene OGI, implying the presence of genetic factors controlling gynoecia development in genetically male D. oleifera. Importantly, discrepancies in the OGI/MeGI module were found in diploid monoecious D. oleifera compared with polyploid monoecious D. kaki, including no Kali insertion on the promoter of OGI, no different abundance of smRNAs targeting MeGI (a counterpart of OGI), and no different expression of MeGI between female and male floral buds. On the contrary, in both single- and co-sexual plants, female function was expressed in the presence of a genome-wide decrease in methylation levels, along with sexually distinct regulatory networks of smRNAs and their targets. Furthermore, a genome-wide association study (GWAS) identified a genomic region and a DUF247 gene cluster strongly associated with the monoecious phenotype and several regions that may contribute to andromonoecy.

CONCLUSIONS

Collectively, our findings demonstrate stable breakdown of the dioecious system in D. oleifera, presumably also a result of genomic features of the Y-linked region.

A multitiered haplotype strategy to enhance phased assembly and fine mapping of a disease resistance locus

Fine mapping of quantitative trait loci (QTL) to dissect the genetic basis of traits of interest is essential to modern breeding practice. Here, we employed a multitiered haplotypic marker system to increase fine mapping accuracy by constructing a chromosome-level, haplotype-resolved parental genome, accurate detection of recombination sites, and allele-specific characterization of the transcriptome. In the first tier of this system, we applied the preexisting panel of 2,000 rhAmpSeq core genome markers that is transferable across the entire Vitis genus and provides a genomic resolution of 200 kb to 1 Mb. The second tier consisted of high-density haplotypic markers generated from Illumina skim sequencing data for samples enriched for relevant recombinations, increasing the potential resolution to hundreds of base pairs. We used this approach to dissect a novel Resistance to Plasmopara viticola-33 (RPV33) locus conferring resistance to grapevine downy mildew, narrowing the candidate region to only 0.46 Mb. In the third tier, we used allele-specific RNA-seq analysis to identify a cluster of 3 putative disease resistance RPP13-like protein 2 genes located tandemly in a nonsyntenic insertion as candidates for the disease resistance trait. In addition, combining the rhAmpSeq core genome haplotype markers and skim sequencing-derived high-density haplotype markers enabled chromosomal-level scaffolding and phasing of the grape Vitis × doaniana 'PI 588149' assembly, initially built solely from Pacific Biosciences (PacBio) high-fidelity (HiFi) reads, leading to the correction of 16 large-scale phasing errors. Our mapping strategy integrates high-density, phased genetic information with individual reference genomes to pinpoint the genetic basis of QTLs and will likely be widely adopted in highly heterozygous species.

Sex-Related Gene Network Revealed by Transcriptome Differentiation of Bisexual and Unisexual Flowers of Orchid Cymbidium tortisepalum

Despite extensive research on orchid reproductive strategies, the genetic studies of sex differentiation in the orchid family are still lacking. In this study, we compared three sexual phenotypes of Cymbidium tortisepalum bisexual flowers as well as female and male unisexual mutants. Through comparative transcriptomes, we analyzed the sex-biased differentially expressed genes (DEGs) and gene co-expression networks of sex organs (gynostemium and ovary) among them, identified the candidate genes of sex differentiation, and validated their expression by qRT-PCR. The C. tortisepalum unisexual mutants with degenerated phenotypes were compared to the bisexual plants with respect to both the flower organs and plant morphologies. Totally, 12,145, 10,789, and 14,447 genes were uniquely expressed in the female, male, and hermaphrodite sex organs, respectively. A total of 4291 sex-biased DEGs were detected among them, with 871, 2867, and 1937 DEGs in the comparisons of bisexual vs. female, bisexual vs. male, and male vs. female flowers, respectively. Two co-expressed network modules, with 81 and 419 genes were tightly correlated with female sexual traits, while two others with 265 and 135 genes were highly correlated with male sexual traits. Two female-biased hub genes (CtSDR3b and CtSDR3b-like) nested in the female modules, the homologs of maize sex determinant tasselseed2, may control the feminization of C. tortisepalum. At the same time, two male-biased hub genes (CtYAB2 and CtYAB5) nested in the male modules, the homologs of grape sex determinant VviYABBY3, may control the androphany of C. tortisepalum. This study discovered the molecular regulation networks and proposed a model for orchid sex differentiation, therefore providing for the first time the genetic basis for the sex separation in the orchid family.

The haplotype-resolved T2T genome of teinturier cultivar Yan73 reveals the genetic basis of anthocyanin biosynthesis in grapes

Teinturier grapes are characterized by the typical accumulation of anthocyanins in grape skin, flesh, and vegetative tissues, endowing them with high utility value in red wine blending and nutrient-enriched foods developing. However, due to the lack of genome information, the mechanism involved in regulating teinturier grape coloring has not yet been elucidated and their genetic utilization research is still insufficient. Here, the cultivar 'Yan73' was used for assembling the telomere-to-telomere (T2T) genome of teinturier grapes by combining the High Fidelity (HiFi)， Hi-C and ultralong Oxford Nanopore Technologies (ONT) reads. Two haplotype genomes were assembled, at the sizes of 501.68 Mb and 493.38 Mb, respectively. In the haplotype 1 genome, the transposable elements (TEs) contained 32.77% of long terminal repeats (LTRs), while in the haplotype 2 genome, 31.53% of LTRs were detected in TEs. Furthermore, obvious inversions were identified in chromosome 18 between the two haplotypes. Transcriptome profiling suggested that the gene expression patterns in 'Cabernet Sauvignon' and 'Yan73' were diverse depending on tissues, developmental stages, and varieties. The transcription program of genes in the anthocyanins biosynthesis pathway between the two cultivars exhibited high similarity in different tissues and developmental stages, whereas the expression levels of numerous genes showed significant differences. Compared with other genes, the expression levels of VvMYBA1 and VvUFGT4 in all samples, VvCHS2 except in young shoots and VvPAL9 except in the E-L23 stage of 'Yan73' were higher than those of 'Cabernet Sauvignon'. Further sequence alignments revealed potential variant gene loci and structure variations of anthocyanins biosynthesis related genes and a 816 bp sequence insertion was found in the promoter of VvMYBA1 of 'Yan73' haplotype 2 genome. The 'Yan73' T2T genome assembly and comparative analysis provided valuable foundations for further revealing the coloring mechanism of teinturier grapes and the genetic improvement of grape coloring traits.

Phylogenomic and syntenic data demonstrate complex evolutionary processes in early radiation of the rosids

Some of the most vexing problems of deep level relationship that remain in angiosperms involve the superrosids. The superrosid clade contains a quarter of all angiosperm species, with 18 orders in three subclades (Vitales, Saxifragales and core rosids) exhibiting remarkable morphological and ecological diversity. To help resolve deep-level relationships, we constructed a high-quality chromosome-level genome assembly for Tiarella polyphylla (Saxifragaceae) thus providing broader genomic representation of Saxifragales. Whole genome microsynteny analysis of superrosids showed that Saxifragales shared more synteny clusters with core rosids than Vitales, further supporting Saxifragales as more closely related with core rosids. To resolve the ordinal phylogeny of superrosids, we screened 122 single copy nuclear genes from genomes of 36 species, representing all 18 superrosid orders. Vitales were recovered as sister to all other superrosids (Saxifragales + core rosids). Our data suggest dramatic differences in relationships compared to earlier studies within core rosids. Fabids should be restricted to the nitrogen-fixing clade, while Picramniales, the Celastrales-Malpighiales (CM) clade, Huerteales, Oxalidales, Sapindales, Malvales and Brassicales formed an "expanded" malvid clade. The Celastrales-Oxalidales-Malpighiales (COM) clade (sensu APG IV) was not monophyletic. Crossosomatales, Geraniales, Myrtales and Zygophyllales did not belong to either of our well-supported malvids or fabids. There is strong discordance between nuclear and plastid phylogenetic hypotheses for superrosid relationships; we show that this is best explained by a combination of incomplete lineage sorting and ancient reticulation.

Vitis labrusca genome assembly reveals diversification between wild and cultivated grapevine genomes

Wild grapevines are important genetic resources in breeding programs to confer adaptive fitness traits and unique fruit characteristics, but the genetics underlying these traits, and their evolutionary origins, are largely unknown. To determine the factors that contributed to grapevine genome diversification, we performed comprehensive intragenomic and intergenomic analyses with three cultivated European (including the PN40024 reference genome) and two wild North American grapevine genomes, including our newly released Vitis labrusca genome. We found the heterozygosity of the cultivated grapevine genomes was twice as high as the wild grapevine genomes studied. Approximately 30% of V. labrusca and 48% of V. vinifera Chardonnay genes were heterozygous or hemizygous and a considerable number of collinear genes between Chardonnay and V. labrusca had different gene zygosity. Our study revealed evidence that supports gene gain-loss events in parental genomes resulted in the inheritance of hemizygous genes in the Chardonnay genome. Thousands of segmental duplications supplied source material for genome-specific genes, further driving diversification of the genomes studied. We found an enrichment of recently duplicated, adaptive genes in similar functional pathways, but differential retention of environment-specific adaptive genes within each genome. For example, large expansions of NLR genes were discovered in the two wild grapevine genomes studied. Our findings support variation in transposable elements contributed to unique traits in grapevines. Our work revealed gene zygosity, segmental duplications, gene gain-and-loss variations, and transposable element polymorphisms can be key driving forces for grapevine genome diversification.

Identification of grapevine clones via high-throughput amplicon sequencing: a proof-of-concept study

Wine cultivars are available to growers in multiple clonal selections with agronomic and enological differences. Phenotypic differences between clones originated from somatic mutations that accrued over thousands of asexual propagation cycles. Genetic diversity between grape cultivars remains unexplored, and tools to discriminate unequivocally clones have been lacking. This study aimed to uncover genetic variations among a group of clonal selections of 4 important Vitis vinifera cultivars: Cabernet sauvignon, Sauvignon blanc, Chardonnay, and Merlot, and use this information to develop genetic markers to discriminate the clones of these cultivars. We sequenced with short-read sequencing technology the genomes of 18 clones, including biological replicates for a total of 46 genomes. Sequences were aligned to their respective cultivar's reference genome for variant calling. We used reference genomes of Cabernet sauvignon, Chardonnay, and Merlot and developed a de novo genome assembly of Sauvignon blanc using long-read sequencing. On average, 4 million variants were detected for each clone, with 74.2% being single nucleotide variants and 25.8% being small insertions or deletions (InDel). The frequency of these variants was consistent across all clones. From these variants, we validated 46 clonal markers using high-throughput amplicon sequencing for 77.7% of the evaluated clones, most of them small InDel. These results represent an advance in grapevine genotyping strategies and will benefit the viticulture industry for the characterization and identification of the plant material.

Chimeras in Merlot grapevine revealed by phased assembly

Chimerism is the phenomenon when several genotypes coexist in a single individual. Used to understand plant ontogenesis they also have been valorised through new cultivar breeding. Viticulture has been taking economic advantage out of chimeras when the variant induced an important modification of wine type such as berry skin colour. Crucial agronomic characters may also be impacted by chimeras that aren't identified yet. Periclinal chimera where the variant has entirely colonised a cell layer is the most stable and can be propagated through cuttings. In grapevine, leaves are derived from both meristem layers, L1 and L2. However, lateral roots are formed from the L2 cell layer only. Thus, comparing DNA sequences of roots and leaves allows chimera detection. In this study we used new generation Hifi long reads sequencing, recent bioinformatics tools and trio-binning with parental sequences to detect periclinal chimeras on 'Merlot' grapevine cultivar. Sequencing of cv. 'Magdeleine Noire des Charentes' and 'Cabernet Franc', the parents of cv. 'Merlot', allowed haplotype resolved assembly. Pseudomolecules were built with a total of 33 to 47 contigs and in few occasions a unique contig for one chromosome. This high resolution allowed haplotype comparison. Annotation was transferred from PN40024 VCost.v3 to all pseudomolecules. After strong selection of variants, 51 and 53 'Merlot' specific periclinal chimeras were found on the Merlot-haplotype-CF and Merlot-haplotype-MG respectively, 9 and 7 been located in a coding region. A subset of positions was analysed using Molecular Inversion Probes (MIPseq) and 69% were unambiguously validated, 25% are doubtful because of technological noise or weak depth and 6% invalidated. These results open new perspectives on chimera detection as an important resource to improve cultivars through clonal selection or breeding.

Adaptive and maladaptive introgression in grapevine domestication

Domesticated grapevines spread to Europe around 3,000 years ago. Previous studies have revealed genomic signals of introgression from wild to cultivated grapes in Europe, but the time, mode, genomic pattern, and biological effects of these introgression events have not been investigated. Here, we studied resequencing data from 345 samples spanning the distributional range of wild (Vitis vinifera ssp. sylvestris) and cultivated (V. vinifera ssp. vinifera) grapes. Based on machine learning-based population genetic analyses, we detected evidence for a single domestication of grapevine, followed by continuous gene flow between European wild grapes (EU) and cultivated grapes over the past ~2,000 y, especially from EU to wine grapes. We also inferred that soft-selective sweeps were the dominant signals of artificial selection. Gene pathways associated with the synthesis of aromatic compounds were enriched in regions that were both selected and introgressed, suggesting EU wild grapes were an important resource for improving the flavor of cultivated grapes. Despite the potential benefits of introgression in grape improvement, the introgressed fragments introduced a higher deleterious burden, with most deleterious SNPs and structural variants hidden in a heterozygous state. Cultivated wine grapes have benefited from adaptive introgression with wild grapes, but introgression has also increased the genetic load. In general, our study of beneficial and harmful effects of introgression is critical for genomic breeding of grapevine to take advantage of wild resources.

Structural variation of a sex-linked region confers monoecy and implicates GATA15 as a master regulator of sex in Salix purpurea

The Salicaceae, including Populus and Salix, are dioecious perennials that utilize different sex determination systems. This family provides a useful system to better understand the evolution of dioecy and sex chromosomes. Here, a rare monoecious genotype of Salix purpurea, 94003, was self- and cross-pollinated and progeny sex ratios were used to test hypotheses on possible mechanisms of sex determination. To delimit genomic regions associated with monoecious expression, the 94003 genome sequence was assembled and DNA- and RNA-Seq of progeny inflorescences was performed. Based on alignments of progeny shotgun DNA sequences to the haplotype-resolved monoecious 94003 genome assembly and reference male and female genomes, a 1.15 Mb sex-linked region on Chr15W was confirmed to be absent in monecious plants. Inheritance of this structural variation is responsible for the loss of a male-suppressing function in what would otherwise be genetic females (ZW), resulting in monoecy (ZW^H or WW^H ), or lethality, if homozygous (W^H W^H ). We present a refined, two-gene sex determination model for Salix purpurea, mediated by ARR17 and GATA15 that is different from the single-gene ARR17-mediated system in the related genus Populus.

The complete reference genome for grapevine (Vitis vinifera L.) genetics and breeding

Grapevine is one of the most economically important crops worldwide. However, the previous versions of the grapevine reference genome tipically consist of thousands of fragments with missing centromeres and telomeres, limiting the accessibility of the repetitive sequences, the centromeric and telomeric regions, and the study of inheritance of important agronomic traits in these regions. Here, we assembled a telomere-to-telomere (T2T) gap-free reference genome for the cultivar PN40024 using PacBio HiFi long reads. The T2T reference genome (PN_T2T) is 69 Mb longer with 9018 more genes identified than the 12X.v0 version. We annotated 67% repetitive sequences, 19 centromeres and 36 telomeres, and incorporated gene annotations of previous versions into the PN_T2T assembly. We detected a total of 377 gene clusters, which showed associations with complex traits, such as aroma and disease resistance. Even though PN40024 derives from nine generations of selfing, we still found nine genomic hotspots of heterozygous sites associated with biological processes, such as the oxidation-reduction process and protein phosphorylation. The fully annotated complete reference genome therefore constitutes an important resource for grapevine genetic studies and breeding programs.

Evolution of sex in crops: recurrent scrap and rebuild

Sexuality is the main strategy for maintaining genetic diversity within a species. In flowering plants (angiosperms), sexuality is derived from ancestral hermaphroditism and multiple sexualities can be expressed in an individual. The mechanisms conferring chromosomal sex determination in plants (or dioecy) have been studied for over a century by both biologists and agricultural scientists, given the importance of this field for crop cultivation and breeding. Despite extensive research, the sex determining gene(s) in plants had not been identified until recently. In this review, we dissect plant sex evolution and determining systems, with a focus on crop species. We introduced classic studies with theoretical, genetic, and cytogenic approaches, as well as more recent research using advanced molecular and genomic techniques. Plants have undergone very frequent transitions into, and out of, dioecy. Although only a few sex determinants have been identified in plants, an integrative viewpoint on their evolutionary trends suggests that recurrent neofunctionalization events are potentially common, in a "scrap and (re)build" cycle. We also discuss the potential association between crop domestication and transitions in sexual systems. We focus on the contribution of duplication events, which are particularly frequent in plant taxa, as a trigger for the creation of new sexual systems.

The chromosome-level genome of Eucommia ulmoides provides insights into sex differentiation and α-linolenic acid biosynthesis

Eucommia ulmoides Oliver is a typical dioecious plant endemic to China that has great medicinal and economic value. Here, we report a high-quality chromosome-level female genome of E. ulmoides obtained by PacBio and Hi-C technologies. The size of the female genome assembly was 1.01 Gb with 17 pseudochromosomes and 31,665 protein coding genes. In addition, Hi-C technology was used to reassemble the male genome released in 2018. The reassembled male genome was 1.24 Gb with the superscaffold N50 (48.30 Mb), which was increased 25.69 times, and the number of predicted genes increased by 11,266. Genome evolution analysis indicated that E. ulmoides has undergone two whole-genome duplication events before the divergence of female and male, including core eudicot γ whole-genome triplication event (γ-WGT) and a recent whole genome duplication (WGD) at approximately 27.3 million years ago (Mya). Based on transcriptome analysis, EuAP3 and EuAG may be the key genes involved in regulating the sex differentiation of E. ulmoides. Pathway analysis showed that the high expression of ω-3 fatty acid desaturase coding gene EU0103017 was an important reason for the high α-linolenic acid content in E. ulmoides. The genome of female and male E. ulmoides presented here is a valuable resource for the molecular biological study of sex differentiation of E. ulmoides and also will provide assistance for the breeding of superior varieties.

Grapevine genome analysis demonstrates the role of gene copy number variation in the formation of monoterpenes

Volatile organic compounds such as terpenes influence the quality parameters of grapevine through their contribution to the flavour and aroma profile of berries. Biosynthesis of volatile organic compounds in grapevine is relatively complex and controlled by multiple genes, the majority of which are unknown or uncharacterised. To identify the genomic regions that associate with modulation of these compounds in grapevine berries, volatile metabolic data generated via GC-MS from a grapevine mapping population was used to identify quantitative trait loci (QTLs). Several significant QTLs were associated with terpenes, and candidate genes were proposed for sesquiterpene and monoterpene biosynthesis. For monoterpenes, loci on chromosomes 12 and 13 were shown to be associated with geraniol and cyclic monoterpene accumulation, respectively. The locus on chromosome 12 was shown to contain a geraniol synthase gene (VvGer), while the locus on chromosome 13 contained an α-terpineol synthase gene (VvTer). Molecular and genomic investigation of VvGer and VvTer revealed that these genes were found in tandemly duplicated clusters, displaying high levels of hemizygosity. Gene copy number analysis further showed that not only did VvTer and VvGer copy numbers vary within the mapping population, but also across recently sequenced Vitis cultivars. Significantly, VvTer copy number correlated with both VvTer gene expression and cyclic monoterpene accumulation in the mapping population. A hypothesis for a hyper-functional VvTer allele linked to increased gene copy number in the mapping population is presented and can potentially lead to selection of cultivars with modulated terpene profiles. The study highlights the impact of VvTPS gene duplication and copy number variation on terpene accumulation in grapevine.

Characterization of Simple Sequence Repeat (SSR) Markers Mined in Whole Grape Genomes

SSR (simple sequence repeat) DNA markers are widely used for genotype DNA identification, QTL mapping, and analyzing genetic biodiversity. However, SSRs in grapes are still in their early stages, with a few primer pairs accessible. With the whole-genome sequencing (WGS) of several grape varieties, characterization of grape SSR changed to be necessary not only to genomics but to also help SSR development and utility. Based on this, we identified the whole-genome SSR of nine grape cultivars (‘PN40024’, ‘Cabernet Sauvignon’, ‘Carménère’, ‘Chardonnay’, ‘Merlot’, ‘Riesling’, ‘Zinfandel’, ‘Shine Muscat’, and ‘Muscat Hamburg’) with whole-genome sequences released publicly and found that there are great differences in the distribution of SSR loci in different varieties. According to the difference in genome size, the number of SSRs ranged from 267,385 (Cabernet Sauvignon) to 627,429 (Carménère), the density of the SSR locus in the genome of nine cultivars was generally 1 per Kb. SSR motif distribution characteristic analysis of these grape cultivars showed that the distribution patterns among grape cultivars were conservative, mainly enriched in A/T. However, there are some differences in motif types (especially tetranucleotides, pentanucleotides, and hexanucleotides), quantity, total length, and average length in different varieties, which might be related to the size of the assembled genome or the specificity of variety domestication. The distribution characteristics of SSRs were revealed by whole-genome analysis of simple repeats of grape varieties. In this study, 32 pairs of primers with lower polymorphism have been screened, which provided an important research foundation for the development of molecular markers of grape variety identification and the construction of linkage maps of important agronomic traits for crop improvement.

Recurrent neo-sex chromosome evolution in kiwifruit

Sex chromosome evolution is thought to be tightly associated with the acquisition and maintenance of sexual dimorphisms. Plant sex chromosomes have evolved independently in many lineages^1,2 and can provide a powerful comparative framework to study this. We assembled and annotated genome sequences of three kiwifruit species (genus Actinidia) and uncovered recurrent sex chromosome turnovers in multiple lineages. Specifically, we observed structural evolution of the neo-Y chromosomes, which was driven via rapid bursts of transposable element insertions. Surprisingly, sexual dimorphisms were conserved in the different species studied, despite the fact that the partially sex-linked genes differ between them. Using gene editing in kiwifruit, we demonstrated that one of the two Y-chromosome-encoded sex-determining genes, Shy Girl, shows pleiotropic effects that can explain the conserved sexual dimorphisms. These plant sex chromosomes therefore maintain sexual dimorphisms through the conservation of a single gene, without a process involving interactions between separate sex-determining genes and genes for sexually dimorphic traits.

Genetic changes in the genus Vitis and the domestication of vine

The genus Vitis belongs to the Vitaceae family and is divided into two subgenera: Muscadinia and Vitis, the main difference between these subgenera being the number of chromosomes. There are many hypotheses about the origin of the genus, which have been formed with archaeological studies and lately with molecular analyses. Even though there is no consensus on the place of origin, these studies have shown that grapes have been used by man since ancient times, starting later on its domestication. Most studies point to the Near East and Greece as the beginning of domestication, current research suggests it took place in parallel in different sites, but in all cases Vitis vinifera (L.) subsp. sylvestris [Vitis vinifera (L.) subsp. sylvestris (Gmelin) Hagi] seems to be the species chosen by our ancestors to give rise to the now known Vitis vinifera (L.) subsp. vinifera [=sativa (Hegi)= caucasica (Vavilov)]. Its evolution and expansion into other territories followed the formation of new empires and their expansion, and this is where the historical importance of this crop lies. In this process, plants with hermaphrodite flowers were preferentially selected, with firmer, sweeter, larger fruits of different colors, thus favoring the selection of genes associated with these traits, also resulting in a change in seed morphology. Currently, genetic improvement programs have made use of wild species for the introgression of disease resistance genes and tolerance to diverse soil and climate environments. In addition, the mapping of genes of interest, both linked to agronomic and fruit quality traits, has allowed the use of molecular markers for assisted selection. Information on the domestication process and genetic resources help to understand the gene pool available for the development of cultivars that respond to producer and consumer requirements.

Chromosome-scale genomics, metabolomics, and transcriptomics provide insight into the synthesis and regulation of phenols in Vitis adenoclada grapes

Vitis adenoclada is a wild grape unique to China. It exhibits well resistance to heat, humidity, fungal disease, drought, and soil infertility. Here, we report the high-quality, chromosome-level genome assembly of GH6 (V. adenoclada). The 498.27 Mb genome contained 221.78 Mb of transposable elements, 28,660 protein-coding genes, and 481.44 Mb of sequences associated with 19 chromosomes. GH6 shares a common ancestor with PN40024 (Vitis vinifera) from approximately 4.26-9.01 million years ago, whose divergence occurred later than Vitis rotundifolia and Vitis riparia. Widely-targeted metabolome and transcriptome analysis revealed that the profiles and metabolism of phenolic compounds in V. adenoclada varieties significantly were differed from other grape varieties. Specifically, V. adenoclada varieties were rich in phenolic acids and flavonols, whereas the flavan-3-ol and anthocyanin content was lower compared with other varieties that have V. vinifera consanguinity in this study. In addition, ferulic acid and stilbenes content were associated with higher expressions of COMT and STSs in V. adenoclada varieties. Furthermore, MYB2, MYB73-1, and MYB73-2 were presumably responsible for the high expression level of COMT in V. adenoclada berries. MYB12 (MYBF1) was positively correlated with PAL, CHS, FLS and UFGT.Meanwhile, MYB4 and MYBC2-L1 may inhibit the synthesis of flavan-3-ols and anthocyanins in two V. adenoclada varieties (YN2 and GH6). The publication of the V. adenoclada grape genome provides a molecular foundation for further revealing its flavor and quality characteristics, is also important for identifying favorable genes of the East Asian species for future breeding.

A multi-omics framework reveals strawberry flavor genes and their regulatory elements

Flavor is essential to consumer preference of foods and is an increasing focus of plant breeding programs. In fruit crops, identifying genes underlying volatile organic compounds has great promise to accelerate flavor improvement, but polyploidy and heterozygosity in many species have slowed progress. Here we use octoploid cultivated strawberry to demonstrate how genomic heterozygosity, transcriptomic intricacy and fruit metabolomic diversity can be treated as strengths and leveraged to uncover fruit flavor genes and their regulatory elements. Multi-omics datasets were generated including an expression quantitative trait loci map with 196 diverse breeding lines, haplotype-phased genomes of a highly-flavored breeding selection, a genome-wide structural variant map using five haplotypes, and volatile genome-wide association study (GWAS) with &gt; 300 individuals. Overlaying regulatory elements, structural variants and GWAS-linked allele-specific expression of numerous genes to variation in volatile compounds important to flavor. In one example, the functional role of anthranilate synthase alpha subunit 1 in methyl anthranilate biosynthesis was supported via fruit transient gene expression assays. These results demonstrate a framework for flavor gene discovery in fruit crops and a pathway to molecular breeding of cultivars with complex and desirable flavor.

HiFi chromosome-scale diploid assemblies of the grape rootstocks 110R, Kober 5BB, and 101-14 Mgt

Cultivated grapevines are commonly grafted on closely related species to cope with specific biotic and abiotic stress conditions. The three North American Vitis species V. riparia, V. rupestris, and V. berlandieri, are the main species used for breeding grape rootstocks. Here, we report the diploid chromosome-scale assembly of three widely used rootstocks derived from these species: Richter 110 (110R), Kober 5BB, and 101–14 Millardet et de Grasset (Mgt). Draft genomes of the three hybrids were assembled using PacBio HiFi sequences at an average coverage of 53.1 X-fold. Using the tool suite HaploSync, we reconstructed the two sets of nineteen chromosome-scale pseudomolecules for each genome with an average haploid genome size of 494.5 Mbp. Residual haplotype switches were resolved using shared-haplotype information. These three reference genomes represent a valuable resource for studying the genetic basis of grape adaption to biotic and abiotic stresses, and designing trait-associated markers for rootstock breeding programs.

Measurement(s)	Genome Assembly Sequence
Technology Type(s)	PacBio Sequel System
Sample Characteristic - Organism	Vitis cinerea var. helleri x Vitis rupestris • Vitis riparia x Vitis rupestris • Vitis cinerea var. helleri x Vitis riparia
Sample Characteristic - Location	State of California

Transcriptomic and metabolomic integration as a resource in grapevine to study fruit metabolite quality traits

Transcriptomics and metabolomics are methodologies being increasingly chosen to perform molecular studies in grapevine (Vitis vinifera L.), focusing either on plant and fruit development or on interaction with abiotic or biotic factors. Currently, the integration of these approaches has become of utmost relevance when studying key plant physiological and metabolic processes. The results from these analyses can undoubtedly be incorporated in breeding programs whereby genes associated with better fruit quality (e.g., those enhancing the accumulation of health-promoting compounds) or with stress resistance (e.g., those regulating beneficial responses to environmental transition) can be used as selection markers in crop improvement programs. Despite the vast amount of data being generated, integrative transcriptome/metabolome meta-analyses (i.e., the joint analysis of several studies) have not yet been fully accomplished in this species, mainly due to particular specificities of metabolomic studies, such as differences in data acquisition (i.e., different compounds being investigated), unappropriated and unstandardized metadata, or simply no deposition of data in public repositories. These meta-analyses require a high computational capacity for data mining a priori, but they also need appropriate tools to explore and visualize the integrated results. This perspective article explores the universe of omics studies conducted in V. vinifera, focusing on fruit-transcriptome and metabolome analyses as leading approaches to understand berry physiology, secondary metabolism, and quality. Moreover, we show how omics data can be integrated in a simple format and offered to the research community as a web resource, giving the chance to inspect potential gene-to-gene and gene-to-metabolite relationships that can later be tested in hypothesis-driven research. In the frame of the activities promoted by the COST Action CA17111 INTEGRAPE, we present the first grapevine transcriptomic and metabolomic integrated database (TransMetaDb) developed within the Vitis Visualization (VitViz) platform (https://tomsbiolab.com/vitviz). This tool also enables the user to conduct and explore meta-analyses utilizing different experiments, therefore hopefully motivating the community to generate Findable, Accessible, Interoperable and Reusable (F.A.I.R.) data to be included in the future.

A recessive mutation in muscadine grapes causes berry color-loss without influencing anthocyanin pathway

Anthocyanins, a major class of flavonoids, are important pigments of grape berries. Despite the recent discovery of the genetic cause underlying the loss of color, the metabolomic and molecular responses are unknown. Anthocyanin quantification among diverse berry color muscadines suggests that all genotypes could produce adequate anthocyanin quantities, irrespective of berry color. Transcriptome profiling of contrasting color muscadine genotypes proposes a potential deficiency that occurs within the anthocyanin transport and/or degradation mechanisms and might cause unpigmented berries. Genome-wide association studies highlighted a region on chromosome-4, comprising several genes encoding glutathione S-transferases involved in anthocyanin transport. Sequence comparison among genotypes reveals the presence of two GST4b alleles that differ by substituting the conserved amino acid residue Pro₁₇₁-to-Leu. Molecular dynamics simulations demonstrate that GST4b2–Leu₁₇₁ encodes an inactive protein due to modifications within the H-binding site. Population genotyping suggests the recessive inheritance of the unpigmented trait with a GST4b2/2 homozygous. A model defining colorless muscadines’ response to the mutation stimulus, avoiding the impact of trapped anthocyanins within the cytoplasm is established.

Transcriptome profiling and mutational analysis suggest a potential deficiency in anthocyanin transport by glutathione S-transferases and/or degradation mechanisms that might cause unpigmented berries.

A convergent mechanism of sex determination in dioecious plants: Distinct sex-determining genes display converged regulation on floral B-class genes

Sex determination in dioecious plants has been broadly and progressively studied with the blooming of genome sequencing and editing techniques. This provides us with a great opportunity to explore the evolution and genetic mechanisms underlining the sex-determining system in dioecious plants. In this study, comprehensively reviewing advances in sex-chromosomes, sex-determining genes, and floral MADS-box genes in dioecious plants, we proposed a convergent model that governs plant dioecy across divergent species using a cascade regulation pathway connecting sex-determining genes and MADS-box genes e.g., B-class genes. We believe that this convergent mechanism of sex determination in dioecious plants will shed light on our understanding of gene regulation and evolution of plant dioecy. Perspectives concerning the evolutionary pathway of plant dioecy are also suggested.

A multi-locus genome-wide association study reveals the genetics underlying muscadine antioxidant in berry skin

Muscadine berries display enhanced nutraceutical value due to the accumulation of distinctive phytochemical constituents with great potential antioxidant activity. Such nutritional and health merits are not only restricted to muscadine, but muscadine berries accumulate higher amounts of bioactive polyphenolics compared with other grape species. For the genetic study of the antioxidant trait in muscadine, a multi-locus genome-wide association study (GWAS) with 350 muscadine genotypes and 1,283 RNase H2 enzyme-dependent amplicon sequencing (rhAmpSeq) markers was performed. Phenotyping was conducted with several antioxidant-related traits, including total phenolic content (TPC), total flavonoid content (TFC), 1,1-diphenyl-2-picrylhydrazyl (DPPH) free radical scavenging activity, and FRAP antioxidant assay in muscadine berry skin. The correlation coefficient analysis revealed that the TPC, and DPPH/FRAP activities were significantly correlated. Through the GWAS analysis, 12 QTNs were identified from the four traits, of which six were pleiotropic QTNs. Two pleiotropic QTNs, chr2_14464718 and chr4_16491374, were commonly identified from the TPC and DPPH/FRAP activities. Co-located genes with the two pleiotropic QTNs were isolated, and two candidate genes were identified with transcriptome analysis. UDP-glycosyltransferase and 4-hydroxy-4-methyl-2-oxoglutarate aldolase were the candidate genes that are positively and negatively correlated to the quantitative property of traits, respectively. These results are the first genetic evidence of the quantitative property of antioxidants in muscadine and provide genetic resources for breeding antioxidant-rich cultivars for both Muscadinia and Euvitis species.

Haplotype-resolved powdery mildew resistance loci reveal the impact of heterozygous structural variation on NLR genes in Muscadinia rotundifolia

Muscadinia rotundifolia cv. Trayshed is a valuable source of resistance to grape powdery mildew. It carries 2 powdery mildew resistance-associated genetic loci, Run1.2 on chromosome 12 and Run2.2 on chromosome 18. The purpose of this study was to identify candidate resistance genes associated with each haplotype of the 2 loci. Both haplotypes of each resistance-associated locus were identified, phased, and reconstructed. Haplotype phasing allowed the identification of several structural variation events between haplotypes of both loci. Combined with a manual refinement of the gene models, we found that the heterozygous structural variants affected the gene content, with some resulting in duplicated or hemizygous nucleotide-binding leucine-rich repeat genes. Heterozygous structural variations were also found to impact the domain composition of some nucleotide-binding leucine-rich repeat proteins. By comparing the nucleotide-binding leucine-rich repeat proteins at Run1.2 and Run2.2 loci, we discovered that the 2 loci include different numbers and classes of nucleotide-binding leucine-rich repeat genes. To identify powdery mildew resistance-associated genes, we performed a gene expression profiling of the nucleotide-binding leucine-rich repeat genes at Run1.2b and Run2.2 loci with or without powdery mildew present. Several nucleotide-binding leucine-rich repeat genes were constitutively expressed, suggesting a role in powdery mildew resistance. These first complete, haplotype-resolved resistance-associated loci and the candidate nucleotide-binding leucine-rich repeat genes identified by this study are new resources that can aid the development of powdery mildew-resistant grape cultivars.

Assembly of complete diploid-phased chromosomes from draft genome sequences

De novo genome assembly is essential for genomic research. High-quality genomes assembled into phased pseudomolecules are challenging to produce and often contain assembly errors because of repeats, heterozygosity, or the chosen assembly strategy. Although algorithms that produce partially phased assemblies exist, haploid draft assemblies that may lack biological information remain favored because they are easier to generate and use. We developed HaploSync, a suite of tools that produces fully phased, chromosome-scale diploid genome assemblies, and performs extensive quality control to limit assembly artifacts. HaploSync scaffolds sequences from a draft diploid assembly into phased pseudomolecules guided by a genetic map and/or the genome of a closely related species. HaploSync generates a report that visualizes the relationships between current and legacy sequences, for both haplotypes, and displays their gene and marker content. This quality control helps the user identify misassemblies and guides Haplosync’s correction of scaffolding errors. Finally, HaploSync fills assembly gaps with unplaced sequences and resolves collapsed homozygous regions. In a series of plant, fungal, and animal kingdom case studies, we demonstrate that HaploSync efficiently increases the assembly contiguity of phased chromosomes, improves completeness by filling gaps, corrects scaffolding, and correctly phases highly heterozygous, complex regions.

Archaeobotanical and chemical investigations on wine amphorae from San Felice Circeo (Italy) shed light on grape beverages at the Roman time

We hereby investigate the pitch used for coating three Roman amphorae from San Felice Circeo (Italy) through a multidisciplinary study. The identification of molecular biomarkers by gas chromatography—mass spectrometry is combined with archaeobotanical evidence of pollen and plant tissues of Vitis flowers. Diterpenic chemical markers together with Pinus pollen and wood revealed Pinaceae tar coating. Aporate 3-zonocolpate pollen, identified as Vitis, together with tartaric, malic and pyruvic acids elucidate the grape-fermented nature of the content. Our conclusions open new consideration on the use of grape derivatives that cannot be supported by traditional analytical methods. Based on the finds of aporate Vitis pollen, found also in local modern and Middle Pleistocene samples, we hypothesize the use of autochthonous vines. The presence of a medicinal wine (historically reported as oenanthium) is also considered. We interrogate Vitis pollen capacity to target grapevine domestication, thereby providing innovative tools to understand such an important process. We anticipate our study to encourage a more systematic multidisciplinary approach regarding the analyses of wine amphorae.

Microstructure of Croatian Wild Grapevine (Vitis vinifera subsp. sylvestris Gmel Hegi) Pollen Grains Revealed by Scanning Electron Microscopy

Wild grapevine (Vitis vinifera subsp. sylvestris Gmel Hegi) is dioecious with male and female plants, whereas domesticated grapevine is mostly hermaphrodite with self-fertile hermaphrodite flowers. The pollen morphology of wild grapevine has been poorly studied. There is no detailed palynological study of V. sylvestris in Croatia and neighboring countries. Here, scanning electron microscopy (SEM) was used to analyze the pollen of V. sylvestris from male and female individuals growing at two natural sites in Croatia. The selective APT3 marker was used to confirm the flower phenotype with the genetic background. SEM analysis showed that the pollen grains of V. sylvestris were isopolar and radially symmetrical, with foveolate perforated ornamentation, regardless of the flower type of the individuals. All male flowers were 3-colporate and prolate in shape, whereas female individuals varied from subprolate to spheroidal and had inaperturate pollen grains. Pollen shape, dimensions and exine ornamentation proved very informative, and here we address the most polymorphic traits in the analyzed V. sylvestris individuals. Principal component analysis (PCA) and clustering based on pollen morphology variables clearly differentiated individuals by their flower type, and no grouping specific to population was observed, pointing to the conserved pollen structure of V. sylvestris. The results indicate the need to continue the palynological study of V. sylvestris and serve as a good phenotypic basis for functional genetic studies on genes involved in pollen morphology and function.

Effects of Different Chemicals on Sexual Regulation in Persimmon (Diospyros kaki Thunb.) Flowers

Research on crop sexuality is important for establishing systems for germplasm innovation and cultivating improved varieties. In this study, androecious persimmon trees were treated with various concentrations of ethrel (100, 500, and 1,000 mg/L) and zeatin (1, 5, and 10 mg/L) to investigate the morphological, physiological, and molecular characteristics of persimmon. Ethrel at 1,000 mg/L and zeatin at 10 mg/L both significantly reduced the stamen length and pollen grain diameter in androecious trees. Ethrel treatment also led to reduced stamen development with degenerated cellular contents; zeatin treatment promoted the development of arrested pistils via maintaining relatively normal mitochondrial morphology. Both treatments altered carbohydrate, amino acid, and endogenous phytohormone contents, as well as genes associated with hormone production and floral organ development. Thereafter, we explored the combined effects of four chemicals, including ethrel and zeatin, as well as zebularine and 5-azacytidine, both of which are DNA methylation inhibitors, on androecious persimmon flower development. Morphological comparisons showed that stamen length, pollen viability, and pollen grain diameter were significantly inhibited after combined treatment. Large numbers of genes involving in carbohydrate metabolic, mitogen-activated protein kinase (MAPK) signaling, and ribosome pathways, and metabolites including uridine monophosphate (UMP) and cyclamic acid were identified in response to the treatment, indicating complex regulatory mechanisms. An association analysis of transcriptomic and metabolomic data indicated that ribosomal genes have distinct effects on UMP and cyclamic acid metabolites, explaining how male floral buds of androecious persimmon trees respond to these exogenous chemicals. These findings extend the knowledge concerning sexual differentiation in persimmon; they also provide a theoretical basis for molecular breeding, high-yield cultivation, and quality improvement in persimmon.

Are plant and animal sex chromosomes really all that different?

Sex chromosomes in plants have often been contrasted with those in animals with the goal of identifying key differences that can be used to elucidate fundamental evolutionary properties. For example, the often homomorphic sex chromosomes in plants have been compared to the highly divergent systems in some animal model systems, such as birds, Drosophila and therian mammals, with many hypotheses offered to explain the apparent dissimilarities, including the younger age of plant sex chromosomes, the lesser prevalence of sexual dimorphism, or the greater extent of haploid selection. Furthermore, many plant sex chromosomes lack complete sex chromosome dosage compensation observed in some animals, including therian mammals, Drosophila, some poeciliids, and Anolis, and plant dosage compensation, where it exists, appears to be incomplete. Even the canonical theoretical models of sex chromosome formation differ somewhat between plants and animals. However, the highly divergent sex chromosomes observed in some animal groups are actually the exception, not the norm, and many animal clades are far more similar to plants in their sex chromosome patterns. This begs the question of how different are plant and animal sex chromosomes, and which of the many unique properties of plants would be expected to affect sex chromosome evolution differently than animals? In fact, plant and animal sex chromosomes exhibit more similarities than differences, and it is not at all clear that they differ in terms of sexual conflict, dosage compensation, or even degree of divergence. Overall, the largest difference between these two groups is the greater potential for haploid selection in plants compared to animals. This may act to accelerate the expansion of the non-recombining region at the same time that it maintains gene function within it. This article is part of the theme issue 'Sex determination and sex chromosome evolution in land plants'.

Recurrent allopolyploidization, Y-chromosome introgression and the evolution of sexual systems in the plant genus Mercurialis

The plant genus Mercurialis includes dioecious, monoecious and androdioecious species (where males coexist with hermaphrodites). Its diversification involved reticulate evolution via hybridization and polyploidization. The Y chromosome of the diploid species Mercurialis annua shows only mild signs of degeneration. We used sequence variation at a Y-linked locus in several species and at multiple autosomal and pseudoautosomal loci to investigate the origin and evolution of the Y chromosome across the genus. Our study provides evidence for further cases of allopolyploid speciation. It also reveals that all lineages with separate sexes (with one possible exception) share the same ancestral Y chromosome. Surprisingly, males in androdioecious populations of hexaploid M. annua carry a Y chromosome that is not derived from either of its two putative progenitor lineages but from a more distantly related perennial dioecious lineage via introgression. These results throw new light on the evolution of sexual systems and polyploidy in Mercurialis and secure it as a promising model for further study of plant sex chromosomes. This article is part of the theme issue 'Sex determination and sex chromosome evolution in land plants'.

Labile sex expression in angiosperm species with sex chromosomes

Here, we review the literature on sexual lability in dioecious angiosperm species with well-studied sex chromosomes. We distinguish three types of departures from strict dioecy, concerning either a minority of flowers in some individuals (leakiness) or the entire individual, which can constantly be bisexual or change sex. We found that for only four of the 22 species studied, reports of lability are lacking. The occurrence of lability is only weakly related to sex chromosome characteristics (number of sex-linked genes, age of the non-recombining region). These results contradict the naive idea that lability is an indication of the absence or the recent evolution of sex chromosomes, and thereby contribute to a growing consensus that sex chromosomes do not necessarily fix sex determination once and for all. We discuss some implications of these findings for the evolution of sex chromosomes, and suggest that more species with well-characterized lability should be studied with genomic data and tools. This article is part of the theme issue 'Sex determination and sex chromosome evolution in land plants'.

The spinach YY genome reveals sex chromosome evolution, domestication, and introgression history of the species

Background

Spinach (Spinacia oleracea L.) is a dioecious species with an XY sex chromosome system, but its Y chromosome has not been fully characterized. Our knowledge about the history of its domestication and improvement remains limited.

Results

A high-quality YY genome of spinach is assembled into 952 Mb in six pseudo-chromosomes. By a combination of genetic mapping, Genome-Wide Association Studies, and genomic analysis, we characterize a 17.42-Mb sex determination region (SDR) on chromosome 1. The sex chromosomes of spinach evolved when an insertion containing sex determination genes occurred, followed by a large genomic inversion about 1.98 Mya. A subsequent burst of SDR-specific repeats (0.1–0.15 Mya) explains the large size of this SDR. We identify a Y-specific gene, NRT1/PTR 6.4 which resides in this insertion, as a strong candidate for the sex determination or differentiation factor. Resequencing of 112 spinach genomes reveals a severe domestication bottleneck approximately 10.87 Kya, which dates the domestication of spinach 7000 years earlier than the archeological record. We demonstrate that a strong selection signal associated with internode elongation and leaf area expansion is associated with domestication of edibility traits in spinach. We find that several strong genomic introgressions from the wild species Spinacia turkestanica and Spinacia tetrandra harbor desirable alleles of genes related to downy mildew resistance, frost resistance, leaf morphology, and flowering-time shift, which likely contribute to spinach improvement.

Conclusions

Analysis of the YY genome uncovers evolutionary forces shaping nascent sex chromosome evolution in spinach. Our findings provide novel insights about the domestication and improvement of spinach.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02633-x.

Berry Anthocyanin, Acid, and Volatile Trait Analyses in a Grapevine-Interspecific F2 Population Using an Integrated GBS and rhAmpSeq Genetic Map

Increased map density and transferability of markers are essential for the genetic analysis of fruit quality and stress tolerance in interspecific grapevine populations. We used 1449 GBS and 2000 rhAmpSeq markers to develop a dense map for an interspecific F₂ population (VRS-F₂) that was derived by selfing a single F₁ from a Vitis riparia x ‘Seyval blanc’ cross. The resultant map contained 2519 markers spanning 1131.3 cM and was highly collinear with the Vitis vinifera ‘PN40024’ genome. Quantitative trait loci (QTL) for berry skin color and flower type were used to validate the map. Four rhAmpSeq transferable markers were identified that can be used in pairs (one pistillate and one hermaphroditic) to predict pistillate and hermaphrodite flower type with ≥99.7% accuracy. Total and individual anthocyanin diglucoside QTL mapped to chromosome 9 near a 5-O-GLUCOSYLTRANSFERASE candidate gene. Malic acid QTL were observed on chromosome 1 and 6 with two MALATE DEHYRDROGENASE CYTOPLASMIC 1 and ALUMINUM-ACTIVATED MALATE TRANSPORTER 2-LIKE (ALMT) candidate genes, respectively. Modeling malic acid identified a potential QTL on chromosome 8 with peak position in proximity of another ALMT. A first-ever reported QTL for the grassy smelling volatile (E)-2-hexenal was found on chromosome 2 with a PHOSPHOLIPID HYDROPEROXIDE GLUTATHIONE PEROXIDASE candidate gene near peak markers.

Reinvention of hermaphroditism via activation of a RADIALIS-like gene in hexaploid persimmon

In flowering plants, different lineages have independently transitioned from the ancestral hermaphroditic state into and out of various sexual systems¹. Polyploidizations are often associated with this plasticity in sexual systems^2,3. Persimmons (the genus Diospyros) have evolved dioecy via lineage-specific palaeoploidizations. More recently, hexaploid D. kaki has established monoecy and also exhibits reversions from male to hermaphrodite flowers in response to natural environmental signals (natural hermaphroditism, NH), or to artificial cytokinin treatment (artificial hermaphroditism, AH). We sought to identify the molecular pathways underlying these polyploid-specific reversions to hermaphroditism. Co-expression network analyses identified regulatory pathways specific to NH or AH transitions. Surprisingly, the two pathways appeared to be antagonistic, with abscisic acid and cytokinin signalling for NH and AH, respectively. Among the genes common to both pathways leading to hermaphroditic flowers, we identified a small-Myb RADIALIS-like gene, named DkRAD, which is specifically activated in hexaploid D. kaki. Consistently, ectopic overexpression of DkRAD in two model plants resulted in hypergrowth of the gynoecium. These results suggest that production of hermaphrodite flowers via polyploidization depends on DkRAD activation, which is not associated with a loss-of-function within the existing sex determination pathway, but rather represents a new path to (or reinvention of) hermaphroditism.

'Nebbiolo' genome assembly allows surveying the occurrence and functional implications of genomic structural variations in grapevines (Vitis vinifera L.)

Background

‘Nebbiolo’ is a grapevine cultivar typical of north-western Italy, appreciated for producing high-quality red wines. Grapevine cultivars are characterized by possessing highly heterozygous genomes, including a great incidence of genomic rearrangements larger than 50 bp, so called structural variations (SVs). Even though abundant, SVs are an under-explored source of genetic variation mainly due to methodological limitations at their detection.

Results

We employed a multiple platform approach to produce long-range genomic data for two different ‘Nebbiolo’ clones, namely: optical mapping, long-reads and linked-reads. We performed a haplotype-resolved de novo assembly for cultivar ‘Nebbiolo’ (clone CVT 71) and used an ab-initio strategy to annotate it. The annotated assembly enhanced our ability to detect SVs, enabling the study of genomic regions not present in the grapevines’ reference genome and accounting for their functional implications. We performed variant calling analyses at three different organizational levels: i) between haplotypes of clone CVT 71 (primary assembly vs haplotigs), ii) between ‘Nebbiolo’ and ‘Cabernet Sauvignon’ assemblies and iii) between clones CVT 71 and CVT 185, representing different ‘Nebbiolo’ biotypes. The cumulative size of non-redundant merged SVs indicated a total of 79.6 Mbp for the first comparison and 136.1 Mbp for the second one, while no SVs were detected for the third comparison. Interestingly, SVs differentiating cultivars and haplotypes affected similar numbers of coding genes.

Conclusions

Our results suggest that SVs accumulation rate and their functional implications in ‘Nebbiolo’ genome are highly-dependent on the organizational level under study. SVs are abundant when comparing ‘Nebbiolo’ to a different cultivar or the two haplotypes of the same individual, while they turned absent between the two analysed clones.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08389-9.

VaAPRT3 Gene is Associated With Sex Determination in Vitis amurensis

In the past decade, progress has been made in sex determination mechanism in Vitis. However, genes responsible for sexual differentiation and its mechanism in V. amurensis remain unknown. Here, we identify a sex determination candidate gene coding adenine phosphoribosyl transferase 3 (VaAPRT3) in V. amurensis. Cloning and sequencing of the VaAPRT3 gene allowed us to develop a molecular marker able to discriminate female individuals from males or hermaphrodites based on a 22-bp InDel. Gene expression and endogenous cytokinin content analysis revealed that the VaAPRT3 gene is involved in sex determination or, to be precise, in female organ differentiation, through regulating cytokinin metabolism in V. amurensis. This study enlarged the understanding of sex determination mechanism in the genus Vitis, and the sex marker could be used as a helpful tool for sexual identification in breeding programs as well as in investigation and collection of V. amurensis germplasms.

Simulated herbivory enhances leaky sex expression in the dioecious herb Mercurialis annua

BACKGROUND AND AIMS

Plant reproductive traits are widely understood to be responsive to the selective pressures exerted by pollinators, but there is also increasing evidence for an important role for antagonists such as herbivores in shaping these traits. Many dioecious species show leaky sex expression, with males and females occasionally producing flowers of the opposite sex. Here, we asked to what extent leakiness in sex expression in Mercurialis annua (Euphorbiaceae) might also be plastically responsive to simulated herbivory. This is important because enhanced leakiness in dioecious populations could lead to a shift in both the mating system and in the conditions for transitions between combined and separate sexes.

METHODS

We examined the effect of simulated herbivory on the sexual expression of males and females of M. annua in two experiments in which different levels of simulated herbivory led to enhanced leakiness in both sexes.

KEY RESULTS

We showed that leaky sex expression in both males and females of the wind-pollinated dioecious herb M. annua is enhanced in response to simulated herbivory, increasing the probability for and the degree of leakiness in both sexes. We also found that leakiness was greater in larger females but not in larger males.

CONCLUSIONS

We discuss hypotheses for a possible functional link between herbivory and leaky sex expression, and consider what simulated herbivory-induced leakiness might imply for the evolutionary ecology of plant reproductive systems, especially the breakdown of dioecy and the evolution of hermaphroditism.

The jojoba genome reveals wide divergence of the sex chromosomes in a dioecious plant

Most flowering plants are hermaphrodites, but around 6% of species are dioecious, having separate male and female plants. Sex chromosomes and some sex-specific genes have been reported in plants, but the genome sequences have not been compared. We now report the genome sequence of male and female jojoba (Simmondsia chinensis) plants, revealing a very large difference in the sex chromosomes. The male genome assembly was 832 Mb and the female 822 Mb. This was explained by the large size differences in the Y chromosome (37.6 Mb) compared with the X chromosome (26.9 Mb). Relative to the X chromosome, the Y chromosome had two large insertions each of more than 5 Mb containing more than 400 genes. Many of the genes in the chromosome-specific regions were novel. These male-specific regions included many flowering-related and stress response genes. Smaller insertions found only in the X chromosome totalled 877 kb. The wide divergence of the sex chromosomes suggests a long period of adaptation to diverging sex-specific roles. Male and female plants may have evolved to accommodate factors such as differing reproductive resource allocation requirements under the stress of the desert environment in which the plants are found. The sex-determining regions accumulate genes beneficial to each sex. This has required the evolution of many more novel sex-specific genes than has been reported for other organisms. This suggest that dioecious plants provide a novel source of genes for manipulation of reproductive performance and environmental adaptation in crops.