The genetic basis of sex determination in grapes

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.

VviPLATZ1 is a major factor that controls female flower morphology determination in grapevine

Plant genetic sex determinants that mediate the transition to dioecy are predicted to be diverse, as this type of mating system independently evolved multiple times in angiosperms. Wild Vitis species are dioecious with individuals producing morphologically distinct female or male flowers; whereas, modern domesticated Vitis vinifera cultivars form hermaphrodite flowers capable of self-pollination. Here, we identify the VviPLATZ1 transcription factor as a key candidate female flower morphology factor that localizes to the Vitis SEX-DETERMINING REGION. The expression pattern of this gene correlates with the formation reflex stamens, a prominent morphological phenotype of female flowers. After generating CRISPR/Cas9 gene-edited alleles in a hermaphrodite genotype, phenotype analysis shows that individual homozygous lines produce flowers with reflex stamens. Taken together, our results demonstrate that loss of VviPLATZ1 function is a major factor that controls female flower morphology in Vitis.

Unlike wild Vitis species, which produce either female or male flowers, modern grapevine cultivars form hermaphrodite flowers for self-pollination. Here, the authors report that the VviPLATZ1 (plant AT-rich sequence-and zinc-binding protein1) transcription factor functions in controlling female flower morphology determination.

AGAMOUS Gene as a New Sex-Identification Marker in Fig (Ficus carica L.) Is More Efficient Than RAN1

Fig is an ancient gynodioecious fruit tree with females for commercial fruit production and hermaphrodites (males) sometimes used as pollen providers. An early sex-identification method would improve breeding efficiency. Three AGAMOUS (AG) genes were recruited from the Ficus carica genome using AG sequences from Ficus microcarpa and Ficus hispida. FcAG was 5230 bp in length, with 7 exons and 6 introns, and a 744-bp coding sequence. The gene was present in both female and male fig genomes, with a 15-bp deletion in the 7th exon. The other two AG genes (FcAG2-Gall_Stamen and FcAG3-Gall_Stamen) were male-specific, without the 15-bp deletion (759-bp coding sequence), and were only expressed in the gall and stamen of the male fig fruit. Using the deletion as the forward primer (AG-Marker), male plants were very efficiently identified by the presence of a 146-bp PCR product. The previously reported fig male and female polymorphism gene RESPONSIVE-TO-ANTAGONIST1 (RAN1) was also cloned and compared between male and female plants. Fifteen SNPs were found in the 3015-bp protein-coding sequence. Among them, 12 SNPs were identified as having sex-differentiating capacity by checking the sequences of 27 known male and 24 known female cultivars. A RAN1-Marker of 608 bp, including 6 SNPs, was designed, and a PCR and sequencing-based method was verified with 352 fig seedlings from two hybrid populations. Our results confirmed that the newly established AG-Marker is as accurate as the RAN1-Marker, and provide new clues to understanding Ficus sex determination.

Horizontal Transfer of LTR Retrotransposons Contributes to the Genome Diversity of Vitis

While horizontally transferred transposable elements (TEs) have been reported in several groups of plants, their importance for genome evolution remains poorly understood. To understand how horizontally transferred TEs contribute to plant genome evolution, we investigated the composition and activity of horizontally transferred TEs in the genomes of four Vitis species. A total of 35 horizontal transfer (HT) events were identified between the four Vitis species and 21 other plant species belonging to 14 different families. We determined the donor and recipient species for 28 of these HTs, with the Vitis species being recipients of 15 of them. As a result of HTs, 8–10 LTR retrotransposon clusters were newly formed in the genomes of the four Vitis species. The activities of the horizontally acquired LTR retrotransposons differed among Vitis species, showing that the consequences of HTs vary during the diversification of the recipient lineage. Our study provides the first evidence that the HT of TEs contributes to the diversification of plant genomes by generating additional TE subfamilies and causing their differential proliferation in host genomes.

A tetraploidization event shaped the Aquilaria sinensis genome and contributed to the ability of sesquiterpenes synthesis

BACKGROUND

Agarwood, generated from the Aquilaria sinensis, has high economic and medicinal value. Although its genome has been sequenced, the ploidy of A. sinensis paleopolyploid remains unclear. Moreover, the expression changes of genes associated with agarwood formation were not analyzed either.

RESULTS

In the present work, we reanalyzed the genome of A. sinensis and found that it experienced a recent tetraploidization event ~ 63-71 million years ago (Mya). The results also demonstrated that the A. sinensis genome had suffered extensive gene deletion or relocation after the tetraploidization event, and exhibited accelerated evolutionary rates. At the same time, an alignment of homologous genes related to different events of polyploidization and speciation were generated as well, which provides an important comparative genomics resource for Thymelaeaceae and related families. Interestingly, the expression changes of genes related to sesquiterpene synthesis in wounded stems of A. sinensis were also observed. Further analysis demonstrated that polyploidization promotes the functional differentiation of the key genes in the sesquiterpene synthesis pathway.

