Molecular and genetic characterization of novel S-RNases from a natural population of Nicotiana alata

Involvement of CgHSFB1 in the regulation of self-incompatibility in 'Shatian' pummelo

As self-incompatibility is a major issue in pummelo breeding and production, its mechanism in citrus was analyzed to improve breeding efficiency and reduce production costs. Rutaceae belongs to S-RNase type of gametophytic self-incompatibility. While the function of S-RNase/SLF and the mechanism of self-incompatibility have been studied extensively, the transcriptional regulation of S-RNase has been less studied. We performed transcriptome sequencing with the styles of 'Shatian' pummelo on the day of anthesis and 1-5 days before anthesis, and found that the transcript level of S-RNase gradually decreased with flower development. By analyzing differentially expressed genes and correlation with the expression trend of S-RNase, we identified a candidate gene, CgHSFB1, and utilized biochemical experiments such as yeast one-hybrid assay, electrophoretic mobility shift assay and dual-luciferase assay, as well as transient transformation of citrus calli and Citrus microcarpa and demonstrated that CgHSFB1 could directly bind to the S1-RNase promoter and repress the expression of S1-RNase, which is involved in the pummelo self-incompatibility response. In contrast, CgHSFB1 did not bind to the promoter of S2-RNase, and there was specificity in the regulation of S-RNase.

Mapping quantitative trait loci associated with self-(in)compatibility in goji berries (Lycium barbarum)

BACKGROUND

Goji (Lycium barbarum L.) is a perennial deciduous shrub widely distributed in arid and semiarid regions of Northwest China. It is highly valued for its medicinal and functional properties. Most goji varieties are naturally self-incompatible, posing challenges in breeding and cultivation. Self-incompatibility is a complex genetic trait, with ongoing debates regarding the number of self-incompatible loci. To date, no genetic mappings has been conducted for S loci or other loci related to self-incompatibility in goji.

RESULTS

We used genome resequencing to create a high-resolution map for detecting de novo single-nucleotide polymorphisms (SNP) in goji. We focused on 229 F1 individuals from self-compatible '13-19' and self-incompatible 'new 9' varieties. Subsequently, we conducted a quantitative trait locus (QTL) analysis on traits associated with self-compatibility in goji berries. The genetic map consisted of 249,327 SNPs distributed across 12 linkage groups (LGs), spanning a total distance of 1243.74 cM, with an average interval of 0.002 cM. Phenotypic data related to self-incompatibility, such as average fruit weight, fruit rate, compatibility index, and comparable compatibility index after self-pollination and geitonogamy, were collected for the years 2021-2022, as well as for an extra year representing the mean data from 2021 to 2022 (2021/22). A total of 43 significant QTL, corresponding to multiple traits were identified, accounting for more than 11% of the observed phenotypic variation. Notably, a specific QTL on chromosome 2 consistently appeared across different years, irrespective of the relationship between self-pollination and geitonogamy. Within the localization interval, 1180 genes were annotated, including Lba02g01102 (annotated as an S-RNase gene), which showed pistil-specific expression. Cloning of S-RNase genes revealed that the parents had two different S-RNase alleles, namely S1S11 and S2S8. S-genotype identification of the F1 population indicated segregation of the four S-alleles from the parents in the offspring, with the type of S-RNase gene significantly associated with self-compatibility.

CONCLUSIONS

In summary, our study provides valuable insights into the genetic mechanism underlying self-compatibility in goji berries. This highlights the importance of further positional cloning investigations and emphasizes the importance of integration of marker-assisted selection in goji breeding programs.

Establishment and application of Agrobacterium-delivered CRISPR/Cas9 system for wild tobacco (Nicotiana alata) genome editing

Clustered regularly interspaced short palindromic repeats (CRISPR)-associated protein 9 (CRISPR-Cas9) system has been widely applied in cultivated crops, but limited in their wild relatives. Nicotiana alata is a typical wild species of genus Nicotiana that is globally distributed as a horticultural plant and well-studied as a self-incompatibility model. It also has valuable genes for disease resistance and ornamental traits. However, it lacks an efficient genetic transformation and genome editing system, which hampers its gene function and breeding research. In this study, we developed an optimized hypocotyl-mediated transformation method for CRISPR-Cas9 delivery. The genetic transformation efficiency was significantly improved from approximately 1% to over 80%. We also applied the CRISPR-Cas9 system to target the phytoene desaturase (NalaPDS) gene in N. alata and obtained edited plants with PDS mutations with over 50% editing efficiency. To generate self-compatible N. alata lines, a polycistronic tRNA-gRNA (PTG) strategy was used to target exonic regions of allelic S-RNase genes and generate targeted knockouts simultaneously. We demonstrated that our system is feasible, stable, and high-efficiency for N. alata genome editing. Our study provides a powerful tool for basic research and genetic improvement of N. alata and an example for other wild tobacco species.

