Impact Factor 3.678

The world's most-cited Plant Sciences journal

Original Research ARTICLE

Front. Plant Sci., 08 August 2016 | https://doi.org/10.3389/fpls.2016.01193

High-Throughput Sequencing Reveals Single Nucleotide Variants in Longer-Kernel Bread Wheat

Feng Chen*, Zibo Zhu, Xiaobian Zhou, Yan Yan, Zhongdong Dong and Dangqun Cui
  • National Key Laboratory of Wheat and Maize Crop, Collaborative Innovation Center of Henan Grain Crops, Agronomy College, Henan Agricultural University, Zhengzhou, China

The transcriptomes of bread wheat Yunong 201 and its ethyl methanesulfonate derivative Yunong 3114 were obtained by next-sequencing technology. Single nucleotide variants (SNVs) in the wheat strains were explored and compared. A total of 5907 and 6287 non-synonymous SNVs were acquired for Yunong 201 and 3114, respectively. A total of 4021 genes with SNVs were obtained. The genes that underwent non-synonymous SNVs were significantly involved in ATP binding, protein phosphorylation, and cellular protein metabolic process. The heat map analysis also indicated that most of these mutant genes were significantly differentially expressed at different developmental stages. The SNVs in these genes possibly contribute to the longer kernel length of Yunong 3114. Our data provide useful information on wheat transcriptome for future studies on wheat functional genomics. This study could also help in illustrating the gene functions of the non-synonymous SNVs of Yunong 201 and 3114.

Introduction

Wheat (Triticum aestivum L.) is one of the three most important cereals (i.e., maize, rice, and wheat), with more than 600 million tons harvested annually (Shewry, 2009). In wheat breeding programs, yield can be divided into three components, namely, spike number per acre, grain number per spike, and thousand-grain weight. Grain weight is mainly determined by grain size in bread wheat. Therefore, kernel size, which is extensively studied in wheat breeding programs, is a key factor affecting wheat yield. Grain size is mainly controlled by heredity. Previously, many genes related to grain size of rice have been successfully cloned (Song et al., 2007; Wang et al., 2008, 2012; Weng et al., 2008; Li et al., 2011). Some quantitative trait loci (QTLs) controlling grain shape and size in bread wheat have been mapped through marker analysis in different wheat populations (Breseghello and Sorrells, 2006; Wu et al., 2015). Many of them were associated with the orthologs of rice grain traits QTLs (Lu et al., 2015; Wu et al., 2015; Zanke et al., 2015). For example, GRAIN SIZE 3 (GS3) is a major QTL for rice grain length and weight (Fan et al., 2006), TaGS-D1 was the syntenic gene in wheat (Zhang et al., 2014). These studies revealed that grain size of bread wheat is regulated via a complex molecular genetic mechanism. However, in bread wheat, the genetic study of grain size has been limited to date because of its large de novo genome. Next-generation sequencing (NGS) technology (Goff et al., 2002) provides a novel method to identify, map, and quantify transcriptomes (Kyndt et al., 2012); this method can also be used for rapid characterization of transcript sequences, gene expression (Wang et al., 2009), and genomic variation in polyploid plants with de novo genomes (Goff et al., 2002; Marioni et al., 2008). The first homolog-specific sequence assembly of wheat transcriptome is based on Roche 454 and Illumina GAIIx (Schreiber et al., 2012). Moreover, the genome of wheat and of its relative were analyzed by recent studies utilizing high-throughput sequencing, which provided references for further study (Brenchley et al., 2012; Ling et al., 2013; The International Wheat Genome Sequencing Consortium [IWGSC], 2014). Subsequently, large amount of advances have been acquired on wheat transcriptome (Duan et al., 2012; Pfeifer et al., 2014; Tanaka et al., 2014; Pingault et al., 2015).

The correlation between the genetic variants and the phenotypes is still a central question for crop improvements. Ethyl methanesulfonate (EMS) induces a large spectrum of mutations, including truncations and missense mutations, thereby allowing to be a readily used chemical in traditional breeding programs because of its flexibility, non-transgenicity, and stable inheritability (McCallum et al., 2000; Henikoff et al., 2004). Furthermore, EMS can create random point mutations at high density in polyploid plants. Many studies has focused on creation of mutations via EMS in specific genes of plants, and all the mutations are G–A or C–T transitions in bread wheat (Feiz et al., 2009; Uauy et al., 2009; Slade et al., 2012; Wang et al., 2014). Through RNAseq, large number of genetic variants across the transcriptomes could be identified. The possible transcriptional mechanism of trait regulation due to genetic variants including single nucleotide polymorphism substitutions has been reported in many crops such as in rice (Mao et al., 2010), but still insufficient in wheat.

A Chinese wheat cultivar Yunong 201, which was released in 2006 (No. Yushenmai2006006), is a high-quality noodle wheat strain that is disease-resistant. Meanwhile, its EMS mutagenesis-derived Yunong 3114 shows longer kernel length and higher production. In the present study, to better understand the genetic basis of kernel size in bread wheat, Illumina (Solexa) sequencing technology was applied in bread wheat Yunong 201 and its EMS mutant line Yunong 3114 to generate their transcriptomes. Gene profiles of Yunong 201 and 3114 were obtained by de novo sequencing. Single nucleotide variants (SNVs) were analyzed. This study provided important information to further understand the transcriptome of hexaploid wheat and determine the wheat-specific genes related to grain size.

Materials and Methods

Plant Materials