CONCLUSIONS

By reanalyzing its genome, we found that the tetraploidization event shaped the A. sinensis genome and contributed to the ability of sesquiterpenes synthesis. We hope that these results will facilitate our understanding of the evolution of A. sinensis and the function of genes involved in agarwood formation.

Introgression among North American wild grapes (Vitis) fuels biotic and abiotic adaptation

BACKGROUND

Introgressive hybridization can reassort genetic variants into beneficial combinations, permitting adaptation to new ecological niches. To evaluate evolutionary patterns and dynamics that contribute to introgression, we investigate six wild Vitis species that are native to the Southwestern United States and useful for breeding grapevine (V. vinifera) rootstocks.

RESULTS

By creating a reference genome assembly from one wild species, V. arizonica, and by resequencing 130 accessions, we focus on identifying putatively introgressed regions (pIRs) between species. We find six species pairs with signals of introgression between them, comprising up to ~ 8% of the extant genome for some pairs. The pIRs tend to be gene poor, located in regions of high recombination and enriched for genes implicated in disease resistance functions. To assess potential pIR function, we explore SNP associations to bioclimatic variables and to bacterial levels after infection with the causative agent of Pierce's disease (Xylella fastidiosa). pIRs are enriched for SNPs associated with both climate and bacterial levels, suggesting that introgression is driven by adaptation to biotic and abiotic stressors.

CONCLUSIONS

Altogether, this study yields insights into the genomic extent of introgression, potential pressures that shape adaptive introgression, and the evolutionary history of economically important wild relatives of a critical crop.

Transcripts switched off at the stop of phloem unloading highlight the energy efficiency of sugar import in the ripening V. vinifera fruit

Transcriptomic changes at the cessation of sugar accumulation in the pericarp of Vitis vinifera were addressed on single berries re-synchronised according to their individual growth patterns. The net rates of water, sugars and K⁺ accumulation inferred from individual growth and solute concentration confirmed that these inflows stopped simultaneously in the ripe berry, while the small amount of malic acid remaining at this stage was still being oxidised at low rate. Re-synchronised individual berries displayed negligible variations in gene expression among triplicates. RNA-seq studies revealed sharp reprogramming of cell-wall enzymes and structural proteins at the stop of phloem unloading, associated with an 80% repression of multiple sugar transporters and aquaporins on the plasma or tonoplast membranes, with the noticeable exception of H⁺/sugar symporters, which were rather weakly and constitutively expressed. This was verified in three genotypes placed in contrasted thermo-hydric conditions. The prevalence of SWEET suggests that electrogenic transporters would play a minor role on the plasma membranes of SE/CC complex and the one of the flesh, while sucrose/H⁺ exchangers dominate on its tonoplast. Cis-regulatory elements present in their promoters allowed to sort these transporters in different groups, also including specific TIPs and PIPs paralogs, and cohorts of cell wall-related genes. Together with simple thermodynamic considerations, these results lead to propose that H⁺/sugar exchangers at the tonoplast, associated with a considerably acidic vacuolar pH, may exhaust cytosolic sugars in the flesh and alleviate the need for supplementary energisation of sugar transport at the plasma membrane.

Rootstock influences the effect of grapevine leafroll-associated viruses on berry development and metabolism via abscisic acid signalling

Grapevine leafroll-associated virus (GLRaV) infections are accompanied by symptoms influenced by host genotype, rootstock, environment, and which individual or combination of GLRaVs is present. Using a dedicated experimental vineyard, we studied the responses to GLRaVs in ripening berries from Cabernet Franc grapevines grafted to different rootstocks and with zero, one, or pairs of leafroll infection(s). RNA sequencing data were mapped to a high-quality Cabernet Franc genome reference assembled to carry out this study and integrated with hormone and metabolite abundance data. This study characterized conserved and condition-dependent responses to GLRaV infection(s). Common responses to GLRaVs were reproduced in two consecutive years and occurred in plants grafted to different rootstocks in more than one infection condition. Though different infections were inconsistently distinguishable from one another, the effects of infections in plants grafted to different rootstocks were distinct at each developmental stage. Conserved responses included the modulation of genes related to pathogen detection, abscisic acid (ABA) signalling, phenylpropanoid biosynthesis, and cytoskeleton remodelling. ABA, ABA glucose ester, ABA and hormone signalling-related gene expression, and the expression of genes in several transcription factor families differentiated the effects of GLRaVs in berries from Cabernet Franc grapevines grafted to different rootstocks. These results support that ABA participates in the shared responses to GLRaV infection and differentiates the responses observed in grapevines grafted to different rootstocks.

Plant genera Cannabis and Humulus share the same pair of well-differentiated sex chromosomes

