Functional gametophytic self-incompatibility in a peripheral population of Solanum peruvianum (Solanaceae)

S-RNase Alleles Associated With Self-Compatibility in the Tomato Clade: Structure, Origins, and Expression Plasticity

The self-incompatibility (SI) system in the Solanaceae is comprised of cytotoxic pistil S-RNases which are countered by S-locus F-box (SLF) resistance factors found in pollen. Under this barrier-resistance architecture, mating system transitions from SI to self-compatibility (SC) typically result from loss-of-function mutations in genes encoding pistil SI factors such as S-RNase. However, the nature of these mutations is often not well characterized. Here we use a combination of S-RNase sequence analysis, transcript profiling, protein expression and reproductive phenotyping to better understand different mechanisms that result in loss of S-RNase function. Our analysis focuses on 12 S-RNase alleles identified in SC species and populations across the tomato clade. In six cases, the reason for gene dysfunction due to mutations is evident. The six other alleles potentially encode functional S-RNase proteins but are typically transcriptionally silenced. We identified three S-RNase alleles which are transcriptionally silenced under some conditions but actively expressed in others. In one case, expression of the S-RNase is associated with SI. In another case, S-RNase expression does not lead to SI, but instead confers a reproductive barrier against pollen tubes from other tomato species. In the third case, expression of S-RNase does not affect self, interspecific or inter-population reproductive barriers. Our results indicate that S-RNase expression is more dynamic than previously thought, and that changes in expression can impact different reproductive barriers within or between natural populations.

Spread of self-compatibility constrained by an intrapopulation crossing barrier

Mating system transitions from self-incompatibility (SI) to self-compatibility (SC) are common in plants. In the absence of high levels of inbreeding depression, SC alleles are predicted to spread due to transmission advantage and reproductive assurance. We characterized mating system and pistil-expressed SI factors in 20 populations of the wild tomato species Solanum habrochaites from the southern half of the species range. We found that a single SI to SC transition is fixed in populations south of the Rio Chillon valley in central Peru. In these populations, SC correlated with the presence of the hab-6 S-haplotype that encodes a low activity S-RNase protein. We identified a single population segregating for SI/SC and hab-6. Intrapopulation crosses showed that hab-6 typically acts in the expected codominant fashion to confer SC. However, we found one specific S-haplotype (hab-10) that consistently rejects pollen of the hab-6 haplotype, and results in SI hab-6/hab-10 heterozygotes. We suggest that the hab-10 haplotype could act as a genetic mechanism to stabilize mixed mating in this population by presenting a disadvantage for the hab-6 haplotype. This barrier may represent a mechanism allowing for the persistence of SI when an SC haplotype appears in or invades a population.

Wolfberry genomes and the evolution of Lycium (Solanaceae)

Wolfberry Lycium, an economically important genus of the Solanaceae family, contains approximately 80 species and shows a fragmented distribution pattern among the Northern and Southern Hemispheres. Although several herbaceous species of Solanaceae have been subjected to genome sequencing, thus far, no genome sequences of woody representatives have been available. Here, we sequenced the genomes of 13 perennial woody species of Lycium, with a focus on Lycium barbarum. Integration with other genomes provides clear evidence supporting a whole-genome triplication (WGT) event shared by all hitherto sequenced solanaceous plants, which occurred shortly after the divergence of Solanaceae and Convolvulaceae. We identified new gene families and gene family expansions and contractions that first appeared in Solanaceae. Based on the identification of self-incompatibility related-gene families, we inferred that hybridization hotspots are enriched for genes that might be functioning in gametophytic self-incompatibility pathways in wolfberry. Extremely low expression of LOCULE NUBER (LC) and COLORLESS NON-RIPENING (CNR) orthologous genes during Lycium fruit development and ripening processes suggests functional diversification of these two genes between Lycium and tomato. The existence of additional flowering locus C-like MADS-box genes might correlate with the perennial flowering cycle of Lycium. Differential gene expression involved in the lignin biosynthetic pathway between Lycium and tomato likely illustrates woody and herbaceous differentiation. We also provide evidence that Lycium migrated from Africa into Asia, and subsequently from Asia into North America. Our results provide functional insights into Solanaceae origins, evolution and diversification.

