Physical Location of New PCR-Based Markers and Powdery Mildew Resistance Gene(s) on Rye (Secale cereale L.) Chromosome 4 Using 4R Dissection Lines.

Rye (Secale cereale L.) 4R chromosome contains elite genes that are applicable for wheat (Triticum aestivum L.) cultivar improvement. PCR-based 4R-specific markers can benefit the detection of elite genes on 4R in wheat backgrounds. In this study, a new fluorescence in situ hybridization (FISH) map of the 4RKu chromosome of rye Kustro has been constructed. A set of 4RKu dissection lines was obtained and 301 new 4RKu-specific markers were developed using specific length amplified fragment sequencing (SLAF-seq) technology. These markers were combined with the 99 4RKu-specific markers previously developed, and were physically mapped to 4RKu chromosome using the new FISH map and the 4RKu dissection lines. A total of 338 of the 400 markers have been successfully mapped to six regions of 4RKu chromosome. Additionally, the powdery mildew resistance gene(s) on the 4RLKu arm was located to the segment between L.4 and L.8, the same region where 115 4RLKu-specific markers were mapped. The markers developed in this study can be used to identify a specific segment of 4R chromatin in wheat backgrounds, help construct a high-density physical map of 4R chromosome, and facilitate the utilization of elite genes on 4R chromosome in wheat breeding programs.

Introduction

Rye (Secale cereale L.) is an important gene source for wheat (Triticum aestivum L.) cultivar improvement. Only the short arm of rye chromosome 1R has been widely used to develop wheat cultivars through 1BL.1RS or 1AL.1RS translocation chromosomes (Berzonsky et al., 1991; Rabinovich, 1998; Kumar et al., 2003; Landjeva et al., 2006). In fact, other rye chromosomes also contain elite genes that can be used for wheat improvement. For example, 4R chromosome contains disease- and insect-resistance genes. Lukaszewski et al. (2001) reported that the 4RL arm of Secale montanum Guss carried a Russian wheat aphid (RWA) resistance gene. New powdery mildew resistance gene(s), which differ from previously reported rye genes, were located on 4R chromosome of rye cultivar German White (An et al., 2013). It has been reported that the 4RL arm of rye Kustro also possesses powdery mildew resistance gene(s) (Fu et al., 2014). Additionally, the 4R chromosome of rye Imperial carries at least two genetic factors that have positive effects on wheat pollination traits (Nguyen et al., 2015). The addition of ‘Kriszta’ chromosome 4R to wheat genome can increase the total protein content (Schneider et al., 2016). However, the utilization of elite genes on 4R chromosome for wheat improvement is difficult because of compensation issues and the very low recombination frequency of the 4R chromosome with its wheat homoeologues (Lukaszewski et al., 2001). Because of these issues, the only available approaches are constructing wheat-rye small-segment translocation lines and cloning rye elite genes. Elite genes can be utilized efficiently when rye chromosomal segments have been transferred into wheat and precisely identified. A well-saturated molecular linkage map can be used for gene tagging. Rye chromosome-specific markers are beneficial to the effective application of rye elite genes in wheat breeding programs, however, very few 4R-specific markers, especially PCR-based and agarose gel electrophoresis-based markers, have been developed so far.

In this study, New PCR-based and 4RKu-specific markers using specific length amplified fragment sequencing (SLAF-seq) technology were developed. Subsequently, these new markers were physically mapped onto six regions (bins) on 4RKu using 4RKu dissection lines in a wheat background.

Materials and Methods

Plant Materials

The octoploid triticale line MK was developed by crossing common wheat T. aestivum L. Mianyang 11 (MY11) with rye S. cereale L. Kustro. Progeny were produced by controlled backcrossing of MK with MY11, followed by self-fertilization. From these progeny, seven wheat-rye monosomic addition lines (MA1RKu-MA7RKu lines), a 4RSKu monotelosomic addition line (MTA4RSKu) and a 4RLKu monotelosomic addition line (MTA4RLKu) were detected (Li et al., 2016; Qiu et al., 2016). Some of the MK seeds were irradiated with fast neutrons at the Institute of Nuclear Physics and Chemistry, China Academy of Engineering Physics, Mianyang, China. The irradiated MK seeds were used as recipients to cross with common wheat T. aestivum L. Chuannong 27 (CN27), and line 12FT2115 was selected from the progeny of this cross combination (Fu et al., 2014). Some of the selfed progeny (seeds) of lines 12FT2115 and MTA4RLKu were irradiated with 60Co-γ rays at the Biotechnology and Nuclear Technology Research Institute, Sichuan Academy of Agricultural Sciences, China. Common wheat T. aestivum L. Chinese Spring (CS) was used as a control.