We recently described, in Cannabis sativa, the oldest sex chromosome system documented so far in plants (12-28 Myr old). Based on the estimated age, we predicted that it should be shared by its sister genus Humulus, which is known also to possess XY chromosomes. Here, we used transcriptome sequencing of an F₁ family of H. lupulus to identify and study the sex chromosomes in this species using the probabilistic method SEX-DETector. We identified 265 sex-linked genes in H. lupulus, which preferentially mapped to the C. sativa X chromosome. Using phylogenies of sex-linked genes, we showed that a region of the sex chromosomes had already stopped recombining in an ancestor of both species. Furthermore, as in C. sativa, Y-linked gene expression reduction is correlated to the position on the X chromosome, and highly Y degenerated genes showed dosage compensation. We report, for the first time in Angiosperms, a sex chromosome system that is shared by two different genera. Thus, recombination suppression started at least 21-25 Myr ago, and then (either gradually or step-wise) spread to a large part of the sex chromosomes (c. 70%), leading to a degenerated Y chromosome.

Genome-Wide Association Studies for Sex Determination and Cross-Compatibility in Water Yam (Dioscorea alata L.)

Yam (Dioscorea spp.) species are predominantly dioecious, with male and female flowers borne on separate individuals. Cross-pollination is, therefore, essential for gene flow among and within yam species to achieve breeding objectives. Understanding genetic mechanisms underlying sex determination and cross-compatibility is crucial for planning a successful hybridization program. This study used the genome-wide association study (GWAS) approach for identifying genomic regions linked to sex and cross-compatibility in water yam (Dioscorea alata L.). We identified 54 markers linked to flower sex determination, among which 53 markers were on chromosome 6 and one on chromosome 11. Our result ascertained that D. alata is characterized by the male heterogametic sex determination system (XX/XY). The cross-compatibility indices, average crossability rate (ACR) and percentage high crossability (PHC), were controlled by loci on chromosomes 1, 6 and 17. Of the significant loci, SNPs located on chromosomes 1 and 17 were the most promising for ACR and PHC, respectively, and should be validated for use in D. alata hybridization activities to predict cross-compatibility success. A total of 61 putative gene/protein families with direct or indirect influence on plant reproduction were annotated in chromosomic regions controlling the target traits. This study provides valuable insights into the genetic control of D. alata sexual reproduction. It opens an avenue for developing genomic tools for predicting hybridization success in water yam breeding programs.

Diploid chromosome-scale assembly of the Muscadinia rotundifolia genome supports chromosome fusion and disease resistance gene expansion during Vitis and Muscadinia divergence

Muscadinia rotundifolia, the muscadine grape, has been cultivated for centuries in the southeastern United States. M. rotundifolia is resistant to many of the pathogens that detrimentally affect Vitis vinifera, the grape species commonly used for winemaking. For this reason, M. rotundifolia is a valuable genetic resource for breeding. Single-molecule real-time reads were combined with optical maps to reconstruct the two haplotypes of each of the 20 M. rotundifolia cv. Trayshed chromosomes. The completeness and accuracy of the assembly were confirmed using a high-density linkage map. Protein-coding genes were annotated using an integrated and comprehensive approach. This included using full-length cDNA sequencing (Iso-Seq) to improve gene structure and hypothetical spliced variant predictions. Our data strongly support that Muscadinia chromosomes 7 and 20 are fused in Vitis and pinpoint the location of the fusion in Cabernet Sauvignon and PN40024 chromosome 7. Disease-related gene numbers in Trayshed and Cabernet Sauvignon were similar, but their clustering locations were different. A dramatic expansion of the Toll/Interleukin-1 Receptor-like Nucleotide-Binding Site Leucine-Rich Repeat (TIR-NBS-LRR) class was detected on Trayshed chromosome 12 at the Resistance to Uncinula necator 1 (RUN1)/Resistance to Plasmopara viticola 1 (RPV1) locus, which confers strong dominant resistance to powdery and downy mildews. A genome browser, annotation, and Blast tool for Trayshed are available at www.grapegenomics.com.

Population genomics of apricots unravels domestication history and adaptive events

Among crop fruit trees, the apricot (Prunus armeniaca) provides an excellent model to study divergence and adaptation processes. Here, we obtain nearly 600 Armeniaca apricot genomes and four high-quality assemblies anchored on genetic maps. Chinese and European apricots form two differentiated gene pools with high genetic diversity, resulting from independent domestication events from distinct wild Central Asian populations, and with subsequent gene flow. A relatively low proportion of the genome is affected by selection. Different genomic regions show footprints of selection in European and Chinese cultivated apricots, despite convergent phenotypic traits, with predicted functions in both groups involved in the perennial life cycle, fruit quality and disease resistance. Selection footprints appear more abundant in European apricots, with a hotspot on chromosome 4, while admixture is more pervasive in Chinese cultivated apricots. Our study provides clues to the biology of selected traits and targets for fruit tree research and breeding.

Three genomes in the algal genus Volvox reveal the fate of a haploid sex-determining region after a transition to homothallism