A Chinese winter wheat cultivar Yunong 201 showing outstanding dry, white noodle quality was treated by 1.0% EMS (0.1 mol/L Na2HPO4 ⋅12H2O, pH 7.0). An elite M2 line was screened from the EMS-mutagenized population encompassing 2000 lines because of its longer kernel length and higher grain weight; this line was self-crossed for three times into a M5 line Yunong 3114. Yunong 201 and its derived line Yunong 3114 were planted and grown at the Zhengzhou Scientific Research and Education Center of Henan Agricultural University during 2011–2012 cropping seasons under non-stressed conditions. Grains were collected at different developmental stages from both lines at 7, 14, 21, 28, and 35 days after flowering, and mature seeds (stored 1 year at 4°C after harvested); pooled samples were obtained by mixing eighteen grains from three spikes of each wheat strain at each developmental stage mentioned above; these samples were used for the generation of small RNA libraries of Yunong 201 and 3114.

RNA Sequencing

The total RNA of each sample was extracted with Trizol Reagent (Invitrogen), according to the manufacturer’s instructions. Full-length cDNAs were synthesized with a TruSeq RNA Sample Preparation Kit (Illumina), according to the manufacturer’s protocol. The cDNA libraries were sequenced on Illumina Hiseq 2000, according to the manufacturer’s instructions. A 10 Gb raw reads were obtained from each cDNA library, and a total of 120 Gb raw reads were generated from the different developmental staged seeds of Yunong 201 and 3114.

Pretreatment of Data and Development of Unigenes

Raw reads were first trimmed by using SeqClean1 86–64, Newbler 2.5.3 (Kumar and Blaxter, 2010), and Lucy 1.20p (Chou and Holmes, 2001) to remove low-quality reads, vector sequences, and reads whose length was less than 50 bp. The trimmed reads were subjected to polygenetic analysis with CD-HIT version 4.0 (Huang et al., 2010) at a sequence identity threshold of 99%. Possible contaminated sequences were removed by aligning to the NT library. The trimmed non-redundant reads were assembled to contigs by Newbler 2.5.3. One or more contigs were assembled to isotigs. A total of 12 groups of isotigs were developed from the small RNA libraries of Yunong 201 and 3114.

Discovery of SNV

All SNVs of Yunong 201 and 3114 were obtained based on the wheat reference sequence mapping (URGI)2 (He et al., 2014). We focused on the sites with non-synonymous mutation. EMS treatment can trigger high-density point mutation on wheat material. We selected the mutation sites occurring in all samples. To explore the different biological mechanisms between Yunong 201 and 3114, three SNV types were significant. The first SNV type was that the mutation sites were same between Yunong 201 and 3114, and the mutation was non-synonymous as compared with reference sequence, but the nucleotides after mutation were different between Yunong 201 and 3114, we considered this mutation significant. The second type was that the nucleotides of Yunong 201 were changed into other nucleotides, but the nucleotides of Yunong 3114 were still the same as the reference sequence. The third SNV indicated that mutation occurred specifically in Yunong 3114, but not in Yunong 201. Thus, we studied the same and different mutation sites between Yunong 201 and Yunong 3114.

Sequence Annotation

Isotigs and singlets sequences were subjected to polygenetic analysis by CD-HIT version 4.0 (Huang et al., 2010) with an identity of 95%. Chromosome mapping and comparison with known genome segments were also performed. Using the BLAST program (Camacho et al., 2009) (a searched threshold of 1e–10), we searched for unigenes in EST libraries of 10 species related to NCBI3 and DFCI4. Nucleotides and protein were separately annotated by comparing the unigenes with those in the NT library (E-value < 1e–5) and the NR library (BLASTX; E-value < 1e–5; similarity of protein >30%). Moreover, we also annotated the sequences based on the wheat sequence database of French National Institute for Agricultural Research5 (INRA). We mainly referenced the sequence of T. aestivum chromosome 3B (ta3b, genomic scaffold, cultivar Chinese Spring) in URGI.

Functional Annotation

We identified the annotated unigene sequences for the possible functions through the Gene Ontology Consortium6. Sequences were mapped to the reference authoritative pathways in Kyoto Encyclopedia of Genes and Genomes (KEGG)7 to determine the active biological pathways in annotated unigene sequences.

Heat Map Analysis

Heat map is shown for preliminary SNV gene expression (logFoldChange) in the mutant strain at different developmental stages with respect to the wild type. Heat map was generated via Microsoft Excel, and logFold-Change values were used. The highest, medium, and lowest values are shown in color red, black, and green, respectively.

Results

Investigation of Agronomic Traits in Yunong 201 and 3114

Wheat strain Yunong 3114 showed longer kernel length, higher grain weights, and higher yield compared with Yunong 201 (data not shown); the former did not also show obvious difference on kernel width (detailed data in Figure 1). Therefore, the higher grain weight of Yunong 3114 resulted from its relatively longer kernel length than that of Yunong 201. The grains with different developmental stages at 7, 14, 21, 28, and 35 days after flowering, and mature seeds (stored 1 year at 4°C after harvested) were collected from both lines; these grains were used to generate six pairs of small RNA libraries of Yunong 201 and 3114 for sequencing. Transcriptomes of both Yunong 201 and 3114 were obtained for the SNVs analysis.

FIGURE 1
www.frontiersin.org

FIGURE 1. Kernel size parameters of Yunong 201 and 3114 during six seed developmental stages. DAF, days after flowering.

Discovery of SNVs in Yunong 201 and 3114

