ORIGINAL RESEARCH article
QTLs Analysis and Validation for Fiber Quality Traits Using Maternal Backcross Population in Upland Cotton
- 1Laboratory of Cotton Genetics, Genomics and Breeding, College of Agronomy and Biotechnology/Beijing Key Laboratory of Crop Genetic Improvement/Key Laboratory of Crop Heterosis and Utilization of Ministry of Education, China Agricultural University, Beijing, China
- 2Department of Cotton Breeding, Institute of Cash Crops, Hubei Academy of Agricultural Sciences, Wuhan, China
Cotton fiber is renewable natural fiber source for textile. Improving fiber quality is an essential goal for cotton breeding project. In present study, F14 recombinant inbred line (RIL) population was backcrossed by the maternal parent to obtain a backcross (BC) population, derived from one Upland cotton hybrid. Three repetitive field trials were performed by randomized complete block design with two replicates in three locations in 2015, together with the BC population, common male parent and the RIL population. Totally, 26 QTLs in BC population explained 5.00–14.17% of phenotype variation (PV) and 37 quantitative trait loci (QTL) were detected in RIL population explaining 5.13–34.00% of PV. Seven common QTLs detected simultaneously in two populations explained PV from 7.69 to 23.05%. A total of 20 QTLs in present study verified the previous results across three environments in 2012. Particularly, qFL-Chr5-2 controlling fiber length on chromosome 5 explained 34.00% of PV, while qFL-Chr5-3 only within a 0.8 cM interval explained 13.93% of PV on average in multiple environments. These stable QTLs explaining great variation offered essential information for marker-assisted selection (MAS) to improve fiber quality traits. Lots of epistasis being detected in both populations acted as one of important genetic compositions of fiber quality traits.
Upland cotton (Gossypium hirsutum L.) is one of the most important fiber crops, which accounts for more than 92% of cultivated cotton worldwide because of its higher yield potential, and adaptation to diverse environments (Zhang T. et al., 2015). But fiber quality of Upland cotton is not as good as that of Sea Island cotton (G. barbadense L.). Because of the fiber quality traits are quantitative traits and highly affect by genetic background and growing environment. Therefore, it is an important target to improve fiber quality in cotton breeding programs.
Recombinant inbred line (RIL) population has been considered as a permanent mapping population, in which the number of homozygous alleles could be used to dissect the additive effect of allele(s). The backcross (BC) population deriving from RIL population by test cross or back cross, is suitable to be used to identify QTLs, to detect the dominant effect, and to carry out multiple trials across multi-environments and repeat observations, similarly to some extent to “immortalized” F2 population (Hua et al., 2002, 2003; Mei et al., 2005; Shang et al., 2016a). Several studies had been performed in other crops in yield and yield-component traits using BCF1 segregation populations based on RIL populations (Xiao et al., 1995; Mei et al., 2005; Larièpe et al., 2012).
Some interspecific BC1F1 populations were developed from G. hirsutum × G. barbadense for fiber quality traits QTL mapping in cotton (Said et al., 2015). For example, seven QTLs related to fiber length and fiber strength were detected using BCF1 segregating population deriving from Guazuncho 2 × VH8-4602 (Lacape et al., 2005), 44 QTLs for fiber quality traits were detected using (CCRI-8 × Pima 90-53) × CCRI-8 BCF1 interspecific population, especially including four QTLs locating on chromosome (Chr) 5 (Yang et al., 2015). A meta-analysis including 728 QTLs for fiber quality traits showed that QTL-rich regions on Chr 5, Chr 7, Chr 14, Chr 19, Chr 21, and Chr 26 were related to fiber length, fiber strength and Micronaire, and 11 QTLs for fiber length and fiber strength were distributed on Chr 21, respectively (Said et al., 2013). The majority of QTLs have also been mapped using intra-species RIL populations of Upland cotton (Wang et al., 2006; Shen et al., 2007; Wan et al., 2007; Wu et al., 2009; Sun et al., 2012; Ning et al., 2014; Tan et al., 2014; Shang et al., 2015, 2016a,b,c; Tang et al., 2015; Zhang J. et al., 2015; Zhang Z. et al., 2015; Jamshed et al., 2016; Li et al., 2016). Recently, intraspecific BCF1 populations in Upland cotton were constructed to map QTLs based on backcross of RIL population to the parental lines, respectively (Shang et al., 2016a; Wang et al., 2016). Shang et al. (2016d) detected 62 stable QTLs related to fiber quality traits in F9 RIL populations and their maternal BCF1 progenies. Wang et al. (2016) used two parental BCF1 populations and mapped four, one and three stable QTLs in more than one environment(s) or population(s) for fiber length, fiber strength, and Micronaire, respectively.
The genome (2.5 Gb) of allotetraploid Upland cotton harbors 26 chromosomes and displays asymmetric evolution between the A subgenome (A01–A13) and the D subgenome (D01–D13) (Zhang T. et al., 2015). The rapid development of genome sequencing technology facilitated the availability of cotton genomic data of diploid species A2 and D5 (Paterson et al., 2012; Wang K. et al., 2012; Li et al., 2014), and tetraploid genomes (AD)1 and (AD)2 (Li et al., 2015; Liu et al., 2015; Yuan et al., 2015; Zhang T. et al., 2015), which provided new insights to the development of molecular markers; from microsatellite (SSR) markers to single nucleotide polymorphism (SNP) markers and other valuable resources in breeding projects. Recently, a genetic map with 2618 SNP loci was constructed by genotyping with cotton 63K SNP array (Hulse-kemp et al., 2015). Nine stable QTLs related to fiber quality traits across two environments were detected, and one fiber length hotspot on Chr 5 carrying three QTLs was observed (Li et al., 2016). Four potential candidate genes for fiber length on Chr Dt7 were found using a total of 81,675 SNP markers genotyping by sequencing through genome-wide association studies (GWAS) (Su et al., 2016). The development of genome re-sequencing researches offered new insights to cotton genome (Fang et al., 2017; Wang M. et al., 2017). It is necessary to verify QTLs accurately under maximum environments and more segregation populations, especially stable QTLs across multiple environments. Because of the importance of cotton fiber quality, it is essential to dissect the genetic mechanism of fiber quality on Upland cotton using diverse design populations, including intraspecific BCF1 segregation population.
In our previous study, segregating populations including F2:3, F2:4, RIL, and BCF1 population derived from the same hybrid “Xinza 1” had been developed and used to dissect the genetic effects for fiber quality traits (Liang et al., 2013; Shang et al., 2015, 2016d) and yield-related traits (Liang et al., 2015; Shang et al., 2016a,b,c). Markers of MUSS193-Gh388 flanking the common qFL-Chr5-2 could be further validated its role in markers assisted selection (MAS), which controlled fiber length on Chr 5 and explained 20.12 and 14.77% of PV across two locations (Liang et al., 2013). In addition, the region of SWU15194-HAU190 on Chr 9 was detected for Micronaire by two times and explained highly to 25.93% of PV; another region of PGML1917-SWU17715 on Chr 5 increased 0.29 mm fiber length with 21.09% of phenotypic variance explained (PVE); Region of NAU4034-SWU17713 on Chr 5 controlled fiber length, fiber strength, and fiber elongation at the same time, in which, a stable QTL namely qFE-Chr5-1 explained 20.41% of PV (Shang et al., 2016d). In order to resolve genetic components related fiber quality traits and reveal the genetic mechanism more details, we further developed BCF1 progenies population based on RIL population by testcrossing with the maternal parent and attempted to investigate QTLs in fiber quality traits in both RIL and BCF1 populations across three environments in present study.
Materials and Methods
Plant Materials and Population Development
The F14 RIL population derived from the Upland cotton hybrid “Xinza 1” (GX 1135 × GX 100-2) was employed, which was developed by single seed descent method (Shang et al., 2016a). Two experimental populations were deduced based on the design: (1) backcross population, comprising 177 BCF1 progenies, was developed by testcrossing 177 F14 RILs with maternal parent GX1135, respectively. (2) RIL population, including 177 lines of F15 generations.
The control set (GX1135, “Xinza 1” F1, GX100-2, and a variety used as competition control) was performed in three field trials. The competition control in Yellow River Region (2015E1 and 2015E2, see details below) was a commercial hybrid “Ruiza 816” F1 by cross of 087 × 884, another in Yangtze River Region (2015E3, see detail below) was a commercial hybrid “Ezamian 10” F1 by cross of Tai 96167 × (GK19 chosen by Tai D-3).
Field Experiment and Trait Evaluation
Three field trials based on RIL population and maternal BC population were conducted in 2015E1, 2015E2, and 2015E3, E1: Quzhou Experimental Station of China Agricultural University, Handan City, Hebei Province (36°78′N, 114°92′E) in Yellow River valley; E2: Guoxin Seed Company Ltd, Cangzhou City, Hebei Province (38°43′N, 116°09′E) in Yellow River valley; E3: Hubei Academy of Agricultural Sciences, Wuhan City, Hubei Province (30°29′N, 114°19′E) in Yangtze valley. In addition, previous study in 2012 performed three same trials in E1 (2012E1) in Handan City, E2 (2012E2) in Cangzhou City, and E4 (2012E4) in Xiangyang City, Hubei Province (30°18′N, 112°15′E) (Shang et al., 2016a,d).
Seedlings were cultured on May 14th and May 15th (E1) and on April 19th, 2015 (E3). The transplanting was on June 12th (20–37°C) and June 13th (20–34°C), 2015 (E1), and on May 9th (17–28°C) and May 10th (15–29°C), 2015 (E3), respectively. Seeds were directly sowed on April 25th and 26th (17–31°C), 2015 (E2). The climate and weather in 2015 (E1, E2, and E3) was regular during cotton growth season, excepting that one trial encountered hailstone one June 11th, 2015 (E2). Since all the experiment materials were under the identical environment so we strengthened field managements in time to recover the seedling development after hailstone disaster (2015E2) so as to obtain the data sets. The flowering time had been postponed 6 weeks in that environment comparing with normal.
The field trials were planted following a randomized complete block design with two replications each trial, which included 904 plots with two rows per plot (22 individuals each). For one replication, 177 BCF1 progenies were inter-planted in the middle of 177 RI Lines as their maternal parents and the common testing-male parent GX1135 (original maternal parent), together with two control sets (GX1135, “Xinza 1” F1, GX100-2, control: “Ruiza816” F1 at E1 and E2, or “Ezamian 10” F1 at E3 and E4). Two-row plots were 90 and 50 cm row spacing alternately for the experiment in E1 and E2, and two-row plots in E3 were 100 and 50 cm row spacing alternately. One plot was on the length of 3.0, 3.0, and 2.4 m respectively at 2015E1, 2015E2, and 2015E3, following 0.7 m pavement apart. Field management followed the local conventional standard field practices.
Twenty-five naturally opening bolls in the middle of plants were hand-harvested for each plot at mature stage in three environments, respectively. Fiber samples were ginned and sampled for measurements of fiber quality traits with HVI 900 instrument (USTER_ HVISPECTRUM, SPINLAB, USA) at Cotton Fiber Quality Inspection and Test Center of Ministry of Agriculture (Anyang, China). The fiber quality traits included fiber length (mm), fiber uniformity (%), fiber strength (cN/tex), fiber elongation, and Micronaire. The total sample numbers were 892, 904, and 884 at 2015E1, 2015E2, and 2015E3, respectively.
Genetic Map and Data Analysis
The genetic map based on RIL population published before by screening polymorphic loci between two parents using 48,836 pairs of SSR primer (Shang et al., 2016d). The source of the SSR markers contained types of BNL, NAU, TM, JESPER, CIR, HAU, CM, MUSS, MUSB, MUCS, SWU, PGML, and CAU from Gossypium, which can be obtained from Cotton Microsatellite Database (http://www.cottonmarker.org) and CottonGen Database (https://www.cottongen.org) (Liang et al., 2013; Shang et al., 2016d). A total of 623 loci distributed on 31 linkage groups and anchored on 26 chromosomes, covering 3889.9 cM (88.20%) of cotton genome with average interval of 6.2 cM. The genotype for each maternal BCF1 was deduced on the basis of the RIL genotype used as the parent for testcross (Shang et al., 2016a,b,c).
Basic statistical analysis was implemented by the software SPSS (Version 19.0, SPSS, Chicago). The software QTL Cartographer (Version 2.5) was used to map single-locus QTL and to estimate the genetic effect (Zeng, 1994; Wang S. et al., 2012). The composite interval mapping (CIM) method was chosen for QTL mapping in the confidence interval of 95%. The precision selected 1.0 cM walk speed. The threshold of LOD 3.0 was used to declare a suggestive QTL after performing a 1,000 permutation test, whereas QTLs in another environment or population with at least LOD 2.5 were considered as common QTL (Liang et al., 2013; Wójcik-Jagła et al., 2013; Shang et al., 2015, 2016d). According to the position linked and shared common markers, QTLs detected in different populations were regarded as common QTL (Shao et al., 2014). The QTL IciMapping 4.1 (www.isbreeding.net) was conducted by the two-locus analysis using inclusive composite interval mapping (ICIM) method. A threshold LOD 2.5 and LOD 5 scores were used to declare significant M-QTL and E-QTLs, respectively. Mapping by IciMapping 4.1, the main-effect QTL (M-QTL) and its environmental interaction (QTL × environment, QE), epistatic QTLs (E-QTLs), and its environmental interactions (QTLs × environment, QQE) were conducted using RIL and BC datasets under three environments.
Trait Performance of Parents and Two Populations
The performances of “original” maternal parent “GX1135” were superior to those of “original” male parent “GX100-2” for fiber length (FL) in three locations, as well as for fiber uniformity (FU), fiber strength (FS), and fiber elongation (FE) in 2015E2 and 2015E3 (Table 1). On the contrary, performance of GX100-2 was superior to performance of GX1135 for Micronaire (FM) in three environments. Phenotype performance of five fiber quality traits varied in three environments such as that fiber quality in 2015E1 was always inferior to fiber quality in 2015E2 and 2015E3 (Table 1). It attributed to the diverse photo-thermal conditions, fertilizer application and growth regulators over the development stages in diverse environments. “Xinza 1” displayed −0.63% mid-parent heterosis (MPH, %) on average for five fiber quality traits (Table 1) and MPH (%) in BC population ranged from 0.08 to 2.91% among 177 BCF1S, in accordance with the previous study (Tables S1, S7) (Shang et al., 2016d). The result demonstrated that there was no significant hybrid vigor for fiber quality traits.
Table 1. Summary statistics on fiber quality trait data of RIL and BC populations in three environments.
In RIL and BC populations, genotype variance and environment variance for five traits showed significant at level of 0.01 except FU in BC population (Table 2). Then, coefficient of variation (CV) ranged from 2.37 to 8.35% for FL, FS and FM in both populations. And, lighter phenotypic variation presented with 0.78–1.63% of CV for FU and FE in both populations. For five fiber quality traits, the majority of mean values in BC population were bigger than mean values in RIL population (Table 2). It indicated that phenotype performance decided by the performance of both parents. Heritability decreased from 0.93 in RIL population to 0.81 in BC population for fiber length and similar tendency between both populations was observed for other four traits. The results indicated that environment affected stronger for phenotype in BC population than that in RIL population. Taken together, wider ranges, bigger CV, and higher heritability for five fiber quality traits were observed in RIL population than that in BC population (Tables 1, 2).
Correlation analysis on five fiber quality traits were performed using the mean phenotypic values both in RIL and BC populations in 2015E1, 2015E2, and 2015E3, respectively (Table 3). For fiber length, the significant and positive correlation existed with fiber strength and fiber elongation, but negatively with Micronaire in three locations in RIL and BC populations. Fiber strength displayed positive correlation with fiber elongation but negative correlation with Micronaire. The correlation showed significant and positive between fiber uniformity and fiber strength in RIL population, as well as between fiber uniformity and fiber elongation. The results showed similar tendency among five fiber quality trait in comparison with previous results (Liang et al., 2013; Shang et al., 2015, 2016d). In BC population, fiber uniformity positively correlated with fiber elongation.
Single Locus QTL Analysis
In three field trials, a total of 55 QTLs were detected in two corresponding populations at 2015E1, 2015E2, and 2015E3, with 5.00–34.00% of phenotypic variance explained (PVE) (Figure 1, Table 4). In our previous study, 71 QTLs were detected in the same maternal BC and same RIL population for fiber quality traits in 2012E1, 2012E2, and 2012E4 (Shang et al., 2016d). Here, 20 QTLs appearing both in 2015 in this study and in 2012 in previous study were showed in Table S2. These QTLs detected repeatedly in both years, were anchored on 21 chromosomes, respectively.
Figure 1. Locations of QTLs for fiber quality traits identified in RIL and BC populations. * and ** (# and ##), marker showed, respectively, segregation distortion significant at P = 0.05 and 0.01 levels; markers with * and ** skewed toward the GX1135 alleles, and markers with # and ## skewed toward the GX100-2 alleles. aFL, fiber length; FU, fiber uniformity; FS, fiber strength; FE, fiber elongation; FM, Micronaire.
For fiber length, six and 9 QTLs were identified in the BC and the RIL populations, respectively (Figure 1, Table 4). Of these, six QTLs showed positive additive effects originated from GX1135 alleles and three QTLs increased fiber length through offering alleles by GX100-2 in RIL population. Five QTLs (45.45%) were distributed on Chr 5. Four stable QTLs, qFL-Chr5-1, qFL-Chr5-2, qFL-Chr5-3, and qFL-Chr5-4, were simultaneously identified in the RIL and BC populations across multiple environments. The qFL-Chr5-2 explained 14.07–34.00% of phenotypic variance (PV, hereinafter same). Nine (82.82%) QTLs, located on Chr 2, Chr 5, Chr 9, and Chr 11, were verified the previous results in 2012 (Table S2). The qFL-Chr5-3 explained highly to 19.63% of PV under multiple environments in both populations, within a 0.8 cM interval flanked by SSR markers Gh388 and SWU17713. In recent study, one fine mapping QTL named qFL-Chr1 was most likely located in a 0.9 cM interval flanked by the SSR markers MUSS422 and CIR018 (Xu et al., 2017). For five stable QTLs detected in both populations in present study, the values of genetic effects reduced in the BC population, while homozygotes became heterozygotes in each BCF1 following backcrossing to GX1135. Ignoring the effect by environment, additive effects of homozygotes contributed more than dominant effects of heterozygotes for fiber length, such as that qFL-Chr5-2 increased 0.61 mm for fiber length in 2015E1 in RIL population whereas it was 0.28 mm in 2015E2 and 0.29 mm in 2012E1 in BC population (Table S2).
Nine QTLs for fiber uniformity were identified. Of those, three QTLs (qFU-Chr5-1, qFU-Chr5-2, qFU-Chr5-3) were anchored on Chr 5, explained 9.90, 14.53, and 7.62% of PV, respectively. One stable QTL (qFU-Chr11-1) was detected at 2015E1 in maternal BC population and at 2015E3 in RIL population.
For fiber strength, a total of 12 QTLs were detected (Table 4). Four QTLs were detected at least in two environments and in 2 years (Table S2). One stable QTL (qFS-Chr5-2) was observed in two environments and both populations, explaining 5.13, 6.00, and 29.65% of PV, respectively. On Chr 5, the flanking markers of qFS-Chr5-1 and qFS-Chr5-2 were same to that of qFL-Chr5-1 and qFL-Chr5-3, respectively, suggesting the pleiotropic effects on improving fiber quality in Upland cotton. One QTL near the two region increased 0.33 mm fiber length and 0.49 cN/tex fiber strength on average. The genetic effects of the QTLs from RI lines were stronger than that from the corresponding BCF1s, whatever QTLs were detected on Chr 5 or on Chr 21. It indicated that additive QTLs played more important role in improving fiber length and fiber strength.
Eight and four QTLs for fiber elongation were identified in the RIL and BC populations, respectively. Two QTLs (qFE-Chr5-1, qFE-Chr5-2) were confirmed by previous result in 2012 under multiple environments and populations. Five QTLs explained PV over 10% involving in qFE-Chr1-1, qFE-Chr4-2, qFE-Chr5-1, qFE-Chr5-2, and qFE-Chr17-1. The phenotypic variation explaining by qFE-Chr5-1 reached to 14.23% and qFE-Chr4-2 explained 12.41% of PV (Table 4). We also observed a pleiotropic region (NAU4034-SWU17713) on Chr 5. Two stable QTLs (qFE-Chr5-2, qFE-Chr24-1) from RIL population were verified in BC population.
A total of 12 QTLs underlying Micronaire were located on 9 chromosomes (Table 4), in which five QTLs were detected repeatedly in 2012 (Table S2), and seven QTLs distributed on At-subgenome and five on Dt-subgenome, respectively. Two stable QTLs (qFM-Chr14-1 and qFM-Chr14-2) were detected across two environments, explained 11.81 and 10.38% of PV on average, respectively. The qFM-Chr9-1 flanking with SWU15194-HAU190 was validated a common QTL explaining highly to 25.93% in 2012E4 in the previous result (Table S2).
On summary, 37 QTLs existed in RIL population with on average 9.43% of PVE; and 26 QTLs in BC population with average 7.92% of PVE (Figure 1, Table 4). The number of QTLs in the RIL population was always more than that in the BC population, which showed similar tendency with the previous conclusion in 2012 (Shang et al., 2016d). In present study, a total of 10 stable QTLs were detected across two populations or in at least two environments. Of these, eight QTLs were detected in both populations: qFL-Chr5-1, qFL-Chr5-2, qFL-Chr5-3, qFL-Chr5-4, qFU-Chr11-1, qFS-Chr5-2, qFE-Chr24-1, and qFM-Chr14-1. A co-location between Gh388 (SWU17717, NAU4034) and TMB1296 (SWU17713) contained four QTLs (qFL-Chr5-3, qFS-Chr5-2, qFE-Chr5-1, qFU-Chr5-2) controlling fiber length, fiber strength, fiber elongation, and fiber uniformity together. The pleiotropic region was consistent with the region of qFL-Chr5-4, qFS-Chr5-2, qFE-Chr5-1 in previous study (Shang et al., 2016d).
QTLs by Environment Interactions and Digenic Interactions
Totally, 102 M-QTLs and environmental interactions (QTL × environment, QE) were identified (Table 5, Tables S3–S6). Of these, 62 and 38 M-QTLs and QEs existed in the RIL population and its BC population, explained 0.13–5.95 and 0.12–7.61% of PV (Tables S3, S5), respectively. Environmental effect prevailed in both populations, but it affected M-QTLs in the BC population in a greater extent than that in RIL population. Then, 97 and 67 E-QTLs and their environmental interactions (digenic interactions × environment, QQE) were respectively identified in RIL and BC populations (Tables S4, S6). The mean number of E-QTL and QQE was 19 with 3.02% of mean PV in RIL population while it was 13 with 3.50% of mean PV in BC population (Table 5). The total number of QTL in both types (62 M-QTLs and 97 E-QTLs) was larger in RIL population than the number of QTL (38 M-QTLs and 67 E-QTL) in BC population. To our surprise, E-QTLs contributed 3.48 and 3.85% of average PV for four traits except fiber uniformity in RIL and BC populations, respectively; whereas M-QTLs explained average PV of 2.41 and 2.42% in RIL and BC populations, respectively. At the same time, QTLs (M-QTLs and E-QTLs) explained 3.04% of PV for four traits on average in both populations whereas environment explained 0.47% of PV on average (Table 5). On the contrary, fiber uniformity displayed vulnerable to environment with 2.89% of average PV.
Table 5. Summary of M-QTLs and E-QTLs by environment interaction controlling fiber quality traits by inclusive composite interval mapping in RIL and maternal BC population.
Here, 18 and seven main-effect QTLs (M-QTLs and QEs) participated in epistatic interactions (E-QTLs and QQE) in RIL and BC populations, respectively (Table 6, Tables S3–S6). Three types of epistasis were checked in RIL population: Type I, both loci were M-QTLs; Type II, either locus between two loci was M-QTL; Type III, both loci were no M-QTLs (Shang et al., 2016d). However, no epistatic QTL with Type I, 7 epistatic QTLs with Type II and 60 epistatic QTLs with Type III were observed in BC population (Table 6). Apparently, Type III was the most popular type of epistatic QTLs in both populations. These results indicated that the epistasis played vital role in improving fiber quality both in RIL and BC populations of Upland cotton. In previous study, epistatic QTLs with significant additive × additive effects were identified for fiber quality traits (Wang B. et al., 2017).
For fiber quality traits, fiber length (FL), fiber strength (FS), and fiber elongation (FE) displayed higher heritability over 0.90 than 0.78 for fiber uniformity (FU) and 0.84 for Micronaire (FM) in RIL population (Table 2). Based on present research and previous results (Shang et al., 2015, 2016c), the performance of “original” maternal parent “GX1135” was superior to those of “original” male parent “GX100-2” for FL, FS, and FE, which was accordance in the result of present study. Taken together, the experiment design was totally finished in six environments across 2 years and four locations (Table S7). Same competition control hybrid (CK) displayed 29.67 and 31.18 mm in 2012E2 and 2015E2, respectively, while “Xinza 1” showed 26.41 and 30.88 mm in 2012E2 and 2015E2, respectively. In addition, 26.96 and 26.53 mm of mean FL values were observed in BC and RIL populations in 2012E2, respectively, while they were 31.01 and 30.85 mm in 2015E2. The results indicated the data validity in present study because of similar variation in “Xinza 1,” the same competition control hybrid and both populations after encountering hailstone disaster in 2015E2.
The identification of stable QTLs across multiple environments and populations plays an essential role in marker-assisted selection (MAS) (Jamshed et al., 2016). In present study, 10 stable QTLs were simultaneously identified in more than one environment or population, of which six stable QTLs focused on Chr 5 (Table 4). Eight QTLs verified each other in the RIL population and its BC progenies as we mentioned above. The results suggested that it is reasonable and effective to map QTLs using the RIL population together with its corresponding BC population across multiple environments and multiple years, which trials were executed in a total of six environments (in present study: 2015E1, 2015E2, 2015E3 and in previous study: 2012E1, 2012E2, 2012E4) by the experiment design using RIL population and its BC population (Shang et al., 2016d).
Twenty QTLs (36.36%) even including seven (70.00%) stable QTLs identified in 2015 were validated the previous result in three environments in 2012 (Shang et al., 2016a; see Table 4, Table S2). Ten of the 20 QTLs were verified each other in two populations. The QTL regions were located on 8 chromosomes. Then, qFL-Chr5-1 was detected simultaneously across five environments and qFM-Chr26-1 across four environments including three environments in 2012 in previous study (Shang et al., 2016a). Two QTLs (qFE-Chr5-1, qFM-Chr9-1) explained 12.86 and 13.83% of phenotypic variations on average, respectively. In 2008 and 2009, three QTLs near MUSS193 or Gh388 markers co-locating on Chr 5, were also detected across multiple environments using F2:3, F2:4 segregating populations derived from the hybrid “Xinza 1” (Liang et al., 2013). The region was same to one QTL (qFL-Chr5-4) in present study, which closely linked with three other QTLs flanking with HAU1603 and SWU17781. Sun et al. (2012) found one cluster underlying MUSS193a harbored four QTLs in RIL population of Upland cotton. The cluster influenced fiber length, fiber elongation, fiber uniformity, and Micronaire. Markers of NAU6240, PGML1671, PGML1917, SWU17715, Gh388 (SWU17717, NAU4034), and SWU17713 (TMB1296) on Chr 5 across multiple environments, multiple populations, and multiple years will be recommended in MAS to improve fiber quality in Upland cotton. Markers flanking with these stable QTLs detected in multiple populations in present study are valuable for future breeding program and facilitate fine mapping, gene cloning, and favorable gene pyramiding project (Shao et al., 2014).
A total of 251 stable mapped QTLs for fiber quality were detected in the past decade in intraspecific populations derived from Upland cotton (Table S8). Due to the different parents and different genetic maps, few common markers and little common QTL for fiber quality traits in the present BCF1 population were found comparing with the previous researches in other labs in Table S8. The qFU-Chr26-2 controlling fiber uniformity in present study shared the common marker BNL2495 underlying qFM12.1 controlling Micronaire in previous study (Zhang et al., 2012). Nine markers flanking QTLs for fiber quality traits in Tang's work were common in our map, but no common QTL was observed in comparison with present study (Tang et al., 2015). The qFS-Chr21-2 explaining 9.06% of PV for fiber strength in present study shared a common SSR marker BNL3171 flanking a stable QTL of qUHM-21-1 controlling fiber length in BC3F2 and BC3F2:4 generation in previous study (Wang B. et al., 2017). Another QTL qFS-Chr1-1 controlling fiber strength in present study also shared SWU10954 which co-segregated with NAU4045 flanking a stable qUHM-15-1 controlling fiber length in previous study (Wang B. et al., 2017). In addition, some QTLs in present study located on same chromosomes in accordance with previous studies, such as 13 QTLs locating on Chr 5 in present study, which offered clues and new information for further exploitation. Yang et al. (2015) detected two clusters on Chr 5 and two QTLs on Chr 5 related to fiber uniformity presented in two environments; Zhang et al. (2005) detected one cluster on Chr 5 controlling LP1 for lint percentage, FS2 for fiber strength, and FE1 for fiber elongation. Wang et al. (2016) identified 6 QTLs including two stable QTLs for fiber length distributing on Chr 2, Chr 5, Chr 10, and Chr 19, which were in the shared linkages / chromosomes comparing with present study. Said et al. (2013) summarized 35 QTLs for fiber quality traits (and six QTLs for yield and its components traits) distributed on Chr 5, and 50 QTLs related fiber quality traits distributed on Chr 14. Wang B. et al. (2017) identified the co-location region of BNL2656 and NAU3498 on Chr 5 containing qUHM-5-1 and qUI-5-1, which improved fiber length and fiber uniformity at the same time. In present study, 11 in 12 QTLs on Chr 5 improved fiber length, fiber uniformity, fiber strength, and fiber elongation (Table S2). Further researches need to be performed to interpret the relationship of fiber related genes with stable QTLs in Upland cotton. Recently, the cotton fiber homolog GhCFE5 on Chr A01 and Chr D05 was found by mapping GhCFE5A and GhCFE5D in BC1 interspecific mapping population, as the same chromosomes in present study (Lv et al., 2016).
The genetic effects were explored in RIL population and its maternal BC progenies. Judging from the overall phenotypic level, RI lines with wider range and maximum values were not obviously worse than some BCF1S, whereas a number of lines in RIL population have better phenotypes than that in BC population (Table 1, Table S9). At the single-locus level, all of six QTLs detected in BC population increased 0.24 mm fiber length on average (Table 4, Table S2). It indicated that the heterozygotes increased fiber length a little than the respective homozygotes (Xiao et al., 1995). For fiber length, qFL-Chr5-3 on Chr 5 increased 0.25 mm in BC population and 0.36 mm in RIL population in 2015E2. At the same time, the QTL increased 0.19 mm in BC population and 0.36 mm in RIL population in 2012E1. The result indicated that the genetic value contributing by one common QTL decreased in BC population in comparison with the value in RIL population.
The epistatic effects and environmental interactions existed commonly for fiber quality traits. Fiber quality traits were sensitive to environment effect (Shao et al., 2014). At two-locus level, there were a number of interactions under environments for fiber quality traits in both populations, but epistasis influenced the performance stronger on average in the maternal BC population than that in RIL population (Table 5). Three types of epistasis combinations were observed (Table 6). Totally, 91 and 81% epistatic interaction of Type III accounted in RIL population and in BC1 population, respectively. The result indicated that the Type III was the common epistasis interaction type in both populations. Another interesting result detected that epistasis is another vital genetic effect affecting fiber quality traits following the previous study (Shang et al., 2016d). Similarly to the previous study, Wang et al. (2006) indicated that both epistasis effect and single-locus effect of QTLs played important genetic roles in cotton fiber quality.
There was no evidence that fiber quality traits had heterosis in “Xinza 1” (Table 1). We also found that fiber quality of BCF1S was not apparently superior to its parents in BC population (Table S1). The maximum of fiber length was 34.05 mm in 2015E3 in RIL population but 33.20 mm in BC population, in agreement with that heterozygotes were not always advantageous for performance (Xiao et al., 1995; Hua et al., 2002). For many crops, such as rice and other cereal crops, yield heterosis was obvious so that people have been interested to exploit and interpret the mechanism of hybrid vigor. Differing from these crops, the fiber quality traits of cotton should have different genetic mechanisms in development stage. In nature, plant species needed to produce more seeds for survival in evolution. Obviously, fiber is not necessary for optimal survival. Maybe this is one of the reasons why negative correlation exists between yield and fiber quality (Yang et al., 2015). Further research will be deduced on both heterosis and development of fiber in Upland cotton in order to improve fiber quality and to ultimate heterosis in breeding program.
LM, YZ, and LS performed the experiments, analyzed the data, and prepared the manuscript. YW maintained the experimental platform and attended bench work. JH conceived the experiments, provided experimental platform, and revised the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We thank Xianda Yi (Hubei Academy of Agricultural Sciences, Hubei Province), Xiaofang Yin, Dongyong Xu, and Huaiyu Lu (Guoxin Seed Company Ltd., Cangzhou, Hebei Province) for their contributions on field experiments and data acquisition. Thanks to Xiaocui Wang (China Agricultural University) for the work in SSR marker evaluation, Kunbo Wang and Fang Liu (Institute of Cotton Research, Chinese Academy of Agricultural Sciences) for providing part of SSR markers, Dr. Zhengsheng Zhang (Southwest University) for providing SWU and ICR SSR primers. This research was supported by the National Key R & D Program for Crop Breeding (2016YFD0101407) to JH.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2017.02168/full#supplementary-material
BC, Backcross population; Chr, Chromosome; CIM, Composite interval mapping; FE, Fiber elongation; FL, Fiber length; FM, Fiber micronaire; FS, Fiber strength; FU, Fiber uniformity; ICIM, Inclusive composite interval mapping; LOD, Likelihood of odd score; QTL, Quantitative trait loci; PV(E), Phenotypic variance (explained); RIL, Recombinant inbred line population.
Fang, L., Wang, Q., Hu, Y., Jia, Y. H., Chen, J. D., Liu, B. L., et al. (2017). Genomic analyses in cotton identify signatures of selection and loci associated with fiber quality and yield traits. Nat. Genet. 19, 1089–1098. doi: 10.1038/ng.3887
Hua, J., Xing, Y., Wu, W., Xu, C., Sun, X., Yu, S., et al. (2003). Single locus heterotic effects and dominance by dominance interactions can adequately explain the genetic basis of heterosis in an elite rice hybrid. Proc. Natl. Acad. Sci. U.S.A. 100, 2574–2579. doi: 10.1073/pnas.0437907100
Hua, J. P., Xing, Y. Z., Xu, C. G., Sun, X. L., Yu, S. B., and Zhang, Q. F. (2002). Genetic dissection of an elite rice hybrid revealed that heterozygotes are not always advantageous for performance. Genetics 162, 1885–1895. Available online at: http://www.genetics.org/content/162/4/1885.article-info
Hulse-kemp, A. M., Lemm, J., Plieske, J., Ashrafi, H., Buyyarapu, R., Fang, D. D., et al. (2015). Development of a 63K SNP array for cotton and high-density mapping of intraspecific and interspecific populations of Gossypium spp. G3 5:1187. doi: 10.1534/g3.115.018416
Islam, M. S., Zeng, L. H., Delhom, C. D., Song, X. L., Kim, H. J., Li, P., et al. (2014). Identification of cotton fiber quality quantitative trait loci using intraspecific crosses derived from two near-isogenic lines differing in fiber bundle strength. Mol. Breeding 34, 373–384. doi: 10.1007/s11032-014-0040-4
Jamshed, M., Jia, F., Gong, J., Palanga, K. K., Shi, Y., Li, J., et al. (2016). Identification of stable quantitative trait loci (QTLs) for fiber quality traits across multiple environments in Gossypium hirsutum recombinant inbred line population. BMC Genomics 17:197. doi: 10.1186/s12864-016-2560-2
Lacape, J. M., Nguyen, T. B., Courtois, B., Belot, J. L., Giband, M., Gourlot, J. P., et al. (2005). QTL analysis of cotton fiber quality using multiple Gossypium hirsutum × Gossypium barbadense backcross generations. Crop Sci. 45, 123–140. doi: 10.2135/cropsci2005.0123a
Larièpe, A., Mangin, B., Jasson, S., Combes, V., Dumas, F., Jamin, P., et al. (2012). The genetic basis of heterosis: multiparental quantitative trait loci mapping reveals contrasted levels of apparent overdominance among traits of agronomical interest in maize (Zea mays L.). Genetics 190, 795–811. doi: 10.1534/genetics.111.133447
Li, C., Dong, Y. T., Zhao, T. L., Li, L., Li, C., Yu, E., et al. (2016). Genome-wide SNP linkage mapping and QTL analysis for fiber quality and yield traits in the Upland cotton recombinant inbred lines population. Front Plant Sci. 7:1356. doi: 10.3389/fpls.2016.01356
Li, F., Fan, G. Y., Lu, C., Xiao, G., Zou, C., Kohel, R. J., et al. (2015). Genome sequence of cultivated Upland cotton (Gossypium hirsutum TM-1) provides insights into genome evolution. Nat. Biotechnol. 33, 524–530. doi: 10.1038/nbt.3208
Liang, Q. Z., Hu, C., Hua, H., Li, Z. H., and Hua, J. P. (2013). Construction of a linkage map and QTL mapping for fiber quality traits in upland cotton (Gossypium hirsutum L.). Chinese Sci. Bull. 58, 3233–3243. doi: 10.1007/s11434-013-5807-1
Liang, Q., Shang, L., Wang, Y., and Hua, J. (2015). Partial dominance, overdominance and epistasis as the genetic basis of heterosis in Upland Cotton. PLoS ONE 10:e0143548. doi: 10.1371/journal.pone.0143548
Liu, X., Zhao, B., Zheng, H. J., Hu, Y., Lu, G., Yang, C. Q., et al. (2015). Gossypium barbadense genome sequence provides insight into the evolution of extra-long staple fiber and specialized metabolites. Sci. Rep. 5:14139. doi: 10.1038/srep14139
Lv, F., Li, P., Zhang, R., Li, N., and Guo, W. (2016). Functional divergence of GhCFE5, homoeologs revealed in cotton fiber and Arabidopsis, root cell development. Plant Cell Rep. 35, 867–881. doi: 10.1007/s00299-015-1928-0
Mei, H. W., Li, Z. K., Shu, Q. Y., Guo, L. B., Wang, Y. P., Yu, X. Q., et al. (2005). Gene actions of QTLs affecting several agronomic traits resolved in a recombinant inbred rice population and two backcross populations. Theor. Appl. Genet. 110, 649–659. doi: 10.1007/s00122-004-1890-7
Ning, Z. Y., Chen, H., Mei, H. X., and Zhang, T. Z. (2014). Molecular tagging of QTLs for fiber quality and yield in the upland cotton cultivar Acala-Prema. Euphytica 195, 143–156. doi: 10.1007/s10681-013-0990-3
Paterson, A. H., Wendel, J. F., Gundlach, H., Guo, H., Jenkins, J., Jin, D., et al. (2012). Repeated polyploidization of Gossypium genomes and the evolution of spinnable cotton fibres. Nature 492, 423–427. doi: 10.1038/nature11798
Said, J. I., Knapka, J. A., Song, M. Z., and Zhang, J. F. (2015). Cotton QTLdb: a cotton QTL database for QTL analysis, visualization, and comparison between Gossypium hirsutum and G. hirsutum × G. barbadense populations. Mol. Genet. Genomics 290, 1615–1625. doi: 10.1007/s00438-015-1021-y
Said, J. I., Lin, Z., Zhang, X., Song, M., and Zhang, J. (2013). A comprehensive meta QTL analysis for fiber quality, yield, yield related and morphological traits, drought tolerance, and disease resistance in tetraploid cotton. BMC Genomics 14:776. doi: 10.1186/1471-2164-14-776
Shang, L. G., Liang, Q. Z., Wang, Y. M., Wang, X. C., Wang, K. B., Abduweli, A., et al. (2015). Identification of stable QTLs controlling fiber traits properties in multi-environment using recombinant inbred lines in Upland cotton (Gossypium hirsutum L.). Euphytica 205, 877–888. doi: 10.1007/s10681-015-1434-z
Shang, L., Liang, Q., Wang, Y., Zhao, Y., Wang, K., and Hua, J. (2016a). Epistasis together with partial dominance, over-dominance and QTL by environment interactions contribute to yield heterosis in upland cotton. Theor. Appl. Genet. 129, 1429–1446. doi: 10.1007/s00122-016-2714-2
Shang, L., Wang, Y., Cai, S., Ma, L., Liu, F., Chen, Z., et al. (2016b). Genetic analysis of Upland cotton dynamic heterosis for boll number per plant at multiple developmental stages. Sci. Rep. 6:35515. doi: 10.1038/srep35515
Shang, L., Wang, Y., Cai, S., Wang, X., Li, Y., Abduweli, A., et al. (2016c). Partial dominance, overdominance, epistasis and QTL by environment interactions contribute to the heterosis in two Upland cotton hybrids. G3 6, 499–507. doi: 10.1534/g3.115.025809
Shang, L., Wang, Y. M., Wang, X. C., Liu, F., Abduweli, A., Cai, S. H., et al. (2016d). Genetic analysis and stable QTL detection on fiber quality traits using two recombinant inbred line populations and their backcross progeny in Upland cotton. G3 6, 2717–2724. doi: 10.1534/g3.116.031302
Shao, Q. S., Zhang, F. J., Tang, S. Y., Liu, Y., Fang, X. M., Liu, D. X., et al. (2014). Identifying QTL for fiber quality traits with three upland cotton (Gossypium hirsutum L.) populations. Euphytica 198, 43–58. doi: 10.1007/s10681-014-1082-8
Shen, X. L., Guo, W. Z., Lu, Q. G., Zhu, X. F., Yuan, Y. L., and Zhang, T. Z. (2007). Genetic mapping of quantitative trait loci for fiber quality and yield trait by RIL approach in Upland cotton. Euphytica 155, 371–380. doi: 10.1007/s10681-006-9338-6
Su, J., Li, L., Pang, C., Wei, H., Wang, C., Song, M., et al. (2016). Two genomic regions associated with fiber quality traits in Chinese upland cotton under apparent breeding selection. Sci. Rep. 6:38496. doi: 10.1038/srep38496
Sun, F. D., Zhang, J. H., Wang, S. F., Gong, W. K., Shi, Y. Z., Liu, A. Y., et al. (2012). QTL mapping for fiber quality traits across multiple generations and environments in upland cotton. Mol. Breeding 30, 569–582. doi: 10.1007/s11032-011-9645-z
Tan, Z. Y., Fang, X. M., Tang, S. Y., Zhang, J., Liu, D. J., Teng, Z. H., et al. (2014). Genetic map and QTL controlling fiber quality traits in upland cotton (Gossypium hirsutum L.). Euphytica 203, 615–628. doi: 10.1007/s10681-014-1288-9
Tang, S. Y., Teng, Z. H., Zhai, T. F., Fang, X. M., Liu, F., Liu, D. J., et al. (2015). Construction of genetic map and QTL analysis of fiber quality traits for Upland cotton (Gossypium hirsutum L.). Euphytica 201, 195–213. doi: 10.1007/s10681-014-1189-y
Wan, Q., Zhang, Z. S., Hu, M. C., Chen, L., Liu, D. J., Chen, X. C., et al. (2007). T1 locus in cotton is the candidate gene affecting lint percentage, fiber quality and spiny bollworm (Earias spp.) resistance. Euphytica 158, 241–247. doi: 10.1007/s10681-007-9446-y
Wang, B., Draye, X., Zhuang, Z., Zhang, Z., Liu, M., Lubbers, E., et al. (2017). QTL analysis of cotton fiber length in advanced backcross populations derived from a cross between Gossypium hirsutum and G. mustelinum. Theor. Appl. Genet. 130, 1297–1308. doi: 10.1007/s00122-017-2889-1
Wang, B. H., Guo, W. Z., Zhu, X. F., Wu, Y. T., Huang, N. T., and Zhang, T. Z., et al. (2006). QTL mapping of fiber quality in an elite hybrid derived-RIL population of upland cotton. Euphytica 152, 367–378. doi: 10.1007/s10681-006-9224-2
Wang, H. T., Huang, C., Zhao, W. X., Dai, B. S., Shen, C., Zhang, B. B., et al. (2016). Identification of QTL for fiber quality and yield trait s using two immortalized backcross populations in Upland cotton. PLoS ONE 11:e0166970. doi: 10.1371/journal.pone.0166970
Wang, M., Tu, L., Lin, M., Lin, Z., Wang, P., Yang, Q., et al. (2017). Asymmetric subgenome selection and cis-regulatory divergence during cotton domestication. Nat. Genet. 49, 579–587. doi: 10.1038/ng.3807
Wang, S., Basten, C. J., and Zeng, Z. B. (2012). Windows QTL Cartographer 2.5. Department of Statistics; North Carolina State University, Raleigh, NC. Available online at: http://statgen.ncsu.edu/qtlcart/WQTLCart.htm
Wójcik-Jagła, M., Rapacz, M., Tyrka, M., Kościelniak, J., Crissy, K., and Zmuda, K. (2013). Comparative QTL analysis of early short-time drought tolerance in polish fodder and malting spring barleys. Theor. Appl. Genet. 126, 3021–3034. doi: 10.1007/s00122-013-2190-x
Wu, J. X., Gutierrez, O. A., Jenkins, J. N., McCarty, J. C., and Zhu, J. (2009). Quantitative analysis and QTL mapping for agronomic and fiber traits in an RI population of Upland cotton. Euphytica 165, 231–245. doi: 10.1007/s10681-008-9748-8
Xu, P., Gao, J., Cao, Z. B., Chee, P. W., Guo, Q., Xu, Z. Z., et al. (2017). Fine mapping and candidate gene analysis of qFL-chr1, a fiber length QTL in cotton. Theor. Appl. Genet. 130, 1309–1319. doi: 10.1007/s00122-017-2890-8
Yang, X. L., Zhou, X. D., Wang, X. F., Li, Z. K., Zhang, Y., Liu, H. G., et al. (2015). Mapping QTL for cotton fiber quality traits using simple sequence repeat markers, conserved intron-scanning primers, and transcript-derived fragments. Euphytica 201, 215–230. doi: 10.1007/s10681-014-1194-1
Yuan, D., Tang, Z., Wang, M., Gao, W., Tu, L., Jin, X., et al. (2015). The genome sequence of Sea-Island cotton (Gossypium barbadense) provides insights into the allopolyploidization and development of superior spinnable fibres. Sci. Rep. 5:17662. doi: 10.1038/srep17662
Zhang, J., Yu, J., Pei, W., Li, X., Said, J., Song, M., et al. (2015). Genetic analysis of Verticillium wilt resistance in a backcross inbred line population and a meta-analysis of quantitative trait loci for disease resistance in cotton. BMC Genomics 16:577. doi: 10.1186/s12864-015-1682-2
Zhang, K., Zhang, J., Ma, J., Tang, S., Liu, D., Teng, Z., et al. (2012). Genetic mapping and quantitative trait locus analysis of fiber quality traits using a three-parent composite population in upland cotton (Gossypium hirsutum L.). Mol. Breed. 29, 335–348. doi: 10.1007/s11032-011-9549-y
Zhang, T., Hu, Y., Jiang, W., Fang, L., Guan, X., Chen, J., et al. (2015). Sequencing of allotetraploid cotton (Gossypium hirsutum L. acc. TM-1) provides a resource for fiber improvement. Nat. Biotechnol. 33, 531–537. doi: 10.1038/nbt.3207
Zhang, Z. S., Xiao, Y. H., Luo, M., Li, X. B., Luo, X. Y., Hou, L., et al. (2005). Construction of a genetic linkage map and QTL analysis of fiber-related traits in upland cotton (Gossypium hirsutum L.). Euphytica 144, 91–99. doi: 10.1007/s10681-005-4629-x
Zhang, Z., Li, J. W., Muhammad, J., Cai, J., Jia, F., Shi, Y. Z., et al. (2015). High resolution consensus mapping of quantitative trait loci for fiber strength, length and Micronaire on chromosome 25 of the Upland cotton (Gossypium hirsutum L.). PLoS ONE 10:e0135430. doi: 10.1371/journal.pone.0135430
Keywords: fiber quality, stable QTL, backcross population, Upland cotton, epistasis detection
Citation: Ma L, Zhao Y, Wang Y, Shang L and Hua J (2017) QTLs Analysis and Validation for Fiber Quality Traits Using Maternal Backcross Population in Upland Cotton. Front. Plant Sci. 8:2168. doi: 10.3389/fpls.2017.02168
Received: 06 August 2017; Accepted: 11 December 2017;
Published: 22 December 2017.
Edited by:Sergio J. Ochatt, UMR 1347 Agroécologie (INRA), France
Reviewed by:Aladdin Hamwieh, International Center for Agricultural Research in the Dry Areas, Egypt
Hamid Khazaei, University of Saskatchewan, Canada
Jinjin Jiang, Yangzhou University, China
Copyright © 2017 Ma, Zhao, Wang, Shang and Hua. 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) or licensor 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: Jinping Hua, email@example.com