Evolutionary transitions between species with separate sexes and species in which individuals have both sex functions have wide-ranging biological implications. It is largely unknown how such transitions occur in systems with haploid male- and female-determining chromosomes in algae and bryophytes. We investigated such a transition in the algal genus Volvox by making whole-genome sequences of two closely related species, one of which is heterothallic (with distinct males and females) and the other homothallic (with only bisexual, self-compatible individuals). The heterothallic species harbors a sex-determining region (SDR), while the homothallic species retains a nearly intact female-derived SDR-like region and separate regions containing key male genes. Thus, an ancestral female has probably become homothallic by acquiring genes that confer male functions.

Transitions between separate sexes (dioecy) and other mating systems are common across eukaryotes. Here, we study a change in a haploid dioecious green algal species with male- and female-determining chromosomes (U and V). The genus Volvox is an oogamous (with large, immotile female gametes and small, motile male gametes) and includes both heterothallic species (with distinct male and female genotypes, associated with a mating-type system that prevents fusion of gametes of the same sex) and homothallic species (bisexual, with the ability to self-fertilize). We date the origin of an expanded sex-determining region (SDR) in Volvox to at least 75 Mya, suggesting that homothallism represents a breakdown of dioecy (heterothallism). We investigated the involvement of the SDR of the U and V chromosomes in this transition. Using de novo whole-genome sequences, we identified a heteromorphic SDR of ca 1 Mbp in male and female genotypes of the heterothallic species Volvox reticuliferus and a homologous region (SDLR) in the closely related homothallic species Volvox africanus, which retained several different hallmark features of an SDR. The V. africanus SDLR includes a large region resembling the female SDR of the presumptive heterothallic ancestor, whereas most genes from the male SDR are absent. However, we found a multicopy array of the male-determining gene, MID, in a different genomic location from the SDLR. Thus, in V. africanus, an ancestrally female genotype may have acquired MID and thereby gained male traits.

Back to the Origins: Background and Perspectives of Grapevine Domestication

Domestication is a process of selection driven by humans, transforming wild progenitors into domesticated crops. The grapevine (Vitis vinifera L.), besides being one of the most extensively cultivated fruit trees in the world, is also a fascinating subject for evolutionary studies. The domestication process started in the Near East and the varieties obtained were successively spread and cultivated in different areas. Whether the domestication occurred only once, or whether successive domestication events occurred independently, is a highly debated mystery. Moreover, introgression events, breeding and intense trade in the Mediterranean basin have followed, in the last thousands of years, obfuscating the genetic relationships. Although a succession of studies has been carried out to explore grapevine origin and different evolution models are proposed, an overview of the topic remains pending. We review here the findings obtained in the main phylogenetic and genomic studies proposed in the last two decades, to clarify the fundamental questions regarding where, when and how many times grapevine domestication took place. Finally, we argue that the realization of the pan-genome of grapes could be a useful resource to discover and track the changes which have occurred in the genomes and to improve our understanding about the domestication.

Epigenetics drive the evolution of sex chromosomes in animals and plants

We review how epigenetics affect sex chromosome evolution in animals and plants. In a few species, sex is determined epigenetically through the action of Y-encoded small RNAs. Epigenetics is also responsible for changing the sex of individuals through time, even in species that carry sex chromosomes, and could favour species adaptation through breeding system plasticity. The Y chromosome accumulates repeats that become epigenetically silenced which leads to an epigenetic conflict with the expression of Y genes and could accelerate Y degeneration. Y heterochromatin can be lost through ageing, which activates transposable elements and lowers male longevity. Y chromosome degeneration has led to the evolution of meiotic sex chromosome inactivation in eutherians (placentals) and marsupials, and dosage compensation mechanisms in animals and plants. X-inactivation convergently evolved in eutherians and marsupials via two independently evolved non-coding RNAs. In Drosophila, male X upregulation by the male specific lethal (MSL) complex can spread to neo-X chromosomes through the transposition of transposable elements that carry an MSL-binding motif. We discuss similarities and possible differences between plants and animals and suggest future directions for this dynamic field of research. This article is part of the theme issue 'How does epigenetics influence the course of evolution?'

Multiple independent recombinations led to hermaphroditism in grapevine

We studied the grape sex-determining region (SDR) in 12 Vitis genomes and demonstrated its conservation across 556 genotypes including 193 accessions from 47 world-wide wild grapevine species and 363 accessions of cultivated grapevine. Although the grape SDR is recombination free in all wild species, we found two distinct hermaphrodite (H) haplotypes (H1 and H2) among the cultivated grapevines, both chimeras of male (M) and female (f) haplotypes. The two independent recombinations carry different genetic signatures which long predate the domestication of grapevine, suggesting independent evolutions of this trait in wild European grapevine gene pools prior to human domestication.