As the wheat reference sequence in URGI (Unité de Recherche Génomique Info)8 was mapped, SNVs in both Yunong 201 and 3114 were detected. A large number of genetic variations were identified across the transcriptomes of both strains, including non-synonymous, synonymous, stop–gain, stop–loss, and stop–stop SNVs (Figure 2). To maintain data credibility, the non-synonymous SNVs that could be simultaneously annotated in reference sequence were analyzed. Non-synonymous SNVs accounted for 49.21% (5904) and 50.32% (6287) in Yunong 201 and 3114 of overall variations occurring in CDS regions, respectively (Table 1). A total of 7432 same mutated positions in Yunong 201 and 3114 were also observed, including 3714 synonymous SNVs and 3718 non-synonymous SNVs. Yunong 3114 had 2569 distinct non-synonymous SNVs, and Yunong 201 had 2186 distinct non-synonymous SNVs as compared to reference sequence (Table 1). Furthermore, one SNV site may induce both synonymous and non-synonymous mutations that were acted as the multi-SNVs. A total of 22 and 18 multi-SNVs in Yunong 201 and 3114 were recorded, respectively.

FIGURE 2
www.frontiersin.org

FIGURE 2. Single nucleotide variants (SNVs) between Yunong 201 and 3114. (A–C) The total SNV, non-synonymous, and synonymous SNVs in Yunong 201 and 3114. (D–F) Three specific kinds of SNVs, namely, stop–gain, stop–loss, and stop–stop.

TABLE 1
www.frontiersin.org

TABLE 1. Single nucleotide variants in Yunong 201 and 3114 based on reference sequence.

Ethyl methanesulfonate preferentially induces G residues alkylation, thereby resulting in many G→A and C→T transitions [CG→TA; (Greene et al., 2003; Till et al., 2007; Tsai et al., 2011; Monson-Miller et al., 2012)]. Thus, in accordance with previous findings, C→T and G→A transitions were the most prevalent patterns of nucleotide substitutions (Figure 3A; Supplementary Figure S1), and the total GC content of Yunong 201 was higher than that of Yunong 3114. Based on the reference wheat sequence, the transitions were higher than transversions both in Yunong 201 and in 3114 (Figures 3B,C). Moreover, the mutations in the first and second codes were significantly more frequent than the third code, both in Yunong 201 and 3114. Our data indicated that the genetic variations were significantly different between the Yunong wheat strains and reference sequence; meanwhile, the difference was smaller between Yunong 201 and 3114, which is mainly because Yunong 3114 was derived from Yunong 201.

FIGURE 3
www.frontiersin.org

FIGURE 3. Patterns of nucleotide substitution for Yunong 201 and 3114 based on reference sequence. (A) Number of genes with different type of nucleotide substitution. (B,C) Percentage (%) of transition and transversion substitutions in Yunong 201 and Yuong 3114. Unpaired t-test was used, P-value < 0.0001.

SNVs and Gene Functional Analysis

Based on the reference sequence, a total of 4021 genes with SNVs were acquired, and most genes have one mutated position (Figure 4). A total of 1147 genes in Yunong 201 had one mutant position where 576 genes underwent a non-synonymous mutation (Figure 4). Meanwhile, a total of 1167 genes had one SNV positions in Yunong 3114, 564 of which had a non-synonymous mutation. We also detected 4755 non-synonymous SNVs and 1168 genes with non-synonymous mutations by comparing Yunong 3114 with Yuong 201 (Table 2). Moreover, genes carrying transistions were more than those undergone to transversions (Figure 5). SNVs were distributed in different scaffolds and contigs of the reference genome. EMS-induced mutagenesis produces high mutant rate in plant. Therefore, we also analyzed the mutagenesis frequencies (Figure 6). However, the mutation sites were distributed atypically across the transcriptome (Figure 6; Supplementary Figure S2). The mutant gene numbers in different mutation patterns were determined, as shown in Supplementary Figure S3.

FIGURE 4
www.frontiersin.org

FIGURE 4. Relationship between SNV and genes in Yunong 201 and 3114 based on the reference sequence. The x axis represents the number of SNV mutations on a single gene sequence. The y axis showed number of genes with SNVs. The number of genes with one SNV ranks as the top.

TABLE 2
www.frontiersin.org

TABLE 2. Single nucleotide variants and gene counts between Yunong 3114 and Yunong 201.

FIGURE 5
www.frontiersin.org

FIGURE 5. Mutation analysis in Yunong 3114 and 201. The x axis indicates the type of mutations; the y axis represents the number of mutant genes or number of SNVs.

FIGURE 6
www.frontiersin.org

FIGURE 6. Single nucleotide variant regions in the gene. A gene was equally divided into 20 regions, the percentage of mutant in each region was statistically calculated. The y axis showed the percentage of mutant; the x axis represented the gene regions.

All the mutant genes with G to A or C to T transitions in the mutant strain compared with Yunong 201 were extracted; function annotation and enrichment analyses were performed. These genes were annotated in ATP binding, zinc ion binding, and nucleotide binding (Figure 7A); they were also mainly enriched in components, such as intracellular membrane-bound and intracellular organelles (Figure 7B); moreover, their functions were enriched in cellular protein metabolic process, modification-dependent macromolecule catabolic process, and cellular protein catabolic process (Figure 7C). The preliminary expression levels of these mutant genes were also analyzed by heat map analysis (Figure 8). For all the six developmental stages, most of these mutant genes were significantly differentially expressed, which indicated that the SNVs on these genes might be associated with the longer kernel length and higher kernel weight of Yunong 3114.

FIGURE 7
www.frontiersin.org

FIGURE 7. Function annotation and enrichment for the genes with G→A and/or C→T mutations. (A) The annotated biological processes; (B) The enrichment of cellular component; (C) The enrichment of biological process.

