Whole genome resequencing in tomato reveals variation associated with introgression and breeding events

Background

One of the goals of genomics is to identify the genetic loci responsible for variation in phenotypic traits. The completion of the tomato genome sequence and recent advances in DNA sequencing technology allow for in-depth characterization of genetic variation present in the tomato genome. Like many self-pollinated crops, cultivated tomato accessions show a low molecular but high phenotypic diversity. Here we describe the whole-genome resequencing of eight accessions (four cherry-type and four large fruited lines) chosen to represent a large range of intra-specific variability and the identification and annotation of novel polymorphisms.

Results

The eight genomes were sequenced using the GAII Illumina platform. Comparison of the sequences with the reference genome yielded more than 4 million single nucleotide polymorphisms (SNPs). This number varied from 80,000 to 1.5 million according to the accessions. Almost 128,000 InDels were detected. The distribution of SNPs and InDels across and within chromosomes was highly heterogeneous revealing introgressions from wild species and the mosaic structure of the genomes of the cherry tomato accessions. In-depth annotation of the polymorphisms identified more than 16,000 unique non-synonymous SNPs. In addition 1,686 putative copy-number variations (CNVs) were identified.

Conclusions

This study represents the first whole genome resequencing experiment in cultivated tomato. Substantial genetic differences exist between the sequenced tomato accessions and the reference sequence. The heterogeneous distribution of the polymorphisms may be related to introgressions that occurred during domestication or breeding. The annotated SNPs, InDels and CNVs identified in this resequencing study will serve as useful genetic tools, and as candidate polymorphisms in the search for phenotype-altering DNA variations.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-791) contains supplementary material, which is available to authorized users.     

The effects of artificial selection on sugar metabolism and transporter genes in grape

Sugar content is a key feature of grape quality. The sugar content of grapes has been significantly improved after nearly a thousand years of artificial selection. However, the mechanism underlying the changes in the grape sugar content during the process of artificial selection remains largely unknown although several genes involved in sugar metabolism and transportation in grape have been identified. In this study, the genomes of 13 wild Vitis species and 14 cultivated Vitis vinifera accessions were resequenced to 2–5 X depth using the Illumina Hiseq2000 platform. Genetic variation of 138 genes involved in sugar biosynthesis and transport was investigated, and 7,690 and 12,717 single nucleotide polymorphisms/insertions and deletions (SNPs/InDel) were identified within the cultivated V. vinifera and wild Vitis species, respectively. The percentages of SNPs/InDels were 0.93 and 1.54 % in cultivated and wild species, respectively, and the wild Vitis species had 1.65-fold more SNPs/InDels than the cultivated V. vinifera. Moreover, the distribution of SNPs/InDels in gene regions was also investigated. Eight genes (HT4, PPFTK4, PPFTK6, PMT3, SPS1, HT8, HT15, SUSy1) showed low level of allelic diversity in cultivated species, suggesting they might have undergone purifying selection during the domestication process of grapes. Our genome DNA resequencing data provided a valuable resource for analyzing the effects of artificial selection on trait-related pathways in grape. The result that eight genes showed lower level of DNA variation in cultivated species than in wild species will be very helpful in understanding sugar accumulation in grapes.     

Genomewide discovery of DNA polymorphisms in rice cultivars with contrasting drought and salinity stress response and their functional relevance

Next‐generation sequencing technologies provide opportunities to understand the genetic basis of phenotypic differences, such as abiotic stress response, even in the closely related cultivars via identification of large number of DNA polymorphisms. We performed whole‐genome resequencing of three rice cultivars with contrasting responses to drought and salinity stress (sensitive IR64, drought‐tolerant Nagina 22 and salinity‐tolerant Pokkali). More than 356 million 90‐bp paired‐end reads were generated, which provided about 85% coverage of the rice genome. Applying stringent parameters, we identified a total of 1 784 583 nonredundant single‐nucleotide polymorphisms (SNPs) and 154 275 InDels between reference (Nipponbare) and the three resequenced cultivars. We detected 401 683 and 662 509 SNPs between IR64 and Pokkali, and IR64 and N22 cultivars, respectively. The distribution of DNA polymorphisms was found to be uneven across and within the rice chromosomes. One‐fourth of the SNPs and InDels were detected in genic regions, and about 3.5% of the total SNPs resulted in nonsynonymous changes. Large‐effect SNPs and InDels, which affect the integrity of the encoded protein, were also identified. Further, we identified DNA polymorphisms present in the differentially expressed genes within the known quantitative trait loci. Among these, a total of 548 SNPs in 232 genes, located in the conserved functional domains, were identified. The data presented in this study provide functional markers and promising target genes for salinity and drought tolerance and present a valuable resource for high‐throughput genotyping and molecular breeding for abiotic stress traits in rice.     