Hermaphroditic (perfect) flowers were a key trait in grapevine domestication, enabling a drastic increase in yields due to the efficiency of self-pollination in the domesticated grapevine (Vitis vinifera L. ssp. vinifera). In contrast, all extant wild Vitis species are dioecious, each plant having only male or female flowers. In this study, we identified the male (M) and female (f) haplotypes of the sex-determining region (SDR) in the wild grapevine species V. cinerea and confirmed the boundaries of the SDR. We also demonstrated that the SDR and its boundaries are precisely conserved across the Vitis genus using shotgun resequencing data of 556 wild and domesticated accessions from North America, East Asia, and Europe. A high linkage disequilibrium was found at the SDR in all wild grape species, while different recombination signatures were observed along the hermaphrodite (H) haplotype of 363 cultivated accessions, revealing two distinct H haplotypes, named H1 and H2. To further examine the H2 haplotype, we sequenced the genome of two grapevine cultivars, 'Riesling' and 'Chardonnay'. By reconstructing the first two H2 haplotypes, we estimated the divergence time between H1 and H2 haplotypes at ∼6 million years ago, which predates the domestication of grapevine (∼8,000 y ago). Our findings emphasize the important role of recombination suppression in maintaining dioecy in wild grape species and lend additional support to the hypothesis that at least two independent recombination events led to the reversion to hermaphroditism in grapevine.

Genome-Wide Identification of Sequence Variations and SSR Marker Development in the Munake Grape Cultivar

The Munake grape cultivar produces uniquely flavored and high-quality fruits. Despite the numerous beneficial agronomic traits of Munake, there are few genetic resources available for this cultivar. To address this knowledge gap, the entire genome was sequenced using whole-genome sequencing approaches and compared with a Vitis vinifera L. reference genome. This study describes more than 3 million single nucleotide polymorphism (SNP), 300,000 insertion and deletion (InDel), 14,000 structural variation (SV), and 80,000 simple sequence repeat (SSR) markers (one SSR per 4.23 kb), as well as their primers. Among the SSRs, 44 SSR primer pairs were randomly selected and validated by polymerase chain reaction (PCR), allowing discrimination between the different Munake cultivar genotypes. The genetic data provided allow a deeper understanding of Munake cultivar genomic sequence and contribute to better knowledge of the genetic basis behind its key agronomic traits.

Plant sex chromosomes defy evolutionary models of expanding recombination suppression and genetic degeneration

Hundreds of land plant lineages have independently evolved separate sexes in either gametophytes (dioicy) or sporophytes (dioecy), but 43% of all dioecious angiosperms are found in just 34 entirely dioecious clades, suggesting that their mode of sex determination evolved a long time ago. Here, we review recent insights on the molecular mechanisms that underlie the evolutionary change from individuals that each produce male and female gametes to individuals specializing in the production of just one type of gamete. The canonical model of sex chromosome evolution in plants predicts that two sex-determining genes will become linked in a sex-determining region (SDR), followed by expanding recombination suppression, chromosome differentiation and, ultimately, degeneration. Experimental work, however, is showing that single genes function as master regulators in model systems, such as the liverwort Marchantia and the angiosperms Diospyros and Populus. In Populus, this type of regulatory function has been demonstrated by genome editing. In other systems, including Actinidia, Asparagus and Vitis, two coinherited factors appear to independently regulate female and male function, yet sex chromosome differentiation has remained low. We discuss the best-understood systems and evolutionary pathways to dioecy, and present a meta-analysis of the sizes and ages of SDRs. We propose that limited sexual conflict explains why most SDRs are small and sex chromosomes remain homomorphic. It appears that models of increasing recombination suppression with age do not apply because selection favours mechanisms in which sex determination depends on minimal differences, keeping it surgically precise.

The Diversity of Plant Sex Chromosomes Highlighted through Advances in Genome Sequencing

For centuries, scientists have been intrigued by the origin of dioecy in plants, characterizing sex-specific development, uncovering cytological differences between the sexes, and developing theoretical models. Through the invention and continued improvements in genomic technologies, we have truly begun to unlock the genetic basis of dioecy in many species. Here we broadly review the advances in research on dioecy and sex chromosomes. We start by first discussing the early works that built the foundation for current studies and the advances in genome sequencing that have facilitated more-recent findings. We next discuss the analyses of sex chromosomes and sex-determination genes uncovered by genome sequencing. We synthesize these results to find some patterns are emerging, such as the role of duplications, the involvement of hormones in sex-determination, and support for the two-locus model for the origin of dioecy. Though across systems, there are also many novel insights into how sex chromosomes evolve, including different sex-determining genes and routes to suppressed recombination. We propose the future of research in plant sex chromosomes should involve interdisciplinary approaches, combining cutting-edge technologies with the classics to unravel the patterns that can be found across the hundreds of independent origins.

The genome of Shanputao (Vitis amurensis) provides a new insight into cold tolerance of grapevine

