Skip to main content
Download PDF
  • Research article
  • Open access
  • Published:

Genome-wide copy number variations in Oryza sativa L.

BMC Genomics volume 14, Article number: 649 (2013) Cite this article

Abstract

Background

Copy number variation (CNV) can lead to intra-specific genome variations. It is not only part of normal genetic variation, but also is the source of phenotypic differences. Rice (Oryza sativa L.) is a model organism with a well-annotated genome, but investigation of CNVs in rice lags behind its mammalian counterparts.

Results

We comprehensively assayed CNVs using high-density array comparative genomic hybridization in a panel of 20 Asian cultivated rice comprising six indica, three aus, two rayada, two aromatic, three tropical japonica, and four temperate japonica varieties. We used a stringent criterion to identify a total of 2886 high-confidence copy number variable regions (CNVRs), which span 10.28 Mb (or 2.69%) of the rice genome, overlapping 1321 genes. These genes were significantly enriched for specific biological functions involved in cell death, protein phosphorylation, and defense response. Transposable elements (TEs) and other repetitive sequences were identified in the majority of CNVRs. Chromosome 11 showed the greatest enrichment for CNVs. Of subspecies-specific CNVRs, 55.75% and 61.96% were observed in only one cultivar of ssp. indica and ssp. japonica, respectively. Some CNVs with high frequency differences among groups resided in genes underlying rice adaptation.

Conclusions

Higher recombination rates and the presence of homologous gene clusters are probably predispositions for generation of the higher number of CNVs on chromosome 11 by non-allelic homologous recombination events. The subspecies-specific variants are enriched for rare alleles, which suggests that CNVs are relatively recent events that have arisen within breeding populations. A number of the CNVs identified in this study are candidates for generation of group-specific phenotypes.

Background