Cytological Techniques and In Situ Hybridization

Non-denaturing fluorescence in situ hybridization (ND-FISH) technology was used to analyze root-tip metaphase cells. Oligonucleotide probes containing Oligo-1162, Oligo-pSc200, Oligo-pSc250, Oligo-pSc119.2-1, and Oligo-pTa535-1 were synthesized following the methods described by Tang et al. (2014) and Fu et al. (2015). ND-FISH was carried out following the procedure described by Fu et al. (2015). Probes Oligo-1162, Oligo-pSc200, Oligo-pSc250 and Oligo-pTa535-1 were 5′-end-labeled with 6-carboxytetramethylrhodamine (Tamra). Probe Oligo-pSc119.2-1 was 5′-end-labeled with 6-carboxyfluorescein (6-FAM). Additionally, a synthetic oligonucleotide probe (AAC)6 was used and was 5′-end-labeled with Cy5. Metaphase chromosomes of the root-tips were prepared following the methods described by Han et al. (2006). Images were made using an epifluorescence Olympus BX51 microscope, which was equipped with a cooled charge-coupled device camera and with the HCIMAGE Live software (version 2.0.1.5). Images were processed using Adobe Photoshop CS 3.0.

Development of PCR-Based Markers

Genomic DNAs of S. cereale L. Kustro and MA4RKu were sequenced using the SLAF-seq technique (Biomarker, Beijing, China). The sequencing procedure followed the methods described by Chen et al. (2013), with some modifications. Genomic DNAs of Kustro and MA4RKu were digested using the restriction enzyme, HaeIII. Subsequently, a Quick Spin column (Qiagen) was used to purify the samples and then run out on a 2% agarose gel. Fragments between 450 to 500 bp were isolated using a Gel Extraction Kit (Qiagen). These isolates were used in a PCR reaction described by Chen et al. (2013). Amplicons with the sizes between 450 to 500 bp were excised and diluted for sequencing, and they were identified, filtered, clustered and corrected following the methods described by Chen et al. (2013). The pair-end reads derived from Kustro and MA4RKu were compared with wheat A genome, D genome and T. aestivum L. Chinese Spring (supported by Biomarker, Beijing, China) sequences using SOAP software (Li et al., 2009). The pair-end reads with low wheat homology were kept. Finally, after comparing specific pair-end reads of Kustro and MA4RKu, the 4RKu specific pair-end reads were obtained. Primers were designed according to the 4RKu specific pair-end reads using the software Primer 3 (version 4.0). For primers, optimal melting temperature and size values were set to 60°C and 20 bases, respectively.

Physical Location of 4R-Specific Markers

Markers, whose products presented in Kustro and MA4RKu, but were absent in CS, MY11, CN27, MA1RKu-MA3RKu and MA5RKu-MA7RKu, were regarded as 4RKu-specific markers. These markers were located to specific regions of the 4RKu chromosome using the 4RKu dissection lines. Additionally, 99 4RKu-specific markers that were previously developed by Qiu et al. (2016) were located on specific regions of the 4RKu chromosome using the 4RKu dissection lines.

PCR Analysis

The PCR amplifications were carried out according to the procedure described by Li et al. (2016). The amplicons were electrophoresed on 2% agarose gels in 1 × TAE buffer. For each of the primer pairs used in this study, PCR reactions were repeated three times.

Powdery Mildew Resistance Test

The resistance of 4RKu dissection lines and parental wheat MY11 and CN27 to powdery mildew was evaluated. Plants were grown in two growing seasons (2015-2016 and 2016-2017) in Qionglai, Sichuan, China. The materials were naturally infected by powdery mildew, and infection types (IT) were scored according to the standard described by Fu et al. (2014).

Results

FISH Map of Chromosome 4R

