IC4R004-Epigenomic-2012-22778444 - Revision history

Xysj1990: /* Research Findings */

2016-06-22T07:55:53Z

‎Research Findings

Xysj1990: /* Research Findings */

2016-06-22T07:55:45Z

‎Research Findings

Xysj1990: /* Research Findings */

2016-06-22T07:55:27Z

‎Research Findings

Xysj1990: /* Research Findings */

2016-06-22T07:55:06Z

‎Research Findings

Xysj1990: /* Research Findings */

2016-06-22T07:52:40Z

‎Research Findings

Xysj1990: /* Research Findings */

2016-06-22T07:51:02Z

‎Research Findings

Xysj1990: /* Research Findings */

2016-06-22T07:50:05Z

‎Research Findings

Xysj1990 at 07:48, 22 June 2016

2016-06-22T07:48:14Z

Xysj1990: /* Plant Materials & Treatment */

2016-06-22T04:48:32Z

‎Plant Materials & Treatment

Xysj1990: /* Plant Materials & Treatment */

2016-06-22T04:48:18Z

‎Plant Materials & Treatment

@@ Line 12: / Line 12: @@
 [[File:IC4R004-Epigenomic-2012-22778444-1.png |center |thumb |927px |'''Figure 1. SNPs between NPB and 93–11. (A) Identification of SNPs required a minimum of three reads per strand in each parental strain. Although there is ambiguity because of cytosine conversion in bisulfite-treated DNA, we can resolve this by considering the sequence of the reverse strand, and, as a result, are able to call all possible types of SNPs. In the example, X indicates the location of the SNP and reads represent, in order, NPB forward strand, NPB reverse strand, 93–11 forward strand, and 93–11 reverse strand. (B) SNP distribution over the genome is plotted on the bottom in windows of 10,000 bases. Each color change represents a new chromosome. The top line represents the density of genes.''']]
 <br>
 * We found one SNP every 253 base pairs between the two parental strains (Fig. 1B), representing a level of divergence that is similar to previous japonica–indica comparisons (13). Not surprisingly, regions enriched with genes had lower densities of SNPs. Of these SNPs, the most frequent were C to T changes, or their complements corresponding to G to A mutations on the opposite strand. Combined, these types of mutations accounted for 73% of the SNPs that were found between varieties (Fig. 1C), which represents a threefold enrichment compared to the expected value if all types of mutations were equally likely. We also observe that methylated cytosines mutate more than three times more frequently than nonmethylated cytosines(Fig. 1D), and they most often mutate to thymines. This is consistent with previous observations that the deamination of methylcytosine, which results in thymine, occurs at higher rates than other spontaneous mutations (14–16).
 <br>

@@ Line 14: / Line 14: @@
 * We found one SNP every 253 base pairs between the two parental strains (Fig. 1B), representing a level of divergence that is similar to previous japonica–indica comparisons (13). Not surprisingly, regions enriched with genes had lower densities of SNPs. Of these SNPs, the most frequent were C to T changes, or their complements corresponding to G to A mutations on the opposite strand. Combined, these types of mutations accounted for 73% of the SNPs that were found between varieties (Fig. 1C), which represents a threefold enrichment compared to the expected value if all types of mutations were equally likely. We also observe that methylated cytosines mutate more than three times more frequently than nonmethylated cytosines(Fig. 1D), and they most often mutate to thymines. This is consistent with previous observations that the deamination of methylcytosine, which results in thymine, occurs at higher rates than other spontaneous mutations (14–16).
 [[File:IC4R004-Epigenomic-2012-22778444-2.png |center |thumb |927px |'''Figure 1. (C) The fraction of the various types of SNPs found in NPB versus 93–11 is shown with the CT and GA SNPs grouped together, and they represent nearly 74% of all SNPs. (D) The relative frequency of each SNP between NPB and 93–11 is shown per 1,000 bases. Methyl-Cs convert to Ts at a rate that is nearly five times that of unmethylated Cs.''']]
 * Next, the researchers compared mutation rates between the parental strains and the hybrids. Not surprisingly, we found that SNPs were nearly absent across generations. If, as is likely, the few SNPs we call are false positives, we estimate the false discovery rate of our SNP calls to be less than 0.0003%.
 * The researchers also found a small number of SNPs between the reference Michigan State University O. sativa japonica genome and our as- sembled version (1∕90651 bp; 2,487 total) (Dataset S2), and considerably more SNPs by comparing our assembled version of the 93–11 genome with the reference Beijing Genome Institute’s O. sativa ssp. indica genome (1∕51 bp), consistent with the fact that the 93–11 assembly is not as complete as that of the NPB cultivar (17). The fact that we only align reads that map to a unique posi- tion in the genome, have at most three mismatches with the reference genome, do not permit insertions or deletions, and re- quire a minimum of three reads on both strands, allowed us to only call bases over about half of the genome in both NPB and 93–11, and thus we are reporting only a subset of all true SNPs.

