Development of Self-Compatible B. rapa by RNAi-Mediated S Locus Gene Silencing

Abstract

The self-incompatibility (SI) system is genetically controlled by a single polymorphic locus known as the S-locus in the Brassicaceae. Pollen rejection occurs when the stigma and pollen share the same S-haplotype. Recognition of S-haplotype specificity has recently been shown to involve at least two S-locus genes, S-receptor kinase (SRK) and S-locus protein 11 or S locus Cysteine-rich (SP11/SCR) protein. Here, we examined the function of S60, one SP11/SCR allele of B. rapa cv. Osome, using a RNAi-mediated gene silencing approach. The transgenic RNAi lines were highly self-compatible, and this trait was stable in subsequent generations, even after crossing with other commercial lines. These findings also suggested that the resultant self-compatibility could be transferred to commercial cultivars with the desired performances in B. rapa.

About the Authors:

Hee-Jeong Jung

Affiliation: Department of Horticulture, Sunchon National University, Suncheon, Republic of Korea

Hyo-Jin Jung

Affiliation: Department of Horticulture, Jeollanam-do Agricultural Research and Extension Services, Najusi, Republic of Korea

Nasar Uddin Ahmed

Affiliation: Department of Horticulture, Sunchon National University, Suncheon, Republic of Korea

Jong-In Park

Affiliation: Department of Horticulture, Sunchon National University, Suncheon, Republic of Korea

Kwon-Kyoo Kang

Affiliation: Department of Horticulture, Hankyong National University, Ansung, Republic of Korea

Yoonkang Hur

Affiliation: Department of Biology, Chungnam National University, Daejeon, Republic of Korea

Yong-Pyo Lim

Affiliation: Department of Horticulture, Chungnam National University, Daejeon, Republic of Korea

Ill-Sup Nou

* E-mail: [email protected]

Affiliation: Department of Horticulture, Sunchon National University, Suncheon, Republic of Korea

Introduction

In flowering plants, self-incompatibility (SI) is a genetic system that promotes outcrossing by rejecting self-related pollen [1]–[3]. In Brassica, recognition of pollen is controlled by S haplotypes (designated S1, S2,... Sn), each of which consists of the pollen determinant gene, the S locus protein 11 (SP11, also called SCR) [4]–[6], and the pistil determinant gene, the S receptor kinase (SRK) [7]–[9]. Genes of SP11/SCR (SP11 hereafter) and SRK are closely linked each other at the S locus [10], and the S locus contain multiple alleles [11]–[14]. The self-incompatibility response occurs when the pollen and pistil share the same S haplotype. SP11 interacts with SRK of the same S haplotype and activates its kinase domain [15]–[17]. This activation is believed to elicit a signaling cascade within the stigmatic papilla cell that leads to rejection of self-pollen [18], [19]. SP11 is expressed in the anther tapetum, a sporophytic tissue. Therefore, the SI phenotype in pollen is determined by the dominance relationships between the two S haplotypes carried by the plant. Based on these relationships, the S haplotypes in Brassica have been broadly classified into two groups: pollen-dominant S haplotypes (class I) and pollen-recessive S haplotypes (class II) [14], [20]–[22]. Pollen-dominant S haplotypes (such as S8, S9, S12 and S52 in Brassica rapa) are generally co-dominant, and are almost always dominant over recessive S haplotypes (such as S44, S60, S40 and S29) [14], [21]. In S heterozygotes with dominant and recessive S haplotypes, the expression of the recessive SP11 allele is silenced as a result of tapetum-specific de novo cytosine methylation in its promoter region immediately before the initiation of SP11 transcription [23], [24]. Additionally, several types of cysteine-rich peptides/polypeptides (CRPs) are expressed specifically in flowers and seeds, where they play reproductive regulatory roles [25]. For instance, SP11, with eight conserved cysteines belonging to a subclass of defensin-like proteins, is involved in the inhibition of self-pollen germination and pollen tube growth [15]. Defensin-like LUREs act as attractants, guiding pollen tubes to the embryo sac [26], [27].