Repetitive DNA sequence (AAC)6, combined with three oligonucleotides Oligo-pSc119.2-1, Oligo-pSc200 and Oligo-pSc250, were used as probes to analyze the root tip metaphase chromosomes of lines MA4RKu, MTA4RSKu, and MTA4RLKu. The signals of probes Oligo-pSc119.2-1, Oligo-pSc200 and Oligo-pSc250 on the telomere of 4RSKu were very strong. Oligo-pSc119.2-1 produced a strong signal at the interstitial site of 4RLKu, and Oligo-pSc200 and Oligo-pSc250 produced clear signals at the sub-telomeric region of 4RLKu (Figure ). Probe (AAC)6 had two signal sites on 4RSKu, one signal site on the pericentromeric region of 4RLKu and one signal site on the telomere of 4RLKu (Figure ). Based on an in situ hybridization map of 4R chromosome constructed by Cuadrado et al. (1995), the signal sites of probes (AAC)6, Oligo-pSc119.2-1, Oligo-pSc200 and Oligo-pSc250 were numbered, and a new FISH map of 4R chromosome was constructed (Figure ).

ND-FISH analysis using probes Oligo-pSc119.2-1 (green), Oligo-pTa535-2 (red), Oligo-pSc200 (red), Oligo-pSc250 (red) and (AAC)6 (white) to identify MA4RKu, MTA4RSKu and MTA4RLKu. (A) Line MA4RKu contains 42 wheat chromosomes and a 4R chromosome. (B) Line MTA4RSKu contains 42 wheat chromosomes and a 4RS arm. (C) Line MTA4RLKu contains 42 wheat chromosomes and a 4RL arm. (D) Cut-pasted 4RKu chromosome, 4RSKu arm and 4RLKu arm. A new FISH map was constructed according to Oligo-pSc119.2-1, Oligo-pTa535-2, Oligo-pSc200, Oligo-pSc250 and (AAC)6 signals. Chromosomes were counterstained with DAPI (blue). Scale bar is 10 μm.

Isolation of 4RKu Dissection Lines

Line 16T75-24, containing a 5BS.5BL-4RL translocation chromosome, was detected from the irradiated seeds of line MTA4RLKu (Figure ). A homozygous 4RS-5DS.5DL translocation line, 16T197-6, was obtained from the selfed progeny of irradiated 12FT2115 (Figure ), and three kinds of broken 5DS-4RS.4RL chromosomes were detected in lines 16T196-22, 16T177-4 and 16T175-1, respectively (Figures ). On the 5BS.5BL-4RL translocation chromosome, the breakpoint of 4RLKu was located at the Oligo-pSc119.2 signal site and the segment from L.4 to L.8 was retained (Figures ). On the 4RS-5DS.5DL chromosomes, the segments from S.4 to S.6 of the 4RKu chromosome were kept (Figures ). Breakpoints on the pair of 5DS-4RS.4RL chromosomes in line 16T196-22 were located at the Oligo-pSc119.2 signal sites of the 4RLKu, and the segments between L.4 and L.8 were lost as were segments between S.4 and S.6 of the short arm (Figures ). The breakpoints on the pair of 5DS-4RS.4RL chromosomes in line 16T177-4 were located on the regions between L.2 and L.4, and the segments between the breakpoints and L.8 were absent as were segments of the short arm from S.4 to S.6 (Figures ). In line 16T175-1, the 5DS-4RS.4RL chromosomes were broken at the centromeres, and this line contained two 5DS-4RS small translocation chromosomes (Figures ). Therefore, the wheat–rye translocation chromosomes in lines 16T75-24, 16T197-6, 16T196-22, 16T177-4 and 16T175-1 composed a set of 4RKu dissection lines, and the 4RKu chromosome was divided into five regions (Figure ).

ND-FISH analysis using probes Oligo-pSc119.2-1 (green), Oligo-pTa535-2 (red), Oligo-pSc200 (red), Oligo-pSc250 (red) and (AAC)6 (white) to identify 4RKu dissection lines. (A) 5BS.5BL-4RL translocation line 16T75-24. (B) 4RS-5DS.5DL translocation line 16T197-6. (C) Line 16T196-22 contains broken 5DS-4RS.4RL chromosomes. (D) Line 16T177-4 contains broken 5DS-4RS.4RL chromosomes. (E) Line 16T175-1 contains 5DS.4RS small translocation chromosomes. (F) Cut-pasted 4RKu dissections and the schematic diagram of FISH pattern of the 4RKu chromosome. The FISH signals are numbered. In the schematic diagram, five regions of the 4RKu chromosome were indicated by Roman numerals, green bands represent Oligo-pSc119.2-1 signals, red bands represent Oligo-pSc200 and Oligo-pSc250 signals, yellow bands represent (AAC)6 signals, and ‘Cen’ represents centromere. Chromosomes were counterstained with DAPI (blue). Scale bar is 10 μm.