@@ Line 9: / Line 9: @@
 ==Research Findings==
 * Using the Michigan State University rice genome version 6.1 annotation (11) as a reference, the researchers mapped the bisulfite-converted reads to the genome using BS Seeker (12). From these alignments we determined SNPs between the NPB and 93–11 parental plants. Although there is ambiguity because of cytosine conversion in bisulfite-treated DNA, we can distinguish conversions from mutations by analyzing the sequence of the complimentary base of each cytosine, and, as a result, were able to call all possible types of SNPs. For our SNP calls we required that we had at least three reads on each strand and over 90% agreement between them before making a call (Fig. 1A). Using these SNP calls we reconstructed two genomes, one for each parent. We then mapped F1 reads against these inferred genomes to identify the parental chromosome from which they were derived.
 [[File:IC4R004-Epigenomic-2012-22778444-1.png |center |thumb |927px |'''Figure 1. SNPs between NPB and 93–11. (A) Identification of SNPs required a minimum of three reads per strand in each parental strain. Although there is ambiguity because of cytosine conversion in bisulfite-treated DNA, we can resolve this by considering the sequence of the reverse strand, and, as a result, are able to call all possible types of SNPs. In the example, X indicates the location of the SNP and reads represent, in order, NPB forward strand, NPB reverse strand, 93–11 forward strand, and 93–11 reverse strand. (B) SNP distribution over the genome is plotted on the bottom in windows of 10,000 bases. Each color change represents a new chromosome. The top line represents the density of genes.''']]
 * We found one SNP every 253 base pairs between the two parental strains (Fig. 1B), representing a level of divergence that is similar to previous japonica–indica comparisons (13). Not surprisingly, regions enriched with genes had lower densities of SNPs. Of these SNPs, the most frequent were C to T changes, or their complements corresponding to G to A mutations on the opposite strand. Combined, these types of mutations accounted for 73% of the SNPs that were found between varieties (Fig. 1C), which represents a threefold enrichment compared to the expected value if all types of mutations were equally likely. We also observe that methylated cytosines mutate more than three times more frequently than nonmethylated cytosines(Fig. 1D), and they most often mutate to thymines. This is consistent with previous observations that the deamination of methylcytosine, which results in thymine, occurs at higher rates than other spontaneous mutations (14–16).
 [[File:IC4R004-Epigenomic-2012-22778444-2.png |center |thumb |927px |'''Figure 1. (C) The fraction of the various types of SNPs found in NPB versus 93–11 is shown with the CT and GA SNPs grouped together, and they represent nearly 74% of all SNPs. (D) The relative frequency of each SNP between NPB and 93–11 is shown per 1,000 bases. Methyl-Cs convert to Ts at a rate that is nearly five times that of unmethylated Cs.''']]

@@ Line 9: / Line 9: @@
 ==Research Findings==
 * Using the Michigan State University rice genome version 6.1 annotation (11) as a reference, the researchers mapped the bisulfite-converted reads to the genome using BS Seeker (12). From these alignments we determined SNPs between the NPB and 93–11 parental plants. Although there is ambiguity because of cytosine conversion in bisulfite-treated DNA, we can distinguish conversions from mutations by analyzing the sequence of the complimentary base of each cytosine, and, as a result, were able to call all possible types of SNPs. For our SNP calls we required that we had at least three reads on each strand and over 90% agreement between them before making a call (Fig. 1A). Using these SNP calls we reconstructed two genomes, one for each parent. We then mapped F1 reads against these inferred genomes to identify the parental chromosome from which they were derived.