SI is one of the most important obstructions to B. rapa seed production and for that self-compatible (SC) cultivar is very important in case of commercial cultivation. Manipulation of S locus genes is one of the most recognized way so far to convert SI into SC in Brassica. RNA interference (RNAi)-mediated suppression of stylar 120 kDa glycoprotein (120K), a non-S-factor, results in breakdown of the capability of the pistil to reject self-pollen, suggesting that it is required for gametophytic self-incompatibility (GSI) function [28]. In plants, RNAi can be used to regulate endogenous genes [29], and by utilizing a partial sequence of an endogenous gene in the inverted repeat regions of the silencing construct, high-level silencing of the target gene expression can also be achieved. The various RNAi techniques each have advantages and disadvantages with respect to how persistent their effect is and the range of plants to which they can be applied. For example, bombardment can be applied to any plant, but produces only transient effects. Alternatively, transformation with ihpRNA-expressing transgenes provides stable and heritable gene silencing. ihpRNA transgenes have been shown to be very effective for a wide range of target genes in various plant species [30]–[35], indicating that the RNAi mechanism is probably conserved in all plant species.

Here, we report for the first time the use of RNAi gene-silencing constructs to achieve silencing of the S locus SP11 gene S60 resulting in the development of a self-compatible B. rapa transgenic line. Furthermore, we demonstrate the stable inheritance of these phenotypes in progeny derived by either selfing or intercrossing and assess the performance of these lines.

Materials and Methods

Plant materials and growth conditions

Brassica rapa cv. Osome plants were grown under sterile conditions on MS medium during transformation and tissue culture. Transgenic B. rapa plants and crossing generations were grown in the greenhouse of the Department of Horticulture, Sunchon National University, Korea, under natural light conditions.

Preparation of RNAi constructs

The coding sequence of S60 containing 285 bp nucleotides was placed upstream and downstream of the Gus gene encoding the β-glucuronidase fragment in opposite directions using the methods described by Chuang and Meyerowitz [36]. This Gus gene with the sense and antisense coding sequence of S60 was placed in the SalI and SacI site between the SSH and nopaline synthase terminator of binary vector pBI101 constructed from pBI121 [37]. The SSH fragment was then removed and the S60 promoter of 489 bp was placed in the KpnI and SalI site. The resulting RNAi construct was denoted as S60-SP11RNAi (Figure 1). This construct was then introduced into Agrobacterium tumefaciens strain EHA105.

[Figure omitted. See PDF.]

Figure 1. Schematic representations of the S60-SP11RNAi vector construct.

A 285 bp cDNA fragment of S60-SP11 was put into the upstream and downstream of the 1023 bp Gus linker of binary vector pBI101 in opposite orientations under control of the S60-SP11 promoter. LB, left border of T-DNA; NPTII, gene for neomycin phosphotransferase for kanamycin resistance; 35S-Pro, 35S promoter of cauliflower mosaic virus; Nos-Pro, promoter of nos (nopaline synthase) gene; Nos-Ter, terminator of nos (nopaline synthase) gene; S60-SP11 pro, S60-SP11 promoter; HPT, the Hygromycin-resistance gene; RB, right border of T-DNA.

https://doi.org/10.1371/journal.pone.0049497.g001

Plant transformation and regeneration

The hypocotyl transformation protocol developed for B. rapa [38], [39] was followed in this study, with some modifications. Briefly, seeds of B. rapa cv. Osome were surface-sterilized by washing with 70% ethanol for 2 min, 1% sodium hypochlorite for 15 min and double distilled water for 3–4 times. Seeds were germinated and grown in 0.1× MS medium in a culture room maintained at 22–24°C with a 16 h light/8 h dark photoperiod at a light intensity of 4500–5500 lux. Hypocotyls were excised from 6 to 7-day-old seedlings, cut into segments 2–4 mm in length, and placed onto MS-1 medium and pre-cultured for 24 h under indirect continuous light. Explants were then immersed in a suspension of 1xl08 bacteria/ml for 30 min with shaking at 40 rpm, then returned to feeder plates. After two days of cocultivation with Agrobacterium, explants were transferred to B5–1 medium supplemented with 500 mg/l carbenicillin and kept at 24°C under continuous light at 7500 lux intensity for 3–7 days, then transferred to B5-BZ shoot regeneration medium. These explants were cultured for seven days in B5–BZ medium supplemented with 500 mg/l carbenicillin, 10 mg/l Hyg, and 10 mg/l AgNO3, followed by 14 days in medium supplemented with 500 mg/l carbenicillin and 20 mg/l Hyg, 14 days in medium supplemented with 500 mg/l carbenicillin and 30 mg/l Hyg and finally, 14 days in medium supplemented with 500 mg/l carbenicillin and 30 mg/l Hyg. The cultures were then transferred onto B5–0 shoot maturation medium supplemented with 500 mg/1 carbenicillin and 50 mg/l Hyg. Two weeks later, shoots were trimmed to contain 2–3 nodes and then placed on B5 root induction medium supplemented with 2 mg/l IBA, 500 mg/l carbenicillin and 50 mg/l Hyg. Roots developed on some of the shoots after two weeks. Shoots that had not rooted were re-cut at the base and placed back onto the medium for another 2–4 weeks.