Genome-wide DNA polymorphisms in elite indica rice inbreds discovered by whole-genome sequencing

Advances in next-generation sequencing technologies have aided discovery of millions of genome-wide DNA polymorphisms, single nucleotide polymorphisms (SNPs) and insertions-deletions (InDels), which are an invaluable resource for marker-assisted breeding. Whole-genome resequencing of six elite indica rice inbreds (three cytoplasmic male sterile and three restorer lines) resulted in the generation of 338?million 75-bp paired-end reads, which provided 85.4% coverage of the Nipponbare genome. A total of 2?819?086 nonredundant DNA polymorphisms including 2?495?052 SNPs, 160?478 insertions and 163?556 deletions were discovered between the inbreds and Nipponbare, providing an average of 6.8 SNPs/kb across the genome. Distribution of SNPs and InDels in the chromosome was nonrandom with SNP-rich and SNP-poor regions being evident across the genome. A contiguous 4.3-Mb region on chromosome 5 with extremely low SNP density was identified. Overall, 83?262 nonsynonymous SNPs spanning 16?379 genes and 3620 nonsynonymous InDels in 2625 genes have been discovered which provide valuable insights into the basis underlying performance of the inbreds and the hybrids between these inbred combinations. SNPs and InDels discovered from this diverse set of indica rice inbreds not only enrich SNP resources for molecular breeding but also enable the study of genome-wide variations on hybrid performance.     

Construction of EMS-Induced Peanut Mutant Libraries and Identification of Pod-Related Traits Mutant Lines

Peanut (Arachis hypogaea L.) is an oil and economic crop of vital importance, and peanut pod is the key organ influencing the yield and processing quality. Hence, the Pod-related traits (PRTs) are considered as important agronomic traits in peanut breeding. To broaden the variability of PRTs in current peanut germplasms, three elite peanut cultivars were used to construct Ethyl methane sulfonate (EMS)-induced mutant libraries in this study. The optimal EMS treatment conditions for the three peanut varieties were determined. It was found that the median lethal dose (LD50) of EMS treatment varied greatly among different genotypes. Finally, the EMS-induced peanut mutant libraries were constructed and a total of 124 mutant lines for PRTs were identified and evaluated. Furthermore, “M-8070”, one of the mutant lines for pod constriction, was re-sequenced via high-throughput sequencing technology. The genome-wide variations between “M-8070” and its wild parent “Fuhua 8” (FH 8) were detected. 2994 EMS-induced single nucleotide polymorphisms (SNPs) and 1188 insertion-deletions (InDels) between “M-8070” and its wild parent were identified. The predominant SNP mutation type was C/G to T/A transitions, while the predominant InDel mutation type was “1-bp”. We analyzed the distribution of identified mutations and annotated their functions. Most of the mutations (91.68% of the SNPs and 77.69% of the InDels) were located in the intergenic region. 72 SNPs were identified in the exonic region, leading to 27 synonymous, 43 non-synonymous and 2 stop-gain variation for gene structure. 13 Indels were identified in the exonic region, leading to 4 frame-shift, 8 non-frame-shift and 1 stop-gain variations of genes. These mutations may lead to the phenotypic variation of “M-8070”. Our study provided valuable resources for peanut improvement and functional genomic research.     

First-generation SNP/InDel markers tagging loci for pathogen resistance in the potato genome