Gene coexpression analysis reveals key pathways and hub genes related to late-acting self-incompatibility in Camellia oleifera

INTRODUCTION

Self-incompatibility (SI) is an important strategy for plants to maintain abundant variation to enhance their adaptability to the environment. Camellia oleifera is one of the most important woody oil plants and is widely cultivated in China. Late acting self-incompatibility (LSI) in C. oleifera results in a relatively poor fruit yield in the natural state, and understanding of the LSI mechanism remains limited.

METHODS

To better understand the molecular expression and gene coexpression network in the LSI reaction in C. oleifera, we conducted self- and cross-pollination experiments at two different flower bud developmental stages (3-4 d before flowering and 1 d before flowering), and cytological observation, fruit setting rate (FSR) investigation and RNA-Seq analysis were performed to investigate the mechanism of the male -female interaction and identify hub genes responsible for the LSI in C. oleifera.

RESULTS

Based on the 21 ovary transcriptomes, a total of 7669 DEGs were identified after filtering out low-expression genes. Weighted gene coexpression network analysis (WGCNA) divided the DEGs into 15 modules. Genes in the blue module (1163 genes) were positively correlated with FSR, and genes in the pink module (339 genes) were negatively correlated with FSR. KEGG analysis indicated that flavonoid biosynthesis, plant MAPK signaling pathways, ubiquitin-mediated proteolysis, and plant-pathogen interaction were the crucial pathways for the LSI reaction. Fifty four transcription factors (TFs) were obtained in the two key modules, and WRKY and MYB were potentially involved in the LSI reaction in C. oleifera. Network establishment indicated that genes encoding G-type lectin S-receptor-like serine (lecRLK), isoflavone 3'-hydroxylase-like (CYP81Q32), cytochrome P450 87A3-like (CYP87A3), and probable calcium-binding protein (CML41) were the hub genes that positively responded to the LSI reaction. The other DEGs inside the two modules, including protein RALF-like 10 (RALF), F-box and pectin acetylesterase (MTERF5), might also play vital roles in the LSI reaction in C. oleifera.

DISCUSSION

Overall, our study provides a meaningful resource for gene network studies of the LSI reaction process and subsequent analyses of pollen-pistil interactions and TF roles in the LSI reaction, and it also provides new insights for exploring the mechanisms of the LSI response.

Allelic diversity of S-RNase alleles in diploid potato species

The S-ribonuclease sequences of 16 S-alleles derived from diploid types of Solanum are presented. A phylogenetic analysis and partial phenotypic analysis support the conclusion that these are functional S-alleles. S-Ribonucleases (S-RNases) control the pistil specificity of the self-incompatibility (SI) response in the genus Solanum and several other members of the Solanaceae. The nucleotide sequences of S-RNases corresponding to a large number of S-alleles or S-haplotypes have been characterised. However, surprisingly, few S-RNase sequences are available for potato species. The identification of new S-alleles in diploid potato species is desirable as these stocks are important sources of traits such as biotic and abiotic resistance. S-RNase sequences are reported here from three distinct diploid types of potato: cultivated Solanum tuberosum Group Phureja, S. tuberosum Group Stenotomum, and the wild species Solanum okadae. Partial S-RNase sequences were obtained from pistil RNA by RT-PCR or 3'RACE (Rapid Amplification of cDNA Ends) using a degenerate primer. Full-length sequences were obtained for two alleles by 5'RACE. Database searches with these sequences identified 16 S-RNases in total, all of which are novel. The sequence analysis revealed all the expected features of functional S-RNases. Phylogenetic analysis with selected published S-RNase and S-like-RNase sequences from the Solanaceae revealed extensive trans-generic evolution of the S-RNases and a clear distinction from S-like-RNases. Pollination tests were used to confirm the self-incompatibility status and cross-compatibility relationships of the S. okadae accessions. All the S. okadae accessions were found to be self-incompatible as expected with crosses amongst them exhibiting both cross-compatibility and semi-compatibility consistent with the S-genotypes determined from the S-RNase sequence data. The progeny analysis of four semi-compatible crosses examined by allele-specific PCR provided further confirmation that these are functional S-RNases.

NnSR1, a class III non-S-RNase specifically induced in Nicotiana alata under phosphate deficiency, is localized in endoplasmic reticulum compartments