Population Genomics of the "Arcanum" Species Group in Wild Tomatoes: Evidence for Separate Origins of Two Self-Compatible Lineages

Given their diverse mating systems and recent divergence, wild tomatoes (Solanum section Lycopersicon) have become an attractive model system to study ecological divergence, the build-up of reproductive barriers, and the causes and consequences of the breakdown of self-incompatibility. Here we report on a lesser-studied group of species known as the "Arcanum" group, comprising the nominal species Solanum arcanum, Solanum chmielewskii, and Solanum neorickii. The latter two taxa are self-compatible but are thought to self-fertilize at different rates, given their distinct manifestations of the morphological "selfing syndrome." Based on experimental crossings and transcriptome sequencing of a total of 39 different genotypes from as many accessions representing each species' geographic range, we provide compelling evidence for deep genealogical divisions within S. arcanum; only the self-incompatible lineage known as "var. marañón" has close genealogical ties to the two self-compatible species. Moreover, there is evidence under multiple inference schemes for different geographic subsets of S. arcanum var. marañón being closest to S. chmielewskii and S. neorickii, respectively. To broadly characterize the population-genomic consequences of these recent mating-system transitions and their associated speciation events, we fit demographic models indicating strong reductions in effective population size, congruent with reduced nucleotide and S-locus diversity in the two independently derived self-compatible species.

The expression pattern of the Pho1a genes encoding plastidic starch phosphorylase correlates with the degradation of starch during fruit ripening in green-fruited and red-fruited tomato species

Genes encoding plastidic starch phosphorylase Pho1a were identified in 10 tomato species (Solanum section Lycopersicon). Pho1a genes showed higher variability in green-fruited than in red-fruited tomato species, but had an extremely low polymorphism level compared with other carbohydrate metabolism genes and an unusually low ratio of intron to exon single nucleotide polymorphisms (SNPs). In red-fruited species, Pho1a was expressed in all analysed tissues, including fruit at different developmental stages, with the highest level in mature green fruit, which is strong sink organ importing sucrose and accumulating starch. In green-fruited species Solanum peruvianum and Solanum arcanum, the Pho1a expression level was similar in mature green and ripe fruit, whereas in Solanum chmielewskii, it was higher in ripe fruit, and in Solanum habrochaites, the dynamics of fruit-specific Pho1a expression was similar to that in red-fruited tomatoes. During fruit development, in red-fruited Solanum lycopersicum, sucrose level was low, the monosaccharide content increased; in green-fruited S. peruvianum, the sucrose concentration increased and those of monosaccharides decreased. In both species, the starch content and Pho1a expression were downregulated. The evolutionary topology based on Pho1a sequences was consistent with the current division of tomatoes into red-fruited and green-fruited species, except for S. habrochaites.

Colonization, Baker's law, and the evolution of gynodioecy in Hawaii: implications from a study of Lycium carolinianum

PREMISE

As Baker's law suggests, the successful colonization of oceanic islands is often associated with uniparental reproduction (self-fertility), but the high incidence of dimorphism (dioecy, gynodioecy) on islands complicates this idea. Lycium carolinianum is widespread, occurring on the North American mainland and the Hawaiian Islands. We examined Baker's ideas for mainland and island populations of L. carolinianum and examined inbreeding depression as a possible contributor to the evolution of gynodioecy on Maui.

METHODS

Controlled crosses were conducted in two mainland populations and two populations in Hawaii. Treatments included self and cross pollination, unmanipulated controls, and autogamy/agamospermy. Alleles from the self-incompatibility S-RNase gene were isolated and compared between mainland and island populations. Given self-compatibility in Hawaii, we germinated seeds from self- and cross- treatments and estimated inbreeding depression using seven traits and a measure of cumulative fitness.

RESULTS

Mainland populations of Lycium carolinianum are predominately self-incompatible with some polymorphism for self-fertility, whereas Hawaiian populations are self-compatible. Concordantly, S-RNase allelic diversity is reduced in Hawaii compared to the mainland. Hawaiian populations also exhibit significant inbreeding depression.

CONCLUSIONS