Vitis amurensis (Shanputao) is the most cold tolerant Vitis species and so is of great interest to grape breeders and producers in areas with low winter temperatures. Here, we report its high-quality, chromosome-level genome assembly based on a combination of sequence data from Illumina and PacBio platforms, BioNano optical mapping and high-throughput chromosome conformation Capture (Hi-C) mapping. The 604.56-Mb genome contains 32 885 protein-coding genes. Shanputao was found to share a common ancestor with PN40024 (V. vinifera) approximately 2.17-2.91 million years ago, and gene expansion observed in Shanputao might contribute to the enhancement of cold tolerance. Transcriptome analysis revealed 17 genes involved in cold signal transduction, suggesting that there was a different response mechanism to chilling temperature and freezing conditions. Furthermore, a genome-wide association study uncovered a phosphoglycerate kinase gene that may contribute to the freezing resistance of buds in the winter. The Shanputao genome sequence not only represents a valuable resource for grape breeders, but also is important for clarifying the molecular mechanisms involved in cold tolerance.

Sharing the Agrarian Knowledge with Archaeology: First Evidence of the Dimorphism of Vitis Pollen from the Middle Bronze Age of N Italy (Terramara Santa Rosa di Poviglio)

The recovery of inaperturate pollen from functionally female flowers in archaeological layers opens the question of a possible pollen-based discrimination between wild and domesticated Vitis vinifera in prehistoric times. Pollen analysis applied to archaeology has not routinely considered the existence of pollen dimorphism in Vitis, a well-known trait in the field of agrarian studies. Therefore, the inaperturate shape of grapevine pollen is ignored by studies on the archaeobotanical history of viticulture. In this paper we investigate pollen morphology of the domesticated and wild subspecies of V. vinifera, and report the first evidence of inaperturate Vitis pollen from an archaeological site. We studied exemplar cases of plants with hermaphroditic flowers, belonging to the subspecies vinifera with fully developed male and female organs, cases of dioecious plants with male or female flowers, belonging to the wild subspecies sylvestris and cases of V. vinifera subsp. vinifera with morphologically hermaphroditic but functionally female flowers. The pollen produced by hermaphroditic and male flowers is usually trizonocolporate; the pollen produced by female flowers is inaperturate. This paper reports on the inaperturate pollen of Vitis found in an archeological site of the Po Plain, Northern Italy. The site dated to the Bronze Age, which is known to have been a critical age for the use of this plant with a transition from wild to domesticated Vitis in central Mediterranean. Can the inaperturate Vitis pollen be a marker of wild Vitis vinifera in prehistoric times? Palynology suggests a possible new investigation strategy on the ancient history of the wild and cultivated grapevine. The pollen dimorphism also implies a different production and dispersal of pollen of the wild and the domesticated subspecies. Grapevine plants are palynologically different from the other Mediterranean “cultural trees”. In fact, Olea, Juglans and Castanea, which are included in the OJC index, have the same pollen morphology and the same pollen dispersal, in wild and domesticated plants. In contrast, the signal of Vitis pollen in past records may be different depending on the hermaphroditic or dioecious subspecies.

The Diversity and Dynamics of Sex Determination in Dioecious Plants

The diversity of inflorescences among flowering plants is captivating. Such charm is not only due to the variety of sizes, shapes, colors, and flowers displayed, but also to the range of reproductive systems. For instance, hermaphrodites occur abundantly throughout the plant kingdom with both stamens and carpels within the same flower. Nevertheless, 10% of flowering plants have separate unisexual flowers, either in different locations of the same individual (monoecy) or on different individuals (dioecy). Despite their rarity, dioecious plants provide an excellent opportunity to investigate the mechanisms involved in sex expression and the evolution of sex-determining regions (SDRs) and sex chromosomes. The SDRs and the evolution of dioecy have been studied in many species ranging from Ginkgo to important fruit crops. Some of these studies, for example in asparagus or kiwifruit, identified two sex-determining genes within the non-recombining SDR and may thus be consistent with the classical model for the evolution of dioecy from hermaphroditism via gynodioecy, that predicts two successive mutations, the first one affecting male and the second one female function, becoming linked in a region of suppressed recombination. On the other hand, aided by genome sequencing and gene editing, single factor sex determination has emerged in other species, such as persimmon or poplar. Despite the diversity of sex-determining mechanisms, a tentative comparative analysis of the known sex-determining genes and candidates in different species suggests that similar genes and pathways may be employed repeatedly for the evolution of dioecy. The cytokinin signaling pathway appears important for sex determination in several species regardless of the underlying genetic system. Additionally, tapetum-related genes often seem to act as male-promoting factors when sex is determined via two genes. We present a unified model that synthesizes the genetic networks of sex determination in monoecious and dioecious plants and will support the generation of hypothesis regarding candidate sex determinants in future studies.

HapSolo: an optimization approach for removing secondary haplotigs during diploid genome assembly and scaffolding

BACKGROUND