A combined strategy of phosphate (Pi) remobilization from internal and external RNA sources seems to be conserved in plants exposed to Pi starvation. Thus far, the only ribonucleases (RNases) reported to be induced in Nicotiana alata undergoing Pi deprivation are extracellular S-like RNase NE and NnSR1. NnSR1 is a class III non S-RNase of unknown subcellular location. Here, we examine the hypothesis that NnSR1 is an intracellular RNase derived from the self-incompatibility system with specific expression in self-incompatible Nicotiana alata. NnSR1 was not induced in self-compatible Nicotiana species exposed to Pi deprivation. NnSR1 conjugated with a fluorescent protein and transiently expressed in Arabidopsis protoplasts and Nicotiana leaves showed that the fusion protein co-localized with an endoplasmic reticulum (ER) marker. Subcellular fractionation by ultracentrifugation of roots exposed to Pi deprivation revealed that the native NnSR1 migrated in parallel with the BiP protein, a typical ER marker. To our knowledge, NnSR1 is the first class III RNase reported to be localized in ER compartments. The induction of NnSR1 was detected earlier than the extracellular RNase NE, suggesting that intracellular RNA may be the first source of Pi used by the cell under Pi stress.

In vitro inhibition of incompatible pollen tubes in Nicotiana alata involves the uncoupling of the F-actin cytoskeleton and the endomembrane trafficking system

In the S-RNase-based self-incompatibility system, subcellular events occurring in the apical region of incompatible pollen tubes during the pollen rejection process are poorly understood. F-actin dynamics and endomembrane trafficking are crucial for polar growth, which is temporally and spatially controlled in the tip region of pollen tubes. Thus, we developed a simple in vitro assay to study the changes in the F-actin cytoskeleton and the endomembrane system at the apical region of incompatible pollen tubes in Nicotiana alata. Growth but not germination of pollen tubes of S c₁₀-, S₇₀-, and S₇₅-haplotypes was selectively inhibited by style extracts carrying the same haplotypes. Pollen F-actin cytoskeleton and endomembrane system, visualized by fluorescent markers, were normal during the initial 60 min of pollen culture in the presence of compatible and incompatible style extracts. Additional culture resulted in complete growth arrest and critical alterations in the integrity of the F-actin cytoskeleton and the endomembrane system of incompatible pollen tubes. The F-actin ring and the V-shaped zone disappeared from the apical region, while distorted F-actin cables and progressive formation of membrane aggregates evolved in the subapical region and the shank. The vacuolar network of incompatible pollen tubes invaded the tip region, but vacuolar membrane integrity remained mostly unaffected. The polar growth machinery of incompatible pollen tubes was uncoupled, as evidenced by the severe disruption of colocalization between the F-actin cytoskeleton and the endomembrane compartments. A model of pollen rejection integrating the main subcellular events occurring in incompatible pollen is discussed.

NnSR1, a class III non-S-RNase constitutively expressed in styles, is induced in roots and stems under phosphate deficiency in Nicotiana alata

BACKGROUND AND AIMS

Non-S-ribonucleases (non-S-RNases) are class III T2 RNases constitutively expressed in styles of species with S-RNase-based self-incompatibility. So far, no function has been attributed to these RNases. The aim of this work is to examine if NnSR1, a non-S-RNase from Nicotiana alata, is induced under conditions of phosphate (Pi) deprivation. The hypothesis is that under Pi-limited conditions, non-S-RNase functions may resemble the role of S-like RNases. To date, the only RNases reported to be induced by Pi deficiency are class I and class II S-like RNases, which are phylogenetically different from the class III clade of RNases.

METHODS

Gene and protein expression of NnSR1 were assayed in plants grown hydroponically with and without Pi, by combining RT-PCR, immunoblot and enzymatic activity approaches.

KEY RESULTS

NnSR1 transcripts were detected in roots 7 d after Pi deprivation and remained stable for several days. Transcript expression was correlated based on Pi availability in the culture medium. Antiserum against a peptide based on a hypervariable domain of NnSR1 recognized NnSR1 in roots and stems but not leaves exposed to Pi shortage. NnSR1 was not detected in culture medium and was pelleted with the microsomal fraction, suggesting that it was membrane-associated or included in large compartments. The anti-NnSR1 inhibited selectively the enzymatic activity of a 31-kDa RNase indicating that NnSR1 was induced in an enzymatically active form.

CONCLUSIONS

The induction of NnSR1 indicates that there is a general recruitment of all classes of T2 RNases in response to Pi shortage. NnSR1 appears to have regained ancestral functions of class III RNases related to strategies to cope with Pi limitation and also possibly with other environmental challenges. This constitutes the first report for a specific function of class III RNases other than S-RNases.

A new dual-specific incompatibility allele revealed by absence of glycosylation in the conserved C2 site of a Solanum chacoense S-RNase