Self-compatibility in Hawaii and individual variation in self-fertility in mainland populations suggests that a colonization filter promoting uniparental reproduction may be acting in this system. Comparison of S-RNase variation suggests a collapse of allelic diversity and heterozygosity at the S-RNase locus in Hawaii, which likely contributed to mate limitation upon arrival to the Pacific. Inbreeding depression coupled with autonomous self-fertilization may have led to the evolution of gynodioecy on Maui.

Overcoming Self-Incompatibility in Diploid Potato Using CRISPR-Cas9

Potato breeding can be redirected to a diploid inbred/F1 hybrid variety breeding strategy if self-compatibility can be introduced into diploid germplasm. However, the majority of diploid potato clones (Solanum spp.) possess gametophytic self-incompatibility that is primarily controlled by a single multiallelic locus called the S-locus which is composed of tightly linked genes, S-RNase (S-locus RNase) and multiple SLFs (S-locus F-box proteins), which are expressed in the style and pollen, respectively. Using S-RNase genes known to function in the Solanaceae gametophytic SI mechanism, we identified S-RNase alleles with flower-specific expression in two diploid self-incompatible potato lines using genome resequencing data. Consistent with the location of the S-locus in potato, we genetically mapped the S-RNase gene using a segregating population to a region of low recombination within the pericentromere of chromosome 1. To generate self-compatible diploid potato lines, a dual single-guide RNA (sgRNA) strategy was used to target conserved exonic regions of the S-RNase gene and generate targeted knockouts (KOs) using a Clustered Regularly Interspaced Short Palindromic Repeats/CRISPR-associated protein 9 (Cas9) approach. Self-compatibility was achieved in nine S-RNase KO T₀ lines which contained bi-allelic and homozygous deletions/insertions in both genotypes, transmitting self compatibility to T₁ progeny. This study demonstrates an efficient approach to achieve stable, consistent self-compatibility through S-RNase KO for use in diploid potato breeding approaches.

TAI vacuolar invertase orthologs: the interspecific variability in tomato plants (Solanum section Lycopersicon)

Understanding the genetic mechanisms underlying carbohydrate metabolism can promote the development of biotechnological advances in fruit plants. The flesh tomato fruit represents an ideal system for examining the role of sucrose cleavage enzymes in fruit development, and wild tomato species differing in storage sugars serve as an excellent research material for this purpose. Plant vacuolar invertase is a key enzyme of sucrose metabolism in the sink organs. In the present study, we identified complete gene sequences encoding the TAI vacuolar invertase in 11 wild and one cultivated tomato accessions of the Solanum section Lycopersicon. The average level of interspecific polymorphism in TAI genes was 8.58%; however, in the green-fruited tomatoes, the TAI genes contained 100 times more SNPs than those in the red-fruited accessions. The TAI proteins demonstrated 8% variability, whereas the red-fruited species had none. A TAI-based phylogenetic tree revealed two main clusters containing self-compatible and self-incompatible species, which concurs with the previous crossability-based division and demonstrates that the TAI genes reflect the evolutionary relationships between the red- and green-fruited tomatoes. Furthermore, we detected differential expression patterns of the TAI genes in the fruits of wild and cultivated tomatoes, which corresponded to sugar composition. The polymorphism analysis of the TAI acid invertases of Solanum section Lycopersicon species will contribute to the understanding of the genetic potential of TAI genes to impact tomato breeding through genetic engineering of the carbohydrate composition in the fruit.

A Comprehensive Study of Molecular Evolution at the Self-Incompatibility Locus of Rosaceae

The family Rosaceae includes a range of important fruit trees, most of which have the S-RNase-based self-incompatibility (SI). Several models have been developed to explain how pollen (SLF) and pistil (S-RNase) components of the S-locus interact. It was discovered in 2010 that additional SLF proteins are involved in pollen specificity, and a Collaborative Non-Self Recognition model has been proposed for SI in Solanaceae; however, the validity of such model remains to be elucidated for other species. The results of this study support the divergent evolution of the S-locus genes from two Rosaceae subfamilies, Prunoideae/Amygdaloideae and Maloideae, The difference identified in the selective pressures between the two lineages provides evidence for positive selection at specific sites in both the S-RNase and the SLF proteins. The evolutionary findings of this study support the role of multiple SLF proteins leading to a Collaborative Non-Self Recognition model for SI in the Maloideae. Furthermore, the identification of the sites responsible for SI specificity determination and the mapping of these sites onto the modelled tertiary structure of ancestor proteins provide useful information for rational functional redesign and protein engineering for the future engineering of new functional alleles providing increased diversity in the SI system in the Maloideae.