Despite marked recent improvements in long-read sequencing technology, the assembly of diploid genomes remains a difficult task. A major obstacle is distinguishing between alternative contigs that represent highly heterozygous regions. If primary and secondary contigs are not properly identified, the primary assembly will overrepresent both the size and complexity of the genome, which complicates downstream analysis such as scaffolding.

RESULTS

Here we illustrate a new method, which we call HapSolo, that identifies secondary contigs and defines a primary assembly based on multiple pairwise contig alignment metrics. HapSolo evaluates candidate primary assemblies using BUSCO scores and then distinguishes among candidate assemblies using a cost function. The cost function can be defined by the user but by default considers the number of missing, duplicated and single BUSCO genes within the assembly. HapSolo performs hill climbing to minimize cost over thousands of candidate assemblies. We illustrate the performance of HapSolo on genome data from three species: the Chardonnay grape (Vitis vinifera), with a genome of 490 Mb, a mosquito (Anopheles funestus; 200 Mb) and the Thorny Skate (Amblyraja radiata; 2650 Mb).

CONCLUSIONS

HapSolo rapidly identified candidate assemblies that yield improvements in assembly metrics, including decreased genome size and improved N50 scores. Contig N50 scores improved by 35%, 9% and 9% for Chardonnay, mosquito and the thorny skate, respectively, relative to unreduced primary assemblies. The benefits of HapSolo were amplified by down-stream analyses, which we illustrated by scaffolding with Hi-C data. We found, for example, that prior to the application of HapSolo, only 52% of the Chardonnay genome was captured in the largest 19 scaffolds, corresponding to the number of chromosomes. After the application of HapSolo, this value increased to ~ 84%. The improvements for the mosquito's largest three scaffolds, representing the number of chromosomes, were from 61 to 86%, and the improvement was even more pronounced for thorny skate. We compared the scaffolding results to assemblies that were based on PurgeDups for identifying secondary contigs, with generally superior results for HapSolo.

Construction of high-density genetic maps defined sex determination region of the Y chromosome in spinach

Spinach (Spinacia olracea L.) is a dioecious leafy vegetable with a highly repetitive genome of around 990 Mb, which is challenging for de-novo genome assembly. In our study, a segregating F₁ (double pseudo-testcross) population from 'Viroflay' × 'Cornell-NO. 9' was used for genetic mapping by resequencing genotyping. In the paternal 'Cornell-NO. 9' map, 212,414 SNPs were mapped, and the total linkage distance was 476.83 cM; the maternal 'Viroflay' map included 29,282 SNPs with 401.28 cM total genetic distance. Both paternal and maternal maps have the expected number of six linkage groups (LGs). A non-recombining region with 5678 SNPs (39 bin markers) co-segregates with sex type which located at 45.2 cM of LG1 in the 'Cornell-NO. 9' map while indicates the sex determination region (SDR). Integration of two maps into a consensus map guided us to anchor additional 1242 contigs to six pseudomolecules from the published reference genome, which improved additional 233 Mb (23.4%) assembly based on spinach estimated genome size. Particularly, the X counterpart of SDR in our assembly is estimated around 18.4 Mb which locates at the largest chromosome, as consensus with sex-biased FISH signals from previous cytogenetics studies. The region is featured by reduced gene density, higher percentage of repetitive sequences, and no recombination. Our linkage maps provide the resource for improving spinach genome de-novo assembly and identification of sex-determining genes in spinach.

The de novo genome assembly of Tapiscia sinensis and the transcriptomic and developmental bases of androdioecy

Tapiscia sinensis (Tapisciaceae) possesses an unusual androdioecious breeding system that has attracted considerable interest from evolutionary biologists. Key aspects of T. sinensis biology, including its biogeography, genomics, and sex-linked genes, are unknown. Here, we report the first de novo assembly of the genome of T. sinensis. The genome size was 410 Mb, with 22,251 predicted genes. Based on whole-genome resequencing of 55 trees from 10 locations, an analysis of population genetic structure indicated that T. sinensis has fragmented into five lineages, with low intrapopulation genetic diversity and little gene flow among populations. By comparing whole-genome scans of male versus hermaphroditic pools, we identified 303 candidate sex-linked genes, 79 of which (25.9%) were located on scaffold 25. A 24-kb region was absent in hermaphroditic individuals, and five genes in that region, TsF-box4, TsF-box10, TsF-box13, TsSUT1, and TsSUT4, showed expression differences between mature male and hermaphroditic flowers. The results of this study shed light on the breeding system evolution and conservation genetics of the Tapisciaceae.

Portuguese wild grapevine genome re-sequencing (Vitis vinifera sylvestris)