A panel of 17 tetraploid and 11 diploid potato genotypes was screened by comparative sequence analysis of polymerase chain reaction (PCR) products for single nucleotide polymorphisms (SNPs) and insertion-deletion polymorphisms (InDels), in regions of the potato genome where genes for qualitative and/or quantitative resistance to different pathogens have been localized. Most SNP and InDel markers were derived from bacterial artificial chromosome (BAC) insertions that contain sequences similar to the family of plant genes for pathogen resistance having nucleotide-binding-site and leucine-rich-repeat domains (NBS-LRR-type genes). Forty-four such NBS-LRR-type genes containing BAC-insertions were mapped to 14 loci, which tag most known resistance quantitative trait loci (QTL) in potato. Resistance QTL not linked to known resistance-gene-like (RGL) sequences were tagged with other markers. In total, 78 genomic DNA fragments with an overall length of 31 kb were comparatively sequenced in the panel of 28 genotypes. 1498 SNPs and 127 InDels were identified, which corresponded, on average, to one SNP every 21 base pairs and one InDel every 243 base pairs. The nucleotide diversity of the tetraploid genotypes (pi = 0.72 x 10(-3)) was lower when compared with diploid genotypes (pi = 2.31 x 10(-3)). RGL sequences showed higher nucleotide diversity when compared with other sequences, suggesting evolution by divergent selection. Information on sequences, sequence similarities, SNPs and InDels is provided in a database that can be queried via the Internet.     

Genome-Wide Identification of SSR and SNP Markers Based on Whole-Genome Re-Sequencing of a Thailand Wild Sacred Lotus (Nelumbo nucifera)

Genomic resources such as single nucleotide polymorphism (SNPs), insertions and deletions (InDels) and SSRs (simple sequence repeats) are essential for crop improvement and better utilization in genetic breeding. However, the resources for the sacred lotus (Nelumbo nucifera Gaertn.) are still limited. In the present study, to dissect large-scale genomic molecular marker resources for sacred lotus, we re-sequenced a Thailand sacred lotus cultivar ‘Chiang Mai wild lotus’ and compared with the reported lotus genome ‘Middle lake wild lotus’. A total of 3,180,059 SNPs, 328, 251 InDels and 14,191 SVs were found between the two genomes. The functional impact analyses of these SNPs indicated that they may be involved in metabolic processes, binding, catalytic activity, etc. Mining the genome sequences for SSRs showed that 191,657 SSRs were identified with a frequency of one SSR per 4.23 kb and 103,656 SSR primer pairs were designed. Furthermore, 14, 502 EST-SSRs were also indentified using the available RNA-seq data in the NCBI. A subset of 150 SSRs (genomic and EST-SSRs) was randomly selected for validation and genetic diversity analysis. The genotypes could be easily distinguished using these SSR markers and the ‘Chiang Mai wild lotus’ was obviously differentiated from the other Chinese accessions. This study provides considerable amounts of genomic resources and markers for the quantitative trait locus (QTL) identification and molecular selection of the species, which could have a potential role in various applications in sacred lotus breeding.     

Genome-wide discovery of DNA polymorphism in Brassica rapa

Single nucleotide polymorphisms (SNPs) and/or insertion/deletions (InDels) are frequent sequence variations in the plant genome, which can be developed as molecular markers for genetic studies on crop improvement. The ongoing Brassica rapa genome sequencing project has generated vast amounts of sequence data useful in genetic research. Here, we report a genome-wide survey of DNA polymorphisms in the B. rapa genome based on the 557 bacterial artificial clone sequences of B. rapa ssp. pekinensis cv. Chiifu. We identified and characterized 21,311 SNPs and 6,753 InDels in the gene space of the B. rapa genome by re-sequencing 1,398 sequence-tagged sites (STSs) in eight genotypes. Comparison of our findings with a B. rapa genetic linkage map confirmed that STS loci were distributed randomly over the B. rapa whole genome. In the 1.4 Mb of aligned sequences, mean nucleotide polymorphism and diversity were θ = 0.00890 and π = 0.00917, respectively. Additionally, the nucleotide diversity in introns was almost three times greater than that in exons, and the frequency of observed InDel was almost 17 times higher in introns than in exons. Information regarding SNPs/InDels obtained here will provide an important resource for genetic studies and breeding programs of B. rapa.     

A catalog of neutral and deleterious polymorphism in yeast