FIGURE 8
www.frontiersin.org

FIGURE 8. Preliminary Expression of genes with mutations of G->A and/or C->T between Yunong 201 and Yunong 3114 in different days. Horizontal axis represent the genes and vertical axis represent the time. The preliminary expression level for a gene in some days was described by the grids, which was calculated by the logFC (FC stands for ‘fold change’). The change of color from green to red represents the expression value from the low to high.

Discussion

Next-generation sequencing has been applied in many species, specifically in crop breeding, to explore genomic variations for investigating population evolutionary history, discovering marker-trait linkages, and genetic cause of phenotype variation. EMS-induced mutagenesis is one of the most effective, reliable, powerful, and frequently used technologies in plant. In the present study, Illumina (Solexa) NGS technology was applied to obtain the transcriptomes of bread wheat strain Yunong 201 and its EMS mutant line Yunong 3114. We comprehensively analyzed the SNVs and their impact on genes in Yunong 3114. Our data demonstrated that transcriptomes could also be used for SNV analysis, and a set of SNVs had significant impact on the expression levels of the mutant genes, potentially associated with the superior traits such as the longer kernel length and higher kernel weight of Yunong 3114. However, further validations are necessary.

Non-synonymous SNVs in Yunong 3114 and 201

Non-synonymous changes directly determine the amino acid change and further affect protein structure, stability, or location. Thus, the identification of non-synonymous sequence variations that can potentially affect gene or protein functionality is the focus of the current study. Stop–gain SNVs may lead to functional consequences due to protein truncation, transcript degradation, and dominant negative influences of protein species (Nagy and Maquat, 1998). The stop–loss non-synonymous SNVs of Yunong 3114 were also significantly less than that of Yunong 201. Current computational methods estimating loss-of-function in genes carrying variants are based on evolutionary conservation and functional redundancy on gene level (Rausell et al., 2014). For all non-synonymous SNVs, [A/G] and [T/C] are the most common base changes both in Yunong 201 and 3114, which agrees with the results of previous studies in legume species (Hand et al., 2008; Deulvot et al., 2010; Gaur et al., 2012; Leonforte et al., 2013). The nucleotide substitution in the first code in Yunong 3114 was significantly more frequent than in Yunong 201. Increased divergence in nucleotide composition shows a corresponding, predictable change in the amino acid compositions of the encoded proteins (Wang et al., 2004). Efficient transcription or mRNA processing is responsible for the high expression of GC-rich genes (Kudla et al., 2006). In early replicating regions, G and C nucleotides are more often misincorporated (Eyre-Walker and Hurst, 2001; Muyle et al., 2011). The distribution of GC content is a remarkable characteristic of genome organization that is often associated with many genomic features, such as meiotic recombination, gene density, and gene length (Mouchiroud et al., 1991; Duret et al., 1995; Duret and Arndt, 2008).

SNVs Impact on Gene Functionality

Most genes had one to three SNV positions both in Yunong 201 and 3114. The genes that underwent non-synonymous SNVs were significantly involved in ATP binding and metal and zinc ion binding; they specifically regulated protein kinase activity, protein phosphorylation, cellular protein metabolic process, and modification-dependent macromolecular catabolic process. Metal ions are crucial elements in enzymatic reaction events in all photosynthetic organisms, such as cyanobacteria, algae, and plants. Metal ions play important roles in maintaining substrate binding in the active site of metalloenzymes and controlling the redox activity of metalloenzymes in enzymatic reaction (Pospisil, 2014). Some oxidoreductases are also important in the metabolic responses. Protein phosphorylation and protein kinase activity play important roles in the regulation of carbon and nitrogen metabolite production in plants (Huber, 2007). Therefore, SNVs in Yunong 3114 may have impact on these metabolic processes, subsequently influencing the metabolism, energy or nutrient accumulation, and grain growth.

We also analyzed and compared the non-synonymous SNVs of the genes associated with these biological processes in both Yunong 201 and 3114. Some SNVs only occurred in Yunong 3114, and these SNVs potentially implicated the relatively superior trait. If the genomic architectures of the traits studied are highly complex, dozens or even hundreds of causative polymorphisms with minor effects may exist across the genome. For example, SNVs on genes of protein phosphatase 5 (PP5), glycerophosphoryl diester phosphodiesterase (GDPD), dephospho-CoA kinase, and alkylated DNA repair protein alkB homolog 1, cytochrome c oxidase subunit 2, and ABA responsive element binding factor may have impact on wheat gene function, thus may result in changes on biological process associated with grain size. It has been reported that Arabidopsis protein phosphatase 5 (AtPP5) performs multi enzymatic activities such as biologically active photoreceptor and protein phosphatase activity (Ryu et al., 2005). The AtPP5 is capable to form a complex with AtHsp90 under heat shock conditions, enhancing thermotolerance in Arabidopsis (Park et al., 2012). Arabidopsis SNC4 encodes a receptor-like kinase with two predicted extracellular GDPD domains and also plays a role as resistance protein (Bi et al., 2010). Another Arabidopsis protein shaven3 also has two tandemly repeated GDPD-like domains which were suggested to be involved in primary cell wall organization (Hayashi et al., 2008). The ABA responsive element binding factor has been found in rice such as OsABF1 (Oryza sativa ABA responsive element binding factor 1) gene which is involved in abiotic stress responses and ABA signaling in rice (Amir Hossain et al., 2010), and OREB1 whose phosphorylation plays important roles in regulating signal integration in the complex stress signaling network of plant cells (Hong et al., 2011). SNVs distributed in these genes might trigger alterations on the secondary structure of proteins that possibly impacts the subsequent functions (such as DNA binding and transcriptional activity) of target genes involved in multiple agronomic traits in wheat. The clues obtained from this study could facilitate the further illustration of gene function with the non-synonymous SNVs between Yunong 201 and 3114.