DNA extraction and PCR analysis

The putative T0, T1 and BC5F2 plants were analyzed by PCR to confirm the presence of transgenes. To accomplish this, plant genomic DNA was isolated from the leaves of each line using the DNeasy Plant Mini Kit (Qiagen, USA). The primers specific for the S60-SP11RNAi cassette (5′- GGC ATA TGA AGC TTG TCG ACA TGA TTT AAC TTT GCA ACAG -3′; 5′- CTG CAG GAG CTC GCG GCC GCA TGA TTT AAC TTT GCA ACAG -3′) and for NPTII (5′- CAA GAT GGA TTG CAC GCA GG -3′; 5′- GAA GAA CTC GTC AAG AAG GCG -3′) were used to identify transgenic plants. PCR reactions were carried out in a 20 l mixture at 94°C for 5 min, then subjected to 35 cycles of amplification at 94°C for 1 min, 55°C for 1 min, and 72°C for 2 min. PCR products were visualized by electrophoresis on 1% agarose gel.

RNA extraction and expression analysis

RNA was extracted from the anthers of T1 and BC5F2 plants using an Rneasy mini kit (Qiagen, USA). RNA was treated with RNase-free DNase (Promega, USA) to remove genomic DNA contaminants. Between 50 and 100 ng of total RNA was used to make single-strand cDNA using SuperScript III reverse transcriptase (Invitrogen, Toronto, Canada) in a 20 µl reaction with oligo (dT18) primers according to the manufacturer's instructions. RT-PCR was conducted using an AMV one step RT-PCR kit (Takara, Japan). Primers specific for S60 (5′- ATG AGA TAT GCT ACT TCT ATA TAT ACA -3′; 5′- TGA TTT AAC TTT GCA ACA GTA GCA -3′) were used for RT-PCR, and actin primers specific for Brassica (5′- ATG GCC GAG GCT GAT GAC AT -3′ and 5′- AGC CTC GGT AAG AAG AAC CG -3′) were used as a control. PCR was conducted using 50 ng of cDNA from the anthers of respective plants as templates in master mixes composed of 20 pmol of each primer, 150 µM of each dNTP, 1.2 U of Taq polymerase, 1x Taq polymerase buffer, and double-distilled H2O diluted to a total volume of 20 µl in 0.5 ml PCR tubes. The samples were then subjected to the following conditions: initial denaturation at 94°C for 5 min, followed by 30 cycles of denaturation at 94°C for 30 s, annealing at 58°C for 30 s and extension at 72°C for 1 min, with a final extension for 5 min at 72°C. PCR products were visualized by electrophoresis on 1% agarose gel.

Real-time quantitative PCR was performed using 1 µl of cDNA in a 25 µl reaction employing iTaq™ SYBR® Green Super-mix with ROX (California, USA). The same primers used for RT-PCR were employed for real-time PCR, while Brassica actin primers (5′- CAA CCA ATC GTC TGT GAC AA -3′; 5′- ATG TCT TGG CCT ACC AAC AA -3′) were used as a control. The conditions for real-time PCR were as follows: initial denaturation for 10 min at 95°C, followed by 40 cycles of 94°C for 30 s, 58°C for 30 s, and 72°C for 45 s. The fluorescence was measured following the last step of each cycle, and three replications were used for each sample. Amplification, detection, and data analysis were conducted using a Rotor-Gene 6000 real-time rotary analyzer (Corbett Life Science, Australia).

Self-compatibility analysis

SC was tested by observing the pollen-tube behavior and fruit set ratio. For observation of pollen-tube behavior, flower buds were covered with a bag one day before anthesis and then pollinated with pollen from the same plant, after which the samples were placed on 1% solid agar plate at room temperature for 24 h. The pollinated pistils were then separated from the buds and fixed in acetic alcohol (ethanol: acetic acid = 3∶1) for 5 h at room temperature. The fixed pistils were hydrolyzed in 1 N NaOH for 2h at 60 °C, after which they were stained with decolorized aniline blue solution (0.01% aniline blue in 2% K3PO4) for 2 h. The stained stigmas were mounted on glass slides with 50% glycerol and observed using a UV fluorescent microscope (Nikon Eclipse 80i, Japan). Three flowers were used from each plant. During flowering, the major inflorescence was bagged for self-pollination and the self-compatibility index (SCI) (self-compatibility index  =  the number of seeds/the number of flowers) was calculated according to Zhang et al. [40].