-[[File:IC4R001-Epigenomic-2008-18263775-1.png |center |thumb |1000px |'''Figure 1. Overview of Epigenetic Modifications of Rice Chromosomes 4 and 10.''']]
+[[File:IC4R004-Epigenomic-2012-22778444-1.png |center |thumb |927px |'''Figure 1. SNPs between NPB and 93–11. (A) Identification of SNPs required a minimum of three reads per strand in each parental strain. Although there is ambiguity because of cytosine conversion in bisulfite-treated DNA, we can resolve this by considering the sequence of the reverse strand, and, as a result, are able to call all possible types of SNPs. In the example, X indicates the location of the SNP and reads represent, in order, NPB forward strand, NPB reverse strand, 93–11 forward strand, and 93–11 reverse strand. (B) SNP distribution over the genome is plotted on the bottom in windows of 10,000 bases. Each color change represents a new chromosome. The top line represents the density of genes.''']]
 * We found one SNP every 253 base pairs between the two parental strains (Fig. 1B), representing a level of divergence that is similar to previous japonica–indica comparisons (13). Not surprisingly, regions enriched with genes had lower densities of SNPs. Of these SNPs, the most frequent were C to T changes, or their complements corresponding to G to A mutations on the opposite strand. Combined, these types of mutations accounted for 73% of the SNPs that were found between varieties (Fig. 1C), which represents a threefold enrichment compared to the expected value if all types of mutations were equally likely. We also observe that methylated cytosines mutate more than three times more frequently than nonmethylated cytosines(Fig. 1D), and they most often mutate to thymines. This is consistent with previous observations that the deamination of methylcytosine, which results in thymine, occurs at higher rates than other spontaneous mutations (14–16).
-[[File:IC4R001-Epigenomic-2008-18263775-1.png |center |thumb |1000px |'''Figure 1. Overview of Epigenetic Modifications of Rice Chromosomes 4 and 10.''']]
+[[File:IC4R004-Epigenomic-2012-22778444-2.png |center |thumb |927px |'''Figure 1. (C) The fraction of the various types of SNPs found in NPB versus 93–11 is shown with the CT and GA SNPs grouped together, and they represent nearly 74% of all SNPs. (D) The relative frequency of each SNP between NPB and 93–11 is shown per 1,000 bases. Methyl-Cs convert to Ts at a rate that is nearly five times that of unmethylated Cs.''']]
 * Next, the researchers compared mutation rates between the parental strains and the hybrids. Not surprisingly, we found that SNPs were nearly absent across generations. If, as is likely, the few SNPs we call are false positives, we estimate the false discovery rate of our SNP calls to be less than 0.0003%.
 * The researchers also found a small number of SNPs between the reference Michigan State University O. sativa japonica genome and our as- sembled version (1∕90651 bp; 2,487 total) (Dataset S2), and considerably more SNPs by comparing our assembled version of the 93–11 genome with the reference Beijing Genome Institute’s O. sativa ssp. indica genome (1∕51 bp), consistent with the fact that the 93–11 assembly is not as complete as that of the NPB cultivar (17). The fact that we only align reads that map to a unique posi- tion in the genome, have at most three mismatches with the reference genome, do not permit insertions or deletions, and re- quire a minimum of three reads on both strands, allowed us to only call bases over about half of the genome in both NPB and 93–11, and thus we are reporting only a subset of all true SNPs.