The abundance and identity of functional variation segregating in natural populations is paramount to dissecting the molecular basis of quantitative traits as well as human genetic diseases. Genome sequencing of multiple organisms of the same species provides an efficient means of cataloging rearrangements, insertion, or deletion polymorphisms (InDels) and single-nucleotide polymorphisms (SNPs). While inbreeding depression and heterosis imply that a substantial amount of polymorphism is deleterious, distinguishing deleterious from neutral polymorphism remains a significant challenge. To identify deleterious and neutral DNA sequence variation within Saccharomyces cerevisiae, we sequenced the genome of a vineyard and oak tree strain and compared them to a reference genome. Among these three strains, 6% of the genome is variable, mostly attributable to variation in genome content that results from large InDels. Out of the 88,000 polymorphisms identified, 93% are SNPs and a small but significant fraction can be attributed to recent interspecific introgression and ectopic gene conversion. In comparison to the reference genome, there is substantial evidence for functional variation in gene content and structure that results from large InDels, frame-shifts, and polymorphic start and stop codons. Comparison of polymorphism to divergence reveals scant evidence for positive selection but an abundance of evidence for deleterious SNPs. We estimate that 12% of coding and 7% of noncoding SNPs are deleterious. Based on divergence among 11 yeast species, we identified 1,666 nonsynonymous SNPs that disrupt conserved amino acids and 1,863 noncoding SNPs that disrupt conserved noncoding motifs. The deleterious coding SNPs include those known to affect quantitative traits, and a subset of the deleterious noncoding SNPs occurs in the promoters of genes that show allele-specific expression, implying that some cis-regulatory SNPs are deleterious. Our results show that the genome sequences of both closely and distantly related species provide a means of identifying deleterious polymorphisms that disrupt functionally conserved coding and noncoding sequences.     

Genome-wide DNA polymorphism and transcriptome analysis of an early-maturing rice mutant

In order to develop a rice population with improved important traits such as flowering time, we developed 2,911 M₂ targeting-induced local lesions in genomes (TILLING) lines by irradiating rice seeds with γ-rays. In all, 15 M₃ lines were obtained from 3 different M₂ lines that exhibited an early-maturing phenotype: these plants matured approximately 25 days faster than wild-type (WT) plants. To identify genome-wide DNA polymorphisms, we performed whole-genome resequencing of both the plant types, i.e., WT and early-maturing TILLING 1 (EMT1), and obtained mapped reads of 118,488,245 bp (99.53 %) and 128,489,860 bp (99.72 %), respectively; Nipponbare was used as the reference genome. We obtained 63,648 and 147,728 single nucleotide polymorphisms (SNPs) and 33,474 and 31,082 insertions and deletions (InDels) for the WT and EMT1, respectively. Interestingly, there was a higher number of SNPs (2.6-fold) and slightly lower number of InDels (0.9-fold) in EMT1 than in WT. The expression of at least 202 structurally altered genes was changed in EMT1, and functional enrichment analysis of these genes revealed that their molecular functions were related to flower development. These results might provide a critical insight into the regulatory pathways of rice flowering.     

Genome-wide searching of single-nucleotide polymorphisms among eight distantly and closely related rice cultivars (Oryza sativa L.) and a wild accession (Oryza rufipogon Griff.).

We searched the genomes of eight rice cultivars (Oryza sativa L. ssp. japonica and ssp. indica) and a wild rice accession (Oryza rufipogon Griffith) for nucleotide polymorphisms, and identified 7805 polymorphic loci, including single-nucleotide polymorphisms (SNPs) and insertions/deletions (InDels), in predicted intergenic regions. Polymorphisms are useful as DNA markers for genetic analysis or positional cloning with segregating populations of crosses. Pairwise comparison between cultivars and a neighbor-joining tree calculated from SNPs agreed very well with relationships between rice strains predicted from pedigree data or calculated with other DNA markers such as p-SINE1 and simple sequence repeats (SSRs), suggesting that whole-genome SNP information can be used for analysis of evolutionary relationships. Using multiple SNPs to identify alleles, we drew a map to illustrate the alleles shared among the eight cultivars and the accession. The map revealed that most of the genome is mono- or di-allelic among japonica cultivars, whereas alleles well conserved among modern japonica paddy rice cultivars were often shared with indica cultivars or wild rice, suggesting that the genome structure of modern cultivars is composed of chromosomal segments from various genetic backgrounds. Use of allele-sharing analysis and association analysis were also tested and are discussed.     