Results and Discussion

Identification of target gene and plasmid construction

The sequences of mature SP11 proteins are highly divergent, except for the presence of conserved cysteines [41]. Class II SP11s show several similarities with class I SP11s such as, they are small, secreted proteins with conserved putative signal peptides, their mature proteins show S-haplotype-specific polymorphisms in spite of their common hydrophilic and basic properties, like other pollen coat proteins (PCPs), they have eight conserved cysteine residues [20]. All eight cysteine residues are conserved in an arrangement that is characteristic of SP11. Conservation of the eight cysteine residues suggests a common three-dimensional protein structure of SP11s that is stabilized by intramolecular disulfide bonds, similar to the defensin family of antimicrobial proteins [42]. Two other residues that are conserved among most of the class I SP11s, a glycine residue between C1 and C2 and an aromatic amino acid residue between C3 and C4 [41], also are conserved among all of the class II SP11s. Hydrophilicity analysis of the class II SP11s does not suggest a hydrophilic (surface-exposed) structure in the C3–C4 region, which has been suggested for class I SP11s [43]. Three other regions, C1–C2, C2–C3, and C5–C6, are highly divergent and contain two or three amino acid residues that are completely variable across the four class II S-haplotypes of B. rapa [20]. B. rapa cv. Osome, a heterozygote of the S52 and S60 haplotype [44], was used for silencing of SP11/SCR gene based on RNA interference. S52 is a class I and S60 is a class II SP11 S-haplotype [9], [14]. Amino acid sequence identities among class II SP11s are 62.3 to 94.6%, rather high compared with those of class I SP11s which ranges from 19.5 to 76.1% [20]. Due to having low diversity and number of alleles of class II SP11s, we considered class II SP11, S60 as a crucial target for conversion of the SI B. rapa cv. Osome to SC. The 285 bp-length sequence of the S60 gene was selected as the target for RNA silencing.

We constructed an RNAi binary vector designated as S60-SP11RNAi (Figure 1) for introduction into SI B. rapa cv. Osome. The target sequence was incorporated into pBI101 in inverted-repeat orientations interrupted by GUS (1023bp). To enhance the efficiency of RNA silencing, strong S60 promoter (−489 to −1) was used to control expression of the transgene, and a plasmid was constructed to specifically suppress expression of S60. It is well known that the 5′-flanking region up to −192 bp SP11 is sufficient to direct gene expression in tapetum and pollen [45]. Moreover, the vector carried plant resistant gene NPTII driven by the Nos promoter and HPT driven by the CaMV 35S promoter.

Transformation and selection of homozygous transgenic plants

Using the hypocotyl as the explants, plasmid pBI101 with a S60-SP11RNAi cassette was transformed into B. rapa mediated by A. tumefaciens EHA105. Subsequently, the NPTII-resistant plantlets were regenerated, and the insertion of silencing fragments was confirmed by genomic DNA PCR (data not shown). Overall, 16 transgenic plants were obtained and grown to generate T1 seeds, after which they were harvested separately. Around, 20 to 25 T1 seeds were harvested from each transgenic line by selfing and then cultured in pots to produce T2 seeds. For genetic segregation of the T2 generation, 15–20 T2 seeds from each T1 line were tested onto 1/10 strength of MS plate with 50 mg/l hygromycin for antibiotic resistance screening. After five weeks of culture, heterozygote and homozygote events were counted on the basis of hygromycin resistance and susceptibility of plants. Among the 16 transformant lines, three lines produced transformants and non–transformants in a 3∶1 (χ2 = 0.003–0.04, P = 0.88–0.93) ratio. This pattern of segregation is only possible when single copy transfer of genes occurs in the transgenic plants. The other transformants were not segregated in such a pattern in the T2 generation, indicating the transfer of more than one copy of the target gene. Therefore, only homozygous transformant lines, self-compatible by RNAi (SR) 6–10, SR11-8 and SR18-7 were used for further study.

Suppression of S60 gene and conversion to self-compatible state

Real-time PCR was conducted to analyze the relative transcription level of the S60 gene among the three aforementioned homozygotic T2 lines. The results revealed much lower transcription levels in the transformants lines than the non-transformant line (Figure 2). Suppression levels varied from almost 2.5 to approximately 20-fold when compared with the untransformed control, with the highest suppression being observed in SR6-10. These three T2 lines were self-pollinated to produce T3 seeds and the compatibility of these lines was also tested by investigating the pollen tube behavior and fruit set. Unlike the self-incompatible tester line (Figure 3A), the S60-SP11RNAi lines showed numerous pollen tubes on stigma papilla cells upon self-pollination (Figure 3B-D). The mean SCI values of the SR6-10, SR11-8 and SR18-7 T2 lines were 6.89, 6.28 and 5.91, respectively, which were comparable to their average seed per open-pollinated flower values (7.95, 7.55 and 7.38). The mean SCI of the tester line was 0.82. Taken together, the pollen tube observations and fruit set data confirmed that the S60-SP11RNAi lines were highly self-compatible.

