Original Research ARTICLE
Development of SSR Markers Based on Transcriptome Sequencing and Association Analysis with Drought Tolerance in Perennial Grass Miscanthus from China
- Department of Grassland Science, Animal Science and Technology College, Sichuan Agricultural University, Sichuan, China
Drought has become a critical environmental stress affecting on plant in temperate area. As one of the promising bio-energy crops to sustainable biomass production, the genus Miscanthus has been widely studied around the world. However, the most widely used hybrid cultivar among this genus, Miscanthus × giganteus is proved poor drought tolerance compared to some parental species. Here we mainly focused on Miscanthus sinensis, which is one of the progenitors of M. × giganteus providing a comparable yield and well abiotic stress tolerance in some places. The main objectives were to characterize the physiological and photosynthetic respond to drought stress and to develop simple sequence repeats (SSRs) markers associated with drought tolerance by transcriptome sequencing within an originally collection of 44 Miscanthus genotypes from southwest China. Significant phenotypic differences were observed among genotypes, and the average of leaf relative water content (RWC) were severely affected by drought stress decreasing from 88.27 to 43.21%, which could well contribute to separating the drought resistant and drought sensitive genotype of Miscanthus. Furthermore, a total of 16,566 gene-associated SSRs markers were identified based on Illumina RNA sequencing under drought conditions, and 93 of them were randomly selected to validate. In total, 70 (75.3%) SSRs were successfully amplified and the generated loci from 30 polymorphic SSRs were used to estimate the genetic differentiation and population structure. Finally, two optimum subgroups of the population were determined by structure analysis and based on association analysis, seven significant associations were identified including two markers with leaf RWC and five markers with photosynthetic traits. With the rich sequencing resources annotation, such associations would serve an efficient tool for Miscanthus drought response mechanism study and facilitate genetic improvement of drought resistant for this species.
As a result of global climate change, drought has become a critical environmental factor negatively affecting on plant growth, development, survival as well as crop yield (Boyer, 1982; Seki et al., 2007; Shinozaki and Yamaguchi-Shinozaki, 2007; Maughan et al., 2012; Malinowska et al., 2017). The genus Miscanthus is a perennial rhizomatous C4 grass which widely recognized as a promising bio-energy crops to sustainable biomass production (Hodkinson et al., 2002; Somerville et al., 2010). Despite the water use efficiency of C4 crops is superior to C3 crops (Gowik and Westhoff, 2011; Ings et al., 2013), water availability still limits the maximum yields achievable. Miscanthus × giganteus, crossing between Miscanthus sacchariflorus and Miscanthus sinensis, is the most widely used hybrid cultivar among this genus in Europe and North America as a biomass feedstock (Greef et al., 1997; Hodkinson et al., 2002; Heaton et al., 2008, 2009; Hastings et al., 2009). However, owing to its natural sterility and a narrow genetic base, it is difficult to improve M. × giganteus with desirable trait through breeding and just one genotype is available for agricultural production in most areas (Nishiwaki et al., 2011; Dwiyanti et al., 2014) Furthermore, poor WUE of M. × giganteus had been reported compared to some parental species (Clifton-Brown and Lewandowski, 2000). As a progenitor of M. × giganteus, new varieties M. sinensis widely adapted to a wide range of geographical conditions (Clifton-Brown and Lewandowski, 2002; Nie et al., 2016) and showed extensive genetic diversity, which is desirable to provide improved drought tolerance and comparable biomass yield. Compared with other C4 grasses like maize (Zea mays L.) and Saccharum (Saccharum officinarum L.), M. sinensis has a relatively high efficient photosynthesis and water use efficiency (Greef et al., 1997; Long and Spence, 2013) owing to its perennial growth habit and extensive root system (Byrt et al., 2011; van der Weijde et al., 2013).
In order to maximize the avoidance competition of arable land with major grain crops (Quinn et al., 2015; van der Weijde et al., 2016), most non-food grasses like Miscanthus were commonly planted in vast marginal land, which always limited water supply, although drought or water deficit always affects crop yield seriously worldwide (Cattivelli et al., 2008). Generally, wild plants showed high tolerant performance against abiotic stresses (Ellis et al., 2000; Bartels, 2001; Mittova et al., 2004), but it may not necessarily favor growth or biomass accumulation under stress rather survival. For M. sinensis, the attractive characteristics is efficient biomass production rather than just survival and yield is always strongly linked to water availability. Therefore, stabilizing crop performance under drought stress, increasing crop potencial production of per unit of water consumption, will become a critical goal for M. sinensis breeding (Ings et al., 2013).
Previous study had been reported that M. × giganteus showed poorly responds to low water availability (Ings et al., 2013). van der Weijde et al. (2016) evaluated the growth and quality of 50 diverse Miscanthus genotypes under drought stress and proved that genotypic variation for drought stress exist in this genus. Malinowska et al. (2017) used a high-throughput phenomics facility to screen 47 Miscanthus genotypes about biomass accumulation and water use under drought stress. All the related researches had offered valuable insights into the varied side of drought-induced responses and the growth as well as yield characteristics of Miscanthus under drought stress. However, despite an increasing concern about the drought resistant of Miscanthus, little is known about the physiological and photosynthetic mechanisms response to drought stress.
In addition, the type and timing of physiological or photosynthetic changes induced upon drought stress always vary in different species and genotypes (Centritto et al., 2009; Costa et al., 2012). The conventional method for screening genotypes with desirable traits is time-consuming and labor intensive (Furbank and Tester, 2011). Thereby, it is necessary to develop efficient methods for screening the favorable traits associated with drought resistance around the wide range Miscanthus population. Molecular markers are essential tools for estimating genetic variation, selecting parental or breeding line and marker-assisted breeding for improvement of crop desirable traits (Hung et al., 2009; Clark et al., 2014; Nie et al., 2014, 2016; Liu et al., 2015). The number of available efficient molecular markers for Miscanthus is still not enough, especially the function markers associated with drought tolerance. SSRs are valuable molecular markers for most plants that genome sequence is not available and have been widely used for molecular MAS of major crop plants. Compared with traditional marker development methods depending on the publicly available cDNA and genomic libraries, recent advances high-throughput next-generation RNA sequencing (NGS) platform provide a high-efficiency tool for developing large amounts of gene-based SSRs markers in most crop species (Hu et al., 2004; Gupta and Prasad, 2009; Varshney, 2009) and proved powerful associated with putative functions involving varies abiotic stress tolerance and important agronomic traits (Kane and Rieseberg, 2007; Yu et al., 2011; Nie et al., 2016).
In this study, the main objective was to characterize the physiological and photosynthetic respond to drought stress with a population of 44 wild collections of Miscanthus. Mapping the drought resistant mechanism of Miscanthus could help us to well predict and assess desirable traits of the crops in extreme environments. Furthermore, the development of SSRs markers associated with drought tolerance by transcriptome sequencing of Miscanthus will facilitate the speed of the screening across the genotypes, which aimed to provide a platform for M. sinensis genetic improvement and molecular marker-assisted breeding.
Materials and Methods
Plant Material and Drought Treatment
A total of 44 diverse Miscanthus genotypes comprising 41 M. sinensis (NO. 1 to NO. 41), 2 Miscanthus floridulus (NO. 48 to NO. 49) and 1 M. × giganteus (NO. 50) species used in this study were mainly collected from the Southwest of China (Supplementary Table S1). Each genotype was propagated to six individuals (three for control and three for treatment) by dividing rhizome and grown in plastic pots (20 cm × 25 cm) with mixture soil (50% loam with 50% fine sandy). All materials were grown in greenhouse with a completely randomized design from April to July in 2014. Before subjected to water stress, all the plants grew with sufficient conditions. The experiment was designed to evaluate the genotypic difference in physiological and photosynthetic aspect responding to water stress. Prior to treat, all the replicates were well watered with the average of 90% SWC measured by Soil Moisture Equipment TDR300 (Santa Barbara, CA, USA). Naturally water stress was applied by stopping water for 30 days. The photosynthetic traits were measured at the end of the drought stress in control and treatment first, and then leaf samples were collected for determination physiological indexes from three replications.
Physiological and Photosynthetic Measurements
The Pn, Gs, Ci, and Tn were determined by photosynthetic system (LI-6400, Li-Cor, Lincoln, NE, USA) with the CO2 concentration was 400 μmol mol-1 and photosynthetic active radiation (PAR) was 1,200 μmol⋅m-2s-1. Gs and Ci were measured under the condition that with a saturated light intensity of 1,000 μmol m-2s-1. Ten randomly selected leaves were measured to determine related traits for each replicate (Zhao et al., 2014).
The RWC of leaves was measured by drying method (Barrs and Weatherley, 1962). Approximate 0.2 g (FW) fresh leaves randomly selected from each pot were wrapped in absorbent paper and put into the 50mL tube with filled water. Tubes were placed in the dark place for 24 h, and then take out the leaves, dry the surface moisture, and weight SW. Followed the leaf samples were placed in the oven with 75°C for 3 days to constant weight, weight DW. The RWC was calculated as follows: RWC (%) = (FW - DW/SW - DW) × 100. For the extraction of MDA, a total of 0.1 g leaf samples were quickly frozen in liquid nitrogen (2 min), sufficiently grind in 2 volume solutions containing 50 mM phosphate buffer (pH = 7.8) and 1% (w/v) polyvinyl polypyrrolidone. The supernatant was then collected after the homogenate centrifuged at 12,000 rpm for 30 min and the content of MDA was examined (Dhindsa and Matowe, 1981).
Analysis of variance of the physiological and photosynthetic index was conducted by using SPSS 17.0 software (IBM, Armonk, NY, USA) and the Pearson correlation analysis of all traits was performed. Effects of both treatment and genotypes were determined by using the model of Least Significant Difference test (LSD).
Identification and Validation of SSRs
A total of six RNA samples from M. sinensis under drought stress were extracted based on the manufacturer’s instructions of the RNeasy Plant Mini Kit (Qiagen, Valencia, CA, USA). RNA purity, concentration, and integrity were checked by using the RNA Nano 6000 Kit for the Agilent 2100 Bioanalyzer 2100 System (Agilent Technologies, USA). The qualified samples were primary stored at –80°C after isolation and quality assessment. Each cDNA library was constructed using 5 μg of total RNA and six cDNA sequencing libraries of M. sinensis were constructed using the NEB Next® Ultra TM RNA Library Prep Kit for Illumina (New England Bio-labs, USA). Initially, the total RNA was treated with RNase-free DNase I (NEB) for 30 min at 37°C and Poly (A) mRNA was isolated using poly-T oligo-linked magnetic beads. After purification, the poly (A)-containing mRNA was fragmented into 200–250 bp pieces using fragment buffer (Ambion), and the first-strand cDNA was synthesized using random hexamer primers and the short fragments as templates. Then the RNase H and DNA polymerase I were applied to the products for second-strand cDNA synthesizing (16°C for 2 h). Finally, before hybridization, NEB Next Adaptor with hairpin loop structure was ligated to the cDNA and the 3′ ends of the DNA fragments were adenylated. Subsequently, the purity and quality of cDNA fragments library were evaluated using the Agilent Bioanalyzer 2100 system. After generating the clusters on a cBot System using the TruSeq PE Cluster kit v3-cBot-HS (Illumina), Illumina sequencing (paired-end technology in Illumina HiSeq 2000 platform) of the six libraries was performed for RNA-seq analysis. The de novo assembly of RNA-seq was conducted using Trinity1 with employing the sequence splicing and redundancy removing. In total, 114,747 assembled unigenes were directly used to identify SSRs based on the software MISA2 (MIcroSAtellite). SSR search criteria were conducted based on perfect mono-, di-, tri-, tetra-, penta-, and hexa-nucleotide motifs minimum number of six, five, five, four, and four repeats, respectively. Primer3 v2.23 was used to design primers in the flanking regions of the SSRs.
A total of 93 SSR markers (Table 1) were randomly selected to validation in 54 plant materials including 44 M. sinensis, 5 M. floridulus, 1 M. sacchariflorus, 1 M. × giganteus, and 3 Hemarthria cultivars (Supplementary Table S1). The genomic DNA was extracted from the young leaves using plant DNA Genomic Extraction Kit (Tiangen Biotech, China). PCR amplifications were performed in 20 μl total volume reactions, including 3 μl genomic DNA (15 ng/μl), 0.4 μl Golden DNA Polymerase, 10 μl 2X Reaction Mix, 0.5 μl of each primer, and 5.6 μl ddH2O. The PCR program was as follows: an initial denaturation at 94°C for 1 min, followed by 5 cycles of 1 min annealing at 60°C, and 1 min extension at 72°C, then followed by 35 cycles of 1 min denaturation at 94°C, 1 min annealing at 55–58°C, and 1 min extension at 72°C. A final 10 min extension at 72°C was used before the PCR reactions were completed. The obtained PCR products were separated on 6% PAGE under 350 volts for 2 h and gel was stained by AgNO3 solutions.
Population Structure, Clustering, and Association Analysis
Population structure analysis was performed using STRUCTRE v2.3.4 software (Pritchard et al., 2000) with the “admixture model”. Briefly, the software default model parameters were set with 10,000 iterations and 10,000 replications of MCMC. 20 independent runs were applied in this study with varying the assumed number of genetic groups (K) from one to six for each run. Maximum likelihood and delta K (ΔK) values were used to determine the optimum number of subgroups (Pritchard et al., 2000; Evanno et al., 2005). In clustering analysis, the UPGMA dendogram were construct by using similarity coefficients based on SAHN module in the NTSYS-pc version 2.10 software.
The association analysis was conducted by using TASSEL 2.1 software to reveal the association between drought related traits and marker alleles (Bradbury et al., 2007). Four models (simple linear, Q, K, and Q+K) were tested by using Quantile–quantile (QQ) plots to identify the best model fitting drought tolerance related traits for association mapping in M. sinensis populations. The significant threshold for selecting the associations between alleles and traits was set at P < 0.01. Phenotypic variation explained (R2) indicated the fixed marker effects.
Results and Discussion
Physiological and Photosynthetic Respond
Relative water content determines plant water status when exposed to drought stress. As a physiological indicators that directly response to the extent of plant water deficiency in drought environments, all plants demonstrated high level of leaf RWC in control conditions with an average RWC of 88.27% while evident effect in drought stress with a decrease to 43.21% (Table 2). The leaf RWC significantly difference among genotypes based on ANOVA analysis under drought stress which ranged from 12.44 to 66.11%, indicating that leaf RWC could well separate the drought resistant and drought sensitive genotype of Miscanthus. The MDA content directly reflected the stability of cellular membrane in the leaves. After a period of drought stress treatment, the MDA content largely increased than control from 11.71 to 54.63 nmol/g FW and showed an extensive variation under drought conditions (19.62–143.78 nmol/g). For the photosynthetic productivity, Pn, Gs, Ci, and Tn were all decreased under drought stress, and significant differences among genotypes were observed under control and drought stress (Table 2).
TABLE 2. The minimum, maximum, mean, standard deviation (SD), and F-value of RWC, MDA, Pn, Gs, Ci, and Tn under control (C) and drought stress (D) of Miscanthus population.
As an undomesticated new bioenergy crop, several characteristics related to drought tolerance successfully distinguish Miscanthus from many other crops. Since high biomass yield is the primary goal by meeting Miscanthus when all above ground biomass is harvested annually, the drought tolerance assessment criterion seems clearly establishment. However, for optimum and stable productivity, at least a 3 years establishment is required after transplanting of Miscanthus, which is a typical perennial grass species (Clifton-Brown and Lewandowski, 2000; Anzoua et al., 2015; Nie et al., 2016). Directly evaluating the biomass productive of Miscanthus in the field is time and labor consuming under drought conditions. Greenhouse experiment is favorable for detecting the mechanism responding to drought stress which can be identified some reliably measured trait to predict future performance for M. sinensis resistant breeding programs.
Leaf RWC reflecting the plant water status, MDA content representing the extent of cellular membrane peroxidation and Pn, Gs, Ci, and Tn comprehensively showing the plant biosynthesis and catabolism process, which provide rapid and easy measurements for screening drought-tolerant plant materials and characterize drought tolerance response mechanism (O’Neill et al., 2006; Jansen et al., 2009; Jiang et al., 2009). Significant correlations were observed between Pn, Gs, Ci, and Tn under well-watered and drought conditions, respectively (Table 3). As expect, Pn significant positively correlated with Gs and Tn, and negatively correlated with Ci under drought stress. Leaf RWC significantly correlated with Pn, Gs, Ci, and Tn under well-watered conditions, but such correlations were not observed under drought stress. MDA had no significant correlation with any other measured traits under control, but significantly correlated with RWC under drought stress. The correlation analysis results showed that there is no significant correlation between the physiological and photosynthetic responds to drought stress, indicating that varied and different regulating mechanism was existed in the panel of Miscanthus accessions.
TABLE 3. Correlation coefficients among RWC, MDA, Pn, Gs, Ci, and Tn under control and drought stress in Miscanthus population‡.
In order to avoid damage caused by drought stress, one hand the plant reduced water transpiration and improved water use efficiency to protect the photosynthetic mechanism through stomatal closure, on the other hand, adjusted osmotic pressure to increase water uptake and maintained the integrity of the cell structure and function (Yu et al., 2013). In this study, a significant correlation between RWC and MDA content were detected under drought condition and the decreasing of RWC and increasing of MDA content significantly different from that of the control, potentially indicating the occurrence of drought stress, suggesting that they are sufficient parameters of drought tolerance in greenhouse experiment. Generally, drought induced ABA synthesis, which causing stomatal closure to prevents water loss through transpiration and play an important role in maintaining leaf water status (Yu et al., 2013). Previously studies had been showed that M. × giganteus exhibited little stomatal regulation ability under mild drought stress compared to M. sinensis (Clifton-Brown and Lewandowski, 2000). However, the significant correlation between RWC and Pn, Gs, Ci, or Tn were just detected in well watered condition in this study, indicating that Miscanthus holding moisture could not through the stomatal closure pathway but antioxidant pathways to maintain the cell structure and function under drought stress.
In addition, different phenotypic responses of Miscanthus to drought stress are important for marker-trait association. The extensive variation observed among the evaluated genotypes regarding physiological and photosynthetic respond indicating that the large range of comprehensive drought tolerance among the evaluated genotypes. The genotypes tested may become important candidates for investigating of mechanisms of plant respond to drought stress and could possibly be used in breeding programs. However, the potential variability in short-term versus long-term drought response should be taken into account in order to further develop new varieties of drought tolerant.
Identification and Validation of SSRs
Simple sequence repeats are considered as one of the most useful molecular markers due to their broad range of applications in genotype identification, genetic map construction, and marker-trait association. Using the Illumina HiSeq 2000 platform, a total of 349,393,396 raw reads were generated (accession numbers: SRP095822 in the NCBI SRA database) and a total of 316,200,846 high-quality clean reads were assembled into 114,747 unigenes after filtering and trimming the raw reads. Using MISA script, the microsatellites were searched across all unigenes and a total of 16,566 SSRs were identified with motif lengths ranging from one to six bp (Table 4). The tri-nucleotide repeats, di-nucleotide repeats and mono-nucleotide repeats were the three most abundant SSRs types in this study and were responsible for 54.63%, 23.84% and 12.81% of the total amount of SSRs, respectively (Table 4). However, the mono-nucleotides were excluded in the subsequent analysis owing to the presence of homopolymer affecting the sequencing quality (Gilles et al., 2011). The number of the given repeat unit of SSRs ranged from 4 to >10 and with the number of repeat units increased, the frequency of the given SSR structure progressively decreased (Table 5). As for the two most abundant repeat motif types (di- and tri-nucleotides), the frequency of AG/TC motif type account for 52.00% in di-nucleotide repeat motifs, and that the frequency of CCG/GGC was the most abundant motif type in the tri-nucleotide, accounting for 41.80% (Table 6).
TABLE 5. Summary information on frequencies of different SSR repeat motif types related to variation of repeat unit numbers in M. sinensis.
The results were very consistent with previous studies that the tri-nucleotide was the most abundance repeat among the types of SSRs (Kaur et al., 2012; Yates et al., 2014; Huang et al., 2016; Tian et al., 2017). In comparison to other nucleotide repeats, the tri-nucleotide expansions always lead to a homopolymeric amino acid run in an unaltered protein due to the retention of the original reading frame in translated regions (Kaur et al., 2012). However, the amplification of repeats other than tri-nucleotide is prevented in coding DNA, resulting in frame-shift mutations (Metzgar et al., 2000), which could be a possible explanation for the abundance of tri-nucleotide repeats. In the other hand, the AG/TC of di-nucleotide SSRs and CCG/GGC of tri-nucleotide SSRs were reported to be the largest number of motif in Hemarthria and grain crops (Temnykh et al., 2000; Kantety et al., 2002; Huang et al., 2016). The inverse GA/CT motif can describe multiple codons when transcribed into mRNA. Furthermore, a high proportion (8 and 10%) of amino acids Ala and Leu were translated from inverse GA/CT motif contributing it widely distribution across the genome (Kantety et al., 2002).
To validate the SSRs markers identified by RNA-seq, a total of 93 SSR primer pairs were randomly selected for validation including different motifs. Validation of these SSRs was based on the different size among 51 different Miscanthus genotypes and 3 Hemarthria species (Supplementary Table S1). In total, 70 (75.3%) SSRs were successfully amplified whereas 30 have polymorphic, including 9 for di-nucleotide, 16 for tri-nucleotide and 5 SSRs with motif lengths greater than 3 (Table 7). A total of 110 polymorphic loci generating from 30 SSRs were determined to be used for cluster analysis. The UPGMA dendrogram constructed by genetic similarity (GS) data obviously showed that three clusters were formed when the GS coefficient was 0.52. The group 1 accessions were mainly Miscanthus sinensis collected from the south of Sichuan and Chongqing which latitude was lower than N 31°. The M. sinensis genotypes of group 2 were primarily collected from the North of Sichuan and other genotypes of M. floridulus, M. sacchariflorus, M. × giganteus were also assigned into this groups. All three Hemarthria cultivars were classed to Group 3 as a foreign species (Figure 1). These results from the cluster analysis suggested that the Miscanthus SSR markers identified in this study have a high level of polymorphism clearly differentiating the Miscanthus species. More importantly, they had been proved having high transferability in the Andropogoneae tribe.
FIGURE 1. Dendrogram of 51 Miscanthus and 3 Hemarthria genotypes based on genetic similarity (GS) coefficient.
Previous studies on Miscanthus marker developing are limited to the number of available for use, and always markers from Tribe Andropogoneae (Poaceae) C4 crops such as sorghum (Sorghum bicolor L. Moench), maize (Zea mays L.), and sugarcane (Saccharum officinarum L.) were used for genetic studies (Ho et al., 2011; Kim et al., 2012; Chae et al., 2014; Yook et al., 2014; Nie et al., 2016). Although these large number of SSRs have been successfully applied to M. sinensis studies and proven with high transferability (Hernandez et al., 2001; Zhao et al., 2011; Xu et al., 2013; Chae et al., 2014; Yook et al., 2014; Nie et al., 2016), the specificity and function of markers seems less than satisfactory. The SSRs developed from Miscanthus transcriptome under drought conditions were aimed to the function regions among expressed unigenes, which could be facilitate better discovering of the marker-trait associations and higher efficiency.
Population Structure and Association Analysis
Population structure of 44 genotypes was analyzed by using STRUCTURE V2.3.3 software. A total number of 110 filtered marker loci (by dropping the loci with minor allele frequency less than 5%) retained for structure and association analysis. Two optimum subgroups were determined by maximum likelihood and delta K (ΔK) values (Figure 2). Accordingly, 24 genotypes were assigned to group 1, 14 genotypes were to group 2, and 6 genotypes were referred not clearly assigned to the specific groups based on a significant membership threshold (Q-value) of 0.60 (Figure 3).
FIGURE 2. Optimal value of K determined. (A), Mean log-likelihoods and their standard deviations from runs assuming different numbers of subpopulations (K). (B) Values of the delta K (ΔK), which tends to peak at the value of K that corresponds to the highest hierarchical level of substructure.
FIGURE 3. Population structure analysis of 44 Miscanthus genotypes from southwest China. Number 1-44 on the x-axis indicate the individual and number on the y-axis show the group membership.
The association analysis was performed by taking Q (population structure) and K (relative kinship) into consideration using GLM and MLM model. Drought resistant related traits measured in this study under drought and control condition were used to test the model of Q, K, Q+K and simple linear model in QQ plots (Figure 4). Generally, in most conditions, the behavior of Q+K model was the best approximation to the excepted cumulative distribution of P-values. By implementing the Q+K models, a significant improvement was observed when compared with the other models. Based on Q+K model, a total of seven loci were detected significantly associate with RWC, Pn, Tn, Gs, and Ci (P < 0.01) under drought stress. Among them, 2 SSRs developing from CL12059.Contig2_All_5665_1 and CL10784.Contig3 _All_5453_1 were significantly associated with RWC, 1 SSRs from CL10592.Contig1_All_5404_1 associate with Pn and Gs simultaneously, 1 SSRs from CL11502.Contig2_All_5573_1 associate with Tn, 1 SSRs from CL10112.Contig3_All_5301_4 associate with Gs and 1 SSRs from CL10164.Contig1_All_5323_2 associate with Ci (Table 8). The contribution of single significant associated markers to the phenotypic variation was varied ranging from 11.89 to 25.30% (Table 8).
FIGURE 4. Quantile–quantile plots for model comparison with physiological and photosynthetic traits under well water and drought treatment. In this figure, black dots represent the standard line means that the predicted value equal to the observed value; yellow dots represent the simple linear model (without population structure and relative kinship); green dots represents the Q model; blue dots represents the K model; and red dots represent the Q+K model. (A) Evaluation of model types using markers for LWC. (B) Evaluation of model types using markers for MDA. (C) Evaluation of model types using markers for Pn. (D) Evaluation of model types using markers for Gs. (E) Evaluation of model types using markers for Ci. (F) Evaluation of model types using markers for Tn.
Molecular markers have been widely used to analysis the genetic relationship of Miscanthus all around its distributed areas (Xu et al., 2013; Clark et al., 2014; Yook et al., 2014; Anzoua et al., 2015; Nie et al., 2016). Although some genetic maps (Kim et al., 2012; Ma et al., 2012; Swaminathan et al., 2012; Liu et al., 2015) and association mapping (Zhao et al., 2011; Slavov et al., 2014; Nie et al., 2016) with high density and resolution have been constructed, mapping of QTLs using the genetic map are limited owing to its absent of genomic information. Zhao et al. (2013) used 23 SSR markers transferred from Brachypodium distachyon to conduct marker-trait association analysis and 115 loci were generated. Among them, nine markers were detected to be significantly (P < 0.01) associated with heading date and biomass yield. Slavov et al. (2014) conducted genome-wide association studies in a 138 M. sinensis population by using single-nucleotide variants markers and a total of 17 significant associations were detected with phenology, morphology, and cell wall composition traits. Nie et al. (2016) used 104 pair of markers to conducted the marker-trait association analysis and among 1059 generated loci, 12 significant associations of biomass yield related traits were identified. In this study, among 110 loci generated from 30 randomly selected SSRs, 7 significant associations were identified with RWC and photosynthetic trait under drought stress indicating that the developed SSR markers have a relative stable linked efficiency to the drought tolerance related traits.
By searching for the function annotation of 7 unigenes accounting for the associated markers in the protein database (Supplementary Table S2), an definite answer was appeared. Especially, unigene CL10592.Contig1_All_5404_1 and CL10112.Contig3_All_5301_4 associated with Pn and Gs were annotated function related to vegetative cell wall protein, while unigene CL12059.Contig2_ All_5665_1 and CL10784.Contig3_All_5453_1 associated with RWC were annotated related to arginine/serine-rich protein and ubiquitin-like-specific protease. With the abundant genetic basis information, the developed markers could be potential candidates for improving Miscanthus selection for enhanced drought tolerance in a breeding program.
XiZ and LH conceived the project and designed the experiments; GN, LT, YZ, and XC performed the experiments; GN, XM, and XuZ analyzed the data; GN, LT, and LP finalized the manuscript.
This work was supported by the Earmarked Fund for the Modern Agro-Industry Technology Research System (#CARS-35-05), the National High-Technology Research and Development Program (863 Program) of China (No. 2012AA101 801-02).
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 sincerely thank Prof. Cuixia Chen from Shandong Agricultural University, Prof. Jianxiong Jiang from Hu’nan Agricultural University, and Chinese Germplasm Seed Bank for providing part of materials used in this study.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2017.00801/full#supplementary-material
ABA, Abscisic acid; ANOVA, Analysis of variance; Ci, Intercellular CO2 concentration; DW, Dry weight; FW, Fresh weight; GLM, General linear model; GO, Gene ontology; Gs, Stomatal conductance; KEGG, Kyoto encyclopedia of genes and genomes; MAS, Marker-assisted selection; MCMC, Markov chain monte carlo; MDA, Malondialdehyde; MLM, Mixed linear model; NGS, Next-generation sequencing; NR, Non-redundant protein sequence; PAGE, Polyacrylamide gel electrophoresis; Pn, Net photosynthetic rate; RWC, Relative water content; SAHN, Sequential agglomerative hierarchical and nested clustering; SSRs, Simple sequence repeats; SW, Saturated weight; SWC, Soil water content; Tn, Transpiration rate; Swiss-Prot, Annotated protein sequence; UPGMA, Unweighted pair group method with arithmetic means; WUE, Water-use efficiency.
Anzoua, K. G., Suzuki, K., Fujita, S., Toma, Y., and Yamada, T. (2015). Evaluation of morphological traits, winter survival and biomass potential in wild Japanese Miscanthus sinensis Anderss. Populations in northern Japan. Grassland Sci. 61, 83–91. doi: 10.1111/grs.12085
Bradbury, P. J., Zhang, Z., Kroon, D. E., Casstevens, T. M., Ramdoss, Y., and Buckler, E. S. (2007). TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics 23, 2633–2635. doi: 10.1093/bioinformatics/btm308
Byrt, C., Grof, C., and Furbank, R. (2011). C4 plants as biofuel feedstocks: optimising biomass production and feedstock quality from a lignocellulosic perspective. J. Integr. Plant Biol. 53, 120–135. doi: 10.1111/j.1744-7909.2010.01023.x
Cattivelli, L., Rizzaa, F., Badeck, F., Mazzucotelli, E., Mastrangelo, A., Francia, E., et al. (2008). Drought tolerance improvement in crop plants: an integrated view from breeding to genomics. Field Crops Res. 105, 1–14. doi: 10.1016/j.fcr
Centritto, M., Lauteri, M., Monteverdi, M. C., and Serraj, R. (2009). Leaf gas exchange, carbon isotope discrimination, and grain yield in contrasting rice genotypes subjected to water deficits during the reproductive stage. J. Exp. Bot. 60, 2325–2339. doi: 10.1093/jxb/erp123
Chae, W. B., Hong, S. J., Gifford, J. M., Rayburn, A. L., Sacks, E. J., and Juvik, J. A. (2014). Plant morphology, genome size, and SSR markers differentiate five distinct taxonomic groups among accessions in the genus Miscanthus. Glob. Change Biol. Bioenergy 6, 646–660. doi: 10.1111/gcbb.12101
Clark, L. V., Brummer, J. E., Glowacka, K., Hall, M. C., Heo, K., Peng, J., et al. (2014). A footprint of past climate change on the diversity and population structure of Miscanthus sinensis. Ann. Bot. 114, 97–107. doi: 10.1093/aob/mcu084
Clifton-Brown, J., and Lewandowski, I. (2000). Water use efficiency and biomass partitioning of three different Miscanthus genotypes with limited and unlimited water supply. Ann. Bot. 86, 191–200. doi: 10.1006/anbo.2000.1183
Clifton-Brown, J. C., and Lewandowski, I. (2002). Screening Miscanthus genotypes in field trials to optimise biomass yield and quality in Southern Germany. Eur. J. Agron. 16, 97–110. doi: 10.1016/S1161-0301(01)00120-4
Costa, J. M., Ortuño, M. F., Lopes, C. M., and Chaves, M. M. (2012). Grapevine varieties exhibiting differences in stomatal response to water deficit. Funct. Plant Biol. 39, 179–189. doi: 10.1071/fp11156
Dwiyanti, M. S., Stewart, J. R., Nishiwaki, A., and Yamada, T. (2014). Natural variation in Miscanthus sinensis seed germination under low temperatures. Grassland Sci. 60, 194–198. doi: 10.1111/Grs.12051
Ellis, R. P., Forster, B. P., Robinson, D., Handley, L. L., Gordon, D. C., Russell, J. R., et al. (2000). Wild barley: a source of genes for crop improvement in the 21st century? J. Exp. Bot. 51, 9–17. doi: 10.1093/jexbot/51.342.9
Evanno, G., Regnaut, S., and 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
Gilles, A., Meglecz, E., Pech, N., Ferreira, S., Malausa, T., and Martin, J. F. (2011). Accuracy and quality assessment of 454 GS-FLX Titanium pyrosequencing. BMC Genomics 12:245. doi: 10.1186/1471-2164-12-245
Greef, J. M., Deuter, M., Jung, C., and Schondelmaier, J. (1997). Genetic diversity of European Miscanthus species revealed by AFLP fingerprinting. Genet. Resour. Crop Evol. 44, 185–195. doi: 10.1023/A:1008693214629
Gupta, S., and Prasad, M. (2009). Development and characterization of genic SSR markers in Medicago truncatula and their transferability in leguminous and non-leguminous species. Genome 52, 761–771. doi: 10.1139/g09-051
Hastings, A., Clifton-Brown, J., Wattenbach, M., Mitchell, P., and Smith, P. (2009). The development of MISCANFOR, a new Miscanthus crop growth model: towards more robust yield predictions under different climatic and soil conditions. Glob. Change Biol. Bioenergy 1, 154–170. doi: 10.1111/j.1757-1707.2009.01007.x
Heaton, E. A., Dohleman, F. G., and Long, S. P. (2009). Seasonal nitrogen dynamics of Miscanthus x giganteus and Panicum virgatum. Glob. Change Biol. Bioenergy 1, 297–307. doi: 10.1111/j.1757-1707.2009.01022.x
Hernandez, P., Dorado, G., Laurie, D. A., Martin, A., and Snape, J. W. (2001). Microsatellites and RFLP probes from maize are efficient sources of molecular markers for the biomass energy crop Miscanthus. Theor. Appl. Genet. 102, 616–622. doi: 10.1007/s001220051688
Ho, C. W., Wu, T. H., Hsu, T. W., Huang, J. C., Huang, C. C., and Chiang, T. Y. (2011). Development of 12 genic microsatellite loci for a biofuel grass, Miscanthus sinensis (Poaceae). Am. J. Bot. 98, e201–e203. doi: 10.3732/ajb.1100071
Hodkinson, T. R., Chase, M. W., Lledo, M. D., Salamin, N., and Renvoize, S. A. (2002). Phylogenetics of Miscanthus, Saccharum and related genera (Saccharinae, Andropogoneae, Poaceae) based on DNA sequences from ITS nuclear ribosomal DNA and plastid trnL intron and trnL-F intergenic spacers. J. Plant Res. 115, 381–392. doi: 10.1007/s10265-002-0049-3
Huang, X., Yan, H. D., Zhang, X. Q., Zhang, J., Frazier, T. P., Huang, D. J., et al. (2016). De novo transcriptome analysis and molecular marker development of two Hemarthria species. Front. Plant Sci. 7:496. doi: 10.3389/fpls.2016.00496
Hung, K. H., Chiang, T. Y., Chiu, C. T., Hsu, T. W., and Ho, C. W. (2009). Isolation and characterization of microsatellite loci from a potential biofuel plant Miscanthus sinensis (Poaceae). Conserv. Genet. 10, 1377–1380. doi: 10.1007/s10592-008-9699-6
Ings, J., Mur, L. A., Robson, P. R., and Bosch, M. (2013). Physiological and growth responses to water deficit in the bioenergy crop Miscanthus x giganteus. Front. Plant Sci. 4:468. doi: 10.3389/fpls.2013.00468
Jansen, M., Gilmer, F., Biskup, B., Nagel, K. A., Rascher, A., Briem, S., et al. (2009). Simultaneous phenotyping of leaf growth and chlorophyll fluorescence via GROWSCREEN FLUORO allows detection of stress tolerance in Arabidopsis thaliana and other rosette plants. Funct. Plant Biol. 36, 902–914. doi: 10.1071/Fpo9095
Jiang, Y., Liu, H., and Cline, V. (2009). Correlations of leaf relative water content, canopy temperature, and spectral reflectance in perennial ryegrass under water deficit conditions. HortScience 44, 459–462.
Kane, N. C., and Rieseberg, L. H. (2007). Selective sweeps reveal candidate genes for adaptation to drought and salt tolerance in common sunflower. Helianthus annuus. Genetics 175, 1823–1834. doi: 10.1534/genetics.106.067728
Kantety, R. V., La Rota, M., Matthews, D. E., and Sorrells, M. E. (2002). Data mining for simple sequence repeats in expressed sequence tags from barley, maize, rice, sorghum and wheat. Plant Mol. Biol. 48, 501–510. doi: 10.1023/a:1014875206165
Kaur, S., Pembleton, L. W., Cogan, N. O., Savin, K. W., Leonforte, T., Paull, J., et al. (2012). Transcriptome sequencing of field pea and faba bean for discovery and validation of SSR genetic markers. BMC Genomics 13:104. doi: 10.1186/1471-2164-13-104
Kim, C., Zhang, D., Auckland, S. A., Rainville, L. K., Jakob, K., Kronmiller, B., et al. (2012). SSR-based genetic maps of Miscanthus sinensis and M. sacchariflorus, and their comparison to sorghum. Theor. Appl. Genet. 124, 1325–1338. doi: 10.1007/s00122-012-1790-1
Liu, S., Clark, L. V., Swaminathan, K., Gifford, J. M., Juvik, J. A., and Sacks, E. J. (2015). High-density genetic map of Miscanthus sinensis reveals inheritance of zebra stripe. GCB Bioenergy 8, 616–630. doi: 10.1111/gcbb.12275
Ma, X. F., Jensen, E., Alexandrov, N., Troukhan, M., Zhang, L., Thomas-Jones, S., et al. (2012). High resolution genetic mapping by genome sequencing reveals genome duplication and tetraploid genetic structure of the diploid Miscanthus sinensis. PLoS ONE 7:e33821. doi: 10.1371/journal.pone.0033821
Malinowska, M., Donnison, I. S., and Robson, P. R. H. (2017). Phenomics analysis of drought responses in Miscanthus collected from different geographical locations. GCB Bioenergy 9, 78–91. doi: 10.1111/gcbb.12350
Maughan, M., Bollero, G., Lee, D. K., Darmody, R., Bonos, S., Cortese, L., et al. (2012). Miscanthus × giganteus productivity: the effects of management in different environments. GCB Bioenergy 4, 253–265. doi: 10.1111/j.1757-1707.2011.01144.x
Mittova, V., Guy, M., Tal, M., and Volokita, M. (2004). Salinity up-regulates the antioxidative system in root mitochondria and peroxisomes of the wild salt-tolerant tomato species Lycopersicon pennellii. J. Exp. Bot. 55, 1105–1113. doi: 10.1093/jxb/erh113
Nie, G., Huang, L., Zhang, X., Taylor, M., Jiang, Y., Yu, X., et al. (2016). Marker-trait association for biomass yield of potential bio-fuel feedstock Miscanthus sinensis from Southwest China. Front. Plant Sci. 7:802. doi: 10.3389/fpls.2016.00802
Nie, G., Zhang, X. Q., Huang, L. K., Xu, W. Z., Wang, J. P., Zhang, Y. W., et al. (2014). Genetic variability and population structure of the potential bioenergy crop Miscanthus sinensis (Poaceae) in Southwest China based on SRAP markers. Molecules 19, 12881–12897. doi: 10.3390/molecules190812881
Nishiwaki, A., Mizuguti, A., Kuwabara, S., Toma, Y., Ishigaki, G., Miyashita, T., et al. (2011). Discovery of natural Miscanthus (Poaceae) triploid plants in sympatric populations of Miscanthus sacchariflorus and Miscanthus sinensis in southern Japan. Am. J. Bot. 98, 154–159. doi: 10.3732/ajb.1000258
O’Neill, P., Shanahan, J., and Schepers, J. (2006). Use of chlorophyll fluorescence assessments to differentiate corn hybrid response to variable water conditions. Crop Sci. 46, 681–687. doi: 10.2135/cropsci2005.06-0170
Quinn, L. D., Straker, K. C., Guo, J., Kim, S., Thapa, S., Kling, G., et al. (2015). Stress-tolerant feedstocks for sustainable bioenergy production on marginal land. Bioenergy Res. 8, 1081–1100. doi: 10.1007/s12155-014-9557-y
Slavov, G. T., Nipper, R., Robson, P., Farrar, K., Allison, G. G., Bosch, M., et al. (2014). Genome-wide association studies and prediction of 17 traits related to phenology, biomass and cell wall composition in the energy grass Miscanthus sinensis. New Phytol. 201, 1227–1239. doi: 10.1111/nph.12621
Swaminathan, K., Chae, W. B., Mitros, T., Varala, K., Xie, L., Barling, A., et al. (2012). A framework genetic map for Miscanthus sinensis from RNAseq-based markers shows recent tetraploidy. BMC Genomics 13:142. doi: 10.1186/1471-2164-13-142
Temnykh, S., Park, W. D., Ayres, N., Cartinhour, S., Hauck, N., Lipovich, L., et al. (2000). Mapping and genome organization of microsatellite sequences in rice (Oryza sativa L.). Theor. Appl. Genet. 100, 697–712. doi: 10.1007/s001220051342
Tian, W., Paudel, D., Vendrame, W., and Wang, J. (2017). Enriching genomic resources and marker development from transcript sequences of Jatropha curcas for microgravity studies. Int. J. Genomics 2017:8614160. doi: 10.1155/2017/8614160
van der Weijde, T., Alvim Kamei, C. L., Torres, A. F., Vermerris, W., Dolstra, O., Visser, R. G., et al. (2013). The potential of C4 grasses for cellulosic biofuel production. Front. Plant Sci. 4:107. doi: 10.3389/fpls.2013.00107
van der Weijde, T., Huxley, L. M., Hawkins, S., Sembiring, E. H., Farrar, K., Dolstra, O., et al. (2016). Impact of drought stress on growth and quality of Miscanthus for biofuel production. GCB Bioenergy 9, 770–782. doi: 10.1111/gcbb.12382
Varshney, R. (2009). “Gene-based marker systems in plants: high throughput approaches for marker discovery and genotyping,” in Molecular Techniques in Crop Improvement, eds S. M. Jain and D. S. Brar (Dordrecht: Springer), 119–142.
Xu, W. Z., Zhang, X. Q., Huang, L. K., Nie, G., and Wang, J. P. (2013). Higher genetic diversity and gene flow in wild populations of Miscanthus sinensis in southwest China. Biochem. Syst. Ecol. 48, 174–181. doi: 10.1016/j.bse.2012.11.024
Yates, S. A., Swain, M. T., Hegarty, M. J., Chernukin, I., Lowe, M., Allison, G. G., et al. (2014). De novo assembly of red clover transcriptome based on RNA-Seq data provides insight into drought response, gene discovery and marker identification. BMC Genomics 15:453. doi: 10.1186/1471-2164-15-453
Yook, M. J., Lim, S. H., Song, J. S., Kim, J. W., Zhang, C. J., Lee, E. J., et al. (2014). Assessment of genetic diversity of Korean Miscanthus using morphological traits and SSR markers. Biomass Bioenergy 66, 81–92. doi: 10.1016/j.biombioe.2014.01.025
Yu, X., Bai, G., Liu, S., Luo, N., Wang, Y., Richmond, D. S., et al. (2013). Association of candidate genes with drought tolerance traits in diverse perennial ryegrass accessions. J. Exp. Bot. 64, 1537–1551. doi: 10.1093/jxb/ert018
Yu, X., Bai, G., Luo, N., Chen, Z., Liu, S., Liu, J., et al. (2011). Association of simple sequence repeat (SSR) markers with submergence tolerance in diverse populations of perennial ryegrass. Plant Sci. 180, 391–398. doi: 10.1016/j.plantsci.2010.10.013
Zhao, H., Wang, B., He, J., Yang, J., Pan, L., Sun, D., et al. (2013). Genetic diversity and population structure of Miscanthus sinensis germplasm in China. PLoS ONE 8:e75672. doi: 10.1371/journal.pone.0075672
Zhao, H., Yu, J., You, F. M., Luo, M., and Peng, J. (2011). Transferability of microsatellite markers from Brachypodium distachyon to Miscanthus sinensis, a potential biomass crop. J. Integr. Plant Biol. 53, 232–245. doi: 10.1111/j.1744-7909.2010.01026.x
Zhao, X. X., Huang, L. K., Zhang, X. Q., Li, Z., and Peng, Y. (2014). Effects of heat acclimation on photosynthesis, antioxidant enzyme activities, and gene expression in orchardgrass under heat stress. Molecules 19, 13564–13576. doi: 10.3390/molecules190913564
Keywords: bio-fuel, RNA sequencing, markers development, association analysis, Miscanthus sinensis
Citation: Nie G, Tang L, Zhang Y, Huang L, Ma X, Cao X, Pan L, Zhang X and Zhang X (2017) Development of SSR Markers Based on Transcriptome Sequencing and Association Analysis with Drought Tolerance in Perennial Grass Miscanthus from China. Front. Plant Sci. 8:801. doi: 10.3389/fpls.2017.00801
Received: 05 March 2017; Accepted: 28 April 2017;
Published: 16 May 2017.
Edited by:Zhulong Chan, Huazhong Agricultural University, China
Copyright © 2017 Nie, Tang, Zhang, Huang, Ma, Cao, Pan, Zhang and Zhang. 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: Xinquan Zhang, email@example.com
†These authors have contributed equally to this work.