Identification of genome-wide variations among three elite restorer lines for hybrid-rice

Rice restorer lines play an important role in three-line hybrid rice production. Previous research based on molecular tagging has suggested that the restorer lines used widely today have narrow genetic backgrounds. However, patterns of genetic variation at a genome-wide scale in these restorer lines remain largely unknown. The present study performed re-sequencing and genome-wide variation analysis of three important representative restorer lines, namely, IR24, MH63, and SH527, using the Solexa sequencing technology. With the genomic sequence of the Indica cultivar 9311 as the reference, the following genetic features were identified: 267,383 single-nucleotide polymorphisms (SNPs), 52,847 insertion/deletion polymorphisms (InDels), and 3,286 structural variations (SVs) in the genome of IR24; 288,764 SNPs, 59,658 InDels, and 3,226 SVs in MH63; and 259,862 SNPs, 55,500 InDels, and 3,127 SVs in SH527. Variations between samples were also determined by comparative analysis of authentic collections of SNPs, InDels, and SVs, and were functionally annotated. Furthermore, variations in several important genes were also surveyed by alignment analysis in these lines. Our results suggest that genetic variations among these lines, although far lower than those reported in the landrace population, are greater than expected, indicating a complicated genetic basis for the phenotypic diversity of the restorer lines. Identification of genome-wide variation and pattern analysis among the restorer lines will facilitate future genetic studies and the molecular improvement of hybrid rice.     

Development of PCR-based SNP markers for rice blast resistance genes at the<Emphasis Type="Italic"> Piz</Emphasis> locus

We assessed the utility of single-nucleotide polymorphisms (SNPs) and small insertion/deletion polymorphisms (InDels) as DNA markers in genetic analysis and breeding of rice. Toward this end, we surveyed SNPs and InDels in the chromosomal region containing the Piz and Piz-t rice blast resistance genes and developed PCR-based markers for typing the SNPs. Analysis of sequences from a blast-susceptible Japanese cultivar and two cultivars each containing one of these genes revealed that SNPs are abundant in the Piz and Piz-t regions (on average, one SNP every 248 bp), but the number of InDels was much lower. The dense distribution of SNPs facilitated the generation of SNP markers in the vicinity of the genes. For typing these SNPs, we used a modified allele-specific PCR method. Of the 49 candidate allele-specific markers, 33 unambiguously and reproducibly discriminated between the two alleles. We used the markers for mapping the Piz and Piz-t genes and evaluating the size of DNA segments introgressed from the Piz donor cultivar in Japanese near-isogenic lines containing Piz. Our findings suggest that, because of its ability to generate numerous markers within a target region and its simplicity in assaying genotypes, SNP genotyping with allele-specific PCR is a valuable tool for gene mapping, map-based cloning, and marker-assisted selection in crops, especially rice.Communicated by D.J. Mackill     

Identification and characterization of nucleotide variations in the genome of Ziziphus jujuba (Rhamnaceae) by next generation sequencing

In this study, single-nucleotide polymorphisms (SNPs) and insertions/deletions (InDels) in the genome of Ziziphus jujuba were identified using sequences generated by the Roche 454 GS-FLX sequencer. A total of, 573,141 reads were produced with an average read length of 360 bp. After quality control, 258,754 of the filtered reads were assembled into 23,864 contigs, and 293,458 remained as singletons. Using the contig assemblies as a reference, 17,160 SNPs and 478 InDels were identified. Among the SNPs, transitions occurred three times more frequently than transversions. In transitions, the number of C/T and G/A transitions was similar. Among the transversions, A/T was the most abundant, and C/G was much rarer than any of the other types of transversions, accounting for only about half the numbers of A/C, A/T and G/T transversions. For the InDels, mononucleotide changes amounted to 64.4 % of the total number of InDels. In general, the frequency of detected InDels decreased as the length of the InDels increased. This study provides valuable marker resources for future genetic studies of Ziziphus spp.     

Whole genome resequencing of tobacco (<Emphasis Type="Italic">Nicotiana tabacum</Emphasis> L.) genotypes and high-throughput SNP discovery