Development of 4R-Specific Markers

Seven wheat–rye monosomic addition lines including MA1RKu, MA2RKu, MA3RKu, MA4RKu, MA5RKu, MA6RKu, and MA7RKu were used to identify 4RKu-specific markers. From the 33,577 4RKu-specific pair-end reads, 780 reads were randomly selected for designing primers. Three hundred and one of the 780 primer pairs amplified specific bands from Kustro and MA4RKu, but not from CS, MY11, CN27, MA1RKu-MA3RKu and MA5RKu-MA7RKu (Figure ). The 301 primer pairs were regarded as 4RKu-specific markers (Supplementary Table S1). Using lines MTA4RSKu and MTA4RLKu, 128 and 173 of the 301 markers were located on the 4RSKu and 4RLKu arms, respectively (Figure and Supplementary Table S2).

Identification of 4RKu-specific, 4RSKu-specific and 4RLKu-specific markers using CS, MY11, CN27, Kustro, MTA4RSKu, MTA4RLKu and MA1RKu–7RKu. (A) Products amplified by primer pair KU-4R.429. (B) Products amplified by primer pair KU-4R.732. (C) Products amplified by primer pair KU-4R.012. (D) Products amplified by primer pair KU.170. Primer pairs KU-4R.429 and KU-4R.732 represent the 4RSKu-specific markers. Primer pairs KU-4R.012 and KU.170 represent the 4RLKu-specific markers. M: DNA marker. CS: Chinese Spring. MY11: Mianyang 11. CN27: Chuannong 27. Kustro: rye kustro. 1R-7R: MA1RKu–7RKu. 4RS: MTA4RSKu. 4RL: MTA4RLKu. Arrows indicate the target bands.

Physical Mapping of 4RKu-Specific Markers

The 301 4RKu-specific markers obtained in this study were physically located on six regions of the 4RKu chromosome using lines 16T75-24, 16T175-1, 16T177-4, 16T196-22 and 16T197-6 (Figure and Supplementary Table S2). Additionally, the 31 4RSKu-specific and 68 4RLKu-specific markers developed by Qiu et al. (2016) were also mapped to the six regions of the 4RKu chromosome (Figure and Supplementary Table S2). Among the 159 (128 + 31) 4RSKu-specific markers, 128 amplified 4RSKu-specific bands from line 16T197-6, but not the other four lines, and 15 amplified 4RSKu-specific bands from lines 16T175-1, 16T177-4, and 16T196-22, but not lines 16T197-6 and 16T75-24 (Figures ). Therefore, the 128 markers were mapped to region I and the 15 markers were mapped to region II (Figure ). There were 16 4RSKu-specific markers that could not amplify their products from any of the five dissection lines, therefore, they could not be mapped to any region of the 4RSKu arm (Supplementary Table S2). For the 241 (173 + 68) 4RLKu-specific markers, 49 markers only amplified their special products from line 16T77-4, carrying the proximal long arm segment from the centromere to L.2, and line 16T196-22, carrying the proximal long arm segment from the centromere to L.4 (Figure and Supplementary Table S2), six markers amplified 4RLKu-specific bands from line 16T196-22 (Figure and Supplementary Table S2), 115 markers only amplified their specific bands from line 16T75-24, which carries the long arm distal segment from L.4 to the telomere (Figure and Supplementary Table S2), and 25 markers amplified 4RLKu-specific bands from both the lines 16T75-24 and 16T196-22 but not the lines 16T175-1, 16T177-4 and 16T197-6 (Figure and Supplementary Table S2). Therefore, the 49, six and 115 markers were mapped to regions III, IV, and V, respectively (Figure ). The 25 markers were located near the site L.4 and across the IV and V regions (Figure ). The remaining 46 4RLKu-specific markers did not amplify their specific band from any of the five deletion lines and they could not be mapped to any region of the 4RLKu arm (Supplementary Table S2).