[Figure omitted. See PDF.]

Figure 2. Real-time PCR expression analysis of the S60-SP11 gene in T2 RNAi plants.

S60: wild type of S60, SR6-10, SR11-8 and SR18-7: RNAi transgenic plants.

https://doi.org/10.1371/journal.pone.0049497.g002

[Figure omitted. See PDF.]

Figure 3. Pollen-tube behaviors of self-pollinated T2 generation plants.

Photographs were obtained by UV fluorescence microscopy. A): WT, B): SR6-10 RNAi plant, C): SR11-8 RNAi plant, D): SR18-7 RNAi plant.

https://doi.org/10.1371/journal.pone.0049497.g003

The SI system is sporophytic in Brassica, and SP11s from class-II S-haplotypes (e.g. S40 and S60) exhibited strictly sporophytic expression patterns, suggesting that the expression of SP11 in the tapetal cell layer was sufficient for SI [20]. Recent biochemical studies revealed that SP11 functions as the sole ligand for its cognate SRK receptor complex. Their interaction induces the autophosphorylation of SRK, which is expected to trigger the signaling cascade that results in the rejection of self-pollen. This so-called ligand receptor complex interaction and receptor activation occurs in an S-haplotype-specific manner, and this specificity is almost certainly the basis for self-pollen recognition [5], [15]. In this study, we silenced the S60 gene in B. rapa to interfere with the so-called ligand receptor interaction and receptor activation and prevent consequent rejection of self-pollen and the resultant S60-SP11RNAi homozygous plants were SC.

Stable inheritance of the trait

To determine the stability of gene silencing in the following generations, three T2 lines that had relatively strong gene silencing based on the reduction of the transcript abundance levels were selected for further analysis. All T3 individuals produced from each T2 line were checked through amplification of NPTII and S60 primers using gDNA extracted from 4-week-old leaves to identify positive and negative progenies and all plants were found to be positive (data not shown). S60 gene expression in these plants was also analyzed through RT-PCR analysis of three plants selected at random from each T3 line using S60 specific primers (Figure 4). The results showed that the gene suppression level was similar to that of T1 transgenic plants when compared to the untransformed control.

[Figure omitted. See PDF.]

Figure 4. Real-time PCR expression analysis of the S60-SP11 gene in T3 generation plants.

WT: wild type of S60, lane 1–3: SR6-10 RNAi plants, lane 4–6: SR11-8 RNAi plants and lane 7–9: SR18-7 RNAi plants.

https://doi.org/10.1371/journal.pone.0049497.g004

Performance evaluation of RNAi transgenic lines

Transgenic SC B. rapa lines were subjected to performance evaluation as seed set experiments under two different conditions; use of pollinating agents (Bees and CO2) and no pollinating agents. It should be noted that the commercial cultivation of different Brassica crops uses bees and CO2 as pollinating agents, and CO2 is not environment friendly. In this evaluation study, three homozygous RNAi transgenic lines and a non-transformant (WT) line were used. For the non-transformant plants, bees and CO2 were used as pollinating agents, while for the transformant plants no bees or CO2 were used (File S1). In each case, the total seed number was counted from four plants and highest average seed set was found in the SR6-10 RNAi transgenic line (Table 1). These findings indicate that the transgenic RNAi lines could perform well in fruit settings without any pollinating agents and could be used for commercial cultivation with great economic benefits and minimal environmental hazards.

[Figure omitted. See PDF.]

Table 1. Comparative seed production in non-transgenic control and S60-SP11RNAi plants of T3 generation.

https://doi.org/10.1371/journal.pone.0049497.t001

Use of RNAi transgenic line as breeding material