Diversification and distinctive structural features of S-RNase alleles in the genus Solanum

The multigenic and multiallelic S-locus in plants is responsible for the gametophytic self-incompatibility system, which is important to prevent the detrimental effects of self-fertilization and inbreeding depression. Several studies have discussed the importance of punctual mutations, recombination, and natural selection in the generation of allelic diversity in the S-locus. However, there has been no wide-ranging study correlating the molecular evolution and structural aspects of the corresponding proteins in Solanum. Therefore, we evaluated the molecular evolution of one gene in this locus and generated a statistically well-supported phylogenetic tree, as well as evidence of positive selection, helping us to understand the diversification of S alleles in Solanum. The three-dimensional structures of some of the proteins corresponding to the major clusters of the phylogenetic tree were constructed and subsequently submitted to molecular dynamics to stabilize the folding and obtain the native structure. The positively selected amino acid residues were predominantly located in the hyper variable regions and on the surface of the protein, which appears to be fundamental for allele specificity. One of the positively selected residues was identified adjacent to a conserved strand that is crucial for enzymatic catalysis. Additionally, we have shown significant differences in the electrostatic potential among the predicted molecular surfaces in S-RNases. The structural results indicate that local changes in the three-dimensional structure are present in some regions of the molecule, although the general structure seems to be conserved. No previous study has described such structural variations in S-RNases.

The timing of molecular and morphological changes underlying reproductive transitions in wild tomatoes (Solanum sect. Lycopersicon)

Molecular mechanisms underlying the transition from genetic self-incompatibility to self-compatibility are well documented, but the evolution of other reproductive trait changes that accompany shifts in reproductive strategy (mating system) remains comparatively under-investigated. A notable exception is the transition from exserted styles to styles with recessed positions relative to the anthers in wild tomatoes (Solanum Section Lycopersicon). This phenotypic change has been previously attributed to a specific mutation in the promoter of a gene that influences style length (style2.1); however, whether this specific regulatory mutation arose concurrently with the transition from long to short styles, and whether it is causally responsible for this phenotypic transition, has been poorly investigated across this group. To address this gap, we assessed 74 accessions (populations) from 13 species for quantitative genetic variation in floral and reproductive traits as well as the presence/absence of deletions at two different locations (StyleD1 and StyleD2) within the regulatory region upstream of style2.1. We confirmed that the putatively causal deletion variant (a 450-bp deletion at StyleD1) arose within self-compatible lineages. However, the variation and history of both StyleD1 and StyleD2 was more complex than previously inferred. In particular, although StyleD1 was statistically associated with differences in style length and stigma exsertion across all species, we found no evidence for this association within two species polymorphic for the StyleD1 mutation. We conclude that the previous association detected between phenotypic and molecular differences is most likely due to a phylogenetic association rather than a causal mechanistic relationship. Phenotypic variation in style length must therefore be due to other unexamined linked variants in the style2.1 regulatory region.

Testing for "snowballing" hybrid incompatibilities in Solanum: impact of ancestral polymorphism and divergence estimates

Two recent high-profile studies offered empirical evidence for a "snowballing" accumulation of postzygotic incompatibilities in Drosophila and Solanum (tomatoes). Here we present a reanalysis of the Solanum data that is motivated by population genetic principles. Specifically, the high levels of intraspecific nucleotide polymorphism in wild tomato species and presumably large effective population size throughout the divergence history of this clade imply that ancestral polymorphism should be taken into account when evaluating sequence divergence between species. Based on our reanalyses of synonymous-site divergence between the four focal Solanum species and a wide range of ancestral polymorphism, we assessed under which conditions the reported accumulation of seed sterility factors supports the snowball effect. Our results highlight the pivotal impact of levels of ancestral polymorphism and alternate divergence values, and they illustrate that robust tests of the snowball effect in Solanum require genome-wide estimates of divergence.