Plant materials
Onions showing unstable MS phenotypes in our previous study (Kim et al. 2019a) were used. Two of them in the OP2 population were crossed with onions showing stable MS phenotypes and containing CMS-T cytoplasm. Two resulting populations were designated as TUMS4 and TUMS9, respectively. Male-fertile (MF) plants in these two populations were self-pollinated or cross-pollinated with MS plants in the same populations to produce large-sized segregating populations. A detail pedigree of onion populations used in this study is depicted in Supplementary Fig. 1.
Seedlings germinated in 128-cell plug trays were transplanted into pots in the greenhouse or fields of Chonnam National University, Gwangju, Republic of Korea (35° 17´ N, 126° 90´ E). Floral induction was initiated naturally during winter seasons. Male-fertility phenotypes were determined by visual examination in June for three years (2018–2020). Statistics of temperatures during flowering seasons are summarized in Supplementary Table 1. Morphologies of anthers and pollen grains were examined using a dissecting microscope (Stemi 2000-C; ZEISS, Oberkochen, Germany) and a scanning electron microscope (SEM, JSM-IT300; JEOL, Tokyo, Japan), respectively.
RNA extraction, RNA-Seq analysis, and real-time RT-PCR
To perform a combined analysis of BSA and RNA-seq, flower buds of 10 MF and 10 MS individuals of the TUMS4 population were pooled, respectively. Total RNAs were extracted from each bulked sample using an RNeasy Plant Mini Kit (QIAGEN, Valencia, CA, USA). In the first RNA-Seq analysis, diverse sizes of flower buds larger than 3 mm were used. Flower buds smaller than 3 mm were included in the second analysis to cover early stages of male gametophyte development. The third analysis was performed using the TUMS4-S1 population after removing heterozygous MF individuals. Genotypes of heterozygous individuals were predicted using tightly linked molecular markers developed in this study.
RNA-Seq analyses were performed by a specialized company (Phyzen Genomics Institute, Seoul, Republic of Korea). Transcriptomes were sequenced using a HiSeq X Ten platform (Illumina, Hayward, CA, USA). Trimmed reads were separately mapped to two reference transcriptomes (Kim et al. 2015b; Sohn et al. 2016). Unmapped reads were pooled and de novo assembled into contigs using Trinity software (Haas et al. 2013). Identification of single nucleotide polymorphisms (SNPs) between bulked RNAs and quantification of transcription levels of contigs were carried out using SAMTools (Li et al. 2009) and RSEM software (Li and Dewey 2011), respectively.
To estimate transcription levels of orf725, real-time RT-PCRs were performed. Total RNAs were extracted from flower buds of five MF and five MS plants of the TUMS4 population using an RNeasy Plant Mini Kit (QIAGEN). Total RNAs extracted in the previous study (Kim et al. 2015a) from MF and MS F2 individuals with male-fertility phenotypes controlled by the Ms locus were also used for comparison of transcription levels. cDNAs were synthesized using a cDNA synthesis kit (SuperScript™ III first-strand synthesis system for RT-PCR, Invitrogen, Carlsbad, CA, USA).
Real-time RT-PCRs were carried out using SYBR® Green Realtime PCR Master Mix (Toyobo Co. Ltd, Osaka, Japan) and a LightCycler®96 Real-Time PCR system (Roche Molecular Systems, Pleasanton, CA, USA) according to manufacturers’ instructions with four technical replicates. cDNAs diluted 100-fold were then used as templates. Onion nad6 gene positioned in the master circle of mitochondrial genome (GenBank accession: KU318712) was used as an internal control. Primer sequences used for real-time RT-PCR are shown in Supplementary Table 2.
DNA extraction, PCR amplification, and high-resolution melting (HRM) marker analysis
Total genomic DNAs of segregating populations were extracted from leaves or flower peduncle tissues using the cetyltrimethylammonium bromide (CTAB) method (Doyle and Doyle 1987). For sequencing of contigs, PCR amplifications were performed in a 25 µL reaction mixture containing 0.1 µg template, 2.5 µL 10x PCR buffer, 0.5 µL forward primer (10 µM), 0.5 µL reverse primer (10 µM), 0.5 µL dNTPs (10 mM each), and 0.25 µL polymerase mix (Advantage 2 Polymerase Mix, Takara Bio, Shiga, Japan). PCR amplification consisted of an initial denaturation step at 95°C for 4 min, 10 cycles at 95°C for 30 s, 65°C (0.8°C decrements in each cycle) for 30 s, and 72°C for 1 min, 35 cycles at 95°C for 30 s, 57°C for 30 s, and 72°C for 1 min, and a final 10 min extension at 72°C. After confirming successful PCR amplification by visualizing on 1.5% agarose gels after ethidium bromide staining, PCR products were purified for sequencing using a QIAquick PCR Purification kit (QIAGEN). Sequencing was performed by a specialized company (Macrogen, Seoul, Republic of Korea). Regarding HRM and cleaved amplified polymorphic sequence (CAPS) markers, detail procedures were described in previous studies (Kim and Kim 2019; Cho et al. 2021). Primer sequences of molecular markers are shown in Supplementary Table 2.