@@ Line 9: / Line 9: @@
 ==Research Findings==
 * Using the Michigan State University rice genome version 6.1 annotation (11) as a reference, the researchers mapped the bisulfite-converted reads to the genome using BS Seeker (12). From these alignments we determined SNPs between the NPB and 93–11 parental plants. Although there is ambiguity because of cytosine conversion in bisulfite-treated DNA, we can distinguish conversions from mutations by analyzing the sequence of the complimentary base of each cytosine, and, as a result, were able to call all possible types of SNPs. For our SNP calls we required that we had at least three reads on each strand and over 90% agreement between them before making a call (Fig. 1A). Using these SNP calls we reconstructed two genomes, one for each parent. We then mapped F1 reads against these inferred genomes to identify the parental chromosome from which they were derived.
 * We found one SNP every 253 base pairs between the two parental strains (Fig. 1B), representing a level of divergence that is similar to previous japonica–indica comparisons (13). Not surprisingly, regions enriched with genes had lower densities of SNPs. Of these SNPs, the most frequent were C to T changes, or their complements corresponding to G to A mutations on the opposite strand. Combined, these types of mutations accounted for 73% of the SNPs that were found between varieties (Fig. 1C), which represents a threefold enrichment compared to the expected value if all types of mutations were equally likely. We also observe that methylated cytosines mutate more than three times more frequently than nonmethylated cytosines(Fig. 1D), and they most often mutate to thymines. This is consistent with previous observations that the deamination of methylcytosine, which results in thymine, occurs at higher rates than other spontaneous mutations (14–16).
 * Next, the researchers compared mutation rates between the parental strains and the hybrids. Not surprisingly, we found that SNPs were nearly absent across generations. If, as is likely, the few SNPs we call are false positives, we estimate the false discovery rate of our SNP calls to be less than 0.0003%.
 * The researchers also found a small number of SNPs between the reference Michigan State University O. sativa japonica genome and our as- sembled version (1∕90651 bp; 2,487 total) (Dataset S2), and considerably more SNPs by comparing our assembled version of the 93–11 genome with the reference Beijing Genome Institute’s O. sativa ssp. indica genome (1∕51 bp), consistent with the fact that the 93–11 assembly is not as complete as that of the NPB cultivar (17). The fact that we only align reads that map to a unique posi- tion in the genome, have at most three mismatches with the reference genome, do not permit insertions or deletions, and re- quire a minimum of three reads on both strands, allowed us to only call bases over about half of the genome in both NPB and 93–11, and thus we are reporting only a subset of all true SNPs.

@@ Line 8: / Line 8: @@
 * BS-seq libraries were made from genomic DNA isolated from leaf tissues of NPB, 93–11, and their two reciprocal hybrid offspring, 93–11 × NPB and NPB × 93–11, by a previously published method using pre-methylated Illumina adapters
 ==Research Findings==
+* Using the Michigan State University rice genome version 6.1 annotation (11) as a reference, the researchers mapped the bisulfite-converted reads to the genome using BS Seeker (12). From these alignments we determined SNPs between the NPB and 93–11 parental plants. Although there is ambiguity because of cytosine conversion in bisulfite-treated DNA, we can distinguish conversions from mutations by analyzing the sequence of the complimentary base of each cytosine, and, as a result, were able to call all possible types of SNPs. For our SNP calls we required that we had at least three reads on each strand and over 90% agreement between them before making a call (Fig. 1A). Using these SNP calls we reconstructed two genomes, one for each parent. We then mapped F1 reads against these inferred genomes to identify the parental chromosome from which they were derived.
 ==Labs working on this Project==

@@ Line 5: / Line 5: @@
 * In this project, the researchers generated integrative maps of whole-genome cytosine methylation profiles [bisulfite sequencing (BS-seq)] and transcriptional profiles (RNA-seq), to characterize two rice subspecies, Oryza sativa spp japonica [Nipponbare (NPB)] and Oryza sativa spp indica (93–11) and their two reciprocal hybrid off- spring. Using a combination of BS-seq, RNA-seq, and siRNA-seq, we were able to generate allele-specific patterns of methylation and transcription in the hybrids, and thus directly measure the degree to which these are altered between the corresponding parental and F1 chromosomes.
 ==Plant Materials & Treatment==
-* Oryza sativa ssp. japonica (NPB), O.sativa ssp. indica (93–11), and their reciprocal cross F1s (NPB × 93–11 and 93–11 × NPB) were used in this study. Rice seeds were surface sterilized with 40% (vol∕vol) sodium hypochlorite solution and transferred on 1∕2 Murashige and Skoog medium. After germination, rice seedlings were trans- planted into soil and grown at 26 °C∕20 °C under a 10-h light/14-h dark cycle in growth chamber. Fully expanded leaves from 6-wk-old plants were collected for library construction.
+* Oryza sativa ssp. japonica (NPB), O.sativa ssp. indica (93–11), and their reciprocal cross F1s (NPB × 93–11 and 93–11 × NPB) were used in this study. Rice seeds were surface sterilized with 40% (vol/vol) sodium hypochlorite solution and transferred on 1∕2 Murashige and Skoog medium. After germination, rice seedlings were trans- planted into soil and grown at 26 °C/20 °C under a 10-h light/14-h dark cycle in growth chamber. Fully expanded leaves from 6-wk-old plants were collected for library construction.
 * BS-seq libraries were made from genomic DNA isolated from leaf tissues of NPB, 93–11, and their two reciprocal hybrid offspring, 93–11 × NPB and NPB × 93–11, by a previously published method using pre-methylated Illumina adapters