The first genome of Vitis vinifera vinifera (PN40024), published in 2007, boosted grapevine related studies. While this reference genome is a suitable tool for the overall studies in the field, it lacks the ability to unveil changes accumulated during V. v. vinifera domestication. The subspecies V. v. sylvestris preserves wild characteristics, making it a good material to provide insights into V. v. vinifera domestication. The difference in the reproductive strategy between both subspecies is one of the characteristics that set them apart. While V. v. vinifera flowers are hermaphrodite, V. v. sylvestris is mostly dioecious. In this paper, we compare the re-sequencing of the genomes from a male and a female individual of the wild sylvestris, against the reference vinifera genome (PN40024). Variant analysis reveals a low number but with high impact modifications in coding regions, essentially non-synonymous single nucleotide polymorphisms and frame shifts caused by insertions and deletions. The sex-locus was manually inspected, and the results obtained are in line with the most recent works related with wild grapevine sex. In this paper we also describe for the first time RNA editing in transcripts of 14 genes in the sex-determining region, including VviYABBY and VviPLATZ.

The wild grape genome sequence provides insights into the transition from dioecy to hermaphroditism during grape domestication

Background

A key step in domestication of the grapevine was the transition from separate sexes (dioecy) in wild Vitis vinifera ssp. sylvestris (V. sylvestris) to hermaphroditism in cultivated Vitis vinifera ssp. sativa (V. vinifera). It is known that V. sylvestris has an XY system and V. vinifera a modified Y haplotype (Yh) and that the sex locus is small, but it has not previously been precisely characterized.

Results

We generate a high-quality de novo reference genome for V. sylvestris, onto which we map whole-genome re-sequencing data of a cross to locate the sex locus. Assembly of the full X, Y, and Yh haplotypes of V. sylvestris and V. vinifera sex locus and examining their gene content and expression profiles during flower development in wild and cultivated accessions show that truncation and deletion of tapetum and pollen development genes on the X haplotype likely causes male sterility, while the upregulation of a Y allele of a cytokinin regulator (APRT3) may cause female sterility. The downregulation of this cytokinin regulator in the Yh haplotype may be sufficient to trigger reversal to hermaphroditism. Molecular dating of X and Y haplotypes is consistent with the sex locus being as old as the Vitis genus, but the mechanism by which recombination was suppressed remains undetermined.

Conclusions

We describe the genomic and evolutionary characterization of the sex locus of cultivated and wild grapevine, providing a coherent model of sex determination in the latter and for transition from dioecy to hermaphroditism during domestication.

Genome Resequencing, Improvement of Variant Calling, and Population Genomic Analyses Provide Insights into the Seedlessness in the Genus Vitis

The seedlessness of grape derived from stenospermocarpy is one of the most prized traits of table or raisin grapes. It is controlled by a complex genetic system containing one dominant gene and multiple recessive genes. Here, we collected dense variation data from high-depth resequencing data of seeded, seedless, and wild relative grape genomes sequenced to > 37x mean depth. Variant calls were made using a modified variant calling pipeline that was suitable for highly diverse interspecific grape accessions. The modified pipeline enabled us to call several million more variants than the commonly recommended pipeline. The quality was validated by Sanger sequencing data and subsequently supported by the genetic population structure and the phylogenetic tree constructed using the obtained variation data, results of which were generally consistent with known pedigree and taxonomic classifications. Variation data enabled us to confirm a dominant gene and identify recessive loci for seedlessness. Incidentally, we found that grape cultivar Rizamat contains an ancestral chromosomal region of the dominant gene in Sultanina, a predominant seedlessness donor cultivar. Furthermore, we predicted new candidate causal genes including Vitvi01g00455, Vitvi08g01528, and Vitvi18g01237 associated with the recessive seedless-regulating loci, which showed high homology with genes that regulate seed development in Arabidopsis. This study provides fundamental insights relevant to variant calling from genome resequencing data of diverse interspecific hybrid germplasms such as grape and will accelerate future efforts aimed at crop improvement.

Default Sex and Single Gene Sex Determination in Dioecious Plants

A well-established hypothesis for the evolution of dioecy involves two genes linked at a sex-determining region (SDR). Recently there has been increased interest in possible single gene sex determination. Work in Populus has finally provided direct experimental evidence for single gene sex determination in plants using CRISPR-Cas9 to knock out a single gene and convert individuals from female to male. In poplar, the feminizing factor popARR17 acts as a "master regulator", analogous to the mammalian masculinizing factor SRY. The production of fully functional males from females by a simple single gene knockout is experimental evidence that an antagonistic male-determining factor does not exist in Populus. Mammals have a "default sex" (female), as do poplar trees (Populus), although the default sex in poplars is male. The occurrence of single gene sex determination with a default sex may be much commoner in plants than hitherto expected, especially when dioecy evolved via monoecy. The master regulator does not even need to be at the SDR (although it may be). In most poplars the feminizing factor popARR17 is not at the SDR, but instead a negative regulator of it. So far there is little information on how high-level regulators are connected to floral phenotype. A model is presented of how sex-determining genes could lead to different floral morphologies via MADS-box floral developmental genes.