Genome-wide single-nucleotide polymorphisms (SNPs) are highly useful in unraveling genetic insights and are essential to accelerate selections for genetic improvement in tobacco. The discovery of genome-wide SNPs in tobacco is very complex due to its high level of repetitive genome and polyploidy. At present, publicly available genomic data on SNPs are very limited, which warrants the need for high-throughput SNPs for application in tobacco breeding. In this research paper, we describe our efforts on SNP discovery by whole genome resequencing of 18 flue-cured Virginia (FCV) tobacco genotypes and annotation of SNPs in the tobacco genome. A large amount of data of about 225 GB per genotype was generated, with an average read depth of 50× using paired-end next-generation sequencing (NGS) with the HiSeq 2500 platform. The discovery of a large number of SNPs and indels was attempted to assist mapping and, thus, the selection processes to develop superior tobacco breeding lines. Discovered SNPs, their functional annotation, mapping to the reference genome, and their relative positioning in the linkage group are discussed in this paper.     

Comparison of whole‐genome (13X) and capture (87X) resequencing methods for SNP and genotype callings

The number of polymorphisms identified with next‐generation sequencing approaches depends directly on the sequencing depth and therefore on the experimental cost. Although higher levels of depth ensure more sensitive and more specific SNP calls, economic constraints limit the increase of depth for whole‐genome resequencing (WGS). For this reason, capture resequencing is used for studies focusing on only some specific regions of the genome. However, several biases in capture resequencing are known to have a negative impact on the sensitivity of SNP detection. Within this framework, the aim of this study was to compare the accuracy of WGS and capture resequencing on SNP detection and genotype calling, which differ in terms of both sequencing depth and biases. Indeed, we have evaluated the SNP calling and genotyping accuracy in a WGS dataset (13X) and in a capture resequencing dataset (87X) performed on 11 individuals. The percentage of SNPs not identified due to a sevenfold sequencing depth decrease was estimated at 7.8% using a down‐sampling procedure on the capture sequencing dataset. A comparison of the 87X capture sequencing dataset with the WGS dataset revealed that capture‐related biases were leading with the loss of 5.2% of SNPs detected with WGS. Nevertheless, when considering the SNPs detected by both approaches, capture sequencing appears to achieve far better SNP genotyping, with about 4.4% of the WGS genotypes that can be considered as erroneous and even 10% focusing on heterozygous genotypes. In conclusion, WGS and capture deep sequencing can be considered equivalent strategies for SNP detection, as the rate of SNPs not identified because of a low sequencing depth in the former is quite similar to SNPs missed because of method biases of the latter. On the other hand, capture deep sequencing clearly appears more adapted for studies requiring great accuracy in genotyping.     

Genome-wide single nucleotide polymorphism and Insertion-Deletion discovery through next-generation sequencing of reduced representation libraries in common bean

Single nucleotide polymorphisms (SNPs) and insertions-deletions (InDels) are valuable molecular markers for genomics and genetics studies and molecular breeding. The advent of next-generation sequencing techniques has enabled researchers to approach high-throughput and cost-effective SNP and InDel discovery on a genomic scale. In this report, 36 common bean genotypes grown in Canada were used to construct reduced representation libraries for next-generation sequencing. Using 76 million sequence reads generated by the Illumina HiSeq 2000 Sequencing System, we identified a total of 43,698 putative SNPs and 1,267 putative InDels. Of the SNPs, 43,504 were bi-allelic and 194 were tri-allelic, and the InDels comprised 574 insertions and 693 deletions. The putative bi-allelic SNPs were distributed across all 11 chromosomes with the highest number of SNPs observed in chromosome 2 (4,788), and the lowest in chromosome 10 (2,941). With the aid of the recent release of the first chromosome-scale version of Phaseolus vulgaris, 24,907 bi-allelic SNPs, 79 tri-allelic SNPs, 315 insertions, and 377 deletions were located in 8,758, 77, 273, and 364 genes, respectively. Among these 24,907 bi-allelic SNPs, 7,168 nonsynonymous bi-allelic SNPs were identified within 36 common bean genotypes that were located in 4,303 genes. A total of 113 putative SNPs were randomly chosen for validation using high-resolution melt analysis. Of the 113 candidate SNPs, 105 (92.9 %) contained the predicted SNPs.