Skip to main content

Breeding Rice lines for physio-functional food through indica ‘Zhaxima’ × japonica ‘Nanjing 46’ haploid technique

Abstract

Resistant starch (RS) encompasses those forms of starch which are not accessible to human digestive enzymes and are fermented in the colons producing short chain fatty acids. The plant materials containing RS are few in the world. In this contribution, the culture ability of callus from anthers of F1 plants from, landraces, ‘Zhaxima’(Oryza sativa var. indica, high-RS rice line with 7.705 ± 0.142, g/100 g) × ‘Nanjing 46’ (Oryza sativa var. japonica, rice variety with RS content (g/100 g) of 0.200 ± 0.001 crosses were studied for obtaining high RS rice plants. The results showed that when M8 basic induction medium was added with 1.5 mg /L 2,4-D、2 mg /LNAA and 0.3 mg /L KT, the inductivity of callus was high as 32.14% for 21 d after pretreatment at 4 °C for 3 d; When MS differentiation basic medium was added with 2 mg /LKT and 3 mg /L ABA, the frequency of regeneration for callus was 50.3% with only a regeneration frequency of 4.55% grown into green seedlings. The RS content in the seeds was between those of the two parents and was partially normally distributed, the highest RS contents of the regenerated plants was as high as 7.66 ± 1.197%. This produced an efficient technology for regenerating stable rice lines with high RS and good eating quality using anthers culture.

Introduction

Increases in the incidence of type-2 diabetes are being observed throughout the world (Zhou et al. 2016). Starch is a major dietary source of carbohydrate. It is composed of two types of molecules, amylose (Am) and amylopectin (Ap) (Jobling 2004). Based on its enzymatic digestion characteristics, starch can be classified into rapidly digestible starch (RDS), slowly digestible starch (SDS), and resistant starch (RS) (Sajilata et al. 2006). RS is a small fraction of starch that is resistant to hydrolysis by exhaustive α-amylase and pullulanase treatment in vitro (Haralampu 2000). RS encompasses those forms of starch, which are not accessible to human digestive enzymes and are fermented in the colon to produce short chain fatty acids (Ashwar et al. 2015). Consumption of foods high in resistant starch (RS) can help to control type-2 diabetes, because its slow digestion and absorption by the small intestine decrease postprandial glucose and insulin responses (Raigond et al. 2015). In addition, RS has similar physiological effects to dietary fiber, affecting body weight and energy balance and increasing lipid excretion to reduce calorie intake and decrease the serum lipid levels (Jacobasch et al. 1999). These unique physical functions of RS have received increasing attention from plant researchers in recent years (Yang et al. 2012; Zhou et al. 2016). Thus, improvement of the amounts and properties of RS in foods is an important goal (Zhou et al. 2016).

As the primary dietary source of carbohydrates in the world, rice (Oryza sativa L.) plays an important role among cereal crops in meeting energy requirements and nutrient intake. However, the RS content in ordinary rice variety is low, generally about 1%, which is not enough to confer the associated health benefits (Frei et al. 2003). In addition, rice with high content of RS also has inferior eating quality, resulting in lower market prices (Raigond et al. 2015). Therefore, many studies have focused on elevating the RS content in rice cultivars via mutation breeding and bioengineering (Hu et al. 2004; Shen et al. 2006; Kubo et al. 2008; Wani et al. 2012). Higher amylose content in starch is generally suggestive of increased RS content, thus producing increased levels of amylose through breeding is a possible method to increase endogenous RS content. Two methods have been used with success: inbreeding of mutants containing genes for high amylose production and the inhibition of starch-branching enzyme (SBE) (Deupuis et al. 2014). Many rice mutants with elevated RS content have been identified, including RS111 (Yang et al. 2006) and ‘Jiangtangdao 1’ (Zhu et al. 2010). ‘Teqing Resistant Starch’ (TRS) is another high-amylose and high RS transgenic line developed by modifying antisense RNA inhibition for SBE in rice (Wei et al. 2010). A putative gene, sbe3-rs of RS was also identified and the codominant cleaved amplified polymorphic sequence (CAPS) marker could be used in marker-assisted breeding to develop rice cultivars with elevated RS, which is otherwise difficult to accurately assess in crop (Yang et al. 2012). Although much progress has been made in obtaining the rice cultivars high in RS, three problems remain to be solved. Firstly, only few rice germplasms high in RS have been reported. Secondly, many rice varieties high in RS belong to the indica subspecies, with low yield, and are very difficult to cross rapidly with japonica rice plants by conventional breeding as this process of breeding needs quite a long time of 5–10 years. Thirdly, the safety of transgenic rice varieties high in RS need to be verified in longer time. Currently, safe, rapid and efficient techniques for breeding rice varieties with both good eating quality and high RS content are lacking.

Doubled haploids (DHs) in plants have complete homozygosity and can be achieved in one generation from hybrid plants. DH production includes two major steps: haploid induction and chromosome doubling. Chromosome doubling of haploid plants has been routinely and successfully performed using colchicine. However, the success and efficiency of haploid induction varies among crop species (Niu et al. 2014). In rice, many theoretical and applied problems remain to be solved, such as the low fertility of indica–japonica F1 hybrids, the lower rate of plant regeneration and the lower seed setting rate of regenerated plants. The factors limit the application of the rice breeding technology by anther culture (Zhao et al. 2012).

In this study, a safe and highly efficient technique for generating new rice lines from the anthers of F1 hybrids of ‘Zhaxima’ × ‘Nanjing 46’ was introduced to obtain rice lines both high in RS and excellent in eating quality. Some genetic characteristics of the regenerated rice plants were also studied to elucidate these traits in the parent indica and japonica plants for tissue culture ability. The genetic characteristics of the partial-japonica type hybrid were helpful to obtain more regeneration of green fertile plants, and the RS content in the seeds of these regeneration plants was all higher than in the male parent with lower content of RS.

Materials and methods

Materials

The indica variety ‘Zhaxima’ from Yunnan was the female parent (stamens removed), and it has high RS content in the seeds and is a landraces (Lin et al. 2013). The japonica variety ‘Nanjing 46’ was male parent, and is known for its good eating quality (Wang et al. 2009). The F1 cross ‘Zhaxima’ × ‘Nanjing 46’ was designed to combine the RS quality of the former with the eating quality of the latter and was produced in the summer of 2012 in Nanjing (118°46′E, 32°24′N 30 m above sea level), China.

Anther culture of the F1 population was carried out as follows. ‘Zhaxima’, ‘Nanjing 46’, and the F1 hybrids were planted in the experimental field of Jiangsu Academy of Agricultural Sciences, Hainan, China, in the spring of 2013. The anthers of the F1 population were collected for anther culturing in the laboratory. Then, anther-culture regenerants were planted in the experimental field of Jiangsu Academy of Agricultural Sciences, Nanjing, China in the summer of 2013.

Regeneration plant callus induction from anthers of F1 hybrids

The media used in the different stages of callus induction were summarized in Table 1. Anthers were cold pretreated in the dark at 4 °C for 1–10 d on induction medium (IM) based on M8 basal medium (Mei et al. 1988). Seven differentiation media (DM) supplemented with different hormone compositions were used in this study (Table 1). Individual dishes were labelled with the anther type and location in the inflorescence (top, middle, base), sealed with Parafilm, and incubated in the dark at 26–28 °C for callus induction. Four different supplemented MS media (Murashige & Skoog 1962) were used for differentiation (Table 1). Calli, at least 2 mm in size from anthers, were transferred to MS basal semisolid medium supplemented with different plant hormone compositions and incubated at 28 °C under a 16-h light photoperiod supplied by cool white fluorescent lamps (66 μmol m− 2 s− 1) for plant regeneration. Tissues producing green buds were transferred to rooting medium (RM) and cultured at 26 ± 2 °C under light for rooting. The plantlets were cultured at 26–28 °C with a 16-h photoperiod at a photosynthetic photon flux density of 200 μmol m− 2 s− 1. Before being transferred to the experimental field, the plantlets were hardened in their tubes by adding sterile water to drown their roots and cultured for 3 d in a culture chamber. When the height of the plants reached 5–6 cm, they were transferred to clay pots (pot with 0.1 m diameter and 0.6 m high) filled with soil, fertilized with slow-release fertilizer, and kept in a growth chamber at 28–30 °C with a 16/8 h (day/night) photoperiod, and finally planted in the field.

Table 1 Plant hormone components in the callus induction, differentiation, and rooting mediaa

Measurements of stomatal guard cell size of the leaf epidermis and plant height of regeneration plants

To detect the ploidy level of regeneration plants, the size of the guard cell hypodermis on the flag leaf was measured at the rice jointing stage (Liang 1979). Clear nail polish was painted on the leaf epidermis, allowed to dry, then pulled off. Subsequently, the stomatal guard cells of epidermis were observed under a microscope. The perimeters of 100 randomly selected guard cells on each leaf blade were measured. Ten blades were observed for each regeneration plant. Plant height was directly measured at the heading stage using a ruler.

Molecular analysis of the regeneration plant subspecies type by insertion–deletion (InDel) markers

‘Nipponbare’ was the japonica control, and ‘Nanjing 16’ was the indica control. Ten seeds of each regeneration plant were randomly selected and germinated at 37 °C. Then, seedlings were grown in an incubator at 28 °C until the leaves reached the heart embryo, one of embryo developmental stages on shape of callus, when 0.1 g of fresh leaves were sampled. The regeneration materials were sampled at the tillering stage. DNA was extracted from all samples by the CTAB method (Murray & Thompson 1980), and the polymorphisms were used to identify the japonica and indica types based on modified protocols of Shen et al. (2003) and Lu et al. (2009). Insertion–deletion (InDel) markers are codominant. Using PCR and electrophoresis of the DNA of indica, ‘Nanjing 16’, and japonica, ‘Nipponbare’ (as the standard), the genotype of each tested cultivated and wild rice sample was determined based on the InDels (Additional file 1: Table S1). Samples consistent with the ‘Nanjing 16’ bands were recorded as the homozygous indica genotype (II), samples matching the banding pattern of ‘Nipponbare’ samples were identified as the homozygous japonica genotype (JJ), and those with both ‘Nanjing 16’ and ‘Nipponbare’ bands were considered to be indica and japonica heterozygous genotype (IJ). The mean gene frequency (F) of all InDel sites was calculated using the following formulas (Shen et al. 2003; Lu et al. 2009),

japonica gene frequency \( Fi=\frac{2{\sum}_1^N Xii+{\sum}_1^N Xij}{2N} \)

indica gene frequency \( Fj=\frac{2{\sum}_1^N Xii+{\sum}_1^N Xij}{2N} \) where Xii is a specific indel locus whose position is the same as a ‘Nanjing 16’ band and denotes an indica homozygous genotype (ii); Xjj is a specific indel site whose position is the same as a ‘Nipponbare’ band and denotes a japonica homozygous genotype (jj); Xij is a specific indel site whose location matches that of both ‘Nanjing 16’ and ‘Nipponbare’ and indicates a heterozygous indica and japonica genotype (ij) and N is the number of indel primer pairs (sites) included (Additional file 2 Table S2).

Wide compatibility variety (WCV) gene sequences have been studied and their functional motifs are known (Ikehashi & Araki 1984; Ji et al. 2005; Ikehashi & Araki 1986; Morinaga & Kuriyama 1958; Qiu et al. 2005;Yanagihara et al. 1995). The S5 WCV gene has been sequenced in indica (accession number EU889295; S5-i), japonica (EU889294; S5-j), and japonica 02428 (EU889293; S5-n). These sequences differ in a 69-bp deletion upstream and a 67-bp deletion downstream of the ATG at the translation initiation point; these deletions lead to loss of function and prevent S5-i and S5-j interaction. According to the deletion sequence of both sides of the gene (Yang et al. 2009), the primers S5136-Forward (5′-ATCAACCCATTTCCTTTCCT-3′) and S5136-Reverse (5′-ATACGCTCGATCGGATTAAC-3′) were designed.

Measurements of RS content in the regeneration plants

RS was measured using the Megazyme RS assay kit (Megazyme, Co. Wicklow, Ireland), which has been widely employed for RS determination in crops (McCleary et al. 2002). The grain sample was treated with 10 mg/mL pancreatic α-amylase and 3 U/mL amyloglucosidase (AMG) enzymes for hydrolysis and solubilization of non-resistant starch. After the enzymatic reaction was terminated by adding 99% ethanol, RS was recovered as a pellet by centrifugation (approx. 3000 g, 10 min). RS in the pellet was dissolved in 2 mol L− 1 KOH before being added into the reaction solution and was repeatedly washed and decanted. Then, starch in the solution was quantitatively hydrolyzed to glucose with AMG. D-glucose was measured with glucose oxidase/peroxidase (GOPOD) at 510 nm wavelength against the reagent blank. All analyses were repeated three times for error control.

Data analysis

Data parameters were calculated as follows:

  • Callus induction rate = (number of calluses) / (number of pollen grains inoculated) × 100%.

  • Differentiation rate = (number of plantlets obtained from callus) / (number of calluses inoculated) × 100%.

  • Plant regeneration rate = (number of plantlets) / (number of inoculated calluses) × 100%.

  • Resistant Starch (g/100 g sample) = ΔE × F × 100/0.1 × 1/1000 × 100/W × 162/180 = ΔE × F/W × 9.27,

  • Where ΔE = absorbance (reaction) read against the reagent.

blank; F = conversion from absorbance to micrograms =100 (mg glucose)/absorbance of 100 mg glucose; 100/0.1 = volume correction (0.1 mL taken from 100 mL); 1/1000 = conversion from micrograms to milligrams; W = dry weight of sample analyzed [= “as is” weight ´ (100-moisture content)/100]; 100/W = factor to present starch as a percentage of sample weight; 162/180 = factor to convert from free glucose, as determined, to anhydroglucose as occurs in starch; 10.3/0.1 = volume correction (0.1 mL taken from 10.3 mL) for samples containing 0–10% RS where the incubation solution is not diluted and the final volume is 10.3 mL (McCleary et al. 2002). One-way analysis of variance of all data was performed using SPSS 17.0 (IBM, Chicago, IL, USA).

Results

Establishment of anther regeneration technology from the hybrid cross “Zhaxima” and “Nanjing 46”

After pretreatment at 4 °C for 0–10 d, the anthers were cultured on M8 induction medium with different hormones. The anthers were induced to form callus during 21 d in the dark. The frequency of callus induction was as high as 32.14% on the IM-3 medium after pretreatment at 4 °C for 3 d (Table 2). The results also showed that anthers from both the basal and middle parts of the panicle had higher callus induction frequency than those from top part (Table 2). The callus was then transferred to DM medium, where it began to form green tissues under light treatment over 14 D. callus with green tissue was cultured for about 30 d on DM-2, with a greening rate of about 4.55% (Fig. 1). The green plantlets could take root on RM for about 20 d. The regeneration plants were placed in the experimental field during the clover heart stage. Notably, the regeneration plants also had a very high rate of albinism. We obtained the seeds of rice plants from anthor culture of F1 from, landraces, ‘Zhaxima’(Oryza sativa var. indica,) × ‘Nanjing 46’ (Oryza sativa var. japonica) crosses that year. In this contribution, the plant regeneration rate from anther culture was as high as 4.55%, but the albino rate of regeneration plants was as high as 47.5%. We further analyzed the morphological and molecular characteristics of the regeneration plants got in this work, in order to provide the basis for improvement of anther culture in rice.

Table 2 Effects of different hormone components and treatment times on callus induction
Fig. 1
figure 1

Effect of phytohormone components on callus differentiation rate. MS was the basic medium (Murashige & Skoog 1962): sucrose 30 g•L−1 and agar 8 g•L−1. DM-1: 2 mg•L−1 KT and 2 mg•L−1 ABA; DM-2: 2 mg•L− 1 KT and 3 mg•L− 1 ABA; DM-3: 0.2 mg•L− 1 2,4-D, 2.5 mg•L− 1 KT, and 2 mg•L− 1 ABA; DM-4: 0.5 mg•L− 1 2,4-D, 2 mg•L− 1 KT, and 2 mg•L− 1 ABA. 2,4-D, 2,4-dichlorophenoxyacetic acid; NAA, 1-naphthylacetic acid; KT, kinetin; ABA, abscisic acid; DM: differentiation medium. Different letters above the bars indicate significant differences (P < 0.05)

Analysis of indica and japonica types of the high regenerated plants

According to 35 InDel sites (Fig. 2), the indica and japonica types of samples (ii, jj and ij) were calculated. Table 3 shows the standard characteristics of indica and japonica types based on the frequency calculations. The japonica genotype frequencies of ‘Zhaxima’ and ‘Nanjing 46’ were 0 and 1, respectively, demonstrated that these parents were typical indica and japonica types, respectively. Type 1 plants had frequencies of 0.37 indica and 0.63 japonica, and were a partial-japonica type. In contrast, type 2 plants had frequencies of 0.63 indica and 0.37 japonica, indicated a partial-indica type. The frequencies of type 3 were 0.41 indica and 0.59 japonica, an intermediate type. We also analyzed the yellow and albino regeneration plants in this study. The indica frequency of yellow plants was 0.44 and their japonica frequency was 0.56, while the frequencies of the albino ones were 0.53 and 0.47, respectively; Both lines were of the intermediate type. These results indicated that japonica plants were better suited to anther culture than indica ones. WCVs with S5-n are able to produce highly fertile hybrids when crossed with both indica and japonica varieties. Using the S5–136 primers detected, the two parents together with their regenerated plants did not contain the wide compatibility gene, S5. We also observed that the seed-setting rate of the F1 generation was very low, which indicated that the genetic obstacles between the indica and japonica without WCVs might be one of the reasons for the low fertility. As the explant, anther from the hybrid with indica japonica incompatibility was one of the main reasons for low regeneration rate. A noteworthy phenomenon was that the higher culture ability of hybrid rice cross indica and japonica seems partial male such as japonica, Nanjing 46 while the equal distribution of indica and japonica characteristics was not suitable for regenerating green plants and setting seeds.

Fig. 2
figure 2

Amplification results of ‘Nanjing 46’, ‘Zhaxima’, F1 hybrids, and regeneration plants using indel primer (R3M23, R4M17, R4M43, and R5M30). In row “a”, lane 1 is Nipponbare, lane 2 is Nanjing 16, lane 3 is Zhaxima, and lane 4 is Nanjing 46. In row “b”, lane 1 is Nipponbare, lane 2 is Nanjing 16, lane 3 is the F1 hybrid, and lane 4 is the regenerated plants

Table 3 Indica (Fi) and japonica (Fj) gene frequencies in rice samplesa

Plant height and seed-set of the regenerants of the green plants

Based on plant height, the 40 regeneration lines could be divided into three types (Table 4, Fig. 3). Type 1 plants grew to 72.1 ± 4.5 cm in plant height, similar to the male parent ‘Nanjing 46’. The stomatal perimeter on the lower epidermis of type 1 plants (59.58 ± 2.20 μm) was larger than those of both parents (‘Nanjing46’: 55.97 ± 1.24 μm; ‘Zhaxima’: 45.98 ± 0.54 μm), and all of those plants set seed. Type 2 plants were 128.8 ± 2.7 cm in plant height and did not set seeds with abnormally hooked anthers and also had a larger stomatal perimeter on the lower epidermis (63.47 ± 3.40 μm). Type 3 plants were the shortest at 58.9 ± 1.4 cm in plant height and also had the smallest lower-epidermal stomatal perimeter (51.9 ± 2.4 μm) of the three type plants. Furthermore,the rice lines in type 3 had significantly smaller spikelet than those of the other two regenerated plants and also did not produce seeds. According to Liang (1979) and Choe et al. (2012), the shortest rice types may be haploid, which are significantly less tall than the high and the intermediate types. Low rate of natural doubling in the regeneration plants was also the main reason for low regeneration rate in this work.

Table 4 External peripheral long of the stoma, stem length and the seed rate in regeneration plants from anthers
Fig. 3
figure 3

Guard cells of ‘Zhaxima’, ‘Nanjing 46’, F1 hybrids, and regenerated plants. a female parent; b male parent; c F1 hybrid; d type 1 plants (height: 72.1 ± 4.5 cm); e type 2 plants (height: 128.8 ± 2.7 cm); f type 3 plants (height: 58.9 ± 1.4 cm)

Resistant starch (RS) content of regenerated rice lines

We obtained seeds from regeneration plants of 12 rice lines and determined their RS contents in their seeds (Table 5). ‘Zhaxima’ and ‘Nanjing 46’ contained 7.705 ± 0.142% and 0.200 ± 0.000% RS, respectively. The RS contents of the regenerated plants lay between those of the two parents, exhibiting the partial normal distribution. The RS contents of the regeneration plants were all higher than that of male parent, Nanjing 46. The rice plants with RS content between 0 and 3% accounted for 83.37% of the regenerated plants, while 16.67% had RS of 3–8%. Furthermore, the amylose starch contents in these regenerated plants were also reduced. Although the regeneration rate was low, the trait of the fertile regeneration plant with higher content of resistant starch has the obvious super mother advantage with good eating quality, which is easier to improve.

Table 5 Amylose starch content and resistant starch (RS) content of regenerated rice plants a

Discussion

Resistant starch (RS) has the potential to protect against diabetes and reduce the incidence of diarrhea, inflammatory bowel disease, colon cancer, and chronic renal and hepatic diseases (Zhou et al. 2016). The cultivation of rice varieties with unique medical value is not only a focal point in breeding but also has particularly received widespread attention among nutritionists. Resistant starch has novel functions similar to dietary fiber that can regulate metabolism. However, breeding rice varieties high in RS by time-consuming conventional methods is inefficient and difficult (Hu et al. 2004), because the RS-rich materials more belong to landraces and their yields are lower. The current focus of breeders is mainly how to obtain the crop high in RS. Some high RS rice varieties, such as ‘RS111’, ‘Zhefu201’, ‘AMF18’, Goami No. 2, Gongmi No. 3, Jiangtangdao 1 and ‘Yitang1’, have been bred using chemical mutagens or aerospace radiation mutagenesis or conventional hybridization (Shen et al. 2006; Kubo et al. 2008; Bai et al. 2012; Lee et al. 2006; Matsumoto et al. 2012). Although these methods have proven successful in crop breeding, it is highly random, not very efficient and also time-consuming (Rahman et al. 2007). With the rapid development of molecular biology, scientist have also tried to breed high RS rice by transgenic techniques (Deupuis et al. 2014) by inhibiting the activity of the SBEs. SBEs are one of the four major enzyme classes involved in starch biosynthesis in plants, and their activities play a crucial role in determining the structure and physical properties of starch granules (Tetlow & Emes 2014). Although inhibiting the activity of the SBE can get higher RS rice lines, modification of SBEs in planta also influences the degradation of starch reserves in developing seeds, thus impacting seedling vigor, this also in turn deteriorates the good eating quality of the rice grain (Nakamura et al. 2010; Sawada et al. 2009; Xia et al. 2011; Nakamura et al. 2012). Zhou’s discovery provides an opportunity to increase RS content of cooked rice, especially in the indica varieties, which predominates in southern Asia (Zhou et al. 2016). However, because of security concerns about transgenic plants and, they are difficult to grow on large tracts of land. In this study, we produced an efficient technology for regenerating stable rice lines with high RS and good eating quality using anthers culture.

Androgenesis (anther or microspore culture) is one of two basic approaches for developing haploids in higher plants. In vitro process, microspore cells with haploid genomes develop into embryo-like structures on culture medium. The embryo-like structures further develop into haploid plantlets (Jauhar et al. 2009). With the rapid development of modern biotechnology, research on androgenesis and haploid breeding in Gramineae is gradually deepening, and some progress has been made in obtaining haploid grasses via anther technology. Although haploids have been successfully developed via in vitro culture of unfertilized ovules and ovaries in more than 20 angiosperm species since 1976 (Wu et al. 2004), anther or microspore culture in crops has experienced obstacles in haploid production, such as high rates of albinism, low response rates of some genotypes, and long periods for the inducing and regenerating processes, which could cause detrimental gametoclonal variation and mixed-ploidy plants (Niu et al. 2014). Cold pretreatment and using anther collected from the top or middle spike position can improve the ability of anther culture (Table 2). There may have been a synergistic effect between 2,4-D concentration and the cold pretreatment, because cold pretreatment is linked to the embryogenic capacity of plant tissues to acquire a specific hormonal status. Tian et al. (2015) reported that pretreatments of low temperature on rice (Oryza sativa L.) anthers changed polysaccharide and protein composition of the anther walls and increased pollen fertility and then callus induction. Furthermore, the molecular characteristics of the two parents are the most important determinants of anther regeneration efficiency. In this study, we selected the female parent, indica ‘Zhaxima’ with high RS content and the male parent, japonica Nanjing 46′ with good eating quality of grain for anther cultural. Our results showed that the higher culture ability of hybrid rice cross indica and japonica seemed partial male such as japonica, Nanjing 46, while the equal distribution of indica and japonica characteristics was not suitable for regenerating green plants and setting seeds.

The ploidy level of microspore-derived regenerant varies among cereals. Microspores when cultured in vitro can spontaneously double the gametophytic chromosome number, but, for unknown reasons, this does not occur in all cases (Jähne & Lörz 1995). In wheat anther culture of a range of genotypes, 20–50% of green regenerants were doubled (Henry & Buyser 1999), while in barley, up to 87% of plants were spontaneous dihaploids (Hoekstra et al. 1993), and in rice up to 72% were reported (Cho & Zapata 1988). Ploidy identification in plant by pressing young root to count chromosome number in cell was a classical method. For this conventional method of ploidy identification of plants using root have certain limitations, because the rice plants from the rice regeneration must be conducted by the grain for the identification of resistant starch. Getting the green plants with seeds is necessary for rice breeding with high content of resistant starch. If the roots of plants were destroyed, the growth of the plants might be affected, conferring the lower seed-setting of rice varieties. We adopted the stomatal perimeter method for identifying ploidy. In this work, based on the observation of stomatal perimeter of the rice plants, combined with other characteristics such as plant height and seed setting rate, we also clearly distinguish the ploidy of the regeneration plants, and this method did not destroy the plants. Our 40 regeneration lines were divided into three types based on plant height, stomatal perimeter on the lower epidermis and setting seed characteristics. Type 2 plants were similar to the female parent ‘Zhaxima’ in plant height and stomatal perimeter on the lower epidermis, but sterile with abnormally hooked anthers. Previous studies have suggested that plants with small stomatal perimeters might be haploid (Liang 1979; Choe et al. 2012), which is similar to our results. The acarpous plants belonged to the type 3 with the smallest stomatal size might be caused by haploidy. Only 30.0% of the plants were dihaploid and produced seeds. There were still 22.5% regeneration plants that were not doubled. In addition to these sterile plants, there were 47.5% sterile plants with similar characteristics of plant height and stomatal sizes as their parents. The results indicated that the reason on low rate of producing RS rice plants from anthers of crosses between indica and japonica might be other intrinsic factors, in addition to the undouble obstacle.

The current focus of breeders is mainly how to obtain the crop high in RS. In this work, we used an efficient and safe regeneration system using F1 hybrids of ‘Zhaxima’ × ‘Nanjing 46’ to get the rice plants with high RS content (7.6 g/100 g dry weight). Based on genetic differentiation of the indica and japonica types using 35 indel sites (Lu et al. 2009), ‘Zhaxima’ was a typical indica variety and ‘Nanjing 46’ was a typical japonica plant. Using the S5–136 primers, we determined that two parents did not contain the wide compatibility gene (S5). Regeneration plants of hybrids of indica and japonica may exhibit morphological traits that differ from those of their parents as a result of random recombination and assortment of chromosomes as well as genetic isolation. Among them, the regeneration plants of type 1 classified in this work can be easy to obtain the fertile green plant, and their subspecies characteristics were partial japonica, like their male parent, Nanjing 46.

Regeneration plants of type 1 all produced seeds, and the RS content of these plants were normally distributed. According to the results of Sun et al. (2012) on the genetic traits of hybrid rice seeds, different combinations of parents produce offspring with different RS contents. The variance analysis of RS content indicated that RS quality might be regulated by different genes and RS level is mainly controlled by additive effects in the hybrid offspring. Thus, the RS contents of offspring are generally distributed between the levels of the parents, also in this experiment they were skewed toward the male parent, just like our results here. The authors did present the amylose content, which is closely related to RS content (Zeng et al. 2016). The high amylose content in rice grain is positively correlated with RS content (Hu et al. 2004). Our results suggest that the high RS in some regenerated plants might be due to Wx-a allele of Wx gene coding for the enzyme for amylose synthesis, which is incorporated from the female indica parent. Through map-based cloning of a RS locus in indica rice, Zhou et al. (2016) have identified a defective soluble starch synthase gene (SSIIIa) responsible for RS production and further showed that RS production is dependent on the high expression of the Waxya (Wxa) allele, which is prevalent in indica varieties. Although the same ssIIIa mutation could be used in japonica rice together with introduction of a Wxa gene, the resulting rice would have higher amylose content than what is normally preferred by consumers of japonica varieties (Zhou et al. 2016). In this study, Nanjing 46 with low amylose and good eating taste was used as male parent to increase RS in japonica rice lines by anther culture, adding a new way of using indica rice to improve japonia. The different alleles of the SSIIa gene are responsible for differences in amylopectin structure between the indica and japonica rice varieties (Umemoto et al. 2002). This suggests that the interactions with different soluble starch synthase genes (SSIIIa and SSIIa) in japonica rice variety (Nanjing 46) and Wxa in indica rice (Zaxima) can be helpful to obtain high RS content and improve the taste quality of rice. Further analysis of the regenerative plants with different partial-japonica or partial-indica genetic background with high RS may provide new clues for improving RS in rice. It is believed that the strategies to increase RS in rice, including indica and japonica, will be developed in the future with the discovery of the molecular basis underlying RS production in rice.

Conclusions

In this study, we established an efficient technology for regenerating stable rice lines high in RS using anthers, which not only shorten the breeding period and improve breeding efficiency, but also avoids having to identify multiple generations. The higher culture ability of hybrid rice cross indica and japonica seemed partial male such as japonica, Nanjing 46, while the equal distribution of indica and japonica characteristics was not suitable for regenerating green plants and setting seeds. Furthermore, the rice materials carrying the genetic components of RS by anther culture will be the basis for further studies on the genetic basis of the control of resistant starch in rice.

Availability of data and materials

All the data and materials were kept in Dr. Li Xia’s lab, Institute of Food Crops, Jiangsu Academy of Agricultural Sciences, Jiangsu High Quality Rice R&D Center, Nanjing Branch, China National Center for Rice Improvement, Nanjing 210014, P.R. China. Data sharing is not applicable to this article as no data sets were analyzed during the current study. If the readers wish to understand these data and material in detail, please contact the corresponding author for data requests.

Abbreviations

2,4-D:

2,4-dichlorophenoxyacetic acid

6-BA:

6-benzylaminopurine

ABA:

Abscisic acid

Am:

Amylose

AMG:

Amyloglucosidase

Ap:

Amylopectin

CAPS:

Codominant cleaved amplified polymorphic sequence

DH:

Doubled haploid

DM:

Differentiation medium

GABA:

r-aminobutyric acid

GI:

Glycemic index

GOPO:

Glucose oxidase/peroxidase

IM:

Induction medium

InDel:

Insertion–deletion

KT:

Kinetin

NAA:

1-naphthaleneacetic acid

RDS:

Rapidly digestible starch

RM:

Rooting medium

RS:

Resistant starch

SBE:

Starch branching enzymes

SDS:

Slowly digestible starch

SM:

Seedling medium

TRS:

Teqing Resistant Starch

WCVs:

Wide compatibility varieties

References

  • Ashwar, B. A., Gani, A., Shah, A., Wani, I. A., & Masoodi, F. A. (2015). Preparation, health benefits and applications of resistant starch-a review. Starch/Stärke, 67, 1–15.

    Article  CAS  Google Scholar 

  • Bai, J. J., Zhu, H. M., Li, D. L., & Piao, Z. Z. (2012). Effect of the high resistant starch rice on blood glucose and lipid metabolism in diabetic GK mice. Chinese Institute of Food Science and Technology, 12, 16–20 (in Chinese with English Abstract).

    CAS  Google Scholar 

  • Cho, M. S., & Zapata, F. J. (1988). Callus formation and plant regeneration in isolated pollen culture of rice (Oryza sativa L. cv. Taipei 309). Plant Science, 58, 239–244.

    Article  Google Scholar 

  • Choe, E., Carbonero, C. H., Mulvaney, K., Rayburn, A. L., & Mumm, R. H. (2012). Improving in vivo maize doubled haploid production efficiency through early detection of false positives. Plant Breeding, 131, 399–401.

    Article  Google Scholar 

  • Chu, C. Y., Sun, C. S., Chu, C. C., Wang, C. C., Pi, F. Y., Yiu, K. C., & Hsu, C. (1975). Establishment of an efficient medium for anther culture of rice through comparative experiments on the nitrogen sources. Sci. Sin. (Peking), 18, 659–668.

    Google Scholar 

  • Deupuis, J. H., Liu, Q., & Yada, R. Y. (2014). Methodologies for increasing the resistant starch content of food starches: A review. Compr Rev Food Sci F, 13, 1219–1234.

    Article  CAS  Google Scholar 

  • Frei, M., Siddhuraju, P., & Becker, K. (2003). Studies on the in vitro starch digestibility and the glycemic index of six different indigenous rice cultivars from the Philippines. Food Chemistry, 83, 395–402.

    Article  CAS  Google Scholar 

  • Haralampu, S. G. (2000). Resistant starch-a review of the physical properties and biological impact of RS3. Carbohyd Poly, 41, 285–292.

    Article  CAS  Google Scholar 

  • Henry, Y., & Buyser, J. D. (1999). Wheat anther culture: Agronomic performance of doubled haploid lines and the release of a new variety ‘florin’. In Y. P. S. Bajaj (Ed.), Biotechnology in agriculture and forestry: Wheat (pp. 285–352). Switzerland: Springer Verlag.

  • Hoekstra, S., van Zijderveld, M. H., Heidekamp, F., & van der Mark, F. (1993). Microspore culture of Hordeum vulgare L.: The influence of density and osmolality. Plant Cell Reports, 12, 661–665.

    Article  CAS  PubMed  Google Scholar 

  • Hu, P. S., Zhao, H. J., Duan, Z. Y., Zhang, Z. L., & Wu, D. X. (2004). Starch digestibility and the estimated glycemic score of different types of rice differing in amylose contents. Journal of Cereal Science, 40, 231–237.

    Article  CAS  Google Scholar 

  • Ikehashi, H., & Araki, H. (1984). Varietal screening for compatibility types revealed in F1 fertility of distant crosses in rice. Jpn J Breed, 34, 304–312.

    Article  Google Scholar 

  • Ikehashi, H., & Araki, H. (1986). Genetics of F1 sterility in remote crosses in rice. International Rice Research Institute. In Rice genetics (pp. 119–130).

    Google Scholar 

  • Jacobasch, G., Schmiedl, D., Kruschewski, M., & Schmehl, K. (1999). Dietary resistant starch and chronic inflammatory bowel diseases. International Journal of Colorectal Disease, 14, 201–211.

    Article  CAS  PubMed  Google Scholar 

  • Jähne, A., & Lörz, H. (1995). Cereal microspore culture. Plant Science, 109, 1–12.

    Article  Google Scholar 

  • Jauhar, P. P., Xu, S. S., & Baenziger, P. S. (2009). Haploidy in cultivated wheats: Induction and utility in basic and applied research. Crop Science, 49, 737–755.

    Article  Google Scholar 

  • Ji, Q., Lu, J., Chao, Q., Gu, M., & Xu, M. (2005). Delimiting a rice wide-compatibility gene S5n to a 50 kb region. Theoretical and Applied Genetics, 111, 1495–1503.

    Article  CAS  PubMed  Google Scholar 

  • Jobling, S. (2004). Improving starch for food and industrial applications. Curr Opin Plant Bio, 7, 210–218.

    Article  CAS  Google Scholar 

  • Kubo, A., Yuguchi, Y., Takemasa, M., Suzuki, S., Satoh, H., & Kitamur, S. (2008). The use of micro-beam X-ray diffraction for the characterization of starch crystal structure in rice mutant kernels of waxy, amylose extender and sugary1. Journal of Cereal Science, 48, 92–97.

    Article  CAS  Google Scholar 

  • Lee, K. W., Song, K. E., Lee, H. S., Kim, Y. K., Lee, S. W., Kim, D. J., Hwang, W. S., Choe, S. J., Kim, Y. S., & Kim, T. Y. (2006). The effects of Goami no. 2 rice, a natural fiber-rich rice, on body weight and lipid metabolism. Obesity (Silver Spring), 14(3), 423–430.

    Article  CAS  Google Scholar 

  • Liang, H. X. (1979). A preliminary study on the identification for the ploidy of the wheat pollen plantlets by the length of stoma. Acta Botanica Yunnan, 1, 138–142 (In Chinese with English Abstract).

    Google Scholar 

  • Lin, J., Sun, B. X., Fang, X. W., Wang, Y. P., Zhang, S. B., & Wang, Y. J. (2013). Starch properties of rice enriched with resistant starch. Acta Agriculturae Boreali-sinica, 28, 58–61 (In Chinese with English Abstract).

    Google Scholar 

  • Lu, B. R., Cai, X. X., & Xin, J. (2009). Efficient indica and japonica rice identification based on the InDel molecular method: Its implication in rice breeding andevolutionary research. Progress in Natural Science, 19, 1241–1252 (In Chinese with English Abstract).

    Article  CAS  Google Scholar 

  • Matsumoto, K., Maekawa, M., Nakaya, M., Takemitsu, H., Satoh, H., & Kitamura, S. (2012). Wx/ae double mutant brown rice prevents the rise in plasma lipid and glucose levels in mice. Bioscience, Biotechnology and Biochemistry, 76, 2112–2117.

    Article  CAS  Google Scholar 

  • McCleary, B. V., McNally, M., & Rossiter, P. (2002). Measurement of resistant starch by enzymatic digestion in starch and selected plant materials: Collaborative study. Journal of AOAC International, 85, 1103–1111.

    CAS  PubMed  Google Scholar 

  • Mei, C. S., Zhang, J. Y., & Wu, G. N. (1988). Improving regeneration rate of anther culture in indica rice (Oryza sativa L. subsp. Indica). Jiangsu Journal of Agricultural. Sciences, 4, 45–48 (In Chinese with English Abstract).

    Google Scholar 

  • Morinaga, T., & Kuriyama, H. (1958). Intermediate type of rice in the subcontinent of India and Java. Jpn J., 7, 253–269.

    Google Scholar 

  • Murashige, T., & Skoog, F. (1962). A revised medium for rapid growth and bioassays with tobacco tissue culture. Plant Physiology, 15, 473–479.

    Article  CAS  Google Scholar 

  • Murray, M. G., & Thompson, W. F. (1980). Rapid isolation of high molecular weight plant DNA. Nucleic Acids Research, 8, 432l–4325l.

    Article  Google Scholar 

  • Nakamura, Y., Ono, M., Utsumi, C., & Steup, M. (2012). Functional interaction between plastidial starch phosphorylase and starch branching enzymes from rice during the synthesis of branched maltodextrins. Plant & Cell Physiology, 53, 869–878.

    Article  CAS  Google Scholar 

  • Nakamura, Y., Utsumi, Y., Sawada, T., Aihara, S., Utsumi, C., Yoshida, M., & Kitamura, S. (2010). Characterization of the reactions of starch branching enzymes from rice endosperm. Plant & Cell Physiology, 51, 776–794.

    Article  CAS  Google Scholar 

  • Niu, Z. X., Jiang, A. X., Hammad, W. B., Oladzadabbasabadi, A., Xu, S. S., Mergoum, M., & Elias, E. M. (2014). Review of doubled haploid production in durum and common wheat through wheat × maize hybridization. Plant Breeding, 133, 313–320.

    Article  CAS  Google Scholar 

  • Qiu, S. Q., Liu, K. D., Jiang, J. X., Song, X., Xu, C. G., Li, X. H., & Zhang, Q. F. (2005). Delimitation of the rice wide compatibility gene S5n to a 40-kb DNA fragment. Theoretical and Applied Genetics, 111, 1080–1086.

    Article  CAS  PubMed  Google Scholar 

  • Rahman, S., Bird, A., Regina, A., Li, Z., & Philippe, R. J. (2007). Resistant starch in cereals: Exploiting genetic engineering and genetic variation. Journal of Cereal Science, 46, 251–260.

    Article  CAS  Google Scholar 

  • Raigond, P., Ezekiel, R., & Raigond, B. (2015). Resistant starch in food: A review. Journal of the Science of Food and Agriculture, 95, 1968–1978.

    Article  CAS  PubMed  Google Scholar 

  • Sajilata, M. G., Singhal, R. S., & Kulkarni, P. R. (2006). Resistant starch–a review. Compr Rev Food Sci F, 2006, 5,1–5,17.

    Google Scholar 

  • Sawada, T., Francisco Jr., P. B., Aihara, S., Utsumi, Y., Yoshida, M., Oyama, Y., Tsuzuki, M., Satoh, H., & Nakamura, Y. (2009). Chlorella starch branching enzyme II (BEII) can complement the function of BEIIb in rice endosperm. Plant & Cell Physiology, 50, 1062–1074.

    Article  CAS  Google Scholar 

  • Shen, W. Q., Shu, X. L., Zhang, L. L., & Wu, D. X. (2006). Development and characteristics of processing functional indica early rice cultivar “ZHEFU 201”. Nuclear Agricultural Sciences, 20, 312–314 (in Chinese with English abstract).

    Google Scholar 

  • Shen, Y. J., Jiang, H., Jin, J. P., Zhang, Z. B., Xi, B., He, Y. Y., Wang, G., Wang, C., Qian, L. L., Li, X., Qing-Bo Yu, Q. B., Liu, H. J., Chen, D. H., Gao, J. H., Huang, H., Shi, T. L., Song, Z. P., Xu, X., Wang, B., Chen, J. K., & Lu, B. R. (2003). Genetic diversity in the northernmost Oryza rufipogon populations estimated by SSR markers. Theor App Gene, 107, 1492–1499.

    Article  Google Scholar 

  • Sun, C. L., Bai, J. J., & Shi, B. (2012). Analysis on combining ability of resistant starch content in rice grain. Chinese Agricultural Science Bulletin, 28, 24–28 (in Chinese with English abstract).

    CAS  Google Scholar 

  • Tetlow, I. J., & Emes, M. J. (2014). A review of starch-branching enzymes and their role in amylopectin biosynthesis. IUBMB Life, 66, 546–558.

    Article  CAS  PubMed  Google Scholar 

  • Tian, Q. Q., Lu, C. M., Li, X., & Fang, X. W. (2015). Low temperature treatments of rice (Oryza sativa L.) anthers changes polysaccharide and protein composition of the anther walls and increases pollen fertility and callus induction. Plant Cell Tiss Organ Cult, 120, 89–98.

    Article  CAS  Google Scholar 

  • Umemoto, T., Yano, M., Sotoh, H., Shomura, A., Nakamuta, Y. (2002). Mapping of a gene responsible for the difference in amyloptctin structure between japonica-type and indica-type rice varieties. Theor Appl Genet, 104(1), 1–8.

    Article  CAS  PubMed  Google Scholar 

  • Wang, C. L., Zhang, Y. D., Zhu, Z., Tao, C., Zhao, L., Lin, J., & Zhou, L. H. (2009). Development of a new japonica rice variety Nanjing 46 with good eating quality by marker assisted selection. Molecular Plant Breeding, 7, 1070–1076 (in Chinese with English abstract).

    CAS  Google Scholar 

  • Wani, A. A., Manzoor, P. S., Shah, A., Schweiggert-Weisz, U., Gul, K., & Wani, I. A. (2012). Rice starch diversity: Effects on structural, morphological, thermal, and physicochemical properties—A review. Compr Rev Food Sci F, 11, 417–436.

    Article  CAS  Google Scholar 

  • Wei, C., Qin, F., Zhu, L., Zhou, W., & Chen, Y. (2010). Microstructure and ultrastructure of high-amylose rice resistant starch granules modified by antisense RNA inhibition of starch branching enzyme. J Agri Food Chem, 58, 1224–1232.

    Article  CAS  Google Scholar 

  • Wu, C. H., Liang, C. Y., & Ye, X. L. (2004). In vitro megagametophyte development via unfertilized ovule culture in Doritis pulcherrima. Acta Botanica Sinica, 46, 839–845 (in Chinese with English abstract).

    CAS  Google Scholar 

  • Xia, H., Yandeau-Nelson, M., Thompson, D. B., & Guiltinan, M. J. (2011). Deficiency of maize starch-branching enzyme I results in altered fine structure, decreased digestibility and reduced coleoptile growth during germination. BMC Plant Biology, 11, 95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • Yanagihara, S., Mccouch, S. R., Ishikawa, K., Ogi, Y., Maruyama, K., & Ikehashi, H. (1995). Molecular analysis of the inheritance of the S-5 locus, conferring wide compatibility in Indica/japonica hybrids of rice (Oryza sativa L.). Theoretical and Applied Genetics, 90, 182–188.

    Article  CAS  PubMed  Google Scholar 

  • Yang, C. Z., Shu, X. L., Zhang, L. L., Wang, X. Y., & Zhao, H. J. (2006). Starch properties of mutant rice high in resistant starch. J Agri Food Chem, 54, 523–528.

    Article  CAS  Google Scholar 

  • Yang, J., Wang, J., Cao, Q., Chen, Z. D., & Zhong, W. G. (2009). Development and application of a functional marker for wide compatibility gene S5-n of Rice. Acta. Agrono Sina, 35, 2000–2007 (in Chinese with English abstract).

    Article  CAS  Google Scholar 

  • Yang, R., Sun, C., Bai, J., Luo, Z., & Shi. (2012). A putative gene sbe3-rs for resistant starch Mutated from SBE3 for starch branching enzyme in Rice (Oryza sativa L.). PLoS One, 7(8), e43026. https://doi.org/10.1371/journal.pone.0043026.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • Zeng, Y. W., Sun, D., Du, J., Pu, X. Y., Yang, S. M., Yang, X. M., Yang, T., & Yang, J. Z. (2016). Identification of QTLs for resistant starch and total alkaloid content in brown and polished rice. Genetics and Molecular Research, 15(3), gmr.15037268.

    Article  CAS  Google Scholar 

  • Zhao, Y. X., Wu, X. X., Li, W., Li, R., Shi, W. J., & Guo, Y. D. (2012). Research progress in androgenesis and haploid breeding of gramineous forage and turfgrass plants. Bulletin of Botany, 47, 87–99 (in Chinese with English abstract).

    Article  CAS  Google Scholar 

  • Zhou, H. G., Wang, L. J., Liu, G. F., Meng, X. B., Jing, Y. H., Shu, X. L., Kong, X. L., Sun, J., Yu, H., Smith, S. M., Wu, D. X., & Li, J. Y. (2016). Critical roles of soluble starch synthase SSIIIa and granule-bound starch synthase waxy in synthesizing resistant starch in rice. PNAS, 45(113), 12844–12849.

    Article  CAS  Google Scholar 

  • Zhu, H. M., Bai, J. J., Wang, H., Li, M. B., & Hao, Z. B. (2010). Starch properties of japonica rice enriched with resistant starch. Chinese Science Bulletin, 26, 108–112 (in Chinese with English abstract).

    Google Scholar 

Download references

Acknowledgements

This work was supported by grants from the national key research and development plan of China (2016YFD0300501-03) and the National Natural Science Foundation of China (31571585). The authors thank the anonymous reviewers and editorial staff for their time and attention.

Funding

This work was supported by grants from the national key research and development plan of China (2016YFD0300501–03) and the National Natural Science Foundation of China (31571585).

Author information

Authors and Affiliations

Authors

Contributions

LX supervised the project, designed and guided the experiment, and wrote the manuscript. TQQ conducted most of the experiments. LCM helped with some experiments and revised the manuscript. FXW guided some experiments. All authors have read and approved the final manuscript.

Corresponding author

Correspondence to X. Li.

Ethics declarations

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1:

Table S1. Information on 35 indel markers (Shen et al., 2003).

Additional file 2:

Table S2. Indel classification standard of dividing the indica and the japonica variety by (Fi) or (Fj) gene frequency.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Tian, Q.Q., Li, X., Lu, C.M. et al. Breeding Rice lines for physio-functional food through indica ‘Zhaxima’ × japonica ‘Nanjing 46’ haploid technique. Food Prod Process and Nutr 1, 12 (2019). https://doi.org/10.1186/s43014-019-0010-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s43014-019-0010-7

Keywords