The stylar determinant of gametophytic self-incompatibility (GSI) in Solanaceae, Rosaceae, and Plantaginaceae is an S-RNase encoded by a multiallelic S-locus. The primary structure of S-RNases shows five conserved (C) and two hypervariable (HV) regions, the latter forming a domain implicated in S-haplotype-specific recognition of the pollen determinant to SI. All S-RNases are glycosylated at a conserved site in the C2 region, although previous studies have shown that N-linked glycans at this position are not required for S-haplotype-specific recognition and pollen rejection. Here the incompatibility phenotype of three constructs derived from an originally monoglycosylated S11-RNase of Solanum chacoense, that were designed to explore the role of the HV domain in determining pollen recognition and the role of the N-linked glycan in the C2 region, is reported. In one series of experiments, a second glycosylation site was introduced in the HVa region to test for inhibition of pollen-specific recognition. This modification does not impede pollen rejection, although analysis shows incomplete glycosylation at the new site in the HVa region. A second construct, designed to permit complete glycosylation at the HVa site by suppression of the conserved site in the C2 region, did increase the degree of site occupancy, but, again, glycosylation was incomplete. Plants expressing this construct rejected S 11 pollen and, surprisingly, also rejected S 13 pollen, thus displaying an unusual dual specificity phenotype. This construct differs from the first by the absence of the conserved C2 glycosylation site, and thus the dual specificity is observed only in the absence of the C2 glycan. A third construct, completely lacking glycosylation sites, conferred an ability to reject only S 11 pollen, disproving the hypothesis that lack of a conserved glycan would confer a universal pollen rejection phenotype to the plant.

Disorganization of F-actin cytoskeleton precedes vacuolar disruption in pollen tubes during the in vivo self-incompatibility response in Nicotiana alata

BACKGROUND AND AIMS

The integrity of actin filaments (F-actin) is essential for pollen-tube growth. In S-RNase-based self-incompatibility (SI), incompatible pollen tubes are inhibited in the style. Consequently, research efforts have focused on the alterations of pollen F-actin cytoskeleton during the SI response. However, so far, these studies were carried out in in vitro-grown pollen tubes. This study aimed to assess the timing of in vivo changes of pollen F-actin cytoskeleton taking place after compatible and incompatible pollinations in Nicotiana alata. To our knowledge, this is the first report of the in vivo F-actin alterations occurring during pollen rejection in the S-RNase-based SI system.

METHODS

The F-actin cytoskeleton and the vacuolar endomembrane system were fluorescently labelled in compatibly and incompatibly pollinated pistils at different times after pollination. The alterations induced by the SI reaction in pollen tubes were visualized by confocal laser scanning microscopy.

KEY RESULTS

Early after pollination, about 70 % of both compatible and incompatible pollen tubes showed an organized pattern of F-actin cables along the main axis of the cell. While in compatible pollinations this percentage was unchanged until pollen tubes reached the ovary, pollen tubes of incompatible pollinations underwent gradual and progressive F-actin disorganization. Colocalization of the F-actin cytoskeleton and the vacuolar endomembrane system, where S-RNases are compartmentalized, revealed that by day 6 after incompatible pollination, when the pollen-tube growth was already arrested, about 80 % of pollen tubes showed disrupted F-actin but a similar percentage had intact vacuolar compartments.

CONCLUSIONS

The results indicate that during the SI response in Nicotiana, disruption of the F-actin cytoskeleton precedes vacuolar membrane breakdown. Thus, incompatible pollen tubes undergo a sequential disorganization process of major subcellular structures. Results also suggest that the large pool of S-RNases released from vacuoles acts late in pollen rejection, after significant subcellular changes in incompatible pollen tubes.

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

The transition from self-incompatibility to self-compatibility is a common transition in angiosperms often reported in populations at the edge of species range limits. Geographically distinct populations of wild tomato species (Solanum section Lycopersicon (Solanaceae)) have been described as polymorphic for mating system with both self-incompatible and self-compatible populations. Using controlled pollinations and sequencing of the S-RNase mating system gene, we test the compatibility status of a population of S. peruvianum located near its southern range limit. Pollinations among plants of known genotypes revealed strong self-incompatibility; fruit set following compatible pollinations was significantly higher than following incompatible pollinations for all tested individuals. Sequencing of the S-RNase gene in parents and progeny arrays was also as predicted under self-incompatibility. Molecular variation at the S-RNase locus revealed a diverse set of alleles, and heterozygosity in over 500 genotyped individuals. We used controlled crosses to test the specificity of sequences recovered in this study; in all cases, results were consistent with a unique allelic specificity for each tested sequence, including two alleles sharing 92% amino-acid similarity. Site-specific patterns of selection at the S-RNase gene indicate positive selection in regions of the gene associated with allelic specificity determination and purifying selection in previously characterized conserved regions. Further, there is broad convergence between the present and previous studies in specific amino-acid positions inferred to be evolving under positive selection.