In recent years, rice genomics has progressed substantially and generated considerable valuable resources, including availability of two independent genome sequences [1, 2] and a composite physical map [3]. These resources provide a foundation for understanding the tremendous genetic diversity that exists in rice. Among organisms for which a high-quality genome sequence from at least one individual is available, such as human, mouse, and Arabidopsis, genome-wide surveys of single nucleotide polymorphisms (SNPs) have captured significant proportions of within-species variation [4–6]. Similarly, the rice research community has migrated to SNPs as the main measure of genetic variation in rice, with initiation of the OryzaSNP project (http://www.OryzaSNP.org), and subsequently ~160,000 nonredundant SNPs distributed across the entire genome of the OryzaSNPset have been discovered [7]. The abundant SNPs can be used to design more targeted SNP assays for immediate use in rice genetics and molecular breeding [8–10], particularly in Genome-wide association study (GWAS) [11–13].

However, genomics research has revealed other forms of genetic variations, such as copy number variation (CNV) in human [14, 15], chimpanzee [16], dog [17, 18], chicken [19], cattle [20], pig [21], rat [22], mouse [23], Drosophila[24], Caenorhabditis elegans[25], yeast [26], Escherichia coli[27], maize [28–30], and Arabidopsis thaliana[31, 32]. CNV is defined as a segment of DNA ≥1 kb that is variable in copy number in comparison with a reference genome [33]. CNV covers more base-pairs [34–38] and has a higher per-locus mutation rate than SNPs [39]. Previous studies indicated that CNVs are in variable linkage disequilibrium (LD) with flanking SNPs [40–42].

High levels of CNV have been found throughout the rice genome [43, 44], and recent study provided insight into the extent of genome-wide structural variations in the important representative restorer lines [45]. Herein we describe a more detailed survey undertaken to detect candidate CNVs in a panel of 20 Asian cultivated rice comprising six indica, three aus, two rayada, two aromatic, three tropical japonica, and four temperate japonica (Additional file 1: Table S1), and examine the genome-wide characteristics of CNVs in subspecies and groups. These resources allowed us to analyze genetic diversity as indicated by CNVs, to evaluate the biological roles of CNVs, and to identify candidate CNVs that are likely to occur independently in subspecies and contribute to the genetic differences among groups.

Results

Detection of CNVs in Asian cultivated rice

We performed array-comparative genomic hybridization (aCGH) covering the entire rice reference genome sequence of ssp. japonica Nipponbare (IRGSP v4.0). DNA samples from Nipponbare, as the reference sample, were fixed in all hybridization experiments. Genomic DNA from 20 accessions of Asian cultivated rice was tested against the Nipponbare reference sample. We used an updated version of a previously described method to identify changes in log2 signal intensity corresponding to copy number gains and losses [43]. Under our conservative calling criteria, self-to-self hybridization showed no detectable false-positives. Figure 1 compares the linearized whole-genome karyogram of IR 64 to that of a self–self (Nipponbare), and shows the whole chromosome 8 in greater detail. Using a set of stringent criteria, a total of 12,224 CNVs in the 20 cultivars were identified, with an average of 611.2 per cultivar, ranging from 235 CNVs in Geumobyeo to 876 CNVs in DV 85 (Additional file 2: Table S2; Figure 2). A higher number of CNV events were detected in ssp. indica (746.8 per cultivar) than in ssp. japonica (500.3 per cultivar) (Table 1). However, among ssp. japonica cultivars, more CNV events were detected in tropical japonica (520.3 per cultivar) than in temperate japonica (320.0 per cultivar) (Additional file 2: Table S2). Analysis of variance (ANOVA) indicated that CNV length did not vary significantly among the cultivars (F = 0.551, P = 0.941), with 98.5% of the CNV lengths within 10 kb.

Figure 1
figure 1

Whole-genome DNA copy number quantification of the IR 64 by aCGH. The linearized whole-genome copy number in shown in (a) Nipponbare (self–self hybridization), (b) IR 64, and (c) on chromosome 8 of IR 64.

Full size image
Figure 2
figure 2

Number of CNVs detected in each of 20 rice cultivars. The numbers of cultivars carrying the CNVs are indicated by different colors.

Full size image
Table 1 Summary of CNVs identified in two rice subspecies
Full size table

Identification and distribution of copy number variable regions

By merging overlapping CNVs identified in all cultivars across the aCGH [15], 2886 high-confidence copy number variable regions (CNVRs) were identified, which covered 10.28 Mb of the rice genome (Additional file 3: Table S3). The CNVRs consisted of 2557 losses and 276 gains in copy number, and 53 with both events. Furthermore, 992 CNVRs (34.4%) were detected in only one cultivar (unique), whereas the remaining CNVRs (65.6%) were detected in two or more cultivars (Additional file 3: Table S3). However, more CNVRs were predicted to be present in ssp. japonica than in ssp. indica (Additional file 4: Table S4), which was inconsistent with the analysis of CNVs. The discrepancy may be attributable to the higher density of CNVs and fewer ssp. indica cultivars sampled.

The CNVRs were distributed throughout all 12 rice chromosomes (Figure 3; Additional file 5: Table S5). The length of CNVRs differed significantly among different chromosomes, and the percentage of entire chromosomes susceptible to CNVRs ranged from 1.07% on chromosome 3 to 5.78% on chromosome 11. Although chromosome 11 spanned only 7.55% of the probes on the whole genome microarray, it included 17.54% of the CNVRs identified and thus showed the greatest enrichment for CNVs with ~2.15-fold the average variable content of the genome. There was little correlation between CNV occurrence and chromosome length (Additional file 5: Table S5), which is consistent with previous studies on the heterogeneous distribution of CNVs throughout the genome [15, 46, 47].

Figure 3
figure 3

Genome-wide distribution of CNVRs. The ideograms depict the chromosomal location of copy number gains (green), losses (red), and gains or losses (blue) relative to Nipponbare in the 12 chromosomes of 20 rice cultivars. The approximate position of each centromere is indicated by a black line on each chromosome.

Full size image

Overlap of CNVs with segmental duplications (SDs) and other repeat contents

When we intersected the 2886 CNVRs with rice SD regions, 10.91% (315/2886) of the CNVRs directly overlapped with SDs. A reason is that only a unique sequence is spotted in the oligonucleotide array. This selectively omits or reduces the probe density in SDs. The regions identified as deletions relative to the reference sample were analyzed for repeat content (including transposable elements; TEs) for each chromosome (Table 2). The most common classes of TEs (long interspersed elements [LINEs], short interspersed elements [SINEs], long terminal repeat [LTR] retrotransposons, and DNA TEs) and other repetitive sequences, such as simple repeat sequences and other unclassified repeats, were identified in the majority of regions (2175/2557) (Additional file 6: Table S6). Of these regions, 639 contained LTR retrotransposons and 1080 contained no DNA TEs. Although differences in the number within each TE class at specific locations varied widely, LTR retrotransposons were the predominant class (45.1% of total TEs) followed by DNA TEs (almost 39.8% of total TEs). In addition, SINEs and LINEs were identified on all chromosomes and only in 6.1% (157) and 9.4% (240), respectively, of all deleted regions. The largest sequence length within the LINE signatures was 5540 bp. Chromosome 2 contained only five LINEs.

Table 2 Transposable elements (TE)/repeat contents of the deleted regions on each rice chromosome
Full size table

Differentiation of subspecies by CNVs

Hierarchical clustering was performed to explore relationships among the 20 Asian cultivated rice based on the CNVs (Figure 4). Unsupervised clustering of the aCGH log2 ratio of all probes in the 2886 CNVRs supported known rice relationships and partially clustered samples consistent with previous studies [48], which indicated that those CNVs shared among samples may be indicators of a common evolutionary history and genetic relationship. However, the two rice subspecies, indica and japonica, did not form separate clusters. This result might be an artifact of the measurements being based on a common reference sample (Nipponbare), or the CNVs may be less phylogenetically reliable than neutral markers by reason of stronger selective constraints [49].

Figure 4
figure 4

Hierarchical clustering of the aCGH log 2 ratio of all probes in 2886 CNVRs. Cultivars are color-coded by group: indica in red, aus in blue, rayada in green, aromatic in orange, tropical japonica in purple, and temperate japonica in black.

Full size image

The approximate number of alleles was estimated using PowerMarker v3.25 software. Nei’s genetic diversity index (He) value for ssp. indica and ssp. japonica at every polymorphic CNVR was 0.329 and 0.281, respectively. This was expected on the basis of previous findings that ssp. indica contains higher genetic diversity than ssp. japonica[50–53]. To better understand how CNV contributes to intra-subspecies diversity, we searched for CNVRs that exhibited high levels of heterozygosity and identified 14 CNVRs with higher diversity in one or both subspecies (He > 0.500). For example, one CNVR (chr10: 19103630–19117640) was identified to be polymorphic in both subspecies. We queried the gene content of this region; interestingly, it contained no known gene. The distribution of Fst across all CNVRs ranged from 0 to 1.000 (Figure 5), with an average Fst value of 0.138, which was slightly lower than estimates reported in previous studies in which global germplasm collections were used in combination with SNP or SSR data [11, 12, 48, 54, 55].

Figure 5
figure 5

Frequency distribution of F st for all CNVRs.

Full size image

By further assessing the frequency of events separately observed in ssp. indica and ssp. japonica, 57.8% (1668/2886) of the variants were observed only in ssp. indica or ssp. japonica, but many of these were unique to one cultivar of ssp. indica or ssp. japonica (373/669 and 619/999, respectively). Only one CNVR (chr1: 35391036–35394841) was observed in all ssp. indica cultivars and no CNVR was observed in all ssp. japonica cultivars (Additional file 7: Table S7). This CNVR also overlapped in a previous aCGH analysis of Guangluai 4 and Nipponbare [43] and contained a subtilase gene Os01g0794800, a member of the superfamily of subtilisin-like serine proteases [56].

To visualize the distribution of both deletions and amplifications within the six groups of cultivars, event frequencies were analyzed by hierarchical clustering (Figure 6). The clustering identified variants that were restricted to certain groups and those that were present in multiple groups. Of 81 CNVRs with high frequency differences among groups (Additional file 8: Table S8), 54.32% corresponded to 51 annotated genes or gene families, of which some are important in rice adaptation, including a NB-ARC domain-containing protein. On the basis of differences in CNV frequency among the groups, we hypothesize that some CNVs arose independently within different groups and contribute to group differences, and therefore are associated with group formation and adaptation.

Figure 6
figure 6

Hierarchical clustering of both CNVRs and groups based on CNVR frequency within groups of rice cultivars. The color indicates the type and frequency of each CNVR: red indicates amplification (a), yellow indicates deletion (b).

Full size image

Gene content of CNVs

A total of 412 genes were identified as entirely affected by CNVs, which included 63 genes that were completely duplicated and 322 that were completely deleted. Five gene-loss events among the 20 Asian cultivated rice were chosen at random for PCR validation, namely Os01g0358700 (similar to the powdery mildew resistance protein PM3b), Os10g0507000 (NB-ARC domain-containing protein), Os11g0668650 (calmodulin-binding protein-like family protein), Os07g0218200 (similar to terpene synthase 7), and Os12g0247700 (similar to jasmonate-induced protein) (Additional file 9: Table S9; Figure 7). Given that the breakpoint definition of CNVRs can be equivocal, some of the gene sequences might partially overlap rather than be encompassed by CNVRs. Nevertheless, 1321 annotated genes wholly or partially spanned 41.6% (1200/2886) of the 2886 CNVRs (Additional file 10: Table S10), of which 46 (3.5%) were non-protein-coding transcripts. Overall we found that 70.0% (924/1321) of the genes identified were gene family members in the rice reference genome.

Figure 7
figure 7

Molecular confirmation of gene loss events detected by aCGH. (a) The gene Os07g0218200 was encompassed entirely by a CNV (chr7: 6548063–6555290). (b) The two internal primers generated a 103 bp fragment; no product was detected by aCGH in the cultivars that contained the putative deletion. Individual lanes represent (left to right): DNA marker, Patnai 23, Tetep, Nanjing 11, Zhongshanbai, IR 64, Kalamkati, DV 85, Chuan 7, Matia Amon, Dal Katra, RTN 904, Xiangmigu, Dom Sofid, Maendalaqili, Kinastano, Honduras, Mansaku, Geumobyeo, Guihuahuang, Xianghe, and Nipponbare.

Full size image

Gene ontology (GO) analysis of the 1321 overlapping annotated genes indicated that the functional classes of genes that were enriched within the 2886 CNVRs were cell death, protein phosphorylation, and defense response genes. We also identified orthologous genes of maize and Arabidopsis based on homolog clustering (Additional file 11: Table S11). The 523 rice-specific genes affected by CNV included 347 for which no additional homologs were found within rice and 176 that were in multigene families. The remaining 798 genes were conserved in maize and/or Arabidopsis.

Discussion

Characterization of CNVs in Asian cultivated rice

Several factors affected estimation of the number of CNVs detected between different rice cultivars against Nipponbare. First, the array designs target only the sequences present in the Nipponbare genome, and therefore we could not detect the sequences present in other rice cultivars but absent from Nipponbare. Second, a portion of the significant differences in hybridization signals may be due to the presence of multiple polymorphisms within the probe sequence between the reference and test sample genome. To reduce errors in CNV detection because of SNPs, we used a secondary manual method to examine all aberrant segments and to remove presumptive false-positives caused by single outlier probes and ensured all probes of each aberrant segment met a log2 ratio cut-off. As a consequence of this approach, our method for detection of CNVs was biased toward detection of small CNVs. Finally, self–self hybridization showed appreciable variability, but a set of conservative calling criteria for the final set of high-confidence CNVs showed no false-positives for self–self control hybridizations. Nonetheless, it is desirable to map the locations of the breakpoints with as high a degree of accuracy as possible.

We compared the 2886 CNVRs detected with 1676 candidate CNVs reported previously using re–sequencing results from 50 accessions of cultivated and wild rice [44]. The majority of the variants identified in these two studies do not overlap. In addition to the differences in detection technology, we suspect that the main reasons for this discordance are (1) differences in the accessions analyzed in the two studies, and (2) genome coverage biases—only unique sequences were spotted in our oligonucleotide array. This approach selectively omits or reduces probe density in CNV-enriched regions, such as SDs and diverse repetitive sequences [57]. In addition, it may also indicate that the CNVs detected using different technological and analytical platforms show large variation in CNV resolution, which may affect CNV sizes.

Previous studies have reported a strong association between CNVs and SDs in human [14, 35, 57–59], chimpanzee [16], mouse [23, 60–62], and cattle [20]. In contrast, when we compared the 2886 CNVRs derived from aCGH data with rice SD regions, only 10.91% (315/2886) of the events overlapped with rice SDs. We suspect that the probe specificity in our arrays might severely interfere with CNV discovery. In the future, CGH arrays with unbiased genome coverage combined with improved CNV calling algorithms could remedy this imperfection.

A notable feature of the rice genome is that LTR retrotransposons were identified in 30.7% of the rice genome [63], which suggests that LTR retrotransposons may be an important source of CNVs in rice. Indeed, we performed a cursory analysis for TE/repeat contents in regions that were identified as being deletions relative to the reference sample. We identified 639 regions that contained LTR retrotransposons, which comprised over 45.1% of the total TEs identified in the deleted regions. Although these results are preliminary and require further analysis, they suggest that TEs and repeats may have played a role in shaping the genomic architecture of CNVs in rice.

Evidence for ancient whole-genome duplications in rice is reported [64, 65]. The most recent segmental duplication event in rice occurred between chromosome 11 and 12 after differentiation of grasses [64]. This fact may explain why CNVR length differed significantly among different chromosomes, with larger CNVRs detected on chromosomes 11 and 12 in the present study. More surprisingly, previous work shows that chromosome 11 is enriched with disease resistance genes [66], which may be attributable to the higher number of CNVs on chromosome 11. Higher recombination rates and the presence of homologous gene clusters are probably predispositions for non-allelic homologous recombination (NAHR) events, which tend to generate CNVs [15, 67].

Differentiation of CNVs between ssp. indica and ssp. japonica

Rice is traditionally classified into two subspecies, ssp. indica and ssp. japonica, on the basis of morphological characteristics, ecological adaptation, crossing ability and geographic origin [68]. In addition, evidence suggests that ssp. indica and ssp. japonica are the products of separate domestication events from the ancestral species, O. rufipogon[69, 70]. We investigated whether any CNVs are associated with ssp. indica– ssp. japonica differentiation. In our study, 57.8% (1668/2886) of variants were observed only in ssp. indica or ssp. japonica. However, many of these ssp. indica- and ssp. japonica-specific events were observed in only one cultivar (373/669 and 619/999, respectively). Therefore, the subspecies-specific variants are enriched for rare alleles and may represent relatively new events that have arisen within breeding populations. This is also consistent with the low proportion of long CNVs identified in our study. We hypothesize that these CNVs are associated with phenotypic diversity among rice cultivars, and further research is important to assess how these variants affect phenotype. We also searched for CNVs potentially associated with ssp. indica– ssp. japonica differentiation on the basis of high frequencies within ssp. indica and ssp. japonica. No variants were present in all ssp. japonica cultivars and not in any ssp. indica cultivar. Only one variant was observed in all ssp. indica cultivars and not in any ssp. japonica cultivar. Thus there was no evidence for strong effects of subspecies differentiation on structural variation. It should be noted that CNVs were documented based on comparison to a reference ssp. japonica genome, and therefore sequences present in the ssp. indica genome but not in the ssp. japonica genome were not detected. Thus, we speculate that subspecies differentiation-associated CNVs would be expected to be present in most ssp. indica cultivars, but in few or no ssp. japonica cultivars.

Insights into Asian cultivated rice relationships based on high-density CNV data

Generally, the population structure of domesticated species is influenced by the evolutionary history of the predomesticated ancestors, as well as by the complexity of breeding practices [55]. Strong selective pressure during the process of rice domestication has led to the formation of population substructure [48, 55, 71]. In some studies, the two subspecies have been further divided into five major groups (indica, aus, tropical japonica, temperate japonica, and aromatic) on the basis of SNP and indel data [55, 72]. In the present study, hierarchical clustering of CNVs indicated that the variation segregating with major rice groups was present in the samples of a rice group. However, smaller groups are adapted to specific ecosystems, which may be recognized as upland, deep water, or floating cultivars [71, 73]. Evidence suggests that rayada is less amenable to ex situ conservation because of their adaptation to deep water. This may be one reason why the two rayada cultivars did not cluster together in the present study.

We conservatively queried CNVRs that have high-confidence frequency differences among groups. The allele frequencies that we observed for structural variants suggested that some variants have been removed entirely from certain groups. The specific regions may be related to diversification within specific rice groups. Further screening is necessary to confirm that these CNVs are truly group-specific. This finding also implies that those loss events detected exclusively in one group will limit the potential for genetic improvement through selection within that group only.

Implications for biological roles for CNVs

It is generally assumed that individuals of the same species have very similar genome contents. Rice is autogamous and has relatively low nucleotide polymorphism rates compared to other crop species. However, the aCGH indicates that considerable genomic variation exists in Asian cultivated rice. Moreover, this variation may include substantial differences in gene content and gene structure, as was observed in the gene deletion analysis of the PCR data in the present study. In addition to 412 genes entirely encompassed by a CNV, some additional genes partially overlapped. This result could be explained by biological and technical factors. In addition to changes in gene dosage, the main mechanisms responsible for the potential effects of CNVs include reshaping of the gene structure and modification of the elements that regulate gene expression [74, 75]. In addition, probe space and stringent criteria for CNV calling may limit coverage of the full length of genes. Notably, the array design is based on the older reference genome (IRGSP v4.0), so additional variant genes would not be identifiable if they reside in regions that were additions to the most recent version of the genome sequence. Collectively, these data suggest that the true number of genes may be substantially larger than the 1321 genes we identified in this study.

High levels of structural variation in plant genomes is related to important quantitative variation [30]. For example, fine mapping, complementation testing, and association analysis of a recently identified quantitative trait locus, qSW5, indicated that a deletion in qSW5 resulted in a significant increase in sink size owing to an increase in cell number in the outer glume of the rice flower [76]. As was demonstrated in maize, CNVs are also important in plant disease responses by directly affecting causative genes [29]. Similarly, it is important to evaluate the effect of such structural variation on phenotypic plasticity in rice cultivars, groups and subspecies. Although there was no obvious phenotypic implication in the present study, a marked enrichment of genes for cell death, protein phosphorylation, and defense response may reflect the importance of CNVs in biological processes. We identified many disease resistance genes within the CNVRs. This is expected because, first, gene families with functions in regulatory processes and signal recognition, such as disease resistance, have higher nonsynonymous-to-synonymous substitution ratios [7, 44]. This result, as for the over-represented set of disease resistance genes, might indicate a relaxation of constraints because of the redundancy expected from the variable number of gene copies [77, 78]. Specifically, sequences encoding leucine-rich repeat and NB-ARC domains are common in plant disease resistance proteins, which are particularly diverse because of pathogen pressure [79–81]. Second, the loss of a single member of a gene family may result in a relatively minor loss of the total function of the gene family, because other family members may genetically cushion the impact. Thus, CNV is likely to contribute to quantitative variation rather than qualitative defects in the complex and highly duplicated plant genome [30]. In addition, the gene content in CNVs is speculated to have contributed to heterosis during domestication [28, 30]. Indeed, high levels of variability in gene content among genotypes will result in hybrids that contain a higher number of genes than either parent. We observed in the present study that many of the rice cultivars were missing unequal numbers of genes relative to Nipponbare. It may be possible to identify a series of recombination events in order to combine all superior alleles for heterosis, which would be valuable in breeding programs.

Conclusions

This study provides insights into the mutational mechanisms and functional effects of CNVs in the rice genome. Our results suggest that many CNVs are generated by NAHR events from higher recombination rates and the presence of homologous gene clusters, which is consistent with previous work that chromosome 11 is enriched with disease resistance genes. We also identify candidate CNVs for involvement in group-specific characteristics. This comprehensive catalogue of CNVs will be useful for future studies to uncover the genetic basis of complex traits in rice.

Methods

Oligonucleotide aCGH construction

A custom 3 × 720 k tiling-path aCGH for whole-genome analysis in Oryza sativa (IRGSP v4.0) was designed and constructed by NimbleGen Systems (http://www.nimblegen.com). Probes were synthesized using an isothermal format, varied in length from 50-mer to 75-mer, and spanned the rice genome with a median spacing of ~500 bp. The arrays were manufactured by maskless array synthesis technology, and the oligonucleotides were synthesized on the arrays by photolithography [82, 83].

Sample processing

We selected 20 rice cultivars to represent all of the major groups of Oryza sativa L. [48, 71]. Their origins and features are summarized in Additional file 1: Table S1. The samples comprised two groups of ssp. indica (indica and aus) and four groups of ssp. japonica (rayada, aromatic, tropical japonica and temperate japonica). We also included the ssp. japonica Nipponbare, which was used to generate the rice reference genome sequence [3]. Genomic DNA was extracted and purified from fresh young leaves using the Plant Genomic DNA Kit (TianGen). DNA quality was assessed by measuring the concentration and purity with a NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies). DNA integrity was assessed by electrophoresis in a 1.0% agarose gel. Standard genomic DNA labeling, hybridization, array scanning and intensity feature extraction were carried out as described previously [23, 84] and performed by CapitalBio Corporation.

Statistical analysis

Before normalization and segmentation analysis, spatial correction was applied, which corrected position-dependent non-uniformity of signals across the array. Locally weighted polynomial regression (LOESS) was used to adjust signal intensities based on x, y feature position [85]. Normalization was then performed using the q-spline method [86], compensating for inherent differences in signal between the two dyes. Segmentation analysis was performed with the segMNT algorithm in NimbleScan 2.5 software [87]. The segments were further filtered to remove all segments in which the average log2 ratio for all probes was > −2 or < 1 to produce a set of stringent segments. Aberrant segments were called deletions if the mean log2 ratio of probes in the segment was ≤ −2 and called amplifications with a mean log2 ratio ≥ 1. We used a secondary manual method to examine all aberrant segments and to remove presumptive false-positives caused by single outlier probes and ensure 100% probes of each aberrant segment meet the log2 ratio cut-off. In most cases these adjustments further decreased the count of the CNVs. We merged overlapping CNV coordinates across hybridizations to form unique CNVRs in accordance with similar criteria as described previously [15]. We queried the probes within the CNVRs against the most recent rice genome data release (IRGSP v5.0). It was possible to identify perfect matches (100% identity and 100% coverage) for 99.98% (715,698/715,851) of the probes for the 12 chromosomes. One hundred and fifty-three probes on the array no longer had perfect sequence matches in the genome, 99.92% of the probes (715,130/715,698) had only a single perfect match and were therefore deemed to be single copies, and 568 probes had more than two perfect matches. The probes were only retained if they had a unique optimal hit (100% sequence identity). The CNVRs that decreased the number of probes to less than five were discarded. Finally, the remaining CNVRs were retained if they did not overlap a large gap in IRGSP v5.0.

We inferred approximate allele frequencies by simplifying CNV phenotypes into three categories: normal, loss, and gain. The frequency of each category was estimated by PowerMarker 3.25 software [88]. The estimated allele frequencies were used to calculate heterozygosity (He) [89] for each subspecies and each polymorphic CNVR. Similarly, for each CNVR we calculated Fst as Fst = 1− Hs/Ht, where Hs and Ht denote average heterozygosity within subspecies and total heterozygosity, respectively.

All probe signals of each CNVR were subjected to unbiased clustering for the 20 rice cultivars using average linkage and correlation (uncentered) as the metric. The frequency of CNVs within each group was analyzed using hierarchical clustering. The single case where we identified loss and gain within the same region could have been treated as a multistate locus, but instead we chose to exclude complex events from this analysis. ANOVA was performed with SPSS 16.0.

Functional characterization of genes affected by CNVs

The gene content of CNVRs was examined using the IRGSP v5.0 gene models (http://rgp.dna.affrc.go.jp/IRGSP/Build5/build5.html). GO annotation from RAP-DB (http://rapdb.dna.affrc.go.jp) of genes that were affected by CNVs was assessed using the Gene Ontology Enrichment Analysis Software Toolkit (GOEAST) [90]. P-values for enrichment were calculated using a hypergeometric test method with false discovery rate correction [91]. Of the 42,081 IRGSP v5.0 gene models used, 15,689 (37.28%) had GO terms from RAP-DB (http://rapdb.dna.affrc.go.jp). Of the 1321 gene models in CNVRs, 561 (42.47%) had GO annotations that were used in our GO enrichment analysis.

Genes belonging to gene families across the genome were separated from a single-copy gene set according to http://green.dna.affrc.go.jp/PGF-DB/. The rice-specific genes and gene families were identified based on homolog clustering with annotated genes of maize and Arabidopsis using the method previously described by McGinnis and Madden [92].

Analysis of TEs and repeats in the deleted regions relative to Nipponbare

The regions identified as deletions relative to the reference sample were analyzed for TEs and repeats. Repeats were detected using REPEATMASKER (http://www.repeatmasker.org). Repeat content was expressed as the percentage of nucleotides masked versus the total.

PCR validation

Gene loss events were assessed by PCR. Amplifications were performed on a 2720 Thermal Cycler (Applied Biosystems) under standard conditions for 30 cycles. The PCR products were electrophoresed on 6% non-denaturing polyacrylamide gel.

Access to data

The full data set from the oligonucleotide aCGH experiments has been submitted to the NCBI Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo) [93] under the accession ID GSE42769.

Abbreviations

CNV:

Copy number variant

CNVR:

Copy number variable region

TE:

Transposable element

SNP:

Single-nucleotide polymorphism

LD:

Linkage disequilibrium

GWAS:

Genome-wide association studies

aCGH:

array comparative genomic hybridization

SD:

Segmental duplication

LINE:

Long interspersed nuclear element

SINE:

Short interspersed nuclear element

LTR:

Long terminal repeat

GO:

Gene ontology

NAHR:

Non-allelic homologous recombination.

References

  1. Goff SA, Ricke D, Lan TH, Presting G, Wang R, Dunn M, Glazebrook J, Sessions A, Oeller P, Varma H, Hadley D, Hutchison D, Martin C, Katagiri F, Lange BM, Moughamer T, Xia Y, Budworth P, Zhong J, Miguel T, Paszkowski U, Zhang S, Colbert M, Sun W, Chen L, Cooper B, Park S, Wood TC, Mao L, Quail P: A Draft Sequence of the Rice Genome (Oryza sativa L. ssp. japonica). Science. 2002, 296: 92-100. 10.1126/science.1068275.

    Article  CAS  PubMed  Google Scholar 

  2. Yu J, Hu S, Wang J, Wong GKS, Li S, Liu B, Deng YJ, Dai L, Zhou Y, Zhang XQ, Cao ML, Liu J, Sun JD, Tang JB, Chen YJ, Huang XB, Lin W, Ye C, Tong W, Cong LJ, Geng JN, Han YJ, Li L, Li W, Hu GQ, Huang XG, Li WJ, Li J, Liu ZW, Li L: A Draft Sequence of the Rice Genome (Oryza sativa L. ssp. indica). Science. 2002, 296: 79-92. 10.1126/science.1068037.

    Article  CAS  PubMed  Google Scholar 

  3. International Rice Genome Sequencing Project: The map-based sequence of the rice genome. Nature. 2005, 436: 793-800. 10.1038/nature03895.

    Article  Google Scholar 

  4. Hinds DA, Stuve LL, Nilsen GB, Halperin E, Eskin E, Ballinger DG, Frazer KA, Cox DR: Whole-genome patterns of common DNA variation in three human populations. Science. 2005, 307: 1072-1079. 10.1126/science.1105436.

    Article  CAS  PubMed  Google Scholar 

  5. Frazer KA, Eskin E, Kang HM, Bogue MA, Hinds DA, Beilharz EJ, Gupta RV, Montgomery J, Morenzoni MM, Nilsen GB, Pethiyagoda CL, Stuve LL, Johnson FM, Daly MJ, Wade CM, Cox DR: A sequence-based variation map of 8.27 million SNPs in inbred mouse strains. Nature. 2007, 448: 1050-1053. 10.1038/nature06067.

    Article  CAS  PubMed  Google Scholar 

  6. Kim S, Plagnol V, Hu TT, Toomajian C, Clark RM, Ossowski S, Ecker JR, Weigel D, Nordborg M: Recombination and linkage disequilibrium in Arabidopsis thaliana. Nat Genet. 2007, 39: 1151-1155. 10.1038/ng2115.

    Article  CAS  PubMed  Google Scholar 

  7. McNally KL, Childs KL, Bohnert R, Davidson RM, Zhao K, Ulat VJ, Zeller G, Clark RM, Hoen DR, Bureau TE, Stokowski R, Ballinger DG, Frazer KA, Cox DR, Padhukasahasram B, Bustamante CD, Weigel D, Mackill DJ, Bruskiewich RM, Rätsch G, Buell CR, Leung H, Leach JE: Genomewide SNP variation reveals relationships among landraces and modern varieties of rice. Proc Natl Acad Sci. 2009, 106: 12273-12278. 10.1073/pnas.0900992106.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. McCouch SR, Zhao K, Wright M, Tung CW, Ebana K, Thomson M, Reynolds A, Wang D, DeClerck G, Ali ML, McClung A, Eizenga G, Bustamante C: Development of genome-wide SNP assays for rice. Breeding Sci. 2010, 60: 524-535. 10.1270/jsbbs.60.524. Carlos Bustamante

    Article  Google Scholar 

  9. Tung CW, Zhao K, Wright MH, Ali ML, Jung J, Tung CW, Zhao KY, Wright MH, Ali ML, Jung J, Kimball J, Tyagi W, Thomson MJ, McNally K, Leung H: Development of a research platform for dissecting phenotype-genotype associations in rice (Oryza spp.). Rice. 2010, 3: 205-217. 10.1007/s12284-010-9056-5.

    Article  Google Scholar 

  10. Ali ML, McClung AM, Jia MH, Kimball JA, McCouch SR, Georgia CE: Rice diversity panel evaluated for genetic and agro-morphological diversity between subpopulations and its geographic distribution. Crop Sci. 2011, 51: 2021-2035. 10.2135/cropsci2010.11.0641.

    Article  Google Scholar 

  11. Huang XH, Wei XH, Sang T, Zhao Q, Feng Q, Zhao Y, Li CY, Zhu CR, Lu TT, Zhang ZW, Li M, Fan DL, Guo YL, Wang AH, Wang L, Deng LW, Li WJ, Lu YQ, Weng QJ, Liu KY, Huang T, Zhou TY, Jing YF, Li W, Lin Z, Buckler ES, Qian Q, Zhang QF, Li JY, Han B: Genome-wide association studies of 14 agronomic traits in rice landraces. Nat Genet. 2010, 42: 961-967. 10.1038/ng.695.

    Article  CAS  PubMed  Google Scholar 

  12. Zhao K, Tung CW, Eizenga GC, Wright MH, Ali ML, Price AH, Norton GJ, Islam MR, Reynolds A, Mezey J, Anna MMC, Bustamante CD, McCouch SR: Genome-wide association mapping reveals rich genetic architecture of complex traits in Oryza sativa. Nat Commun. 2011, 2: 467-

    Article  PubMed Central  PubMed  Google Scholar 

  13. Huang XH, Zhao Y, Wei XH, Li C, Wang A, Zhao Q, Li WJ, Guo YL, Deng LW, Zhu CR, Fan DL, Lu YQ, Weng QJ, Liu KY, Zhou TY, Jing YF, Si LZ, Dong GJ, Huang T, Lu TT, Feng Q, Qian Q, Li JY, Han B: Genome-wide association study of flowering time and grain yield traits in a worldwide collection of rice germplasm. Nat Genet. 2012, 44: 32-39.

    Article  Google Scholar 

  14. Iafrate AJ, Feuk L, Rivera MN, Listewnik ML, Donahoe PK, Qi Y, Scherer SW, Lee C: Detection of large-scale variation in the human genome. Nat Genet. 2004, 36: 949-951. 10.1038/ng1416.

    Article  CAS  PubMed  Google Scholar 

  15. Redon R, Ishikawa S, Fitch KR, Feuk L, Perry GH, Andrews TD, Fiegler H, Shapero MH, Carson AR, Chen W, Cho EK, Dallaire S, Freeman JL, González JR, Gratacòs M, Huang J, Kalaitzopoulos D, Komura D, MacDonald JR, Marshall CR, Mei R, Montgomery L, Nishimura K, Okamura K, Shen F, Somerville MJ, Tchinda J, Valsesia A, Woodwark C, Yang F: Global variation in copy number in the human genome. Nature. 2006, 444: 444-454. 10.1038/nature05329.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Perry GH, Yang F, Marques-Bonet T, Murphy C, Fitzgerald T, Lee AS, Hyland C, Stone AC, Hurles ME, Chris TS, Eichler EE, Carter NP, Lee C, Redon R: Copy number variation and evolution in humans and chimpanzees. Genome Res. 2008, 18: 1689-1710.

    Article  Google Scholar 

  17. Chen WK, Swartz JD, Rush LJ, Alvarez CE: Mapping DNA structural variation in dogs. Genome Res. 2009, 19: 500-509.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  18. Nicholas TJ, Baker C, Eichler EE, Akey JM: A high-resolution integrated map of copy number polymorphisms within and between breeds of the modern domesticated dog. BMC Genomics. 2011, 12: 414-10.1186/1471-2164-12-414.

    Article  PubMed Central  PubMed  Google Scholar 

  19. Wang X, Nahashon S, Feaster TK, Bohannon-Stewart A, Adefope N: An initial map of chromosomal segmental copy number variations in the chicken. BMC Genomics. 2010, 11: 351-10.1186/1471-2164-11-351.

    Article  PubMed Central  PubMed  Google Scholar 

  20. Liu GE, Hou Y, Zhu B, Cardone MF, Jiang L, Cellamare A, Mitra A, Alexander LJ, Coutinho LL, Dell’Aquila ME, Gasbarre LC, Lacalandra G, Li RW, Matukumalli LK, Nonneman D, Regitano LC, Smith TP, Song J, Sonstegard TS, Van Tassell CP, Ventura M, Eichler EE, McDaneld TG, Keele JW: Analysis of copy number variations among diverse cattle breeds. Genome Res. 2010, 20: 693-703. 10.1101/gr.105403.110.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  21. Fadista J, Nygaard M, Holm LE, Thomsen B, Bendixen C: A snapshot of CNVs in the pig genome. PLoS ONE. 2008, 3: e3916-10.1371/journal.pone.0003916.

    Article  PubMed Central  PubMed  Google Scholar 

  22. Guryev V, Saar K, Adamovic T, Verheul M, van Heesch SAAC, Cook S, Pravenec M, Aitman T, Jacob H, Shull JD, Hubner N, Cuppen E: Distribution and functional impact of DNA copy number variation in the rat. Nat Genet. 2008, 40: 538-545. 10.1038/ng.141.

    Article  CAS  PubMed  Google Scholar 

  23. Graubert TA, Cahan P, Edwin D, Selzer RR, Richmond TA, Eis PS, Shannon WD, Li X, McLeod HL, Cheverud JM, Ley TJ: A high-resolution map of segmental DNA copy number variation in the mouse genome. PLoS Genet. 2007, 3: e3-10.1371/journal.pgen.0030003.

    Article  PubMed Central  PubMed  Google Scholar 

  24. Emerson JJ, Cardoso-Moreira M, Borevitz JO, Long M: Natural selection shapes genome-wide patterns of copy-number polymorphism in Drosophila melanogaster. Science. 2008, 320: 1629-1631. 10.1126/science.1158078.

    Article  CAS  PubMed  Google Scholar 

  25. Maydan JS, Lorch A, Edgley ML, Flibotte S, Moerman DG: Copy number variation in the genomes of twelve natural isolates of Caenorhabditis elegans. BMC Genomics. 2010, 11: 62-10.1186/1471-2164-11-62.

    Article  PubMed Central  PubMed  Google Scholar 

  26. Infante JJ, Dombek KM, Rebordinos L, Cantoral JM, Young ET: Genome-wide amplifications caused by chromosomal rearrangements play a major role in the adaptive evolution of natural yeast. Genetics. 2003, 165: 1745-1759.

    PubMed Central  CAS  PubMed  Google Scholar 

  27. Skvortsov D, Abdueva D, Stitzer ME, Finkel SE, Tavaré S: Using expression arrays for copy number detection: an example from E. coli. BMC Bioinforma. 2007, 8: 203-10.1186/1471-2105-8-203.

    Article  Google Scholar 

  28. Springer NM, Ying K, Fu Y, Ji T, Yeh CT, Jia Y, Wu W, Richmond T, Kitzman J, Rosenbaum H, Iniguez AL, Barbazuk WB, Jeddeloh JA, Nettleton D, Schnable PS: Maize inbreds exhibit high levels of copy number variation (CNV) and presence/absence variation (PAV) in genome content. PLoS Genet. 2009, 5: e1000734-10.1371/journal.pgen.1000734.

    Article  PubMed Central  PubMed  Google Scholar 

  29. Beló A, Beatty MK, Hondred D, Fengler KA, Li B, Rafalski A: Allelic genome structural variations in maize detected by array comparative genome hybridization. Theor Appl Genet. 2009, 120: 355-367.

    Article  PubMed  Google Scholar 

  30. Swanson-Wagner RA, Eichten SR, Kumari S, Tiffin P, Stein JC, Ware D, Springer NM: Pervasive gene content variation and copy number variation in maize and its undomesticated progenitor. Genome Res. 2010, 20: 1689-1699. 10.1101/gr.109165.110.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. DeBolt S: Copy number variation shapes genome diversity in Arabidopsis over immediate family generational scales. Genome Biol Evol. 2010, 2: 441-453. 10.1093/gbe/evq033.

    Article  PubMed Central  PubMed  Google Scholar 

  32. Santuari L, Pradervand S, Amiguet-Vercher AM, Thomas J, Dorcey E, Harshman K, Xenarios I, Juenger TE, Hardtke CS: Substantial deletion overlap among divergent Arabidopsis genomes revealed by intersection of short reads and tiling arrays. Genome Biol. 2010, 11: R4-10.1186/gb-2010-11-1-r4.

    Article  PubMed Central  PubMed  Google Scholar 

  33. Feuk L, Carson AR, Scherer SW: Structural variation in the human genome. Nat Rev Genet. 2006, 7: 85-97.

    Article  CAS  PubMed  Google Scholar 

  34. Levy S, Sutton G, Ng PC, Feuk L, Halpern AL, Walenz BP, Axelrod N, Huang JQ, Kirkness EF, Denisov G, Lin Y, MacDonald JR, Pang AWC, Shago M, Stockwell TB, Tsiamouri A, Bafna V, Bansal V, Kravitz SA, Busam DA, Beeson KY, McIntosh TC, Remington KA, Abril F, Gill J, Borman J, Rogers YH, Frazier ME, Scherer SW, Strausberg RL: The diploid genome sequence of an individual human. PLoS Biol. 2007, 5: e254-10.1371/journal.pbio.0050254.

    Article  PubMed Central  PubMed  Google Scholar 

  35. Wong KK, DeLeeuw RJ, Dosanjh NS, Kimm LR, Cheng Z, Horsmac DE, MacAulayb C, Ngd RT, Browne CJ, Eichlerf EE, Lama WL: A comprehensive analysis of common copy-number variations in the human genome. Am J Hum Genet. 2007, 80: 91-104. 10.1086/510560.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  36. Zogopoulos G, Ha KC, Naqib F, Moore S, Kim H, Zogopoulos G, Ha KCH, Naqib F, Moore S, Kim H, Montpetit A, Robidoux F, Laflamme P, Cotterchio M, Greenwood C, Scherer SW, Zanke B, Hudson TJ, Bader GD, Gallinger S: Germ-line DNA copy number variation frequencies in a large north american population. Hum Genet. 2007, 122: 345-353. 10.1007/s00439-007-0404-5.

    Article  CAS  PubMed  Google Scholar 

  37. Perry GH, Ben-Dor A, Tsalenko A, Sampas N, Rodriguez-Revenga L, Tran CW, Scheffer A, Steinfeld I, Tsang P, Yamada NA, ParkH AS, Kim JI, Seo JS, Yakhini Z, Laderman S, Bruhn L, Lee C: The fine-scale and complex architecture of human copy-number variation. Am J Hum Genet. 2008, 82: 685-695. 10.1016/j.ajhg.2007.12.010.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  38. Shaikh TH, Gai X, Perin JC, Glessner JT, Xie H, Murphy K, Hara RO, Casalunovo T, Conlin LK, Arcy MD, Frackelton EC, Geiger EA, Englert CH, Imielinski M, Kim CE, Medne L, Annaiah K, Bradfield JP, Dabaghyan E, Eckert A, Onyiah CC, Ostapenko S, Otieno FG, Santa E, Shaner JL, Skraban R, Smith RM, Elia J, Goldmuntz E, Spinner NB: High-resolution mapping and analysis of copy number variations in the human genome: a data resource for clinical and research applications. Genome Res. 2009, 19: 1682-1690. 10.1101/gr.083501.108.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  39. Lupski JR: Genomic rearrangements and sporadic disease. Nat Genet. 2007, 39: S43-S47. 10.1038/ng2084.

    Article  CAS  PubMed  Google Scholar 

  40. Hinds DA, Kloek AP, Jen M, Chen X, Frazer KA: Common deletions and SNPs are in linkage disequilibrium in the human genome. Nat Genet. 2006, 38: 82-85. 10.1038/ng1695.

    Article  CAS  PubMed  Google Scholar 

  41. Locke DP, Sharp AJ, McCarroll SA, McGtath SD, Newman TL, Cheng Z, Schwartz S, Albertson DG, Pinkel D, Altshuler DM, Eichler EE: Linkage disequilibrium and heritability of copy-number polymorphisms within duplicated regions of the human genome. Am J Hum Genet. 2006, 79: 275-290. 10.1086/505653.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  42. McCarroll SA, Hadnott TN, Perry GH, Sabeti PC, Zody MC, Barrett JC, Dallaire S, Gabriel SB, Lee C, Daly MJ, Altshuler DM: The International HapMap Consortium Common deletion polymorphisms in the human genome. Nat Genet. 2006, 38: 86-92. 10.1038/ng1696.

    Article  CAS  PubMed  Google Scholar 

  43. Yu P, Wang C, Xu Q, Feng Y, Yuan X, Yu H, Wang Y, Tang S, Wei X: Detection of copy number variations in rice using array-based comparative genomic hybridization. BMC Genomics. 2011, 12: 372-10.1186/1471-2164-12-372.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  44. Xu X, Liu X, Ge S, Jensen JD, Hu F, Li X, Dong Y, Gutenkunst RN, Fang L, Huang L, Li J, He W, Zhang G, Zheng X, Zhang F, Li Y, Yu C, Kristiansen K, Zhang X, Wang J, Susan MW, Nielsen MR, Wang J, Wang W: Resequencing 50 accessions of cultivated and wild rice yields markers for identifying agronomically important genes. Nature biotechnology. 2012, 30: 105-111.

    Article  CAS  Google Scholar 

  45. Li SC, Wang SQ, Deng QM, Zheng AP, Zhu J, Liu HN, Wang LX, Gao FY, Zou T, Huang B, Cao XM, Xu LZ, Yu C, Ai P, Li P: Identification of genome-wide variations among three elite restorer lines for hybrid-rice. PloS ONE. 2012, 7: e30952-10.1371/journal.pone.0030952.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  46. Fadista J, Thomsen B, Holm LE, Bendixen C: Copy number variation in the bovine genome. BMC Genomics. 2010, 11: 284-10.1186/1471-2164-11-284.

    Article  PubMed Central  PubMed  Google Scholar 

  47. Lee AS, Gutierrez-Arcelus M, Perry GH, Vallender EJ, Johnson WE, Miller GM, Korbel JO, Lee C: Analysis of copy number variation in the rhesus macaque genome identifies candidate loci for evolutionary and human disease studies. Hum Mol Genet. 2008, 17: 1127-1136. 10.1093/hmg/ddn002.

    Article  CAS  PubMed  Google Scholar 

  48. Garris AJ, Tai TH, Coburn J, Kresovich S, McCouch S: Genetic Structure and Diversity in Oryza sativa L. Genetics. 2005, 169: 1631-1638.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  49. Henrichsen CN, Vinckenbosch N, Zollner S, Chaignat E, Pradervand S, Schütz F, Ruedi M, Kaessmann H, Reymond A: Segmental copy number variation shapes tissue transcriptomes. Nat Genet. 2009, 41: 424-429. 10.1038/ng.345.

    Article  CAS  PubMed  Google Scholar 

  50. Sun CQ, Wang XK, Li ZC, Yoshimura A, Lwata N: Comparison of the genetic diversity of common wild rice (Oryza rufipogon Griff.) and cultivated rice (O. sativa L.) using RFLP markers. Theor Appl Genet. 2001, 102: 157-162. 10.1007/s001220051631.

    Article  CAS  Google Scholar 

  51. Gao LZ, Zhang CH, Chang LP, Jia JZ, Qiu ZE, Dong YS: Microsatellite diversity within Oryza sativa with emphasis on indica-japonica divergence. Genet Res. 2005, 85: 1-14. 10.1017/S0016672304007293.

    Article  CAS  PubMed  Google Scholar 

  52. Thomson MJ, Septiningsih EM, Suwardjo F, Santoso TJ, Silitonga TS, McCouch SR: Genetic diversity analysis of traditional and improved Indonesian rice (Oryza sativa L.) germplasm using microsatellite markers. Theor Appl Genet. 2007, 114: 559-568. 10.1007/s00122-006-0457-1.

    Article  CAS  PubMed  Google Scholar 

  53. Gao LZ, Innan H: Nonindependent domestication of the two Rice subspecies, Oryza sativa ssp. indica and ssp. japonica, demonstrated by multilocus microsatellites. Genetics. 2008, 179: 965-976. 10.1534/genetics.106.068072.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  54. Agrama HA, Yan WG, Jia M, Fjellstrom R, McClung AM: Genetic structure associated with diversity and geographic distribution in the USDA rice world collection. Nat Sci. 2010, 2: 247-291.

    Google Scholar 

  55. Zhao K, Wright M, Kimball J, Eizenga G, McClung A, Kovach M, Tyagi W, Ali ML, Tung CW, Reynolds A, Bustamante CD, McCouch SR: Genomic diversity and introgression in O. sativa reveal the impact of domestication and breeding on the rice genome. PLoS ONE. 2010, 5: e10780-10.1371/journal.pone.0010780.

    Article  PubMed Central  PubMed  Google Scholar 

  56. Siezen RJ, Leunissen JAM: Subtilases: The superfamily of subtilisin-like serine proteases. Protein Sci. 1997, 6: 501-523.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  57. Sharp AJ, Locke DP, McGrath SD, Cheng Z, Bailey JA, Vallente RU, Pertz LM, Clark RA, Schwartz S, Seqraves R, Oseroff W, Albertson DG, Pinkel D, Eichler EE: Segmental duplications and copy-number variation in the human genome. Am J Hum Genet. 2005, 77: 78-88. 10.1086/431652.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  58. Sebat J, Lakshmi B, Troge J, Alexander J, Young J, Lundin P, Månér S, Massa H, Walker M, Chi M, Navin N, Lucito R, Healy J, Hicks J, Ye K, Reiner A, Gilliam TC, Trask B, Patterson N, Zetterberg A, Wigler M: Large-scale copy number polymorphism in the human genome. Science. 2004, 305: 525-528. 10.1126/science.1098918.

    Article  CAS  PubMed  Google Scholar 

  59. Tuzun E, Sharp AJ, Bailey JA, Kaul R, Morrison VA, Pertz LM, Haugen E, Hayden H, Albertson D, Pinkel D, Olson MV, Eichler EE: Fine-scale structural variation of the human genome. Nat Genet. 2005, 37: 727-732. 10.1038/ng1562.

    Article  CAS  PubMed  Google Scholar 

  60. Li J, Jiang T, Mao JH, Balmain A, Peterson L, Harris C, Rao PH, Havlak P, Gibbs R, Cai WW: Genomic segmental polymorphisms in inbred mouse strains. Nat Genet. 2004, 36: 952-954. 10.1038/ng1417.

    Article  CAS  PubMed  Google Scholar 

  61. Cutler G, Marshall LA, Chin N, Baribault H, Kassner PD: Significant gene content variation characterizes the genomes of inbred mouse strains. Genome Res. 2007, 12: 1743-1754.

    Article  Google Scholar 

  62. Watkins-Chow DE, Pavan WJ: Genomic copy number and expression variation within the C57BL/6J inbred mouse strain. Genome Res. 2008, 18: 60-66.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  63. Hurwitz BL, Kudrna D, Yu Y, Sebastian A, Zuccolo A, Jackson SA, Ware D, Wing RA, Stein L: Rice structural variation: a comparative analysis of structural variation between rice and three of its closest relatives in the genus Oryza. Plant J. 2010, 63: 990-1003. 10.1111/j.1365-313X.2010.04293.x.

    Article  CAS  PubMed  Google Scholar 

  64. Yu J, Wang J, Lin W, Li S, Li H, Zhou J, Ni P, Dong W, Hu S, Zeng C, Zhang J, Zhang Y, Li R, Xu Z, Li S, Li X, Zheng H, Cong L, Lin L, Yin J, Geng J, Li G, Shi J, Liu J, Lv H, Li J, Wang J, Deng Y, Ran L, Shi X: The genomes of Oryza sativa: a history of duplications. PLoS Biol. 2005, 3: e38-10.1371/journal.pbio.0030038.

    Article  PubMed Central  PubMed  Google Scholar 

  65. Paterson AH, Chapman BA, Kissinger JC, Bowers JE, Feltus FA, Estill JC: Many gene and domain families have convergent fates following independent whole-genome duplication events in Arabidopsis, Oryza, Saccharomyces and Tetraodon. Trends Genet. 2006, 22: 597-602. 10.1016/j.tig.2006.09.003.

    Article  CAS  PubMed  Google Scholar 

  66. The Rice Chromosomes 11 and 12 Sequencing Consortia: The sequence of rice chromosomes 11 and 12, rich in disease resistance genes and recent gene duplications. BMC Biol. 2005, 3: 20-37.

    Article  PubMed Central  Google Scholar 

  67. Conrad DF, Hurles ME: The population genetics of structural variation. Nat Genet. 2007, 39: S30-S36. 10.1038/ng2042.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  68. Oka HI: Indica-Japonica differentiation of rice cultivars. Origin of Cultivated Rice. 1988, Tokyo: Elsevier Science/Japan Scientific Societies Press, 141-179.

    Google Scholar 

  69. Second G: Origin of the genic diversity of cultivated rice (Oryza spp.): study of the polymorphism scored at 40 isozyme loci. Jpn J Genet. 1982, 57: 25-57. 10.1266/jjg.57.25.

    Article  Google Scholar 

  70. Cheng CY, Motohashi R, Tsuchimoto S, Fukuta Y, Ohtsubo H, Ohtsubo E: Polyphyletic origin of cultivated rice: based on the interspersion pattern of SINEs. Mol Biol Evol. 2003, 20: 67-75. 10.1093/molbev/msg004.

    Article  CAS  PubMed  Google Scholar 

  71. Glaszmann JC: Isozymes and classification of asian rice varieties. Theor Appl Genet. 1987, 74: 21-30. 10.1007/BF00290078.

    Article  CAS  PubMed  Google Scholar 

  72. Caicedo AL, Williamson SH, Hernandez RD, Boyko A, Fledel-Alon A, York TL, Polato NR, Olsen KM, Nielsen R, McCouch SR, Bustamante CD, Purugganan MD: Genome-wide patterns of nucleotide polymorphism in domesticated rice. PLoS Genet. 2007, 3: e163-10.1371/journal.pgen.0030163.

    Article  PubMed Central  Google Scholar 

  73. Khush G: Origin, dispersal, cultivation and variation of rice. Plant Mol Biol. 1997, 35: 25-34. 10.1023/A:1005810616885.

    Article  CAS  PubMed  Google Scholar 

  74. Henrichsen CN, Chaignat E, Reymond A: Copy number variants, diseases and gene expression. Hum Mol Genet. 2009, 18: R1-R8. 10.1093/hmg/ddp011.

    Article  CAS  PubMed  Google Scholar 

  75. Zhang F, Gu W, Hurles ME, Lupski JR: Copy number variation in human health, disease, and evolution. Annu Rev Genomics Hum Genet. 2009, 10: 451-481. 10.1146/annurev.genom.9.081307.164217.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  76. Shomura A, Izawa T, Ebana K, Ebitani T, Kanegae H, Konishi S, Yano M: Deletion in a gene associated with grain size increased yields during rice domestication. Nat Genet. 2008, 40: 1023-1028. 10.1038/ng.169.

    Article  CAS  PubMed  Google Scholar 

  77. Ohno S: Evolution by gene duplication. 1970, New York: Springer-Verlag

    Book  Google Scholar 

  78. Kondrashov FA, Kondrashov AS: Role of selection in fixation of gene duplications. J Theor Biol. 2006, 239: 141-151. 10.1016/j.jtbi.2005.08.033.

    Article  CAS  PubMed  Google Scholar 

  79. Grant MR, McDowell JM, Sharpe AG, Zabala MT, Lydiate DJ, Dangl JL: Independent deletions of a pathogen-resistance gene in Brassica and Arabidopsis. Proc Natl Acad Sci USA. 1998, 95: 15843-15848. 10.1073/pnas.95.26.15843.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  80. Shen J, Araki H, Chen L, Chen J, Tian D: Unique evolutionary mechanism in R genes under the presence/absence polymorphism in Arabidopsis thaliana. Genetics. 2006, 172: 1243-1250.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  81. Clark RM, Schweikert G, Toomajian C, Ossowski S, Zeller G, Shinn P, Warthmann N, Hu TT, Fu G, Hinds DA, Chen H, Frazer KA, Huson DH, Scholkopf B, Nordborg M, Ratsch G, Ecker JR, Weigel D: Common sequence polymorphisms shaping genetic diversity in Arabidopsis thaliana. Science. 2007, 317: 338-342. 10.1126/science.1138632.

    Article  CAS  PubMed  Google Scholar 

  82. Singh-Gasson S, Green RD, Yue Y, Nelson C, Blattner F, Sussman MR, Cerrina F: Maskless fabrication of light-directed oligonucleotide microarrays using a digital micromirror array. Nat Biotech. 1999, 17: 974-978. 10.1038/13664.

    Article  CAS  Google Scholar 

  83. Nuwaysir EF, Huang W, Albert TJ, Singh J, Nuwaysir K, Pitas A, Richmond T, Gorski T, Berg JP, Ballin J, McCormick M, Norton J, Pollock T, Sumwalt T, Butcher L, Porter D, Molla M, Hall C, Blattner F, Sussman MR, Wallace RL, Cerrina F, Green RD: Gene expression analysis using oligonucleotide arrays produced by maskless photolithography. Genome Res. 2002, 12: 1749-1755. 10.1101/gr.362402.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  84. Selzer RR, Richmond TA, Pofahl NJ, Green RD, Eis PS, Nair P, Brothman AR, Stallings R: Analysis of chromosome breakpoints in neuroblastoma at sub-kilobase resolution using fine-tiling oligonucleotide array CGH. Gene Chromosomes Canc. 2005, 44: 305-319. 10.1002/gcc.20243.

    Article  CAS  Google Scholar 

  85. Smyth GK, Speed TP: Normalization of cDNA microarray data. Methods. 2003, 31: 265-273. 10.1016/S1046-2023(03)00155-5.

    Article  CAS  PubMed  Google Scholar 

  86. Workman C, Jensen LJ, Jarmer H, Berka R, Gautier L, Nielsen HB, Saxild HH, Nielsen C, Brunak S, Knudsen S: A new non-linear normalization method for reducing variability in DNA microarray experiments. Genome Biol. 2002, 3: 1-16.

    Article  Google Scholar 

  87. Molla M: Novel uses for machine learning and other computational methods for the design and interpretation of genetic microarrays. PhD thesis. 2007, Madison: USA: University of Wisconsin-Madison, Computer Science Department

    Google Scholar 

  88. Liu K, Muse SV: Integrated analysis environment for genetic marker data. Bioinformatics. 2005, 21: 2128-2129. 10.1093/bioinformatics/bti282.

    Article  CAS  PubMed  Google Scholar 

  89. Nei M: Molecular Evolutionary Genetic. 1987, New York: Columbia University Press, 190-191.

    Google Scholar 

  90. Zheng Q, Wang XJ: GOEAST: a web-based software toolkit for Gene Ontology enrichment analysis. Nucleic Acids Res. 2008, 36: W358-W363. 10.1093/nar/gkn276.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  91. Benjamini Y, Yekutieli D: The control of the false discovery rate in multiple testing under dependency. Ann Stat. 2001, 29: 1165-1188. 10.1214/aos/1013699998.

    Article  Google Scholar 

  92. McGinnis S, Madden TL: BLAST: at the core of a powerful and diverse set of sequence analysis tools. Nucleic Acids Res. 2004, 32: W20-W25. 10.1093/nar/gkh435.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  93. Barrett T, Troup DB, Wilhite SE, Ledoux P, Rudnev D, Evangelista C, Kim IF, Soboleva A, Tomashevsky M, Marshall KA, Phillippy KH, Sherman PM, Muertter RN, Edgar R: NCBI GEO: archive for high-throughput functional genomic data. Nucleic Acids Res. 2009, 37: D885-D890. 10.1093/nar/gkn764.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgments

We thank Zhenglin Du (Beijing Institute of Genomics, Chinese Academy of Sciences) for computational assistance. We also thank Bin Han (Institute of Plant Physiology and Ecology, Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences) for critical reading of the manuscript. This work was supported by the Crop Germplasm Resources Protection Project of MOA, China (NB2012-2130135), and the National Key Technology R&D Program in China (2013BAD01B02-14).

Author information

Authors and Affiliations

  1. State Key Laboratory of Rice Biology, China National Rice Research Institute, Hangzhou, China

    Ping Yu, Cai-Hong Wang, Qun Xu, Yue Feng, Xiao-Ping Yuan, Han-Yong Yu, Yi-Ping Wang, Sheng-Xiang Tang & Xing-Hua Wei

Authors
  1. Ping Yu
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Cai-Hong Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Qun Xu
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Yue Feng
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Xiao-Ping Yuan
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Han-Yong Yu
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Yi-Ping Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Sheng-Xiang Tang
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Xing-Hua Wei
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Xing-Hua Wei.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

PY and XHW conceived and designed the experiments. YF, XPY, HYY and YPW performed the experiments. PY, CHW, and QX contributed to the interpretation of the data. PY, XHW, and SXT drafted the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

Additional file 1: Table S1: Selection of rice samples. (XLS 26 KB)

Additional file 2: Table S2: Sample level CNV calls from aCGH. (XLS 2 MB)

Additional file 3: Table S3: IRGSP v5.0 CNVRs and their frequencies. (XLS 364 KB)

Additional file 4: Table S4: CNVRs in ssp. indica and ssp. japonica.(XLS 517 KB)

Additional file 5: Table S5: CNVR coverage in all chromosomes. (XLS 32 KB)

12864_2013_5351_MOESM6_ESM.xls

Additional file 6: Table S6: Locations and classifications of repeats found in the deletions relative to the reference sample. (XLS 550 KB)

Additional file 7: Table S7: CNV frequency difference between subspecies. (XLS 399 KB)

Additional file 8: Table S8: High confidence CNVR frequency differences among groups. (XLS 46 KB)

Additional file 9: Table S9: Primers used for PCR validation of gene loss events. (XLS 29 KB)

Additional file 10: Table S10: Gene content of 2886 CNVRs. (XLS 248 KB)

Additional file 11: Table S11: Conservation analysis of 1321 genes in Arabidopsis and maize. (XLS 219 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Yu, P., Wang, CH., Xu, Q. et al. Genome-wide copy number variations in Oryza sativa L.. BMC Genomics 14, 649 (2013). https://doi.org/10.1186/1471-2164-14-649

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-14-649

Keywords

BMC Genomics

ISSN: 1471-2164

Contact us