- Yunnan Key Laboratory for Rice Genetic Improvement, Food Crops Research Institute, Yunnan Academy of Agricultural Sciences, Kunming, China
Rice improvement depends on the availability of genetic variation, and AA genome Oryza species are the natural reservoir of favorable alleles that are useful for rice breeding. To systematically evaluate and utilize potentially valuable traits of new QTLs or genes for the Asian cultivated rice improvement from all AA genome Oryza species, 6,372 agronomic trait introgression lines (ILs) from BC2 to BC6 were screened and raised based on the variations in agronomic traits by crossing 170 accessions of 7 AA genome species and 160 upland rice accessions of O. sativa as the donor parents, with three elite cultivars of O. sativa, Dianjingyou 1 (a japonica variety), Yundao 1 (a japonica variety), and RD23 (an indica variety) as the recurrent parents, respectively. The agronomic traits, such as spreading panicle, erect panicle, dense panicle, lax panicle, awn, prostrate growth, plant height, pericarp color, kernel color, glabrous hull, grain size, 1,000-grain weight, drought resistance and aerobic adaption, and blast resistance, were derived from more than one species. Further, 1,401 agronomic trait ILs in the Dianjingyou 1 background were genotyped using 168 SSR markers distributed on the whole genome. A total of twenty-two novel allelic variations were identified to be highly related to the traits of grain length (GL) and grain width (GW), respectively. In addition, allelic variations for the same locus were detected from the different donor species, which suggest that these QTLs or genes were conserved and the different haplotypes of a QTL (gene) were valuable resources for broadening the genetic basis in Asian cultivated rice. Thus, this agronomic trait introgression library from multiple species and accessions provided a powerful resource for future rice improvement and genetic dissection of agronomic traits.
Introduction
Rice is one of the most important staple crops for almost half of the world’s population. The Food and Agriculture Organization of the United Nations predicts that rice yield will have to be increased 50 to 70% by 2050 to meet human’s demands, which increases that rice yield is still central for maintaining global food security (Ray et al., 2013). Whereas rice yield potential has been stagnant since the introduction of semidwarf gene into cultivated rice and the utilization of heterosis (Virmani et al., 1982; Monna et al., 2002; Sasaki et al., 2002), the narrow genetic basis that results from the overuse of few parental materials and the lack of favorable variations led to yield bottleneck in rice breeding (Tanksley and Mccouch, 1997).
Genus Oryza contains twenty-two wild species and two cultivated rice species that represent 11 genomes: AA, BB, CC, BBCC, CCDD, EE, FF, GG, HHJJ, HHKK, and KKLL (Khush, 1997). Among these, six wild rice species (O. nivara, O. rufipogon, O. barthii, O. glumaepatula, O. longistaminata, and O. meridionalis) and two cultivated species (O. sativa and O. glaberrima) were classified into the AA genome. Asian cultivated rice (O. sativa L.) was domesticated from wild species O. rufipogon thousands of years ago (Khush, 1997; Huang et al., 2012). Previous reports indicated that 40% of the alleles of O. rufipogon was lost during the domestication from common wild rice to the cultivated rice (Sun et al., 2002), and only 10–20% of the genetic diversity in O. rufipogon and O. nivara was retained in two subspecies of the cultivated rice (Zhu et al., 2007). Since sharing the same AA genome, O. glaberrima and the six wild rice species are the most accessible gene pool for rice improvement (Ren et al., 2003). Thus, the exploitation and utilization of the useful alleles of AA genome species may overcome yield plateaus of O. sativa (Xiao et al., 1998). However, it is difficult to utilize the natural genetic diversity because of reproductive isolation, linkage drag, and background noise. Moreover, many important agronomic traits that include yield are controlled by quantitative trait loci (QTL) with smaller effect, which can be influenced by the environment. It is difficult to understand the QTL-controlled agronomic traits because of their complex inheritance and the genetic background noise.
Introgression lines are genetic resource in which the whole genome of a donor genotype (DG) is represented by the different segments in the genetic background of elite varieties. Genetic background noise of ILs can be eliminated significantly, which can be evaluated for any traits’ improvement over the recurrent parents for rice breeding, also for QTL mapping and gene discovering as a single Mendelian factor; in addition, potential favorable genes hidden in the background of related species could be expressed in the genetic background of cultivated rice (Ballini et al., 2007; Eizenga et al., 2009; Bian et al., 2010; Rama et al., 2015; Jin et al., 2016; Yang et al., 2016; Bhatia et al., 2017; Bhatia et al., 2018; Yamagata et al., 2019). Thus, ILs that eliminate hybrid sterility, linkage drags, and background noise are one of the most important genetic resources for QTL mapping, gene identification, and discovery and rapid utilization for commercial breeding. Though a series of introgression lines developed by the genome-wide marker selections were obtained from the intersubspecific crosses between indica and japonica varieties and from the interspecific crosses between Asian cultivated rice and wild relatives of Oryza sativa (Bhatia et al., 2017; Divya et al., 2019), some lines showed the remarkable phenotype, whereas others did not exhibit obvious agronomic traits, which were difficult to be used for QTL identification, gene cloning, and breeding improvement. Establishing the introgression lines based on agronomic trait selection might be time-consuming and less laborious strategy.
AA genome species distributed in the natural and wild environment, which contains amount of useful allelic genes for improving rice yield and resistance to biotic and abiotic stresses (Khush, 1997). Different AA genome wild species and native varieties with unique characteristics and ecological adaptability represented the independent center of genetic diversity in rice. Comprehensively and systemically developing the ILs from all the AA genome species in different elite cultivar varieties background will help us to achieve sustainable yield improvement, diverse requirements for quality, and broad-spectrum resistance so as to meet the demand of the modern breeding program.
In this study, to explore and utilize wild relatives in rice improvement, we systematically introduced foreign segments from eight different AA genome species (O. longistaminata, O. barthii, O. glumaepatula, O. meridinalis, O. nivara, O. rufipogon, O. glaberrima, and upland rice of O. sativa) into three elite, highly productive O. sativa varieties (Dianjingyou 1, Yundao 1 and RD23). A total of six thousand three hundred and seventy-two agronomic ILs in three different backgrounds were screened and developed based on the repeated evaluation and selection of agronomic traits. One thousand four hundred and one of 6,372 agronomic ILs in the Dianjingyou 1 background were used to analyze genotype and discover novel alleles for grain size. Thus, this agronomic introgression library provided a powerful resource for future rice improvement and genetic dissection of agronomic traits.
Materials and Methods
Plant Materials
The plant materials included 1 accession of O. longistaminata, 13 accessions of O. barthii, 6 accessions of O. glumaepatula, 8 accessions of O. meridionalis, 19 accessions of O. rufipogon, 20 accessions of O. nivara, 103 accessions of O. glaberrima, and 160 upland rice varieties of O. sativa (Supplementary Table 1). Three elite varieties, Dianjingyou 1 (a japonica variety), Yundao 1 (a japonica variety), and RD23 (an indica variety), were used as the recurrent parents.
A total of three hundred and twenty-nine accessions of AA genome species as the donor parents except for O. longistaminata were crossed with Dianjingyou 1 as the recurrent parent. A total of two hundred and twenty-six accessions as the donor parents, except for O. longistaminata and O. glaberrima, were used to cross with the recurrent parent Yundao 1. All the F1 plants were used as female parents to backcross to their respective recurrent parents to produce BC1F1 generation. More than 200 BC1F1 seeds were generated for each of the combinations. The moderate heading date of individuals was selected to backcross with the recurrent parents, and about 200 BC2F1 seeds were obtained. From each of the BC2F1 progeny, individuals that showed a significant agronomic difference from the recurrent parents were selected for further backcrossing or selfing. After 2–6 times backcrossing and 2–7 times selfing, the progeny with stable and different target traits from their recurrent parents was developed as agronomic ILs.
The F1 plants were obtained by embryo rescue technique from the cross between 1 accession of O. longistaminata as the donor parent and an indica variety RD23 as the recurrent parent, and crossing and selfing from BC1F1 generation were performed as above mentioned procedure.
All materials were grown at the Sanya Breeding Station, Sanya (18.24° N, 109.50° E), Hainan province, China. Ten individuals per row were planted at a spacing of 20 cm × 25 cm. All materials were grown and managed according to the local protocol.
Agronomic Trait Evaluation
A randomized complete block design was carried out with three replications for agronomic trait evaluation under two different environments (E1: winter and dry season, December to April 2008–2009; E2: summer and rainy season, July to November 2009), respectively. Each line was planted in three rows with 10 individuals per row. The five plants in the middle of each row were used for scoring traits. The recurrent parents, Dianjingyou 1, Yundao 1, and RD23, were used as controls in the experiment, respectively.
Prostrate growth habit was observed for the tiller angle in three main stages, which includes booting stage, heading stage, and grain filling stage. That tiller angle in ILs was larger than that in recurrent parent, which was regarded as the prostrate growth.
Primary branches at the base of panicle of the lines extend outward were regarded as the spreading panicle. Erect panicle or drooping panicle was evaluated according to the angle between the lines that connecting panicle pedestal with panicle tip and the elongation line of stem; spikelet numbers were measured as the total number of spikelets of the whole plant divided by its total number of panicles. Dense panicle was scored by the ratio of spikelet numbers to panicle length.
Tiller number was recorded from five random plants; plant height was measured from the ground level to the tip of the tallest panicle.
To measure the grain size, grains were selected from primary panicle and stored at room temperature for at least 3 months before testing. Twenty grains were used to measure grain length (GL), grain width (GW), and the ratio of grain length to grain width (RLW) from each plant. Photographs of grains per individual were taken using stereomicroscope, and then, grain size was measured by software Image J. The average value of 20 grains was used as phenotypic data. The weight of one thousand grains was measured by weighting fertile, fully mature grains from five panicles.
Aerobic adaptation was evaluated by biomass, yield, harvest index, heading date, and plant height in both aerobic and irrigated environments. Drought tolerance was assessed by the same traits as the aerobic adaptation in both upland and irrigated environments. For the aerobic and upland treatments, we used direct sowing with 4 seeds per hole and retained one seedling at the three-leaf stage. There is some difference in water management, and the rainfall provided the essential water for plant growth without the extra irrigation under aerobic treatment, whereas mobile sprinkler irrigation facilities were used to maintain a humid soil environment at the sow, tiller, and heading stage under upland treatment. For the irrigating treatment, sowing and transplanting single seedlings were done, and the field was managed according to the local standard practices.
To evaluate blast resistance, introgression lines were inoculated with Magnaporthe oryzae for 3 weeks after sowing by spraying with conidial suspension. After 7 days, lesion types on rice leaves were observed and scored according to a standard reference scale based on a dominant lesion type (Xu P. et al., 2015).
For a simply inherited trait, awn, pericarp color, and kernel color were observed directly in the field.
DNA Extraction and PCR Protocol
The experimental procedure for DNA extraction was performed as previously described (Edwards et al., 1991); A total of 168 SSR markers were selected from the Gramene database1 or previously published polymorphic SSR markers within the Oryza AA genome species (McCouch et al., 2002; Orjuela et al., 2010). PCR was performed as follows: a total volume of 10 μl containing 10 ng template DNA, 1 × buffer, 0.2 μM of each primer, 50 μM of dNTPs, and 0.5 unit of Taq polymerase (Tiangen Company, Beijing, China). The reaction mixture was incubated at 94°C for an initial 4 min, followed by 30 cycles of 94°C 30 s, 55°C 30 s, and 72°C 30 s, and a final extension step of 5 min at 72°C. PCR products were separated on 8% non-denaturing polyacrylamide gel and detected using the silver staining method.
Determination of the Length of a Substituted Segment in Introgression Lines
The substituted segment was counted based on the SSR markers distributed on twelve chromosomes (Paterson et al., 1988; Young and Tanksley, 1989; McCouch et al., 2002). Intervals between two markers homozygous for the DG were regarded as 100% introgression segment, and a chromosome segment flanked by one marker of the DG and one marker of the recurrent type (DR) was considered as 50% introgression segment, whereas intervals between two markers homozygous for the recurrent genotype (RG) represented the background genotype. Thus, the physical distance of both DD and half of DR was used to estimate the length of introgression segments. The expected introgression length of the genome is divided by the total genome size to yield the expected proportion of introgression.
Exploring Loci for Grain Size Based on the Introgression Lines
Statistical analyses were performed on the SAS software package. The linkage between loci and grain size was scored by binomial distribution based on the genotype and phenotype between ILs and the recurrent parent. The genotypes of ILs that showed a significant difference from the recurrent parent in grain size were used to perform QTL analysis, if the rate of DG at some loci was significantly higher than that of the theoretical prediction, this locus could be linked with the grain size. Significant level was determined by the comparison between the ILs and the recurrent parents. Dunnett’s t-test at p < 0.0001 was set to decrease the false probability (Eshed and Zamir, 1995; Xu et al., 2014).
Results
Agronomic Trait Diversity in Introgression Library From AA Genome Donors
To systematically explore potentially valuable genes hidden in the AA genome wild relatives and two cultivated species, 67 accessions of AA genome wild rice, 103 accessions of O. glaberrima, and 160 upland rice varieties of O. sativa as the donors were used to raise the agronomic introgression line library. Of these accessions, all accessions except for O. longistaminata were used for generating ILs in the Dianjingyou 1 background, 233 donors except for the accessions of O. glaberrima and O. longistaminata were used for developing ILs in the background of Yundao 1, and 1 accession of O. longistaminata was used for raising ILs in the RD23 background. A total of 6,372 introgression lines with multiple donors showed a remarkable difference in the agronomic traits, which includes spreading panicle, erect panicle, dense panicle, lax panicle, awn, prostrate growth, plant height, pericarp color, kernel color, glabrous hull, grain size, 1,000-grain weight, drought resistance, and aerobic adaption, compared with their recurrent parents (Figures 1A–D). Among these, 74, 61, 179, 824, 135, 251, and 1561 ILs that show distinguished traits in the Dianjingyou 1 background were selected from the donors of O. barthii, O. glumaepatula, O. meridionalis, O. rufipogon, O. nivara, O. glaberrima, and upland rice of O. sativa, respectively (Figures 1A,C and Supplementary Tables 2, 3). Additionally, 244, 85, 547, 714, 858, and 825 ILs that exhibit different agronomic traits in the Yundao 1 background were developed from the donors of O. barthii, O. glumaepatula, O. meridionalis, O. rufipogon, O. nivara, and upland rice of O. sativa, respectively (Figure 1A and Supplementary Table 4). A total of two hundred and sixty-five ILs were derived from the cross between 1 accession of O. longistaminata as the donor and an indica variety RD23 as the recurrent parent (Figure 1D and Supplementary Table 5). Thus, the agronomic introgression library derived from the multiple donors of AA genome species in the different backgrounds showed the abundant genetic variations in the agronomic traits.
Figure 1. The summary of introgression libraries with the donor of 8 AA genome species in three different genetic backgrounds. (A) A number of ILs from different donors in the Dianjingyou 1, Yundao 1, and RD23 background. (B) A number of ILs derived from O. barthii, O. glumaepatula, O. meridionalis, O. nivara, O. rufipogon, and upland rice, respectively, showed agronomic traits distinguishing from their recurrent parents Dianjingyou 1 and Yundao 1. (C) Agronomic traits of ILs derived from crosses between O. glaberrima and Dianjingyou 1. (D) Agronomic traits of ILs derived from cross between O. longistaminata and RD23.
For the same donor parent, phenotype variations for the agronomic traits varied with the genetic background. The numbers of ILs that show erect panicle, dense panicle, lax panicle, awn, plant height, pericarp color, 1,000-grain weight, drought-resistance, and aerobic adaption in Yundao 1 background was more than those of Dianjingyou 1 background, whereas the number of ILs that exhibit spreading panicle, prostrate growth, kernel color, glabrous hull, GL, and GW in Yundao 1 background was less than those of Dianjingyou 1 background (Figure 1B). It suggested that target trait expression was depended on the background of recurrent parent difference to a certain degree. Developing an introgression library in the different backgrounds will be beneficial to express hidden genes in the donor and discover more genetic variations for further study.
Characteristics of Chromosome Substituted Segments in the Introgression Library
A total of 168 SSR markers distributed on 12 chromosomes were selected to genotype agronomic introgression library in Dianjingyou 1 background (Supplementary Figure 1). The length of the interval between two markers ranged from 0.2 to 5.5 Mb, with an average of 2.22 Mb on the rice physical map (Table 1 and Supplementary Figure 1). The polymorphism rate displayed from 82.74 to 98.43% between seven AA genome species and Dianjingyou 1 (Table 1).
A total of one thousand four hundred and one IL in the Dianjingyou 1 background were used to detect the characteristics of chromosome segments from seven AA genome species, which include 29 ILs from O. barthii, 30 ILs from O. glumaepatula, 76 ILs from O. meridionalis, 380 ILs from O. nivara, 74 ILs from O. rufipogon, 81 ILs from O. glaberrima, and 731 ILs from upland rice of O. sativa (Supplementary Table 6). In the 29 ILs from the donor of O. barthii, the length of introgression segments ranged from 2.66 to 28.98 Mb, averaging 6.99 Mb (Table 2). Different coverage rate was observed in a different chromosome. Chromosomes 3, 8, and 9 had 100% coverage rate, whereas chromosome 11 only had 39.14% coverage rate (Supplementary Figure 2 and Supplementary Table 7).
A total of thirty ILs were obtained from an interspecific backcross between the cultivated rice O. sativa Dianjingyou 1 and the wild relative O. glumaepatula, and the length of introgression segments ranged from 660 kb to 25.6 Mb, averaging 5.35 Mb (Table 2). The ILs covered 73.11% of O. glumaepatula genome in the Dianjingyou 1 background (Supplementary Figure 3 and Supplementary Table 8).
In the Dianjingyou 1/O. meridionalis introgression library, the length of introgression segments was detected from 0.27 to 23.77 Mb, averaging 5.83 Mb (Table 2). The donor introgressions covered 89.17% O. meridionalis genome. Chromosomes 3, 6, and 8 exhibited complete coverage, whereas chromosome 11 showed the least coverage rate of 60.71% (Supplementary Figure 4 and Supplementary Table 9).
A total of three hundred and eight ILs were developed with the donor of O. nivara, and the average length of introgression segments was 5.84 Mb (Table 2). Average coverage rate per chromosome was about 97.17% (Supplementary Figure 5 and Supplementary Table 10).
In the introgression library with the donor of O. rufipogon, the length of introgression segments varied from 280 kb to 23.46 Mb, averaging 6.37 Mb, and the average number of segments per chromosome was 52.08 (Table 2). Chromosomes 3, 8, and 9 exhibited 100% coverage rate, whereas chromosome 2 showed the least coverage rate of 74.43% (Supplementary Figure 6 and Supplementary Table 11).
We developed 81 ILs with the O. glaberrima species as the donor, and the length of introgression segments varied from 270 kb to 24.17 Mb, averaging 5.80 Mb (Table 2). Chromosomes 3, 6, 9, and 11 exhibited complete coverage, whereas chromosome 10 showed the least coverage rate of 78.65% (Supplementary Figure 7 and Supplementary Table 12).
In the introgression library derived from the donor of upland rice in O. sativa, the average length of introgression segments was 6.67 Mb (Table 2). All the chromosomes except for chromosome 10 showed 100% upland rice genome coverage (Supplementary Figure 8 and Supplementary Table 13).
Taken together, ILs covered 81.19%, 73.11%, 89.17%, 97.17%, 89.19%, 90.70%, and 99.10% of O. barthii, O. glumaepatula, O. meridionalis, O. nivara, O. rufipogon, O. glaberrima, and upland rice genome information, respectively, which suggests that systematic and comprehensive agronomic IL library with the donor of AA genome species was developed by agronomic trait selection.
Detection of Potential Allelic Variations for Grain Size in the Introgression Library
Seed size plays an important role in rice yield (Xing and Zhang, 2010). Grain size not only determines rice appearance, but also affects milling, cooking, and eating quality of rice (Fan et al., 2006). Significant variations were observed for GL, GW, and the RLW in the introgression library with multiple donors in the background of Dianjingyou 1. Some ILs for GL, GW, and RLW were found to be significantly superior to the recurrent parent Dianjingyou 1. For GL, 133 and 125 ILs were found to be significantly longer than the Dianjingyou 1 in two seasons, respectively. For GW, 412 and 508 ILs were observed to be significantly wider than the recurrent parent in two environments, respectively. For the RLW, 277 and 178 ILs were found to be higher than the Dianjingyou 1 in different seasons, respectively (Figure 2). In addition, the same traits in the different environments showed a highly significant correlation (Supplementary Table 14). These results suggested that abundant genetic variations for grain size existed in the wild and cultivated accessions of rice.
Figure 2. Frequency distribution of grain size traits in introgression lines in the Dianjingyou 1 background. Arrows indicated the mean values for recurrent parent Dianjingyou 1.
To explore favorable allelic variation for grain size, QTLs were detected based on the genotype and phenotype data. A total of forty-one loci linking with GL, forty-four loci linking with GW, and thirty-two loci linking with RLW were identified. It indicated that abundant gene pool for grain size existed in the AA genome species (Figures 3–5). Among these, 26 loci for GL were detected from multiple donors, and 12, 11, 2, and 1 loci were found from the donors of two species, three species, four species, and six species, respectively. It suggested that the same locus that contributes to GL is a potential allelic variation from different donors. Moreover, 4 loci from the different donors were only responsible for long grain, 12 loci derived from the multiple donors only contributed to short grain, and 10 loci from the different species controlled both long grain and short grain. Moreover, 22 novel allelic variations from multiple donors contributed to GL (Figure 3).
Figure 3. Allelic variation linking with GL was detected from the different donors. Asterisk indicated published genes.
Figure 4. Allelic variation linking with GW was detected from the different donors. Asterisk indicated published genes.
A total of 27 loci for GW were examined from the multiple donors, which include 13 loci from the only two species, 9 loci from three species, 3 loci from four species, and 1 locus from five species (Figure 4). Moreover, 12 loci from the different donors were only responsible for wide grain, 4 loci from multiple species only led to thin grain, and 10 loci from the different species controlled both wide and thin grains. In addition, 22 novel allelic variations for GW were found in agronomic IL library (Figure 4).
Nineteen loci for RLW from multiple donors were explored on 12 chromosomes, which include 12, 6, and 1 locus detected simultaneously in two, three, and four donor species, respectively (Figure 5).
These results indicated that detection of favorable genes using multiple donors could help us to find the novel allelic variations. The allelic genes were detected in the different donors, which suggest that some loci for grain size were conserved in Genus Oryza. Some loci controlled the opposite phenotype, long grain vs. short grain, wide grain vs. thin grain, validating these loci’s functions in forward and reverse direction and also suggesting that the loci functioned divergence in the different donors. Taken together, these results would provide the information that the loci for grain size from the different donors were the same or different haplotypes; it also indicated that agronomic IL library with the donor of 7 AA genome species was an excellent resource and tool to discover favorable allelic variations and new QTLs or genes for rice improvement.
Discussion
Agronomic Introgression Line Library for All AA Genome Species Is an Important Stock for Breeding Improvement in Rice
The geographical distribution of wild relative species with the AA genome is in a wide range of environments, O. nivara and O. rufipogon mainly in Asia, O. glaberrima, O. barthii, and O. longistaminata in Africa, O. meridionalis in Australia, O. glumaepatula in Latin America, and wild relative species evolve a large extend on genetic differentiation and morphological variation under the different ecological environments (Vaughan et al., 2005). Constructing introgression lines are a feasible and effective approach to transferring the favorable genes from wild relative species to the cultivated varieties for improving elite cultivars. By now, more than 40 sets of ILs were developed using the AA genome wild species as the donor parents (Chen et al., 2006; Tian et al., 2006; McCouch et al., 2007; Rangel et al., 2008; Hao et al., 2009; Ali et al., 2010; Gutierrez et al., 2010; Ramos et al., 2016; Bhatia et al., 2017), but ILs were almost derived from a single accession of AA genome species in a single background, which leads to the lack of systematic utilization of favorable genes. One of the challenges of constructing interspecific introgression lines was to overcome interspecific hybrid sterility. In this study, we selected the typical 330 accessions of AA genome species distributed in the different geographical region as the donor parents to raise agronomic IL library. We observed that pollen fertility of F1 varied from 1.92% to 93.19% dependent on the different accessions of O. nivara and O. rufipogon, whereas all the crosses with the accessions of O. barthii, O. glumaepatula, and O. meridionalis showed complete pollen sterility in the F1 combinations (data not shown). When the japonica varieties Dianjingyou 1 and Yundao 1 used as the recurrent parents were crossed with the accession of O. longistaminata, the crossing was failed despite many efforts. Only the cross using an indica variety RD23 as the recurrent parent and O. longistaminata as the donor was obtained by embryo rescue. Fortunately, the female gametes from the interspecific hybrids were partially fertile, and some hybridization seeds in the different combinations could be harvested by backcrossing F1 as the female parent with O. sativa as the male parent. Finally, agronomic IL library that contains 6,372 lines was developed based on the agronomic trait selection, and agronomic ILs showed a significant difference from the recurrent parents, which include spreading panicle, erect panicle, dense panicle, lax panicle, awn, prostrate growth, plant height, pericarp color, kernel color, glabrous hull, grain size, 1,000-grain weight, drought resistance and aerobic adaption, and blast resistance (Figure 1), and this IL library help us to understand the genetic base of agronomic traits and explore the favorable genes or allelic variations, and also breeding improvement systematically and comprehensively.
Variation for agronomic traits in this IL library shows the importance of AA genome species for further breeding improvement in rice. For example, 18 ILs in Dianjingyou 1 background and 57 ILs in Yundao 1 background showed dense panicle (Supplementary Tables 2–4), 137 ILs registered significant improvement over the recurrent parents in grain weight (data not shown), and these ILs derived from multiple donors could contribute to the variations for yield.
The upland rice is a predominant ecotype adapted to aerobic and rain-fed conditions in the mountainous areas that have high genetic variability in the characteristics of morphology and physiology, such as glabrous hull and aerobic adaptation (Bridhikitti and Overcamp, 2012; Sandhu et al., 2013). In this study, 160 accessions of upland rice that represent an abundant genetic diversity were used to raise introgression lines. Among those agronomic ILs, 93 and 18 ILs in the Diangjingyou 1 and Yundao 1 background exhibited the glabrous hull phenotype. Those resources were not only good for breeding varieties suitable for agricultural operation, but also helpful to understand the genetic mechanism of glabrous hull development. A total of 125 ILs in aerobic adaptation were superior to the recurrent parents, which are important and environmental-friendly breeding materials to meet the need of aerobic rice development.
The effect of alleles on the agronomic traits varied with genetic background. The different recurrent parents helped us to find background-dependent useful traits or stable traits in the different backgrounds. In this study, with the donor accessions of O. barthii, O. glumaepatula, O. meridionalis, O. nivara, O. rufipogon, and upland rice of O. sativa, 62, 37, 147, 470, 122, and 710 ILs that conferred agronomic traits were found in the genetic background of the Dianjingyou 1 and Yundao 1, which suggested that the genes for agronomic traits had a stable effect on the different genetic backgrounds (Supplementary Tables 2–4). In addition, we found ILs that showed aerobic adaptation with the donors of O. barthii, O. glumaepatula, O. meridionalis, O. nivara, O. rufipogon were detected in the Yundao 1 background, rather than Dianjingyou 1 background (Supplementary Tables 2–4). The results would provide the theoretical guidance of the relationship between the traits and genetic background in rice breeding.
Exploration of Natural Allelic Variations Using Agronomic Introgression Line Library
Genetic diversity and allele were lost during the domestication from the wild species of rice to the cultivated rice (Sun et al., 2002), whereas narrow genetic basis led to the yield bottleneck of Asian cultivated rice. In the recent years, mining and utilization of useful allele variation have made great progress in rice breeding. For example, the allelic variation in the Wx gene and SSSI was proved to contribute greatly to the differences in rice eating and cooking qualities (ECQs) in the two subspecies (Li et al., 2018). Allelic variation at the E1/Ghd7 locus allowed an expansion of the rice cultivation area through adjusting heading date (Saito et al., 2019). The allele types of BPH9 conferred varying levels of resistance to different biotypes of BPH and enabled rice to combat planthopper variation (Zhao et al., 2016). The allelic variation at the rice blast resistance (R) Pid3 locus was analyzed based on the 3K RGP sequencing data, and different strategies were developed to apply the functional Pid3 alleles to indica and japonica rice breeding (Lv et al., 2017). In this study, one locus for GL and one locus for GW were explored from the six and five different donor species, respectively. Two loci for GL, three loci for GW, and one locus for the RLW were detected from the donors of four species, respectively (Figures 3–5). Additionally, many of published genes for grain size were found based on the agronomic introgression library analysis, such as GW2 (Song et al., 2007), GL2 (Hu et al., 2015), PGL2 (Heang and Sassa, 2012b), PGL1 (Heang and Sassa, 2012a), GL3.2/CYP78A13 (Xu F. et al., 2015), GS3 (Fan et al., 2006; Takano-Kai et al., 2009;Mao et al., 2010), qGL3-1 (Qi et al., 2012), qGL3.3 (Hu et al., 2018), GS6 (Sun et al., 2013), TGW6 (Ishimaru et al., 2013), GL7 (Wang et al., 2015), and OsSPL13 (Si et al., 2016). In addition, 22 loci might be new QTLs or genes controlled for GL and GW from the different AA genome donors. Accordingly, agronomic introgression library with multiple donors from different relatives of Asian cultivated rice is a powerful resource platform to discover novel and functional allelic variations for agronomic traits. Mining the natural functional variations in the useful genes derived from the multiple donors and combing the different alleles through diversification could be useful for an accurate rice breeding program.
Therefore, agronomic IL libraries derived from the multiple donors have some advantages: (1) an abundant genetic variations were introgressed into the cultivated rice genome; (2) target genes or QTLs for the same phenotype could be validated by the different donors, and it will provide the information that these target genes or QTLs could be the same haplotype; (3) the genes or QTLs responsible for the opposite phenotypes, for example, long-grain size and short-grain size, could also be confirmed using the different populations from multiple donors, and it could be the different haplotypes. Therefore, this agronomic IL library will help us to improve rice breeding and interesting gene discovery and utilization.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
YZ and DT draft the manuscript. DT designed the research. JZ, PX, WD, and XD developed the introgression lines. YZ, YYa, YYu, JL, and QP participated the genotype and phenotype evaluation. YZ performed the data analysis. All authors reviewed and approved the final manuscript.
Funding
This research was supported by the National Natural Science Foundation of China (Grant Nos. 31991221, 31660380, and 31201196, and 32160489), the Yunnan Provincial Science and Technology Department, China (2018FG001-086, 202101AS070036, and 202101AT070193), the Yunnan Provincial Government (YNWR-QNBJ-2018-359), and the Yunnan Seed Industrialization Laboratory Program.
Conflict of Interest
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.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.856514/full#supplementary-material
Supplementary Figure 1 | Physical map of the 168 SSR markers in rice used for this study.
Supplementary Figure 2 | Chromosome segments distribution and coverage degree of introgression library from the donor of O. barthii. Each block at the row represented introgression from a total of 29 ILs on the target chromosome, regardless of introgression segments in other chromosomes, and each column represented a molecular marker locus.
Supplementary Figure 3 | Chromosome segments distribution and coverage degree of introgression library from the donor of O. glumaepatula. Each block at the row represented introgression from a total of 30 ILs on the target chromosome, regardless of introgression segments in other chromosomes, and each column represented a molecular marker locus.
Supplementary Figure 4 | Chromosome segments distribution and coverage degree of introgression library from the donor of O. meridionalis. Each block at the row represented introgression from a total of 76 ILs on the target chromosome, regardless of introgression segments in other chromosomes, and each column represented a molecular marker locus.
Supplementary Figure 5 | Chromosome segments distribution and coverage degree of introgression library from the donor of O. nivara. Each block at the row represented introgression from a total of 380 ILs on the target chromosome, regardless of introgression segments in other chromosomes, and each column represented a molecular marker locus.
Supplementary Figure 6 | Chromosome segments distribution and coverage degree of introgression library from the donor of O. rufipogon. Each block at the row represented introgression from a total of 74 ILs on the target chromosome, regardless of introgression segments in other chromosomes, and each column represented a molecular marker locus.
Supplementary Figure 7 | Chromosome segments distribution and coverage degree of introgression library from the donor of O. glaberrima. Each block at the row represented introgression from a total of 81 ILs on the target chromosome, regardless of introgression segments in other chromosomes, and each column represented a molecular marker locus.
Supplementary Figure 8 | Chromosome segments distribution and coverage degree of introgression library from the donor of upland rice. Each block at the row represented introgression from a total of 731 ILs on the target chromosome, regardless of introgression segments in other chromosomes, and each column represented a molecular marker locus.
Footnotes
References
Ali, M., Sanchez, P., Yu, S., Lorieux, M., and Eizenga, G. (2010). Chromosome segment substitution lines: a powerful tool for the introgression of valuable genes from Oryza wild species into cultivated rice (O. sativa). Rice 3, 218–234.
Ballini, E., Berruyer, R., Morel, J., Lebrun, M., Notteghem, J., and Tharreau, D. (2007). Modern elite rice varieties of the ‘Green revolution’ have retained a large introgression from wild rice around the Pi33 rice blast resistance locus. New Phytol. 175, 340–350. doi: 10.1111/j.1469-8137.2007.02105.x
Bhatia, D., Joshi, S., Das, A., Vikal, Y., Gurpreet, K., Neelam, et al. (2017). Introgression of yield component traits in rice (Oryza sativa ssp. indica) through interspecific hybridization. Crop Sci. 57, 1–17.
Bhatia, D., Wing, R. A., Yu, Y., Chougule, K., Kudrna, D., Lee, S., et al. (2018). Genotyping by sequencing of rice interspecific backcross inbred lines identifies QTLs for grain weight and grain length. Euphytica 214:41.
Bian, J., Jiang, L., Liu, L., Wei, X., Xiao, Y., Zhang, L., et al. (2010). Construction of a new set of rice chromosome segment substitution lines and identification of grain weight and related traits QTLs. Breed Sci. 60, 305–313.
Bridhikitti, A., and Overcamp, T. (2012). Estimation of Southeast Asian rice paddy areas with different ecosystems from moderate-resolution satellite imagery. Agric. Ecosyst. Environ. 146, 113–120.
Chen, J., Hafeez, U., Chen, D., Liu, G., Zheng, K., and Zhuang, J. (2006). Development of chromosomal segment substitution lines from a backcross recombinant inbred population of interspecific rice cross. Rice Sci. 13, 15–21.
Divya, B., Malathi, S., Sukumar, M., and Sarla, N. (2019). Development and use of chromosome segment substitution lines as a genetic resource for crop improvement. Theor. Appl. Genet. 132, 1–25. doi: 10.1007/s00122-018-3219-y
Edwards, K., Johnstone, C., and Thompson, C. (1991). A simple and rapid method for the preparation of plant genomic DNA for PCR analysis. Nucleic Acids Res. 19, 1349. doi: 10.1093/nar/19.6.1349
Eizenga, G., Agrama, H., Lee, F., and Jia, Y. (2009). Exploring genetic diversity and potential novel disease resistance genes in a collection of rice (Oryza spp.) wild relatives. Genet. Resour. Crop Evol. 56, 65–76.
Eshed, Y., and Zamir, D. (1995). An introgression line population of Lycopersicon pennellii in the cultivated tomato enables the identification and fine mapping of yield-associated QTL. Genetics 141, 1147–1162. doi: 10.1093/genetics/141.3.1147
Fan, C., Xing, Y., Mao, H., Lu, T., Han, B., Xu, C., et al. (2006). GS3 a major QTL for grain length and weight and minor QTL for grain width and thickness in rice encodes a putative transmembrane protein. Theor. Appl. Genet. 112, 1164–1171. doi: 10.1007/s00122-006-0218-1
Gutierrez, A., Carabali, S., Giraldo, O., Martinez, C., Correa, F., Prado, G., et al. (2010). Identification of a rice stripe necrosis virus resistance locus and yield component QTLs using Oryza sativa x O. glaberrima introgression lines. BMC Plant Biol. 10:6. doi: 10.1186/1471-2229-10-6
Hao, W., Zhu, M., Gao, J., Sun, S., and Lin, H. (2009). Identification of quantitative trait loci for rice quality in a population of chromosome segment substitution lines. J. Integr. Plant Biol. 51, 500–512. doi: 10.1111/j.1744-7909.2009.00822.x
Heang, D., and Sassa, H. (2012a). Antagonistic actions of HLH/bHLH proteins are involved in grain length and weight in rice. PLoS One 7:e31325. doi: 10.1371/journal.pone.0031325
Heang, D., and Sassa, H. (2012b). An atypical bHLH protein encoded by positive regulator of grain length 2 is involved in controlling grain length and weight of rice through interaction with a typical bHLH protein APG. Breed Sci. 62, 133–141. doi: 10.1270/jsbbs.62.133
Hu, J., Wang, Y., Fang, Y., Zeng, L., Xu, J., and Yu, H. (2015). A Rare allele of GS2 enhances grain size and grain yield in rice. Mol. Plant 8, 1455–1465. doi: 10.1016/j.molp.2015.07.002
Hu, Z., Lu, S., Wang, M., He, H., Sun, L., and Wang, H. (2018). A novel QTL qTGW3 encodes the GSK3/SHAGGY-Like kinase OsGSK5/OsSK41 that interacts with OsARF4 to negatively regulate grain size and weight in rice. Mol. Plant 11, 736–749. doi: 10.1016/j.molp.2018.03.005
Huang, X., Kurata, N., Wei, X., Wang, Z., Wang, A., and Zhao, Q. (2012). A map of rice genome variation reveals the origin of cultivated rice. Nature 490, 497–501. doi: 10.1038/nature11532
Ishimaru, K., Hirotsu, N., Madoka, Y., Murakami, N., Hara, N., and Onodera, H. (2013). Loss of function of the IAA-glucose hydrolase gene TGW6 enhances rice grain weight and increases yield. Nat. Genet. 45, 707–711.
Jin, J., Hua, L., Zhu, Z., Tan, L., Zhao, X., and Zhang, W. (2016). GAD1 encodes a secreted peptide that regulates grain number grain length and awn development in rice domestication. Plant Cell 28, 2453–2463. doi: 10.1105/tpc.16.00379
Li, Q., Liu, X., Zhang, C., Jiang, L., Jiang, M., and Zhong, M. (2018). Rice soluble starch synthase I: allelic variation, expression, function, and interaction with Waxy. Front. Plant Sci. 9:1591.
Lv, Q., Huang, Z., Xu, X., Tang, L., Liu, H., and Wang, C. (2017). Allelic variation of the rice blast resistance gene Pid3 in cultivated rice worldwide. Sci. Rep. 7:10362. doi: 10.1038/s41598-017-10617-2
Mao, H., Sun, S., Yao, J., Wang, C., Yu, S., Xu, C., et al. (2010). Linking differential domain functions of the GS3 protein to natural variation of grain size in rice. Proc. Natl. Acad. Sci. U S A. 107, 19579–19584. doi: 10.1073/pnas.1014419107
McCouch, S., Sweeney, M., Li, J., Jiang, H., Thomson, M., and Septiningsih, E. (2007). Through the genetic bottleneck: O. rufipogon as a source of trait-enhancing alleles for O. sativa. Euphytica 154, 317–339.
McCouch, S., Teytelman, L., Xu, Y., Lobos, K., Clare, K., and Walton, M. (2002). Development and mapping of 2240 new SSR markers for rice (Oryza sativa L.) (supplement). DNA Res. 9, 257–279. doi: 10.1093/dnares/9.6.257
Monna, L., Kitazawa, N., Yoshino, R., Suzuki, J., Masuda, H., and Maehara, Y. (2002). Positional cloning of rice semidwarfing gene sd-1: rice “green revolution gene” encodes a mutant enzyme involved in gibberellin synthesis. DNA Res. 9, 11–17. doi: 10.1093/dnares/9.1.11
Orjuela, J., Garavito, A., Bouniol, M., Arbelaez, J., Moreno, L., and Kimball, J. (2010). A universal core genetic map for rice. Theor. Appl. Genet. 120, 563–572. doi: 10.1007/s00122-009-1176-1
Paterson, A., Lander, E., Hewitt, J., Peterson, S., Lincoln, S., and Tanksley, S. (1988). Resolution of quantitative traits into Mendelian factors by using a complete linkage map of restriction fragment length polymorphisms. Nature 335, 721–726. doi: 10.1038/335721a0
Qi, P., Lin, Y., Song, X., Shen, J., Huang, W., Shan, J., et al. (2012). The novel quantitative trait locus GL3.1 controls rice grain size and yield by regulating Cyclin-T1;3. Cell Res. 22, 1666–1680. doi: 10.1038/cr.2012.151
Rama, S., Singh, K., Umakanth, B., Vishalakshi, B., Renuka, P., and Vijay, S. (2015). Development and identification of novel rice blast resistant sources and their characterization using molecular markers. Rice Sci. 22, 300–308.
Ramos, J., Furuta, T., Kanako, U., Niwa, C., Rosalyn, B., and Angeles, S. (2016). Development of chromosome segment substitution lines (CSSLs) of Oryza longistaminata A., Chev., Röhr in the background of the elite japonica rice cultivar Taichung 65 and their evaluation for yield traits. Euphytica 210, 151–163.
Rangel, P., Brondani, R., Rangel, P., and Brondani, C. (2008). Agronomic and molecular characterization of introgression lines from the interspecific cross Oryza sativa (BG90-2) x Oryza glumaepatula (RS-16). Genet. Mol. Res. 7, 184–195. doi: 10.4238/vol7-1gmr406
Ray, D., Mueller, N., West, P., and Foley, J. (2013). Yield trends are insufficient to double global crop production by 2050. PLoS One 8:e66428. doi: 10.1371/journal.pone.0066428
Ren, F., Lu, B., Li, S., Huang, J., and Zhu, Y. (2003). A comparative study of genetic relationships among the AA-genome Oryza species using RAPD and SSR markers. Theor. Appl. Genet. 108, 113–120. doi: 10.1007/s00122-003-1414-x
Saito, H., Okumoto, Y., Tsukiyama, T., Xu, C., Teraishi, M., and Tanisaka, T. (2019). Allelic differentiation at the E1/Ghd7 Locus has allowed expansion of rice cultivation area. Plants 8:550. doi: 10.3390/plants8120550
Sandhu, N., Jain, S., Kumar, A., Mehla, B. S., and Jain, R. (2013). Genetic variation, linkage mapping of QTL and correlation studies for yield, root, and agronomic traits for aerobic adaptation. BMC Genet. 14:104. doi: 10.1186/1471-2156-14-104
Sasaki, A., Ashikari, M., Ueguchi-Tanaka, M., Itoh, H., Nishimura, A., and Swapan, D. (2002). Green revolution: a mutant gibberellin-synthesis gene in rice. Nature 416, 701–702.
Si, L., Chen, J., Huang, X., Gong, H., Luo, J., and Hou, Q. (2016). OsSPL13 controls grain size in cultivated rice. Nat. Genet. 48, 447–456. doi: 10.1038/ng.3518
Song, X., Huang, W., Shi, M., Zhu, M., and Lin, H. (2007). A QTL for rice grain width and weight encodes a previously unknown RING-type E3 ubiquitin ligase. Nat. Genet. 39, 623–630. doi: 10.1038/ng2014
Sun, L., Li, X., Fu, Y., Zhu, Z., Tan, L., and Liu, F. (2013). GS6, A member of the GRAS gene family, negatively regulates grain size in rice. J. Integr. Plant Biol. 55, 38–949. doi: 10.1111/jipb.12062
Sun, Q., Wang, K., Yoshimura, A., and Doi, K. (2002). Genetic differentiation for nuclear mitochondrial and chloroplast genomes in common wild rice (Oryza rufipogon Griff.) and cultivated rice (Oryza sativa L.). Theor. Appl. Genet. 104, 1335–1345. doi: 10.1007/s00122-002-0878-4
Takano-Kai, N., Jiang, H., Kubo, T., Sweeney, M., Matsumoto, T., and Kanamori, H. (2009). Evolutionary history of GS3 a gene conferring grain length in rice. Genetics 182, 1323–1334. doi: 10.1534/genetics.109.103002
Tanksley, S., and Mccouch, S. (1997). Seed banks and molecular maps: unlocking genetic potential from the wild. Science 277, 1063–1066. doi: 10.1126/science.277.5329.1063
Tian, F., Li, D., Fu, Q., Zhu, Z., Fu, Y., Wang, X., et al. (2006). Construction of introgression lines carrying wild rice (Oryza rufipogon Griff.) segments in cultivated rice (Oryza sativa L.) background and characterization of introgressed segments associated with yield-related traits. Theor. Appl. Genet. 112, 570–580. doi: 10.1007/s00122-005-0165-2
Vaughan, D., Kadowaki, K., Kaga, A., and Tomooka, N. (2005). On the phlylogeny and biogeography of the genus Oryza. Breed. Sci. 55, 113–122. doi: 10.1016/j.simyco.2014.10.002
Virmani, S., Aquino, R., and Khush, G. (1982). Heterosis breeding in rice (Oryza sativa L.). Theor. Appl. Genet. 63, 373–380.
Wang, Y., Xiong, G., Hu, J., Jiang, L., Yu, H., and Xu, J. (2015). Copy number variation at the GL7 locus contributes to grain size diversity in rice. Nat. Genet. 47, 944–948. doi: 10.1038/ng.3346
Xiao, J., Li, J., Grandillo, S., Ahn, S., Yuan, L., Tanksley, S., et al. (1998). Identification of trait-improving quantitative trait loci alleles from a wild rice relative. Oryza rufipogon. Genetics 150, 899–909. doi: 10.1093/genetics/150.2.899
Xing, Y., and Zhang, Q. (2010). Genetic and molecular bases of rice yield. Annu. Rev. Plant Biol. 61, 421–442. doi: 10.1146/annurev-arplant-042809-112209
Xu, F., Fang, J., Ou, S., Gao, S., Zhang, F., and Du, L. (2015). Variations in CYP78A13 coding region influence grain size and yield in rice. Plant Cell Environ. 38, 800–811. doi: 10.1111/pce.12452
Xu, P., Zhou, J., Li, J., Hu, F., Deng, X., and Feng, S. (2014). Mapping three new interspecific hybrid sterile loci between Oryza sativa and O. glaberrima. Breed. Sci. 63, 476–482. doi: 10.1270/jsbbs.63.476
Xu, P., Zhou, J., Li, J., Zhang, Y., Hu, F., and Liu, S. (2015). Identification and mapping of a novel blast resistance gene pi57(t) in Oryza longistaminata. Euphytica 205, 95–102. doi: 10.1007/s10681-015-1402-7
Yamagata, Y., Win, K., Miyazaki, Y., Ogata, C., Yasui, H., and Yoshimura, A. (2019). Development of introgression lines of AA genome Oryza species O. glaberrima, O. rufipogon and O. nivara in the genetic background of O. sativa L. cv Taichung 65. Breed. Sci. 69, 359–363. doi: 10.1270/jsbbs.19002
Yang, D., Ye, X., Zheng, X., Cheng, C., Ye, N., and Huang, F. (2016). Development and evaluation of chromosome segment substitution lines carrying overlapping chromosome segments of the whole wild rice genome. Front. Plant Sci. 7:1737. doi: 10.3389/fpls.2016.01737
Young, N., and Tanksley, S. (1989). Restriction fragment length polymorphism maps and the concept of graphical genotypes. Theor. Appl. Genet. 77, 95–101.
Zhao, Y., Huang, J., Wang, Z., Jing, S., Wang, Y., and Ouyang, Y. (2016). Allelic diversity in an NLR gene BPH9 enables rice to combat planthopper variation. Proc. Natl. Acad. Sci. U S A 113, 12850–12855.
Keywords: rice, AA genome, introgression line, grain size, allelic variation
Citation: Zhang Y, Zhou J, Xu P, Li J, Deng X, Deng W, Yang Y, Yu Y, Pu Q and Tao D (2022) A Genetic Resource for Rice Improvement: Introgression Library of Agronomic Traits for All AA Genome Oryza Species. Front. Plant Sci. 13:856514. doi: 10.3389/fpls.2022.856514
Received: 17 January 2022; Accepted: 07 February 2022;
Published: 24 March 2022.
Edited by:
Mallikarjuna Swamy, International Rice Research Institute, PhilippinesReviewed by:
Hu Zhao, Huazhong Agricultural University, ChinaSung-Ryul Kim, International Rice Research Institute, Philippines
Copyright © 2022 Zhang, Zhou, Xu, Li, Deng, Deng, Yang, Yu, Pu and Tao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Dayun Tao, dGFvZHkxMkBhbGl5dW4uY29t
†These authors have contributed equally to this work