Conclusion

The current study analyzed the SNVs on transcriptome levels of bread wheat Yunong 201 and its EMS mutagenesis strain Yunong 3114. Our data provided useful information on wheat transcriptome by presenting new resources that could be used in future studies on wheat functional genomics. In addition, detecting new genetic variation for wheat and further information on gene functions that eventually generate new phenotypic variation provide insight into wheat breeding program.

Author Contributions

FC designed the project. FC, ZZ, and XZ performed RNA-seq experiments. FC wrote the paper. FC, ZD, YY, and DC performed the computational analyses.

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.

Acknowledgments

This project was funded by the 973 projects (2014CB160303 and 2014CB138105), Program for New Century Excellent Talents in University (NCET-13-0776) of China, and Program for Innovative Research Team (in Science and Technology) in University of Henan Province (14IRTSTHN010).

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2016.01193

FIGURE S1 | Patterns of nucleotide substitution in the SNVs of Yunong 201 and 3114 based on the reference genome. (A) Nucleotide substitutions of the SNVs for Yunong 201; (B) Nucleotide substitutions of the SNVs for Yunong 3114.

FIGURE S2 | Different SNVs of 18 functional genes between Yunong 201 and 3114. Boxes with letters show the mutant positions.

FIGURE S3 | Mutant gene numbers in different mutation patterns in both Yunong 201 and 3114.

Footnotes

  1. ^ http://sourceforge.net/projects/seqclean/
  2. ^ https://urgi.versailles.inra.fr/blast/blast.php
  3. ^ http://www.ncbi.nlm.nih.gov/
  4. ^ ftp://occams.dfci.harvard.edu/pub/bio/tgi/data
  5. ^ https://urgi.versailles.inra.fr/blast/blast.php
  6. ^ http://www.geneontology.org/
  7. ^ http://www.genome.jp/kegg/
  8. ^ https://urgi.versailles.inra.fr/

References

Amir Hossain, M., Lee, Y., Cho, J. I., Ahn, C. H., Lee, S. K., Jeon, J. S., et al. (2010). The bZIP transcription factor OsABF1 is an ABA responsive element binding factor that enhances abiotic stress signaling in rice. Plant Mol. Biol. 72, 557–566. doi: 10.1007/s11103-009-9592-9

CrossRef Full Text | Google Scholar

Bi, D., Cheng, Y. T., Li, X., and Zhang, Y. (2010). Activation of plant immune responses by a gain-of-function mutation in an atypical receptor-like kinase. Plant Physiol. 153, 1771–1779. doi: 10.1104/pp.110.158501

CrossRef Full Text | Google Scholar

Brenchley, R., Spannagl, M., Pfeifer, M., Barker, G. L., D’Amore, R., Allen, A. M., et al. (2012). Analysis of the bread wheat genome using whole-genome shotgun sequencing. Nature 491, 705–710. doi: 10.1038/nature11650

CrossRef Full Text | Google Scholar

Breseghello, F., and Sorrells, M. E. (2006). Association mapping of kernel size and milling quality in wheat (Triticum aestivum L.) cultivars. Genetics 172, 1165–1177. doi: 10.1534/genetics.105.044586

CrossRef Full Text | Google Scholar

Camacho, C., Coulouris, G., Avagyan, V., Ma, N., Papadopoulos, J., Bealer, K., et al. (2009). BLAST+: architecture and applications. BMC Bioinformatics 10:421. doi: 10.1186/1471-2105-10-421

CrossRef Full Text | Google Scholar

Chou, H. H., and Holmes, M. H. (2001). DNA sequence quality trimming and vector removal. Bioinformatics 17, 1093–1104. doi: 10.1093/bioinformatics/17.12.1093

CrossRef Full Text | Google Scholar

Deulvot, C., Charrel, H., Marty, A., Jacquin, F., Donnadieu, C., Lejeune-Henaut, I., et al. (2010). Highly-multiplexed SNP genotyping for genetic mapping and germplasm diversity studies in pea. BMC Genomics 11:468. doi: 10.1186/1471-2164-11-468

CrossRef Full Text | Google Scholar

Duan, J., Xia, C., Zhao, G., Jia, J., and Kong, X. (2012). Optimizing de novo common wheat transcriptome assembly using short-read RNA-Seq data. BMC Genomics 13:392. doi: 10.1186/1471-2164-13-392

CrossRef Full Text | Google Scholar

Duret, L., and Arndt, P. F. (2008). The impact of recombination on nucleotide substitutions in the human genome. PLoS Genet. 4:e1000071. doi: 10.1371/journal.pgen.1000071

CrossRef Full Text | Google Scholar

Duret, L., Mouchiroud, D., and Gautier, C. (1995). Statistical analysis of vertebrate sequences reveals that long genes are scarce in GC-rich isochores. J. Mol. Evol. 40, 308–317. doi: 10.1007/BF00163235

CrossRef Full Text | Google Scholar

Eyre-Walker, A., and Hurst, L. D. (2001). The evolution of isochores. Nat. Rev. Genet. 2, 549–555. doi: 10.1038/35080577

CrossRef Full Text | Google Scholar