Physically localized 4RKu-specific markers using 4RKu dissection lines. (A) Products amplified by primer pair KU.737. (B) Products amplified by primer pair KU-4R.115. (C) Products amplified by primer pair KU-4R.149. (D) Products amplified by primer pair KU-4R.563. (E) Products amplified by primer pair KU-4R.004. (F) Products amplified by primer pair KU.1058. Primer pairs KU.737, KU-4R.115, KU-4R.149, KU-4R.563, KU-4R.004 and KU.1058 respectively represent the markers that mapped to the six regions including I, II, III, IV, V, and the region cross IV and V. M: DNA marker. 4R: MA4RKu. 4RS: MTA4RSKu. 4RL: MTA4RLKu. Arrows indicate the target bands.

Physical map of 338 4RKu-specific markers in six regions on 4RKu chromosome based on PCR amplification using 4RKu dissection lines. The six regions are divided by the five dark lines. In each region, the names of the 4RKu-specific markers are listed on the right. The values on the left indicate numbered sites, which were determined by the FISH signals. The Roman numerals on the left indicates the five regions of the 4RKu chromosome divided by FISH signals. In the schematic diagram, green bands represent Oligo-pSc119.2-1 signals, red bands represent Oligo-pSc200 and Oligo-pSc250 signals, yellow bands represent (AAC)6 signals, and ‘Cen’ represents centromere.

Physical Location of the Powdery Mildew Resistance Gene(s) on 4RL

The resistance of the lines MA4RKu, MTA4RLKu, 16T75-24, 16T177-4, 16T196-22, and parental wheat MY11 and CN27 to powdery mildew was tested in the field. Results showed that MA4RKu, MTA4RLKu and 16T75-24 were highly resistant to powdery mildew (IT = 1), and 16T177-4, 16T196-22, MY11 and CN27 were highly susceptible (IT = 4) (Figure ). Therefore, the powdery mildew resistance gene(s) on 4RLKu was localized to the region V (segment between L.4 and L.8).

Powdery mildew resistance testing. MY11, CN27, line 16T177-4 and line 16T196-22 are highly susceptible to powdery mildew. Lines MA4RKu, MTA4RLKu, 16T75-24 are highly resistant to powdery mildew.

Discussion

4R-Specific Markers

Isozyme, protein, RFLP, DArT, and SSR markers have been used to develop genetic maps of 4R chromosome (Benito et al., 1994; Korzun et al., 1998; Saal and Wricke, 1999; Korzun et al., 2001; Ma et al., 2001; Milczarski et al., 2007; Gustafson et al., 2009; Milczarski et al., 2011; Milczarski et al., 2016). However, rare specific markers can be used to detect 4R chromosome in wheat backgrounds. Three Secale cereale inter-microsatellite (SCIM) markers were located on 4R chromosomes (Camacho et al., 2005). A wheat simple sequence repeat (SSR) marker, Xgwm260, can produce a 4R-specific band (Fu et al., 2010). Two 4R-specific markers were developed using EcoO109I primers (Tomita and Seno, 2012). Six 4R-specific markers derived from expressed sequence tags (ESTs) were obtained (Xu et al., 2012). Li et al. (2013) developed eight 4R-specific markers using the polymerase chain reaction (PCR)-based landmark unique gene (PLUG) system. Qiu et al. (2016) developed 101 4RKu-specific markers using SLAF-seq technology. Almost all of the 4R-specific markers mentioned above are PCR-based markers, and are easier to perform than RFLP, AFLP, and DArT markers. Additionally, these markers not only can be used to distinguish 4R chromatin in wheat backgrounds, but also can be used to construct a map of 4R chromosome. Lukaszewski et al. (2001) attempted to map the genetic position of the RWA resistance locus on 4RL using 4RL arms, which were derived from two different rye lines. However, this attempt was unsuccessful, partially because of the absence of cytologically identifiable recombinants (Lukaszewski et al., 2001). This case indicates the current lack of and the immediate need for additional 4R-specific markers. In this study, 301 new PCR-based and 4RKu-specific markers were developed, adding to the specific markers for distinguishing 4R chromatin in wheat backgrounds and 4R chromosome mapping.

Physical Location of 4RKu-Specific Markers Using Dissection Lines

