Dan Liu1†Yanjiao Cui2†Jihong He3†Suying Li2†*Qiang Li2Dan Liang1Jianhe Wang1Xiaowei Shi1Conglei Wang1Kongjun Dong3Tianpeng Liu3Lei Zhang3Ruiyu Ren3Tianyu Yang3*Gang Feng1*
Zhengli Liu2*- 1Tianjin Crop Research Institute, Tianjin Academy of Agricultural Sciences, Tianjin, China
- 2Department of Life Sciences, Tangshan Normal University, Tangshan, China
- 3Crop Research Institute, Gansu Academy of Agricultural Sciences, Lanzhou, China
Due to the maternal inheritance of cytoplasm, using foxtail millet [Setaria italica (L.) P. Beauv.] male sterile lines with a single cytoplasmic source as the female parent will inevitably lead to a narrow source of cytoplasm in hybrids, which may make them vulnerable to infection by cytoplasm-specific pathogens, ultimately leading to destructive yield losses. To assess cytoplasmic genetic diversity in plants, molecular markers derived from chloroplast DNA (cpDNA) have been used. However, such markers have not yet been applied to foxtail millet. In this study, we designed and screened nine pairs of polymorphic foxtail millet-specific primers based on its completely sequenced cpDNA. Using these primers, we analyzed the genetic diversity and cytoplasmic types of 130 elite foxtail millet parental lines collected in China. Our results revealed that the cytoplasmic genetic diversity of these accessions was low and needs to be increased. The parental lines were divided into four cytoplasmic types according to population structure analysis and a female parent-derivative evolutionary graph, indicating that the cytoplasmic types of elite foxtail millet lines were rather limited. A principal component analysis (PCA) plot was linked with the geographic and ecological distribution of accessions for each cytoplasmic type, as well as their basal maternal parents. Collectively, our results suggest that enriching cytoplasmic sources through the use of accessions from diverse ecological regions and other countries as the female parent may improve foxtail millet breeding programs, and prevent infection by cytoplasm-specific pathogens.
Introduction
Heterosis refers to the superior performance of hybrids compared with their parental lines and has been extensively exploited to increase grain yield in many crops for more than a century, including maize (Zea mays L.), sorghum (Sorghum bicolor L.), and rice (Oryza sativa L.) (Guo et al., 2006; Springer and Stupar, 2007; Baranwal et al., 2012). Foxtail millet [Setaria italica (L.) P. Beauv.] originated in China, and has been a traditional cereal food crop in China since ancient times. Recently, researchers have developed foxtail millet hybrid cultivars based on the heterosis theory and successfully increased yields (Huang et al., 2018; Li et al., 2018). During the development of foxtail millet hybrids, a male sterile line and a restorer line are used. China is the center of diversity for foxtail millet landraces, and to date, 27,059 accessions have been recorded in the Chinese National GeneBank (CNGB) (Diao and Jia, 2017). Despite this, male sterile lines are represented by limited genetic diversity. Pedigree analysis indicates that almost all Chinese spring foxtail millet male sterile lines are derived from Chang10A, with cytoplasm from Qinyuanmujizui (Wang et al., 1998; Liu et al., 2011). Similarly, most Chinese summer foxtail millet male sterile lines are derived from Huangmi1A, which contains cytoplasm from Dahuanggu (Liu et al., 1996; Liu et al., 2006).
The wide use of a single source of cytoplasm presents considerable risks to agricultural production. For example, in the 1970s in the United States, the widespread use of maize Texas male-sterile cytoplasm (cms-T) rendered hybrid cultivars susceptible to the fungal pathogen Helminthosporium maydis, which caused an outbreak of southern leaf blight and resulted in a serious yield loss (Hooker et al., 1970; Scheifele et al., 1970; Smith et al., 1970). In 1984, the widespread use of Xianyou 2 in China, a wild abortion type rice male sterile cytoplasm, resulted in a major epidemic of Pyricularia oryzae, representing one of the most devastating rice fungal diseases (Liu et al., 1992; Yuan, 1997). Currently, the primary male sterile cytoplasms of foxtail millet cultivars in China are Qinyuanmujizui and Dahuanggu. Due to the maternal inheritance of cytoplasm, hybrid cultivars produced by using these male sterile lines still carry the same source of cytoplasm. Long-term and large-scale use of a single source of cytoplasm in foxtail millet male sterile lines may eventually cause infections by certain pathogens and ultimately lead to destructive losses in foxtail millet hybrid production. Thus, it is extremely important to enrich the cytoplasm resources in male sterile lines, in particular through a better understanding of the genetic diversity and classification of Chinese foxtail millet cytoplasms.
In recent years, molecular markers have been used to assess the genetic diversity of foxtail millet. To investigate the genetic diversity and population structure, Liu et al. (2011) first screened a collection of 128 foxtail millet germplasm and elite breeding lines from three major ecological areas of China for foxtail millet growth using 79 SSR markers. Subsequently, researchers conducted population diversity and structure analyses of Chinese foxtail millet landraces and cultivars using microsatellite markers (Wang et al., 2012a; He et al., 2015; Jia et al., 2015). However, these studies mainly focused to nuclear genomes. As such, there is a scarcity of reports on the genetic diversity of foxtail millet cytoplasm.
Both types of cytoplasmic genomes, mitochondrial DNA (mtDNA) and chloroplast DNA (cpDNA), represent important genetic material in addition to the nuclear genome of plants. Their distinct features, including small genome sizes, high conservation, and the absence of sexual recombination make mtDNA and cpDNA useful tools for determining phylogenetic relationships and studying plant populations (Avise, 2000; Provan et al., 2001; Norouzi et al., 2012; Cheng et al., 2015; Skuza et al., 2019). Recently, molecular markers derived from cpDNA and mtDNA have been designed and used for research on cytoplasmic variation in soybean (Glycine max) (Shimamoto et al., 2000), rice (Rajendrakumar et al., 2007), Lolium species (McGrath et al., 2007), and Brassica species (Wang et al., 2012b; Zamani-Nour et al., 2013; Shu et al., 2016). However, genetic diversity analysis and classification of cytoplasm types in foxtail millet accessions with molecular markers based on cpDNA or mtDNA have rarely been reported. In a previous study, we analyzed the genetic diversity of the cytoplasm of 111 elite foxtail millet breeding lines using 23 mtDNA consensus primer pairs, and constructed a phylogenetic tree to classify these accessions (Liu et al., 2014). However, the genetic diversity analysis and classification of elite parental lines have not yet been investigated, but is crucial for creating male sterile lines that carry new cytoplasm types.
In this study, the cytoplasmic genetic diversity of 130 elite foxtail millet parental lines was assessed based on cpDNA polymorphism. We used nine foxtail millet-specific primer pairs, which were developed based on the completely sequenced chloroplast genome. To classify the cytoplasmic types of these elite parental lines, we used population structure analysis and female parent-derivative evolutionary graphing. We aimed to understand the relationship between the geographic and ecological distribution of the tested accessions and their basal maternal sources by using principal component analysis (PCA). Our results provide insight for enriching the cytoplasmic types of foxtail millet hybrids in China and around the world, and also provide a reference to broaden the cytoplasmic types of conventional cultivars and hybrids of other crops.
Materials and Methods
Plant Materials
A total of 130 representative Chinese elite foxtail millet parental lines were collected. These materials were sourced from diverse geographic origins that spanned three ecotypes (north summer millet, northwest spring millet, and northeast spring millet). These accessions were primarily from eight foxtail millet growing provinces and cities (1 accession in Beijing, 107 accessions in Hebei, 3 accessions in Henan, 12 accessions in Shanxi, 1 accession in Shandong, 2 accessions in Jilin, 3 accessions in Liaoning, and 1 accession in Shaanxi), where the majority of Chinese foxtail millet breeding programs are located. Among the 130 elite parental lines, 51 accessions were developed and released after 2014. Information on the collected samples is summarized in Supplemental Table S1.
Deoxyribonucleic Acid Extraction, Primer Definition, and Selection
The cpDNA of four foxtail millet accessions (Gu56A, Gu572A, Datong28, and Datong29lv) was extracted and sequenced as described in a previous study (unpublished data, the chloroplast genome sequences have been deposited to GenBank). Genomic DNA used for DNA polymorphism detection was extracted using the DNAquick Plant System DP321 (Tiangen Biotech, Beijing, China) following the manufacturer’s instructions.
Single nucleotide polymorphisms (SNPs) and insertion-deletion polymorphisms (InDels) in four cpDNAs (Gu56A, Gu572A, Datong28, and Datong29lv) were detected by comparing the aligned sequences in DNAMAN software to identify those conducive for genotyping. Twenty-seven primer pairs were designed to amplify these regions and were selected to screen the genomes of 130 foxtail millet accessions. A total of nine primers with high levels of polymorphisms were used for further analysis. Primer sequences, annealing temperatures for PCR amplification, and amplification product sizes are summarized in Table 1.
Table 1 Information on the nine chloroplast markers used in this study.
A set of 42 universal cpDNA primers (https://www.ncbi.nlm.nih.gov/probe/?term=chloroplast+ssr) and 34 mtDNA consensus primers (Duminil et al., 2002) were selected to amplify the target sequences of genomes from the 130 selected foxtail millet accessions. Ten cpDNA primers and 17 mtDNA primers were used in this study according to their high levels of polymorphisms across these accessions. Primer sequences, annealing temperatures for PCR amplification, and amplification product sizes are summarized in Supplemental Tables S2 and S3.
Polymerase Chain Reaction Amplification and Gel Electrophoresis
PCR amplification was conducted in 10 µl reaction mixtures, which contained 20 ng template DNA, 10 pmol of each forward and reverse primers, and the 2× Taq PCR Master Mix (Tiangen Biotech, Beijing, China). The reaction was carried out with a thermocycler (Biometra, Göttingen, Germany) using the following conditions: an initial denaturation step of 95°C for 5 min, followed by 35 cycles of 95°C for 30 s, an appropriate annealing temperature of different primers (Table 1 and Supplemental Tables S2 and S3) for 30 s, 72°C for 2 min, and a final elongation step of 72°C for 10 min.
Amplified products were separated on a 1.5% agarose gel in 1×TAE buffer along with the DNA ladder DL2000 (Sangon, Shanghai, China) as a size marker. Next, the PCR products were visualized and photographed under a Gel Dox XR+ (Bio-Rad, Madison, CA, USA) photography system using ultraviolet light. Data were scored manually for band presence or absence and entered into a matrix for further analysis.
Data Analysis
Summary statistics, including the number of alleles per locus, major allele frequency, gene diversity, and polymorphism information content (PIC) values were determined using PowerMarker v3.25 (Liu and Muse, 2005). To detect population genetic structure and assign individuals to subpopulations, the program STRUCTURE v2.3.4 (Pritchard et al., 2000; Falush et al., 2003) was used, which employs a Bayesian clustering approach. Ten independent runs for cluster (K) values, which ranged from 2 to 10, were performed after a burn-in period of 5 × 105 Markov Chain Monte Carlo steps, followed by 1 × 105 replicates using the admixture model. The output was exported into Structure Harvester (Earl and vonHoldt, 2012) to determine the most likely number of K clusters (K = 4 was optimum for this analysis) using Evanno’s ΔK method (Evanno et al., 2005). Results from 10 independent STRUCTURE runs for the most likely K were assessed with the software CLUMPP (Jakobsson and Rosenberg, 2007) and plotted using the DISTRUCT program (Rosenberg, 2010).
Arlequin 3.11 (Excoffier et al., 2005) was used to conduct the analysis of molecular variance (AMOVA) and to calculate the pairwise FST values. A PCA was performed with SPSS 17.0 for Windows (SPSS Inc., Chicago, USA) using standard statistical methods.
Results
Designing of Chloroplast Deoxyribonucleic Acid Primer Pairs
Prior to this study, cpDNAs of four different foxtail millet accessions (Gu56A, Gu572A, Datong28, and Datong29lv) were extracted and sequenced. In this study, 27 primer pairs were designed and used for PCR amplification using the total DNA from 130 accessions. Primer pairs were designed based on the SNPs and insertion-deletion polymorphisms (InDels) in the foxtail millet cpDNAs. Nine pairs showed polymorphisms and were selected according to preliminary screening, with the percentage of polymorphism accounting for 33.33%. These primer pairs were all located in the large single copy (LSC) region of the foxtail millet chloroplast genome, and amplified coding (exon 44%) and mixed regions [exon + intron 12%, Intergenic Spacers (IGS) + gene 44%].
After genotyping the 130 foxtail millet accessions using the nine chloroplast DNA markers, a total of 20 alleles were detected. There was an average of two alleles per locus and the average major allele frequency was 0.70. The PIC values ranged from 0.03 to 0.61, with an average of 0.30. The mean gene diversity was 0.37 (Table 2). These results indicated that the cytoplasmic genetic diversity of foxtail millet accessions assessed here was low.
Table 2 Data summary for 130 foxtail millet accessions based on nine chloroplast DNA markers.
Population Structure Analysis
An analysis of the population structure of the 130 foxtail millet accessions showed that the most appropriate grouping was four subpopulations with a ΔK peak (K = 4 using the program STRUCTURE; Figure 1A). Thus, the group of foxtail millet accessions were divided into four subpopulations termed G1, G2, G3, and G4 (Figure 1B). Among the four subpopulations, the level of genetic diversity within G2 was the highest (0.28), followed by G1 (0.26), G4 (0.25), and G3 (0.23). Additionally, G2 also had the most alleles scored (Table 3).
Figure 1 Population structure analysis for 130 accessions of foxtail millet. (A) Delta K values for different numbers of populations assumed (K) in the structure analysis. (B) Classification of 130 accessions into four subpopulations according to preset K value using STRUCTURE program. The distribution of the accessions to different subpopulations is indicated by color (G1: red, G2: blue; G3: green; G4: purple).
Table 3 Summary statistics for subpopulations detected by structure analysis based on nine chloroplast DNA markers.
An analysis of molecular variance (AMOVA) was performed to estimate the genetic variance among and within populations. The test showed that the genetic variation among individuals within subpopulations was 65.20%, higher than that between subpopulations, which was 34.80% (Table 4). The pairwise FST analysis indicated that the genetic distance between G2 and G3 was the smallest (0.25646), and the genetic distance between G2 and G4 was the largest (0.42528) (Table 5).
Table 4 The analysis of molecular variance (AMOVA) among and within four populations of 130 foxtail millet accessions identified by structure analysis.
Table 5 Pairwise FST among four populations of 130 foxtail millet accessions identified by structure analysis.
Construction of a Female Parent-Derivative Evolutionary Graph and Classification of Cytoplasm Types
We constructed a female parent-derivative evolutionary graph of the 130 accessions according to their maternal genealogy, and clustered them into four cytoplasmic types by combination with the results of population structure analysis (Figure 2). Type 1 (G1) contained 30 accessions; 8 were Dahuanggu derivatives, 3 were Mihuanggu derivatives, and 19 were from unknown maternal sources. As Dahuanggu and Mihuanggu both originated from Henan Province in China and are classified into the G1 group, we speculated that Dahuanggu and Mihuanggu were the same accession. Thus, the G1 type likely has the same cytoplasm as Dahuanggu, and can be designated as a Dahuanggu type. A total of 35 accessions were grouped in type 2 (G2), including 11 Riben60ri derivatives, two Zhengai2 derivatives, one Yingsuigu derivative, one Heizhigu derivative, and 20 accessions with unknown pedigrees. Since Zhengai2 is a derivative of Riben60ri (Wei et al., 1998), and the majority of the accessions in the G2 group (86.7%) were from Riben60ri, their cytoplasm type is considered as Riben60ri type. Type 3 (G3) consisted of 31 accessions; 12 were Moligu derivatives, and the remaining 19 accessions were from unknown maternal sources. The G3 group fully matches the Moligu cytoplasm type and was designated as a Moligu type. The remaining 34 accessions were categorized into type 4 (G4). Type 4 was a mixed subpopulation, including 4 Huangruangu derivatives, 2 Qinyuanmujizui derivatives, 1 Xiannong3 derivative, 1 Changsuihuang derivative, and 26 accessions having unknown pedigrees. Since half the accessions in the G4 group (50%) were from Huangruangu, their cytoplasm type is defined as Van Huangruangu type.
Figure 2 Classification of the cytoplasm types of the 130 foxtail millet accessions based on the female parent-derivative evolutionary graph and population structure analysis. The colors represent population membership identified using STRUCTURE in Figure 1.
Geographic and Ecological Distribution of Different Cytoplasmic Types
When we linked sample position on the two-dimensional PCA plot to the geographic distribution on the world map, we found that the accessions of same cytoplasmic type were closely related genetically, however, their geographic and ecological distributions were not the same (Figure 3). Among the 30 accessions of the Dahuanggu type (G1), 27 (90.01%) accessions were from central and south Hebei Province, and one each (3.33%) was from Henan Province, Liaoning Province, and the northwest Hebei Province. Except for one accession from the northwest Hebei Province (3.33%) belonging to the northwest spring millet ecotype, and one accession from Liaoning Province (3.33%) belonging to northeast spring millet ecotype, all other accessions (93.34%) were cultivated in the north China summer millet region. Among the Riben60ri type (G2), 29 (82.86%) accessions were from central and south Hebei Province, two (5.70%) were from Henan Province, and one each (2.86%) was from Jilin, Liaoning, Shanxi, and Shandong Provinces. These accessions mainly belong to three different ecological ecotypes: the north summer millet ecotype (91.43%), the northwest spring millet ecotype (2.86%), and the northeast China spring millet ecotype (5.71%).
Figure 3 Geographical and genetic distributions of 130 accessions of foxtail millet. The lower plot is from the principal component analysis (PCA). The lines link the position of accessions in the PCA graph and their geographic origin on the map. The colors indicate population membership identified using STRUCTURE in Figure 1.
Among the 31 accessions in the Moligu type (G3), 26 (83.87%) were from Beijing and central and south Hebei Province, and belonged to the north China summer millet ecotype. Four (12.90%) accessions were from Shanxi and Shaanxi Province and were grown in the northwest China spring millet region. One (3.23%) accession from Liaoning Province was distributed in the northeast China spring millet region. Among the Van Huangruangu type (G4), 25 (73.53%) accessions from central and south Hebei Province were grown in the north China summer millet region, eight accessions (23.53%) from Shanxi Province were distributed in the northwest China spring millet region, and one (2.94%) accession belonging to the northeast China spring millet ecotype was from Jilin Province. It is clear that the majority of the elite foxtail millet parental lines of each cytoplasmic type belong to the north China summer millet ecotype, although they have different geographic distributions.
We also focused on the distribution of basal maternal sources of the 130 tested accessions. The basal maternal parent of accessions in the Dahuanggu type (G1), from Henan Province, belonged to the north China summer millet region. Riben60ri, which was the basal maternal source of accessions in the Riben60ri type (G2), was originally introduced from Japan and is widely used in the north China summer millet region. The ecological distribution of Dahuanggu and Riben60ri coincided with the distribution of most of the accessions in the Dahuanggu (G1) and Riben60ri (G2) types. By contrast, accessions in the Moligu (G3) and Van Huangruangu (G4) types could be traced to the maternal parents Moligu and Huangruangu, respectively, both of which were from Shanxi Province and belonged to the northwest China spring millet ecotype. This was inconsistent with the ecological distribution of most of their derivatives.
Discussion
Plant organelle genomes (e.g., cpDNA and mtDNA) have numerous important applications in phylogenetic and population studies. However, the sequence information of organelle genomes for many plant species remains unclear, and cpDNA or mtDNA of high purity and yield are usually difficult to extract due to nuclear DNA interference. In a previous study, we developed a new protocol and successfully isolated high-quality cpDNAs for whole-genome sequencing from four different foxtail millet accessions (Gu56A, Gu572A, Datong28, and Datong29lv) (unpublished data). According to the sequence alignment, 45 SNPs and 9 InDels were identified, which provided the most direct sequence information for the development of molecular markers based on cpDNA. This information also ensured that the PCR products that were amplified with these markers were from the foxtail millet chloroplast genomes.
To our knowledge, this is the first study that utilized foxtail millet-specific chloroplast DNA primer pairs. Moreover, through screening, nine primers with polymorphisms were selected from 27 primers and were used for cytoplasmic variation analysis. The match rate between the population structure analysis and the female parent-derivative evolutionary graph was 86.67% (Figures 1 and 2). These results revealed that the nine cpDNA markers that we developed were sufficient for the classification of the cytoplasmic types. Our markers may also help create new cytoplasm sources for male sterile lines and enrich the cytoplasmic types of foxtail millet hybrids.
We also analyzed the foxtail millet cytoplasmic genetic diversity with 10 universal cpDNA and 17 mtDNA primer pairs with polymorphisms. All 130 foxtail millet accessions were classified into 3 cytoplasmic types, and the coincidence rates between the population structure and the female parent-derivative evolutionary graph were 35.56 and 42.22%, respectively (Supplemental Figures S1 and S2 and Supplemental Tables S4 and S5). These rates were much lower than the rate when specific chloroplast DNA primers were used. The cpDNA universal primers were developed based on complete chloroplast genomes of rice, Arabidopsis thaliana, and olive (Olea europaea L.) (https://www.ncbi.nlm.nih.gov/probe/?term=chloroplast+ssr). The set of mtDNA consensus primers were designed based on the fully sequenced mitochondrial genomes of Arabidopsis thaliana and Beta vulgaris (Duminil et al., 2002). Thus, it is possible that some DNA fragments that were amplified by these universal primers may not be from the cytoplasmic genome of foxtail millet. Furthermore, the genetic diversity was less informative than that revealed using specific cpDNA markers based on chloroplast genome sequences (Chaïr et al., 2005; Scarcelli et al., 2011). These results further suggest that the specific cpDNA primers we reported were more suitable for genetic variation, population analysis, and the classification of cytoplasmic types of foxtail millet.
Nuclear DNA microsatellite markers have been previously used in genetic diversity and population structure analysis for many crops, including 79 and 77 markers for foxtail millet (Liu et al., 2011; Wang et al., 2012a), 10 markers for wild rice (Oryza rufipogon Griff.) (Zhou et al., 2003), 96 markers for maize (Vigouroux et al., 2008), and 95 markers for sorghum (Wang et al., 2009). Along with the sequencing of the chloroplast genome, molecular markers based on cpDNA have been applied to population genetics and evolutionary studies in plants (Provan et al., 2001). In this study, nine chloroplast DNA markers specific to foxtail millet were used for population structure analysis. This number is similar to other studies in which chloroplast microsatellite markers were developed, such as Canary Island pine (Pinus canariensis C. Sm.) (six markers) (Gomez et al., 2003), wild soybean (Glycine soja Sieb. Et Zucc) (five markers) (He et al., 2012), and oilseed rape (Brassica napus L.) (14 markers) (Zamani-Nour et al., 2013). The size of circular DNA in the chloroplast ranges from 120 to 160 kilobases (kb) in most plants; this is much less than that the nuclear genome (Sugiura, 1992; Sugiura, 1995; Jansen et al., 2005). Thus, the number of cpDNA markers used in our research is sufficient for examining cytoplasmic genetic variation and population structure of foxtail millet accessions.
In this study, population structure analysis revealed the genetic relationships among different foxtail millet parental lines. Based on our results, we divided the 130 tested accessions into four groups. Most accessions in each group had highly similar genetic background (Figure 1). This resulted from the maternal inheritance mode of cpDNA (Shaver et al., 2006; Ravi et al., 2008) and further demonstrated that our population structure determination was reliable.
We detected a high level of polymorphism in the 9 chloroplast DNA markers across 130 foxtail millet accessions (Table 2). The mean gene diversity (0.37) and mean PIC value (0.30) were greater than in a previous report (0.29 and 0.23, respectively), which were identified using 23 consensus mtDNA primers across 111 accessions (Liu et al., 2014). This might be due to the fact that we used a more diverse selection of samples containing some accessions developed and released in the last few years (Supplemental Table S1), as indicated in our PCA analysis results (Figure 3). Before 2010, breeders primarily used locally adapted material as the female parent in order to develop new foxtail millet cultivars. This led to a narrow source of cytoplasm, which is illustrated in a previous study by Liu et al. (2014).
We found that the cytoplasmic types of Chinese foxtail millet are gradually diversifying. For example, among the 31 accessions in the Moligu type (G3), 83.87% accessions were grown in the north China summer millet region, while their basal maternal source, Moligu, was from Shanxi Province, which belongs to the northwest China spring millet region (Figure 3). This indicated that a greater number of foxtail millet male sterile lines were being used as the female parent cross-ecological regions in recent years. However, on the whole, the cytoplasm of Chinese foxtail millet is still of limited genetic variation, and the utilization of the foxtail millet germplasm is still primarily limited to the ecological distribution they belonged to originally. For example, the majority of accessions in Dahuanggu (G1) and Riben60ri (G2) types (93.34 and 91.43%, respectively) are cultivated in the north China summer millet region, and their basal maternal sources, Dahuanggu and Riben60ri, are also from (or widely applied) in the same ecological region (Figure 3). Thus, to avoid risks and losses from using of single source of cytoplasm in foxtail millet male sterile lines, finding and developing diverse cytoplasm types by using female parents introduced from different ecological regions and foreign countries is still key for improving hybrid cultivars.
It is worth noting that the Van Huangruangu type (G4) was a mixed subpopulation, which contained 25% Qinyuanmujizui derivatives, 12.5% Xiannong3 derivatives, 12.5% Changsuihuang derivatives, and 50% Huangruangu derivatives (Figure 2). This suggested that the accessions in the Van Huangruangu type (G4) may be derived from more than one single cytoplasmic type. Thus, the Van Huangruangu type (G4) may be a source for elite parental lines that contains different cytoplasmic types, and may therefore be used in foxtail millet breeding.
Data Availability Statement
The datasets generated for this study can be found in the GenBank under the following accession numbers: 56A (MK348603), 572A (MK348609), Datong28 (MK348605), and Datong29lv (MK348604).
Author Contributions
ZL and SL contributed to the study conception and design. TY, XS, CW, TL, and QL collected the foxtail millet materials. DLiu designed the cpDNA primers. SL, DLiu and QL genotyped the accessions, and determined the summary statistics. DLiu, DLia, and JW performed the population structure analysis, conducted the analysis of molecular variance, and calculated the pairwise FST values. JH and YC constructed the female parent-derivative evolutionary graph. The principal component analysis was performed by GF and DLiu. KD, RR, and LZ coordinated the experiments and data analysis. The first draft of the manuscript was written by YC. ZL and YC revised the article. All authors read and approved the final manuscript.
Funding
This work was supported by National Natural Science Foundation of China [grant number 31560094].
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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2019.01198/full#supplementary-material
References
Avise, J. C. (2000). Phylogeography: the history and formation of species. Cambridge, MA: Harvard University Press.
Baranwal, V. K., Mikkilineni, V., Zehr, U. B., Tyagi, A. K., Kapoor, S. (2012). Heterosis: emerging ideas about hybrid vigour. J. Exp. Bot. 63, 6309–6314. doi: 10.1093/jxb/errs321931
Chaïr, H., Perrier, X., Agbangla, C., Marchand, J. L., Dainou, O., Noyer, J. L. (2005). Use of cpSSRs for the characterisation of yam phylogeny in Benin. Genome 48, 674–684. doi: 10.1139/g05-018
Cheng, B., Zheng, Y., Sun, Q. (2015). Genetic diversity and population structure of Taxus cuspidata in the Changbai Mountains assessed by chloroplast DNA sequences and microsatellite markers. Biochem. Syst. Ecol. 63, 157–164. doi: 10.1016/j.bse.2015.10.009
Diao, X., Jia, G. (2017). “Foxtail millet germplasm and inheritance of morphological characteristics” in Genetics and genomics of Setaria. Eds. Doust, A., Diao, X. (Cham: Springer International Publishing), 73–92. doi: 10.1007/978-3-319-45105-3_5
Duminil, J., Pemonge, M. H., Petit, R. J. (2002). A set of 35 consensus primer pairs amplifying genes and introns of plant mitochondrial DNA. Mol. Ecol. Notes 2, 428–430. doi: 10.1046/j.1471-8286.2002.00263.x
Earl, D. A., vonHoldt, B. M. (2012). Structure harvester: a website and program for visualizing structure output and implementing the Evanno method. Conserv. Genet. Resour. 4, 359–361. doi: 10.1007/s12686-011-9548-7
Evanno, G., Regnaut, S., Goudet, J. (2005). Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol. Ecol. 14, 2611–2620. doi: 10.1111/j.1365-294X.2005.02553.x
Excoffier, L., Laval, G., Schneider, S. (2005). Arlequin (version 3.0): An integrated software package for population genetics data analysis. Evol. Bioinform. 1, 47–50. doi: 10.1177/117693430500100003
Falush, D., Stephens, M., Pritchard, J. K. (2003). Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics 164 (4), 1567–1587.
Gomez, A., Gonzalez-Martinez, S. C., Collada, C., Climent, J., Gil, L. (2003). Complex population genetic structure in the endemic Canary Island pine revealed using chloroplast microsatellite markers. Theor. Appl. Genet. 107, 1123–1131. doi: 10.1007/s00122-003-1320-2
Guo, M., Rupe, M. A., Yang, X., Crasta, O., Zinselmeier, C., Smith, O. S. (2006). Genome-wide transcript analysis of maize hybrids: allelic additive gene expression and yield heterosis. Theor. Appl. Genet. 113, 831–845. doi: 10.1007/S00122-006-0335-X
He, S., Wang, Y., Volis, S., Li, D., Yi, T. (2012). Genetic diversity and population structure: implications for conservation of wild soybean (Glycine soja Sieb. et Zucc) based on nuclear and chloroplast microsatellite variation. Int. J. Mol. Sci. 13, 12608–12628. doi: 10.3390/ijms131012608
He, S., Yang, Y., Morrell, P. L., Yi, T. (2015). Nucleotide sequence diversity and linkage disequilibrium of four nuclear loci in foxtail millet (Setaria italica). PloS One 10, e0137088. doi: 10.1371/journal.pone.0137088
Hooker, A. L., Smith, D. R., Lim, S. M., Beckett, J. B. (1970). Reaction of corn seedlings with male-sterile cytoplasm to Helminthosporium maydis. Plant Dis. Rep. 54, 708–712.
Huang, X., Huang, M., Liu, H., Zhao, C., Wang, J. (2018). Effects of annual precipitation and population density on tiller-earing and yield of Zhangzagu 5 under film mulching and hole sowing. Crops 4, 106–113. doi: 10.16035/j.issn.1001-7283.2018.04.018
Jakobsson, M., Rosenberg, N. A. (2007). CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics 23, 1801–1806. doi: 10.1093/bioinformatics/btm233
Jansen, R. K., Raubeson, L. A., Boore, J. L., Pamphilis, C. W., Chumley, T. W., Haberle, R. C. (2005). Methods for obtaining and analyzing whole chloroplast genome sequences. Methods Enzymol. 395, 348. doi: 10.1016/S0076-6879(05)95020-9
Jia, G., Liu, X., Schnable, J. C., Niu, Z., Wang, C., Li, Y. (2015). Microsatellite variations of elite Setaria varieties released during last six decades in China. PloS One 10, e0125688. doi: 10.1371/journal.pone.0125688
Li, S., Liu, D., Li, Q., Chen, C., Liu, Z. (2018). Breeding of advantage heterosis foxtail millet hybrid Jizagu5 with characteristic of herbicide-resistant and the research of simplified cultivation technique in seedlings. J. Tangshan Normal Univ. 40, 54–58. doi: 10.3969/j.issn.1009-9115.2018.03.014
Liu, K., Muse, S. V. (2005). PowerMarker: an integrated analysis environment for genetic marker analysis. Bioinformatics 21, 2128–2129. doi: 10.1093/bioinformatics/bti282
Liu, K., Wang, L., Wei, J., Zhu, X., Wu, Q. (1992). Reaction of rice male sterile cytoplasm of wild abortion type to the infection of Pyricularia oryzae. Sci. Agric. Sin. 25, 92.
Liu, Z., Cheng, R., Li, X. (1996). The pedigree analysis and evaluation of north China summer millets. Crops 5, 24.
Liu, Z., Cheng, R., Zhang, F., Xia, X., Shi, Z., Hou, S. (2006). Millet variety in boreali-sinica summer millets region and its pedigree evolution and analysis on genetic foundation. Acta Agric. Boreali-Sin. Suppl. 21, 103–109. doi: 10.3321/j.issn:1000-7091.2006.z2.025
Liu, Z., Bai, G., Zhang, D., Zhu, C., Xia, X., Cheng, R. (2011). Genetic diversity and population structure of elite foxtail millet [Setaria italica (L.) P. Beauv.] germplasm in China. Crop Sci. 51, 1655. doi: 10.2135/cropsci2010.11.0643
Liu, Z., Zhang, T., Li, C., Bai, G. (2014). Genetic diversity and classification of cytoplasm of Chinese elite foxtail millet [Setaria italica (L.) P. Beauv.]. Germplasm Crop Sci. 54, 659. doi: 10.2135/cropsci2012.11.0646
McGrath, S., Hodkinson, T. R., Barth, S. (2007). Extremely high cytoplasmic diversity in natural and breeding populations of Lolium (Poaceae). Heredity 99, 531–544. doi: 10.1038/sj.hdy.6801030
Norouzi, M., Talebi, M., Sayed-Tabatabaei, B.-E. (2012). Chloroplast microsatellite diversity and population genetic structure of Iranian pomegranate (Punica granatum L.) genotypes. Sci. Hortic. 137, 114–120. doi: 10.1016/j.scienta.2012.01.034
Pritchard, J. K., Stephens, M., Donnelly, P. (2000). Inference of population structure using multilocus genotype data. Genetics 155 (2), 945–959.
Provan, J., Powell, W., Hollingsworth, P. M. (2001). Chloroplast microsatellites: new tools for studies in plant ecology and evolution. Trends Ecol. Evol. 16, 142–147. doi: 10.1016/S0169-5347(00)02097-8
Rajendrakumar, P., Biswal, A. K., Balachandran, S. M., Ramesha, M. S., Viraktamath, B. C., Sundaram, R. M. (2007). A mitochondrial repeat specific marker for distinguishing wild abortive type cytoplasmic male sterile rice lines from their cognate isogenic maintainer lines. Crop Sci. 47, 207–211. doi: 10.2135/cropsci2006.06.0365
Ravi, V., Khurana, J. P., Tyagi, A. K., Khurana, P. (2008). An update on chloroplast genomes. Plant Syst. Evol. 271, 101–122. doi: 10.1007/s00606-007-0608-0
Rosenberg, N. A. (2010). DISTRUCT: a program for the graphical display of population structure. Mol. Ecol. Notes 4, 137–138. doi: 10.1046/j.1471-8286.2003.00566.x
Scarcelli, N., Barnaud, A., Eiserhardt, W., Treier, U. A., Seveno, M., d’Anfray, A. (2011). A set of 100 chloroplast DNA primer pairs to study population genetics and phylogeny in monocotyledons. PloS One 6, e19954. doi: 10.1371/journal.pone.0019954
Scheifele, G. L., Whitehead, W., Rowe, C. (1970). Increased susceptibility to Southern leaf spot (Helminthosporium maydis) in inbred lines and hybrids of maize with Texas male-sterile cytoplasm. Plant Dis. Rep. 54, 501–503.
Shaver, J. M., Oldenburg, D. J., Bendich, A. J. (2006). Changes in chloroplast DNA during development in tobacco, Medicago truncatula, pea, and maize. Planta 224, 72–82. doi: 10.1007/s00425-005-0195-7
Shimamoto, Y., Abe, J., Gao, Z., Gai, J., Thseng, F.-S. (2000). Characterizing the cytoplasmic diversity and phyletic relationship of Chinese landraces of soybean, Glycine max, based on RFLPs of chloroplast and mitochondrial DNA. Genet. Resour. Crop Ev. 47, 611–617. doi: 10.1023/a:1026538907387
Shu, J., Liu, Y., Li, Z., Zhang, L., Fang, Z., Yang, L. (2016). Detection of the diversity of cytoplasmic male sterility sources in Broccoli (Brassica Oleracea var. Italica) using mitochondrial markers. Front. Plant Sci. 7, 927. doi: 10.3389/fpls.2016.00927
Skuza, L., Szućko, I., Filip, E., Strzała, T. (2019). Genetic diversity and relationship between cultivated, weedy and wild rye species as revealed by chloroplast and mitochondrial DNA non-coding regions analysis. PloS One 14, e0213023. doi: 10.1371/journal.pone.0213023
Smith, D. R., Hooker, A. L., Lim, S. M. (1970). Physiologic races of Helminthosporium maydis. Plant Dis. Rep. 54, 819–822.
Springer, N. M., Stupar, R. M. (2007). Allelic variation and heterosis in maize: How do two halves make more than a whole? Genome Res. 17, 264–275. doi: 10.1101/gr.5347007
Sugiura, M. (1992). The chloroplast genome. Plant Mol. Biol. 19 (1), 149–168. doi: 10.1007/BF00015612
Vigouroux, Y., Glaubitz, J. C., Matsuoka, Y., Goodman, M. M., Sanchez, G, J., Doebley, J. (2008). Population structure and genetic diversity of New World maize races assessed by DNA microsatellites. Am. J. Bot. 95, 1240–1253. doi: 10.3732/ajb.0800097
Wang, Y., Li, H., Wang, G., Tian, G. (1998). Breeding of foxtail millet highly sterile line “Chang10A”. Gansu Agr. Sci. Techno. 12, 12–13.
Wang, M. L., Zhu, C., Barkley, N. A., Chen, Z., Erpelding, J. E., Murray, S. C. (2009). Genetic diversity and population structure analysis of accessions in the US historic sweet sorghum collection. Theor. Appl. Genet. 120, 13–23. doi: 10.1007/s00122-009-1155-6
Wang, C., Jia, G., Zhi, H., Niu, Z., Chai, Y., Li, W. (2012a). Genetic diversity and population structure of Chinese foxtail millet [Setaria italica (L.) Beauv.] landraces. G3. Genes Genomes Genet. 2, 769–777. doi: 10.1534/g3.112.002907
Wang, Q., Zhang, Y., Fang, Z., Liu, Y., Yang, L., Zhuang, M. (2012). Chloroplast and mitochondrial SSR help to distinguish allo-cytoplasmic male sterile types in cabbage (Brassica oleracea L. var. capitata). Mol. Breed. 30, 709–716. doi: 10.1007/s11032-011-9656-9
Wei, L., Meng, Z., Qi, M., Ding, Y. (1998). Innovation of foxtail millet dwarfing germplasm resource Zhengai2. Crops 5, 28.
Yuan, L. (1997). Current status and developing prospects in two-line hybrid rice research in China. Res. Agric. Modern. 1, 1–3.
Zamani-Nour, S., Clemens, R., Möllers, C. (2013). Cytoplasmic diversity of Brassica napus L., Brassica oleracea L. and Brassica rapa L. as determined by chloroplast microsatellite markers. Genet. Resour. Crop Ev. 60, 953–965. doi: 10.1007/s10722-012-9891-x
Keywords: foxtail millet, elite parental lines, genetic diversity, cytoplasmic type, chloroplast deoxyribonucleic acid
Citation: Liu D, Cui Y, He J, Li S, Li Q, Liang D, Wang J, Shi X, Wang C, Dong K, Liu T, Zhang L, Ren R, Yang T, Feng G and Liu Z (2019) Genetic Diversity and Classification of the Cytoplasm of Chinese Elite Foxtail Millet [Setaria italica (L.) P. Beauv.] Parental Lines Revealed by Chloroplast Deoxyribonucleic Acid Variation. Front. Genet. 10:1198. doi: 10.3389/fgene.2019.01198
Received: 12 August 2019; Accepted: 29 October 2019;
Published: 22 November 2019.
Edited by:
Longjiang Fan, Zhejiang University, ChinaReviewed by:
Jiasheng Wu, Zhejiang Agriculture and Forestry University, ChinaGuanglin He, Sichuan University, China
Copyright © 2019 Liu, Cui, He, Li, Li, Liang, Wang, Shi, Wang, Dong, Liu, Zhang, Ren, Yang, Feng and Liu. 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: Suying Li lisuying65@126.com; Tianyu Yang 13519638111@163.com; Gang Feng fg7195@aliyun.com; Zhengli Liu, liuzhengli65@126.com
†These authors have contributed equally to this work