Fan, C., Xing, Y., Mao, H., Lu, T., Han, B., Xu, C., et al. (2006). GS3, a major QTL for grain length and weight and minor QTL for grain width and thickness in rice, encodes a putative transmembrane protein. Theor. Appl. Genet. 112, 1164–1171. doi: 10.1007/s00122-006-0218-1

CrossRef Full Text | Google Scholar

Feiz, L., Martin, J. M., and Giroux, M. J. (2009). Creation and functional analysis of new Puroindoline alleles in Triticum aestivum. Theor. Appl. Genet. 118, 247–257. doi: 10.1007/s00122-008-0893-1

CrossRef Full Text | Google Scholar

Gaur, R., Azam, S., Jeena, G., Khan, A. W., Choudhary, S., Jain, M., et al. (2012). High-throughput SNP discovery and genotyping for constructing a saturated linkage map of chickpea (Cicer arietinum L.). DNA Res. 19, 357–373. doi: 10.1093/dnares/dss018

CrossRef Full Text | Google Scholar

Goff, S. A., Ricke, D., Lan, T. H., Presting, G., Wang, R., Dunn, M., et al. (2002). A draft sequence of the rice genome (Oryza sativa L. ssp. japonica). Science 296, 92–100. doi: 10.1126/science.1068275

CrossRef Full Text | Google Scholar

Greene, E. A., Codomo, C. A., Taylor, N. E., Henikoff, J. G., Till, B. J., Reynolds, S. H., et al. (2003). Spectrum of chemically induced mutations from a large-scale reverse-genetic screen in Arabidopsis. Genetics 164, 731–740.

Google Scholar

Hand, M. L., Ponting, R. C., Drayton, M. C., Lawless, K. A., Cogan, N. O., Charles Brummer, E., et al. (2008). Identification of homologous, homoeologous and paralogous sequence variants in an outbreeding allopolyploid species based on comparison with progenitor taxa. Mol. Genet. Genomics 280, 293–304. doi: 10.1007/s00438-008-0365-y

CrossRef Full Text | Google Scholar

Hayashi, S., Ishii, T., Matsunaga, T., Tominaga, R., Kuromori, T., Wada, T., et al. (2008). The glycerophosphoryl diester phosphodiesterase-like proteins SHV3 and its homologs play important roles in cell wall organization. Plant Cell Physiol. 49, 1522–1535. doi: 10.1093/pcp/pcn120

CrossRef Full Text | Google Scholar

He, Q., He, Q., Liu, X., Wei, Y., Shen, S., Hu, X., et al. (2014). Genome-wide prediction of cancer driver genes based on SNP and cancer SNV data. Am. J. Cancer Res. 4, 394–410.

Google Scholar

Henikoff, S., Till, B. J., and Comai, L. (2004). TILLING. Traditional mutagenesis meets functional genomics. Plant Physiol. 135, 630–636. doi: 10.1104/pp.104.041061

CrossRef Full Text | Google Scholar

Hong, J. Y., Chae, M. J., Lee, I. S., Lee, Y. N., Nam, M. H., Kim, D. Y., et al. (2011). Phosphorylation-mediated regulation of a rice ABA responsive element binding factor. Phytochemistry 72, 27–36. doi: 10.1016/j.phytochem.2010.10.005

CrossRef Full Text | Google Scholar

Huang, Y., Niu, B., Gao, Y., Fu, L., and Li, W. (2010). CD-HIT Suite: a web server for clustering and comparing biological sequences. Bioinformatics 26, 680–682. doi: 10.1093/bioinformatics/btq003

CrossRef Full Text | Google Scholar

Huber, S. C. (2007). Exploring the role of protein phosphorylation in plants: from signalling to metabolism. Biochem. Soc. Trans. 35(Pt 1), 28–32. doi: 10.1042/BST0350028

CrossRef Full Text | Google Scholar

Kudla, G., Lipinski, L., Caffin, F., Helwak, A., and Zylicz, M. (2006). High guanine and cytosine content increases mRNA levels in mammalian cells. PLoS Biol. 4:e180. doi: 10.1371/journal.pbio.0040180

CrossRef Full Text | Google Scholar

Kumar, S., and Blaxter, M. L. (2010). Comparing de novo assemblers for 454 transcriptome data. BMC Genomics 11:571. doi: 10.1186/1471-2164-11-571

CrossRef Full Text | Google Scholar

Kyndt, T., Denil, S., Haegeman, A., Trooskens, G., De Meyer, T., Van Criekinge, W., et al. (2012). Transcriptome analysis of rice mature root tissue and root tips in early development by massive parallel sequencing. J. Exp. Bot. 63, 2141–2157. doi: 10.1093/jxb/err435

CrossRef Full Text | Google Scholar

Leonforte, A., Sudheesh, S., Cogan, N. O., Salisbury, P. A., Nicolas, M. E., Materne, M., et al. (2013). SNP marker discovery, linkage map construction and identification of QTLs for enhanced salinity tolerance in field pea (Pisum sativum L.). BMC Plant Biol. 13:161. doi: 10.1186/1471-2229-13-161

CrossRef Full Text | Google Scholar

Li, Y., Fan, C., Xing, Y., Jiang, Y., Luo, L., Sun, L., et al. (2011). Natural variation in GS5 plays an important role in regulating grain size and yield in rice. Nat. Genet. 43, 1266–1269. doi: 10.1038/ng.977

CrossRef Full Text | Google Scholar

Ling, H. Q., Zhao, S., Liu, D., Wang, J., Sun, H., Zhang, C., et al. (2013). Draft genome of the wheat A-genome progenitor Triticum urartu. Nature 496, 87–90. doi: 10.1038/nature11997