It has been reported that genetic maps based on recombination rates can not represent the actual physical location of genes and molecular markers on chromosomes (DeScenzo and Wise, 1996). Physical maps complement genetic maps. For wheat and its relatives, chromosomal dissection or deletion lines are useful for determining the physical location of genes and molecular markers (Tsuchida et al., 2008). An array of chromosome deletion stocks have been used to construct physical maps of common wheat (Endo and Gill, 1996). Many molecular markers were physically located on specific segments of barley chromosomes (Ashida et al., 2007; Sakata et al., 2010; Ishihara et al., 2014). Pu et al. (2015) obtained a series of structural aberrations of Thinopyrum bessarabicum chromosome 4J and used them to physically map 101 4J-specific markers and the blue-grained gene Bathb. Agropyron cristatum chromosome 6P-specific STS markers were physically located in 14 regions of this chromosome using 6P deletion lines (Song et al., 2016). Rye chromosome dissection lines were also used to physically map rye-specific markers that are mainly restricted to chromosomes 1R, 2R, and 6R (Kofler et al., 2008; Tsuchida et al., 2008; Gyawali et al., 2009, 2010; Zhuang et al., 2015; Li et al., 2016). In this study, a new FISH map of the 4RKu chromosome was constructed using probes (AAC)6, Oligo-pSc119.2-1, Oligo-pSc200 and Oligo-pSc250. The 4RKu dissection lines combined with the new FISH map of the 4RKu chromosome were used to physically map 400 4RKu-specific markers to five regions (Figure ). Physically locating these markers is beneficial for selecting introgressed 4R small segments in wheat backgrounds to create translocation lines.

Variation of 4RKu-Specific Markers

During the development of rye-specific markers, some rye chromosome-specific markers that were identified using wheat–rye addition lines were often not able to be located to rye chromosome arms or given chromosome segments. For example, some 6R-specific bands could not be amplified from either 6RS or 6RL telosomic addition lines (Xu et al., 2012). Two of the 101 4RKu-specific markers could not be physically mapped to 4RSKu and 4RLKu arms (Qiu et al., 2016). It is possible that this was caused by variations in the structures of rye chromosome arms during the procedure of development of wheat–rye telosomic addition lines (Qiu et al., 2016). In this study, there were 16 4RSKu-specific and 46 4RLKu-specific markers could not be mapped to any regions of 4RKu chromosome using the 4RKu dissection lines. This indicated that the 4RKu chromosome segments in these dissection lines were altered during irradiation.

Localization of Powdery Mildew Resistance Gene(s) on 4RL

It has already been reported that 4R chromosome of rye cultivar German White carried powdery mildew resistance gene (s) (An et al., 2013). Subsequently, it was reported that 4R chromosome of rye Kustro also possesses powdery mildew resistance gene(s), which was mapped on the long arm of 4RKu (4RLKu) (Fu et al., 2014). In this study, the powdery mildew resistance gene(s) on 4RL of Kustro was physically mapped to the V region (segment between L.4 and L.8). Therefore, the powdery mildew resistance gene(s) on 4RLKu has been located to a more explicit segment on 4RLKu arm, and 115 PCR-based and 4RLKu-specific markers that were located to this region has also been developed. These 115 markers will be helpful for localizing the powdery mildew resistance gene(s) on 4RLKu arm in wheat breeding programs. However, more 4R dissection lines are needed to localize the powdery mildew resistance gene(s) to a smaller region on 4RL arm.

Conclusion

A new FISH map of 4RKu chromosome has been constructed. A set of 4RKu dissection lines was obtained and 301 new 4RKu-specific markers were developed. The 301 markers were combined with the 99 4RKu-specific markers developed previously, and were physically mapped to 4RKu chromosome using the new FISH map of 4RKu chromosome and the 4RKu dissection lines. 338 of the 400 markers have been successfully mapped to six regions of 4RKu chromosome. In addition, the powdery mildew resistance gene(s) on 4RLKu arm has been located to a clearly defined region and 115 markers will be helpful for the further localizing the powdery mildew resistance gene(s). The markers developed in this study have enriched the collection of markers that can specifically identify the 4R chromatin in wheat backgrounds and can be used to construct high-density map of 4R chromosome.

Author Contributions

SF and ZT designed the study, analyzed the data and wrote the manuscript. QD, YW, and LQ designed the primers and performed experiments. TR and ZL performed experiments.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.