SI is generally the main obstruction for B. rapa seed production and for that SC cultivar is crucial issue in case of commercial seed production. In this study, we developed a SC B. rapa line and then transferred this trait in a strictly SI commercial variety, B. rapa ‘Seoulbechhu’. We utilized the backcross method to transfer the S60-SP11RNAi gene from the SR6-10 RNAi transgenic homozygous line into this variety and obtained 19 lines in the BC5F2 generation. The insertion of silencing fragments was confirmed by genomic DNA PCR of these 19 lines using hygromycin resistant gene primers (data not shown). We then re-checked the inserts using class I and class II SRK universal primers to identify the homozygous lines for the S60 gene introduced into these 19 lines and found no class I insert in lines no. 7, 10, 12 and 17 (File S2), indicating that these four lines are homozygous for class II S60-SP11RNAi gene. Because pollen-dominant class I S haplotypes (such as S8, S9, S12 and S52 in Brassica rapa) are generally co-dominant, they are almost always dominant over recessive class II S haplotypes (such as S44, S60, S40 and S29) [14], [21]. Upon expression analysis, only these four lines showed expression of the S60-SP11RNAi gene at different levels, while no expression was observed in the other lines (File S3). The expression level of the S60-SP11RNAi gene in these four homozygous lines was compared with that of non-transformed S60 containing plants. Line 10 showed the highest level of suppression when compared to non-transformed S60 bearing plants, followed by lines 7, 17 and 12 (Figure 5). The fruit setting of these four lines was also very high in this generation.

Taken together, it can be concluded that the transgenic RNAi lines were fully self-compatible, and that this trait was stable in subsequent generations, even after crossing with other commercial lines with higher performance than the non-transgenic lines.

[Figure omitted. See PDF.]

Figure 5. Expression analysis of S60-SP11RNAi gene in BC5F2 generation plants using S60-SP11 specific primers and Br-Actin primers were used as control.

A) RT-PCR and B) Real-time PCR expression analysis. WT: wild type of S60; 7, 10, 12 and 17: Homozygous lines of BC5F2 generation after backcrossing between SR6-10 transgenic plant and B. rapa ‘Seoulbechhu’ plant.

https://doi.org/10.1371/journal.pone.0049497.g005

Supporting Information

[Figure omitted. See PDF.]

File S1.

Fruit set analysis in S60-SP11RNAi plants of T3 generation. A) Non-transgenic control (with CO2 & Bees), B) SR6-10 lines (without CO2 & Bees), C) SR11-8 lines (without CO2 & Bees), D) SR18-7 lines (without CO2 & Bees).

https://doi.org/10.1371/journal.pone.0049497.s001

(TIF)

File S2.

Genomic DNA PCR analysis of 19 BC5F2 (backcrossing between SR6-10 transgenic plant and B. rapa ‘Seoulbechhu’ plant) generation plants using SRK, A) class I and B) class II universal primers.

https://doi.org/10.1371/journal.pone.0049497.s002

(TIF)

File S3.

RT-PCR expression analysis of S60-SP11RNAi gene in 19 BC5F2 generation plants using S60-SP11 specific primers and Br-Actin primers were used as control. WT: wild type of S60, 1–19: lines of BC5F2 generation.

https://doi.org/10.1371/journal.pone.0049497.s003

(TIF)

Author Contributions

Conceived and designed the experiments: Hee-Jeong Jung Hyo-Jin Jung NUA, JIP, ISN. Performed the experiments: Hee-Jeong Jung Hyo-Jin Jung NUA. Analyzed the data: Hee-Jeong Jung Hyo-Jin Jung NUA JIP KKK YKH YPL ISN. Contributed reagents/materials/analysis tools: Hee-Jeong Jung Hyo-Jin Jung NUA JIP KKK YKH YPL ISN. Wrote the paper: NUA JIP.

Citation: Jung H-J, Jung H-J, Ahmed NU, Park J-I, Kang K-K, Hur Y, et al. (2012) Development of Self-Compatible B. rapa by RNAi-Mediated S Locus Gene Silencing. PLoS ONE 7(11): e49497. https://doi.org/10.1371/journal.pone.0049497

References

1. McCubbin AG, Kao TH (2000) Molecular recognition and response in pollen and pistil interactions. Annu Rev Cell Dev Biol 16: 333–364.

2. Watanabe M, Suzuki G, Takayama S (2008) Chapter 7 "Milestones identifying self-incompatibility genes in Brassica species: From old stories to new findings", In "Self-incompatibility in Flowering Plants - Evolution, Diversity, and Mechanisms" (Ed.: V. E. Franklin-Tong), Springer, pp151–172.

3. Suwabe K, Suzuki G, Watanabe M (2010) Achievement of genetics in plant reproduction research: past decade for coming decade. Genes Genet Syst 85: 297–310.

4. Suzuki G, Kai N, Hirose T, Fukui K, Nishio T, et al. (1999) Genomic organization of the S locus: identification and characterization of genes in SLG/SRK region of S9 haplotype of Brassica campestris (syn. rapa). Genetics 153: 391–400.