CrossRef Full Text | Google Scholar

Lu, J., Chang, C., Zhang, H.-P., Wang, S.-X., Sun, G., Xiao, S.-H., et al. (2015). Identification of a novel allele of TaCKX6a02 associated with grain size, filling rate and weight of common wheat. PLoS ONE 10:e0144765. doi: 10.1371/journal.pone.0144765

CrossRef Full Text | Google Scholar

Mao, H., Sun, S., Yao, J., Wang, C., Yu, S., Xu, C., et al. (2010). Linking differential domain functions of the GS3 protein to natural variation of grain size in rice. Proc. Natl. Acad. Sci. U.S.A. 107, 19579–19584. doi: 10.1073/pnas.1014419107

CrossRef Full Text | Google Scholar

Marioni, J. C., Mason, C. E., Mane, S. M., Stephens, M., and Gilad, Y. (2008). RNA-seq: an assessment of technical reproducibility and comparison with gene expression arrays. Genome Res. 18, 1509–1517. doi: 10.1101/gr.079558.108

CrossRef Full Text | Google Scholar

McCallum, C. M., Comai, L., Greene, E. A., and Henikoff, S. (2000). Targeting induced local lesions IN genomes (TILLING) for plant functional genomics. Plant Physiol. 123, 439–442. doi: 10.1104/pp.123.2.439

CrossRef Full Text | Google Scholar

Monson-Miller, J., Sanchez-Mendez, D. C., Fass, J., Henry, I. M., Tai, T. H., and Comai, L. (2012). Reference genome-independent assessment of mutation density using restriction enzyme-phased sequencing. BMC Genomics 13:72. doi: 10.1186/1471-2164-13-72

CrossRef Full Text | Google Scholar

Mouchiroud, D., D’Onofrio, G., Aissani, B., Macaya, G., Gautier, C., and Bernardi, G. (1991). The distribution of genes in the human genome. Gene 100, 181–187. doi: 10.1016/0378-1119(91)90364-H

CrossRef Full Text | Google Scholar

Muyle, A., Serres-Giardi, L., Ressayre, A., Escobar, J., and Glemin, S. (2011). GC-biased gene conversion and selection affect GC content in the Oryza genus (rice). Mol. Biol. Evol. 28, 2695–2706. doi: 10.1093/molbev/msr104

CrossRef Full Text | Google Scholar

Nagy, E., and Maquat, L. E. (1998). A rule for termination-codon position within intron-containing genes: when nonsense affects RNA abundance. Trends Biochem. Sci. 23, 198–199. doi: 10.1016/S0968-0004(98)01208-0

CrossRef Full Text | Google Scholar

Park, J. H., Kim, W. Y., Chae, H. B., Kim, M. G., and Lee, S. Y. (2012). Serine/threonine protein phosphatase 5 (PP5) interacts with substrate under heat stress conditions and forms protein complex in Arabidopsis. Plant Signal. Behav. 7, 535–538. doi: 10.4161/psb.19699

CrossRef Full Text | Google Scholar

Pfeifer, M., Kugler, K. G., Sandve, S. R., Zhan, B., Rudi, H., Hvidsten, T. R., et al. (2014). Genome interplay in the grain transcriptome of hexaploid bread wheat. Science 345, 1250091. doi: 10.1126/science.1250091

CrossRef Full Text | Google Scholar

Pingault, L., Choulet, F., Alberti, A., Glover, N., Wincker, P., Feuillet, C., et al. (2015). Deep transcriptome sequencing provides new insights into the structural and functional organization of the wheat genome. Genome Biol. 16, 29. doi: 10.1186/s13059-015-0601-9

CrossRef Full Text | Google Scholar

Pospisil, P. (2014). The role of metals in production and scavenging of reactive oxygen species in photosystem II. Plant Cell Physiol. 55, 1224–1232. doi: 10.1093/pcp/pcu053

CrossRef Full Text | Google Scholar

Rausell, A., Mohammadi, P., McLaren, P. J., Bartha, I., Xenarios, I., Fellay, J., et al. (2014). Analysis of stop-gain and frameshift variants in human innate immunity genes. PLoS Comput. Biol. 10:e1003757. doi: 10.1371/journal.pcbi.1003757

CrossRef Full Text | Google Scholar

Ryu, J. S., Kim, J. I., Kunkel, T., Kim, B. C., Cho, D. S., Hong, S. H., et al. (2005). Phytochrome-specific type 5 phosphatase controls light signal flux by enhancing phytochrome stability and affinity for a signal transducer. Cell 120, 395–406. doi: 10.1016/j.cell.2004.12.019

CrossRef Full Text | Google Scholar

Schreiber, A. W., Hayden, M. J., Forrest, K. L., Kong, S. L., Langridge, P., and Baumann, U. (2012). Transcriptome-scale homoeolog-specific transcript assemblies of bread wheat. BMC Genomics 13:492. doi: 10.1186/1471-2164-13-492

CrossRef Full Text | Google Scholar

Shewry, P. R. (2009). Wheat. J. Exp. Bot. 60, 1537–1553. doi: 10.1093/jxb/erp058

CrossRef Full Text | Google Scholar

Slade, A. J., McGuire, C., Loeffler, D., Mullenberg, J., Skinner, W., Fazio, G., et al. (2012). Development of high amylose wheat through TILLING. BMC Plant Biol. 12:69. doi: 10.1186/1471-2229-12-69

CrossRef Full Text | Google Scholar