5. Takayama S, Shiba H, Iwano M, Shimosato H, Che FS, et al. (2000) The pollen determinant of self-incompatibility in Brassica campestris. P?Natl Acad Sci, USA 97: 1920–1925.

6. Schopfer CR, Nasrallah ME, Nasrallah JB (1999) The male determinant of self-incompatibility in Brassica. Science 286: 1697–1700.

7. Stein JC, Howlett B, Boyes DC, Nasrallah ME, Nasrallah JB (1991) Molecular cloning of a putative receptor protein kinase gene encoded at the self-incompatibility locus of Brassica oleracea. Proc Natl Acad Sci, USA 88: 8816–20.

8. Watanabe M, Takasaki T, Toriyama K, Yamakawa S, Isogai A, et al. (1994) A high degree of homology exists between the protein encoded by SLG and the S receptor domain encoded by SRK in self-incompatible Brassica campestris L. Plant Cell Physiol 35: 1221–1229.

9. Takasaki T, Hatakeyama K, Suzuki G, Watanabe M, Isogai A, et al. (2000) The S receptor kinase determines self-incompatibility in Brassica stigma. Nature 403: 913–916.

10. Suzuki G, Watanabe M, Nishio T (2000) Physical distances between S-locus genes in various S haplotypes of Brassica rapa and B. oleracea. Theor Appl Genet 101: 80–85.

11. Nou IS, Watanabe M, Isogai A, Shiozawa H, Suzuki A, et al. (1991) Variation of S alleles and S glycoproteins in a naturalized population of self incompatible Brassica campestris L. Jpn J Genet 66: :227–239.

12. Nou IS, Watanabe M, Isogai A, Hinata K (1993) Comparison of S alleles and S glycoproteins between two wild populations of Brassica campestris in Turkey and Japan. Sex Plant Reprod. 67: :79–86.

13. Nou IS, Watanabe M, Isuzugawa K, Isogai A, Hinata K (1993) Isolation of S allele from a wild population of Brassica campestris L. at Balcesme, Turkey and their characterization by S glycoprotein. Sex Plant Reprod 6: 71–78.

14. Hatakeyama K, Watanabe M, Takasaki T, Ojima K, Hinata K (1998) Dominance relationships between S-alleles in self-incompatible Brassica campestris L. Heredity 80: 241–247.

15. Takayama S, Shimosato H, Shiba H, Funato M, Che FS, et al. (2001) Direct ligand-receptor complex interaction controls Brassica self-incompatibility. Nature 413: 534.

16. Kachroo A, Schopfer CR, Nasrallah ME, Nasrallah JB (2001) Allele-specific receptor-ligand interactions in Brassica self-incompatibility. Science 293: 1824–1826.

17. Shimosato H, Shiba H, Iwano M, Entani T, Che FS, et al. (2007) Characterization of the SP11/SCR high-affinity binding site involved in self/nonself recognition in Brassica self-incompatibility. Plant Cell 19: 107–117.

18. Murase K, Shiba H, Iwano M, Che FS, Watanabe M, et al. (2004) A membrane-anchored protein kinase involved in Brassica self-incompatibility signaling. Science 303: 1516–1519.

19. Kakita M, Murase K, Iwano M, Matsumoto T, Watanabe M, et al. (2007) Two distinct forms of M-locus protein kinase localize to the plasmamembrane and interact directly with S-locus receptor kinase to transduce self-incompatibility signaling in Brassica rapa. Plant Cell 19: 3961–3973.

20. Shiba H, Iwano M, Entani T, Ishimoto K, Shimosato H, et al. (2002) The dominance of alleles controlling self-incompatibility in Brassica pollen is regulated at the RNA level. Plant Cell 14: 491–504.

21. Kakizaki T, Takada Y, Ito A, Suzuki G, Shiba H, et al. (2003) Linear dominance relationship among four class-II S haplotypes in pollen is determined by the expression of SP11 in Brassica self-incompatibility. Plant Cell Physiol 44: 70–75.

22. Nasrallah JB, Nasrallah ME (1993) Pollen–stigma signaling in the sporophytic self-incompatibility response. Plant Cell 5: 1325–1335.

23. Shiba H, Kakizaki T, Iwano M, Tarutani Y, Watanabe M, et al. (2006) Dominance relationships between self-incompatibility alleles controlled by DNA methylation. Nature Genet 38: 297–299.

24. Tarutani Y, Shiba H, Iwano M, Kakizaki T, Suzuki G, et al. (2010) Trans-acting small RNA determines dominance relationships in Brassica self-incompatibility. Nature 466: 983–986.