Song, X. J., Huang, W., Shi, M., Zhu, M. Z., and Lin, H. X. (2007). A QTL for rice grain width and weight encodes a previously unknown RING-type E3 ubiquitin ligase. Nat. Genet. 39, 623–630. doi: 10.1038/ng2014

CrossRef Full Text | Google Scholar

Tanaka, T., Kobayashi, F., Joshi, G. P., Onuki, R., Sakai, H., Kanamori, H., et al. (2014). Next-generation survey sequencing and the molecular organization of wheat chromosome 6B. DNA Res. 21, 103–114. doi: 10.1093/dnares/dst041.

CrossRef Full Text | Google Scholar

The International Wheat Genome Sequencing Consortium [IWGSC] (2014). A chromosome-based draft sequence of the hexaploid bread wheat (Triticum aestivum) genome. Science 345:1251788. doi: 10.1126/science.1251788

CrossRef Full Text | Google Scholar

Till, B. J., Cooper, J., Tai, T. H., Colowit, P., Greene, E. A., Henikoff, S., et al. (2007). Discovery of chemically induced mutations in rice by TILLING. BMC Plant Biol. 7:19. doi: 10.1186/1471-2229-7-19

CrossRef Full Text | Google Scholar

Tsai, H., Howell, T., Nitcher, R., Missirian, V., Watson, B., Ngo, K. J., et al. (2011). Discovery of rare mutations in populations: TILLING by sequencing. Plant Physiol. 156, 1257–1268. doi: 10.1104/pp.110.169748

CrossRef Full Text | Google Scholar

Uauy, C., Paraiso, F., Colasuonno, P., Tran, R. K., Tsai, H., Berardi, S., et al. (2009). A modified TILLING approach to detect induced mutations in tetraploid and hexaploid wheat. BMC Plant Biol. 9:115. doi: 10.1186/1471-2229-9-115

CrossRef Full Text | Google Scholar

Wang, E., Wang, J., Zhu, X., Hao, W., Wang, L., Li, Q., et al. (2008). Control of rice grain-filling and yield by a gene with a potential signature of domestication. Nat. Genet. 40, 1370–1374. doi: 10.1038/ng.220

CrossRef Full Text | Google Scholar

Wang, H. C., Singer, G. A., and Hickey, D. A. (2004). Mutational bias affects protein evolution in flowering plants. Mol. Biol. Evol. 21, 90–96. doi: 10.1093/molbev/msh003

CrossRef Full Text | Google Scholar

Wang, S., Wu, K., Yuan, Q., Liu, X., Liu, Z., Lin, X., et al. (2012). Control of grain size, shape and quality by OsSPL16 in rice. Nat. Genet. 44, 950–954. doi: 10.1038/ng.2327

CrossRef Full Text | Google Scholar

Wang, Y., Cheng, X., Shan, Q., Zhang, Y., Liu, J., Gao, C., et al. (2014). Simultaneous editing of three homoeoalleles in hexaploid bread wheat confers heritable resistance to powdery mildew. Nat. Biotechnol. 32, 947–951. doi: 10.1038/nbt.2969

CrossRef Full Text | Google Scholar

Wang, Z., Gerstein, M., and Snyder, M. (2009). RNA-Seq: a revolutionary tool for transcriptomics. Nat. Rev. Genet. 10, 57–63. doi: 10.1038/nrg2484

CrossRef Full Text | Google Scholar

Weng, J., Gu, S., Wan, X., Gao, H., Guo, T., Su, N., et al. (2008). Isolation and initial characterization of GW5, a major QTL associated with rice grain width and weight. Cell Res. 18, 1199–1209. doi: 10.1038/cr.2008.307

CrossRef Full Text | Google Scholar

Wu, Q. H., Chen, Y. X., Zhou, S. H., Fu, L., Chen, J. J., Xiao, Y., et al. (2015). High-density genetic linkage map construction and QTL mapping of grain shape and size in the wheat population Yanda1817 x Beinong6. PLoS ONE 10:e0118144. doi: 10.1371/journal.pone.0118144

CrossRef Full Text | Google Scholar

Zanke, C. D., Ling, J., Plieske, J., Kollers, S., Ebmeyer, E., Korzun, V., et al. (2015). Analysis of main effect QTL for thousand grain weight in European winter wheat (Triticum aestivum L.) by genome-wide association mapping. Front. Plant Sci. 6:644. doi: 10.3389/fpls.2015.00644

CrossRef Full Text | Google Scholar

Zhang, Y., Liu, J., Xia, X., and He, Z. (2014). TaGS-D1, an ortholog of rice OsGS3, is associated with grain weight and grain length in common wheat. Mol. Breed. 34, 1097–1107. doi: 10.1007/s11032-014-0102-7

CrossRef Full Text | Google Scholar

Keywords: wheat, ethyl methanesulfonate (EMS), mutation, next-generation sequencing, transcriptome, single nucleotide variations (SNVs)

Citation: Chen F, Zhu Z, Zhou X, Yan Y, Dong Z and Cui D (2016) High-Throughput Sequencing Reveals Single Nucleotide Variants in Longer-Kernel Bread Wheat. Front. Plant Sci. 7:1193. doi: 10.3389/fpls.2016.01193

Received: 24 February 2016; Accepted: 25 July 2016;
Published: 08 August 2016.

Edited by:

Chengdao Li, Murdoch University, Australia

Reviewed by:

Andrea Mazzucato, Tuscia University, Italy
Zhongfu Ni, China Agricultural University, China

Copyright © 2016 Chen, Zhu, Zhou, Yan, Dong and Cui. 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: Feng Chen, chf0088@163.com