25. Fukuda H, Higashiyama T (2012) Diverse Functions of Plant Peptides: Entering a New Phase. Plant Cell Physiol 52: 1–4.

26. Higashiyama T (2010) Peptide Signaling in Pollen–Pistil Interactions. Plant Cell Physiol 51: 177–189.

27. Goto H, Okuda S, Mizukami A, Mori H, Sasaki N, et al. (2011) Chemical Visualization of an Attractant Peptide, LURE. Plant Cell Physiol 52: 49–58.

28. Hancock CN, Kent L, McClure B (2005) The 120 kDa glycoprotein is required for S-specific pollen rejection in Nicotiana. The Plant Journal 43: 716–723.

29. Hunter C, Poethig R (2003) Missing links: miRNAs and plant development. Curr Opin Genet Dev 13: 372–378.

30. Waterhouse P, Helliwell C (2003) Exploring plant genomes by RNA-induced gene silencing. Nature Genet 4: 29–38.

31. Wesley S, Helliwell C, Smith N, Wang MB, Rouse DT, et al. (2001) Construct design for efficient, effective and high-throughput gene silencing in plants. Plant J 27: 581–590.

32. Urayama S, Moriyama H, Aoki N, Nakazawa Y, Okada R, et al. (2010) Knock-down of OsDCL2 in Rice Negatively Affects Maintenance of the Endogenous dsRNA Virus, Oryza sativa Endornavirus. Plant Cell Physiol 51: 58–67.

33. Yadav SR, Khanday I, Majhi BB, Veluthambi K, Vijayraghavan U (2011) Auxin-Responsive OsMGH3, a Common Downstream Target of OsMADS1 and OsMADS6, Controls Rice Floret Fertility. Plant Cell Physiol 52: 2123–2135.

34. Wakasa Y, Hirano K, Urisu A, Matsuda T, Takaiwa F (2011) Generation of Transgenic Rice Lines with Reduced Contents of Multiple Potential Allergens Using a Null Mutant in Combination with an RNA Silencing Method. Plant Cell Physiol 52: 2190–2199.

35. Nagamine A, Matsusaka H, Ushijima T, Kawagoe Y, Ogawa M, et al. (2011) A Role for the Cysteine-Rich 10 kDa Prolamin in Protein Body I Formation in Rice. Plant Cell Physiol 52: 1003–1016.

36. Chuang CF, Meyerowitz EM (2000) Specific and heritable genetic interference by double-stranded RNA in Arabidopsis thaliana. Proc Natl Acad Sci, USA 97: 4985–4990.

37. Jefferson RA (1987) Assaying chimeric genes in plants: the GUS gene fusion system. Plant Mol Biol Rep 5: 387–405.

38. Radke SE, Turner JC, Facciotti D (1992) Transformation and regeneration of Brassica rapa using Agrobacterium tumefaciens. Plant Cell Reports 11(10): 499–505.

39. Zhang FL, Takahata Y, Watanabe M, Xu JB (2000) Agrobacterium mediated transformation of cotyledonary explants of Chinese cabbage (Brassica campestris L. ssp. pekinensis). Plant Cell Rep 19: 569–575.

40. Zhang XG, Ma CZ, Fu TD, Li YY, Wang TH, et al. (2008) Development of SCAR markers linked to self-incompatibility in Brassica napus L. Mol Breeding 21: 305–315.

41. Watanabe M, Ito A, Takada Y, Ninomiya C, Kakizaki T, et al. (2000) Highly divergent sequences of the pollen self-incompatibility (S) gene in class-I S haplotypes of Brassica campestris (syn. rapa) L. FEBS Lett 473: 139–144.

42. Broekaert WF, Terras FRG, Cammue BPA, Osborn RW (1995) Plant defensins: Novel antimicrobial peptides as components of the host defense system. Plant Physiol 108: 1353–1358.

43. Schopfer CR, Nasrallah JB (2000) Self-incompatibility:Prospects for a novel putative peptide-signaling molecule. Plant Physiol 124: 935–939.

44. Takasaki T, Hatakeyama K, Ojima K, Watanabe M, Toriyama K, et al. (1997) Factors influencing Agrobacterium-mediated transformation of Brassica rapa. L Breed Sci 47: 127–134.

45. Shiba H, Takayama S, Iwano M, Shimosato H, Funato M, et al. (2001) A Pollen Coat Protein, SP11/SCR, Determines the Pollen S-Specificity in the Self-Incompatibility of Brassica Species. Plant Physiology 125: 2095–2103.